Skip to content

Commit

Permalink
separate 1 and 2 jets in et/mt
Browse files Browse the repository at this point in the history
  • Loading branch information
ccaillol committed Nov 25, 2019
1 parent be719b0 commit 0c0ec41
Show file tree
Hide file tree
Showing 17 changed files with 3,313 additions and 246 deletions.
4 changes: 3 additions & 1 deletion FFcode/bin/OSSScorrection_et.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -376,6 +376,7 @@ int main(int argc, char** argv) {
TFile frawff("uncorrected_fakefactors_et.root");
TF1* ff_qcd_0jet=(TF1*) frawff.Get("rawFF_et_qcd_0jetSSloose");
TF1* ff_qcd_1jet=(TF1*) frawff.Get("rawFF_et_qcd_1jetSSloose");
TF1* ff_qcd_2jet=(TF1*) frawff.Get("rawFF_et_qcd_2jetSSloose");

TFile fmvisclosure ("FF_corrections_1.root");
TF1* mvisclosure_qcdloose=(TF1*) fmvisclosure.Get("closure_mvis_et_qcdloose");
Expand Down Expand Up @@ -634,7 +635,8 @@ int main(int argc, char** argv) {
//if (myvar>300) myvar=299;

float ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_0jet)*get_mvis_closure((myele+mytau).M(),mvisclosure_qcdloose);
if (njets>0) ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_1jet)*get_mvis_closure((myele+mytau).M(),mvisclosure_qcdloose);
if (njets==1) ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_1jet)*get_mvis_closure((myele+mytau).M(),mvisclosure_qcdloose);
else if (njets>1) ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_2jet)*get_mvis_closure((myele+mytau).M(),mvisclosure_qcdloose);

if (!is_includedInEmbedded){

Expand Down
4 changes: 3 additions & 1 deletion FFcode/bin/OSSScorrection_mt.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -415,6 +415,7 @@ int main(int argc, char** argv) {
TFile frawff("uncorrected_fakefactors_mt.root");
TF1* ff_qcd_0jet=(TF1*) frawff.Get("rawFF_mt_qcd_0jetSSloose");
TF1* ff_qcd_1jet=(TF1*) frawff.Get("rawFF_mt_qcd_1jetSSloose");
TF1* ff_qcd_2jet=(TF1*) frawff.Get("rawFF_mt_qcd_2jetSSloose");

TFile fmvisclosure ("FF_corrections_1.root");
TF1* mvisclosure_qcdloose=(TF1*) fmvisclosure.Get("closure_mvis_mt_qcdloose");
Expand Down Expand Up @@ -676,7 +677,8 @@ int main(int argc, char** argv) {
//if (myvar>300) myvar=299;

float ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_0jet)*get_mvis_closure((mymu+mytau).M(),mvisclosure_qcdloose);
if (njets>0) ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_1jet)*get_mvis_closure((mymu+mytau).M(),mvisclosure_qcdloose);
if (njets==1) ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_1jet)*get_mvis_closure((mymu+mytau).M(),mvisclosure_qcdloose);
if (njets>1) ff_qcd=get_raw_FF(mytau.Pt(),ff_qcd_2jet)*get_mvis_closure((mymu+mytau).M(),mvisclosure_qcdloose);

