#ifndef __CINT__ #include "RooGlobalFunc.h" #endif #include "RooWorkspace.h" #include "RooRealVar.h" #include "RooAbsPdf.h" #include "RooGenericPdf.h" #include "RooArgList.h" #include "RooDataSet.h" #include "RooHist.h" #include "RooDataHist.h" #include "RooHistPdf.h" #include "RooPlot.h" #include "RooGlobalFunc.h" #include "RooGaussian.h" #include "RooConstVar.h" #include "RooTruthModel.h" #include "RooFormulaVar.h" #include "RooRealSumPdf.h" #include "RooPolyVar.h" #include "RooProduct.h" #include "RooFitResult.h" #include "TH1.h" #include "TCanvas.h" #include "TAxis.h" #include "RooStats/HistFactory/PiecewiseInterpolation.h" #include "RooStats/ModelConfig.h" #include "RooStats/BernsteinCorrection.h" using namespace RooFit; using namespace RooStats; void histo_maker() { //Open data and Monte Carlo ntuple TChain* dataSVSV=new TChain("Tuple"); dataSVSV->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/SVSV/Down/SVSV_data.root"); dataSVSV->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/SVSV/Up/SVSV_data.root"); TChain* chain_Hbb=new TChain("Tuple"); chain_Hbb->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/Hbb.root"); TChain* chain_Zbb=new TChain("Tuple"); chain_Zbb->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/Zbb.root"); TChain* chain_Zcc=new TChain("Tuple"); chain_Zcc->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/Zcc.root"); TChain* chain_Hcc=new TChain("Tuple"); chain_Hcc->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/Hcc.root"); TChain* chain_QCD_cc=new TChain("Tuple"); //for checks chain_QCD_cc->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVjet_Data/cc.root"); TChain* chain_QCD_cc_svsv=new TChain("Tuple"); //for checks chain_QCD_cc_svsv->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/cc.root"); chain_QCD_cc_svsv->Add("/home/user/Scrivania/Tesi_Hbb_Hcc/Transfer function/SVSV_MC/bb.root"); //Define interesting variables Double_t jet0_pe, jet0_px,jet0_py,jet0_pz,jet0_pt,jet0_eta; Double_t jet0_svtag, jet0_efb,jet0_efc; Float_t jet0_bdt0,jet0_bdt1,jet0_pb,jet0_pc,jet0_pq; Double_t jet1_pe,jet1_px,jet1_py, jet1_pz,jet1_pt,jet1_eta; Double_t jet1_svtag, jet1_efb,jet1_efc; Float_t jet1_bdt0,jet1_bdt1,jet1_pb,jet1_pc,jet1_pq; Double_t dijet_m,dijet_dphi,weight; //Open branch address chain_Hbb->SetBranchAddress("jet0_pb", &jet0_pb); chain_Hbb->SetBranchAddress("jet0_pc", &jet0_pc); chain_Hbb->SetBranchAddress("jet0_pq", &jet0_pq); chain_Hbb->SetBranchAddress("jet1_pb", &jet1_pb); chain_Hbb->SetBranchAddress("jet1_pc", &jet1_pc); chain_Hbb->SetBranchAddress("jet1_pq", &jet1_pq); chain_Hbb->SetBranchAddress("dijet_m", &dijet_m); chain_Hcc->SetBranchAddress("jet0_pb", &jet0_pb); chain_Hcc->SetBranchAddress("jet0_pc", &jet0_pc); chain_Hcc->SetBranchAddress("jet0_pq", &jet0_pq); chain_Hcc->SetBranchAddress("jet1_pb", &jet1_pb); chain_Hcc->SetBranchAddress("jet1_pc", &jet1_pc); chain_Hcc->SetBranchAddress("jet1_pq", &jet1_pq); chain_Hcc->SetBranchAddress("dijet_m", &dijet_m); chain_Zbb->SetBranchAddress("jet0_pb", &jet0_pb); chain_Zbb->SetBranchAddress("jet0_pc", &jet0_pc); chain_Zbb->SetBranchAddress("jet0_pq", &jet0_pq); chain_Zbb->SetBranchAddress("jet1_pb", &jet1_pb); chain_Zbb->SetBranchAddress("jet1_pc", &jet1_pc); chain_Zbb->SetBranchAddress("jet1_pq", &jet1_pq); chain_Zbb->SetBranchAddress("dijet_m", &dijet_m); chain_Zcc->SetBranchAddress("jet0_pb", &jet0_pb); chain_Zcc->SetBranchAddress("jet0_pc", &jet0_pc); chain_Zcc->SetBranchAddress("jet0_pq", &jet0_pq); chain_Zcc->SetBranchAddress("jet1_pb", &jet1_pb); chain_Zcc->SetBranchAddress("jet1_pc", &jet1_pc); chain_Zcc->SetBranchAddress("jet1_pq", &jet1_pq); chain_Zcc->SetBranchAddress("dijet_m", &dijet_m); dataSVSV->SetBranchAddress("jet0_pb", &jet0_pb); dataSVSV->SetBranchAddress("jet0_pc", &jet0_pc); dataSVSV->SetBranchAddress("jet0_pq", &jet0_pq); dataSVSV->SetBranchAddress("jet1_pb", &jet1_pb); dataSVSV->SetBranchAddress("jet1_pc", &jet1_pc); dataSVSV->SetBranchAddress("jet1_pq", &jet1_pq); dataSVSV->SetBranchAddress("dijet_m", &dijet_m); chain_QCD_cc->SetBranchAddress("weight", &weight); chain_QCD_cc->SetBranchAddress("jet0_pb", &jet0_pb); chain_QCD_cc->SetBranchAddress("jet0_pc", &jet0_pc); chain_QCD_cc->SetBranchAddress("jet0_pq", &jet0_pq); chain_QCD_cc->SetBranchAddress("jet1_pb", &jet1_pb); chain_QCD_cc->SetBranchAddress("jet1_pc", &jet1_pc); chain_QCD_cc->SetBranchAddress("jet1_pq", &jet1_pq); chain_QCD_cc->SetBranchAddress("dijet_m", &dijet_m); chain_QCD_cc->SetBranchAddress("jet1_svtag", &jet1_svtag); chain_QCD_cc->SetBranchAddress("jet0_svtag", &jet0_svtag); chain_QCD_cc_svsv->SetBranchAddress("weight", &weight); chain_QCD_cc_svsv->SetBranchAddress("jet0_pb", &jet0_pb); chain_QCD_cc_svsv->SetBranchAddress("jet0_pc", &jet0_pc); chain_QCD_cc_svsv->SetBranchAddress("jet0_pq", &jet0_pq); chain_QCD_cc_svsv->SetBranchAddress("jet1_pb", &jet1_pb); chain_QCD_cc_svsv->SetBranchAddress("jet1_pc", &jet1_pc); chain_QCD_cc_svsv->SetBranchAddress("jet1_pq", &jet1_pq); chain_QCD_cc_svsv->SetBranchAddress("dijet_m", &dijet_m); //define histo ranges and n bins int x= 38; //n bins int min_bin=60; //min inv mass range int max_bin=250; //max inv mass range //define SR and CR in Monte Carlo for svjet //TH1F* MCbb_SR_svjet=new TH1F("MCbb_SR_svjet","MCbb_SR_svjet",x,min_bin,max_bin); //TH1F* MCbb_CR_svjet=new TH1F("MCbb_CR_svjet","MCbb_CR_svjet",x,min_bin,max_bin); TH1F* MCbb_SR_svsv=new TH1F("MCbb_SR_svsv","MCbb_SR_svsv",x,min_bin,max_bin); TH1F* MCbb_CR_svsv=new TH1F("MCbb_CR_svsv","MCbb_CR_svsv",x,min_bin,max_bin); TH1F* TF_MC=new TH1F("TF_MC","TF_MC",x,min_bin,max_bin); TH1F* TF_ratio=new TH1F("TF_ratio","TF_ratio",x,min_bin,max_bin); TH1F* TF_data=new TH1F("TF_data","TF_data",x,min_bin,max_bin); TH1F* NON_GUARDARE_data_SR=new TH1F("NON_GUARDARE_data_SR","NON_GUARDARE_data_SR",x,min_bin,max_bin); TH1F* data_SR=new TH1F("data_SR","data_SR",x,min_bin,max_bin); TH1F* data_CR=new TH1F("data_CR","data_CR",x,min_bin,max_bin); TH1F* Higgsbb_SR=new TH1F("Higgsbb_SR","Higgsbb_SR",x,min_bin,max_bin); TH1F* Zbb_SR=new TH1F("Zbb_SR","Zbb_SR",x,min_bin,max_bin); TH1F* Higgscc_SR=new TH1F("Higgscc_SR","Higgscc_SR",x,min_bin,max_bin); TH1F* signal=new TH1F("signal","signal",x,min_bin,max_bin); TH1F* Zcc_SR=new TH1F("Zcc_SR","Zcc_SR",x,min_bin,max_bin); TH1F* Data_SR_correction=new TH1F("Data_SR_correction","Data_SR_correction",x,min_bin,max_bin); TH1F* Data_CR_correction=new TH1F("Data_CR_correction","Data_CR_correction",x,min_bin,max_bin); //TH1F* TF_MC_from_funct=new TH1F("TF_MC_from_funct","TF_MC_from_funct",x,min_bin,max_bin); int nEntries_data = (int)dataSVSV->GetEntries(); //int nEntries_QCD = (int)chain_QCD_cc->GetEntries(); int nEntries_QCD_svsv = (int)chain_QCD_cc_svsv->GetEntries(); int nEntries_Hbb = (int)chain_Hbb->GetEntries(); int nEntries_Hcc = (int)chain_Hcc->GetEntries(); int nEntries_Zbb = (int)chain_Zbb->GetEntries(); int nEntries_Zcc = (int)chain_Zcc->GetEntries(); //FOR LOOPS ON DATA AND MONTE CARLO int count_number_SVSV_SR=0; cout << nEntries_data << endl; for (int i=0;i<1000000; i++) //nEntries_data { dataSVSV->GetEntry(i); if ( jet0_pc>0.15 && jet1_pc>0.15) {if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);}continue;} } for (int i=1000000;iGetEntry(i); //Fill Signal region SVSV data if ( jet0_pc>0.15 && jet1_pc>0.15) {if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);data_SR->Fill(dijet_m*0.001);count_number_SVSV_SR++;}continue;} if ((jet1_pc>0.15 && jet0_pc<0.02)||(jet0_pc>0.15 && jet1_pc<0.02)) {if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);}} } //TF READY //MONTE CARLO DISTRIBUTIONS SIGNAL REGION for (int i=0;iGetEntry(i); if ( jet0_pc>0.15 && jet1_pc>0.15 ) { if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);} } } for (int i=0;iGetEntry(i); if ( jet0_pc>0.15 && jet1_pc>0.15 ) { if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);signal->Fill(dijet_m*0.001);} } } for (int i=0;iGetEntry(i); if ( jet0_pc>0.15 && jet1_pc>0.15 ) { if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);} } } for (int i=0;iGetEntry(i); if ( jet0_pc>0.15 && jet1_pc>0.15) { if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001);} } } for (int y=0; yGetEntry(y); if ( jet0_pc>0.15 && jet1_pc>0.15) { if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001,weight);} } } for (int y=0; yGetEntry(y); if ( jet0_pc>0.15 && jet1_pc>0.15 ) {continue;} if ((jet1_pc>0.15 && jet0_pc<0.02)||(jet0_pc>0.15 && jet1_pc<0.02)) { if(dijet_m*0.001>min_bin && dijet_m*0.001Fill(dijet_m*0.001,weight);} } } //TF_MC->Draw(); double luminosity=1576.86;//1644.16;//1644.16; //pb^-1 1576.86 double acceptance_bb=0.051;//0.051; double acceptance_cc=0.064; double cross_section_cc=1.647;//29.6; //pb^-1 double cross_section_bb=31.8; //pb^-1 double efficiency_sel_bb=0.11; double efficiency_sel_cc=0.014; double N_H_bb_ev=acceptance_bb*luminosity*cross_section_bb*efficiency_sel_bb; double N_H_cc_ev=acceptance_cc*luminosity*cross_section_cc*efficiency_sel_cc; double acceptance_Zbb=0.161;//0.051; double acceptance_Zcc=0.213; double cross_section_Zcc=4154;//29.6; //pb^-1 double cross_section_Zbb=5307.0; //pb^-1 double efficiency_sel_Zbb=0.14; double efficiency_sel_Zcc=0.011; double N_Z_bb_ev=acceptance_Zbb*luminosity*cross_section_Zbb*efficiency_sel_Zbb; double N_Z_cc_ev=acceptance_Zcc*luminosity*cross_section_Zcc*efficiency_sel_Zcc; cout << "Number of data in Control region " << data_CR->Integral(0,320000) << endl; cout << "Number of data in Signal region " << count_number_SVSV_SR << endl; cout << "Number of Hbb in Signal region " << Higgsbb_SR->GetEntries() << endl; cout << "Number of Hcc in Signal region " << Higgscc_SR->GetEntries() << endl; cout << "Number of Zbb in Signal region " << Zbb_SR->GetEntries() << endl; cout << "Number of Zcc in Signal region " << Zcc_SR->GetEntries() << endl; cout << endl; cout << endl; cout << "SM number of Higgs in bb " <GetEntries()/nEntries_Hbb<< endl; cout << "Efficiency Hcc in Signal region " << Higgscc_SR->GetEntries()/nEntries_Hcc << endl; cout << "Efficiency Zbb in Signal region " << Zbb_SR->GetEntries()/nEntries_Zbb<< endl; cout << "Efficiency Zcc in Signal region " << Zcc_SR->GetEntries()/nEntries_Zcc << endl; cout << endl; cout << endl; cout << "CONTEGGIO NUMERO EVENTI STANDARD MODEL " << endl; cout << "Number of Hbb in Signal region " << (Higgsbb_SR->GetEntries()/nEntries_Hbb)*N_H_bb_ev << endl; cout << "Number of Hcc in Signal region " << (Higgscc_SR->GetEntries()/nEntries_Hcc)*N_H_cc_ev << endl; cout << "Number of Zbb in Signal region " << (Zbb_SR->GetEntries()/nEntries_Zbb)*N_Z_bb_ev << endl; cout << "Number of Zcc in Signal region " << (Zcc_SR->GetEntries()/nEntries_Zcc)*N_Z_cc_ev << endl; Data_SR_correction->Sumw2(kTRUE); Data_SR_correction->Divide(data_CR); //MCbb_SR_svjet->Sumw2(kTRUE); //MCbb_SR_svjet->Divide(MCbb_CR_svjet); MCbb_SR_svsv->Sumw2(kTRUE); MCbb_SR_svsv->Divide(MCbb_CR_svsv); cout << "Ratio Data_SR_correction " << Data_SR_correction->GetBinContent(1)<Sumw2(kTRUE); cout << "error Data_SR_correction " << Data_SR_correction->GetBinError(1) << endl; cout << "Ratio MCbb_SR_svsv " << MCbb_SR_svsv->GetBinContent(1)<Sumw2(kTRUE); cout << "error MCbb_SR_svsv " << MCbb_SR_svsv->GetBinError(1) << endl; Data_SR_correction->Sumw2(kTRUE); Data_SR_correction->Divide(MCbb_SR_svsv); cout << "Ratio Data_SR_correction" << Data_SR_correction->GetBinContent(1)<Sumw2(kTRUE); cout << "error Data_SR_correction" << Data_SR_correction->GetBinError(1) << endl; RooRealVar y("y","y",min_bin,max_bin); RooDataHist Zbb("Zbb","Zbb",y,Import(*Zbb_SR)); RooDataHist Zcc("Zcc","Zcc",y,Import(*Zcc_SR)); RooDataHist Hbb("Hbb","Hbb",y,Import(*Higgsbb_SR)); RooDataHist Hcc("Hcc","Hcc",y,Import(*Higgscc_SR)); RooDataHist DataCR("DataCR","DataCR",y,Import(*data_CR)); RooHistPdf Hbb_pdf("Hbb_pdf","Hbb_pdf",y,Hbb,0); RooHistPdf Zbb_pdf("Zbb_pdf","Zbb_pdf",y,Zbb,0); RooHistPdf Hcc_pdf("Hcc_pdf","Hcc_pdf",y,Hcc,0); RooHistPdf Zcc_pdf("Zcc_pdf","Zcc_pdf",y,Zcc,0); RooHistPdf DataCR_pdf("DataCR_pdf","DataCR_pdf",y,DataCR); int N_data=2036834-9376-11986-1.8-14; RooRealVar n_bkg_SR("n_bkg_SR", "n_bkg_SR",N_data,0,5000000 );//,0,20000000000 RooRealVar n_Zbb("n_Zbb", "n_Zbb",9376); // ,0,20000000 RooRealVar n_Hbb("n_Hbb", "n_Hbb",14); // ,-100000,100000 RooRealVar n_Zcc("n_Zcc", "n_Zcc",11986); RooRealVar n_Hcc("n_Hcc", "n_Hcc",1.8); //CORRECTION FUNCTION +1 sigma RooRealVar x0("x0","x0",1.03665); //,0.,10. RooRealVar x1("x1","x1", -0.00656052); //,-0.1,0. RooRealVar x2("x2","x2",2.77415e-05); //,1.e-03 RooArgSet set1(y,x0,x1,x2); RooGenericPdf TF_func("TF_func","TF_func","x0+y*x1+y*y*x2", RooArgSet(set1)); RooRealVar x0_1("x0_1","x0_1",4.98278); RooRealVar x1_1("x1_1","x1_1",4.72565 ); RooRealVar x2_1("x2_1","x2_1",8.11861 ); RooRealVar x3_1("x3_1","x3_1",7.00158); RooRealVar x4_1("x4_1","x4_1",6.85781); RooBernstein bg_bern_after_fit_1("bg_bern_after_fit_1","bg_bern_after_fit_1",y,RooArgList(x0_1,x1_1,x2_1,x3_1,x4_1)); y.setBins(38); RooProdPdf bkgProd_1("bkgProd_1","bkgProd_1*TF_func",bg_bern_after_fit_1,TF_func); RooBinSamplingPdf binSampler_after_fit_1("Name", "Plot Label", y, bkgProd_1); RooProdPdf bkgProd_fin_1("bkgProd_fin_1","bkgProd_fin_1*DataCR_pdf",binSampler_after_fit_1,DataCR_pdf); RooAddPdf model_0("model_0","model_0",RooArgList(bkgProd_fin_1,Hbb_pdf,Zcc_pdf,Zbb_pdf,Hcc_pdf),RooArgList(n_bkg_SR,n_Hbb,n_Zcc,n_Zbb,n_Hcc)); RooDataHist *mod_final_1=model_0.generateBinned(y,50000000); mod_final_1->SetName("mod_final_1"); RooHistPdf mod_final_1_pdf("mod_final_1_pdf","mod_final_1_pdf",y,*mod_final_1,0); RooRealVar x0_m1("x0_m1","x0_m1",4.98278-0.574746); RooRealVar x1_m1("x1_m1","x1_m1", 4.72565 - 0.554297 ); RooRealVar x2_m1("x2_m1","x2_m1",8.11861 - 0.970454); RooRealVar x3_m1("x3_m1","x3_m1",7.00158 - 0.924254); RooRealVar x4_m1("x4_m1","x4_m1",6.85781 - 0.817415); RooBernstein bg_bern_after_fit_m1("bg_bern_after_fit_m1","bg_bern_after_fit_m1",y,RooArgList(x0_m1,x1_m1,x2_m1,x3_m1,x4_m1)); y.setBins(38); RooProdPdf bkgProd_m1("bkgProd_m1","bkgProd_m1*TF_func",bg_bern_after_fit_m1,TF_func); RooBinSamplingPdf binSampler_after_fit_m1("Name", "Plot Label", y, bkgProd_m1); RooProdPdf bkgProd_fin_m1("bkgProd_fin_m1","bkgProd_fin_m1*DataCR_pdf",binSampler_after_fit_m1,DataCR_pdf); RooAddPdf model_m1("model_m1","model_m1",RooArgList(bkgProd_fin_m1,Hbb_pdf,Zcc_pdf,Zbb_pdf,Hcc_pdf),RooArgList(n_bkg_SR,n_Hbb,n_Zcc,n_Zbb,n_Hcc)); RooDataHist *mod_final_m1=model_m1.generateBinned(y,50000000); mod_final_m1->SetName("mod_final_m1"); RooHistPdf mod_final_m1_pdf("mod_final_m1_pdf","mod_final_m1_pdf",y,*mod_final_m1,0); RooRealVar x0_p1("x0_p1","x0_p1",4.98278+0.574746); RooRealVar x1_p1("x1_p1","x1_p1", 4.72565 + 0.554297 ); RooRealVar x2_p1("x2_p1","x2_p1",8.11861 + 0.970454); RooRealVar x3_p1("x3_p1","x3_p1",7.00158 + 0.924254); RooRealVar x4_p1("x4_p1","x4_p1",6.85781 + 0.817415); RooBernstein bg_bern_after_fit_p1("bg_bern_after_fit_p1","bg_bern_after_fit_p1",y,RooArgList(x0_p1,x1_p1,x2_p1,x3_p1,x4_p1)); y.setBins(38); RooProdPdf bkgProd_p1("bkgProd_p1","bkgProd_p1*TF_func",bg_bern_after_fit_p1,TF_func); RooBinSamplingPdf binSampler_after_fit_p1("Name", "Plot Label", y, bkgProd_p1); RooProdPdf bkgProd_fin_p1("bkgProd_fin_p1","bkgProd_fin_p1*DataCR_pdf",binSampler_after_fit_p1,DataCR_pdf); RooAddPdf model_p1("model_p1","model_p1",RooArgList(bkgProd_fin_p1,Hbb_pdf,Zcc_pdf,Zbb_pdf,Hcc_pdf),RooArgList(n_bkg_SR,n_Hbb,n_Zcc,n_Zbb,n_Hcc)); RooDataHist *mod_final_p1=model_p1.generateBinned(y,50000000); mod_final_p1->SetName("mod_final_p1"); RooHistPdf mod_final_p1_pdf("mod_final_p1_pdf","mod_final_p1_pdf",y,*mod_final_p1,0); TH1 *HighHistForSyst4 = mod_final_p1->createHistogram("HighHistForSyst4",y,Binning(38)); TH1 *LowHistForSyst4 = mod_final_m1->createHistogram("LowHistForSyst4",y,Binning(38)); TH1 *bkg_model_full = mod_final_1->createHistogram("bkg_model_full",y,Binning(38)); HighHistForSyst4->Scale(2036832./HighHistForSyst4->Integral()); LowHistForSyst4->Scale(2036832./LowHistForSyst4->Integral()); bkg_model_full->Scale(2036832./bkg_model_full->Integral()); HighHistForSyst4->Sumw2(kTRUE); LowHistForSyst4->Sumw2(kTRUE); bkg_model_full->Sumw2(kTRUE); Higgsbb_SR->Scale(n_Hbb.getVal()/Higgsbb_SR->Integral()); Zbb_SR->Scale(n_Zbb.getVal()/Zbb_SR->Integral()); signal->Scale(n_Hcc.getVal()/Higgscc_SR->Integral()); Zcc_SR->Scale(n_Zcc.getVal()/Zcc_SR->Integral()); TFile* SVSV_analysis= new TFile("analysis.root","RECREATE"); signal->Write(); bkg_model_full->Write(); HighHistForSyst4->Write(); LowHistForSyst4->Write(); TCanvas *c1 = new TCanvas("c1"); RooPlot * modello_plot =y.frame(Title("modello_plot")); //modello->plotOn(modello_plot) ; mod_final_p1->Draw() ; /* Number of data in Control region 3.09383e+06 Number of data in Signal region 19673732 Number of Hbb in Signal region 21605 Number of Hcc in Signal region 987 Number of Zbb in Signal region 26982 Number of Zcc in Signal region 365 Efficiency data in Signal region 19673732 Efficiency Hbb in Signal region 0.940575 Efficiency Hcc in Signal region 0.36994 Efficiency Zbb in Signal region 0.953899 Efficiency Zcc in Signal region 0.337338 CONTEGGIO NUMERO EVENTI STANDARD MODEL Number of Hbb in Signal region 246.287 Number of Hcc in Signal region 0.768334 Number of Zbb in Signal region 179928 Number of Zcc in Signal region 5177.23 */ }