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 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);
}
