From b1292bcaec580c2889529a841d30c496bfc303ab Mon Sep 17 00:00:00 2001 From: David Date: Thu, 4 Feb 2016 09:34:56 -0600 Subject: [PATCH 1/3] changes to Macro.cc testDir and include file structure. --- plugins/Macros.cc | 152 +++++++++++++++++++++++++--------------------- 1 file changed, 83 insertions(+), 69 deletions(-) diff --git a/plugins/Macros.cc b/plugins/Macros.cc index 3667d73..6279088 100644 --- a/plugins/Macros.cc +++ b/plugins/Macros.cc @@ -4,10 +4,25 @@ #include #include #include +/* #include "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/plugins/Sample.cc" #include "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/plugins/CutClass.cc" #include "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/interface/TreeReader.h" #include "TGraphAsymmErrors.h" +*/ +#include "Sample.cc" +#include "CutClass.cc" +#include "../interface/TreeReader.h" +#include "TGraphAsymmErrors.h" + +std::string testDir = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_10_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/"; +std::string bgTestDir = testDir; +std::string x53String = "X53X53"; +/* +std::string testDir = "/uscms/home/dzou/CMSSW_7_3_0/src/X53ThirteenTeVAnalysisCode/test/"; //global variable for test directory used for samples, etc. +std::string testDir = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_10_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/"; +std::string x53String = "X53"; +*/ std::vector getVariableVec(){ @@ -34,7 +49,7 @@ std::vector getVariableVec(){ Variable* ak4jet1phi = new Variable("AK4Jet1Phi",20,-3.5,3.5,"Leading Jet #phi","N_{Events}"); vVar.push_back(ak4jet1phi); - Variable* ak4jet2pt = new Variable("AK4Jet2Pt",15,0,600,"subLeading Jet p_{T} (GeV)","N_{Events}"); +< Variable* ak4jet2pt = new Variable("AK4Jet2Pt",15,0,600,"subLeading Jet p_{T} (GeV)","N_{Events}"); vVar.push_back(ak4jet2pt); Variable* ak4jet2eta = new Variable("AK4Jet2Eta",25,-5,5,"subLeading Jet #eta","N_{Events}"); vVar.push_back(ak4jet2eta); @@ -180,8 +195,7 @@ std::vector getNPHistos(int nMu){ Sample* getDataSample(std::string cut,std::string elID, std::string muID){ - - std::string filename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/Data_Mu"+muID+"_El"+elID+".root"; + std::string filename = testDir.c_str()+muID+"_El"+elID+".root"; TFile* tfile = new TFile(filename.c_str()); Sample* dataSample = new Sample("Data",tfile,1,1,cut,kBlack,1); return dataSample; @@ -215,29 +229,29 @@ std::vector getSigSampleVecForPlots(std::string cut, float lumi, std::s std::vector vSigSamples; - std::string lh700name ="/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m700LH_Mu"+muID+"_El"+elID+".root"; + std::string lh700name =testDir.c_str()+"X53X53m700LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Lfile = new TFile(lh700name.c_str()); Sample* x53x53m700L = new Sample(vSigNames.at(0),x53x53700Lfile,vWeights.at(0),vXsec.at(0),cut,kGreen,2); vSigSamples.push_back(x53x53m700L); - std::string rh700name = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m700RH_Mu"+muID+"_El"+elID+".root"; + std::string rh700name = testDir.c_str()+"X53X53m700RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Rfile = new TFile(rh700name.c_str()); Sample* x53x53m700R = new Sample(vSigNames.at(1),x53x53700Rfile,vWeights.at(1),vXsec.at(1),cut,kGreen,1); vSigSamples.push_back(x53x53m700R); - /* std::string lh1000name = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1000LH_Mu"+muID+"_El"+elID+".root"; + /* std::string lh1000name = testDir.c_str()+"X53X53m1000LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Lfile = new TFile(lh1000name.c_str()); Sample* x53x53m1000L = new Sample(vSigNames.at(2),x53x531000Lfile,vWeights.at(2),vXsec.at(2),cut,kCyan,2); vSigSamples.push_back(x53x53m1000L); - std::string rh1000name = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1000RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1000name = testDir.c_str()+"X53X53m1000RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Rfile = new TFile(rh1000name.c_str()); Sample* x53x53m1000R = new Sample(vSigNames.at(3),x53x531000Rfile,vWeights.at(3),vXsec.at(3),cut,kCyan,1); vSigSamples.push_back(x53x53m1000R);*/ - /* TFile* x53x531300Lfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1300LH.root"); + /* TFile* x53x531300Lfile = new TFile(testDir.c_str()+"X53X53m1300LH.root"); Sample* x53x53m1300L = new Sample(vSigNames.at(4),x53x531300Lfile,vWeights.at(4),vXsec.at(4),cut,kBlue,2); - TFile* x53x531300Rfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1300RH.root"); + TFile* x53x531300Rfile = new TFile(testDir.c_str()+"X53X53m1300RH.root"); Sample* x53x53m1300R = new Sample(vSigNames.at(5),x53x531300Rfile,vWeights.at(5),vXsec.at(5),cut,kBlue,1); */ return vSigSamples; @@ -279,82 +293,82 @@ std::vector getSigSampleVecForTable(std::string cut, float lumi, std::s } std::vector vSigSamples; - std::string lh700 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m700LH_Mu"+muID+"_El"+elID+".root"; + std::string lh700 = testDir.c_str()+"X53X53m700LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Lfile = new TFile(lh700.c_str()); Sample* x53x53m700L = new Sample(vSigNames.at(0),x53x53700Lfile,vWeights.at(0),vXsec.at(0),cut,kGreen,2); vSigSamples.push_back(x53x53m700L); - std::string rh700 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m700RH_Mu"+muID+"_El"+elID+".root"; + std::string rh700 = testDir.c_str()+"X53X53m700RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Rfile = new TFile(rh700.c_str()); Sample* x53x53m700R = new Sample(vSigNames.at(1),x53x53700Rfile,vWeights.at(1),vXsec.at(1),cut,kGreen,1); vSigSamples.push_back(x53x53m700R); - std::string lh800 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m800LH_Mu"+muID+"_El"+elID+".root"; + std::string lh800 = testDir.c_str()+"X53X53m800LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Lfile = new TFile(lh800.c_str()); Sample* x53x53m800L = new Sample(vSigNames.at(2),x53x53800Lfile,vWeights.at(2),vXsec.at(2),cut,kBlue,2); vSigSamples.push_back(x53x53m800L); - std::string rh800 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m800RH_Mu"+muID+"_El"+elID+".root"; + std::string rh800 = testDir.c_str()+"X53X53m800RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Rfile = new TFile(rh800.c_str()); Sample* x53x53m800R = new Sample(vSigNames.at(3),x53x53800Rfile,vWeights.at(3),vXsec.at(3),cut,kBlue,1); vSigSamples.push_back(x53x53m800R); - std::string lh900 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m900LH_Mu"+muID+"_El"+elID+".root"; + std::string lh900 = testDir.c_str()+"X53X53m900LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Lfile = new TFile(lh900.c_str()); Sample* x53x53m900L = new Sample(vSigNames.at(4),x53x53900Lfile,vWeights.at(4),vXsec.at(4),cut,kGreen,2); //vSigSamples.push_back(x53x53m900L); - std::string rh900 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m900RH_Mu"+muID+"_El"+elID+".root"; + std::string rh900 = testDir.c_str()+"X53X53m900RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Rfile = new TFile(rh900.c_str()); Sample* x53x53m900R = new Sample(vSigNames.at(5),x53x53900Rfile,vWeights.at(5),vXsec.at(5),cut,kGreen,1); vSigSamples.push_back(x53x53m900R); - std::string lh1000 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1000LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1000 = testDir.c_str()+"X53X53m1000LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Lfile = new TFile(lh1000.c_str()); Sample* x53x53m1000L = new Sample(vSigNames.at(6),x53x531000Lfile,vWeights.at(6),vXsec.at(6),cut,kCyan,2); //vSigSamples.push_back(x53x53m1000L); - std::string rh1000 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1000RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1000 = testDir.c_str()+"X53X53m1000RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Rfile = new TFile(rh1000.c_str()); Sample* x53x53m1000R = new Sample(vSigNames.at(7),x53x531000Rfile,vWeights.at(7),vXsec.at(7),cut,kCyan,1); vSigSamples.push_back(x53x53m1000R); - std::string lh1100 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1100LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1100 = testDir.c_str()+"X53X53m1100LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Lfile = new TFile(lh1100.c_str()); Sample* x53x53m1100L = new Sample(vSigNames.at(8),x53x531100Lfile,vWeights.at(8),vXsec.at(8),cut,kBlue,2); vSigSamples.push_back(x53x53m1100L); - std::string rh1100 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1100RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1100 = testDir.c_str()+"X53X53m1100RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Rfile = new TFile(rh1100.c_str()); Sample* x53x53m1100R = new Sample(vSigNames.at(9),x53x531100Rfile,vWeights.at(9),vXsec.at(9),cut,kBlue,1); vSigSamples.push_back(x53x53m1100R); - std::string lh1300 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1300LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1300 = testDir.c_str()+"X53X53m1300LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Lfile = new TFile(lh1300.c_str()); Sample* x53x53m1300L = new Sample(vSigNames.at(10),x53x531300Lfile,vWeights.at(10),vXsec.at(10),cut,kBlue,2); vSigSamples.push_back(x53x53m1300L); - std::string rh1300 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1300RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1300 = testDir.c_str()+"X53X53m1300RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Rfile = new TFile(rh1300.c_str()); Sample* x53x53m1300R = new Sample(vSigNames.at(11),x53x531300Rfile,vWeights.at(11),vXsec.at(11),cut,kBlue,1); //vSigSamples.push_back(x53x53m1300R); - std::string lh1400 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1400LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1400 = testDir.c_str()+"X53X53m1400LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Lfile = new TFile(lh1400.c_str()); Sample* x53x53m1400L = new Sample(vSigNames.at(12),x53x531400Lfile,vWeights.at(12),vXsec.at(12),cut,kBlue,2); vSigSamples.push_back(x53x53m1400L); - std::string rh1400 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1400RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1400 = testDir.c_str()+"X53X53m1400RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Rfile = new TFile(rh1400.c_str()); Sample* x53x53m1400R = new Sample(vSigNames.at(13),x53x531400Rfile,vWeights.at(13),vXsec.at(13),cut,kBlue,1); vSigSamples.push_back(x53x53m1400R); - std::string lh1500 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1500LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1500 = testDir.c_str()+"X53X53m1500LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Lfile = new TFile(lh1500.c_str()); Sample* x53x53m1500L = new Sample(vSigNames.at(14),x53x531500Lfile,vWeights.at(14),vXsec.at(14),cut,kBlue,2); vSigSamples.push_back(x53x53m1500L); - std::string rh1500 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1500RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1500 = testDir.c_str()+"X53X53m1500RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Rfile = new TFile(rh1500.c_str()); Sample* x53x53m1500R = new Sample(vSigNames.at(15),x53x531500Rfile,vWeights.at(15),vXsec.at(15),cut,kBlue,1); vSigSamples.push_back(x53x53m1500R); - /* TFile* x53x531600Lfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1600LH_Mu"+muID+"_El"+elID+".root"); + /* TFile* x53x531600Lfile = new TFile(testDir.c_str()+"X53X53m1600LH_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600L = new Sample(vSigNames.at(16),x53x531600Lfile,vWeights.at(16),vXsec.at(16),cut,kBlue,2); vSigSamples.push_back(x53x53m1600L); - TFile* x53x531600Rfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1600RH_Mu"+muID+"_El"+elID+".root"); + TFile* x53x531600Rfile = new TFile(testDir.c_str()+"X53X53m1600RH_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600R = new Sample(vSigNames.at(17),x53x531600Rfile,vWeights.at(17),vXsec.at(17),cut,kBlue,1); vSigSamples.push_back(x53x53m1600R);*/ @@ -401,91 +415,91 @@ std::vector getInclusiveSigSampleVecForTable(std::string cut, float lum } std::vector vSigSamples; - std::string lh700 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m700LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh700 = testDir.c_str()+"X53X53m700LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Lfile = new TFile(lh700.c_str()); Sample* x53x53m700L = new Sample(vSigNames.at(0),x53x53700Lfile,vWeights.at(0),vXsec.at(0),cut,kGreen,2); vSigSamples.push_back(x53x53m700L); - std::string rh700 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m700RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh700 = testDir.c_str()+"X53X53m700RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Rfile = new TFile(rh700.c_str()); Sample* x53x53m700R = new Sample(vSigNames.at(1),x53x53700Rfile,vWeights.at(1),vXsec.at(1),cut,kGreen,1); vSigSamples.push_back(x53x53m700R); - std::string lh800 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m800LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh800 = testDir.c_str()+"X53X53m800LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Lfile = new TFile(lh800.c_str()); Sample* x53x53m800L = new Sample(vSigNames.at(2),x53x53800Lfile,vWeights.at(2),vXsec.at(2),cut,kBlue,2); vSigSamples.push_back(x53x53m800L); - std::string rh800 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m800RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh800 = testDir.c_str()+"X53X53m800RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Rfile = new TFile(rh800.c_str()); Sample* x53x53m800R = new Sample(vSigNames.at(3),x53x53800Rfile,vWeights.at(3),vXsec.at(3),cut,kBlue,1); vSigSamples.push_back(x53x53m800R); - std::string lh900 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m900LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh900 = testDir.c_str()+"X53X53m900LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Lfile = new TFile(lh900.c_str()); Sample* x53x53m900L = new Sample(vSigNames.at(4),x53x53900Lfile,vWeights.at(4),vXsec.at(4),cut,kGreen,2); vSigSamples.push_back(x53x53m900L); - std::string rh900 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m900RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh900 = testDir.c_str()+"X53X53m900RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Rfile = new TFile(rh900.c_str()); Sample* x53x53m900R = new Sample(vSigNames.at(5),x53x53900Rfile,vWeights.at(5),vXsec.at(5),cut,kGreen,1); vSigSamples.push_back(x53x53m900R); - std::string lh1000 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1000LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1000 = testDir.c_str()+"X53X53m1000LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Lfile = new TFile(lh1000.c_str()); Sample* x53x53m1000L = new Sample(vSigNames.at(6),x53x531000Lfile,vWeights.at(6),vXsec.at(6),cut,kCyan,2); vSigSamples.push_back(x53x53m1000L); - std::string rh1000 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1000RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1000 = testDir.c_str()+"X53X53m1000RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Rfile = new TFile(rh1000.c_str()); Sample* x53x53m1000R = new Sample(vSigNames.at(7),x53x531000Rfile,vWeights.at(7),vXsec.at(7),cut,kCyan,1); vSigSamples.push_back(x53x53m1000R); - std::string lh1100 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1100LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1100 = testDir.c_str()+"X53X53m1100LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Lfile = new TFile(lh1100.c_str()); Sample* x53x53m1100L = new Sample(vSigNames.at(8),x53x531100Lfile,vWeights.at(8),vXsec.at(8),cut,kBlue,2); vSigSamples.push_back(x53x53m1100L); - std::string rh1100 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1100RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1100 = testDir.c_str()+"X53X53m1100RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Rfile = new TFile(rh1100.c_str()); Sample* x53x53m1100R = new Sample(vSigNames.at(9),x53x531100Rfile,vWeights.at(9),vXsec.at(9),cut,kBlue,1); vSigSamples.push_back(x53x53m1100R); - std::string lh1200 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1200LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1200 = testDir.c_str()+"X53X53m1200LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531200Lfile = new TFile(lh1200.c_str()); Sample* x53x53m1200L = new Sample(vSigNames.at(10),x53x531200Lfile,vWeights.at(10),vXsec.at(10),cut,kBlue,2); vSigSamples.push_back(x53x53m1200L); - std::string rh1200 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1200RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1200 = testDir.c_str()+"X53X53m1200RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531200Rfile = new TFile(rh1200.c_str()); Sample* x53x53m1200R = new Sample(vSigNames.at(11),x53x531200Rfile,vWeights.at(11),vXsec.at(11),cut,kBlue,1); vSigSamples.push_back(x53x53m1200R); - std::string lh1300 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1300LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1300 = testDir.c_str()+"X53X53m1300LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Lfile = new TFile(lh1300.c_str()); Sample* x53x53m1300L = new Sample(vSigNames.at(12),x53x531300Lfile,vWeights.at(12),vXsec.at(12),cut,kBlue,2); vSigSamples.push_back(x53x53m1300L); - std::string rh1300 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1300RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1300 = testDir.c_str()+"X53X53m1300RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Rfile = new TFile(rh1300.c_str()); Sample* x53x53m1300R = new Sample(vSigNames.at(13),x53x531300Rfile,vWeights.at(13),vXsec.at(13),cut,kBlue,1); vSigSamples.push_back(x53x53m1300R); - std::string lh1400 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1400LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1400 = testDir.c_str()+"X53X53m1400LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Lfile = new TFile(lh1400.c_str()); Sample* x53x53m1400L = new Sample(vSigNames.at(14),x53x531400Lfile,vWeights.at(14),vXsec.at(14),cut,kBlue,2); vSigSamples.push_back(x53x53m1400L); - std::string rh1400 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1400RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1400 = testDir.c_str()+"X53X53m1400RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Rfile = new TFile(rh1400.c_str()); Sample* x53x53m1400R = new Sample(vSigNames.at(15),x53x531400Rfile,vWeights.at(15),vXsec.at(15),cut,kBlue,1); vSigSamples.push_back(x53x53m1400R); - std::string lh1500 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1500LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1500 = testDir.c_str()+"X53X53m1500LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Lfile = new TFile(lh1500.c_str()); Sample* x53x53m1500L = new Sample(vSigNames.at(16),x53x531500Lfile,vWeights.at(16),vXsec.at(16),cut,kBlue,2); vSigSamples.push_back(x53x53m1500L); - std::string rh1500 = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1500RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1500 = testDir.c_str()+"X53X53m1500RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Rfile = new TFile(rh1500.c_str()); Sample* x53x53m1500R = new Sample(vSigNames.at(17),x53x531500Rfile,vWeights.at(17),vXsec.at(17),cut,kBlue,1); vSigSamples.push_back(x53x53m1500R); - /* TFile* x53x531600Lfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1600LH_Inc_Mu"+muID+"_El"+elID+".root"); + /* TFile* x53x531600Lfile = new TFile(testDir.c_str()+"X53X53m1600LH_Inc_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600L = new Sample(vSigNames.at(16),x53x531600Lfile,vWeights.at(16),vXsec.at(16),cut,kBlue,2); vSigSamples.push_back(x53x53m1600L); - TFile* x53x531600Rfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/X53X53m1600RH_Inc_Mu"+muID+"_El"+elID+".root"); + TFile* x53x531600Rfile = new TFile(testDir.c_str()+"X53X53m1600RH_Inc_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600R = new Sample(vSigNames.at(17),x53x531600Rfile,vWeights.at(17),vXsec.at(17),cut,kBlue,1); vSigSamples.push_back(x53x53m1600R);*/ @@ -512,23 +526,23 @@ std::vector getSamplesForClosureTest(std::string cut, float lumi, std:: std::vector vSample; - std::string tt = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTJets_Mu"+muID+"_El"+elID+".root"; + std::string tt = testDir.c_str()+"TTJets_Mu"+muID+"_El"+elID+".root"; TFile* ttfile = new TFile(tt.c_str()); Sample* ttSample = new Sample(vBkgNames.at(0),ttfile, vWeights.at(0),vXsec.at(0),cut,kRed+2); vSample.push_back(ttSample); - std::string nptt = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; + std::string nptt = testDir.c_str()+"NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; TFile* npttfile = new TFile(nptt.c_str()); Sample* npttSample = new Sample(vBkgNames.at(1),npttfile, vWeights.at(1),vXsec.at(1),cut,kBlue+2); vSample.push_back(npttSample); - std::string ttb = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTbar-powheg_Mu"+muID+"_El"+elID+".root"; + std::string ttb = testDir.c_str()+"TTbar-powheg_Mu"+muID+"_El"+elID+".root"; TFile* ttbfile = new TFile(ttb.c_str()); Sample* ttbSample = new Sample(vBkgNames.at(2),ttbfile, vWeights.at(2),vXsec.at(2),cut,kGreen+2); vSample.push_back(ttbSample); - std::string npttb = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/NonPromptTTbar-powheg_Mu"+muID+"_El"+elID+".root"; + std::string npttb = testDir.c_str()+"NonPromptTTbar-powheg_Mu"+muID+"_El"+elID+".root"; TFile* npttbfile = new TFile(npttb.c_str()); Sample* npttbSample = new Sample(vBkgNames.at(3),npttbfile, vWeights.at(3),vXsec.at(3),cut,kBlack); vSample.push_back(npttbSample); @@ -571,83 +585,83 @@ std::vector getBkgSampleVec(std::string cut, float lumi, std::string el //now make samples and add to vector std::vector vSample; - //TFile* ttfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTbar.root"); + //TFile* ttfile = new TFile(testDir.c_str()+"TTbar.root"); //Sample* ttSample = new Sample(vBkgNames.at(0),ttfile, vWeights.at(0),vXsec.at(0),cut,kRed+2); //vSample.push_back(ttSample); - std::string ttZfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTZ_Mu"+muID+"_El"+elID+".root"; + std::string ttZfilename = bgTestDir.c_str()+"TTZ_Mu"+muID+"_El"+elID+".root"; TFile* ttZfile = new TFile(ttZfilename.c_str()); Sample* ttZSample = new Sample(vBkgNames.at(0),ttZfile, vWeights.at(0),vXsec.at(0),cut,kRed); vSample.push_back(ttZSample); - std::string ttWfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTW_Mu"+muID+"_El"+elID+".root"; + std::string ttWfilename = bgTestDir.c_str()+"TTW_Mu"+muID+"_El"+elID+".root"; TFile* ttwfile = new TFile(ttWfilename.c_str()); Sample* ttwSample = new Sample(vBkgNames.at(1),ttwfile, vWeights.at(1),vXsec.at(1),cut,kYellow-2); vSample.push_back(ttwSample); - std::string tthfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTH_Mu"+muID+"_El"+elID+".root"; + std::string tthfilename = bgTestDir.c_str()+"TTH_Mu"+muID+"_El"+elID+".root"; TFile* tthfile = new TFile(tthfilename.c_str()); Sample* tthSample = new Sample(vBkgNames.at(2),tthfile, vWeights.at(2),vXsec.at(2),cut,kYellow); vSample.push_back(tthSample); - std::string ttttfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/TTTT_Mu"+muID+"_El"+elID+".root"; + std::string ttttfilename = bgTestDir.c_str()+"TTTT_Mu"+muID+"_El"+elID+".root"; TFile* ttttfile = new TFile(ttttfilename.c_str()); Sample* ttttSample = new Sample(vBkgNames.at(3),ttttfile, vWeights.at(3),vXsec.at(3),cut,kRed+2); vSample.push_back(ttttSample); - std::string wzfilename="/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/WZ_Mu"+muID+"_El"+elID+".root"; + std::string wzfilename=bgTestDir.c_str()+"WZ_Mu"+muID+"_El"+elID+".root"; TFile* wzfile = new TFile(wzfilename.c_str()); Sample* wzSample = new Sample(vBkgNames.at(4),wzfile, vWeights.at(4),vXsec.at(4),cut,kBlue-3); vSample.push_back(wzSample); - //TFile* wjfile = new TFile("/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/WJets_Mu"+muID+"_El"+elID+".root"); + //TFile* wjfile = new TFile(bgTestDir.c_str()+"WJets_Mu"+muID+"_El"+elID+".root"); //Sample* wjSample = new Sample(vBkgNames.at(4),wjfile, vWeights.at(4),vXsec.at(4),cut,kGreen+2); //vSample.push_back(wjSample); - /*std::string dyfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/DYJets_Mu"+muID+"_El"+elID+".root"; + /*std::string dyfilename = bgTestDir.c_str()+"DYJets_Mu"+muID+"_El"+elID+".root"; TFile* dyjfile = new TFile(dyfilename.c_str()); Sample* dyjSample = new Sample(vBkgNames.at(3),dyjfile, vWeights.at(3),vXsec.at(3),cut,kMagenta+2); vSample.push_back(dyjSample);*/ - std::string zzfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/ZZ_Mu"+muID+"_El"+elID+".root"; + std::string zzfilename = bgTestDir.c_str()+"ZZ_Mu"+muID+"_El"+elID+".root"; TFile* zzjfile = new TFile(zzfilename.c_str()); Sample* zzjSample = new Sample(vBkgNames.at(5),zzjfile, vWeights.at(5),vXsec.at(5),cut,kOrange+1); //std::cout<<"weight for ZZ is: "< getBkgSampleVec(std::string cut, float lumi, std::string el //********** Nonprompt *************** - /*std::string npTTfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; + /*std::string npTTfilename = bgTestDir.c_str()+"NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; TFile* npttfile = new TFile(nptTfilename.c_str()); Sample* npttSample = new Sample(vBkgNames.at(-1),npttfile,vWeights.at(-1),vXsec.at(01),cut,kBlue); vSample.push_back(npttSample)*/ - std::string npfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/NonPromptData_Mu"+muID+"_El"+elID+".root"; + std::string npfilename = bgTestDir.c_str()+"NonPromptData_Mu"+muID+"_El"+elID+".root"; TFile* npfile = new TFile(npfilename.c_str()); Sample* npSample = new Sample(vBkgNames.at(12),npfile,vWeights.at(13),vXsec.at(13),cut,kGray); vSample.push_back(npSample); //********ChargeMisID********** - std::string cmidfilename = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/ChargeMisID_Mu"+muID+"_El"+elID+".root"; + std::string cmidfilename = bgTestDir.c_str()+"ChargeMisID_Mu"+muID+"_El"+elID+".root"; TFile* cmidfile = new TFile(cmidfilename.c_str()); Sample* cmidSample = new Sample(vBkgNames.at(13),cmidfile,vWeights.at(13),vXsec.at(13),cut,kAzure+6); //force charge misID to start here since only at this point do we filter events vSample.push_back(cmidSample); From e9021fd38f04ecc85a604f8726fb7bdb3c8aca57 Mon Sep 17 00:00:00 2001 From: David Date: Thu, 4 Feb 2016 13:06:28 -0600 Subject: [PATCH 2/3] Add Dave's directory and ljmet files for single production. --- plugins/Macros.cc | 296 ++++++++++++++++++++++++++------------------ test/X53Analyzer.cc | 30 ++++- 2 files changed, 202 insertions(+), 124 deletions(-) diff --git a/plugins/Macros.cc b/plugins/Macros.cc index 6279088..4af2d21 100644 --- a/plugins/Macros.cc +++ b/plugins/Macros.cc @@ -15,13 +15,18 @@ #include "../interface/TreeReader.h" #include "TGraphAsymmErrors.h" +//global variable for test directory used for samples, etc. + std::string testDir = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_10_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/"; std::string bgTestDir = testDir; std::string x53String = "X53X53"; + + +//std::string testDir = "/uscms/home/dzou/CMSSW_7_3_0/src/X53ThirteenTeVAnalysisCode/test/"; /* -std::string testDir = "/uscms/home/dzou/CMSSW_7_3_0/src/X53ThirteenTeVAnalysisCode/test/"; //global variable for test directory used for samples, etc. -std::string testDir = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_10_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/"; -std::string x53String = "X53"; +std::string testDir = "/uscms_data/d3/dzou/using_git/T53/ljmet/CMSSW_7_4_15_patch1/src/dzoucms-fork/X53ThirteenTeVAnalysisCode/test"; +std::string bgTestDir = "/uscms_data/d3/clint/using_git/T53/ljmet/CMSSW_7_4_10_patch1/src/AnalysisCode/X53ThirteenTeVAnalysisCode/test/"; +std::string x53String = "X53T"; */ std::vector getVariableVec(){ @@ -49,7 +54,7 @@ std::vector getVariableVec(){ Variable* ak4jet1phi = new Variable("AK4Jet1Phi",20,-3.5,3.5,"Leading Jet #phi","N_{Events}"); vVar.push_back(ak4jet1phi); -< Variable* ak4jet2pt = new Variable("AK4Jet2Pt",15,0,600,"subLeading Jet p_{T} (GeV)","N_{Events}"); + Variable* ak4jet2pt = new Variable("AK4Jet2Pt",15,0,600,"subLeading Jet p_{T} (GeV)","N_{Events}"); vVar.push_back(ak4jet2pt); Variable* ak4jet2eta = new Variable("AK4Jet2Eta",25,-5,5,"subLeading Jet #eta","N_{Events}"); vVar.push_back(ak4jet2eta); @@ -195,25 +200,27 @@ std::vector getNPHistos(int nMu){ Sample* getDataSample(std::string cut,std::string elID, std::string muID){ - std::string filename = testDir.c_str()+muID+"_El"+elID+".root"; + std::string filename = testDir+muID+"_El"+elID+".root"; TFile* tfile = new TFile(filename.c_str()); Sample* dataSample = new Sample("Data",tfile,1,1,cut,kBlack,1); return dataSample; } std::vector getSigSampleVecForPlots(std::string cut, float lumi, std::string elID, std::string muID){ + //this method seems like it is no longer used and has been replaced by getSigSampleVecForTable and getInclusiveSigSampleVecForTable //br is OR of either side decaying to ssdl. BR(ssdl) for one side = BR(W->enu) OR BR(W->munu) OR BR(W->taunu) **2 where the square comes from AND requiring both to decay leptonically - float BRssdl= 2*( pow((.1063 + .1071 + .1138),2)); + float BRssdl= 2*( pow((.1063 + .1071 + .1138),2)); // for pair production + //float BRssdl= ( pow((.1063 + .1071 + .1138),2)); // for single production //make names vector //make x-sec vector - NEED TO FIX THESE VALUES************************************************ //make vector for number of events ran to get weight std::vector vNEvts; std::vector vSigNames; std::vector vXsec; - vSigNames.push_back("LH_X53X53_M-700"); vXsec.push_back(0.442*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-700"); vXsec.push_back(0.442*BRssdl);vNEvts.push_back(20000.); - //vSigNames.push_back("LH_X53X53_M-1000"); vXsec.push_back(0.0427*BRssdl);vNEvts.push_back(20000.); - //vSigNames.push_back("RH_X53X53_M-1000"); vXsec.push_back(0.0427*BRssdl);vNEvts.push_back(19400.); + vSigNames.push_back("LH_"+x53String+"_M-700"); vXsec.push_back(0.442*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-700"); vXsec.push_back(0.442*BRssdl);vNEvts.push_back(20000.); + //vSigNames.push_back("LH_"+x53String+"_M-1000"); vXsec.push_back(0.0427*BRssdl);vNEvts.push_back(20000.); + //vSigNames.push_back("RH_"+x53String+"_M-1000"); vXsec.push_back(0.0427*BRssdl);vNEvts.push_back(19400.); //FOR RUNNING ON INCLUSIVE DON'T APPLY BR TO SSDL!! //vXsec.push_back(0.442);vXsec.push_back(0.442);vXsec.push_back(0.0427);vXsec.push_back(0.0427);vXsec.push_back(0.00618);vXsec.push_back(0.00618); @@ -229,29 +236,29 @@ std::vector getSigSampleVecForPlots(std::string cut, float lumi, std::s std::vector vSigSamples; - std::string lh700name =testDir.c_str()+"X53X53m700LH_Mu"+muID+"_El"+elID+".root"; + std::string lh700name =testDir+x53String+"m700LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Lfile = new TFile(lh700name.c_str()); Sample* x53x53m700L = new Sample(vSigNames.at(0),x53x53700Lfile,vWeights.at(0),vXsec.at(0),cut,kGreen,2); vSigSamples.push_back(x53x53m700L); - std::string rh700name = testDir.c_str()+"X53X53m700RH_Mu"+muID+"_El"+elID+".root"; + std::string rh700name = testDir+x53String+"m700RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Rfile = new TFile(rh700name.c_str()); Sample* x53x53m700R = new Sample(vSigNames.at(1),x53x53700Rfile,vWeights.at(1),vXsec.at(1),cut,kGreen,1); vSigSamples.push_back(x53x53m700R); - /* std::string lh1000name = testDir.c_str()+"X53X53m1000LH_Mu"+muID+"_El"+elID+".root"; + /* std::string lh1000name = testDir+x53String+"m1000LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Lfile = new TFile(lh1000name.c_str()); Sample* x53x53m1000L = new Sample(vSigNames.at(2),x53x531000Lfile,vWeights.at(2),vXsec.at(2),cut,kCyan,2); vSigSamples.push_back(x53x53m1000L); - std::string rh1000name = testDir.c_str()+"X53X53m1000RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1000name = testDir+x53String+"m1000RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Rfile = new TFile(rh1000name.c_str()); Sample* x53x53m1000R = new Sample(vSigNames.at(3),x53x531000Rfile,vWeights.at(3),vXsec.at(3),cut,kCyan,1); vSigSamples.push_back(x53x53m1000R);*/ - /* TFile* x53x531300Lfile = new TFile(testDir.c_str()+"X53X53m1300LH.root"); + /* TFile* x53x531300Lfile = new TFile(testDir+x53String+"m1300LH.root"); Sample* x53x53m1300L = new Sample(vSigNames.at(4),x53x531300Lfile,vWeights.at(4),vXsec.at(4),cut,kBlue,2); - TFile* x53x531300Rfile = new TFile(testDir.c_str()+"X53X53m1300RH.root"); + TFile* x53x531300Rfile = new TFile(testDir+x53String+"m1300RH.root"); Sample* x53x53m1300R = new Sample(vSigNames.at(5),x53x531300Rfile,vWeights.at(5),vXsec.at(5),cut,kBlue,1); */ return vSigSamples; @@ -263,28 +270,51 @@ std::vector getSigSampleVecForTable(std::string cut, float lumi, std::s std::vector vSigNames; std::vector vXsec; //br is OR of either side decaying to ssdl. BR(ssdl) for one side = BR(W->enu) OR BR(W->munu) **2 where the square comes from AND requiring both to decay leptonically - float BRssdl= 2*( pow((.1063 + .1071 + .1138),2)); std::vector vNEvts; - vSigNames.push_back("LH_X53X53_M-700"); vXsec.push_back(0.442*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-700"); vXsec.push_back(0.442*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("LH_X53X53_M-800"); vXsec.push_back(.190*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-800"); vXsec.push_back(.190*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("LH_X53X53_M-900"); vXsec.push_back(.0877*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-900"); vXsec.push_back(.0877*BRssdl); vNEvts.push_back(19800.); - vSigNames.push_back("LH_X53X53_M-1000"); vXsec.push_back(0.0427*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-1000"); vXsec.push_back(0.0427*BRssdl); vNEvts.push_back(19400.); - vSigNames.push_back("LH_X53X53_M-1100"); vXsec.push_back(0.0217*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-1100"); vXsec.push_back(0.0217*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("LH_X53X53_M-1300"); vXsec.push_back(0.00618*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-1300"); vXsec.push_back(0.00618*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("LH_X53X53_M-1400"); vXsec.push_back(0.00342*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-1400"); vXsec.push_back(0.00342*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("LH_X53X53_M-1500"); vXsec.push_back(0.00193*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-1500"); vXsec.push_back(0.00193*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("LH_X53X53_M-1600"); vXsec.push_back(0.00111*BRssdl); vNEvts.push_back(20000.); - vSigNames.push_back("RH_X53X53_M-1600"); vXsec.push_back(0.00111*BRssdl); vNEvts.push_back(20000.); - - + if ( x53String.compare("X53X53") == 0 ) { // pair production + float BRssdl= 2*( pow((.1063 + .1071 + .1138),2)); + vSigNames.push_back("LH_"+x53String+"_M-700"); vXsec.push_back(0.442*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-700"); vXsec.push_back(0.442*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("LH_"+x53String+"_M-800"); vXsec.push_back(.190*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-800"); vXsec.push_back(.190*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("LH_"+x53String+"_M-900"); vXsec.push_back(.0877*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-900"); vXsec.push_back(.0877*BRssdl); vNEvts.push_back(19800.); + vSigNames.push_back("LH_"+x53String+"_M-1000"); vXsec.push_back(0.0427*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-1000"); vXsec.push_back(0.0427*BRssdl); vNEvts.push_back(19400.); + vSigNames.push_back("LH_"+x53String+"_M-1100"); vXsec.push_back(0.0217*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-1100"); vXsec.push_back(0.0217*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("LH_"+x53String+"_M-1300"); vXsec.push_back(0.00618*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-1300"); vXsec.push_back(0.00618*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("LH_"+x53String+"_M-1400"); vXsec.push_back(0.00342*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-1400"); vXsec.push_back(0.00342*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("LH_"+x53String+"_M-1500"); vXsec.push_back(0.00193*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-1500"); vXsec.push_back(0.00193*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("LH_"+x53String+"_M-1600"); vXsec.push_back(0.00111*BRssdl); vNEvts.push_back(20000.); + vSigNames.push_back("RH_"+x53String+"_M-1600"); vXsec.push_back(0.00111*BRssdl); vNEvts.push_back(20000.); + } + else { // single production + float BRssdl= ( pow((.1063 + .1071 + .1138),2)); + vSigNames.push_back("LH_"+x53String+"_M-700"); vXsec.push_back(0.745*BRssdl); vNEvts.push_back(49900.); + vSigNames.push_back("RH_"+x53String+"_M-700"); vXsec.push_back(0.745*BRssdl); vNEvts.push_back(49800.); + vSigNames.push_back("LH_"+x53String+"_M-800"); vXsec.push_back(.532*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("RH_"+x53String+"_M-800"); vXsec.push_back(.532*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("LH_"+x53String+"_M-900"); vXsec.push_back(.388*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("RH_"+x53String+"_M-900"); vXsec.push_back(.388*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("LH_"+x53String+"_M-1000"); vXsec.push_back(0.285*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("RH_"+x53String+"_M-1000"); vXsec.push_back(0.285*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("LH_"+x53String+"_M-1100"); vXsec.push_back(0.212*BRssdl); vNEvts.push_back(299600.);// + vSigNames.push_back("RH_"+x53String+"_M-1100"); vXsec.push_back(0.212*BRssdl); vNEvts.push_back(299800.);// + vSigNames.push_back("LH_"+x53String+"_M-1200"); vXsec.push_back(0.159*BRssdl); vNEvts.push_back(295400.);// + vSigNames.push_back("RH_"+x53String+"_M-1200"); vXsec.push_back(0.159*BRssdl); vNEvts.push_back(300000.);// + vSigNames.push_back("LH_"+x53String+"_M-1300"); vXsec.push_back(0.120*BRssdl); vNEvts.push_back(50000.); + vSigNames.push_back("RH_"+x53String+"_M-1300"); vXsec.push_back(0.120*BRssdl); vNEvts.push_back(49400.); + vSigNames.push_back("LH_"+x53String+"_M-1400"); vXsec.push_back(0.0917*BRssdl); vNEvts.push_back(300000.);// + vSigNames.push_back("RH_"+x53String+"_M-1400"); vXsec.push_back(0.0917*BRssdl); vNEvts.push_back(298400.);// + vSigNames.push_back("LH_"+x53String+"_M-1500"); vXsec.push_back(0.0706*BRssdl); vNEvts.push_back(298400.);// + vSigNames.push_back("RH_"+x53String+"_M-1500"); vXsec.push_back(0.0706*BRssdl); vNEvts.push_back(300000.);// + vSigNames.push_back("LH_"+x53String+"_M-1600"); vXsec.push_back(0.0541*BRssdl); vNEvts.push_back(300000.);// + vSigNames.push_back("RH_"+x53String+"_M-1600"); vXsec.push_back(0.0541*BRssdl); vNEvts.push_back(300000.);// + } //vector to hold weights std::vector vWeights; @@ -293,82 +323,82 @@ std::vector getSigSampleVecForTable(std::string cut, float lumi, std::s } std::vector vSigSamples; - std::string lh700 = testDir.c_str()+"X53X53m700LH_Mu"+muID+"_El"+elID+".root"; + std::string lh700 = testDir+x53String+"m700LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Lfile = new TFile(lh700.c_str()); Sample* x53x53m700L = new Sample(vSigNames.at(0),x53x53700Lfile,vWeights.at(0),vXsec.at(0),cut,kGreen,2); vSigSamples.push_back(x53x53m700L); - std::string rh700 = testDir.c_str()+"X53X53m700RH_Mu"+muID+"_El"+elID+".root"; + std::string rh700 = testDir+x53String+"m700RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Rfile = new TFile(rh700.c_str()); Sample* x53x53m700R = new Sample(vSigNames.at(1),x53x53700Rfile,vWeights.at(1),vXsec.at(1),cut,kGreen,1); vSigSamples.push_back(x53x53m700R); - std::string lh800 = testDir.c_str()+"X53X53m800LH_Mu"+muID+"_El"+elID+".root"; + std::string lh800 = testDir+x53String+"m800LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Lfile = new TFile(lh800.c_str()); Sample* x53x53m800L = new Sample(vSigNames.at(2),x53x53800Lfile,vWeights.at(2),vXsec.at(2),cut,kBlue,2); vSigSamples.push_back(x53x53m800L); - std::string rh800 = testDir.c_str()+"X53X53m800RH_Mu"+muID+"_El"+elID+".root"; + std::string rh800 = testDir+x53String+"m800RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Rfile = new TFile(rh800.c_str()); Sample* x53x53m800R = new Sample(vSigNames.at(3),x53x53800Rfile,vWeights.at(3),vXsec.at(3),cut,kBlue,1); vSigSamples.push_back(x53x53m800R); - std::string lh900 = testDir.c_str()+"X53X53m900LH_Mu"+muID+"_El"+elID+".root"; + std::string lh900 = testDir+x53String+"m900LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Lfile = new TFile(lh900.c_str()); Sample* x53x53m900L = new Sample(vSigNames.at(4),x53x53900Lfile,vWeights.at(4),vXsec.at(4),cut,kGreen,2); - //vSigSamples.push_back(x53x53m900L); - std::string rh900 = testDir.c_str()+"X53X53m900RH_Mu"+muID+"_El"+elID+".root"; + vSigSamples.push_back(x53x53m900L); + std::string rh900 = testDir+x53String+"m900RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Rfile = new TFile(rh900.c_str()); Sample* x53x53m900R = new Sample(vSigNames.at(5),x53x53900Rfile,vWeights.at(5),vXsec.at(5),cut,kGreen,1); vSigSamples.push_back(x53x53m900R); - std::string lh1000 = testDir.c_str()+"X53X53m1000LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1000 = testDir+x53String+"m1000LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Lfile = new TFile(lh1000.c_str()); Sample* x53x53m1000L = new Sample(vSigNames.at(6),x53x531000Lfile,vWeights.at(6),vXsec.at(6),cut,kCyan,2); - //vSigSamples.push_back(x53x53m1000L); - std::string rh1000 = testDir.c_str()+"X53X53m1000RH_Mu"+muID+"_El"+elID+".root"; + vSigSamples.push_back(x53x53m1000L); + std::string rh1000 = testDir+x53String+"m1000RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Rfile = new TFile(rh1000.c_str()); Sample* x53x53m1000R = new Sample(vSigNames.at(7),x53x531000Rfile,vWeights.at(7),vXsec.at(7),cut,kCyan,1); vSigSamples.push_back(x53x53m1000R); - std::string lh1100 = testDir.c_str()+"X53X53m1100LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1100 = testDir+x53String+"m1100LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Lfile = new TFile(lh1100.c_str()); Sample* x53x53m1100L = new Sample(vSigNames.at(8),x53x531100Lfile,vWeights.at(8),vXsec.at(8),cut,kBlue,2); vSigSamples.push_back(x53x53m1100L); - std::string rh1100 = testDir.c_str()+"X53X53m1100RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1100 = testDir+x53String+"m1100RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Rfile = new TFile(rh1100.c_str()); Sample* x53x53m1100R = new Sample(vSigNames.at(9),x53x531100Rfile,vWeights.at(9),vXsec.at(9),cut,kBlue,1); vSigSamples.push_back(x53x53m1100R); - std::string lh1300 = testDir.c_str()+"X53X53m1300LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1300 = testDir+x53String+"m1300LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Lfile = new TFile(lh1300.c_str()); Sample* x53x53m1300L = new Sample(vSigNames.at(10),x53x531300Lfile,vWeights.at(10),vXsec.at(10),cut,kBlue,2); vSigSamples.push_back(x53x53m1300L); - std::string rh1300 = testDir.c_str()+"X53X53m1300RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1300 = testDir+x53String+"m1300RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Rfile = new TFile(rh1300.c_str()); Sample* x53x53m1300R = new Sample(vSigNames.at(11),x53x531300Rfile,vWeights.at(11),vXsec.at(11),cut,kBlue,1); - //vSigSamples.push_back(x53x53m1300R); + vSigSamples.push_back(x53x53m1300R); - std::string lh1400 = testDir.c_str()+"X53X53m1400LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1400 = testDir+x53String+"m1400LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Lfile = new TFile(lh1400.c_str()); Sample* x53x53m1400L = new Sample(vSigNames.at(12),x53x531400Lfile,vWeights.at(12),vXsec.at(12),cut,kBlue,2); vSigSamples.push_back(x53x53m1400L); - std::string rh1400 = testDir.c_str()+"X53X53m1400RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1400 = testDir+x53String+"m1400RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Rfile = new TFile(rh1400.c_str()); Sample* x53x53m1400R = new Sample(vSigNames.at(13),x53x531400Rfile,vWeights.at(13),vXsec.at(13),cut,kBlue,1); vSigSamples.push_back(x53x53m1400R); - std::string lh1500 = testDir.c_str()+"X53X53m1500LH_Mu"+muID+"_El"+elID+".root"; + std::string lh1500 = testDir+x53String+"m1500LH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Lfile = new TFile(lh1500.c_str()); Sample* x53x53m1500L = new Sample(vSigNames.at(14),x53x531500Lfile,vWeights.at(14),vXsec.at(14),cut,kBlue,2); vSigSamples.push_back(x53x53m1500L); - std::string rh1500 = testDir.c_str()+"X53X53m1500RH_Mu"+muID+"_El"+elID+".root"; + std::string rh1500 = testDir+x53String+"m1500RH_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Rfile = new TFile(rh1500.c_str()); Sample* x53x53m1500R = new Sample(vSigNames.at(15),x53x531500Rfile,vWeights.at(15),vXsec.at(15),cut,kBlue,1); vSigSamples.push_back(x53x53m1500R); - /* TFile* x53x531600Lfile = new TFile(testDir.c_str()+"X53X53m1600LH_Mu"+muID+"_El"+elID+".root"); + /* TFile* x53x531600Lfile = new TFile(testDir+x53String+"m1600LH_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600L = new Sample(vSigNames.at(16),x53x531600Lfile,vWeights.at(16),vXsec.at(16),cut,kBlue,2); vSigSamples.push_back(x53x53m1600L); - TFile* x53x531600Rfile = new TFile(testDir.c_str()+"X53X53m1600RH_Mu"+muID+"_El"+elID+".root"); + TFile* x53x531600Rfile = new TFile(testDir+x53String+"m1600RH_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600R = new Sample(vSigNames.at(17),x53x531600Rfile,vWeights.at(17),vXsec.at(17),cut,kBlue,1); vSigSamples.push_back(x53x53m1600R);*/ @@ -385,28 +415,50 @@ std::vector getInclusiveSigSampleVecForTable(std::string cut, float lum //br is OR of either side decaying to ssdl. BR(ssdl) for one side = BR(W->enu) OR BR(W->munu) **2 where the square comes from AND requiring both to decay leptonically //float BRssdl= 2*( pow((.1063 + .1071 + .1138),2)); INCLUSIVE SO NO BRANCHING RATIO std::vector vNEvts; - vSigNames.push_back("LH_X53X53_M-700"); vXsec.push_back(0.442); vNEvts.push_back(300000.); - vSigNames.push_back("RH_X53X53_M-700"); vXsec.push_back(0.442); vNEvts.push_back(297400.); - vSigNames.push_back("LH_X53X53_M-800"); vXsec.push_back(.190); vNEvts.push_back(295600.); - vSigNames.push_back("RH_X53X53_M-800"); vXsec.push_back(.190); vNEvts.push_back(299600.); - vSigNames.push_back("LH_X53X53_M-900"); vXsec.push_back(.0877); vNEvts.push_back(300000.); - vSigNames.push_back("RH_X53X53_M-900"); vXsec.push_back(.0877); vNEvts.push_back(299800.); - vSigNames.push_back("LH_X53X53_M-1000"); vXsec.push_back(0.0427); vNEvts.push_back(293600.); - vSigNames.push_back("RH_X53X53_M-1000"); vXsec.push_back(0.0427); vNEvts.push_back(299000.); - vSigNames.push_back("LH_X53X53_M-1100"); vXsec.push_back(0.0217); vNEvts.push_back(299600.); - vSigNames.push_back("RH_X53X53_M-1100"); vXsec.push_back(0.0217); vNEvts.push_back(299800.); - vSigNames.push_back("LH_X53X53_M-1200"); vXsec.push_back(0.0114); vNEvts.push_back(295400.); - vSigNames.push_back("RH_X53X53_M-1200"); vXsec.push_back(0.0114); vNEvts.push_back(300000.); - vSigNames.push_back("LH_X53X53_M-1300"); vXsec.push_back(0.00618); vNEvts.push_back(300000.); - vSigNames.push_back("RH_X53X53_M-1300"); vXsec.push_back(0.00618); vNEvts.push_back(293600.); - vSigNames.push_back("LH_X53X53_M-1400"); vXsec.push_back(0.00342); vNEvts.push_back(300000.); - vSigNames.push_back("RH_X53X53_M-1400"); vXsec.push_back(0.00342); vNEvts.push_back(298400.); - vSigNames.push_back("LH_X53X53_M-1500"); vXsec.push_back(0.00193); vNEvts.push_back(298400.); - vSigNames.push_back("RH_X53X53_M-1500"); vXsec.push_back(0.00193); vNEvts.push_back(300000.); - vSigNames.push_back("LH_X53X53_M-1600"); vXsec.push_back(0.00111); vNEvts.push_back(300000.); - vSigNames.push_back("RH_X53X53_M-1600"); vXsec.push_back(0.00111); vNEvts.push_back(300000.); - - + if ( x53String.compare("X53X53") == 0 ) { //for pair production + vSigNames.push_back("LH_"+x53String+"_M-700"); vXsec.push_back(0.442); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-700"); vXsec.push_back(0.442); vNEvts.push_back(297400.); + vSigNames.push_back("LH_"+x53String+"_M-800"); vXsec.push_back(.190); vNEvts.push_back(295600.); + vSigNames.push_back("RH_"+x53String+"_M-800"); vXsec.push_back(.190); vNEvts.push_back(299600.); + vSigNames.push_back("LH_"+x53String+"_M-900"); vXsec.push_back(.0877); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-900"); vXsec.push_back(.0877); vNEvts.push_back(299800.); + vSigNames.push_back("LH_"+x53String+"_M-1000"); vXsec.push_back(0.0427); vNEvts.push_back(293600.); + vSigNames.push_back("RH_"+x53String+"_M-1000"); vXsec.push_back(0.0427); vNEvts.push_back(299000.); + vSigNames.push_back("LH_"+x53String+"_M-1100"); vXsec.push_back(0.0217); vNEvts.push_back(299600.); + vSigNames.push_back("RH_"+x53String+"_M-1100"); vXsec.push_back(0.0217); vNEvts.push_back(299800.); + vSigNames.push_back("LH_"+x53String+"_M-1200"); vXsec.push_back(0.0114); vNEvts.push_back(295400.); + vSigNames.push_back("RH_"+x53String+"_M-1200"); vXsec.push_back(0.0114); vNEvts.push_back(300000.); + vSigNames.push_back("LH_"+x53String+"_M-1300"); vXsec.push_back(0.00618); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-1300"); vXsec.push_back(0.00618); vNEvts.push_back(293600.); + vSigNames.push_back("LH_"+x53String+"_M-1400"); vXsec.push_back(0.00342); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-1400"); vXsec.push_back(0.00342); vNEvts.push_back(298400.); + vSigNames.push_back("LH_"+x53String+"_M-1500"); vXsec.push_back(0.00193); vNEvts.push_back(298400.); + vSigNames.push_back("RH_"+x53String+"_M-1500"); vXsec.push_back(0.00193); vNEvts.push_back(300000.); + vSigNames.push_back("LH_"+x53String+"_M-1600"); vXsec.push_back(0.00111); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-1600"); vXsec.push_back(0.00111); vNEvts.push_back(300000.); + } + else { + vSigNames.push_back("LH_"+x53String+"_M-700"); vXsec.push_back(0.745); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-700"); vXsec.push_back(0.745); vNEvts.push_back(297400.); + vSigNames.push_back("LH_"+x53String+"_M-800"); vXsec.push_back(.532); vNEvts.push_back(295600.); + vSigNames.push_back("RH_"+x53String+"_M-800"); vXsec.push_back(.532); vNEvts.push_back(299600.); + vSigNames.push_back("LH_"+x53String+"_M-900"); vXsec.push_back(.388); vNEvts.push_back(50000.); + vSigNames.push_back("RH_"+x53String+"_M-900"); vXsec.push_back(.388); vNEvts.push_back(299800.); + vSigNames.push_back("LH_"+x53String+"_M-1000"); vXsec.push_back(0.285); vNEvts.push_back(293600.); + vSigNames.push_back("RH_"+x53String+"_M-1000"); vXsec.push_back(0.285); vNEvts.push_back(299000.); + vSigNames.push_back("LH_"+x53String+"_M-1100"); vXsec.push_back(0.212); vNEvts.push_back(299600.); + vSigNames.push_back("RH_"+x53String+"_M-1100"); vXsec.push_back(0.212); vNEvts.push_back(299800.); + vSigNames.push_back("LH_"+x53String+"_M-1200"); vXsec.push_back(0.159); vNEvts.push_back(295400.); + vSigNames.push_back("RH_"+x53String+"_M-1200"); vXsec.push_back(0.159); vNEvts.push_back(300000.); + vSigNames.push_back("LH_"+x53String+"_M-1300"); vXsec.push_back(0.120); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-1300"); vXsec.push_back(0.120); vNEvts.push_back(293600.); + vSigNames.push_back("LH_"+x53String+"_M-1400"); vXsec.push_back(0.0917); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-1400"); vXsec.push_back(0.0917); vNEvts.push_back(298400.); + vSigNames.push_back("LH_"+x53String+"_M-1500"); vXsec.push_back(0.0706); vNEvts.push_back(298400.); + vSigNames.push_back("RH_"+x53String+"_M-1500"); vXsec.push_back(0.0706); vNEvts.push_back(300000.); + vSigNames.push_back("LH_"+x53String+"_M-1600"); vXsec.push_back(0.0541); vNEvts.push_back(300000.); + vSigNames.push_back("RH_"+x53String+"_M-1600"); vXsec.push_back(0.0541); vNEvts.push_back(300000.); + } //vector to hold weights std::vector vWeights; @@ -415,91 +467,91 @@ std::vector getInclusiveSigSampleVecForTable(std::string cut, float lum } std::vector vSigSamples; - std::string lh700 = testDir.c_str()+"X53X53m700LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh700 = testDir+x53String+"m700LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Lfile = new TFile(lh700.c_str()); Sample* x53x53m700L = new Sample(vSigNames.at(0),x53x53700Lfile,vWeights.at(0),vXsec.at(0),cut,kGreen,2); vSigSamples.push_back(x53x53m700L); - std::string rh700 = testDir.c_str()+"X53X53m700RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh700 = testDir+x53String+"m700RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53700Rfile = new TFile(rh700.c_str()); Sample* x53x53m700R = new Sample(vSigNames.at(1),x53x53700Rfile,vWeights.at(1),vXsec.at(1),cut,kGreen,1); vSigSamples.push_back(x53x53m700R); - std::string lh800 = testDir.c_str()+"X53X53m800LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh800 = testDir+x53String+"m800LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Lfile = new TFile(lh800.c_str()); Sample* x53x53m800L = new Sample(vSigNames.at(2),x53x53800Lfile,vWeights.at(2),vXsec.at(2),cut,kBlue,2); vSigSamples.push_back(x53x53m800L); - std::string rh800 = testDir.c_str()+"X53X53m800RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh800 = testDir+x53String+"m800RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53800Rfile = new TFile(rh800.c_str()); Sample* x53x53m800R = new Sample(vSigNames.at(3),x53x53800Rfile,vWeights.at(3),vXsec.at(3),cut,kBlue,1); vSigSamples.push_back(x53x53m800R); - std::string lh900 = testDir.c_str()+"X53X53m900LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh900 = testDir+x53String+"m900LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Lfile = new TFile(lh900.c_str()); Sample* x53x53m900L = new Sample(vSigNames.at(4),x53x53900Lfile,vWeights.at(4),vXsec.at(4),cut,kGreen,2); vSigSamples.push_back(x53x53m900L); - std::string rh900 = testDir.c_str()+"X53X53m900RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh900 = testDir+x53String+"m900RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x53900Rfile = new TFile(rh900.c_str()); Sample* x53x53m900R = new Sample(vSigNames.at(5),x53x53900Rfile,vWeights.at(5),vXsec.at(5),cut,kGreen,1); vSigSamples.push_back(x53x53m900R); - std::string lh1000 = testDir.c_str()+"X53X53m1000LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1000 = testDir+x53String+"m1000LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Lfile = new TFile(lh1000.c_str()); Sample* x53x53m1000L = new Sample(vSigNames.at(6),x53x531000Lfile,vWeights.at(6),vXsec.at(6),cut,kCyan,2); vSigSamples.push_back(x53x53m1000L); - std::string rh1000 = testDir.c_str()+"X53X53m1000RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1000 = testDir+x53String+"m1000RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531000Rfile = new TFile(rh1000.c_str()); Sample* x53x53m1000R = new Sample(vSigNames.at(7),x53x531000Rfile,vWeights.at(7),vXsec.at(7),cut,kCyan,1); vSigSamples.push_back(x53x53m1000R); - std::string lh1100 = testDir.c_str()+"X53X53m1100LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1100 = testDir+x53String+"m1100LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Lfile = new TFile(lh1100.c_str()); Sample* x53x53m1100L = new Sample(vSigNames.at(8),x53x531100Lfile,vWeights.at(8),vXsec.at(8),cut,kBlue,2); vSigSamples.push_back(x53x53m1100L); - std::string rh1100 = testDir.c_str()+"X53X53m1100RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1100 = testDir+x53String+"m1100RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531100Rfile = new TFile(rh1100.c_str()); Sample* x53x53m1100R = new Sample(vSigNames.at(9),x53x531100Rfile,vWeights.at(9),vXsec.at(9),cut,kBlue,1); vSigSamples.push_back(x53x53m1100R); - std::string lh1200 = testDir.c_str()+"X53X53m1200LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1200 = testDir+x53String+"m1200LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531200Lfile = new TFile(lh1200.c_str()); Sample* x53x53m1200L = new Sample(vSigNames.at(10),x53x531200Lfile,vWeights.at(10),vXsec.at(10),cut,kBlue,2); vSigSamples.push_back(x53x53m1200L); - std::string rh1200 = testDir.c_str()+"X53X53m1200RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1200 = testDir+x53String+"m1200RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531200Rfile = new TFile(rh1200.c_str()); Sample* x53x53m1200R = new Sample(vSigNames.at(11),x53x531200Rfile,vWeights.at(11),vXsec.at(11),cut,kBlue,1); vSigSamples.push_back(x53x53m1200R); - std::string lh1300 = testDir.c_str()+"X53X53m1300LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1300 = testDir+x53String+"m1300LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Lfile = new TFile(lh1300.c_str()); Sample* x53x53m1300L = new Sample(vSigNames.at(12),x53x531300Lfile,vWeights.at(12),vXsec.at(12),cut,kBlue,2); vSigSamples.push_back(x53x53m1300L); - std::string rh1300 = testDir.c_str()+"X53X53m1300RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1300 = testDir+x53String+"m1300RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531300Rfile = new TFile(rh1300.c_str()); Sample* x53x53m1300R = new Sample(vSigNames.at(13),x53x531300Rfile,vWeights.at(13),vXsec.at(13),cut,kBlue,1); vSigSamples.push_back(x53x53m1300R); - std::string lh1400 = testDir.c_str()+"X53X53m1400LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1400 = testDir+x53String+"m1400LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Lfile = new TFile(lh1400.c_str()); Sample* x53x53m1400L = new Sample(vSigNames.at(14),x53x531400Lfile,vWeights.at(14),vXsec.at(14),cut,kBlue,2); vSigSamples.push_back(x53x53m1400L); - std::string rh1400 = testDir.c_str()+"X53X53m1400RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1400 = testDir+x53String+"m1400RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531400Rfile = new TFile(rh1400.c_str()); Sample* x53x53m1400R = new Sample(vSigNames.at(15),x53x531400Rfile,vWeights.at(15),vXsec.at(15),cut,kBlue,1); vSigSamples.push_back(x53x53m1400R); - std::string lh1500 = testDir.c_str()+"X53X53m1500LH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string lh1500 = testDir+x53String+"m1500LH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Lfile = new TFile(lh1500.c_str()); Sample* x53x53m1500L = new Sample(vSigNames.at(16),x53x531500Lfile,vWeights.at(16),vXsec.at(16),cut,kBlue,2); vSigSamples.push_back(x53x53m1500L); - std::string rh1500 = testDir.c_str()+"X53X53m1500RH_Inc_Mu"+muID+"_El"+elID+".root"; + std::string rh1500 = testDir+x53String+"m1500RH_Inc_Mu"+muID+"_El"+elID+".root"; TFile* x53x531500Rfile = new TFile(rh1500.c_str()); Sample* x53x53m1500R = new Sample(vSigNames.at(17),x53x531500Rfile,vWeights.at(17),vXsec.at(17),cut,kBlue,1); vSigSamples.push_back(x53x53m1500R); - /* TFile* x53x531600Lfile = new TFile(testDir.c_str()+"X53X53m1600LH_Inc_Mu"+muID+"_El"+elID+".root"); + /* TFile* x53x531600Lfile = new TFile(testDir+x53String+"m1600LH_Inc_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600L = new Sample(vSigNames.at(16),x53x531600Lfile,vWeights.at(16),vXsec.at(16),cut,kBlue,2); vSigSamples.push_back(x53x53m1600L); - TFile* x53x531600Rfile = new TFile(testDir.c_str()+"X53X53m1600RH_Inc_Mu"+muID+"_El"+elID+".root"); + TFile* x53x531600Rfile = new TFile(testDir+x53String+"m1600RH_Inc_Mu"+muID+"_El"+elID+".root"); Sample* x53x53m1600R = new Sample(vSigNames.at(17),x53x531600Rfile,vWeights.at(17),vXsec.at(17),cut,kBlue,1); vSigSamples.push_back(x53x53m1600R);*/ @@ -526,23 +578,23 @@ std::vector getSamplesForClosureTest(std::string cut, float lumi, std:: std::vector vSample; - std::string tt = testDir.c_str()+"TTJets_Mu"+muID+"_El"+elID+".root"; + std::string tt = testDir+"TTJets_Mu"+muID+"_El"+elID+".root"; TFile* ttfile = new TFile(tt.c_str()); Sample* ttSample = new Sample(vBkgNames.at(0),ttfile, vWeights.at(0),vXsec.at(0),cut,kRed+2); vSample.push_back(ttSample); - std::string nptt = testDir.c_str()+"NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; + std::string nptt = testDir+"NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; TFile* npttfile = new TFile(nptt.c_str()); Sample* npttSample = new Sample(vBkgNames.at(1),npttfile, vWeights.at(1),vXsec.at(1),cut,kBlue+2); vSample.push_back(npttSample); - std::string ttb = testDir.c_str()+"TTbar-powheg_Mu"+muID+"_El"+elID+".root"; + std::string ttb = testDir+"TTbar-powheg_Mu"+muID+"_El"+elID+".root"; TFile* ttbfile = new TFile(ttb.c_str()); Sample* ttbSample = new Sample(vBkgNames.at(2),ttbfile, vWeights.at(2),vXsec.at(2),cut,kGreen+2); vSample.push_back(ttbSample); - std::string npttb = testDir.c_str()+"NonPromptTTbar-powheg_Mu"+muID+"_El"+elID+".root"; + std::string npttb = testDir+"NonPromptTTbar-powheg_Mu"+muID+"_El"+elID+".root"; TFile* npttbfile = new TFile(npttb.c_str()); Sample* npttbSample = new Sample(vBkgNames.at(3),npttbfile, vWeights.at(3),vXsec.at(3),cut,kBlack); vSample.push_back(npttbSample); @@ -585,83 +637,83 @@ std::vector getBkgSampleVec(std::string cut, float lumi, std::string el //now make samples and add to vector std::vector vSample; - //TFile* ttfile = new TFile(testDir.c_str()+"TTbar.root"); + //TFile* ttfile = new TFile(testDir+"TTbar.root"); //Sample* ttSample = new Sample(vBkgNames.at(0),ttfile, vWeights.at(0),vXsec.at(0),cut,kRed+2); //vSample.push_back(ttSample); - std::string ttZfilename = bgTestDir.c_str()+"TTZ_Mu"+muID+"_El"+elID+".root"; + std::string ttZfilename = bgTestDir+"TTZ_Mu"+muID+"_El"+elID+".root"; TFile* ttZfile = new TFile(ttZfilename.c_str()); Sample* ttZSample = new Sample(vBkgNames.at(0),ttZfile, vWeights.at(0),vXsec.at(0),cut,kRed); vSample.push_back(ttZSample); - std::string ttWfilename = bgTestDir.c_str()+"TTW_Mu"+muID+"_El"+elID+".root"; + std::string ttWfilename = bgTestDir+"TTW_Mu"+muID+"_El"+elID+".root"; TFile* ttwfile = new TFile(ttWfilename.c_str()); Sample* ttwSample = new Sample(vBkgNames.at(1),ttwfile, vWeights.at(1),vXsec.at(1),cut,kYellow-2); vSample.push_back(ttwSample); - std::string tthfilename = bgTestDir.c_str()+"TTH_Mu"+muID+"_El"+elID+".root"; + std::string tthfilename = bgTestDir+"TTH_Mu"+muID+"_El"+elID+".root"; TFile* tthfile = new TFile(tthfilename.c_str()); Sample* tthSample = new Sample(vBkgNames.at(2),tthfile, vWeights.at(2),vXsec.at(2),cut,kYellow); vSample.push_back(tthSample); - std::string ttttfilename = bgTestDir.c_str()+"TTTT_Mu"+muID+"_El"+elID+".root"; + std::string ttttfilename = bgTestDir+"TTTT_Mu"+muID+"_El"+elID+".root"; TFile* ttttfile = new TFile(ttttfilename.c_str()); Sample* ttttSample = new Sample(vBkgNames.at(3),ttttfile, vWeights.at(3),vXsec.at(3),cut,kRed+2); vSample.push_back(ttttSample); - std::string wzfilename=bgTestDir.c_str()+"WZ_Mu"+muID+"_El"+elID+".root"; + std::string wzfilename=bgTestDir+"WZ_Mu"+muID+"_El"+elID+".root"; TFile* wzfile = new TFile(wzfilename.c_str()); Sample* wzSample = new Sample(vBkgNames.at(4),wzfile, vWeights.at(4),vXsec.at(4),cut,kBlue-3); vSample.push_back(wzSample); - //TFile* wjfile = new TFile(bgTestDir.c_str()+"WJets_Mu"+muID+"_El"+elID+".root"); + //TFile* wjfile = new TFile(bgTestDir+"WJets_Mu"+muID+"_El"+elID+".root"); //Sample* wjSample = new Sample(vBkgNames.at(4),wjfile, vWeights.at(4),vXsec.at(4),cut,kGreen+2); //vSample.push_back(wjSample); - /*std::string dyfilename = bgTestDir.c_str()+"DYJets_Mu"+muID+"_El"+elID+".root"; + /*std::string dyfilename = bgTestDir+"DYJets_Mu"+muID+"_El"+elID+".root"; TFile* dyjfile = new TFile(dyfilename.c_str()); Sample* dyjSample = new Sample(vBkgNames.at(3),dyjfile, vWeights.at(3),vXsec.at(3),cut,kMagenta+2); vSample.push_back(dyjSample);*/ - std::string zzfilename = bgTestDir.c_str()+"ZZ_Mu"+muID+"_El"+elID+".root"; + std::string zzfilename = bgTestDir+"ZZ_Mu"+muID+"_El"+elID+".root"; TFile* zzjfile = new TFile(zzfilename.c_str()); Sample* zzjSample = new Sample(vBkgNames.at(5),zzjfile, vWeights.at(5),vXsec.at(5),cut,kOrange+1); //std::cout<<"weight for ZZ is: "< getBkgSampleVec(std::string cut, float lumi, std::string el //********** Nonprompt *************** - /*std::string npTTfilename = bgTestDir.c_str()+"NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; + /*std::string npTTfilename = bgTestDir+"NonPromptTTJets_Mu"+muID+"_El"+elID+".root"; TFile* npttfile = new TFile(nptTfilename.c_str()); Sample* npttSample = new Sample(vBkgNames.at(-1),npttfile,vWeights.at(-1),vXsec.at(01),cut,kBlue); vSample.push_back(npttSample)*/ - std::string npfilename = bgTestDir.c_str()+"NonPromptData_Mu"+muID+"_El"+elID+".root"; + std::string npfilename = bgTestDir+"NonPromptData_Mu"+muID+"_El"+elID+".root"; TFile* npfile = new TFile(npfilename.c_str()); Sample* npSample = new Sample(vBkgNames.at(12),npfile,vWeights.at(13),vXsec.at(13),cut,kGray); vSample.push_back(npSample); //********ChargeMisID********** - std::string cmidfilename = bgTestDir.c_str()+"ChargeMisID_Mu"+muID+"_El"+elID+".root"; + std::string cmidfilename = bgTestDir+"ChargeMisID_Mu"+muID+"_El"+elID+".root"; TFile* cmidfile = new TFile(cmidfilename.c_str()); Sample* cmidSample = new Sample(vBkgNames.at(13),cmidfile,vWeights.at(13),vXsec.at(13),cut,kAzure+6); //force charge misID to start here since only at this point do we filter events vSample.push_back(cmidSample); diff --git a/test/X53Analyzer.cc b/test/X53Analyzer.cc index 2161fe7..cf662b1 100644 --- a/test/X53Analyzer.cc +++ b/test/X53Analyzer.cc @@ -124,7 +124,29 @@ int main(int argc, char* argv[]){ sig_samples["X53X53m1500LH_Inc"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53X53m1500LH_Inc.root"; sig_samples["X53X53m1600LH_Inc"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53X53m1600LH_Inc.root"; - + // Dave ljmet trees PHYS14 + sig_samples["X53Tm700RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-700_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm800RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-800_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm900RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-900_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1000RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1000_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1100RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1100_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1200RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1200_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1300RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1300_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1400RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1400_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1500RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1500_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1600RH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1600_right_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm700LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-700_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm800LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-800_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm900LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-900_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1000LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1000_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1100LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1100_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1200LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1200_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1300LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1300_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1400LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1400_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1500LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1500_left_Tune4C_13TeV-madgraph-tauola.root"; + sig_samples["X53Tm1600LH"]="/eos/uscms/store/user/lpctlbsm/dzou/PHYS14/Single_SSDL/PU20/ljmet_trees/ljmet_tree_T53TTo2L2Nu_M-1600_left_Tune4C_13TeV-madgraph-tauola.root"; + + /* sig_samples["X53Tm700RH"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53TToAll_M-700_right.root"; sig_samples["X53Tm800RH"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53TToAll_M-800_right.root"; sig_samples["X53Tm900RH"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53TToAll_M-900_right.root"; @@ -145,6 +167,7 @@ int main(int argc, char* argv[]){ sig_samples["X53Tm1400LH"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53TToAll_M-1400_left.root"; sig_samples["X53Tm1500LH"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53TToAll_M-1500_left.root"; sig_samples["X53Tm1600LH"]="/eos/uscms/store/user/lpctlbsm/clint/Spring15/25ns/Approval/ljmet_trees/ljmet_X53TToAll_M-1600_left.root"; + */ //data samples data_samples["Data"]="/eos/uscms/store/user/lpctlbsm/clint/Run2015D/Nov17/ljmet_trees/ljmet_Data_All.root"; @@ -475,12 +498,15 @@ int main(int argc, char* argv[]){ //require OR of triggers but make each channel consistent + // dzou - comment out trigger requirements for use with PHYS14 ljmet - trigger menus have changed + bool skip = false; + /* bool skip = true; if(mumu && tr->HLT_Mu27TkMu8) skip =false; if(elmu && (tr->HLT_Mu8Ele17 || tr->HLT_Mu17Ele12)) skip = false; if(elel && tr->HLT_Ele17Ele12) skip = false; if(skip) continue; - +*/ //now skip if not the channel from the corresponding dataset if((argv1=="DataMuMu" || argv1=="NonPromptDataMuMu") && !mumu) continue; From a3300bb3b7e76f73ec5120a7747a879e08a2c59a Mon Sep 17 00:00:00 2001 From: David Date: Fri, 26 Aug 2016 10:37:54 -0500 Subject: [PATCH 3/3] Fixed typo from previous pull request --- plugins/Macros.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/plugins/Macros.cc b/plugins/Macros.cc index 4af2d21..4e51f6e 100644 --- a/plugins/Macros.cc +++ b/plugins/Macros.cc @@ -200,7 +200,7 @@ std::vector getNPHistos(int nMu){ Sample* getDataSample(std::string cut,std::string elID, std::string muID){ - std::string filename = testDir+muID+"_El"+elID+".root"; + std::string filename = testDir+"Data_Mu"+muID+"_El"+elID+".root"; TFile* tfile = new TFile(filename.c_str()); Sample* dataSample = new Sample("Data",tfile,1,1,cut,kBlack,1); return dataSample;