void PlotAll_DATA_PREAPP()
{
gROOT->LoadMacro("invar_yield_ana_v9.C");
char dir[100];
char dir_corr[100];
bool GEN = false;
bool CORRECT = false;
bool REBIN = true;
bool MC = false;
bool multcrct =false;
bool seccrct = false;
bool cross = false;
bool loose = false;
bool jet_based_correction = true;
bool entire_range = true;
bool zerobin = false;
bool onetwothreebin = false;
bool evteffcorr = false;
bool full_eta =true;
bool oneoverpt=true;
bool minpt=true;
double eta_max;
float pt_max;
float ymin, ymax;
float ymin_r, ymax_r;
if(entire_range) {
if(MC) pt_max = 55, ymin = 3e-13, ymax = 9e1, ymin_r = 0.30, ymax_r = 1.65;
else pt_max = 100, ymin = 3e-14, ymax = 9e1, ymin_r = 0.6, ymax_r = 1.65;
} else {
if(MC) pt_max = 10, ymin = 5e-7, ymax = 2e1, ymin_r = 0.7, ymax_r = 1.4;
//else pt_max = 10, ymin = 5e-7, ymax = 5e1, ymin_r = 0.97, ymax_r = 1.03;
else pt_max = 10, ymin = 5e-7, ymax = 5e1, ymin_r = 0.5, ymax_r = 1.5;
}
if(full_eta) eta_max = 2.4;
else eta_max = 1.0;
char file1[100], file2[100], file3[100], file4[100], file5[100], file6[100], file7[100];
sprintf(file1,"../root_files/MB-C10-A20RR-TRKANASKIM-MERGED.root");
//sprintf(file1,"../root_files/MB-C10-MERGED.root");
//sprintf(file1,"../root_files/MB-C10-A20PR_proc0531.root");
//sprintf(file2,"../root_files/TrkHistMC_june04_10Mv2.root");
//sprintf(file3,"../root_files/TrkHistMC_june04_10Mv2.root");
sprintf(file2,"../root_files/TrkHistMC_june04_QCD_10Mv6.root");
sprintf(file3,"../root_files/TrkHistMC_june04_QCD_10Mv6.root");
sprintf(file4,"../root_files/TrkHistMCv10_QCD_Pt15v3.root");
sprintf(file5,"../root_files/TrkHistMCv11_qcdPt30.root");
sprintf(file6,"../root_files/TrkHistMCv11_qcdPt80.root");
sprintf(file7,"../root_files/TrkHistMCv11_qcdPt170v2.root");
MC = false, GEN = false, CORRECT = true, multcrct = true, seccrct = true, evteffcorr = true, zerobin = true, onetwothreebin = true;
sprintf(dir,"trackAna_STD");
sprintf(dir_corr,"trkEffAnalyzer");
invar_yield_ana_v9_data spec_rec_tight =
invar_yield_ana_v9_graph(file1,file2,file3,file4,file5,file6,file7,
dir,dir_corr,GEN,CORRECT,multcrct,seccrct,MC,jet_based_correction,
evteffcorr,zerobin,onetwothreebin,cross,oneoverpt,0,eta_max);
//sprintf(file1,"../root_files/TrkHistMC_june04_10Mv2.root");
//sprintf(file2,"../root_files/TrkHistMC_june04_10Mv2.root");
//sprintf(file3,"../root_files/TrkHistMC_june04_10Mv2.root");
sprintf(file1,"../root_files/TrkHistMC_june04_QCD_10Mv6.root");
sprintf(file2,"../root_files/TrkHistMC_june04_QCD_10Mv6.root");
sprintf(file3,"../root_files/TrkHistMC_june04_QCD_10Mv6.root");
sprintf(file4,"../root_files/TrkHistMCv10_QCD_Pt15v3.root");
sprintf(file5,"../root_files/TrkHistMCv11_qcdPt30.root");
sprintf(file6,"../root_files/TrkHistMCv11_qcdPt80.root");
sprintf(file7,"../root_files/TrkHistMCv11_qcdPt170v2.root");
// GEN MC
MC = true, GEN = true, CORRECT = false, multcrct = false, seccrct = false, evteffcorr = false, zerobin = false, onetwothreebin = false;
sprintf(dir,"preTrackAna");
sprintf(dir_corr,"trkEffAnalyzer");
invar_yield_ana_v9_data spec_rec_tight_MC =
invar_yield_ana_v9_graph(file1,file2,file3,file4,file5,file6,file7,
dir,dir_corr,GEN,CORRECT,multcrct,seccrct,MC,jet_based_correction,
evteffcorr,zerobin,onetwothreebin,cross,oneoverpt,0,eta_max);
//--------------------------- CMS measurement ---------------------------
TGraphErrors* cms_7000GeV = (TGraphErrors*) CMS_7TEV(1);
TF1 *f2 = new TF1("f2","[0]*(1+(sqrt(0.1396**2+x**2)-0.1396)/([1]*[2]))**(-[2])",0,6);
f2->SetParameters(2*TMath::Pi(),0.13,7.7);
cms_7000GeV->Fit(f2,"RN");
f2->SetLineColor(kRed);
//--------------------------------------------------- PYTHIA truth GEN------------------------
gROOT->LoadMacro("invar_yield_ana_GEN_v9.C");
char file1_GEN[100], file2_GEN[100], file3_GEN[100], file4_GEN[100], file5_GEN[100], file6_GEN[100], file7_GEN[100];
char file8_GEN[100], file9_GEN[100], file10_GEN[100];
sprintf(file1_GEN,"../root_files/TrkHistMC_june04_qcdPt0To15v3cFix.root");
sprintf(file2_GEN,"../root_files/TrkHistMC_june04_qcdPt15to20v2.root");
sprintf(file3_GEN,"../root_files/TrkHistMC_june04_qcdPt20to30v2.root");
sprintf(file4_GEN,"../root_files/TrkHistMC_june04_qcdPt30to50v2.root");
sprintf(file5_GEN,"../root_files/TrkHistMC_june04_qcdPt50to80v2.root");
sprintf(file6_GEN,"../root_files/TrkHistMC_june04_qcdPt80to120v2.root");
sprintf(file7_GEN,"../root_files/TrkHistMC_june04_qcdPt120to170v2.root");
sprintf(file8_GEN,"../root_files/TrkHistMC_june04_qcdPt170to230v2.root");
sprintf(file9_GEN,"../root_files/TrkHistMC_june04_qcdPt230to300v2.root");
sprintf(file10_GEN,"../root_files/TrkHistMC_june04_qcdPt300to380v2.root");
invar_yield_ana_GEN_v9_data gen_nsd_spect =
invar_yield_ana_GEN_v9_graph(file1_GEN,file2_GEN,file3_GEN,file4_GEN,
file5_GEN,file6_GEN,file7_GEN,file8_GEN,
file9_GEN,file10_GEN,true,117.,0,eta_max);
//--------------------------------------------------- PYTHIA truth GEN------------------------
double intLum;
if(REBIN) {
TH1D *dndpt_rec_tight = spec_rec_tight.hRInvX;
intLum = spec_rec_tight.integratedLum;
dndpt_rec_tight->Scale(1./intLum);
TH1D *dndpt_rec_tight_MC = spec_rec_tight_MC.hRInvX;
intLum = spec_rec_tight_MC.integratedLum;
dndpt_rec_tight_MC->Scale(1./intLum);
TH1D *dndpt_gen_nsd = gen_nsd_spect.hRInvX;
TGraphErrors *dndpt_gen_nsd_tg = gen_nsd_spect.RInvX;
} else {
}
// Basic canvas and dummy histogram
TCanvas *call = new TCanvas("call","call",510,670);
call->cd();
//if(entire_range) call->SetLogx();
// ------------------------------- Pad settings
pp1 = new TPad("p1","p1",0,0.34,1,1,0,0,0);
pp1->SetBottomMargin(0.0);
pp1->SetTopMargin(0.05*(1/0.72));
pp1->Draw();
pp1->cd();
pp1->SetNumber(1);
call->cd();
pp1_1 = new TPad("p1_1","p1_1",0,0.0,1,0.34,0,0,0);
pp1_1->SetTopMargin(0);
pp1_1->SetBottomMargin(0.14*(1/0.34));
pp1_1->Draw();
pp1_1->cd();
pp1_1->SetNumber(2);
// ----- pad 1
pp1->cd();
pp1->SetLogy();
if(entire_range) pp1->SetLogx();
Char_t xTitle[100],yTitle[100];
if(!cross) sprintf(yTitle,"Ed^{3}N/d^{3}p [GeV^{-2}c^{3}]");
else sprintf(yTitle,"Ed^{3}#sigma/d^{3}p [mb GeV^{-2}c^{3}]");
sprintf(xTitle,"p_{T} [GeV/c]");
TH1D *dum = GetDummyHist(pt_max,ymin,ymax,xTitle,yTitle);
dum->Draw();
dum->GetXaxis()->SetNdivisions(908);
dum->GetYaxis()->SetNdivisions(407);
TPad *call_pd = (TPad*)gPad;
call_pd->cd();
gPad->SetLogy();
if(minpt) {
dndpt_rec_tight->GetXaxis()->SetRange(3,dndpt_rec_tight->GetXaxis()->GetLast());
dndpt_rec_tight_MC->GetXaxis()->SetRange(3,dndpt_rec_tight_MC->GetXaxis()->GetLast());
dndpt_gen_nsd->GetXaxis()->SetRange(3,dndpt_gen_nsd->GetXaxis()->GetLast());
}
//plotting
th1Style1(cms_7000GeV,13,30,1.0,13,1,1,1);
//th1Style1(dndpt_gen_nsd,96,20,1,96,2,1,2);
th1Style1(dndpt_gen_nsd_tg,96,20,1,96,2,1,2);
//th1Style1(dndpt_rec_tight_MC,2,20,1.0,2,1.5,1,3);
th1Style1(dndpt_rec_tight,1,20,1.0,1,1.5,1,1);
if(entire_range) TLegend *leg2 = new TLegend(0.21,0.09,0.58,0.47);
else TLegend *leg2 = new TLegend(0.37,0.47,0.74,0.85);
leg2->SetMargin(0.3);
leg2->SetBorderSize(0);
leg2->SetFillColor(kWhite);
leg2->SetFillStyle(0);
leg2->SetTextSize(0.039);
leg2->SetMargin(0.3);
if(full_eta) leg2->SetHeader("pp #rightarrow h+X, 7 TeV DATA, |#eta|<2.4");
else leg2->SetHeader("pp #rightarrow h+X, 7 TeV DATA, |#eta|<1.0");
leg2->AddEntry(dndpt_rec_tight,"RECO (TRK+EVT corrected)","pl");
leg2->AddEntry(cms_7000GeV,"CMS 7 TeV |#eta|<2.4","pl");
leg2->AddEntry(dndpt_gen_nsd_tg,"PYTHIA D6T","l");
leg2->Draw();
TLatex * tex;
tex = new TLatex(0.69,0.88,"CMS Preliminary");
tex->SetTextSize(0.04);
tex->SetLineWidth(2);
tex->SetNDC();
tex->Draw();
// ------ pad2
pp1_1->cd();
if(entire_range) pp1_1->SetLogx();
sprintf(yTitle,"Ratio");
//TH1D *dum2 = GetDummyHist(-0.1,2.3,xTitle,yTitle);
TH1D *dum2 = GetDummyHist(pt_max,ymin_r,ymax_r,xTitle,yTitle);
dum2->GetYaxis()->SetLabelSize(0.07);
dum2->GetYaxis()->SetTitleSize(0.10);
dum2->GetYaxis()->SetTitleOffset(0.75);
dum2->GetXaxis()->SetLabelSize(0.09);
dum2->GetXaxis()->SetLabelOffset(0.05);
dum2->GetXaxis()->SetTitleSize(0.11);
dum2->GetXaxis()->SetTitleOffset(1.47);
dum2->GetYaxis()->SetNdivisions(405);
dum2->GetYaxis()-> SetDecimals(true);
dum2->Draw();
TLine *line = new TLine(0,1,pt_max,1);
line->SetLineStyle(2);
line->SetLineWidth(2);
line->Draw();
TH1D *dndpt_rec_tight_dum = (TH1D*) dndpt_rec_tight->Clone("dndpt_rec_tight_dum");
TH1D *dndpt_rec_tight_dum2 = (TH1D*) dndpt_rec_tight->Clone("dndpt_rec_tight_dum2");
TH1D *dndpt_rec_tight_MC_dum = (TH1D*) dndpt_rec_tight_MC->Clone("dndpt_rec_tight_MC_dum");
TH1D *dndpt_gen_nsd_dum = (TH1D*) dndpt_gen_nsd->Clone("dndpt_gen_nsd_dum");
TH1D *dndpt1_dum_div_cms7000 = (TH1D*) ratio_hist_to_func(dndpt_rec_tight_dum,
f2->GetParameter(0),
f2->GetParameter(1),
f2->GetParameter(2));
dndpt_gen_nsd_dum->Sumw2();
dndpt_rec_tight_dum->Sumw2();
dndpt_rec_tight_MC_dum->Sumw2();
dndpt_rec_tight_dum->Divide(dndpt_gen_nsd_dum);
//dndpt_rec_tight_dum2->Divide(dndpt_rec_tight_MC_dum);
if(minpt) {
dndpt_rec_tight_dum->GetXaxis()->SetRange(3,dndpt_rec_tight_dum->GetXaxis()->GetLast());
//dndpt_rec_tight_dum2->GetXaxis()->SetRange(3,dndpt_rec_tight_dum2->GetXaxis()->GetLast());
}
th1Style1(dndpt_rec_tight_dum,96,20,1,96,2,1,3);
//th1Style1(dndpt_rec_tight_dum2,2,24,1.0,2,1.5,1,1);
th1Style1(dndpt1_dum_div_cms7000,13,30,1.0,13,1.5,1,1);
// Print the values!!
//PrintXsection(dndpt_rec_tight,117.);
PrintXsection(dndpt_rec_tight,1.);
if(entire_range) {
if(full_eta) printCanvases(call,"spectra_DATA_PREAPP_fullEta",2,1);
else printCanvases(call,"spectra_DATA_PREAPP",2,1);
} else {
if(full_eta) printCanvases(call,"spectra_DATA_PREAPP_narrow_fullEta",1,1);
else printCanvases(call,"spectra_DATA_PREAPP_narrow",1,1);
}
}
void PlotPASFigFourAllSysErrorV2(bool save=false){
bool lowpT = false;
bool atlas = false;
float pt_min, pt_max;
float ymin, ymax;
float ymin_r, ymax_r;
//57.38
//double sigma_nsd = 57.38;
//double sigma_nsd = 71.3;
//double sigma_nsd = 57.38; // M6RR
double sigma_nsd = 56.47; // M6RR + PR9
if(lowpT){
pt_min = 0.5, pt_max = 7;
ymin = 6E-6, ymax = 5E+1;
ymin_r = 0.7, ymax_r = 1.3;
}else{
//pt_min = 0.4, pt_max = 365;
pt_min = 0.4, pt_max = 220;
ymin = 2E-18, ymax = 5E+1;
//ymin_r = 0.24, ymax_r = 1.56;
ymin_r = 0.44, ymax_r = 1.56;
}
setSysErrorFunction();
// file loading -----------------------------------------
file1 = (TFile*) loadFile(file1,
"./../rootOutput_postApp_3rd/feb072011/MergedAll_MBproc0710_Jet15Uproc0111_Jet50Uproc0111_eta2.4_feb07011_V2.root");
//"./../rootOutput_postApp_3rd/feb072011/MergedAll_MBproc0710_Jet15Uproc0111_Jet50Uproc0111_eta1.0_feb07011_V2.root");
//"./../rootOutput_postApp_3rd/feb072011/MergedAll_MBproc0710_Jet15Uproc0111_Jet50Uproc0111_eta2.4_feb07011_V2.root");
//"./../rootOutput_postApp_3rd/feb072011/MergedAll_MBproc0710_Jet15Uproc0111_Jet50Uproc0111_eta1.0_feb07011.root");
//"./../rootOutput_postApp_2nd/jan122011/MergedAll_MBproc0710_Jet15Uproc0111_Jet50Uproc0111_eta2.4_jan122011.root");
//"./../rootOutput_postApp_2nd/dec222010/MergedAll_MBproc0710_Jet15Uproc0710_Jet50Uproc1022_eta2.4_dec222010_scen_v2.root");
//"./../rootOutput_postApp_2nd/dec222010/MergedAll_MBproc0710_Jet15Uproc0710_Jet50Uproc1022_eta2.4_dec222010.root");
//"./../rootOutput_postApp_2nd/dec222010/MergedAll_MBproc0710_Jet15Uproc0710_Jet50Uproc1022_eta1.0_dec222010.root");
//"./../rootOutput_postApp_2nd/dec172010/MergedAll_MBproc0710_Jet15Uproc0710_Jet50Uproc1022_eta1.0_dec172010.root");
//"../corrections/rootOutput_aug3/CORR_TEST1_MB-C10-PR9-M6RR-JTMskim_Merged_0to2.4_proc0710.root"); // this has a shorter name because it's merged files
//"./rootOutput/CORR_All_Merged_witScale.root");
hdndpt1 = (TH1D*) file1->Get("corrTypeTwo/hInvCorrTwoLev2R");
TH1D *hdndpt1_temp = 0;
hdndpt1_temp = (TH1D*) hdndpt1->Clone("hdndpt1_temp");
hdndpt1_temp = (TH1D*) resetSysStatError(hdndpt1_temp);
tasydndpt1 = (TGraphAsymmErrors*) TGraphAsyErrIt(hdndpt1_temp);
//tdndpt1 = (TGraphErrors*) TgraphIt(hdndpt1_temp, true);
//tdndpt1 = (TGraphErrors*) file1->Get("corrTypeTwo/tInvCorrTwoLev2R");
if(hdndpt1==0) {
Error("[HistLoading]","could not get histogram");
return;
}
// D6T
file2 = (TFile*) loadFile(file2,
"./../rootOutput_postApp_2nd/dec222010_gen/GEN_TrkHistGEN_7TeV_dec14_qcdNSD_D6T_Pt0to20_eta_0.0to2.4_jet_0.0to2400.0_preTrackAna_varBin1_rebin4.root");
//"./../rootOutput_postApp_2nd/dec222010_gen/GEN_TrkHistGEN_7TeV_dec14_qcdNSD_D6T_Pt0to20_eta_0.0to1.0_jet_0.0to2400.0_preTrackAna_varBin1_rebin4.root");
//"../corrections//rootOutput/GEN_TrkHistGEN_july28_d6tv2_qcdPt0to20V2_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july19_qcdPt0to15_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
hdndpt2 = (TH1D*) file2->Get("corrTypeTwoGEN/hInvCorrR");
tdndpt2 = (TGraphErrors*) file2->Get("corrTypeTwoGEN/tInvCorr");
if(hdndpt2==0 || tdndpt2==0) {
Error("[HistLoading]","could not get histogram");
return;
}
// P0
file3 = (TFile*) loadFile(file3,
"./../rootOutput_postApp_2nd/dec222010_gen/GEN_TrkHistGEN_7TeV_dec14_qcdNSD_P0_Pt0to20_eta_0.0to2.4_jet_0.0to2400.0_preTrackAna_varBin1_rebin4.root");
//"./../rootOutput_postApp_2nd/dec222010_gen/GEN_TrkHistGEN_7TeV_dec14_qcdNSD_P0_Pt0to20_eta_0.0to1.0_jet_0.0to2400.0_preTrackAna_varBin1_rebin4.root");
//"../corrections/rootOutput/GEN_TrkHistGEN_july28_p0_qcdPt0to20V2_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july19_p0_qcdPt0to15_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july19_p0v2_qcdPt0to15V2_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july19_p0v3_qcdPt0to15_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
hdndpt3 = (TH1D*) file3->Get("corrTypeTwoGEN/hInvCorrR");
tdndpt3 = (TGraphErrors*) file3->Get("corrTypeTwoGEN/tInvCorr");
if(hdndpt3==0 || tdndpt3==0) {
Error("[HistLoading]","could not get histogram");
return;
}
// PYTHIA 8
file4 = (TFile*) loadFile(file4,
"./../rootOutput_postApp_2nd/dec222010_gen/GEN_TrkHistGEN_7TeV_dec14_qcdNSD_PY8_Pt0to30_eta_0.0to2.4_jet_0.0to2400.0_preTrackAna_varBin1_rebin4.root");
//"../corrections/rootOutput/GEN_TrkHistGEN_july28_py8v3_qcdPt0to20V3_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july28_proq20_qcdPt0to20V2_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july19_atlas_qcdPt0to15_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
hdndpt4 = (TH1D*) file4->Get("corrTypeTwoGEN/hInvCorrR");
tdndpt4 = (TGraphErrors*) file4->Get("corrTypeTwoGEN/tInvCorr");
if(hdndpt4==0 || tdndpt4==0) {
Error("[HistLoading]","could not get histogram");
return;
}
// X1 -> ProQ 20
file5 = (TFile*) loadFile(file5,
"./../rootOutput_postApp_2nd/dec222010_gen/GEN_TrkHistGEN_7TeV_dec14_qcdNSD_ProQ20_Pt0to20_eta_0.0to2.4_jet_0.0to2400.0_preTrackAna_varBin1_rebin4.root");
//"../corrections/rootOutput/GEN_TrkHistGEN_july28_proq20_qcdPt0to20V2_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july28_py8v3_qcdPt0to20V3_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
//"./rootOutput/GEN_TrkHistGEN_july19_x1_qcdPt0to15_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
hdndpt5 = (TH1D*) file5->Get("corrTypeTwoGEN/hInvCorrR");
tdndpt5 = (TGraphErrors*) file5->Get("corrTypeTwoGEN/tInvCorr");
if(hdndpt5==0 || tdndpt5==0) {
Error("[HistLoading]","could not get histogram");
return;
}
// Atlas
if(atlas){
file6 = (TFile*) loadFile(file6,
"../corrections/rootOutput/GEN_TrkHistGEN_july28_atlas_qcdPt0to20_eta_0.0to2.4_jet_0.0to2000.0_preTrackAna_varBin1_rebin4.root");
hdndpt6 = (TH1D*) file6->Get("corrTypeTwoGEN/hInvCorrR");
tdndpt6 = (TGraphErrors*) file6->Get("corrTypeTwoGEN/tInvCorr");
if(hdndpt6==0 || tdndpt6==0) {
Error("[HistLoading]","could not get histogram");
return;
}
}
// Normalization
//float N_evt = 1.13643e+07;
//hdndpt1->Scale(sigma_nsd);
//--------------------------- CMS measurement ---------------------------
cms_7000GeV = (TGraphErrors*) CMS_7TEV(1);
//---------------------------- fits --------------------------------------
TF1 *f2 = new TF1("fitTsallis","[0]*(1+(sqrt(0.1396**2+x**2)-0.1396)/([1]*[2]))**(-[2])",0.5,200);
f2->SetParameters(1.82045e+01,1.53540e-01,7.06942e+00); // fit the fit (used in paper)
f2->SetLineColor(2), f2->SetLineWidth(2);
// Basic canvas and dummy histogram
TCanvas *call = new TCanvas("call","call",510,670);
setPad(call);
// ----- pad 1
pp1->cd();
pp1->SetLogy();
pp1->SetLogx();
dum1 = GetDummyHist(pt_min,pt_max,ymin,ymax,"p_{T} [GeV/c]","Ed^{3}N/dp^{3} [GeV^{-2}c^{3}]",false);
dum1->Draw("");
dum1->GetXaxis()->SetNdivisions(908);
dum1->GetYaxis()->SetNdivisions(407);
int fillsty = 3001;
th1Style2(tdndpt2,2,20,1.0,2,2,1,2,true);
th1Style2(tdndpt3,4,20,1.0,4,2,2,2,false);
//th1Style2(tdndpt4,6,20,1.0,6,2,9,2,false);
th1Style2(tdndpt4,6,20,1.0,6,2,4,2,false);
//th1Style2(tdndpt5,19,20,1.0,19,2,3,3,false);
th1Style2(tdndpt5,20,20,1.0,20,2,9,3,false);
if(atlas) th1Style2(tdndpt6,20,20,1.0,20,2,4,3,false);
//th1Style1(hdndpt2,1,20,1.0,1,1.5,1,1);
tasydndpt1->SetMarkerSize(1.0); tasydndpt1->SetFillStyle(fillsty), tasydndpt1->SetFillColor(kGray);
tasydndpt1->Draw("2");
th1Style1(hdndpt1,1,20,1.0,1,1.5,1,1); // draw the measurement on top of everything
//TLegend *leg1 = new TLegend(0.37,0.47,0.74,0.85);
TLegend *leg1 = new TLegend(0.22,0.09,0.59,0.47);
leg1->SetMargin(0.3), leg1->SetBorderSize(0);
leg1->SetFillColor(kWhite), leg1->SetFillStyle(0);
leg1->SetTextSize(0.039), leg1->SetMargin(0.3);
leg1->SetHeader(" #sqrt{s} = 7 TeV, |#eta|<2.4");
leg1->AddEntry(tasydndpt1,"CMS #scale[0.7]{#int}Ldt = 2.96 pb^{-1}","plf");
leg1->AddEntry(tdndpt2,"PYTHIA D6T","l");
leg1->AddEntry(tdndpt3,"PYTHIA Perugia0","l");
leg1->AddEntry(tdndpt5,"PYTHIA ProQ20","l");
leg1->AddEntry(tdndpt4,"PYTHIA 8","l");
if(atlas) leg1->AddEntry(tdndpt6,"PYTHIA Atlas","l");
leg1->Draw();
//putCMSPrel();
//putIntLum();
putB();
// ----- pad 2
pp1_1->cd();
pp1_1->SetLogx();
dum1_1 = GetDummyHist(pt_min,pt_max,ymin_r,ymax_r,"p_{T} [GeV/c]","DATA / MC",false);
setLowerPad(dum1_1);
dum1_1->Draw("");
TLine *line = new TLine(0,1,pt_max,1);
line->SetLineStyle(2);
line->SetLineWidth(2);
line->Draw();
hdndpt1_2_cln = (TH1D*) hdndpt1->Clone("hdndpt1_2_cln");
hdndpt1_3_cln = (TH1D*) hdndpt1->Clone("hdndpt1_3_cln");
hdndpt1_4_cln = (TH1D*) hdndpt1->Clone("hdndpt1_4_cln");
hdndpt1_5_cln = (TH1D*) hdndpt1->Clone("hdndpt1_5_cln");
hdndpt2_cln = (TH1D*) hdndpt2->Clone("hdndpt2_cln");
hdndpt3_cln = (TH1D*) hdndpt3->Clone("hdndpt3_cln");
hdndpt4_cln = (TH1D*) hdndpt4->Clone("hdndpt4_cln");
hdndpt5_cln = (TH1D*) hdndpt5->Clone("hdndpt5_cln");
hdndpt1_2_cln->Divide(hdndpt2_cln);
hdndpt1_3_cln->Divide(hdndpt3_cln);
hdndpt1_4_cln->Divide(hdndpt4_cln);
hdndpt1_5_cln->Divide(hdndpt5_cln);
if(atlas){
hdndpt1_6_cln = (TH1D*) hdndpt1->Clone("hdndpt1_6_cln");
hdndpt6_cln = (TH1D*) hdndpt6->Clone("hdndpt6_cln");
hdndpt1_6_cln->Divide(hdndpt6_cln);
}
// Draw sys + stat combined error
createSysStatErrorHist(hdndpt1);
tsyserr = (TGraphErrors*) TgraphIt(hsyserr);
tsyserr->SetFillColor(kGray);
tsyserr->SetFillStyle(fillsty);
tsyserr->Draw("3");
//hsyserr->Draw("pzsame");
//th1Style1(cms_7000GeV_div,13,30,1.0,13,1,1,1);
th1Style1(hdndpt1_2_cln,2,20,1.0,2,2.0,1,3);
th1Style1(hdndpt1_3_cln,4,20,1.0,4,2.0,2,3);
//th1Style1(hdndpt1_4_cln,6,20,1.0,6,2.0,9,3);
//th1Style1(hdndpt1_5_cln,19,20,1.0,19,2.0,3,3);
th1Style1(hdndpt1_4_cln,6,20,1.0,6,2.0,4,3);
th1Style1(hdndpt1_5_cln,20,20,1.0,20,2.0,9,3);
if(atlas) th1Style1(hdndpt1_6_cln,20,20,1.0,20,2.0,5,3);
// print Xsection
PrintXsection(hdndpt1,1.);
//PrintXsection(hdndpt1,1.);
// save canvas
if(save){
if(atlas) printCanvases(call, "spectra_inv_cms_vs_pythia_all_v3",0);
else printCanvases(call, "spectra_inv_cms_vs_pythia_v3",0);
}
}
