| 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() |
| 128 |
|
|
| 129 |
|
fillTriggerBits( GetHLTTable(), fTrigMask, fTriggerBits ); |
| 130 |
|
|
| 131 |
+ |
PFnoPUflag.clear(); |
| 132 |
|
makePFnoPUArray(fPFCandidates, PFnoPUflag, fPrimVerts ); |
| 133 |
|
|
| 134 |
|
if(fPFMet->At(0)->Et() >25) return; |
| 142 |
|
if (passMuonDenominatorCuts(mu, fPrimVerts->At(0), denominatorType[denominatorTypeIndex])) |
| 143 |
|
nMuons[denominatorTypeIndex]++; |
| 144 |
|
} |
| 145 |
< |
ControlFlags ctrl; |
| 143 |
< |
ctrl.era = 2012; |
| 145 |
> |
|
| 146 |
|
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... |
| 147 |
|
passMuonDenominatorCuts(mu, fPrimVerts->At(0), denominatorType[1]) |
| 148 |
|
&& |
| 149 |
< |
muon2012CutBasedIDTight(ctrl,mu,fPrimVerts->At(0),fPFCandidates,fPileupEnergyDensity,mithep::MuonTools::kMuEAData2012) |
| 149 |
> |
muonPOG2012CutBasedIDTight(mu,fPrimVerts->At(0),fPFCandidates,fPileupEnergyDensity) |
| 150 |
|
) |
| 151 |
|
goodMuons.push_back(mu); |
| 152 |
|
|
| 153 |
|
} |
| 154 |
|
for(Int_t i=0; i<fElectrons->GetEntries(); i++) { |
| 155 |
|
const Electron *ele = fElectrons->At(i); |
| 156 |
< |
ControlFlags ctrl; |
| 155 |
< |
ctrl.era = 2012; |
| 156 |
< |
if(electron2012CutBasedIDMedium(ctrl,ele, fPrimVerts->At(0), fPFCandidates, fConversions, fPileupEnergyDensity,mithep::ElectronTools::kEleEAData2012)) goodElectrons.push_back(ele); |
| 156 |
> |
if(electronPOG2012CutBasedIDMedium(ele, fPrimVerts->At(0), fPFCandidates, fConversions, fPileupEnergyDensity)) goodElectrons.push_back(ele); |
| 157 |
|
|
| 158 |
|
} |
| 159 |
|
|
| 293 |
|
if(!muonIDPFSelection(ctrl, mu, fPrimVerts->At(0), fPFCandidates ).looseID()) passNumerator = kFALSE; |
| 294 |
|
if(fabs(mu->Ip3dPVSignificance()) > 4) passNumerator = kFALSE; |
| 295 |
|
if(!muonReferenceIsoSelection(ctrl,mu,fPrimVerts->At(0),fPFCandidates,fPileupEnergyDensity,eraMu,photonsToVeto).passLooseIso()) passNumerator = kFALSE; |
| 296 |
< |
|
| 297 |
< |
if (passNumerator) |
| 298 |
< |
fake_muon_pass[denominatorTypeIndex] = passNumerator ? 1 : 0; |
| 296 |
> |
|
| 297 |
> |
fake_muon_pass[denominatorTypeIndex] = passNumerator ? 1 : 0; |
| 298 |
|
|
| 299 |
|
fOutputTrees[denominatorTypeIndex]->Fill(); |
| 300 |
|
|
