//--------------------------------------------------------------------------------------------------
void executeComputeElectronFakeRateMod(const char *fileset = "",
const char *dataset = "s8-ttbar-id9",
const char *book = "mit/filler/006",
const char *catalogDir = "/home/mitprod/catalog",
const char *outputName = "JetChargeCorrelationStudy",
int nEvents = -1)
{
//------------------------------------------------------------------------------------------------
// some global setups
//------------------------------------------------------------------------------------------------
using namespace mithep;
gDebugMask = Debug::kAnalysis;
gDebugLevel = 3;
//------------------------------------------------------------------------------------------------
// generator information
//------------------------------------------------------------------------------------------------
GeneratorMod *generatorMod = new GeneratorMod;
generatorMod->SetFillHist(kTRUE);
//------------------------------------------------------------------------------------------------
// HLT information
//------------------------------------------------------------------------------------------------
HLTMod *hltmod = new HLTMod;
hltmod->SetPrintTable(kTRUE);
hltmod->AddTrigger("HLT_IsoEle15_L1I");
hltmod->AddTrigger("HLT_Ele15_SW_L1R");
hltmod->AddTrigger("HLT_IsoEle15_LW_L1I");
hltmod->AddTrigger("HLT_IsoEle18_L1R");
hltmod->AddTrigger("HLT_Ele15_LW_L1R");
hltmod->AddTrigger("HLT_LooseIsoEle15_LW_L1R");
hltmod->AddTrigger("HLT_IsoMu15");
hltmod->AddTrigger("HLT_Mu15");
hltmod->SetTrigObjsName("myhltobjs");
//------------------------------------------------------------------------------------------------
// object id and cleaning sequence
//------------------------------------------------------------------------------------------------
MuonIDMod *muonID = new MuonIDMod;
muonID->SetIDType("Loose");
muonID->SetIsoType("TrackCaloSliding");
ElectronIDMod *electronID = new ElectronIDMod;
electronID->SetIDType(TString("Tight"));
electronID->SetIsoType(TString("TrackJuraSliding"));
PhotonIDMod *photonID = new PhotonIDMod;
photonID->SetIDType(TString("Tight"));
photonID->SetPtMin(20.0);
JetIDMod *jetID = new JetIDMod;
jetID->SetInputName(TString("SisCone5Jets"));
jetID->SetUseCorrection(kTRUE);
jetID->SetPtCut(35.0);
ElectronCleaningMod *electronCleaning = new ElectronCleaningMod;
PhotonCleaningMod *photonCleaning = new PhotonCleaningMod;
JetCleaningMod *jetCleaning = new JetCleaningMod;
//------------------------------------------------------------------------------------------------
// merge modules
//------------------------------------------------------------------------------------------------
//MergeLeptonsMod *mergeLeptonsMod = new MergeLeptonsMod;
//------------------------------------------------------------------------------------------------
// analyses modules
//------------------------------------------------------------------------------------------------
JetChargeCorrelationStudyMod *analysisMod = new JetChargeCorrelationStudyMod;
analysisMod->SetVetoGeneratorLeptons(kTRUE);
//------------------------------------------------------------------------------------------------
// making analysis chain
//------------------------------------------------------------------------------------------------
hltmod->Add(muonID);
generatorMod->Add(muonID);
muonID->Add(electronID);
electronID->Add(photonID);
photonID->Add(jetID);
jetID->Add(electronCleaning);
electronCleaning->Add(photonCleaning);
photonCleaning->Add(jetCleaning);
jetCleaning->Add(analysisMod);
//------------------------------------------------------------------------------------------------
// setup analysis
//------------------------------------------------------------------------------------------------
Analysis *ana = new Analysis;
ana->SetUseHLT(kFALSE);
ana->SetKeepHierarchy(kTRUE);
if (nEvents >= 0)
ana->SetProcessNEvents(nEvents);
//ana->SetSuperModule(muonID);
ana->SetSuperModule(generatorMod);
//------------------------------------------------------------------------------------------------
// organize input
//------------------------------------------------------------------------------------------------
printf("\nRely on Catalog: %s\n",catalogDir);
printf(" -> Book: %s Dataset: %s Fileset: %s <-\n\n",book,dataset,fileset);
Catalog *c = new Catalog(catalogDir);
Dataset *d = c->FindDataset(book,dataset,fileset);
ana->AddDataset(d);
//------------------------------------------------------------------------------------------------
// organize output
//------------------------------------------------------------------------------------------------
TString rootFile = TString(outputName);
rootFile += TString("_") + TString(dataset);
if (TString(fileset) != TString(""))
rootFile += TString("_") + TString(fileset);
rootFile += TString(".root");
printf("\nRoot output: %s\n\n",rootFile.Data());
ana->SetOutputName(rootFile.Data());
//------------------------------------------------------------------------------------------------
// run the analysis after successful initialisation
//------------------------------------------------------------------------------------------------
ana->Run(!gROOT->IsBatch());
return;
}
//--------------------------------------------------------------------------------------------------
void executeComputeElectronFakeRateMod(const char *fileset = "",
const char *skim = "noskim",
const char *dataset = "s8-ttbar-id9",
const char *book = "mit/filler/006",
const char *catalogDir = "/home/mitprod/catalog",
const char *outputName = "ComputeElectronFakeRateMod",
Int_t sampletype = 0,
int nEvents = -1)
{
//------------------------------------------------------------------------------------------------
// some global setups
//------------------------------------------------------------------------------------------------
using namespace mithep;
gDebugMask = Debug::kAnalysis;
gDebugLevel = 3;
string mctype = "";
Bool_t VetoTriggerJet = false;
Bool_t VetoGenLeptons = true;
if (sampletype == 10) {
VetoTriggerJet = true;
VetoGenLeptons = false;
}
if (sampletype == 20) {
VetoGenLeptons = false;
}
if (sampletype == 3 || string(dataset) == "s8-wjets-mg-id11" )
mctype = "kMCTypeVLightJets";
if (sampletype == 4 || string(dataset) == "s8-vqq-mg-id11")
mctype = "kMCTypeVQQ";
if (sampletype == 5 || string(dataset) == "f8-wc-mg-id11")
mctype = "kMCTypeWc";
//------------------------------------------------------------------------------------------------
// generator information
//------------------------------------------------------------------------------------------------
GeneratorMod *generatorMod = new GeneratorMod;
//generatorMod->SetFillHist(kTRUE);
PartonFlavorHistoryMod *partonFlavorHistoryMod = new PartonFlavorHistoryMod;
partonFlavorHistoryMod->SetMCSampleType(mctype.c_str());
partonFlavorHistoryMod->SetApplyPartonFlavorFilter(kTRUE);
//------------------------------------------------------------------------------------------------
// HLT information
//------------------------------------------------------------------------------------------------
// HLTMod *hltmod = new HLTMod;
// hltmod->SetPrintTable(kTRUE);
// hltmod->AddTrigger("HLT_IsoEle15_L1I");
// hltmod->AddTrigger("HLT_Ele15_SW_L1R");
// hltmod->AddTrigger("HLT_IsoEle15_LW_L1I");
// hltmod->AddTrigger("HLT_IsoEle18_L1R");
// hltmod->AddTrigger("HLT_Ele15_LW_L1R");
// hltmod->AddTrigger("HLT_LooseIsoEle15_LW_L1R");
// hltmod->AddTrigger("HLT_IsoMu15");
// hltmod->AddTrigger("HLT_Mu15");
// hltmod->SetTrigObjsName("myhltobjs");
HLTMod *hltmod_Jet30 = new HLTMod;
hltmod_Jet30->SetName("HLTMod_Jet30");
hltmod_Jet30->SetPrintTable(kFALSE);
hltmod_Jet30->AddTrigger("HLT_Jet30");
hltmod_Jet30->SetTrigObjsName("myhltobjs_Jet30");
HLTMod *hltmod_Jet50 = new HLTMod;
hltmod_Jet50->SetName("HLTMod_Jet50");
hltmod_Jet50->SetPrintTable(kFALSE);
hltmod_Jet50->AddTrigger("HLT_Jet50");
hltmod_Jet50->SetTrigObjsName("myhltobjs_Jet50");
HLTMod *hltmod_Jet80 = new HLTMod;
hltmod_Jet80->SetName("HLTMod_Jet80");
hltmod_Jet80->SetPrintTable(kFALSE);
hltmod_Jet80->AddTrigger("HLT_Jet80");
hltmod_Jet80->SetTrigObjsName("myhltobjs_Jet80");
HLTMod *hltmod_Jet110 = new HLTMod;
hltmod_Jet110->SetName("HLTMod_Jet110");
hltmod_Jet110->SetPrintTable(kFALSE);
hltmod_Jet110->AddTrigger("HLT_Jet110");
hltmod_Jet110->SetTrigObjsName("myhltobjs_Jet110");
HLTMod *hltmod_Photon15 = new HLTMod;
hltmod_Photon15->SetName("HLTMod_Photon15");
hltmod_Photon15->SetPrintTable(kFALSE);
hltmod_Photon15->AddTrigger("HLT_Photon15_L1R");
hltmod_Photon15->SetTrigObjsName("myhltobjs_Photon15");
HLTMod *hltmod_Photon25 = new HLTMod;
hltmod_Photon25->SetName("HLTMod_Photon25");
hltmod_Photon25->SetPrintTable(kFALSE);
hltmod_Photon25->AddTrigger("HLT_Photon25_L1R");
hltmod_Photon25->SetTrigObjsName("myhltobjs_Photon25");
//------------------------------------------------------------------------------------------------
// publisher Mod
//------------------------------------------------------------------------------------------------
PublisherMod<CaloJet,Jet> *pubJet = new PublisherMod<CaloJet,Jet>("JetPub");
pubJet->SetInputName("SisCone5Jets");
pubJet->SetOutputName("PubSisCone5Jets");
PublisherMod<CaloMet,Met> *pubMet = new PublisherMod<CaloMet,Met>("MetPub");
pubMet->SetInputName("SisCone5Met");
pubMet->SetOutputName(ModNames::gkCleanCaloMetName);
//------------------------------------------------------------------------------------------------
// object id and cleaning sequence
//------------------------------------------------------------------------------------------------
MuonIDMod *muonID = new MuonIDMod;
muonID->SetIDType("Loose");
muonID->SetIsoType("TrackCaloSliding");
ElectronIDMod *electronID = new ElectronIDMod;
electronID->SetIDType(TString("Tight"));
electronID->SetIsoType(TString("TrackJuraSliding"));
electronID->SetApplyConversionFilter(kTRUE);
PhotonIDMod *photonID = new PhotonIDMod;
photonID->SetIDType(TString("Custom"));
photonID->SetIsoType(TString("CombinedIso"));
photonID->SetPtMin(20.0);
JetIDMod *jetID = new JetIDMod;
jetID->SetInputName(pubJet->GetOutputName());
jetID->SetUseCorrection(kTRUE);
jetID->SetPtCut(35.0);
ElectronCleaningMod *electronCleaning = new ElectronCleaningMod;
PhotonCleaningMod *photonCleaning = new PhotonCleaningMod;
JetCleaningMod *jetCleaning = new JetCleaningMod;
MergeLeptonsMod *mergeLeptonsMod = new MergeLeptonsMod;
//------------------------------------------------------------------------------------------------
// generate fakeable objects
//------------------------------------------------------------------------------------------------
GenerateFakeableObjectsMod *generateFakeableObjectsMod = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod->SetName("GenerateFakeableObjectsMod");
generateFakeableObjectsMod->SetElectronFOType("GsfPlusSC");
generateFakeableObjectsMod->SetMuonFOType("IsoTrack");
generateFakeableObjectsMod->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects");
generateFakeableObjectsMod->SetMuonFakeableObjectsName("IsoTrackMuonFakeableObjects");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Reco = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Reco->SetName("GenerateFakeableObjectsMod_Reco");
generateFakeableObjectsMod_Reco->SetElectronFOType("Reco");
generateFakeableObjectsMod_Reco->SetMuonFOType("Global");
generateFakeableObjectsMod_Reco->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Reco->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Reco->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Reco->SetElectronFakeableObjectsName("RecoElectronFakeableObjects");
generateFakeableObjectsMod_Reco->SetMuonFakeableObjectsName("GlobalMuonFakeableObjects");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Jet30 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Jet30->SetName("GenerateFakeableObjectsMod_Jet30");
generateFakeableObjectsMod_Jet30->SetElectronFOType("GsfPlusSC");
generateFakeableObjectsMod_Jet30->SetMuonFOType("IsoTrack");
generateFakeableObjectsMod_Jet30->SetTriggerName("HLTMod_Jet30");
generateFakeableObjectsMod_Jet30->SetTriggerObjectsName("myhltobjs_Jet30");
generateFakeableObjectsMod_Jet30->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Jet30->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Jet30->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Jet30->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet30");
generateFakeableObjectsMod_Jet30->SetMuonFakeableObjectsName("IsoTrackMuonFakeableObjects_Jet30");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Jet50 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Jet50->SetName("GenerateFakeableObjectsMod_Jet50");
generateFakeableObjectsMod_Jet50->SetElectronFOType("GsfPlusSC");
generateFakeableObjectsMod_Jet50->SetMuonFOType("IsoTrack");
generateFakeableObjectsMod_Jet50->SetTriggerName("HLTMod_Jet50");
generateFakeableObjectsMod_Jet50->SetTriggerObjectsName("myhltobjs_Jet50");
generateFakeableObjectsMod_Jet50->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Jet50->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Jet50->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Jet50->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet50");
generateFakeableObjectsMod_Jet50->SetMuonFakeableObjectsName("IsoTrackMuonFakeableObjects_Jet50");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Jet80 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Jet80->SetName("GenerateFakeableObjectsMod_Jet80");
generateFakeableObjectsMod_Jet80->SetElectronFOType("GsfPlusSC");
generateFakeableObjectsMod_Jet80->SetMuonFOType("IsoTrack");
generateFakeableObjectsMod_Jet80->SetTriggerName("HLTMod_Jet80");
generateFakeableObjectsMod_Jet80->SetTriggerObjectsName("myhltobjs_Jet80");
generateFakeableObjectsMod_Jet80->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Jet80->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Jet80->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Jet80->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet80");
generateFakeableObjectsMod_Jet80->SetMuonFakeableObjectsName("IsoTrackMuonFakeableObjects_Jet80");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Jet110 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Jet110->SetName("GenerateFakeableObjectsMod_Jet110");
generateFakeableObjectsMod_Jet110->SetElectronFOType("GsfPlusSC");
generateFakeableObjectsMod_Jet110->SetMuonFOType("IsoTrack");
generateFakeableObjectsMod_Jet110->SetTriggerName("HLTMod_Jet110");
generateFakeableObjectsMod_Jet110->SetTriggerObjectsName("myhltobjs_Jet110");
generateFakeableObjectsMod_Jet110->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Jet110->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Jet110->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Jet110->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet110");
generateFakeableObjectsMod_Jet110->SetMuonFakeableObjectsName("IsoTrackMuonFakeableObjects_Jet110");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Reco_Jet30 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Reco_Jet30->SetName("GenerateFakeableObjectsMod_Reco_Jet30");
generateFakeableObjectsMod_Reco_Jet30->SetElectronFOType("Reco");
generateFakeableObjectsMod_Reco_Jet30->SetMuonFOType("Global");
generateFakeableObjectsMod_Reco_Jet30->SetTriggerName("HLTMod_Jet30");
generateFakeableObjectsMod_Reco_Jet30->SetTriggerObjectsName("myhltobjs_Jet30");
generateFakeableObjectsMod_Reco_Jet30->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Reco_Jet30->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Reco_Jet30->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Reco_Jet30->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet30");
generateFakeableObjectsMod_Reco_Jet30->SetMuonFakeableObjectsName("GlobalMuonFakeableObjects_Jet30");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Reco_Jet50 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Reco_Jet50->SetName("GenerateFakeableObjectsMod_Reco_Jet50");
generateFakeableObjectsMod_Reco_Jet50->SetElectronFOType("Reco");
generateFakeableObjectsMod_Reco_Jet50->SetMuonFOType("Global");
generateFakeableObjectsMod_Reco_Jet50->SetTriggerName("HLTMod_Jet50");
generateFakeableObjectsMod_Reco_Jet50->SetTriggerObjectsName("myhltobjs_Jet50");
generateFakeableObjectsMod_Reco_Jet50->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Reco_Jet50->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Reco_Jet50->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Reco_Jet50->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet50");
generateFakeableObjectsMod_Reco_Jet50->SetMuonFakeableObjectsName("GlobalMuonFakeableObjects_Jet50");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Reco_Jet80 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Reco_Jet80->SetName("GenerateFakeableObjectsMod_Reco_Jet80");
generateFakeableObjectsMod_Reco_Jet80->SetElectronFOType("Reco");
generateFakeableObjectsMod_Reco_Jet80->SetMuonFOType("Global");
generateFakeableObjectsMod_Reco_Jet80->SetTriggerName("HLTMod_Jet80");
generateFakeableObjectsMod_Reco_Jet80->SetTriggerObjectsName("myhltobjs_Jet80");
generateFakeableObjectsMod_Reco_Jet80->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Reco_Jet80->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Reco_Jet80->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Reco_Jet80->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet80");
generateFakeableObjectsMod_Reco_Jet80->SetMuonFakeableObjectsName("GlobalMuonFakeableObjects_Jet80");
GenerateFakeableObjectsMod *generateFakeableObjectsMod_Reco_Jet110 = new GenerateFakeableObjectsMod;
generateFakeableObjectsMod_Reco_Jet110->SetName("GenerateFakeableObjectsMod_Reco_Jet110");
generateFakeableObjectsMod_Reco_Jet110->SetElectronFOType("Reco");
generateFakeableObjectsMod_Reco_Jet110->SetMuonFOType("Global");
generateFakeableObjectsMod_Reco_Jet110->SetTriggerName("HLTMod_Jet110");
generateFakeableObjectsMod_Reco_Jet110->SetTriggerObjectsName("myhltobjs_Jet110");
generateFakeableObjectsMod_Reco_Jet110->SetVetoTriggerJet(VetoTriggerJet);
generateFakeableObjectsMod_Reco_Jet110->SetVetoGenLeptons(VetoGenLeptons);
generateFakeableObjectsMod_Reco_Jet110->SetVetoCleanLeptons(kFALSE);
generateFakeableObjectsMod_Reco_Jet110->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet110");
generateFakeableObjectsMod_Reco_Jet110->SetMuonFakeableObjectsName("GlobalMuonFakeableObjects_Jet110");
//------------------------------------------------------------------------------------------------
// analyses modules
//------------------------------------------------------------------------------------------------
ComputeElectronFakeRateMod *analysisMod = new ComputeElectronFakeRateMod;
analysisMod->SetName("ComputeElectronFakeRateMod");
analysisMod->SetVetoTriggerJet(VetoTriggerJet);
analysisMod->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects");
analysisMod->SetMetName(pubMet->GetOutputName());
analysisMod->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco = new ComputeElectronFakeRateMod;
analysisMod_Reco->SetName("ComputeElectronFakeRateMod_Reco");
analysisMod_Reco->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco->SetElectronFakeableObjectsName("RecoElectronFakeableObjects");
analysisMod_Reco->SetMetName(pubMet->GetOutputName());
analysisMod_Reco->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Jet30 = new ComputeElectronFakeRateMod;
analysisMod_Jet30->SetName("ComputeElectronFakeRateMod_Jet30");
analysisMod_Jet30->SetTriggerName("HLT_Jet30");
analysisMod_Jet30->SetTriggerObjectsName("myhltobjs_Jet30");
analysisMod_Jet30->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Jet30->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Jet30->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet30");
analysisMod_Jet30->SetMetName(pubMet->GetOutputName());
analysisMod_Jet30->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Jet50 = new ComputeElectronFakeRateMod;
analysisMod_Jet50->SetName("ComputeElectronFakeRateMod_Jet50");
analysisMod_Jet50->SetTriggerName("HLT_Jet50");
analysisMod_Jet50->SetTriggerObjectsName("myhltobjs_Jet50");
analysisMod_Jet50->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Jet50->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Jet50->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet50");
analysisMod_Jet50->SetMetName(pubMet->GetOutputName());
analysisMod_Jet50->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Jet80 = new ComputeElectronFakeRateMod;
analysisMod_Jet80->SetName("ComputeElectronFakeRateMod_Jet80");
analysisMod_Jet80->SetTriggerName("HLT_Jet80");
analysisMod_Jet80->SetTriggerObjectsName("myhltobjs_Jet80");
analysisMod_Jet80->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Jet80->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Jet80->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet80");
analysisMod_Jet80->SetMetName(pubMet->GetOutputName());
analysisMod_Jet80->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Jet110 = new ComputeElectronFakeRateMod;
analysisMod_Jet110->SetName("ComputeElectronFakeRateMod_Jet110");
analysisMod_Jet110->SetTriggerName("HLT_Jet110");
analysisMod_Jet110->SetTriggerObjectsName("myhltobjs_Jet110");
analysisMod_Jet110->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Jet110->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Jet110->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects_Jet110");
analysisMod_Jet110->SetMetName(pubMet->GetOutputName());
analysisMod_Jet110->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco_Jet30 = new ComputeElectronFakeRateMod;
analysisMod_Reco_Jet30->SetName("ComputeElectronFakeRateMod_Reco_Jet30");
analysisMod_Reco_Jet30->SetTriggerName("HLT_Jet30");
analysisMod_Reco_Jet30->SetTriggerObjectsName("myhltobjs_Jet30");
analysisMod_Reco_Jet30->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco_Jet30->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco_Jet30->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet30");
analysisMod_Reco_Jet30->SetMetName(pubMet->GetOutputName());
analysisMod_Reco_Jet30->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco_Jet50 = new ComputeElectronFakeRateMod;
analysisMod_Reco_Jet50->SetName("ComputeElectronFakeRateMod_Reco_Jet50");
analysisMod_Reco_Jet50->SetTriggerName("HLT_Jet50");
analysisMod_Reco_Jet50->SetTriggerObjectsName("myhltobjs_Jet50");
analysisMod_Reco_Jet50->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco_Jet50->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco_Jet50->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet50");
analysisMod_Reco_Jet50->SetMetName(pubMet->GetOutputName());
analysisMod_Reco_Jet50->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco_Jet80 = new ComputeElectronFakeRateMod;
analysisMod_Reco_Jet80->SetName("ComputeElectronFakeRateMod_Reco_Jet80");
analysisMod_Reco_Jet80->SetTriggerName("HLT_Jet80");
analysisMod_Reco_Jet80->SetTriggerObjectsName("myhltobjs_Jet80");
analysisMod_Reco_Jet80->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco_Jet80->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco_Jet80->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet80");
analysisMod_Reco_Jet80->SetMetName(pubMet->GetOutputName());
analysisMod_Reco_Jet80->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco_Jet110 = new ComputeElectronFakeRateMod;
analysisMod_Reco_Jet110->SetName("ComputeElectronFakeRateMod_Reco_Jet110");
analysisMod_Reco_Jet110->SetTriggerName("HLT_Jet110");
analysisMod_Reco_Jet110->SetTriggerObjectsName("myhltobjs_Jet110");
analysisMod_Reco_Jet110->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco_Jet110->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco_Jet110->SetElectronFakeableObjectsName("RecoElectronFakeableObjects_Jet110");
analysisMod_Reco_Jet110->SetMetName(pubMet->GetOutputName());
analysisMod_Reco_Jet110->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Photon15 = new ComputeElectronFakeRateMod;
analysisMod_Photon15->SetName("ComputeElectronFakeRateMod_Photon15");
analysisMod_Photon15->SetTriggerName("HLT_Photon15_L1R");
analysisMod_Photon15->SetTriggerObjectsName("myhltobjs_Photon15");
analysisMod_Photon15->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Photon15->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Photon15->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects");
analysisMod_Photon15->SetMetName(pubMet->GetOutputName());
analysisMod_Photon15->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Photon25 = new ComputeElectronFakeRateMod;
analysisMod_Photon25->SetName("ComputeElectronFakeRateMod_Photon25");
analysisMod_Photon25->SetTriggerName("HLT_Photon25_L1R");
analysisMod_Photon25->SetTriggerObjectsName("myhltobjs_Photon25");
analysisMod_Photon25->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Photon25->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Photon25->SetElectronFakeableObjectsName("GsfPlusSCElectronFakeableObjects");
analysisMod_Photon25->SetMetName(pubMet->GetOutputName());
analysisMod_Photon25->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco_Photon15 = new ComputeElectronFakeRateMod;
analysisMod_Reco_Photon15->SetName("ComputeElectronFakeRateMod_Reco_Photon15");
analysisMod_Reco_Photon15->SetTriggerName("HLT_Photon15_L1R");
analysisMod_Reco_Photon15->SetTriggerObjectsName("myhltobjs_Photon15");
analysisMod_Reco_Photon15->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco_Photon15->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco_Photon15->SetElectronFakeableObjectsName("RecoElectronFakeableObjects");
analysisMod_Reco_Photon15->SetMetName(pubMet->GetOutputName());
analysisMod_Reco_Photon15->SetSampleType(sampletype);
ComputeElectronFakeRateMod *analysisMod_Reco_Photon25 = new ComputeElectronFakeRateMod;
analysisMod_Reco_Photon25->SetName("ComputeElectronFakeRateMod_Reco_Photon25");
analysisMod_Reco_Photon25->SetTriggerName("HLT_Photon25_L1R");
analysisMod_Reco_Photon25->SetTriggerObjectsName("myhltobjs_Photon25");
analysisMod_Reco_Photon25->SetVetoTriggerJet(VetoTriggerJet);
analysisMod_Reco_Photon25->SetVetoGeneratorLeptons(VetoGenLeptons);
analysisMod_Reco_Photon25->SetElectronFakeableObjectsName("RecoElectronFakeableObjects");
analysisMod_Reco_Photon25->SetMetName(pubMet->GetOutputName());
analysisMod_Reco_Photon25->SetSampleType(sampletype);
//------------------------------------------------------------------------------------------------
// making analysis chain
//------------------------------------------------------------------------------------------------
if (sampletype == 3 ||sampletype == 4 || sampletype == 5
|| string(dataset) == "s8-wjets-mg-id11" || string(dataset) == "s8-vqq-mg-id11"
) {
generatorMod->Add(partonFlavorHistoryMod);
partonFlavorHistoryMod->Add(pubJet);
} else {
generatorMod->Add(pubJet);
}
pubJet->Add(pubMet);
pubMet->Add(muonID);
muonID->Add(electronID);
electronID->Add(photonID);
photonID->Add(jetID);
jetID->Add(electronCleaning);
electronCleaning->Add(photonCleaning);
photonCleaning->Add(jetCleaning);
jetCleaning->Add(mergeLeptonsMod);
//For W+Jet samples
if (sampletype == 3 || sampletype == 4 || sampletype == 5 || sampletype == 1 ) {
mergeLeptonsMod->Add(generateFakeableObjectsMod);
generateFakeableObjectsMod->Add(analysisMod);
mergeLeptonsMod->Add(generateFakeableObjectsMod_Reco);
generateFakeableObjectsMod_Reco->Add(analysisMod_Reco);
} else if (sampletype == 20) {
//For photon trigger samples
mergeLeptonsMod->Add(hltmod_Photon15);
hltmod_Photon15->Add(generateFakeableObjectsMod);
generateFakeableObjectsMod->Add(analysisMod_Photon15);
hltmod_Photon15->Add(generateFakeableObjectsMod_Reco);
generateFakeableObjectsMod_Reco->Add(analysisMod_Reco_Photon15);
// mergeLeptonsMod->Add(hltmod_Photon25);
// hltmod_Photon25->Add(generateFakeableObjectsMod);
// generateFakeableObjectsMod->Add(analysisMod_Photon25);
// hltmod_Photon25->Add(generateFakeableObjectsMod_Reco);
// generateFakeableObjectsMod_Reco->Add(analysisMod_Reco_Photon25);
} else if (sampletype == 10) {
//For Jet Trigger samples
mergeLeptonsMod->Add(hltmod_Jet30);
hltmod_Jet30->Add(generateFakeableObjectsMod_Jet30);
generateFakeableObjectsMod_Jet30->Add(analysisMod_Jet30);
hltmod_Jet30->Add(generateFakeableObjectsMod_Reco_Jet30);
generateFakeableObjectsMod_Reco_Jet30->Add(analysisMod_Reco_Jet30);
mergeLeptonsMod->Add(hltmod_Jet50);
hltmod_Jet50->Add(generateFakeableObjectsMod_Jet50);
generateFakeableObjectsMod_Jet50->Add(analysisMod_Jet50);
hltmod_Jet50->Add(generateFakeableObjectsMod_Reco_Jet50);
generateFakeableObjectsMod_Reco_Jet50->Add(analysisMod_Reco_Jet50);
mergeLeptonsMod->Add(hltmod_Jet80);
hltmod_Jet80->Add(generateFakeableObjectsMod_Jet80);
generateFakeableObjectsMod_Jet80->Add(analysisMod_Jet80);
hltmod_Jet80->Add(generateFakeableObjectsMod_Reco_Jet80);
generateFakeableObjectsMod_Reco_Jet80->Add(analysisMod_Reco_Jet80);
mergeLeptonsMod->Add(hltmod_Jet110);
hltmod_Jet110->Add(generateFakeableObjectsMod_Jet110);
generateFakeableObjectsMod_Jet110->Add(analysisMod_Jet110);
hltmod_Jet110->Add(generateFakeableObjectsMod_Reco_Jet110);
generateFakeableObjectsMod_Reco_Jet110->Add(analysisMod_Reco_Jet110);
}
//------------------------------------------------------------------------------------------------
// setup analysis
//------------------------------------------------------------------------------------------------
Analysis *ana = new Analysis;
ana->SetUseHLT(kTRUE);
ana->SetKeepHierarchy(kTRUE);
if (nEvents >= 0)
ana->SetProcessNEvents(nEvents);
ana->SetSuperModule(generatorMod);
//ana->SetSkipFirstNEvents(47500);
//------------------------------------------------------------------------------------------------
// organize input
//------------------------------------------------------------------------------------------------
printf("\nRely on Catalog: %s\n",catalogDir);
printf(" -> Book: %s Dataset: %s Skim: %s Fileset: %s <-\n\n",book,dataset,skim,fileset);
Catalog *c = new Catalog(catalogDir);
TString skimdataset = TString(dataset)+TString("/") +TString(skim);
Dataset *d = NULL;
if (TString(skim).CompareTo("noskim") == 0)
d = c->FindDataset(book,dataset,fileset);
else
d = c->FindDataset(book,skimdataset.Data(),fileset);
ana->AddDataset(d);
//------------------------------------------------------------------------------------------------
// organize output
//------------------------------------------------------------------------------------------------
TString rootFile = TString(outputName);
rootFile += TString("_") + TString(dataset) + TString("_") + TString(skim);
if (TString(fileset) != TString(""))
rootFile += TString("_") + TString(fileset);
rootFile += TString(".root");
printf("\nRoot output: %s\n\n",rootFile.Data());
ana->SetOutputName(rootFile.Data());
//------------------------------------------------------------------------------------------------
// run the analysis after successful initialisation
//------------------------------------------------------------------------------------------------
ana->Run(!gROOT->IsBatch());
if (sampletype == 10) {
cout << "HLTMod Jet30: Accepted " << hltmod_Jet30->GetNAccepted()
<< " of " << hltmod_Jet30->GetNEvents() << " events" << endl;
cout << "HLTMod Jet50: Accepted " << hltmod_Jet50->GetNAccepted()
<< " of " << hltmod_Jet50->GetNEvents() << " events" << endl;
cout << "HLTMod Jet80: Accepted " << hltmod_Jet80->GetNAccepted()
<< " of " << hltmod_Jet80->GetNEvents() << " events" << endl;
cout << "HLTMod Jet110: Accepted " << hltmod_Jet110->GetNAccepted()
<< " of " << hltmod_Jet110->GetNEvents() << " events" << endl;
}
if (sampletype == 20) {
cout << "HLTMod Photon15: Accepted " << hltmod_Photon15->GetNAccepted()
<< " of " << hltmod_Photon15->GetNEvents() << " events" << endl;
cout << "HLTMod Photon25: Accepted " << hltmod_Photon25->GetNAccepted()
<< " of " << hltmod_Photon25->GetNEvents() << " events" << endl;
}
return;
}
//--------------------------------------------------------------------------------------------------
void runShowerCorrections(const char *fileset = "0000",
const char *skim = "noskim",
// const char *dataset = "s12-zllm50-v7n",
// const char *dataset = "r10b-pho-d22",
const char *dataset = "r12a-pho-j22-v1",
// const char *book = "local/filefi/031",
// const char *book = "t2mit/filefi/030",
const char *book = "t2mit/filefi/029",
const char *catalogDir = "/home/cmsprod/catalog",
const char *outputName = "hgg",
int nEvents = 10000)
{
//------------------------------------------------------------------------------------------------
// some parameters get passed through the environment
//------------------------------------------------------------------------------------------------
char json[1024], overlap[1024];
float overlapCut = -1;
if (gSystem->Getenv("MIT_PROD_JSON"))
sprintf(json, "%s",gSystem->Getenv("MIT_PROD_JSON"));
else {
printf(" JSON file was not properly defined. EXIT!\n");
return;
}
//TString jsonFile = TString("/home/cmsprod/json/") + TString(json);
TString jsonFile = TString("/home/fabstoec/cms/json/") + TString(json);
Bool_t isData = ( (jsonFile.CompareTo("/home/fabstoec/cms/json/~") != 0) );
if (gSystem->Getenv("MIT_PROD_OVERLAP")) {
sprintf(overlap,"%s",gSystem->Getenv("MIT_PROD_OVERLAP"));
if (EOF == sscanf(overlap,"%f",&overlapCut)) {
printf(" Overlap was not properly defined. EXIT!\n");
return;
}
}
else {
printf(" OVERLAP file was not properly defined. EXIT!\n");
return;
}
//------------------------------------------------------------------------------------------------
// some global setups
//------------------------------------------------------------------------------------------------
using namespace mithep;
gDebugMask = Debug::kGeneral;
gDebugLevel = 3;
//------------------------------------------------------------------------------------------------
// set up information
//------------------------------------------------------------------------------------------------
RunLumiSelectionMod *runLumiSel = new RunLumiSelectionMod;
runLumiSel->SetAcceptMC(kTRUE); // Monte Carlo events are always accepted
// only select on run- and lumisection numbers when valid json file present
if ((jsonFile.CompareTo("/home/fabstoec/cms/json/~") != 0) &&
(jsonFile.CompareTo("/home/fabstoec/cms/json/-") != 0) ) {
runLumiSel->AddJSONFile(jsonFile.Data());
}
if ((jsonFile.CompareTo("/home/cmsprod/json/-") == 0) ) {
printf("\n WARNING -- Looking at data without JSON file: always accept.\n\n");
runLumiSel->SetAbortIfNotAccepted(kFALSE); // accept all events if there is no valid JSON file
}
printf("\n Run lumi worked. \n\n");
//------------------------------------------------------------------------------------------------
// HLT information
//------------------------------------------------------------------------------------------------
HLTMod *hltModE = new HLTMod("HLTModE");
/// HLT paths for tag and probe in 2012
hltModE->AddTrigger("HLT_Ele17_CaloIdVT_CaloIsoVT_TrkIdT_TrkIsoVT_Ele8_Mass50_v*");
hltModE->AddTrigger("HLT_Ele20_CaloIdVT_CaloIsoVT_TrkIdT_TrkIsoVT_SC4_Mass50_v*");
hltModE->AddTrigger("HLT_Ele32_CaloIdT_CaloIsoT_TrkIdT_TrkIsoT_SC17_Mass50_v*");
hltModE->SetTrigObjsName("MyHltElecObjs");
hltModE->SetAbortIfNotAccepted(kFALSE);
//------------------------------------------------------------------------------------------------
// select events with a good primary vertex
//------------------------------------------------------------------------------------------------
GoodPVFilterMod *goodPVFilterMod = new GoodPVFilterMod;
goodPVFilterMod->SetMinVertexNTracks(0);
goodPVFilterMod->SetMinNDof (5);
goodPVFilterMod->SetMaxAbsZ (24.0);
goodPVFilterMod->SetMaxRho (2.0);
goodPVFilterMod->SetIsMC (!isData);
// our MIT Photons ID with iverted veto - we are looking at electrons
PhotonIDMod *photidpreselinvert = new PhotonIDMod("PhotonIDModPreselInvert");
photidpreselinvert->SetPtMin(15.);
photidpreselinvert->SetOutputName("GoodPhotonsPreselidInvert");
photidpreselinvert->SetIDType("MITPFSelectionNoEcal");
//photidpreselinvert->SetIDType("TrivialSelection");
photidpreselinvert->SetIdMVAType("2013FinalIdMVA");
photidpreselinvert->SetApplyElectronVeto(kFALSE);
photidpreselinvert->SetInvertElectronVeto(kTRUE);
photidpreselinvert->SetIsData(isData);
//photidpreselinvert->SetApplyFiduciality(kFALSE);
// some default VBTF Electron ID
ElectronIDMod *elecId = new ElectronIDMod;
elecId->SetIDType ("VBTFWorkingPoint80Id");
elecId->SetIsoType ("TrackJuraSliding");
elecId->SetApplyConversionFilterType1(kFALSE);
elecId->SetApplyConversionFilterType2(kTRUE);
elecId->SetChargeFilter (kFALSE);
elecId->SetApplyD0Cut (kTRUE);
elecId->SetNExpectedHitsInnerCut (0);
elecId->SetVertexName (goodPVFilterMod->GetOutputName());
elecId->SetPtMin (30.);
PhotonMCDataSFMod *phSFMod = new PhotonMCDataSFMod;
phSFMod-> SetIsData(isData);
phSFMod-> SetInputPhotonsName(photidpreselinvert->GetOutputName());
phSFMod-> SetPhotonsFromBranch(false);
phSFMod-> SetPVName(goodPVFilterMod->GetOutputName());
phSFMod-> SetPVFromBranch(false);
phSFMod-> SetElectronName(elecId->GetOutputName());
phSFMod-> SetEleFromBranch(false);
phSFMod-> AddTriggerLegs("hltEle32CaloIdLCaloIsoVLSC17PixelMatchFilter"/*,"hltEle32CaloIdLCaloIsoVLSC17HEDoubleFilter"*/);
phSFMod-> AddTriggerLegs("hltEle17CaloIdVTCaloIsoVTTrkIdTTrkIsoVTEle8TrackIsoFilter");
phSFMod-> AddTriggerLegs("hltEle20CaloIdVTCaloIsoVTTrkIdTTrkIsoVTSC4TrackIsoFilter");
phSFMod-> AddTriggerLegs("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsoFilter");
phSFMod-> SetTrigObjsName(hltModE->GetOutputName());
phSFMod-> SetTupleName("MITtuple");
//------------------------------------------------------------------------------------------------
// making analysis chain
//------------------------------------------------------------------------------------------------
// this is how it always starts
runLumiSel ->Add(hltModE);
hltModE ->Add(goodPVFilterMod);
goodPVFilterMod ->Add(elecId);
// ... the MIT flow ...
elecId ->Add(photidpreselinvert);
photidpreselinvert->Add(phSFMod);
//------------------------------------------------------------------------------------------------
// setup analysis
//------------------------------------------------------------------------------------------------
Analysis *ana = new Analysis;
ana->SetUseHLT(kTRUE);
ana->SetKeepHierarchy(kTRUE);
ana->SetSuperModule(runLumiSel);
ana->SetPrintScale(100);
if (nEvents >= 0)
ana->SetProcessNEvents(nEvents);
//------------------------------------------------------------------------------------------------
// organize input
//------------------------------------------------------------------------------------------------
// Catalog *c = new Catalog(catalogDir);
// TString skimdataset = TString(dataset)+TString("/") +TString(skim);
// Dataset *d = NULL;
// if (TString(skim).CompareTo("noskim") == 0)
// d = c->FindDataset(book,dataset,fileset);
// else
// d = c->FindDataset(book,skimdataset.Data(),fileset);
// ana->AddDataset(d);
Catalog *c = new Catalog(catalogDir);
TString skimdataset = TString(dataset)+TString("/") +TString(skim);
Dataset *d = NULL;
TString bookstr = book;
Bool_t caching = kTRUE;
//if (TString(dataset).Contains("s11-h")) bookstr.ReplaceAll("local","t2mit");
if (TString(skim).CompareTo("noskim") == 0)
d = c->FindDataset(bookstr,dataset,fileset,caching);
else
d = c->FindDataset(bookstr,skimdataset.Data(),fileset,caching);
ana->AddDataset(d);
//------------------------------------------------------------------------------------------------
// organize output
//------------------------------------------------------------------------------------------------
TString rootFile = TString(outputName);
rootFile += TString("_") + TString(dataset) + TString("_") + TString(skim);
if (TString(fileset) != TString(""))
rootFile += TString("_") + TString(fileset);
rootFile += TString(".root");
ana->SetOutputName(rootFile.Data());
// ana->SetCacheSize(64*1024*1024);
ana->SetCacheSize(0);
//------------------------------------------------------------------------------------------------
// Say what we are doing
//------------------------------------------------------------------------------------------------
printf("\n==== PARAMETER SUMMARY FOR THIS JOB ====\n");
printf("\n JSON file: %s\n and overlap cut: %f (%s)\n",jsonFile.Data(),overlapCut,overlap);
printf("\n Rely on Catalog: %s\n",catalogDir);
printf(" -> Book: %s Dataset: %s Skim: %s Fileset: %s <-\n",book,dataset,skim,fileset);
printf("\n Root output: %s\n\n",rootFile.Data());
printf("\n========================================\n");
//------------------------------------------------------------------------------------------------
// run the analysis after successful initialisation
//------------------------------------------------------------------------------------------------
ana->Run(!gROOT->IsBatch());
return;
}
