30 |
|
|
31 |
|
void mithep::MuonFakeRateMod::SlaveBegin() |
32 |
|
{ |
33 |
+ |
|
34 |
|
denominatorType.push_back(1); |
35 |
|
denominatorType.push_back(2); |
36 |
|
|
86 |
|
} |
87 |
|
|
88 |
|
rlrm.AddJSONFile("/afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/certification/Collisions12/8TeV/Prompt/Cert_190456-194479_8TeV_PromptReco_Collisions12_JSON.txt"); |
89 |
+ |
|
90 |
|
} |
91 |
|
|
92 |
|
void mithep::MuonFakeRateMod::SlaveTerminate() |
141 |
|
if (passMuonDenominatorCuts(mu, fPrimVerts->At(0), denominatorType[denominatorTypeIndex])) |
142 |
|
nMuons[denominatorTypeIndex]++; |
143 |
|
} |
144 |
< |
ControlFlags ctrl; |
143 |
< |
ctrl.era = 2012; |
144 |
> |
|
145 |
|
if ( //maybe you should do a loop to make sure the muon passes all of the denominators? right now this is ok because 2 is tighter... |
146 |
|
passMuonDenominatorCuts(mu, fPrimVerts->At(0), denominatorType[1]) |
147 |
|
&& |
148 |
< |
muon2012CutBasedIDTight(ctrl,mu,fPrimVerts->At(0),fPFCandidates,fPileupEnergyDensity,mithep::MuonTools::kMuEAData2012) |
148 |
> |
muonPOG2012CutBasedIDTight(mu,fPrimVerts->At(0),fPFCandidates,fPileupEnergyDensity) |
149 |
|
) |
150 |
|
goodMuons.push_back(mu); |
151 |
|
|
152 |
|
} |
153 |
|
for(Int_t i=0; i<fElectrons->GetEntries(); i++) { |
154 |
|
const Electron *ele = fElectrons->At(i); |
155 |
< |
ControlFlags ctrl; |
155 |
< |
ctrl.era = 2012; |
156 |
< |
if(electron2012CutBasedIDMedium(ctrl,ele, fPrimVerts->At(0), fPFCandidates, fConversions, fPileupEnergyDensity,mithep::ElectronTools::kEleEAData2012)) goodElectrons.push_back(ele); |
155 |
> |
if(electronPOG2012CutBasedIDMedium(ele, fPrimVerts->At(0), fPFCandidates, fConversions, fPileupEnergyDensity)) goodElectrons.push_back(ele); |
156 |
|
|
157 |
|
} |
158 |
|
|
292 |
|
if(!muonIDPFSelection(ctrl, mu, fPrimVerts->At(0), fPFCandidates ).looseID()) passNumerator = kFALSE; |
293 |
|
if(fabs(mu->Ip3dPVSignificance()) > 4) passNumerator = kFALSE; |
294 |
|
if(!muonReferenceIsoSelection(ctrl,mu,fPrimVerts->At(0),fPFCandidates,fPileupEnergyDensity,eraMu,photonsToVeto).passLooseIso()) passNumerator = kFALSE; |
295 |
< |
|
296 |
< |
if (passNumerator) |
298 |
< |
fake_muon_pass[denominatorTypeIndex] = passNumerator ? 1 : 0; |
295 |
> |
|
296 |
> |
fake_muon_pass[denominatorTypeIndex] = passNumerator ? 1 : 0; |
297 |
|
|
298 |
|
fOutputTrees[denominatorTypeIndex]->Fill(); |
299 |
|
|