void doAll_zmet_looper(bool skipFWLite = true)
{
// Load various tools
gROOT->ProcessLine(Form(".x setup.C(%d)", skipFWLite));
// Load FWLite
gSystem->Load("Tools/MiniFWLite/libMiniFWLite.so");
// Load and compile the looping code
gSystem->CompileMacro("zmet_looper.C","++k", "libzmet_looper");
zmet_looper* looper = new zmet_looper();
//use OS/SS baseline selection as documented in:
//http://www.t2.ucsd.edu/tastwiki/bin/view/CMS/SusyStudies3x
looper->set_susybaseline(0);
//make baby ntuple
looper->set_createTree(1);
//use bitmask selection
looper->set_useBitMask(0);
// K-factors
// these have been k-factors NLO/LO before
// now using them as sample normalizations to NLO
// these two are taken from Ceballos's pdf.
// It looks like the top x-section is for mtop = 175 GeV
float kttdil = 1.; // 375pb, 127000 events processed
float kttotr = 1.; // 375pb, 127000 events processed
float kWW = 1.;
float kWZ = 1.;
float kZZ = 1.;
float kWjets = 1.; // 11850 pb, 980000 events processed
float kWcharm = 1.1;
float kZjets = 1.;
float kDYee = 1.; // 1230 pb, 970360 events processed
float kDYmm = 1.; // 1230 pb, 970360 events processed
float kDYtautau = 1.; // 1230 pb, 970360 events processed
float kppMuX = 1.; // xsec/nevents
float kEM = 1.;
float ktW = 1.; // the evtScale is all negative for some reason
float kVQQ = 1.;
float kLM0 = 1.;
float kLM1 = 1.;
float kLM2 = 1.;
float kLM3 = 1.;
float kLM4 = 1.;
float kLM5 = 1.;
float kLM6 = 1.;
float kLM7 = 1.;
float kLM8 = 1.;
float kLM9 = 1.;
float kLM10 = 1.;
float kLM11 = 1.;
float kLM12 = 1.;
float kLMscan = 1.;
float kML1 = 1.;
float kML2 = 1.;
float kML3 = 1.;
float kML4 = 1.;
float kML5 = 1.;
float kML6 = 1.;
float kML7 = 1.;
float kML8 = 1.;
// Prescales
int prettdil = 1;
int prettotr = 1;
int preWW = 1;
int preWZ = 1;
int preZZ = 1;
int preWjets = 1;
int preWcharm = 1;
int preZjets = 1;
int preDYee = 1;
int preDYmm = 1;
int preDYtautau = 1;
int preppMuX = 1;
int preEM = 1;
int pretW = 1;
int preVQQ = 1;
int preLM0 = 1;
int preLM1 = 1;
int preLM2 = 1;
int preLM3 = 1;
int preLM4 = 1;
int preLM5 = 1;
int preLM6 = 1;
int preLM7 = 1;
int preLM8 = 1;
int preLM9 = 1;
int preLM10 = 1;
int preLM11 = 1;
int preLM12 = 1;
int preML1 = 1;
int preML2 = 1;
int preML3 = 1;
int preML4 = 1;
int preML5 = 1;
int preML6 = 1;
int preML7 = 1;
int preML8 = 1;
int preLMscan = 1;
//Flags for files to run over
bool runttdil = 1;
bool runttotr = 1;
bool runWW = 1;
bool runWZ = 1;
bool runZZ = 1;
bool runWjets = 1;
bool runWcharm = 0;
bool runZjets = 1;
bool runDYee = 0;
bool runDYmm = 0;
bool runDYtautau = 0;
bool runppMuX = 0;
bool runEM = 0;
bool runtW = 1;
bool runVQQ = 0;
bool runLM0 = 1;
bool runLM1 = 1;
bool runLM2 = 1;
bool runLM3 = 1;
bool runLM4 = 1;
bool runLM5 = 1;
bool runLM6 = 1;
bool runLM7 = 1;
bool runLM8 = 1;
bool runLM9 = 1;
bool runLM10 = 1;
bool runLM11 = 1;
bool runLM12 = 1;
bool runML1 = 0;
bool runML2 = 0;
bool runML3 = 0;
bool runML4 = 0;
bool runML5 = 0;
bool runML6 = 0;
bool runML7 = 0;
bool runML8 = 0;
bool runLMscan = 0;
// //Flags for files to run over
// bool runttdil = 0;
// bool runttotr = 0;
// bool runWW = 0;
// bool runWZ = 0;
// bool runZZ = 0;
// bool runWjets = 0;
// bool runWcharm = 0;
// bool runZjets = 0;
// bool runDYee = 0;
// bool runDYmm = 0;
// bool runDYtautau = 0;
// bool runppMuX = 0;
// bool runEM = 0;
// bool runtW = 0;
// bool runVQQ = 0;
// bool runLM0 = 0;
// bool runLM1 = 0;
// bool runLM2 = 0;
// bool runLM3 = 0;
// bool runLM4 = 0;
// bool runLM5 = 0;
// bool runLM6 = 0;
// bool runLM7 = 0;
// bool runLM8 = 0;
// bool runLM9 = 0;
// bool runLM10 = 0;
// bool runLM11 = 0;
// bool runLM12 = 0;
// bool runML1 = 1;
// bool runML2 = 1;
// bool runML3 = 1;
// bool runML4 = 1;
// bool runML5 = 1;
// bool runML6 = 1;
// bool runML7 = 1;
// bool runML8 = 1;
// bool runLMscan = 0;
TChain* chZjets = new TChain("Events");
if(runZjets){
pickSkimIfExists(chZjets,
"/tas/cms2/ZJets-madgraph_Spring10-START3X_V26_S09-v1/V03-04-08/merged*root",
"Zjets");
}
TChain* chtopdil = new TChain("Events");
if (runttdil) {
pickSkimIfExists(chtopdil,
"/tas/cms2/TTbarJets-madgraph_Spring10-START3X_V26_S09-v1/V03-04-07/merged*.root",
"TTJets");
}
TChain* chtopotr = new TChain("Events");
if(runttotr){
pickSkimIfExists(chtopotr,
"/tas/cms2/TTbarJets-madgraph_Spring10-START3X_V26_S09-v1/V03-04-07/merged*.root",
"TTJets");
}
TChain* chww = new TChain("Events");
if(runWW){
pickSkimIfExists(chww,
"/tas/cms2/WW_Spring10-START3X_V26_S09-v1_DiLep/V03-04-08/merged*root",
"WW");
}
TChain* chWZ = new TChain("Events");
if(runWZ){
pickSkimIfExists(chWZ,
"/tas/cms2/WZ_Spring10-START3X_V26_S09-v1/V03-04-08/merged*root",
"WZ_incl"); // can try WZ_3l-Pythia
}
TChain* chZZ = new TChain("Events");
if(runZZ){
pickSkimIfExists(chZZ,
"/tas/cms2/ZZ_Spring10-START3X_V26_S09-v1_DiLep/V03-04-08/merged*.root",
"ZZ");
}
TChain* chWjets = new TChain("Events");
if(runWjets){
pickSkimIfExists(chWjets,
"/tas/cms2/WJets-madgraph_Spring10-START3X_V26_S09-v1_SingleLep/V03-04-08/merged*.root",
"WJets");
}
//SAMPLE NOT YET AVAILABLE
// TChain* chWcharm = new TChain("Events");
// if(runWcharm){
// pickSkimIfExists(chWcharm,
// "data/Wc-madgraph_Fall08_IDEAL_V11_redigi_v1/merged*.root",
// "Wc");
// }
TChain* chDYtautau = new TChain("Events");
if(runDYtautau){
pickSkimIfExists(chDYtautau,
"/tas/cms2/Ztautau_Spring10-START3X_V26_S09-v1/V03-04-08-01/merged*.root",
"DYtautau");
}
// the low-mass splice has no choice other than the skim
// pickSkimIfExists(chDYtautau, "data/Ztautau_M20_Summer08_IDEAL_V11_redigi_v1/merged*.root", "Ztautau_M20");
TChain* chDYee = new TChain("Events");
if(runDYee){
pickSkimIfExists(chDYee,
"/tas/cms2/Zee_Spring10-START3X_V26_S09-v1/V03-04-08-01/merged*.root",
"DYee");
}
// the low-mass splice has no choice other than the skim
// pickSkimIfExists(chDYee, "data/Zee_M20_Summer08_IDEAL_V11_redigi_v1/merged*.root", "Zee_M20");
if(runDYmm){
TChain* chDYmm = new TChain("Events");
pickSkimIfExists(chDYmm,
"/tas/cms2/Zmumu_Spring10-START3X_V26_S09-v1/V03-04-08-01/merged*.root",
"DYmm");
}
// the low-mass splice has no choice other than the skim
// pickSkimIfExists(chDYmm, "data/Zmumu_M20_Summer08_IDEAL_V11_redigi_v1/merged*.root", "Zmumu_M20");
// ppMuX
// Spring10 sample not yet available!!!
// TChain* chppMuX = new TChain("Events");
// if (runppMuX) {
// pickSkimIfExists(chppMuX,
// "data3x/InclusiveMu15_Summer09-MC_31X_V3_7TeV-v1_dilepfilt/merged*.root",
// "InclusiveMuPt15");
// // can try InclusiveMu5Pt50 .. figure out how to merge later
// }
// ppEM
// Spring10 sample not yet available!!!
// TChain* chEM = new TChain("Events");
// if (runEM) {
// pickSkimIfExists(chEM,
// "data3x/QCD_EMEnriched_Pt20to30_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
// "_skimSimple2020");
// pickSkimIfExists(chEM,
// "data3x/QCD_EMEnriched_Pt30to80_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
// "_skimSimple2020");
// pickSkimIfExists(chEM,
// "data3x/QCD_EMEnriched_Pt80to170_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
// "_skimSimple2020");
// pickSkimIfExists(chEM,
// "data3x/QCD_BCtoE_Pt20to30_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
// "_skimSimple2020");
// pickSkimIfExists(chEM,
// "data3x/QCD_BCtoE_Pt30to80_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
// "_skimSimple2020");
// pickSkimIfExists(chEM,
// "data3x/QCD_BCtoE_Pt80to170_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
// "_skimSimple2020");
// }
// tW
TChain* chtW = new TChain("Events");
if (runtW) {
// pickSkimIfExists(chtW,
// //TRANSFER!!!
// "SingleTop_sChannel");
// pickSkimIfExists(chtW,
// //TRANSFER!!!
// "SingleTop_tChannel");
pickSkimIfExists(chtW,
"/tas/cms2/SingleTop_tWChannel-madgraph_Spring10-START3X_V26_S09-v1/V03-04-07/merged*.root",
"SingleTop_tWChannel");
}
//SAMPLE NOT YET AVAILABLE
// // VQQ
// TChain* chVQQ = new TChain("Events");
// if (runVQQ) {
// pickSkimIfExists(chVQQ,
// "data/VQQ-madgraph_Summer08_IDEAL_V11_redigi_v2/merged*.root",
// "VQQ");
// }
// LM points currently being processed!!!!
// LM0
TChain *chLM0 = new TChain("Events");
if (runLM0) {
pickSkimIfExists(chLM0,
"data3x/LM0_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM0");
}
// LM1
TChain *chLM1 = new TChain("Events");
if (runLM1) {
pickSkimIfExists(chLM1,
"data3x/LM1_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM1");
}
// LM2
TChain *chLM2 = new TChain("Events");
if (runLM2) {
pickSkimIfExists(chLM2,
"data3x/LM2_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM2");
}
// LM3
TChain *chLM3 = new TChain("Events");
if (runLM3) {
pickSkimIfExists(chLM3,
"data3x/LM3_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM3");
}
// LM4
TChain *chLM4 = new TChain("Events");
if (runLM4) {
pickSkimIfExists(chLM4,
"data3x/LM4_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM4");
}
// LM5
TChain *chLM5 = new TChain("Events");
if (runLM5) {
pickSkimIfExists(chLM5,
"data3x/LM5_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM5");
}
// LM6
TChain *chLM6 = new TChain("Events");
if (runLM6) {
pickSkimIfExists(chLM6,
"data3x/LM6_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM6");
}
// LM7
TChain *chLM7 = new TChain("Events");
if (runLM7) {
pickSkimIfExists(chLM7,
"data3x/LM7_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM7");
}
// LM8
TChain *chLM8 = new TChain("Events");
if (runLM8) {
pickSkimIfExists(chLM8,
"data3x/LM8_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM8");
}
// LM9
TChain *chLM9 = new TChain("Events");
if (runLM9) {
pickSkimIfExists(chLM9,
"data3x/LM9_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM9");
}
// LM10
TChain *chLM10 = new TChain("Events");
if (runLM10) {
pickSkimIfExists(chLM10,
"data3x/LM10_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM10");
}
// LM11
TChain *chLM11 = new TChain("Events");
if (runLM11) {
pickSkimIfExists(chLM11,
"data3x/LM11_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM11");
}
// LM12
TChain *chLM12 = new TChain("Events");
if (runLM12) {
pickSkimIfExists(chLM12,
"data3x/LM12_Summer09-MC_31X_V3_7TeV-v1/V03-00-35/merged*.root",
"SUSY_LM12");
}
// ML1
TChain *chML1 = new TChain("Events");
if (runML1) {
pickSkimIfExists(chML1,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML01_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML1");
}
// ML2
TChain *chML2 = new TChain("Events");
if (runML2) {
pickSkimIfExists(chML2,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML02_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML2");
}
// ML3
TChain *chML3 = new TChain("Events");
if (runML3) {
pickSkimIfExists(chML3,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML03_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML3");
}
// ML4
TChain *chML4 = new TChain("Events");
if (runML4) {
pickSkimIfExists(chML4,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML04_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML4");
}
// ML5
TChain *chML5 = new TChain("Events");
if (runML5) {
pickSkimIfExists(chML5,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML05_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML5");
}
// ML6
TChain *chML6 = new TChain("Events");
if (runML6) {
pickSkimIfExists(chML6,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML06_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML6");
}
// ML7
TChain *chML7 = new TChain("Events");
if (runML7) {
pickSkimIfExists(chML7,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML07_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML7");
}
// ML8
TChain *chML8 = new TChain("Events");
if (runML8) {
pickSkimIfExists(chML8,
"/tas/cms2/PhysicsProcess_PYTHIA6_SUSY_GMSM_SC_ML08_7TeV_v0/V03-04-13-01-gmsb/merged*root",
"SUSY_ML8");
}
// LMscan
TChain *chLMscan = new TChain("Events");
if (runLMscan) {
pickSkimIfExists(chLMscan,
"data3x/TANB3_CMSW336FASTv3/V03-00-37/merged*.root",
"LMscan");
}
char jetTypeStrings[2][5] = {"JPT", "calo"};
char metTypeStrings[3][8] = {"tcmet", "muon", "muonjes"};
bool doFakeApp = false;
char* zvetoStrings[4] = {"", "_allzveto", "_nozveto","_selectz"};
// Process files one at a time, and color them as needed
for (int jetTypeIdx = 1; jetTypeIdx < 2; ++jetTypeIdx)
{
for (int metTypeIdx = 0; metTypeIdx < 1; ++metTypeIdx)
{
for (int zvetoIdx = 0; zvetoIdx < 1; ++zvetoIdx)
{
zmet_looper::JetTypeEnum jetType(jetTypeIdx);
zmet_looper::MetTypeEnum metType(metTypeIdx);
zmet_looper::ZVetoEnum zveto(zvetoIdx);
if (runZjets) {
cout << "Processing Zjets" << endl;
looper->ScanChain(chZjets,"Zjets", kZjets, preZjets, jetType, metType, zveto, doFakeApp);
cout << "Done processing Zjets" << endl;
hist::color("Zjets", kBlack);
}
if (runttdil) {
cout << "Processing ttbar dileptonic.. " << endl;
looper->ScanChain(chtopdil,"ttdil", kttdil, prettdil, jetType, metType, zveto, doFakeApp);
cout << "Done processing ttbar dileptonic.. " << endl;
hist::color("ttdil", kYellow);
}
if (runttotr) {
cout << "Processing ttbar no-dileptons.. " << endl;
looper->ScanChain(chtopotr,"ttotr", kttotr, prettotr, jetType, metType, zveto, doFakeApp);
cout << "Done processing ttbar no-dileptons.. " << endl;
hist::color("ttotr", 30);
}
if (runWW) {
cout << "Processing WW.." << endl;
looper->ScanChain(chww,"ww", kWW, preWW, jetType, metType, zveto, doFakeApp);
cout << "Done processing WW.." << endl;
hist::color("ww", kRed);
}
if (runWZ) {
cout << "Processing WZ.." << endl;
looper->ScanChain(chWZ,"wz", kWZ, preWZ, jetType, metType, zveto, doFakeApp);
cout << "Done processing WZ.." << endl;
hist::color("wz", kBlue);
}
if (runZZ) {
cout << "Processing ZZ.." << endl;
looper->ScanChain(chZZ,"zz", kZZ, preZZ, jetType, metType, zveto, doFakeApp);
cout << "Done processing ZZ.." << endl;
hist::color("zz", kGreen);
}
if (runWjets) {
cout << "Processing Wjets.." << endl;
looper->ScanChain(chWjets,"wjets", kWjets, preWjets, jetType, metType, zveto, doFakeApp);
cout << "Done processing Wjets.." << endl;
hist::color("wjets", 40);
}
// if (runWcharm) { //SAMPLE NOT YET AVAILABLE
// cout << "Processing Wcharm.." << endl;
// looper->ScanChain(chWcharm, "wcharm", kWcharm, preWcharm, jetType, metType, zveto, doFakeApp);
// cout << "Done processing Wcharm.." << endl;
// hist::color("wcharm", 50);
// }
if (runDYtautau) {
cout << "Processing DY->tautau" << endl;
looper->ScanChain(chDYtautau,"DYtautau", kDYtautau, preDYtautau, jetType, metType, zveto, doFakeApp);
cout << "Done processing DY->tautau" << endl;
hist::color("DYtautau", kBlack);
}
if (runDYee) {
cout << "Processing DY->ee" << endl;
looper->ScanChain(chDYee,"DYee", kDYee, preDYee, jetType, metType, zveto, doFakeApp);
cout << "Done rocessing DY->ee" << endl;
hist::color("DYee", kMagenta);
}
if (runDYmm) {
cout << "Processing DY->mm" << endl;
looper->ScanChain(chDYmm,"DYmm", kDYmm, preDYmm, jetType, metType, zveto, doFakeApp);
cout << "Done processing DY->mm" << endl;
hist::color("DYmm", kCyan);
}
// if (runppMuX) {
// cout << "Processing ppMuX" << endl;
// looper->ScanChain(chppMuX,"ppMuX", kppMuX, preppMuX, jetType, metType, zveto, doFakeApp);
// cout << "Done processing ppMuX" << endl;
// hist::color("ppMuX", 51);
// }
// if (runEM) {
// cout << "Processing EM" << endl;
// looper->ScanChain(chEM,"EM", kEM, preEM, jetType, metType, zveto, doFakeApp);
// cout << "Done processing EM" << endl;
// hist::color("EM", 49);
// }
if (runtW) {
cout << "Processing tW" << endl;
looper->ScanChain(chtW,"tW", ktW, pretW, jetType, metType, zveto, doFakeApp);
cout << "Done processing tW" << endl;
hist::color("tW", 63);
}
// if (runVQQ) { //SAMPLE NOT YET AVAILABLE
// cout << "Processing VQQ" << endl;
// looper->ScanChain(chVQQ,"VQQ", kVQQ, preVQQ, jetType, metType, zveto, doFakeApp);
// cout << "Done processing VQQ" << endl;
// hist::color("VQQ", 45);
// }
if (runLM0) {
cout << "Processing LM0" << endl;
looper->ScanChain(chLM0, "LM0", kLM0, preLM0, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM0" << endl;
hist::color("LM0", kOrange);
}
if (runLM1) {
cout << "Processing LM1" << endl;
looper->ScanChain(chLM1, "LM1", kLM1, preLM1, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM1" << endl;
hist::color("LM1", kOrange+1);
}
if (runLM2) {
cout << "Processing LM2" << endl;
looper->ScanChain(chLM2, "LM2", kLM2, preLM2, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM2" << endl;
hist::color("LM2", kOrange+2);
}
if (runLM3) {
cout << "Processing LM3" << endl;
looper->ScanChain(chLM3, "LM3", kLM3, preLM3, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM3" << endl;
hist::color("LM3", kOrange+3);
}
if (runLM4) {
cout << "Processing LM4" << endl;
looper->ScanChain(chLM4, "LM4", kLM4, preLM4, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM4" << endl;
hist::color("LM4", kOrange+4);
}
if (runLM5) {
cout << "Processing LM5" << endl;
looper->ScanChain(chLM5, "LM5", kLM5, preLM5, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM5" << endl;
hist::color("LM5", kOrange+5);
}
if (runLM6) {
cout << "Processing LM6" << endl;
looper->ScanChain(chLM6, "LM6", kLM6, preLM6, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM6" << endl;
hist::color("LM6", kOrange+6);
}
if (runLM7) {
cout << "Processing LM7" << endl;
looper->ScanChain(chLM7, "LM7", kLM7, preLM7, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM7" << endl;
hist::color("LM7", kOrange+7);
}
if (runLM8) {
cout << "Processing LM8" << endl;
looper->ScanChain(chLM8, "LM8", kLM8, preLM8, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM8" << endl;
hist::color("LM8", kOrange+8);
}
if (runLM9) {
cout << "Processing LM9" << endl;
looper->ScanChain(chLM9, "LM9", kLM9, preLM9, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM9" << endl;
hist::color("LM9", kOrange+9);
}
if (runLM10) {
cout << "Processing LM10" << endl;
looper->ScanChain(chLM10, "LM10", kLM10, preLM10, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM10" << endl;
hist::color("LM10", kOrange+10);
}
if (runLM11) {
cout << "Processing LM11" << endl;
looper->ScanChain(chLM11, "LM11", kLM11, preLM11, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM11" << endl;
hist::color("LM11", kOrange-7);
}
if (runLM12) {
cout << "Processing LM12" << endl;
looper->ScanChain(chLM12, "LM12", kLM12, preLM12, jetType, metType, zveto, doFakeApp);
cout << "Done processing LM12" << endl;
hist::color("LM12", kOrange-7);
}
if (runML1) {
cout << "Processing ML1" << endl;
looper->ScanChain(chML1, "ML1", kML1, preML1, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML1" << endl;
}
if (runML2) {
cout << "Processing ML2" << endl;
looper->ScanChain(chML2, "ML2", kML2, preML2, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML2" << endl;
}
if (runML3) {
cout << "Processing ML3" << endl;
looper->ScanChain(chML3, "ML3", kML3, preML3, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML3" << endl;
}
if (runML4) {
cout << "Processing ML4" << endl;
looper->ScanChain(chML4, "ML4", kML4, preML4, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML4" << endl;
}
if (runML5) {
cout << "Processing ML5" << endl;
looper->ScanChain(chML5, "ML5", kML5, preML5, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML5" << endl;
}
if (runML6) {
cout << "Processing ML6" << endl;
looper->ScanChain(chML6, "ML6", kML6, preML6, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML6" << endl;
}
if (runML7) {
cout << "Processing ML7" << endl;
looper->ScanChain(chML7, "ML7", kML7, preML7, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML7" << endl;
}
if (runML8) {
cout << "Processing ML8" << endl;
looper->ScanChain(chML8, "ML8", kML8, preML8, jetType, metType, zveto, doFakeApp);
cout << "Done processing ML8" << endl;
}
if (runLMscan) {
cout << "Processing LMscan" << endl;
looper->ScanChain(chLMscan, "LMscan", kLMscan, preLMscan, jetType, metType, zveto, doFakeApp);
cout << "Done processing LMscan" << endl;
hist::color("LMscan", kOrange-7);
}
// save all the histograms
if(doFakeApp) {
const char* outFile = Form("root/victory_baseline_metgt50_sumjetptgt200_cand01_%s_%s%s_FakeApp_3x.root",
jetTypeStrings[jetTypeIdx], metTypeStrings[metTypeIdx],zvetoStrings[zvetoIdx]);
}
else {
const char* outFile = Form("root/victory_baseline_metgt50_sumjetptgt200_cand01_%s_%s%s_ML_3x.root",
jetTypeStrings[jetTypeIdx], metTypeStrings[metTypeIdx],zvetoStrings[zvetoIdx]);
}
//const char* outFile = Form("victory_baseline_genmetgt50_nosumjetptcut_%s_%s_pleasework_varbins.root",
//jetTypeStrings[jetTypeIdx], metTypeStrings[metTypeIdx]);
TDirectory *rootdir = gDirectory->GetDirectory("Rint:");
rootdir->cd();
saveHist(outFile);
deleteHistos();
}//zvetoIdx
} // metTypeIdx
} // jetTypeIdx
gSystem->Exit(0);
}
void doMassPoint(float analysis, Option option, RunEra runEra, int mHiggs[20],
double DYBkgScaleFactorWWPreselection[3], double DYBkgScaleFactorWWPreselectionKappa[3],
double DYBkgScaleFactorHiggsSelection[3][20], double DYBkgScaleFactorHiggsSelectionKappa[3][20],
double DYBkgScaleFactorHiggsSelectionMVA[3][20], double DYBkgScaleFactorHiggsSelectionKappaMVA[3][20])
{
float lumi = 19467; // in unit of 1/pb
//
// set up the looper
//
SmurfDYLooper *looper = new SmurfDYLooper(analysis, option, runEra);
//looper->setGoodRunList("../runlists/runlist_1092.txt");
looper->setLumiScale(lumi, lumi);
//
// set up samples
//
int ana = (int)analysis;
SmurfSample *sample_data = new SmurfSample(option, DATA, kBlack, "Data", analysis);
SmurfSample *sample_dyll = new SmurfSample(option, ZLL, kBlue, "DYLL", analysis);
SmurfSample *sample_vv = new SmurfSample(option, VV, kGreen, "VV", analysis);
// examples of using the unskimed files
char *dataDir = "/smurf/jaehyeok/data/Run2012_Summer12_SmurfV9_53X/mitf-alljets_19p5ifb_new/dyskim/"; // TAS
//char *dataDir = "/nfs-7/userdata/jaehyeok/smurfntuples/mitf-alljets/dyskim/"; // UAF
sample_data->add(Form("%s/data.root", dataDir));
sample_dyll->add(Form("%s/dyll.root", dataDir));
sample_vv->add(Form("%s/wz.root", dataDir));
sample_vv->add(Form("%s/zz.root", dataDir));
//
// do the looping
//
looper->loop(sample_data);
looper->loop(sample_dyll);
looper->loop(sample_vv);
//
// save histograms
//
std::string outFile = Form("dyhistos_ww_%i.root", int(analysis));
if ( option == WW_OPT_SMURFXSECSEL ) {
outFile = Form("dyhistos_ww_%i_wwxsec.root", int(analysis));
}
if ( option == HWW_OPT_SMURFCUTSEL ) {
outFile = Form("dyhistos_ww_%i_cut.root", int(analysis));
}
/*
if ( option == HWW_OPT_SMURFMVASEL ) {
outFile = Form("dyhistos_ww_%i_mva.root", int(analysis));
std::cout << "outputFile Name = " << outFile.c_str() << "\n";
}
*/
saveHist(outFile.c_str());
deleteHistos();
std::string debugFileName = Form("dyest_mH%i_%.0fpb.txt", int(analysis), lumi);
if ( option == WW_OPT_SMURFXSECSEL)
debugFileName = Form("dyest_mH%i_%.0fpb_wwxsec.txt", int(analysis), lumi);
if ( option == HWW_OPT_SMURFCUTSEL)
debugFileName = Form("dyest_mH%i_%.0fpb_cut.txt", int(analysis), lumi);
/*
if ( option == HWW_OPT_SMURFMVASEL)
debugFileName = Form("dyest_mH%i_%.0fpb_mva.txt", int(analysis), lumi);
*/
FILE *debugtext = fopen(debugFileName.c_str(), "w");
//
// fill the R
//
std::string ratioFileName = Form("Routin_mH%i_%.0fpb.root", int(analysis), lumi);
if ( option == WW_OPT_SMURFXSECSEL) {
ratioFileName = Form("Routin_mH%i_%.0fpb_wwxsec.root", int(analysis), lumi);
}
if ( option == HWW_OPT_SMURFCUTSEL) {
ratioFileName = Form("Routin_mH%i_%.0fpb_cut.root", int(analysis), lumi);
}
/*
if ( option == HWW_OPT_SMURFMVASEL) {
ratioFileName = Form("Routin_mH%i_%.0fpb_mva.root", int(analysis), lumi);
}
*/
fillRoutin(outFile.c_str(), ratioFileName.c_str(), debugtext);
dyest(analysis, option, outFile.c_str(), ratioFileName.c_str(), debugtext, mHiggs,
DYBkgScaleFactorWWPreselection, DYBkgScaleFactorWWPreselectionKappa,
DYBkgScaleFactorHiggsSelection, DYBkgScaleFactorHiggsSelectionKappa,
DYBkgScaleFactorHiggsSelectionMVA, DYBkgScaleFactorHiggsSelectionKappaMVA );
fclose(debugtext);
//
// clean up
//
delete looper;
delete sample_data;
delete sample_dyll;
delete sample_vv;
}
void doAll(bool skipFWLite = true)
{
//---------------------------------------------------------------
// choose version, output will be written to output/[version]
//---------------------------------------------------------------
const char* version = "V00-01-06";
const char* jsonfile = "jsons/Cert_EPSFINAL_May10ReReco_v2_PromptReco_160404_167913_JSON_goodruns.txt";
const bool useMCSkims = true;
cout << "Version : " << version << endl;
cout << "json : " << jsonfile << endl;
//Load CORE stuff
gROOT->ProcessLine(".L ../CORE/CMS2.cc+");
gROOT->ProcessLine(".L ../CORE/utilities.cc+");
gROOT->ProcessLine(".L ../CORE/trackSelections.cc+");
gROOT->ProcessLine(".L ../CORE/eventSelections.cc+");
gROOT->ProcessLine(".L ../CORE/MITConversionUtilities.cc+");
gROOT->ProcessLine(".L ../CORE/muonSelections.cc+");
gROOT->ProcessLine(".L ../CORE/electronSelectionsParameters.cc+");
gROOT->ProcessLine(".L ../CORE/electronSelections.cc+");
gROOT->ProcessLine(".L ../CORE/metSelections.cc+");
gROOT->ProcessLine(".L ../CORE/SimpleFakeRate.cc+");
gROOT->ProcessLine(".L ../CORE/mcSelections.cc+");
gROOT->ProcessLine(".L ../CORE/MT2/MT2.cc+");
gROOT->ProcessLine(".L ../CORE/triggerUtils.cc+");
gROOT->ProcessLine(".L ../CORE/susySelections.cc+");
gROOT->ProcessLine(".L ../CORE/mcSUSYkfactor.cc+");
gROOT->ProcessLine(".L ../CORE/triggerSuperModel.cc+");
gROOT->ProcessLine(".L ../CORE/triggerUtils.cc+");
//gROOT->ProcessLine(".L ../CORE/jetSelections.cc+");
gROOT->ProcessLine(".L ../CORE/ttbarSelections.cc+");
// Load various tools
gROOT->ProcessLine(Form(".x setup.C(%d)", skipFWLite));
//gROOT->ProcessLine(".L ../CORE/topmass/ttdilepsolve.cpp+");
// Load FWLite
gSystem->Load("../Tools/MiniFWLite/libMiniFWLite.so");
// Load and compile the looping code
gSystem->CompileMacro("singleLeptonLooper.C","++k", "libsingleLeptonLooper");
singleLeptonLooper* looper = new singleLeptonLooper();
//set looper parameters
looper->set_susybaseline(0);
//make baby ntuple
looper->set_createTree(1);
//use bitmask selection
looper->set_useBitMask(0);
//set version
looper->set_version(version);
//set json
looper->set_json( jsonfile );
// k-factors
float kttall = 1.;
float kqcd = 1.;
float kWjets = 1.;
float kDYtot = 1.;
float ktW = 1.;
// prescales
int preqcd = 1;
int prettall = 1;
int preWjets = 1;
int preDYtot = 1;
int pretW = 1;
// flags for files to run over
bool rundata = 0;
bool runttall = 1;
bool runWjets = 0;
bool runQCD = 0;
bool runtW = 0;
bool runDYtot = 0;
bool runT2tt = 0;
bool runT2tt_few = 0;
bool runT2bw = 0;
if( useMCSkims ) cout << "Using MC skims" << endl;
else cout << "Using full MC samples" << endl;
//----------------------------------------
// QCD
//----------------------------------------
TChain* chQCD = new TChain("Events");
if(runQCD){
string skimdir = "/home/users/benhoob/filters/output/";
//SingleLeptonSkim ntuples
if( useMCSkims ){
pickSkimIfExists(chQCD, skimdir+"QCD_Pt-15to30_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/SingleLeptonSkim/merged*root");
pickSkimIfExists(chQCD, skimdir+"QCD_Pt-30to50_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/SingleLeptonSkim/merged*root");
pickSkimIfExists(chQCD, skimdir+"QCD_Pt-50to80_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/SingleLeptonSkim/merged*root");
pickSkimIfExists(chQCD, skimdir+"QCD_Pt-80to120_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/SingleLeptonSkim/merged*root");
pickSkimIfExists(chQCD, skimdir+"QCD_Pt-120to170_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/SingleLeptonSkim/merged*root");
pickSkimIfExists(chQCD, skimdir+"QCD_Pt-170to300_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/SingleLeptonSkim/merged*root");
}
//full cms2 ntuples
else{
pickSkimIfExists(chQCD, "/nfs-7/userdata/cms2/QCD_Pt-15to30_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/merged*root");
pickSkimIfExists(chQCD, "/nfs-7/userdata/cms2/QCD_Pt-30to50_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/merged*root");
pickSkimIfExists(chQCD, "/nfs-7/userdata/cms2/QCD_Pt-50to80_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/merged*root");
pickSkimIfExists(chQCD, "/nfs-7/userdata/cms2/QCD_Pt-80to120_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/merged*root");
pickSkimIfExists(chQCD, "/nfs-7/userdata/cms2/QCD_Pt-120to170_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/merged*root");
pickSkimIfExists(chQCD, "/nfs-7/userdata/cms2/QCD_Pt-170to300_TuneZ2_7TeV_pythia6_Summer11-PU_S4_START42_V11-v1/V04-02-31/merged*root");
}
}
//----------------------------------------
// ttbar
//----------------------------------------
TChain* chtopall = new TChain("Events");
if (runttall) {
pickSkimIfExists(chtopall,"/nfs-7/userdata/cms2/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29_singleLepton/merged*root");
//pickSkimIfExists(chtopall,"/nfs-7/userdata/cms2/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29_singleLepton/merged_ntuple.root");
}
//----------------------------------------
// W+jets
//----------------------------------------
TChain* chWjets = new TChain("Events");
if(runWjets){
//SingleLeptonSkim ntuples
if( useMCSkims ){
pickSkimIfExists(chWjets,"/hadoop/cms/store/group/snt/papers2011/Summer11MC/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/SingleLeptonAndTwoJets/merged_ntuple_*.root");
//pickSkimIfExists(chWjets,"/nfs-7/userdata/cms2/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/SingleLeptonAndJets/merged*root");//l+3j filtered
//pickSkimIfExists(chWjets,"/hadoop/cms/store/user/imacneill/Summer11MC/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/SingleLeptonAndJets/merged*root");
}
//full cms2 ntuples
else{
pickSkimIfExists(chWjets,"/nfs-7/userdata/cms2/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29_singleLepton/merged*root");
}
}
//----------------------------------------
// DY
//----------------------------------------
TChain* chDYtot = new TChain("Events");
if(runDYtot){
cout << "ERROR! NEED TO ADD SINGLE LEPTON SKIM OF DY SAMPLES!!!" << endl;
exit(0);
}
//----------------------------------------
// single top
//----------------------------------------
TChain* chtW = new TChain("Events");
if (runtW) {
pickSkimIfExists(chtW,"/nfs-7/userdata/benhoob/cms2/T_TuneZ2_s-channel_7TeV-powheg-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/merged*root");
pickSkimIfExists(chtW,"/nfs-7/userdata/benhoob/cms2/Tbar_TuneZ2_s-channel_7TeV-powheg-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/merged*root");
pickSkimIfExists(chtW,"/nfs-7/userdata/benhoob/cms2/T_TuneZ2_t-channel_7TeV-powheg-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/merged*root");
pickSkimIfExists(chtW,"/nfs-7/userdata/benhoob/cms2/Tbar_TuneZ2_t-channel_7TeV-powheg-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/merged*root");
pickSkimIfExists(chtW,"/nfs-7/userdata/benhoob/cms2/T_TuneZ2_tW-channel-DR_7TeV-powheg-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/merged*root");
pickSkimIfExists(chtW,"/nfs-7/userdata/benhoob/cms2/Tbar_TuneZ2_tW-channel-DR_7TeV-powheg-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29/merged*root");
}
//----------------------------------------
// T2tt (a few sample points)
//----------------------------------------
TChain *chT2tt_few = new TChain("Events");
if (runT2tt_few) {
// these are the T2tt files with m(stop),m(LSP) = 400,100 OR 500,100 OR 500,300
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_37.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_46.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_66.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_71.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_84.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_85.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_88.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_147.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_149.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_151.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_191.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_193.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_212.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_213.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_222.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_225.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_232.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_268.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_328.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_390.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_395.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_396.root");
// these are the T2tt files with m(stop),m(LSP) = 250,50 OR 300,50 OR 400,50 OR 400,200
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_1.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_5.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_6.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_16.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_62.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_63.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_70.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_91.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_105.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_149.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_150.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_160.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_171.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_243.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_246.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_264.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_274.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_312.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_315.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_325.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_326.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_333.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_335.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_349.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_351.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_369.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_388.root");
// pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_394.root");
//350/100 (from ATLAS paper)
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_122.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_274.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_286.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_326.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_347.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_364.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_391.root");
pickSkimIfExists(chT2tt_few,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged_ntuple_392.root");
}
//----------------------------------------
// T2tt (a few sample points)
//----------------------------------------
TChain *chT2tt = new TChain("Events");
if (runT2tt) {
pickSkimIfExists(chT2tt,"/nfs-7/userdata/cms2/SMS-T2tt_Mstop-225to1200_mLSP-50to1025_7TeV-Pythia6Z_Summer11-PU_START42_V11_FastSim-v1/V04-02-20-04/merged*root");
}
//----------------------------------------
// T2bw
//----------------------------------------
TChain *chT2bw = new TChain("Events");
if (runT2bw) {
// pickSkimIfExists(chT2bw,"/nfs-7a/userdata/cms2/SMS-T2bw_x-0p25to0p75_mStop-50to850_mLSP-50to800_7TeV-Pythia6Z_Summer11-PU_START42_V11_FSIM-v1/VB04-02-29_Fastsim/merged_ntuple_98.root");
// pickSkimIfExists(chT2bw,"/nfs-7a/userdata/cms2/SMS-T2bw_x-0p25to0p75_mStop-50to850_mLSP-50to800_7TeV-Pythia6Z_Summer11-PU_START42_V11_FSIM-v1/VB04-02-29_Fastsim/merged_ntuple_99.root");
// pickSkimIfExists(chT2bw,"/nfs-7a/userdata/cms2/SMS-T2bw_x-0p25to0p75_mStop-50to850_mLSP-50to800_7TeV-Pythia6Z_Summer11-PU_START42_V11_FSIM-v1/VB04-02-29_Fastsim/merged_ntuple_100.root");
// pickSkimIfExists(chT2bw,"/nfs-7a/userdata/cms2/SMS-T2bw_x-0p25to0p75_mStop-50to850_mLSP-50to800_7TeV-Pythia6Z_Summer11-PU_START42_V11_FSIM-v1/VB04-02-29_Fastsim/merged_ntuple_101.root");
pickSkimIfExists(chT2bw,"/nfs-7a/userdata/cms2/SMS-T2bw_x-0p25to0p75_mStop-50to850_mLSP-50to800_7TeV-Pythia6Z_Summer11-PU_START42_V11_FSIM-v1/VB04-02-29_Fastsim/merged*root");
}
//----------------------------------------
// DATA
//----------------------------------------
TChain* chdata = new TChain("Events");
if(rundata){
cout << "adding ElectronHad and MuHad or SingleMu data" << endl;
//---------------------------
// May10 rereco
//---------------------------
//pickSkimIfExists(chdata,"/hadoop/cms/store/user/yanjuntu/CMSSW_4_2_7_patch1_V04-02-33/MuHad_Run2011A-May10ReReco-v1_AOD/CMSSW_4_2_7_patch1_V04-02-33_merged/V04-02-33/merged*root");
//pickSkimIfExists(chdata,"/hadoop/cms/store/user/yanjuntu/CMSSW_4_2_7_patch1_V04-02-33/ElectronHad_Run2011A-May10ReReco-v1_AOD/CMSSW_4_2_7_patch1_V04-02-33_merged/V04-02-33/merged*root");
// pickSkimIfExists(chdata,"/nfs-7/userdata/cms2/ElectronHad_Run2011A-May10ReReco-v1_AOD/V04-02-33/SingleLeptonAndJets/merged*root");//l+3j filtered
// pickSkimIfExists(chdata,"/nfs-7/userdata/cms2/MuHad_Run2011A-May10ReReco-v1_AOD/V04-02-33/SingleLeptonAndJets/merged*root");//l+3j filtered
pickSkimIfExists(chdata,"/hadoop/cms/store/user/vimartin/SingleLeptonAndTwoJets/SingleMu_Run2011A-May10ReReco-v1_AOD/V04-02-33/SingleLeptonAndTwoJets/merged*root");//l+2j filtered
pickSkimIfExists(chdata,"/hadoop/cms/store/user/vimartin/SingleLeptonAndTwoJets/ElectronHad_Run2011A-May10ReReco-v1_AOD/V04-02-33/SingleLeptonAndTwoJets/merged*root");//l+2j filtered
//---------------------------
// prompt reco v4
//---------------------------
//pickSkimIfExists(chdata,"/hadoop/cms/store/user/yanjuntu/CMSSW_4_2_7_patch1_V04-02-33/ElectronHad_Run2011A-PromptReco-v4_AOD/CMSSW_4_2_7_patch1_V04-02-33_merged/V04-02-33/merged*root");
//pickSkimIfExists(chdata,"/hadoop/cms/store/user/yanjuntu/CMSSW_4_2_7_patch1_V04-02-33/MuHad_Run2011A-PromptReco-v4_AOD/CMSSW_4_2_7_patch1_V04-02-33_merged/V04-02-33/merged*root");
// pickSkimIfExists(chdata,"/nfs-7/userdata/cms2/ElectronHad_Run2011A-PromptReco-v4_AOD/V04-02-33/SingleLeptonAndJets/merged*root");//l+3j filtered
// pickSkimIfExists(chdata,"/nfs-7/userdata/cms2/MuHad_Run2011A-PromptReco-v4_AOD/V04-02-33/SingleLeptonAndJets/merged*root");//l+3j filtered
pickSkimIfExists(chdata,"/hadoop/cms/store/user/vimartin/SingleLeptonAndTwoJets/SingleMu_Run2011A-PromptReco-v4_AOD/V04-02-33/SingleLeptonAndTwoJets/merged*root");//l+2j filtered
pickSkimIfExists(chdata,"/hadoop/cms/store/user/vimartin/SingleLeptonAndTwoJets/ElectronHad_Run2011A-PromptReco-v4_AOD/V04-02-33/SingleLeptonAndTwoJets/merged*root");//l+2j filtered
}
//--------------------------------
//set luminosity to scale to
//--------------------------------
float lumi = 1.0;
bool calculateTCMET = false; //redo tcmet calculation on the fly
char* jetTypeStrings[3] = {"JPT", "calo","pfjet"};
char* metTypeStrings[4] = {"tcmet", "muon", "muonjes","pfmet"};
char* zvetoStrings[4] = {"", "_allzveto", "_nozveto","_selectz"};
char* frmodeStrings[2] = {"QCDType","WjetsType"}; //e_qcd = 0, e_wjets
bool doFakeApp = false;
singleLeptonLooper::TrigEnum trig;
for (int jetTypeIdx = 2; jetTypeIdx < 3; ++jetTypeIdx)
{
for (int metTypeIdx = 3; metTypeIdx < 4; ++metTypeIdx)
{
for (int zvetoIdx = 0; zvetoIdx < 1; ++zvetoIdx)
{
for (int frmodeIdx = 0; frmodeIdx < (2-(1*!doFakeApp)); ++frmodeIdx)
//for (int frmodeIdx = 0; frmodeIdx < 1; ++frmodeIdx)
{
singleLeptonLooper::JetTypeEnum jetType(jetTypeIdx);
singleLeptonLooper::MetTypeEnum metType(metTypeIdx);
singleLeptonLooper::ZVetoEnum zveto(zvetoIdx);
singleLeptonLooper::FREnum frmode(frmodeIdx);
if( doFakeApp ){
if( frmodeIdx == 0 ) cout << "Doing double fake estimate" << endl;
if( frmodeIdx == 1 ) cout << "Doing single fake estimate" << endl;
}
//--------------------------------------------------------------------
if (runttall) {
cout << "Processing ttbar all.. " << endl;
looper->ScanChain(chtopall,"ttall", kttall, prettall, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing ttbar all.. " << endl;
}
//--------------------------------------------------------------------
if (rundata) {
cout << "Processing data" << endl;
looper->ScanChain(chdata,"data", 1, 1, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing data" << endl;
}
//--------------------------------------------------------------------
if (runDYtot) {
cout << "Processing DY->all" << endl;
looper->ScanChain(chDYtot,"DYtot", kDYtot, preDYtot, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done rocessing DY->ee" << endl;
}
//--------------------------------------------------------------------
if (runQCD) {
cout << "Processing QCD.. " << endl;
looper->ScanChain(chQCD,"qcd", kqcd, preqcd, lumi, jetType, metType, zveto,frmode, doFakeApp, calculateTCMET);
cout << "Done processing QCD.. " << endl;
}
//--------------------------------------------------------------------
if (runWjets) {
cout << "Processing Wjets.." << endl;
looper->ScanChain(chWjets,"wjets", kWjets, preWjets, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing Wjets.." << endl;
}
//--------------------------------------------------------------------
if (runtW) {
cout << "Processing tW" << endl;
looper->ScanChain(chtW,"tW", ktW, pretW, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing tW" << endl;
}
//--------------------------------------------------------------------
if (runT2tt) {
cout << "Processing T2tt" << endl;
looper->ScanChain(chT2tt, "T2tt", 1, 1, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing T2tt" << endl;
}
//--------------------------------------------------------------------
if (runT2tt_few) {
cout << "Processing T2tt_few" << endl;
looper->ScanChain(chT2tt_few, "T2tt_few", 1, 1, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing T2tt_few" << endl;
}
//--------------------------------------------------------------------
if (runT2bw) {
cout << "Processing T2bw all.. " << endl;
looper->ScanChain(chT2bw,"T2bw", 1, 1, lumi, jetType, metType, zveto, frmode, doFakeApp, calculateTCMET);
cout << "Done processing T2bw all.. " << endl;
}
// save all the histograms
const char* outFile;
if(doFakeApp) {
outFile = Form("../output/%s/singleLepton_%s_%s%s_%s_FakeApp.root", version,
jetTypeStrings[jetTypeIdx], metTypeStrings[metTypeIdx],zvetoStrings[zvetoIdx],frmodeStrings[frmode]);
}
else {
outFile = Form("../output/%s/singleLepton_%s_%s%s.root", version,
jetTypeStrings[jetTypeIdx], metTypeStrings[metTypeIdx],zvetoStrings[zvetoIdx]);
}
//const char* outFile = Form("victory_baseline_genmetgt50_nosumjetptcut_%s_%s_pleasework_varbins.root",
//jetTypeStrings[jetTypeIdx], metTypeStrings[metTypeIdx]);
TDirectory *rootdir = gDirectory->GetDirectory("Rint:");
rootdir->cd();
saveHist(outFile);
deleteHistos();
} // frmodeIdx
}//zvetoIdx
} // metTypeIdx
} // jetTypeIdx
//}
gSystem->Exit(0);
}
void doAll(bool skipFWLite = true)
{
//---------------------------------------------------------------
// choose version, output will be written to output/[version]
//---------------------------------------------------------------
const char* version = "V00-00-10";
const char* jsonfile = "jsons/Cert_190456-202305_8TeV_PromptReco_Collisions12_JSON_goodruns.txt";
const bool useMCSkims = true;
cout << "Version : " << version << endl;
cout << "json : " << jsonfile << endl;
// Load various tools
gROOT->ProcessLine(Form(".x setup.C(%d)", skipFWLite));
// Load FWLite
gSystem->Load("../Tools/MiniFWLite/libMiniFWLite.so");
// Load and compile the looping code
gSystem->CompileMacro("looper.C","++k", "liblooper");
looper* looper = new looper();
//set looper parameters
looper->set_susybaseline(0);
//make baby ntuple
looper->set_createTree(1);
//use bitmask selection
looper->set_useBitMask(0);
//set version
looper->set_version(version);
//set json
looper->set_json( jsonfile );
//----------------------------------------
// flags for files to run over
//----------------------------------------
bool runData = 0;
bool runElData = 0;
bool runMuData = 0;
bool runSingleElData = 0;
bool runSingleMuData = 0;
bool runElHadData = 0;
bool runDoubleElData = 0;
bool runDoubleMuData = 0;
bool runMuHadData = 0;
bool runPhotonData = 0;
bool runData2012A = 0;
bool runData2012B = 1;
bool runData2012C = 1;
//----------------------------------------------------------------------------------------------------------
// muon trigger skim: 2012A
//----------------------------------------------------------------------------------------------------------
TChain* chData2012A = new TChain("Events");
if( runData2012A ){
cout << "Adding all 2012A data muon trigger skim" << endl;
pickSkimIfExists(chData2012A,"/home/users/benhoob/CondorSkim/MuonTrigger/MuHad_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged_ntuple_999999_12_2_skim.root");
//pickSkimIfExists(chData,"/home/users/benhoob/CondorSkim/MuonTrigger/MuHad_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------------------------
// muon trigger skim: 2012B
//----------------------------------------------------------------------------------------------------------
TChain* chData2012B = new TChain("Events");
if( runData2012B ){
cout << "Adding all 2012B data muon trigger skim" << endl;
//pickSkimIfExists(chData2012B,"/home/users/benhoob/CondorSkim/MuonTrigger/MuHad_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged_ntuple_999999_21_6_skim.root");
pickSkimIfExists(chData2012B,"/home/users/benhoob/CondorSkim/MuonTrigger/MuHad_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------------------------
// muon trigger skim: 2012C
//----------------------------------------------------------------------------------------------------------
TChain* chData2012C = new TChain("Events");
if( runData2012C ){
cout << "Adding all 2012C data muon trigger skim" << endl;
//pickSkimIfExists(chData2012C,"/home/users/benhoob/CondorSkim/MuonTrigger/MuHad_Run2012C-PromptReco-v2_AOD/V05-03-13/merged_ntuple_199572_0_skim.root");
pickSkimIfExists(chData2012C,"/home/users/benhoob/CondorSkim/MuonTrigger/MuHad_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------------------------
// all data
//----------------------------------------------------------------------------------------------------------
TChain* chData = new TChain("Events");
if( runData ){
cout << "Adding all data" << endl;
pickSkimIfExists(chData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleMu_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/MuEG_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/SingleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/macneill/CMSSW_5_2_3_patch3_V05-02-07/ElectronHad_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/macneill/CMSSW_5_2_3_patch3_V05-02-07/Photon_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/MuHad_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// all electron data
//----------------------------------------------------------------------------------------------------------
TChain* chElData = new TChain("Events");
if( runElData ){
cout << "Adding all electron data" << endl;
pickSkimIfExists(chElData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chElData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/SingleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chElData,"/hadoop/cms/store/user/macneill/CMSSW_5_2_3_patch3_V05-02-07/ElectronHad_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// all muon data
//----------------------------------------------------------------------------------------------------------
TChain* chMuData = new TChain("Events");
if( runMuData ){
cout << "Adding all muon data" << endl;
//pickSkimIfExists(chMuData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store_data_Run2012A_SingleMu_AOD_PromptReco-v1_000_190_663_125F6FAA-C382-E111-8D4A-003048F1110E.root");
//pickSkimIfExists(chMuData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store_data_Run2012A_SingleMu_AOD_PromptReco-v1_000_190_663_BE881788-D182-E111-86AB-001D09F2462D.root");
// pickSkimIfExists(chMuData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/MuEG_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chMuData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleMu_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chMuData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/MuHad_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chMuData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// single electron data
//----------------------------------------------------------------------------------------------------------
TChain* chSingleElData = new TChain("Events");
if( runSingleElData ){
cout << "Adding single electron data" << endl;
pickSkimIfExists(chSingleElData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/SingleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// single muon data
//----------------------------------------------------------------------------------------------------------
TChain* chSingleMuData = new TChain("Events");
if( runSingleMuData ){
cout << "Adding single muon data" << endl;
//pickSkimIfExists(chSingleMuData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/*root");
//pickSkimIfExists(chSingleMuData,"/tas/benhoob/testFiles/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store_data_Run2012A_SingleMu_AOD_PromptReco-v1_000_191_718_1031CF12-DC8B-E111-9983-003048F117EC.root");
//pickSkimIfExists(chSingleMuData,"/tas/benhoob/testFiles/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store*root");
pickSkimIfExists(chSingleMuData,"/tas/benhoob/testFiles/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store*191_718*root");
//pickSkimIfExists(chSingleMuData,"/tas/benhoob/testFiles/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store*191_830*root");
//pickSkimIfExists(chSingleMuData,"/tas/benhoob/testFiles/SingleMu_Run2012A-PromptReco-v1_AOD/unmerged/store_data_Run2012A_SingleMu_AOD_PromptReco-v1_000_191_830_5E887D82-A98C-E111-9455-003048F117F6.root");
}
//----------------------------------------------------------------------------------------------------------
// ElHad data
//----------------------------------------------------------------------------------------------------------
TChain* chElHadData = new TChain("Events");
if( runElHadData ){
cout << "Adding electronhad data" << endl;
//pickSkimIfExists(chElHadData,"/hadoop/cms/store/user/cwelke/CMSSW_5_2_3_patch3_V05-02-07/ElectronHad_Run2012A-PromptReco-v1_AOD/unmerged/store_data_Run2012A_ElectronHad_AOD_PromptReco-v1_000_191_718_C448B8E1-028C-E111-9E99-003048F110BE.root");
//pickSkimIfExists(chElHadData,"/hadoop/cms/store/user/cwelke/CMSSW_5_2_3_patch3_V05-02-07/ElectronHad_Run2012A-PromptReco-v1_AOD/unmerged/*191_718*root");
pickSkimIfExists(chElHadData,"/hadoop/cms/store/user/cwelke/CMSSW_5_2_3_patch3_V05-02-07/ElectronHad_Run2012A-PromptReco-v1_AOD/unmerged/*191_830*root");
//pickSkimIfExists(chElHadData,"/hadoop/cms/store/user/macneill/CMSSW_5_2_3_patch3_V05-02-07/ElectronHad_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// doubleElectron data
//----------------------------------------------------------------------------------------------------------
TChain* chDoubleElData = new TChain("Events");
if( runDoubleElData ){
cout << "Adding double electron data" << endl;
pickSkimIfExists(chDoubleElData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// doubleMu data
//----------------------------------------------------------------------------------------------------------
TChain* chDoubleMuData = new TChain("Events");
if( runDoubleMuData ){
cout << "Adding double muon data" << endl;
pickSkimIfExists(chDoubleMuData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleMu_Run2012A-PromptReco-v1_AOD/unmerged/*root");
//pickSkimIfExists(chDoubleMuData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleMu_Run2012A-PromptReco-v1_AOD/unmerged/store_data_Run2012A_DoubleMu_AOD_PromptReco-v1_000_191_090_7C2A77D9-2087-E111-877A-003048D2BC4C.root");
}
//----------------------------------------------------------------------------------------------------------
// MuHad data
//----------------------------------------------------------------------------------------------------------
TChain* chMuHadData = new TChain("Events");
if( runMuHadData ){
cout << "Adding muhad data" << endl;
pickSkimIfExists(chMuHadData,"/hadoop/cms/store/user/jaehyeok/CMSSW_5_2_3_patch3_V05-02-07/MuHad_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// photon data
//----------------------------------------------------------------------------------------------------------
TChain* chPhotonData = new TChain("Events");
if( runPhotonData ){
cout << "Adding all photon data" << endl;
pickSkimIfExists(chPhotonData,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch3_V05-02-07/DoubleElectron_Run2012A-PromptReco-v1_AOD/unmerged/*root");
pickSkimIfExists(chPhotonData,"/hadoop/cms/store/user/macneill/CMSSW_5_2_3_patch3_V05-02-07/Photon_Run2012A-PromptReco-v1_AOD/unmerged/*root");
}
//----------------------------------------------------------------------------------------------------------
// run on samples
//----------------------------------------------------------------------------------------------------------
if (runData2012A) {
cout << "Processing all 2012A muon trigger data" << endl;
looper->ScanChain(chData2012A,"data2012A_mutrig");
cout << "Done processing all 2012A muon trigger data" << endl;
}
if (runData2012B) {
cout << "Processing all 2012B muon trigger data" << endl;
looper->ScanChain(chData2012B,"data2012B_mutrig");
cout << "Done processing all 2012B muon trigger data" << endl;
}
if (runData2012C) {
cout << "Processing all 2012C muon trigger data" << endl;
looper->ScanChain(chData2012C,"data2012C_mutrig");
cout << "Done processing all 2012C muon trigger data" << endl;
}
if (runData) {
cout << "Processing all data" << endl;
looper->ScanChain(chData,"allData");
cout << "Done processing all data" << endl;
}
if (runDoubleMuData) {
cout << "Processing double muon data" << endl;
looper->ScanChain(chDoubleMuData,"doubleMuData");
cout << "Done processing double muon data" << endl;
}
if (runSingleMuData) {
cout << "Processing single muon data" << endl;
looper->ScanChain(chSingleMuData,"singleMuData");
cout << "Done processing single muon data" << endl;
}
if (runSingleElData) {
cout << "Processing single electron data" << endl;
looper->ScanChain(chSingleElData,"singleElData");
cout << "Done processing single electron data" << endl;
}
if (runElHadData) {
cout << "Processing elhad data" << endl;
looper->ScanChain(chElHadData,"elHadData");
cout << "Done processing elhad data" << endl;
}
if (runDoubleElData) {
cout << "Processing double electron data" << endl;
looper->ScanChain(chDoubleElData,"doubleElData");
cout << "Done processing double electron data" << endl;
}
if (runMuHadData) {
cout << "Processing muhad data" << endl;
looper->ScanChain(chMuHadData,"muHadData");
cout << "Done processing muhad data" << endl;
}
if (runElData) {
cout << "Processing all electron data" << endl;
looper->ScanChain(chElData,"elData");
cout << "Done processing all electron data" << endl;
}
if (runMuData) {
cout << "Processing all muon data" << endl;
looper->ScanChain(chMuData,"muData");
cout << "Done processing all muon data" << endl;
}
if (runPhotonData) {
cout << "Processing all photon data" << endl;
looper->ScanChain(chPhotonData,"photonData");
cout << "Done processing all photon data" << endl;
}
//----------------------------------------
// save all the histograms
//----------------------------------------
const char* outFile = Form("../output/%s/singleLepton.root", version );
TDirectory *rootdir = gDirectory->GetDirectory("Rint:");
rootdir->cd();
saveHist(outFile);
deleteHistos();
gSystem->Exit(0);
}
void doMassPoint(float analysis, Option option, RunEra runEra)
{
bool blind = false;
//
// set up the looper
//
SmurfLooper *looper = new SmurfLooper(analysis, option, runEra);
//looper->setGoodRunList("../runlists/runlist_1092.txt");
const float lumi = 19467;
looper->setLumiScale(lumi, lumi);
//
// set up samples
//
int ana = (int)analysis;
SmurfSample *sample_data = new SmurfSample(option, DATA, kBlack, "Data", analysis);
// higgs samples
SmurfSample *sample_gghww = new SmurfSample(option, GGHWW, kCyan, "ggH", analysis);
SmurfSample *sample_qqhww = new SmurfSample(option, QQHWW, kCyan, "qqH", analysis);
SmurfSample *sample_zhww = new SmurfSample(option, ZHWW, kCyan, "ZH", analysis);
SmurfSample *sample_whww = new SmurfSample(option, WHWW, kCyan, "WH", analysis);
// backgrounds
SmurfSample *sample_qqww = new SmurfSample(option, QQWW, kYellow+2, "qqWW", analysis);
SmurfSample *sample_ggww = new SmurfSample(option, GGWW, kYellow+2, "ggWW", analysis);
SmurfSample *sample_vv = new SmurfSample(option, VV, kGreen, "VV", analysis);
SmurfSample *sample_top = new SmurfSample(option, TOP, kMagenta, "Top", analysis);
// note this considers only the Drell Yan part
SmurfSample *sample_dyll = new SmurfSample(option, ZLL, kBlue, "DYLL", analysis);
// NOTE. this zjets is a historial nameing..
// this corresponds to the WZ/ZZ where both leptons are from the same Z
//SmurfSample *sample_zjets = new SmurfSample(option, ZJETS, kBlue, "Zjets", analysis);
SmurfSample *sample_wjets = new SmurfSample(option, WJETSDATA, kCyan, "Wjets", analysis);
SmurfSample *sample_wjetsEle = new SmurfSample(option, WJETSELEDATA, kCyan, "WjetsE", analysis);
SmurfSample *sample_wjetsMu = new SmurfSample(option, WJETSMUDATA, kCyan, "WjetsM", analysis);
SmurfSample *sample_wgamma = new SmurfSample(option, WGAMMA, kCyan, "Wgamma", analysis);
SmurfSample *sample_wgammanorm = new SmurfSample(option, WGAMMANORM, kCyan, "Wgammanorm", analysis);
SmurfSample *sample_wg3l = new SmurfSample(option, WG3L, kCyan, "Wg3l", analysis);
SmurfSample *sample_ztt = new SmurfSample(option, ZTT, kBlue+2, "Ztt", analysis);
// Below are needed for the central shape
//SmurfSample *sample_dyll_loosemet = new SmurfSample(option, ZLLLOOSEMET, kBlue, "DYLLLooseMET", analysis);
// Below are needed for the shape variation studies
SmurfSample *sample_top_var = new SmurfSample(option, TOPALTER, kMagenta, "TopVar", analysis);
SmurfSample *sample_wwmcnlo = new SmurfSample(option, WWMCNLO, kBlue+2, "WWMCNLO", analysis);
SmurfSample *sample_wwmcnlo_up = new SmurfSample(option, WWMCNLOUP, kBlue+2, "WWMCNLOUp", analysis);
SmurfSample *sample_wwmcnlo_down = new SmurfSample(option, WWMCNLODOWN, kBlue+2, "WWMCNLODown", analysis);
//SmurfSample *sample_wjets_mc = new SmurfSample(option, WJETS, kCyan, "WjetsMC", analysis);
//SmurfSample *sample_wjets_mc_loose = new SmurfSample(option, WJETSMCLOOSE, kCyan, "WjetsMCLoose", analysis);
//SmurfSample *sample_dyll_data = new SmurfSample(option, ZLLDATA, kBlue, "DYLLDATA", analysis);
SmurfSample *sample_gghww_ref = new SmurfSample(option, GGHWWREF, kCyan, "ggHRef", analysis);
SmurfSample *sample_gghww_jhu = new SmurfSample(option, GGHWWJHU, kCyan, "ggHJHU", analysis);
bool skimwithmva = true;
// Below is the setup for running the smurf ntuples at the ww-preselection skim
if ( skimwithmva) {
char *dataDir = "/smurf/jaehyeok/data/Run2012_Summer12_SmurfV9_53X/mitf-alljets_19p5ifb_new/WW/"; // TAS
//char *dataDir = "/nfs-7/userdata/jaehyeok/smurfntuples/mitf-alljets/WW/"; // UAF
// set up the mva histogram ranges
unsigned int n = 20;
float min = -1.0;
float max = 1.0;
if ( option == HWW_OPT_SMURFMESEL) min = 0.0;
if ( analysis > 0. ) {
// higgs
for (int njet=0; njet < 2; njet++) {
if ( option == XWW_OPT_MT2DMLL_JCP )
sample_gghww->add(Form("%s/mva/%i/xww2p%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
else
sample_gghww->add(Form("%s/mva/%i/hww%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
}
sample_gghww->addShapeVariation2D(LEPEFFVAR, "histo_ggH_CMS_hww_MVALepEffBounding", false);
sample_gghww->addShapeVariation2D(LEPRESVAR, "histo_ggH_CMS_hww_MVALepResBounding", true);
sample_gghww->addShapeVariation2D(METVAR, "histo_ggH_CMS_hww_MVAMETResBounding", true);
sample_gghww->addShapeVariation2D(JETRESVAR, "histo_ggH_CMS_hww_MVAJESBounding", false);
for (int njet=0; njet < 2; njet++)
sample_qqhww->add(Form("%s/mva/%i/hww%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
sample_qqhww->addShapeVariation2D(LEPEFFVAR, "histo_qqH_CMS_hww_MVALepEffBounding", false);
sample_qqhww->addShapeVariation2D(LEPRESVAR, "histo_qqH_CMS_hww_MVALepResBounding", true);
sample_qqhww->addShapeVariation2D(METVAR, "histo_qqH_CMS_hww_MVAMETResBounding", true);
sample_qqhww->addShapeVariation2D(JETRESVAR, "histo_qqH_CMS_hww_MVAJESBounding", false);
for (int njet=0; njet < 2; njet++)
sample_zhww->add(Form("%s/mva/%i/hww%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
sample_zhww->addShapeVariation2D(LEPEFFVAR, "histo_ZH_CMS_hww_MVALepEffBounding", false);
sample_zhww->addShapeVariation2D(LEPRESVAR, "histo_ZH_CMS_hww_MVALepResBounding", true);
sample_zhww->addShapeVariation2D(METVAR, "histo_ZH_CMS_hww_MVAMETResBounding", true);
sample_zhww->addShapeVariation2D(JETRESVAR, "histo_ZH_CMS_hww_MVAJESBounding", false);
for (int njet=0; njet < 2; njet++)
sample_whww->add(Form("%s/mva/%i/hww%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
sample_whww->addShapeVariation2D(LEPEFFVAR, "histo_WH_CMS_hww_MVALepEffBounding", false);
sample_whww->addShapeVariation2D(LEPRESVAR, "histo_WH_CMS_hww_MVALepResBounding", true);
sample_whww->addShapeVariation2D(METVAR, "histo_WH_CMS_hww_MVAMETResBounding", true);
sample_whww->addShapeVariation2D(JETRESVAR, "histo_WH_CMS_hww_MVAJESBounding", false);
}
// bkgs
// qqWW
for (int njet=0; njet < 2; njet++)
sample_qqww->add(Form("%s/mva/%i/qqww_%ij.root", dataDir, int(analysis), njet));
if ( analysis < 200)
sample_qqww->addShapeVariation2D(LEPEFFVAR, "histo_qqWW_CMS_hww_MVALepEffBounding", true);
else
sample_qqww->addShapeVariation2D(LEPEFFVAR, "histo_qqWW_CMS_hww_MVALepEffBounding", false);
sample_qqww->addShapeVariation2D(LEPRESVAR, "histo_qqWW_CMS_hww_MVALepResBounding", true);
sample_qqww->addShapeVariation2D(METVAR, "histo_qqWW_CMS_hww_MVAMETResBounding", true);
sample_qqww->addShapeVariation2D(JETRESVAR, "histo_qqWW_CMS_hww_MVAJESBounding", false);
// ggWW
for (int njet=0; njet < 2; njet++)
sample_ggww->add(Form("%s/mva/%i/ggww_%ij.root", dataDir, int(analysis), njet));
if ( analysis < 200)
sample_ggww->addShapeVariation2D(LEPEFFVAR, "histo_ggWW_CMS_hww_MVALepEffBounding", true);
else
sample_ggww->addShapeVariation2D(LEPEFFVAR, "histo_ggWW_CMS_hww_MVALepEffBounding", false);
sample_ggww->addShapeVariation2D(LEPRESVAR, "histo_ggWW_CMS_hww_MVALepResBounding", true);
sample_ggww->addShapeVariation2D(METVAR, "histo_ggWW_CMS_hww_MVAMETResBounding", true);
sample_ggww->addShapeVariation2D(JETRESVAR, "histo_ggWW_CMS_hww_MVAJESBounding", false);
// Drell-Yan
for (int njet=0; njet < 2; njet++) {
sample_dyll->add(Form("%s/mva/%i/dyll_%ij.root", dataDir, int(analysis), njet));
//sample_dyll_loosemet->add(Form("%s/mva/%i/dyll_%ij.root", dataDir, int(analysis), njet));
//ample_zjets->add(Form("%s/mva/%i/dyll_%ij.root", dataDir, int(analysis), njet));
}
// Top
for (int njet=0; njet < 2; njet++) {
sample_top->add(Form("%s/mva/%i/ttbar_powheg_%ij.root", dataDir, int(analysis), njet));
sample_top->add(Form("%s/mva/%i/tw_%ij.root", dataDir, int(analysis), njet));
}
sample_top->addShapeVariation2D(LEPRESVAR, "histo_Top_CMS_hww_MVALepResBounding", true);
sample_top->addShapeVariation2D(METVAR, "histo_Top_CMS_hww_MVAMETResBounding", true);
sample_top->addShapeVariation2D(JETRESVAR, "histo_Top_CMS_hww_MVAJESBounding", false);
// VV
for (int njet=0; njet < 2; njet++) {
sample_vv->add(Form("%s/mva/%i/wz_%ij.root", dataDir, int(analysis), njet));
sample_vv->add(Form("%s/mva/%i/zz_%ij.root", dataDir, int(analysis), njet));
sample_vv->add(Form("%s/mva/%i/www_%ij.root", dataDir, int(analysis), njet));
}
sample_vv->addShapeVariation2D(LEPEFFVAR, "histo_VV_CMS_hww_MVALepEffBounding", false);
sample_vv->addShapeVariation2D(LEPRESVAR, "histo_VV_CMS_hww_MVALepResBounding", true);
sample_vv->addShapeVariation2D(METVAR, "histo_VV_CMS_hww_MVAMETResBounding", true);
sample_vv->addShapeVariation2D(JETRESVAR, "histo_VV_CMS_hww_MVAJESBounding", false);
// Wgamma : l + gamma sample
for (int njet=0; njet < 2; njet++) {
sample_wgamma->add(Form("%s/mva/%i/wgammafo_%ij.root", dataDir, int(analysis), njet));
sample_wgamma->add(Form("%s/mva/%i/zgammafo_%ij.root", dataDir, int(analysis), njet));
}
// Wgamma Normalization
for (int njet=0; njet < 2; njet++) {
sample_wgammanorm->add(Form("%s/mva/%i/wgamma_%ij.root", dataDir, int(analysis), njet));
sample_wgammanorm->add(Form("%s/mva/%i/zgamma_%ij.root", dataDir, int(analysis), njet));
}
// Wg3l ( Wgamma* )
for (int njet=0; njet < 2; njet++) {
sample_wg3l->add(Form("%s/mva/%i/wglll_%ij.root", dataDir, int(analysis), njet));
}
// Ztt
for (int njet=0; njet < 2; njet++) {
sample_ztt->add(Form("%s/mva/%i/data_ztt_%ij.root", dataDir, int(analysis), njet));
}
// Wjets
for (int njet=0; njet < 2; njet++) {
sample_wjetsEle->add(Form("%s/mva/%i/data_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/qqww_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/ggww_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/ttbar_powheg_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/tw_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/wz_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/zz_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/wgamma_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/zgamma_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/wglll_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsEle->add(Form("%s/mva/%i/dyll_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/data_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/qqww_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/ggww_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/ttbar_powheg_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/tw_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/wz_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/zz_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/wgamma_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/zgamma_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/wglll_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_wjetsMu->add(Form("%s/mva/%i/dyll_PassFail_%ij.root", dataDir, int(analysis), njet));
}
sample_wjetsEle->addShapeVariation2D(WJETSELESHAPEVAR, "histo_WjetsE_CMS_hww_MVAWEBounding", true);
sample_wjetsMu->addShapeVariation2D(WJETSMUSHAPEVAR, "histo_WjetsM_CMS_hww_MVAWMBounding", true);
// data
if ( !blind) {
for (int njet=0; njet < 2; njet++)
sample_data->add(Form("%s/mva/%i/data_%ij.root", dataDir, int(analysis), njet));
}
// for Shape variations
// wjets
for (int njet=0; njet < 2; njet++) {
//sample_wjets_mc->add(Form("%s/mva/%i/wjets_%ij.root", dataDir, int(analysis), njet));
//sample_wjets_mc_loose->add(Form("%s/mva/%i/wjets_PassFail_%ij.root", dataDir, int(analysis), njet));
sample_top_var->add(Form("%s/mva/%i/ttbar_%ij.root", dataDir, int(analysis), njet));
sample_top_var->add(Form("%s/mva/%i/tw_%ij.root", dataDir, int(analysis), njet));
sample_wwmcnlo->add(Form("%s/mva/%i/wwmcnlo_%ij.root", dataDir, int(analysis), njet));
sample_wwmcnlo_up->add(Form("%s/mva/%i/wwmcnloup_%ij.root", dataDir, int(analysis), njet));
sample_wwmcnlo_down->add(Form("%s/mva/%i/wwmcnlodown_%ij.root", dataDir, int(analysis), njet));
//sample_dyll_data->add(Form("%s/mva/%i/data_%ij.root", dataDir, int(analysis), njet));
//sample_dyll_data->add(Form("%s/mva/%i/wz_%ij.root", dataDir, int(analysis), njet));
//sample_dyll_data->add(Form("%s/mva/%i/zz_%ij.root", dataDir, int(analysis), njet));
sample_gghww_ref->add(Form("%s/mva/%i/hww%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
sample_gghww_jhu->add(Form("%s/mva/%i/xww2p%i_%ij.root", dataDir, int(analysis), int(analysis), njet));
}
}
//
// do the looping
//
looper->loop(sample_data);
looper->loop(sample_gghww);
looper->loop(sample_qqhww);
looper->loop(sample_whww);
looper->loop(sample_zhww);
looper->loop(sample_qqww);
looper->loop(sample_ggww);
looper->loop(sample_vv);
looper->loop(sample_top);
looper->loop(sample_dyll);
//looper->loop(sample_dyll_loosemet);
//looper->loop(sample_zjets);
looper->loop(sample_wjetsEle);
looper->loop(sample_wjetsMu);
looper->loop(sample_wgammanorm);
looper->loop(sample_wgamma);
looper->loop(sample_wg3l);
looper->loop(sample_ztt);
// for shape syst.
//looper->loop(sample_wjets_mc_loose);
//looper->loop(sample_wjets_mc);
looper->loop(sample_top_var);
looper->loop(sample_wwmcnlo);
looper->loop(sample_wwmcnlo_up);
looper->loop(sample_wwmcnlo_down);
looper->loop(sample_dyll_data);
looper->loop(sample_gghww_ref);
looper->loop(sample_gghww_jhu);
//
// make tables
//
std::vector<SmurfSample*> samplesToTabulate;
samplesToTabulate.push_back(sample_data);
samplesToTabulate.push_back(sample_gghww);
samplesToTabulate.push_back(sample_qqhww);
samplesToTabulate.push_back(sample_whww);
samplesToTabulate.push_back(sample_zhww);
samplesToTabulate.push_back(sample_qqww);
samplesToTabulate.push_back(sample_ggww);
samplesToTabulate.push_back(sample_vv);
samplesToTabulate.push_back(sample_top);
samplesToTabulate.push_back(sample_dyll);
samplesToTabulate.push_back(sample_zjets);
samplesToTabulate.push_back(sample_wjetsEle);
samplesToTabulate.push_back(sample_wjetsMu);
samplesToTabulate.push_back(sample_wgammanorm);
samplesToTabulate.push_back(sample_wgamma);
samplesToTabulate.push_back(sample_wg3l);
samplesToTabulate.push_back(sample_ztt);
printResultsTable(samplesToTabulate, option, true);
//
// for shape variations
//
//samplesToTabulate.push_back(sample_dyll_loosemet);
//samplesToTabulate.push_back(sample_wjets_mc);
//samplesToTabulate.push_back(sample_wjets_mc_loose);
samplesToTabulate.push_back(sample_top_var);
samplesToTabulate.push_back(sample_wwmcnlo);
samplesToTabulate.push_back(sample_wwmcnlo_up);
samplesToTabulate.push_back(sample_wwmcnlo_down);
samplesToTabulate.push_back(sample_dyll_data);
samplesToTabulate.push_back(sample_gghww_ref);
samplesToTabulate.push_back(sample_gghww_jhu);
//
// make cards
//
printf("\n\n[doAllHWW::doMassPoint] Writing cards\n");
const std::string cardDir = "../cards/hwwjcp_19p5fb/";
ShapeVar_t mva_option = (1ll<<STATVAR) | (1ll<<TOPSHAPEVAR) | (1ll<<WWSHAPEVAR) | (1ll<<DYSHAPEVAR)
| (1ll<<LEPEFFVAR) | (1ll<<METVAR) | (1ll<<LEPRESVAR) | (1ll<<JETRESVAR) | (1ll<<WJETSELESHAPEVAR) | (1ll<<WJETSMUSHAPEVAR);
for (int jetbin = 0; jetbin < 2; ++jetbin) {
printCard(samplesToTabulate, option, jetbin, analysis, cardDir, ((1<<1)|(1<<2)), mva_option, runEra);
print2DShapeHistograms(samplesToTabulate, option, jetbin, analysis, cardDir, ((1<<1)|(1<<2)), runEra);
}
//
// save histograms
//
const std::string outFile = Form("histos_hww_analysis%i_%i.root", option, int(analysis));
saveHist(outFile.c_str());
deleteHistos();
//
// clean up
//
delete looper;
delete sample_data;
delete sample_gghww;
delete sample_qqhww;
delete sample_zhww;
delete sample_whww;
delete sample_qqww;
delete sample_ggww;
delete sample_vv;
delete sample_top;
delete sample_dyll;
delete sample_ztt;
//delete sample_zjets;
delete sample_wjets;
delete sample_wjetsEle;
delete sample_wjetsMu;
delete sample_wgamma;
delete sample_wgammanorm;
delete sample_wg3l;
//delete sample_dyll_loosemet;
delete sample_top_var;
delete sample_wwmcnlo;
delete sample_wwmcnlo_up;
delete sample_wwmcnlo_down;
//delete sample_wjets_mc;
//delete sample_wjets_mc_loose;
//delete sample_dyll_data;
delete sample_gghww_ref;
delete sample_gghww_jhu;
}
void doAll() {
//
// the looper
//
gSystem->Load("libTree.so");
gSystem->Load("libPhysics.so");
gSystem->Load("libEG.so");
gSystem->Load("libMathCore.so");
gSystem->Load("/tas/dlevans/HWW2012/CMSSW_5_2_3/src/LHAPDF-5.8.92b/lib/libLHAPDF.so");
gSystem->Load("libCMS2NtupleMacrosLooper.so");
//
// create looper
//
// create a looper for a sample that was generated by using
// CTEQ6LL. Check if this is the case for your sample!
// generally, LO samples are generated with CTEQ6LL and
// NLO samples are generated with CT10. Note that in the
// CT10 case the central subset is at index 5, not 0.
MyScanChain *looper_cteq6ll = new MyScanChain("cteq6ll.LHpdf", 0);
//
// run all pdf sets
//
std::vector<std::string> pdfSets;
// CT10
pdfSets.push_back("CT10");
pdfSets.push_back("CT10as");
// MSTW
pdfSets.push_back("MSTW2008nlo68cl");
pdfSets.push_back("MSTW2008nlo68cl_asmz+68cl");
pdfSets.push_back("MSTW2008nlo68cl_asmz+68clhalf");
pdfSets.push_back("MSTW2008nlo68cl_asmz-68cl");
pdfSets.push_back("MSTW2008nlo68cl_asmz-68clhalf");
// NNPDF
pdfSets.push_back("NNPDF20_as_0116_100");
pdfSets.push_back("NNPDF20_as_0117_100");
pdfSets.push_back("NNPDF20_as_0118_100");
pdfSets.push_back("NNPDF20_100");
pdfSets.push_back("NNPDF20_as_0120_100");
pdfSets.push_back("NNPDF20_as_0121_100");
pdfSets.push_back("NNPDF20_as_0122_100");
// data sample
TChain *chain_ttbar = new TChain("Events");
chain_ttbar->Add("/tas/dlevans_data/MCNtupling/CMSSW/CMSSW_5_3_2_patch4_V05-03-23/crab/makecms2ntuple/post_processed_ntuple.root");
// do gensets
// the variation of the genset with respect to itself
// is by definition 1.0, but we'll need this number later
// stored just like all the others
looper_cteq6ll->ScanChain("ttbar", chain_ttbar, "cteq6ll");
// do other sets
for (unsigned int i = 0; i < pdfSets.size(); ++i) {
std::cout << "===== Doing =====> " << pdfSets[i] << std::endl;
looper_cteq6ll->ScanChain("ttbar", chain_ttbar, pdfSets[i]);
}
//
// write histograms
//
const std::string outFile = "results.root";
saveHist(outFile.c_str());
deleteHistos();
//
// tidy up
//
delete looper_cteq6ll;
delete chain_ttbar;
}