if (!is_includedInEmbedded){

Expand Down
144 changes: 132 additions & 12 deletions FFcode/bin/RawFF_et.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -330,23 +330,31 @@ int main(int argc, char** argv) {

TH1F* h0LT_qcd_iso = new TH1F ("h0LT_qcd_iso","h0LT_qcd_iso",binnum_mtt0,bins_mtt0); h0LT_qcd_iso->Sumw2();
TH1F* h1LT_qcd_iso = new TH1F ("h1LT_qcd_iso","h1LT_qcd_iso",binnum_mtt0,bins_mtt0); h1LT_qcd_iso->Sumw2();
TH1F* h2LT_qcd_iso = new TH1F ("h2LT_qcd_iso","h2LT_qcd_iso",binnum_mtt0,bins_mtt0); h2LT_qcd_iso->Sumw2();
TH1F* h0LT_qcd_anti = new TH1F ("h0LT_qcd_anti","h0LT_qcd_anti",binnum_mtt0,bins_mtt0); h0LT_qcd_anti->Sumw2();
TH1F* h1LT_qcd_anti = new TH1F ("h1LT_qcd_anti","h1LT_qcd_anti",binnum_mtt0,bins_mtt0); h1LT_qcd_anti->Sumw2();
TH1F* h2LT_qcd_anti = new TH1F ("h2LT_qcd_anti","h2LT_qcd_anti",binnum_mtt0,bins_mtt0); h2LT_qcd_anti->Sumw2();
TH1F* h0LT_w_iso = new TH1F ("h0LT_w_iso","h0LT_w_iso",binnum_mtt0,bins_mtt0); h0LT_w_iso->Sumw2();
TH1F* h1LT_w_iso = new TH1F ("h1LT_w_iso","h1LT_w_iso",binnum_mtt0,bins_mtt0); h1LT_w_iso->Sumw2();
TH1F* h2LT_w_iso = new TH1F ("h2LT_w_iso","h2LT_w_iso",binnum_mtt0,bins_mtt0); h2LT_w_iso->Sumw2();
TH1F* h0LT_w_anti = new TH1F ("h0LT_w_anti","h0LT_w_anti",binnum_mtt0,bins_mtt0); h0LT_w_anti->Sumw2();
TH1F* h1LT_w_anti = new TH1F ("h1LT_w_anti","h1LT_w_anti",binnum_mtt0,bins_mtt0); h1LT_w_anti->Sumw2();
TH1F* h2LT_w_anti = new TH1F ("h2LT_w_anti","h2LT_w_anti",binnum_mtt0,bins_mtt0); h2LT_w_anti->Sumw2();
TH1F* h0LT_tt_iso = new TH1F ("h0LT_tt_iso","h0LT_tt_iso",binnum_mtt0,bins_mtt0); h0LT_tt_iso->Sumw2();
TH1F* h0LT_tt_anti = new TH1F ("h0LT_tt_anti","h0LT_tt_anti",binnum_mtt0,bins_mtt0); h0LT_tt_anti->Sumw2();

TH1F* h0J_qcd_iso = new TH1F ("h0J_qcd_iso","h0J_qcd_iso",binnum_mtt0,bins_mtt0); h0J_qcd_iso->Sumw2();
TH1F* h1J_qcd_iso = new TH1F ("h1J_qcd_iso","h1J_qcd_iso",binnum_mtt0,bins_mtt0); h1J_qcd_iso->Sumw2();
TH1F* h2J_qcd_iso = new TH1F ("h2J_qcd_iso","h2J_qcd_iso",binnum_mtt0,bins_mtt0); h2J_qcd_iso->Sumw2();
TH1F* h0J_qcd_anti = new TH1F ("h0J_qcd_anti","h0J_qcd_anti",binnum_mtt0,bins_mtt0); h0J_qcd_anti->Sumw2();
TH1F* h1J_qcd_anti = new TH1F ("h1J_qcd_anti","h1J_qcd_anti",binnum_mtt0,bins_mtt0); h1J_qcd_anti->Sumw2();
TH1F* h2J_qcd_anti = new TH1F ("h2J_qcd_anti","h2J_qcd_anti",binnum_mtt0,bins_mtt0); h2J_qcd_anti->Sumw2();
TH1F* h0J_w_iso = new TH1F ("h0J_w_iso","h0J_w_iso",binnum_mtt0,bins_mtt0); h0J_w_iso->Sumw2();
TH1F* h1J_w_iso = new TH1F ("h1J_w_iso","h1J_w_iso",binnum_mtt0,bins_mtt0); h1J_w_iso->Sumw2();
TH1F* h2J_w_iso = new TH1F ("h2J_w_iso","h2J_w_iso",binnum_mtt0,bins_mtt0); h2J_w_iso->Sumw2();
TH1F* h0J_w_anti = new TH1F ("h0J_w_anti","h0J_w_anti",binnum_mtt0,bins_mtt0); h0J_w_anti->Sumw2();
TH1F* h1J_w_anti = new TH1F ("h1J_w_anti","h1J_w_anti",binnum_mtt0,bins_mtt0); h1J_w_anti->Sumw2();
TH1F* h2J_w_anti = new TH1F ("h2J_w_anti","h2J_w_anti",binnum_mtt0,bins_mtt0); h2J_w_anti->Sumw2();
TH1F* h0J_tt_iso = new TH1F ("h0J_tt_iso","h0J_tt_iso",binnum_mtt0,bins_mtt0); h0J_tt_iso->Sumw2();
TH1F* h0J_tt_anti = new TH1F ("h0J_tt_anti","h0J_tt_anti",binnum_mtt0,bins_mtt0); h0J_tt_anti->Sumw2();

Expand All @@ -358,6 +366,10 @@ int main(int argc, char** argv) {
TH1F* h1SSlooseLT_qcd_anti = new TH1F ("h1SSlooseLT_qcd_anti","h1SSlooseLT_qcd_anti",binnum_mtt0,bins_mtt0); h1SSlooseLT_qcd_anti->Sumw2();
TH1F* h1SSlooseJ_qcd_iso = new TH1F ("h1SSlooseJ_qcd_iso","h1SSlooseJ_qcd_iso",binnum_mtt0,bins_mtt0); h1SSlooseJ_qcd_iso->Sumw2();
TH1F* h1SSlooseJ_qcd_anti = new TH1F ("h1SSlooseJ_qcd_anti","h1SSlooseJ_qcd_anti",binnum_mtt0,bins_mtt0); h1SSlooseJ_qcd_anti->Sumw2();
TH1F* h2SSlooseLT_qcd_iso = new TH1F ("h2SSlooseLT_qcd_iso","h2SSlooseLT_qcd_iso",binnum_mtt0,bins_mtt0); h2SSlooseLT_qcd_iso->Sumw2();
TH1F* h2SSlooseLT_qcd_anti = new TH1F ("h2SSlooseLT_qcd_anti","h2SSlooseLT_qcd_anti",binnum_mtt0,bins_mtt0); h2SSlooseLT_qcd_anti->Sumw2();
TH1F* h2SSlooseJ_qcd_iso = new TH1F ("h2SSlooseJ_qcd_iso","h2SSlooseJ_qcd_iso",binnum_mtt0,bins_mtt0); h2SSlooseJ_qcd_iso->Sumw2();
TH1F* h2SSlooseJ_qcd_anti = new TH1F ("h2SSlooseJ_qcd_anti","h2SSlooseJ_qcd_anti",binnum_mtt0,bins_mtt0); h2SSlooseJ_qcd_anti->Sumw2();

string datapath = string(std::getenv("CMSSW_BASE"))+"/src/ComputeFF2018/FFcode/data/";

Expand Down Expand Up @@ -676,28 +688,40 @@ int main(int argc, char** argv) {
h0LT_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets==0 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h0LT_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>0 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
if (signalRegion && njets==1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h1LT_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>0 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
if (antiisoRegion && njets==1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h1LT_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h2LT_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h2LT_qcd_anti->Fill(myvar,aweight*weight2);

if (signalRegion && njets==0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h0SSlooseLT_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets==0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h0SSlooseLT_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
if (signalRegion && njets==1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h1SSlooseLT_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
if (antiisoRegion && njets==1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h1SSlooseLT_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h2SSlooseLT_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h2SSlooseLT_qcd_anti->Fill(myvar,aweight*weight2);

if (signalRegion && njets==0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h0LT_w_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets==0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h0LT_w_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
if (signalRegion && njets==1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h1LT_w_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
if (antiisoRegion && njets==1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h1LT_w_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h2LT_w_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h2LT_w_anti->Fill(myvar,aweight*weight2);

if (signalRegion && njets>=0 && iso_1<0.15 && nbtag>0 && mt<50 && q_1*q_2<0)
h0LT_tt_iso->Fill(myvar,aweight*weight2);
Expand All @@ -709,28 +733,40 @@ int main(int argc, char** argv) {
h0J_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets==0 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h0J_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>0 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
if (signalRegion && njets==1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h1J_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>0 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
if (antiisoRegion && njets==1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h1J_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h2J_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>1 && iso_1>0.02 && iso_1<0.15 && nbtag==0 && mt<50 && q_1*q_2>0)
h2J_qcd_anti->Fill(myvar,aweight*weight2);

if (signalRegion && njets==0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h0SSlooseJ_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets==0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h0SSlooseJ_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
if (signalRegion && njets==1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h1SSlooseJ_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>0 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
if (antiisoRegion && njets==1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h1SSlooseJ_qcd_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h2SSlooseJ_qcd_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>1 && iso_1>0.15 && iso_1<0.25 && nbtag==0 && mt<50 && q_1*q_2>0)
h2SSlooseJ_qcd_anti->Fill(myvar,aweight*weight2);

if (signalRegion && njets==0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h0J_w_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets==0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h0J_w_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
if (signalRegion && njets==1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h1J_w_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>0 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
if (antiisoRegion && njets==1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h1J_w_anti->Fill(myvar,aweight*weight2);
if (signalRegion && njets>1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h2J_w_iso->Fill(myvar,aweight*weight2);
if (antiisoRegion && njets>1 && iso_1<0.15 && nbtag==0 && mt>70 && q_1*q_2<0)
h2J_w_anti->Fill(myvar,aweight*weight2);

if (signalRegion && njets>=0 && iso_1<0.15 && nbtag>0 && mt<50 && q_1*q_2<0)
h0J_tt_iso->Fill(myvar,aweight*weight2);
Expand Down Expand Up @@ -804,6 +840,34 @@ int main(int argc, char** argv) {
h1J_qcd_anti->Write();
}

TDirectory *d2_qcd_iso =fout->mkdir("et_2jet_qcd_iso");
d2_qcd_iso->cd();
if (sample=="data_obs" or sample=="W"){
h2LT_qcd_iso->SetName(name.c_str());
h2LT_qcd_iso->Add(h2J_qcd_iso);
h2LT_qcd_iso->Write();
}
else{
h2LT_qcd_iso->SetName(name.c_str()+postfixLT);
h2LT_qcd_iso->Write();
h2J_qcd_iso->SetName(name.c_str()+postfixJ);
h2J_qcd_iso->Write();
}

TDirectory *d2_qcd_anti =fout->mkdir("et_2jet_qcd_anti");
d2_qcd_anti->cd();
if (sample=="data_obs" or sample=="W"){
h2LT_qcd_anti->SetName(name.c_str());
h2LT_qcd_anti->Add(h2J_qcd_anti);
h2LT_qcd_anti->Write();
}
else{
h2LT_qcd_anti->SetName(name.c_str()+postfixLT);
h2LT_qcd_anti->Write();
h2J_qcd_anti->SetName(name.c_str()+postfixJ);
h2J_qcd_anti->Write();
}

TDirectory *d0SSloose_qcd_iso =fout->mkdir("et_0SSloose_qcd_iso");
d0SSloose_qcd_iso->cd();
if (sample=="data_obs" or sample=="W"){
Expand Down Expand Up @@ -860,6 +924,34 @@ int main(int argc, char** argv) {
h1SSlooseJ_qcd_anti->Write();
}

TDirectory *d2SSloose_qcd_iso =fout->mkdir("et_2SSloose_qcd_iso");
d2SSloose_qcd_iso->cd();
if (sample=="data_obs" or sample=="W"){
h2SSlooseLT_qcd_iso->SetName(name.c_str());
h2SSlooseLT_qcd_iso->Add(h2SSlooseJ_qcd_iso);
h2SSlooseLT_qcd_iso->Write();
}
else{
h2SSlooseLT_qcd_iso->SetName(name.c_str()+postfixLT);
h2SSlooseLT_qcd_iso->Write();
h2SSlooseJ_qcd_iso->SetName(name.c_str()+postfixJ);
h2SSlooseJ_qcd_iso->Write();
}

TDirectory *d2SSloose_qcd_anti =fout->mkdir("et_2SSloose_qcd_anti");
d2SSloose_qcd_anti->cd();
if (sample=="data_obs" or sample=="W"){
h2SSlooseLT_qcd_anti->SetName(name.c_str());
h2SSlooseLT_qcd_anti->Add(h2SSlooseJ_qcd_anti);
h2SSlooseLT_qcd_anti->Write();
}
else{
h2SSlooseLT_qcd_anti->SetName(name.c_str()+postfixLT);
h2SSlooseLT_qcd_anti->Write();
h2SSlooseJ_qcd_anti->SetName(name.c_str()+postfixJ);
h2SSlooseJ_qcd_anti->Write();
}

TDirectory *d0_w_iso =fout->mkdir("et_0jet_w_iso");
d0_w_iso->cd();
if (sample=="data_obs" or sample=="W"){
Expand Down Expand Up @@ -916,6 +1008,34 @@ int main(int argc, char** argv) {
h1J_w_anti->Write();
}

TDirectory *d2_w_iso =fout->mkdir("et_2jet_w_iso");
d2_w_iso->cd();
if (sample=="data_obs" or sample=="W"){
h2LT_w_iso->SetName(name.c_str());
h2LT_w_iso->Add(h2J_w_iso);
h2LT_w_iso->Write();
}
else{
h2LT_w_iso->SetName(name.c_str()+postfixLT);
h2LT_w_iso->Write();
h2J_w_iso->SetName(name.c_str()+postfixJ);
h2J_w_iso->Write();
}

TDirectory *d2_w_anti =fout->mkdir("et_2jet_w_anti");
d2_w_anti->cd();
if (sample=="data_obs" or sample=="W"){
h2LT_w_anti->SetName(name.c_str());
h2LT_w_anti->Add(h2J_w_anti);
h2LT_w_anti->Write();
}
else{
h2LT_w_anti->SetName(name.c_str()+postfixLT);
h2LT_w_anti->Write();
h2J_w_anti->SetName(name.c_str()+postfixJ);
h2J_w_anti->Write();
}

TDirectory *d0_tt_iso =fout->mkdir("et_0jet_tt_iso");
d0_tt_iso->cd();
if (sample=="data_obs" or sample=="W"){
Expand Down
Loading

0 comments on commit 0c0ec41

Please sign in to comment.