void EstimateFakes(unsigned int bitmask, bool skipFWLite = false){
//cut <-> bit mask
//ID cuts -> 2**0 (1)
//Isolation cuts -> 2**1 (2) (default is both legs are isolated. Using relative isolation, TRK+CALO)
// 2**8+2**1 (Require one-only hyp lepton to be isolated. In emu, that will mean that the muon
// will have rel iso > 0.92 and the el will have 0.6 < relIso < 0.92)
// 2**9+2**1 (Require one-only hyp lepton to be isolated. In emu, that will mean that the el
// will have rel iso > 0.92 and the mu will have 0.6 < relIso < 0.92)
// 2**9+2**8+2**1 (require both leptons to have 0.6 < relIso < 0.92)
//DileptonMassVeto -> 2**2 (4)
//METcut -> 2**3 (8)
//nJets -> 2**4 (16)
//extra MuTag -> 2**5 (32)
//METveto -> 2**6 (64)
//Extra MuTag (pt>5) -> 2**7 (128)
//looseDileptonSelection, TTdil note 2008
// -> 2**10 (1024)
//fullMultipleHypsOnly -> 2**11 (2048) !!!!! Not implemented, so does nothing right now !!!!!!!
//applyZWindow cut -> 2**12 (4096)
//Opp. Sign Selection -> 2**13 (8192)
//fillMaxWeightDilOnly -> 2**14 (16284
//leptonIsolationDilSelectionTTDil08 -> 2**15 (uses trk and calo isolation seperately, reltrkIso > 0.9,
// relCaloIso > 0.9. NO OTHER CUTS BUT ISOLATION APPLIED)
//looseDilSelectionNoIsoTTDil08 -> 2**16 (basic muon preselection cuts, no isolation)
//lepton20Eta2p4DilSelection -> 2**17 (only pt and eta cuts applied to leptons)
//metBaselineSelectionTTDil108 -> 2**18 (calls passPatMet_OF20_SF30 -> corrMET > 20 (emu),
// corrMET > 30 (ee, mumu))
//dilepMassVetoCutTTDil08 -> 2**19 pretty much what it sounds like
//applyTriggersMu9orLisoE15 -> 2**20 mm -> HLT_Mu9
// ee -> HLT_Ele15_SW_L1R
// em -> HLT_Mu9 or HLT_Ele15_SW_L1R
//applyTriggersTTDil08JanTrial -> 2**21 mm -> HLT_Mu15_L1Mu7 || HLT_DoubleMu3
// ee -> HLT_IsoEle18_L1R || HLT_DoubleIsoEle12_L1R
// em -> HLT_IsoEle18_L1R || HLT_Mu15_L1Mu7 || HLT_IsoEle10_Mu10_L1R
// dilepAdditionalMassVetoCutTTDil08 -> 2**22
// corJES10ptDn -> 2**23 rescale JES by 10% up
// corJES10ptUp -> 2**24 rescale JES by 10% up
// invertZMassVeto -> 2**28 only use events in Z-veto region
// Load various tools
//gROOT->ProcessLine(Form(".x setup.C(%d)",skipFWLite));
gROOT->ProcessLine(".L setup.C");
gROOT->ProcessLine("setup(true)");
// Load and compile the looping code
gSystem->CompileMacro("ttDilCounts_looper.C", "++k", "libttDilCounts_looper");
// Flag for jet selection
// true = hyp_jet selection (15 GeV uncorrectedm eta<3)
// false = 30 GeV corrected, eta<2.4
bool oldjet=false;
// 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 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;
// Prescales
int prettdil = 1;
int prettotr = 1;
int preWW = 1;
int preWZ = 1;
int preZZ = 1;
int preWjets = 1;
int preWcharm = 1;
int preDYee = 1;
int preDYmm = 1;
int preDYtautau = 1;
int preppMuX = 1;
int preEM = 1;
int pretW = 1;
int preVQQ = 1;
// Flags for files to run over
bool runttdil = true;
bool runttotr = true;
bool runWW = true;
bool runWZ = true;
bool runZZ = true;
bool runWjets = true;
bool runWcharm = false;
bool runDYee = true;
bool runDYmm = true;
bool runDYtautau = true;
bool runppMuX = true;
bool runEM = true;
bool runtW = true;
bool runVQQ = true;
TChain* chtopdil = new TChain("Events");
pickSkimIfExists(chtopdil, "data/TTJets-madgraph_Fall08_IDEAL_V9_v2/merged*.root", "_skimSimple2020anydil");
TChain* chtopotr = new TChain("Events");
pickSkimIfExists(chtopotr, "data/TTJets-madgraph_Fall08_IDEAL_V9_v2/merged*.root", "_skimSimple2020nodil");
TChain* chww = new TChain("Events");
// pickSkimIfExists(chww, "data/WW_2l_Summer08_IDEAL_V9_v2/merged*.root", "");
pickSkimIfExists(chww, "data/WW_Summer08_IDEAL_V9_v1/merged*.root", "");
TChain* chWZ = new TChain("Events");
// pickSkimIfExists(chWZ, "data/WZ_3l_Summer08_IDEAL_V9_v2/merged*.root", ""); // can try WZ_3l-Pythia
pickSkimIfExists(chWZ, "data/WZ_incl_Summer08_IDEAL_V9_v2/merged*.root", ""); // can try WZ_3l-Pythia
TChain* chZZ = new TChain("Events");
// pickSkimIfExists(chZZ, "data/ZZ_2l2n_Summer08_IDEAL_V9_v2/merged*.root", "");
pickSkimIfExists(chZZ, "data/ZZ_Summer08_IDEAL_V9_v1/merged*.root", "");
TChain* chWjets = new TChain("Events");
pickSkimIfExists(chWjets, "data/WJets-madgraph_Fall08_IDEAL_V9_v1/merged*.root", "");
TChain* chWcharm = new TChain("Events");
pickSkimIfExists(chWcharm, "data/Wc-madgraph_Fall08_IDEAL_V9_reco-v1/merged*.root", "");
TChain* chDYtautau = new TChain("Events");
pickSkimIfExists(chDYtautau, "data/ZJets-madgraph_Fall08_IDEAL_V9_reco-v2/merged*.root", "_skimSimple2020tautau");
//the low-mass splice has no choice other than the skim
pickSkimIfExists(chDYtautau, "data/Ztautau_M20_Summer08_IDEAL_V9_v1/merged*.root_skimSimple2020_20m50", "");
TChain* chDYee = new TChain("Events");
pickSkimIfExists(chDYee, "data/ZJets-madgraph_Fall08_IDEAL_V9_reco-v2/merged*.root", "_skimSimple2020ee");
//the low-mass splice has no choice other than the skim
pickSkimIfExists(chDYee, "data/Zee_M20_Summer08_IDEAL_V9_reco-v3/merged*.root_skimSimple2020_20m50", "");
TChain* chDYmm = new TChain("Events");
pickSkimIfExists(chDYmm, "data/ZJets-madgraph_Fall08_IDEAL_V9_reco-v2/merged*.root", "_skimSimple2020mm");
//the low-mass splice has no choice other than the skim
pickSkimIfExists(chDYmm, "data/Zmumu_M20_Summer08_IDEAL_V9_reco-v2/merged*.root_skimSimple2020_20m50", "");
//ppMuX
TChain* chppMuX = new TChain("Events");
if (runppMuX) {
pickSkimIfExists(chppMuX, "data/InclusiveMuPt15/merged*.root", "_skimSimple2020");
//can try InclusiveMu5Pt50 .. figure out how to merge later
}
//ppEM
TChain* chEM = new TChain("Events");
if (runEM) {
pickSkimIfExists(chEM, "data/QCD_EMenriched_Pt20to30/merged*.root", "_skimSimple2020");
pickSkimIfExists(chEM, "data/QCD_EMenriched_Pt30to80/merged*.root", "_skimSimple2020");
pickSkimIfExists(chEM, "data/QCD_EMenriched_Pt80to170/merged*.root", "_skimSimple2020");
pickSkimIfExists(chEM, "data/QCD_BCtoE_Pt20to30/merged*.root", "_skimSimple2020");
pickSkimIfExists(chEM, "data/QCD_BCtoE_Pt30to80/merged*.root", "_skimSimple2020");
pickSkimIfExists(chEM, "data/QCD_BCtoE_Pt80to170/merged*.root", "_skimSimple2020");
}
//tW
TChain* chtW = new TChain("Events");
if (runtW) {
pickSkimIfExists(chtW, "data/SingleTop_sChannel-madgraph-LHE/merged*.root", "");
pickSkimIfExists(chtW, "data/SingleTop_tChannel-madgraph-LHE/merged*.root", "");
pickSkimIfExists(chtW, "data/SingleTop_tWChannel-madgraph-LHE/merged*.root", "");
}
//VQQ
TChain* chVQQ = new TChain("Events");
if (runVQQ) {
pickSkimIfExists(chVQQ, "data/VQQ-madgraph_Fall08_IDEAL_V9_v1/merged*.root", "");
}
// Define colors numbers:
gStyle->SetPalette(1);
enum EColor { kWhite, kBlack, kRed, kGreen, kBlue, kYellow, kMagenta, kCyan };
ttDilCounts_looper* looper = new ttDilCounts_looper();
// Process files one at a time, and color them as needed
if (runttdil) {
cout << "Processing ttbar dileptonic.. "<<endl;
looper->ScanChain(chtopdil,"ttdil", kttdil, prettdil, oldjet, bitmask);
cout << "Done Processing ttbar dileptonic.. "<<endl;
hist::color("ttdil", kYellow);
}
if (runttotr) {
cout << "Processing ttbar no-dileptons.. "<<endl;
looper->ScanChain(chtopotr,"ttotr", kttotr, prettotr, oldjet, bitmask);
hist::color("ttotr", 30);
}
if (runWW) {
cout << "Processing WW.."<<endl;
looper->ScanChain(chww,"ww", kWW, preWW, oldjet, bitmask);
hist::color("ww", kRed);
}
if (runWZ) {
cout << "Processing WZ.."<<endl;
looper->ScanChain(chWZ,"wz", kWZ, preWZ, oldjet, bitmask);
hist::color("wz", kBlue);
}
if (runZZ) {
cout << "Processing ZZ.."<<endl;
looper->ScanChain(chZZ,"zz", kZZ, preZZ, oldjet, bitmask);
hist::color("zz", kGreen);
}
if (runWjets) {
cout << "Processing Wjets.."<<endl;
looper->ScanChain(chWjets,"wjets", kWjets, preWjets, oldjet, bitmask);
hist::color("wjets", 40);
}
if (runDYtautau) {
cout << "Processing DY->tautau" << endl;
looper->ScanChain(chDYtautau,"DYtautau", kDYtautau, preDYtautau, oldjet, bitmask);
hist::color("DYtautau", kBlack);
}
if (runDYee) {
cout << "Processing DY->ee" << endl;
looper->ScanChain(chDYee,"DYee", kDYee, preDYee, oldjet, bitmask);
hist::color("DYee", kMagenta);
}
if (runDYmm) {
cout << "Processing DY->mm" << endl;
looper->ScanChain(chDYmm,"DYmm", kDYmm, preDYmm, oldjet, bitmask);
hist::color("DYmm", kCyan);
}
if (runppMuX) {
cout << "Processing ppMuX"<<endl;
looper->ScanChain(chppMuX,"ppMuX", kppMuX, preppMuX, oldjet, bitmask);
hist::color("ppMuX", 51);
}
if (runEM) {
cout << "Processing EM"<<endl;
looper->ScanChain(chEM,"EM", kEM, preEM, oldjet, bitmask);
hist::color("EM", 49);
}
if (runtW) {
cout << "Processing tW"<<endl;
looper->ScanChain(chtW,"tW", ktW, pretW, oldjet, bitmask);
hist::color("tW", 63);
}
if (runVQQ) {
cout << "Processing VQQ"<<endl;
looper->ScanChain(chVQQ,"VQQ", kVQQ, preVQQ, oldjet, bitmask);
hist::color("VQQ", 45);
}
//save all the histograms
//const char* outFile = Form("myHist_FRmu_emu_%d_%s.root", bitmask, looper->compactConfig.c_str());
//const char* outFile = Form("myHist_FRmu_05_05_emu_%d_%s.root", bitmask, looper->compactConfig.c_str());
const char* outFile = Form("myHist_FRwLowMassDY_%d_%s.root", bitmask, looper->compactConfig.c_str());
hist::saveHist(outFile);
hist::deleteHistos();
cout << "Finished with bitmask "<<bitmask << endl;
//gSystem->Exit(0);
}
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 runZLooper(char* prefix , bool isData = true, float kFactor = 1.){
TChain* ch = new TChain("Events");
//----------------------------------------------------------------------------------------
if( strcmp( prefix , "data_53X_2012A" ) == 0 ){
//pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged_ntuple_999999_1_skim.root");
//pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged_ntuple_999999_2_5_skim.root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012A-recover-06Aug2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012A-recover-06Aug2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012A-recover-06Aug2012-v1_AOD/V05-03-18_slim/merged*root");
}
else if( strcmp( prefix , "data_53X_2012B" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012B-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012B-13Jul2012-v4_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012B-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
}
else if( strcmp( prefix , "data_53X_2012C" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012C-24Aug2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012C-24Aug2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012C-24Aug2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012C-PromptReco-v2_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012C-PromptReco-v2_AOD/V05-03-18_slim/merged*root");
}
else if( strcmp( prefix , "data_53X_2012D" ) == 0 ){
// pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged_ntuple_205718_0_skim.root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
}
else if( strcmp( prefix , "data_53X_2012D_ee" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_slim/DoubleElectron_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
}
else if( strcmp( prefix , "data_53X_2012D_mm" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_slim/DoubleMu_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
}
else if( strcmp( prefix , "data_53X_2012D_em" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_slim/MuEG_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*root");
//pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple_1.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_small_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*2.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_part1_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*0.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*1.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*2.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*3.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_part2_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*4.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*5.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*6.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_part3_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*7.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*8.root");
pickSkimIfExists(ch,"Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*9.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttbar_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*root");
// pickSkimIfExists(ch,"datalinks/Summer12_53X_MC/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-18_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz2l2q_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/ZZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz2l2nu_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/ZZJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v3/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz4l_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wz2l2q_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/WZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wz3lnu_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttw_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/TTWJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttz_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/TTZJets_8TeV-madgraph_v2_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "tbz_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/TBZToLL_4F_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7C-v1/V05-03-20/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ww_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "t_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/T_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
pickSkimIfExists(ch,"Summer12_53X_MC/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "vvv_53X_slim" ) == 0 ){
pickSkimIfExists(ch,"Summer12_53X_MC/ZZZNoGstarJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged*root");
pickSkimIfExists(ch,"Summer12_53X_MC/WWWJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged*root");
pickSkimIfExists(ch,"Summer12_53X_MC/WZZNoGstarJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged*root");
pickSkimIfExists(ch,"Summer12_53X_MC/WWZNoGstarJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13_slim/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data_2012A_53X" ) == 0 ){
//pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged_ntuple_999999_6_9_skim.root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data_2012B_53X" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012B-13Jul2012-v4_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data_2012C_53X" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data_allTrig_53X" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleMu_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleElectron_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleMu_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012B-13Jul2012-v4_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleElectron_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleMu_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data_allTrig_2012C_53X" ) == 0 ){
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
//pickSkimIfExists(ch,"ZMet2012_pt2010/SingleElectron_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
//pickSkimIfExists(ch,"ZMet2012_pt2010/SingleMu_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleElectron_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/SingleMu_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "dataskim" ) == 0 ){
// 2012A
pickSkimIfExists(ch,"/nfs-6/userdata/benhoob/ZMet2012/DoubleElectron_Run2012A-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"/nfs-6/userdata/benhoob/ZMet2012/DoubleMu_Run2012A-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"/nfs-6/userdata/benhoob/ZMet2012/MuEG_Run2012A-PromptReco-v1_AOD/V05-02-27/merged*root");
// 2012B
pickSkimIfExists(ch,"/nfs-6/userdata/benhoob/ZMet2012/DoubleElectron_Run2012B-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"/nfs-6/userdata/benhoob/ZMet2012/DoubleMu_Run2012B-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"/nfs-6/userdata/benhoob/ZMet2012/MuEG_Run2012B-PromptReco-v1_AOD/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "dataskim2010" ) == 0 ){
// 2012A
//pickSkimIfExists(ch,"ZMet2012_pt2010/DoubleElectron_Run2012A-PromptReco-v1_AOD/V05-02-27/merged_ntuple_193621_0_skim.root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012A-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012A-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012A-PromptReco-v1_AOD/V05-02-27/merged*root");
// 2012B
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012B-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012B-PromptReco-v1_AOD/V05-02-27/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012B-PromptReco-v1_AOD/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data2012c" ) == 0 ){
// pickSkimIfExists(ch,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_3_2_patch4_V05-03-09/DoubleElectron_Run2012C-PromptReco-v2_AOD/merged/merged*root");
// pickSkimIfExists(ch,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_3_2_patch4_V05-03-09/DoubleMu_Run2012C-PromptReco-v2_AOD/merged/merged*root");
// pickSkimIfExists(ch,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_3_2_patch4_V05-03-09/MuEG_Run2012C-PromptReco-v2_AOD/merged/merged*root");
// 2012B
pickSkimIfExists(ch,"ZMet2012_pt2010/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-09/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/DoubleMu_Run2012C-PromptReco-v2_AOD/V05-03-09/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010/MuEG_Run2012C-PromptReco-v2_AOD/V05-03-09/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets" ) == 0 ){
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12-PU_S7_START52_V9-v2/V05-02-27/merged*1.root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12-PU_S7_START52_V9-v2/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/SingleOrDiLepton/merged*root");
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/SingleOrDiLepton/merged_ntuple_1.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "t_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/T_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wz2l2q_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/WZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ww_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz2l2q_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/ZZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wjets_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/user/benhoob/ZMet2010_pt2010/WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_full_53X" ) == 0 ){
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/SingleOrDiLepton/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/user/benhoob/ZMet2012_pt2010/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_MET50_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/user/benhoob/MET50/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "data2012cv2" ) == 0 ){
// 2012C v2
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/DoubleMu_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"ZMet2012_pt2010_nfs/MuEG_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zjets_10to50" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/DYJetsToLL_M-10To50filter_8TeV-madgraph_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "testfilter_newJEC" ) == 0 ){
//pickSkimIfExists(ch,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch4_V05-02-27/DoubleElectron_Run2012A-PromptReco-v1_AOD/merged/merged_ntuple_193334_0.root");
//pickSkimIfExists(ch,"/hadoop/cms/store/user/yanjuntu/CMSSW_5_2_3_patch4_V05-02-27/DoubleMu_Run2012A-PromptReco-v1_AOD/merged/merged_ntuple_193334_0.root");
pickSkimIfExists(ch,"/home/users/benhoob/filters/output/CMSSW_5_2_3_patch4_V05-02-27/DoubleMu_Run2012A-PromptReco-v1_AOD/merged/ZMet2012/merged_ntuple.root");
pickSkimIfExists(ch,"/home/users/benhoob/filters/output/CMSSW_5_2_3_patch4_V05-02-27/DoubleElectron_Run2012A-PromptReco-v1_AOD/merged/ZMet2012/merged_ntuple.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttbar" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/TTJets_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged*root");
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/TTJets_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged_ntuple.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttbar_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged_ntuple.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttbar_massiveb" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S6_START52_V9-v1/V05-02-28/SingleOrDiLepton/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz" ) == 0 ){
pickSkimIfExists(ch,"/nfs-6/userdata/cms2/ZZJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v3/V05-02-27/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v3/V05-02-27/merged*root");
pickSkimIfExists(ch,"/nfs-6/userdata/cms2/ZZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v3/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ww" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged*root");
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged_ntuple.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wz" ) == 0 ){
pickSkimIfExists(ch,"/nfs-6/userdata/cms2/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v2/V05-02-27/merged*root");
pickSkimIfExists(ch,"/nfs-6/userdata/cms2/WZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wz_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/ZZJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v3/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "zz4l_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttW_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/TTWJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "ttZ_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/TTZJets_8TeV-madgraph_v2_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "VVV_53X" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/WZZNoGstarJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/WWZNoGstarJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/WWWJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/ZZZNoGstarJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1/V05-03-13/merged_ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "wzsms" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/SMS-TChiwz_Mchargino-100to500_mLSP-0to400_8TeV-Pythia6Z_Summer12-START52_V9_FSIM-v1/V05-03-13/merged*root");
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/SMS-TChiwz_Mchargino-100to500_mLSP-0to400_8TeV-Pythia6Z_Summer12-START52_V9_FSIM-v1/V05-03-13/merged_ntuple_1.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "gmsb" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/Higgsino_StoreResults-Higgsino-ebe1b6443ab75fb2cd06c2581e9a7621/V05-02-28/merged*root");
//pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/Higgsino_StoreResults-Higgsino-ebe1b6443ab75fb2cd06c2581e9a7621/V05-02-28/merged_ntuple.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "gmsb_526" ) == 0 ){
//pickSkimIfExists(ch,"/hadoop/cms/store/user/algeorge/CMS2_V05-03-13/Higgsino_mu110-410_8TeV_v1_algeorge-Higgsino_mu110-410_8TeV_v1-800d6ac945c15e195cbd98c0883f9d09/ntuple_8_7_8k5.root");
pickSkimIfExists(ch,"/hadoop/cms/store/user/algeorge/CMS2_V05-03-13/Higgsino_mu110-410_8TeV_v1_algeorge-Higgsino_mu110-410_8TeV_v1-800d6ac945c15e195cbd98c0883f9d09/ntuple*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "gmsb_526_v2" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/user/algeorge/CMS2_V05-03-13/Higgsino_mu110-410_8TeV_v1_algeorge-Higgsino_mu110-410_8TeV_v1-800d6ac945c15e195cbd98c0883f9d09_v2/ntuple_*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "RelValZEE" ) == 0 ){
pickSkimIfExists(ch,"/tas/benhoob/home/ntupling/CMSSW_5_3_2_patch4/src/CMS2/NtupleMaker/test/RelValZEE_53X.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "RelValZMM" ) == 0 ){
pickSkimIfExists(ch,"/tas/benhoob/home/ntupling/CMSSW_5_3_2_patch4/src/CMS2/NtupleMaker/test/RelValZMM_53X.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "t" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/T_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged*root");
pickSkimIfExists(ch,"/hadoop/cms/store/group/snt/papers2012/Summer12MC/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_Summer12-PU_S7_START52_V9-v1/V05-02-27/merged*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "DoubleElectron_199752" ) == 0 ){
pickSkimIfExists(ch,"/tas/benhoob/testFiles/DoubleElectron_Run2012C-PromptReco-v2/V05-03-13/run199752/DoubleElectron_53X_199752.root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "DoubleMu_199752" ) == 0 ){
pickSkimIfExists(ch,"/tas/benhoob/testFiles/DoubleMu_Run2012C-PromptReco-v2/V05-03-13/run199752/ntuple_*root");
}
//----------------------------------------------------------------------------------------
else if( strcmp( prefix , "MuEG_199752" ) == 0 ){
pickSkimIfExists(ch,"/tas/benhoob/testFiles/MuEG_Run2012C-PromptReco-v2/V05-03-13/run199752/ntuple_*root");
}
//----------------------------------------------------------------------------------------
else{
cout << "ERROR: cannot find sample " << prefix << endl;
exit(0);
}
//----------------------------------------------------------------------------------------
cout << endl << endl;
cout << "Checking for corrupt files" << endl;
cout << "Entries in chain: " << ch->GetEntries() << endl;
ch->Draw("evt_run");
cout << endl << endl;
bool calculateTCMET = false; //recalculate tcmet on-the-fly?
Z_looper* myLooper = new Z_looper();
cout << "Running on sample " << prefix << endl;
myLooper->ScanChain(ch, prefix, isData, calculateTCMET, -1 ,kFactor);
}
void processBaby( TString outfileid = "tt_test", TString infile = "/hadoop/cms/store/group/snt/papers2012/Summer12_53X_MC/TTJets_SemiLeptMGDecays_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V7C-v1/V05-03-25/merged_ntuple_450.root" )
{
//---------------------------------------------------------------
// choose version, output will be written to output/[version]
//---------------------------------------------------------------
const char* version = "V00-00-02";
const char* jsonfile = "jsons/Cert_198050-207279_8TeV_19p47ifb_Collisions12_JSON_goodruns.txt";
const bool useMCSkims = true;
cout << "Version : " << version << endl;
cout << "json : " << jsonfile << endl;
// Load Everything
gSystem->Load("libTree.so");
gSystem->Load("libPhysics.so");
gSystem->Load("libEG.so");
gSystem->Load("libMathCore.so");
gSystem->Load("libMiniFWLite.so");
gSystem->Load("libsingleLeptonCORE.so");
gSystem->Load("libsingleLeptonLooper.so");
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 );
//batch mode
// looper->SetBatchMode(true);
TChain *chain = new TChain("Events");
pickSkimIfExists(chain, infile.Data());
//-------------------------------------
//set name to get comprehensible output
//-------------------------------------
char* sample;
//MC
if (infile.Contains("TTJets_MassiveBinDECAY_TuneZ2star_8TeV")) sample = Form("ttall_%s", outfileid.Data());
else if (infile.Contains("TTJets_HadronicMGDecays_8TeV-madgraph-tauola")) sample = Form("ttfake_lmgtau_%s", outfileid.Data());
else if (infile.Contains("TTJets_FullLeptMGDecays_8TeV-madgraph-tauola")) sample = Form("ttdl_lmgtau_%s", outfileid.Data());
else if (infile.Contains("TTJets_SemiLeptMGDecays_8TeV-madgraph-tauola")) sample = Form("ttsl_lmgtau_%s", outfileid.Data());
else if (infile.Contains("TTJets_HadronicMGDecays_8TeV-madgraph")) sample = Form("ttfake_lmg_%s", outfileid.Data());
else if (infile.Contains("TTJets_FullLeptMGDecays_8TeV-madgraph")) sample = Form("ttdl_lmg_%s", outfileid.Data());
else if (infile.Contains("TTJets_SemiLeptMGDecays_8TeV-madgraph")) sample = Form("ttsl_lmg_%s", outfileid.Data());
else if (infile.Contains("WbbJetsToLNu")) sample = Form("wbbjets_%s", outfileid.Data());
else if (infile.Contains("WJetsToLNu")) sample = Form("wjets_%s", outfileid.Data());
else if (infile.Contains("W1JetsToLNu")) sample = Form("w1jets_%s", outfileid.Data());
else if (infile.Contains("W2JetsToLNu")) sample = Form("w2jets_%s", outfileid.Data());
else if (infile.Contains("W3JetsToLNu")) sample = Form("w3jets_%s", outfileid.Data());
else if (infile.Contains("W4JetsToLNu")) sample = Form("w4jets_%s", outfileid.Data());
else if (infile.Contains("DYJetsToLL")) sample = Form("DYtot_%s", outfileid.Data());
else if (infile.Contains("DY1JetsToLL")) sample = Form("DY1Jtot_%s", outfileid.Data());
else if (infile.Contains("DY2JetsToLL")) sample = Form("DY2Jtot_%s", outfileid.Data());
else if (infile.Contains("DY3JetsToLL")) sample = Form("DY3Jtot_%s", outfileid.Data());
else if (infile.Contains("DY4JetsToLL")) sample = Form("DY4Jtot_%s", outfileid.Data());
else if (infile.Contains("TBarToDilepton_tW-channel")) sample = Form("tbartWchanb_2l_%s", outfileid.Data());
else if (infile.Contains("TBarToLeptons_s-channel")) sample = Form("tbarschan_l_%s", outfileid.Data());
else if (infile.Contains("TBarToLeptons_t-channel")) sample = Form("tbartchan_l_%s", outfileid.Data());
else if (infile.Contains("TBarToThadWlep_tW-channel")) sample = Form("tbartWchan_thWl_%s", outfileid.Data());
else if (infile.Contains("TBarToTlepWhad_tW-channel")) sample = Form("tbartWchan_tlWh_%s", outfileid.Data());
else if (infile.Contains("TToDilepton_tW-channel")) sample = Form("ttWchan_2l_%s", outfileid.Data());
else if (infile.Contains("TToLeptons_s-channel")) sample = Form("tschan_l_%s", outfileid.Data());
else if (infile.Contains("TToLeptons_t-channel")) sample = Form("ttchan_l_%s", outfileid.Data());
else if (infile.Contains("TToThadWlep_tW-channel")) sample = Form("ttWchan_thWl_%s", outfileid.Data());
else if (infile.Contains("TToTlepWhad_tW-channel")) sample = Form("ttWchan_tlWh_%s", outfileid.Data());
else if (infile.Contains("T_s-channel")) sample = Form("tschan_%s", outfileid.Data());
else if (infile.Contains("Tbar_s-channel")) sample = Form("tbarschan_%s", outfileid.Data());
else if (infile.Contains("T_t-channel")) sample = Form("ttchan_%s", outfileid.Data());
else if (infile.Contains("Tbar_t-channel")) sample = Form("tbartchan_%s", outfileid.Data());
else if (infile.Contains("T_tW-channel")) sample = Form("ttWchan_%s", outfileid.Data());
else if (infile.Contains("Tbar_tW-channel")) sample = Form("tbartWchan_%s", outfileid.Data());
else if (infile.Contains("WWJetsTo2L2Nu")) sample = Form("ww2l2nujets_%s", outfileid.Data());
else if (infile.Contains("ZZJetsTo4L")) sample = Form("zz4ljets_%s", outfileid.Data());
else if (infile.Contains("ZZJetsTo2L2Nu")) sample = Form("zz2l2nujets_%s", outfileid.Data());
else if (infile.Contains("ZZJetsTo2L2Q")) sample = Form("zz2l2qjets_%s", outfileid.Data());
else if (infile.Contains("WZJetsTo3LNu")) sample = Form("wz3lnujets_%s", outfileid.Data());
else if (infile.Contains("WZJetsTo2L2Q")) sample = Form("wz2l2qjets_%s", outfileid.Data());
else if (infile.Contains("WGstarToLNu2E")) sample = Form("wglnu2ejets_%s", outfileid.Data());
else if (infile.Contains("WGstarToLNu2Mu")) sample = Form("wglnu2mujets_%s", outfileid.Data());
else if (infile.Contains("WGstarToLNu2Tau")) sample = Form("wglnu2taujets_%s", outfileid.Data());
else if (infile.Contains("ZZZNoGstarJets")) sample = Form("zzzjets_%s", outfileid.Data());
else if (infile.Contains("WZZNoGstarJets")) sample = Form("wzzjets_%s", outfileid.Data());
else if (infile.Contains("WWZNoGstarJets")) sample = Form("wwzjets_%s", outfileid.Data());
else if (infile.Contains("WWWJets")) sample = Form("wwwjets_%s", outfileid.Data());
else if (infile.Contains("WWGJets")) sample = Form("wwgjets_%s", outfileid.Data());
else if (infile.Contains("TBZ")) sample = Form("tbz_%s", outfileid.Data());
else if (infile.Contains("TTZJets")) sample = Form("ttzjets_%s", outfileid.Data());
else if (infile.Contains("TTWJets")) sample = Form("ttwjets_%s", outfileid.Data());
else if (infile.Contains("TTGJets")) sample = Form("ttgjets_%s", outfileid.Data());
else if (infile.Contains("TTWWJets")) sample = Form("ttwwjets_%s", outfileid.Data());
else if (infile.Contains("TTJets_scaleup_TuneZ2star_8TeV")) sample = Form("tt_scaleup_%s", outfileid.Data());
else if (infile.Contains("TTJets_scaledown_TuneZ2star_8TeV")) sample = Form("tt_scaledw_%s", outfileid.Data());
else if (infile.Contains("TTJets_matchingup_TuneZ2star_8TeV")) sample = Form("tt_matchup_%s", outfileid.Data());
else if (infile.Contains("TTJets_matchingdown_TuneZ2star_8TeV")) sample = Form("tt_matchdw_%s", outfileid.Data());
else if (infile.Contains("TTJets_mass178_5_TuneZ2star_8TeV")) sample = Form("tt_massup_%s", outfileid.Data());
else if (infile.Contains("TTJets_mass166_5_TuneZ2star_8TeV")) sample = Form("tt_massdw_%s", outfileid.Data());
else if (infile.Contains("TT_CT10_TuneZ2star_8TeV-powheg")) sample = Form("tt_powheg_%s", outfileid.Data());
else if (infile.Contains("TT_8TeV-mcatnlo")) sample = Form("tt_mcatnlo_%s", outfileid.Data());
else if (infile.Contains("TTJets_CT10_8TeV-sherpa")) sample = Form("tt_sherpa_%s", outfileid.Data());
else if (infile.Contains("WH_WToLNu_HToBB")) sample = Form("whbb_%s", outfileid.Data());
else if (infile.Contains("SMS-T2tt")) sample = Form("T2tt_%s", outfileid.Data());
else if (infile.Contains("SMS-T2bw")) sample = Form("T2bw_%s", outfileid.Data());
else if (infile.Contains("SMS-TChiWH")) sample = Form("TChiWH_%s", outfileid.Data());
else if (infile.Contains("TChiwh_125_1")) sample = Form("TChiwh_125_1_%s", outfileid.Data());
else if (infile.Contains("TChiwh_150_1")) sample = Form("TChiwh_150_1_%s", outfileid.Data());
else if (infile.Contains("TChiwh_200_1")) sample = Form("TChiwh_200_1_%s", outfileid.Data());
else if (infile.Contains("TChiwh_250_1")) sample = Form("TChiwh_250_1_%s", outfileid.Data());
else if (infile.Contains("TChiwh_250_25")) sample = Form("TChiwh_250_25_%s", outfileid.Data());
else if (infile.Contains("TChiwh_250_50")) sample = Form("TChiwh_250_50_%s", outfileid.Data());
else if (infile.Contains("TChiwh_300_1")) sample = Form("TChiwh_300_1_%s", outfileid.Data());
else if (infile.Contains("TChiwh_350_1")) sample = Form("TChiwh_350_1_%s", outfileid.Data());
else if (infile.Contains("TChiwh_400_1")) sample = Form("TChiwh_400_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino130_bino1")) sample = Form("Wino_130_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino150_bino1")) sample = Form("Wino_150_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino175_bino1")) sample = Form("Wino_175_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino200_bino1")) sample = Form("Wino_200_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino250_bino1")) sample = Form("Wino_250_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino300_bino1")) sample = Form("Wino_300_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino350_bino1")) sample = Form("Wino_350_1_%s", outfileid.Data());
else if (infile.Contains("WinoNLSP_chargino400_bino1")) sample = Form("Wino_400_1_%s", outfileid.Data());
//Data
//single mu-had
else if (infile.Contains("MuHad_Run2012A-recover-06Aug2012-v1_AOD")) sample = Form("MuHad2012A_recover06Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("MuHad_Run2012A-13Jul2012-v1_AOD")) sample = Form("MuHad2012A_13Jul2012v1V532_%s", outfileid.Data());
//single muon
else if (infile.Contains("SingleMu_Run2012A-recover-06Aug2012-v1_AOD")) sample = Form("SingleMu2012A_recover06Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2012A-13Jul2012-v1_AOD")) sample = Form("SingleMu2012A_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2012B-13Jul2012-v1_AOD")) sample = Form("SingleMu2012B_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2012C-24Aug2012-v1_AOD")) sample = Form("SingleMu2012C_24Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2012C-PromptReco-v2_AOD")) sample = Form("SingleMu2012C_PromptRecov2V532_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2012D-PromptReco-v1_AOD")) sample = Form("SingleMu2012D_PromptRecov1V532_%s", outfileid.Data());
//single electron
else if (infile.Contains("SingleElectron_Run2012A-recover-06Aug2012-v1_AOD")) sample = Form("SingleElectron2012A_recover06Aug2012V532_%s", outfileid.Data());
else if (infile.Contains("SingleElectron_Run2012A-13Jul2012-v1_AOD")) sample = Form("SingleElectron2012A_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleElectron_Run2012B-13Jul2012-v1_AOD")) sample = Form("SingleElectron2012B_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleElectron_Run2012C-24Aug2012-v1_AOD")) sample = Form("SingleElectron2012C_24Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("SingleElectron_Run2012C-PromptReco-v2_AOD")) sample = Form("SingleElectron2012C_PromptRecov2V532_%s", outfileid.Data());
else if (infile.Contains("SingleElectron_Run2012D-PromptReco-v1_AOD")) sample = Form("SingleElectron2012D_PromptRecov1V532_%s", outfileid.Data());
//dimuon
else if (infile.Contains("DoubleMu_Run2012A-recover-06Aug2012-v1_AOD")) sample = Form("DoubleMu2012A_recover06Aug2012V532_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2012A-13Jul2012-v1_AOD")) sample = Form("DoubleMu2012A_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2012B-13Jul2012-v4_AOD")) sample = Form("DoubleMu2012B_13Jul2012v4V532_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2012C-24Aug2012-v1_AOD")) sample = Form("DoubleMu2012C_24Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2012C-PromptReco-v2_AOD")) sample = Form("DoubleMu2012C_PromptRecov2V532_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2012D-PromptReco-v1_AOD")) sample = Form("DoubleMu2012D_PromptRecov1V532_%s", outfileid.Data());
//electron+muon
else if (infile.Contains("MuEG_Run2012A-recover-06Aug2012-v1_AOD")) sample = Form("MuEG2012A_recover06Aug2012V532_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2012A-13Jul2012-v1_AOD")) sample = Form("MuEG2012A_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2012B-13Jul2012-v1_AOD")) sample = Form("MuEG2012B_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2012C-24Aug2012-v1_AOD")) sample = Form("MuEG2012C_24Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2012C-PromptReco-v2_AOD")) sample = Form("MuEG2012C_PromptRecov2V532_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2012D-PromptReco-v1_AOD")) sample = Form("MuEG2012D_PromptRecov1V532_%s", outfileid.Data());
//dielectron
else if (infile.Contains("DoubleElectron_Run2012A-recover-06Aug2012-v1_AOD")) sample = Form("DoubleElectron2012A_recover06Aug2012V532_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2012A-13Jul2012-v1_AOD")) sample = Form("DoubleElectron2012A_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2012B-13Jul2012-v1_AOD")) sample = Form("DoubleElectron2012B_13Jul2012v1V532_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2012C-24Aug2012-v1_AOD")) sample = Form("DoubleElectron2012C_24Aug2012v1V532_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2012C-PromptReco-v2_AOD")) sample = Form("DoubleElectron2012C_PromptRecov2V532_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2012D-PromptReco-v1_AOD")) sample = Form("DoubleElectron2012D_PromptRecov1V532_%s", outfileid.Data());
//otherwise
else sample = Form("boiade_%s", outfileid.Data());
//old names
// else if (infile.Contains("TTTo2L2Nu2B_7TeV-powheg-pythia6")) sample = Form("tt_notauola_%s", outfileid.Data());
// else if (infile.Contains("WW_TuneZ2")) sample = Form("ww_%s", outfileid.Data());
// else if (infile.Contains("WZ_TuneZ2")) sample = Form("wz_%s", outfileid.Data());
// else if (infile.Contains("ZZ_TuneZ2")) sample = Form("zz_%s", outfileid.Data());
// else if (infile.Contains("QCD_Pt-15to20_MuPt5Enriched")) sample = Form("MuQCD15to20_%s", outfileid.Data());
// else if (infile.Contains("QCD_Pt-20to30_MuPt5Enriched")) sample = Form("MuQCD20to30_%s", outfileid.Data());
// else if (infile.Contains("QCD_Pt-30to50_MuPt5Enriched")) sample = Form("MuQCD30to50_%s", outfileid.Data());
// else if (infile.Contains("QCD_Pt-50to80_MuPt5Enriched")) sample = Form("MuQCD50to80_%s", outfileid.Data());
// else if (infile.Contains("QCD_Pt-80to120_MuPt5Enriched")) sample = Form("MuQCD80to120_%s", outfileid.Data());
//
cout<<"sample is "<<sample<<endl;
//--------------------------------
//set luminosity to scale to
//--------------------------------
float kfactor = 1.0;
int prescale = 1;
float lumi = 1.0;
looper->ScanChain(chain, sample, kfactor, prescale, lumi);
// 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);
}
void runMakePhotonBabies(char* prefix , bool isData = true, float kFactor = 1.){
TChain* ch = new TChain("Events");
//-----------------------------------------------------------------------------------
if( strcmp( prefix , "data_53X_2012A" ) == 0 ){
//pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged_ntuple_999999_9_1_skim.root");
pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012A-recover-06Aug2012-v1_AOD/V05-03-18_slim/merged*root");
}
//-----------------------------------------------------------------------------------
else if( strcmp( prefix , "data_53X_2012B" ) == 0 ){
pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012B-13Jul2012-v1_AOD/V05-03-18_slim/merged*root");
}
//-----------------------------------------------------------------------------------
else if( strcmp( prefix , "data_53X_2012C" ) == 0 ){
pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012C-24Aug2012-v1_AOD/V05-03-18_slim/merged*root");
pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-18_slim/merged*root");
}
//-----------------------------------------------------------------------------------
else if( strcmp( prefix , "data_53X_2012D" ) == 0 ){
pickSkimIfExists(ch,"PhotonTriggerSkim_slim/DoubleElectron_Run2012D-PromptReco-v1_AOD/V05-03-18_slim/merged*root");
}
//-----------------------------------------------------------------------------------
else if( strcmp( prefix , "Photon" ) == 0 ){
pickSkimIfExists(ch,"/hadoop/cms/store/user/cwelke/CMSSW_5_2_3_patch4_V05-02-27/Photon_Run2012A-PromptReco-v1_AOD/unmerged/store*root");
pickSkimIfExists(ch,"/hadoop/cms/store/user/cwelke/CMSSW_5_2_3_patch4_V05-02-27/SinglePhoton_Run2012B-PromptReco-v1_AOD/unmerged/store*root");
}
//-----------------------------------------------------------------------------------
else if( strcmp( prefix , "DoubleElectron" ) == 0 ){
pickSkimIfExists(ch,"PhotonTriggerSkim/DoubleElectron_Run2012A-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"PhotonTriggerSkim/DoubleElectron_Run2012A-recover-06Aug2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"PhotonTriggerSkim/DoubleElectron_Run2012B-13Jul2012-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"PhotonTriggerSkim/DoubleElectron_Run2012C-PromptReco-v1_AOD/V05-03-13/merged*root");
pickSkimIfExists(ch,"PhotonTriggerSkim/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
}
//-----------------------------------------------------------------------------------
else if( strcmp( prefix , "DoubleElectron_2012Cv2" ) == 0 ){
pickSkimIfExists(ch,"PhotonTriggerSkim/DoubleElectron_Run2012C-PromptReco-v2_AOD/V05-03-13/merged*root");
}
//-----------------------------------------------------------------------------------
else{
cout << "ERROR: cannot find sample " << prefix << endl;
exit(0);
}
//-----------------------------------------------------------------------------------
bool calculateTCMET = false; //recalculate tcmet on-the-fly?
cout << endl;
cout << "Checking for corrupt files" << endl;
cout << "Entries " << ch->GetEntries() << endl;
ch->Draw("evt_run");
cout << endl;
makePhotonBabies* myLooper = new makePhotonBabies();
cout << "Running on sample " << prefix << endl;
myLooper->ScanChain(ch, prefix, isData, calculateTCMET, -1 ,kFactor);
}
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 processBaby( TString outfileid = "merged_ntuple_35", TString infile = "/nfs-7/userdata/cms2/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11-PU_S4_START42_V11-v1/V04-02-29_singleLepton/merged_ntuple_35.root" )
//void processBaby( TString outfileid = "merged_ntuple_180296_0", TString infile = "/hadoop/cms/store/user/yanjuntu/CMSSW_4_2_7_patch1_V04-02-34/DoubleMu_Run2011B-PromptReco-v1_AOD/CMSSW_4_2_7_patch1_V04-02-34_merged/V04-02-34//merged_ntuple_180296_0.root")
//void processBaby( TString outfileid = "merged_ntuple", TString infile = "/hadoop/cms/store/user/jaehyeok/CMSSW_4_2_7_patch1_V04-02-33/SingleMu_Run2011A-PromptReco-v6_AOD/CMSSW_4_2_7_patch1_V04-02-33_merged/V04-02-33//merged_ntuple_172620_0.root" )
{
//---------------------------------------------------------------
// choose version, output will be written to output/[version]
//---------------------------------------------------------------
const char* version = "V00-03-04";
const char* jsonfile = "jsons/Cert_160404-180252_7TeV_mergePromptMay10Aug5_JSON_goodruns.txt";
const bool useMCSkims = true;
cout << "Version : " << version << endl;
cout << "json : " << jsonfile << endl;
// Load Everything
gSystem->Load("libTree.so");
gSystem->Load("libPhysics.so");
gSystem->Load("libEG.so");
gSystem->Load("libMathCore.so");
gSystem->Load("libMiniFWLite.so");
gSystem->Load("libsingleLeptonCORE.so");
gSystem->Load("libsingleLeptonLooper.so");
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 );
//batch mode
// looper->SetBatchMode(true);
TChain *chain = new TChain("Events");
pickSkimIfExists(chain, infile.Data());
//-------------------------------------
//set name to get comprehensible output
//-------------------------------------
char* sample;
//MC
if (infile.Contains("TTJets_TuneZ2_7TeV")) sample = Form("ttall_%s", outfileid.Data());
else if (infile.Contains("TTjets_TuneZ2_scaleup_7TeV")) sample = Form("tt_scaleup_%s", outfileid.Data());
else if (infile.Contains("TTjets_TuneZ2_scaledown_7TeV")) sample = Form("tt_scaledw_%s", outfileid.Data());
else if (infile.Contains("TTjets_TuneZ2_matchingup_7TeV")) sample = Form("tt_matchup_%s", outfileid.Data());
else if (infile.Contains("TTjets_TuneZ2_matchingdown_7TeV")) sample = Form("tt_matchdw_%s", outfileid.Data());
else if (infile.Contains("TTjets_TuneZ2_mass178_5_7TeV")) sample = Form("tt_massup_%s", outfileid.Data());
else if (infile.Contains("TTjets_TuneZ2_mass166_5_7TeV")) sample = Form("tt_massdw_%s", outfileid.Data());
else if (infile.Contains("TT_TuneZ2_7TeV-pythia6-tauola")) sample = Form("tt_pythia_%s", outfileid.Data());
else if (infile.Contains("TT_TuneZ2_7TeV-mcatnlo")) sample = Form("tt_mcatnlo_%s", outfileid.Data());
else if (infile.Contains("TT_TuneZ2_7TeV-powheg-tauola")) sample = Form("tt_powheg_%s", outfileid.Data());
else if (infile.Contains("TTTo2L2Nu2B_7TeV-powheg-pythia6")) sample = Form("tt_notauola_%s", outfileid.Data());
else if (infile.Contains("WJetsToLNu")) sample = Form("wjets_%s", outfileid.Data());
else if (infile.Contains("DYJetsToLL")) sample = Form("DYtot_%s", outfileid.Data());
else if (infile.Contains("T_TuneZ2_s-channel")) sample = Form("tschan_%s", outfileid.Data());
else if (infile.Contains("Tbar_TuneZ2_s-channel")) sample = Form("tbarschan_%s", outfileid.Data());
else if (infile.Contains("T_TuneZ2_t-channel")) sample = Form("ttchan_%s", outfileid.Data());
else if (infile.Contains("Tbar_TuneZ2_t-channel")) sample = Form("tbartchan_%s", outfileid.Data());
else if (infile.Contains("T_TuneZ2_tW-channel")) sample = Form("ttWchan_%s", outfileid.Data());
else if (infile.Contains("Tbar_TuneZ2_tW-channel")) sample = Form("tbartWchan_%s", outfileid.Data());
else if (infile.Contains("WW_TuneZ2")) sample = Form("ww_%s", outfileid.Data());
else if (infile.Contains("WZ_TuneZ2")) sample = Form("wz_%s", outfileid.Data());
else if (infile.Contains("ZZ_TuneZ2")) sample = Form("zz_%s", outfileid.Data());
else if (infile.Contains("QCD_Pt-15to20_MuPt5Enriched")) sample = Form("MuQCD15to20_%s", outfileid.Data());
else if (infile.Contains("QCD_Pt-20to30_MuPt5Enriched")) sample = Form("MuQCD20to30_%s", outfileid.Data());
else if (infile.Contains("QCD_Pt-30to50_MuPt5Enriched")) sample = Form("MuQCD30to50_%s", outfileid.Data());
else if (infile.Contains("QCD_Pt-50to80_MuPt5Enriched")) sample = Form("MuQCD50to80_%s", outfileid.Data());
else if (infile.Contains("QCD_Pt-80to120_MuPt5Enriched")) sample = Form("MuQCD80to120_%s", outfileid.Data());
else if (infile.Contains("SMS-T2tt")) sample = Form("T2tt_%s", outfileid.Data());
else if (infile.Contains("SMS-T2bw")) sample = Form("T2bw_%s", outfileid.Data());
else if (infile.Contains("WWJetsTo2L2Nu")) sample = Form("ww2l2nujets_%s", outfileid.Data());
//Data
//single muon
else if (infile.Contains("SingleMu_Run2011A-PromptReco-v4_AOD") && infile.Contains("V04-02-33")) sample = Form("SingleMu2011A_PromptRecov4V33_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2011A-PromptReco-v6_AOD") && infile.Contains("V04-02-33")) sample = Form("SingleMu2011A_PromptRecov6V33_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-33")) sample = Form("SingleMu2011B_PromptRecov1V33_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-34")) sample = Form("SingleMu2011B_PromptRecov1V34_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2011A-May10ReReco-v1_AOD") && infile.Contains("V04-02-33")) sample = Form("SingleMu2011A-May10ReRecov1V33_%s", outfileid.Data());
else if (infile.Contains("SingleMu_Run2011A-05Aug2011-v1_AOD") && infile.Contains("V04-02-33")) sample = Form("SingleMu2011A-05Aug2011v1V33_%s", outfileid.Data());
//dimuon
else if (infile.Contains("DoubleMu_Run2011A-PromptReco-v4_AOD") && infile.Contains("V04-02-20")) sample = Form("DoubleMu2011A_PromptRecov4V20_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2011A-PromptReco-v6_AOD") && infile.Contains("V04-02-30")) sample = Form("DoubleMu2011A_PromptRecov6V30_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-30")) sample = Form("DoubleMu2011B_PromptRecov1V30_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-34")) sample = Form("DoubleMu2011B_PromptRecov1V34_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2011A-May10ReReco-v1_AOD") && infile.Contains("V04-02-20")) sample = Form("DoubleMu2011A-May10ReRecov1V20_%s", outfileid.Data());
else if (infile.Contains("DoubleMu_Run2011A-05Aug2011-v1_AOD") && infile.Contains("V04-02-30")) sample = Form("DoubleMu2011A-05Aug2011v1V30_%s", outfileid.Data());
//electron+muon
else if (infile.Contains("MuEG_Run2011A-PromptReco-v4_AOD") && infile.Contains("V04-02-20")) sample = Form("MuEG2011A_PromptRecov4V20_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2011A-PromptReco-v6_AOD") && infile.Contains("V04-02-30")) sample = Form("MuEG2011A_PromptRecov6V30_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-30")) sample = Form("MuEG2011B_PromptRecov1V30_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-34")) sample = Form("MuEG2011B_PromptRecov1V34_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2011A-May10ReReco-v1_AOD") && infile.Contains("V04-02-20")) sample = Form("MuEG2011A-May10ReRecov1V20_%s", outfileid.Data());
else if (infile.Contains("MuEG_Run2011A-05Aug2011-v1_AOD") && infile.Contains("V04-02-30")) sample = Form("MuEG2011A-05Aug2011v1V30_%s", outfileid.Data());
//dielectron
else if (infile.Contains("DoubleElectron_Run2011A-PromptReco-v4_AOD") && infile.Contains("V04-02-20")) sample = Form("DoubleElectron2011A_PromptRecov4V20_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2011A-PromptReco-v6_AOD") && infile.Contains("V04-02-30")) sample = Form("DoubleElectron2011A_PromptRecov6V30_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-30")) sample = Form("DoubleElectron2011B_PromptRecov1V30_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2011B-PromptReco-v1_AOD") && infile.Contains("V04-02-34")) sample = Form("DoubleElectron2011B_PromptRecov1V34_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2011A-May10ReReco-v1_AOD") && infile.Contains("V04-02-20")) sample = Form("DoubleElectron2011A-May10ReRecov1V20_%s", outfileid.Data());
else if (infile.Contains("DoubleElectron_Run2011A-05Aug2011-v1_AOD") && infile.Contains("V04-02-30")) sample = Form("DoubleElectron2011A-05Aug2011v1V30_%s", outfileid.Data());
//otherwise
else sample = Form("boiade_%s", outfileid.Data());
cout<<"sample is "<<sample<<endl;
//--------------------------------
//set luminosity to scale to
//--------------------------------
float kfactor = 1.0;
int prescale = 1;
float lumi = 1.0;
looper->ScanChain(chain, sample, kfactor, prescale, lumi);
// gSystem->Exit(0);
}
