ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitHzz4l/NonMCBackground/src/SelectionEMU.cc
(Generate patch)

Comparing UserCode/MitHzz4l/NonMCBackground/src/SelectionEMU.cc (file contents):
Revision 1.5 by dkralph, Tue Jun 12 22:39:31 2012 UTC vs.
Revision 1.8 by dkralph, Thu Jun 21 21:25:38 2012 UTC

# Line 106 | Line 106 | EventData apply_HZZ4L_EMU_selection(Cont
106      {
107        const mithep::Muon *mu = (mithep::Muon*)((*muonArr)[i]);      
108        
109 –      SelectionStatus musel;
109  
110 <      musel |= (*MuonPreSelector)(ctrl,mu,vtx,pfCandidates);
111 <      if( !(musel.getStatus() & SelectionStatus::PRESELECTION) ) continue;
110 >      SelectionStatus denomSel;
111 >      denomSel |= muonDenominatorSelection(ctrl,mu,vtx,pfCandidates);
112 >      if( !(denomSel.getStatus() & SelectionStatus::PRESELECTION) ) continue;
113  
114 +      SelectionStatus musel;
115 +      musel |= (*MuonPreSelector)(ctrl,mu,vtx,pfCandidates);
116        musel |= (*MuonIDSelector)(ctrl,mu,vtx,pfCandidates );
117 <
116 <      if( !(ctrl.doFSR) ) {
117 <        musel |=  (*MuonIsoSelector)(ctrl,mu,vtx,pfCandidates,puEnergyDensity,eraMu,photonsToVeto);
118 <      }
117 >      musel |= (*MuonIsoSelector)(ctrl,mu,vtx,pfCandidates,puEnergyDensity,eraMu,photonsToVeto);
118  
119        SimpleLepton tmplep;
120        float pt = mu->Pt();
# Line 144 | Line 143 | EventData apply_HZZ4L_EMU_selection(Cont
143        tmplep.isLoose = musel.loose();
144        tmplep.status  = musel;    
145        tmplep.fsrRecoveryAttempted = false;
146 <      tmplep.tightCutsApplied   = false;
146 >      SelectionStatus tmpstat = PassWwMuonSel(mu,vtx,pfCandidates);
147 >      tmplep.tightCutsApplied = tmpstat.tight();
148        lepvec.push_back(tmplep);
149        if( ctrl.debug ) cout << endl;
150      }
# Line 155 | Line 155 | EventData apply_HZZ4L_EMU_selection(Cont
155        {
156          const mithep::Electron *ele = (mithep::Electron*)((*electronArr)[i]);
157  
158 <        SelectionStatus elesel;
158 >        SelectionStatus denomSel;
159 >        denomSel |= electronDenominatorSelection(ctrl,ele,vtx);
160 >        if( !(denomSel.getStatus() & SelectionStatus::PRESELECTION) ) continue;
161          
162 +        SelectionStatus elesel;
163          elesel |= (*ElectronPreSelector)(ctrl,ele,vtx);
161 –        if( !(elesel.getStatus() & SelectionStatus::PRESELECTION) ) continue;
162 –        
164          elesel |= (*ElectronIDSelector)(ctrl,ele,vtx);
165 <        
165 <        if( !(ctrl.doFSR) ) {
166 <          elesel |= (*ElectronIsoSelector)(ctrl,ele,vtx,pfCandidates,puEnergyDensity,eraEle,photonsToVeto);
167 <        }
165 >        elesel |= (*ElectronIsoSelector)(ctrl,ele,vtx,pfCandidates,puEnergyDensity,eraEle,photonsToVeto);
166          
167          SimpleLepton tmplep;
168          float pt = ele->Pt();
# Line 215 | Line 213 | EventData apply_HZZ4L_EMU_selection(Cont
213        for (vector<SimpleLepton>::iterator it2=lepvec.begin(); it2 != lepvec.end(); it2++ ) {
214          if ( it2 == it1 )                       continue;
215          if ( abs(it2->type) != 13 )             continue;
216 <        if( !(it2->status.looseIDAndPre()) )    continue;
216 >        // if( !(it2->status.looseIDAndPre()) )    continue;
217          TVector3 mvec = it2->vec.Vect();
218          
219          if ( evec.DrEtaPhi(mvec) < 0.05 ) {
# Line 260 | Line 258 | EventData apply_HZZ4L_EMU_selection(Cont
258      //******************************************************************************
259      // W + (OF SS lepton) Selection
260      //******************************************************************************
261 <    if(ctrl.fakeScheme.Contains("emu-")) {
262 <      if(has_ssof_lepton(ctrl)) {
263 <        ret.status.setStatus(SelectionStatus::EVTPASS);
264 <        ret.Z1leptons.push_back(passingLeptons[0]);
265 <        ret.Z1leptons.push_back(passingLeptons[0]);
266 <        ret.Z2leptons.push_back(passingLeptons[0]);
267 <        ret.Z2leptons.push_back(passingLeptons[0]);
268 <      } else {
271 <        ret.status.setStatus(SelectionStatus::FAIL);
272 <      }
273 <      return ret;
261 >    if(has_ssof_lepton(ctrl)) {
262 >      ret.status.setStatus(SelectionStatus::EVTPASS);
263 >      ret.Z1leptons.push_back(passingLeptons[0]);
264 >      ret.Z1leptons.push_back(passingLeptons[0]);
265 >      ret.Z2leptons.push_back(passingLeptons[0]);
266 >      ret.Z2leptons.push_back(passingLeptons[0]);
267 >    } else {
268 >      ret.status.setStatus(SelectionStatus::FAIL);
269      }
270 +    return ret;
271   }
276 –
277 – // //----------------------------------------------------------------------------
278 – // //
279 – // // Get primary vertices
280 – // // Assumes primary vertices are ordered by sum-pT^2 (as should be in CMSSW)
281 – // // NOTE: if no PV is found from fitting tracks, the beamspot is used
282 – // //
283 – // //----------------------------------------------------------------------------
284 – // bool setPV(ControlFlags ctrl,
285 – //         const mithep::Array<mithep::Vertex> * vtxArr,
286 – //         const mithep::Vertex* &vtx)
287 – // //----------------------------------------------------------------------------
288 – // {
289 –
290 – //   const mithep::Vertex *bestPV = 0;    
291 – //   float best_sumpt=-1;
292 –
293 – //   // good PV requirements
294 – //   const UInt_t   fMinNTracksFit = 0;
295 – //   const Double_t fMinNdof       = 4;
296 – //   const Double_t fMaxAbsZ       = 24;
297 – //   const Double_t fMaxRho        = 2;
298 –  
299 – //   for(int i=0; i<vtxArr->GetEntries(); ++i) {
300 – //     const mithep::Vertex *pv = (mithep::Vertex*)(vtxArr->At(i));
301 – //     if( ctrl.debug ) cout << "vertex :: " << i << "\tntrks: " << pv->NTracks() << endl;
302 –    
303 – //     // Select best PV for corrected d0; if no PV passing cuts, the first PV in the collection will be used
304 – //     if(!pv->IsValid())                                continue;
305 – //     if(pv->NTracksFit()       < fMinNTracksFit)       continue;
306 – //     if(pv->Ndof()              < fMinNdof)         continue;
307 – //     if(fabs(pv->Z()) > fMaxAbsZ)                   continue;
308 – //     if(pv->Position().Rho()   > fMaxRho)           continue;    
309 –    
310 – //     // take the first ...
311 – //     bestPV = pv;
312 – //     break;
313 –
314 – //     // this never reached ...    
315 – //     float tmp_sumpt=0;
316 – //     for( int t=0; t<pv->NTracks(); t++ )
317 – //       tmp_sumpt += pv->Trk(t)->Pt();
318 –    
319 – //     if( tmp_sumpt > best_sumpt  ) {
320 – //       bestPV = pv;
321 – //       best_sumpt = tmp_sumpt;
322 – //       if( ctrl.debug) cout << "new PV set, pt : " << best_sumpt << endl;
323 – //     }
324 – //   }
325 –        
326 – //   // sync
327 – //   if(!bestPV)
328 – //     return false;
329 – //   else {    
330 – //     vtx = bestPV;
331 – //     return true;
332 – //   }  
333 – // };
272   //----------------------------------------------------------------------------------------
273   bool has_ssof_lepton(ControlFlags &ctrl)
274   {
# Line 361 | Line 299 | bool has_ssof_lepton(ControlFlags &ctrl)
299  
300    return has_ssof;
301   }
364 – // //----------------------------------------------------------------------------------------
365 – // void getEATargets(ControlFlags &ctrl, mithep::MuonTools::EMuonEffectiveAreaTarget &eraMu, mithep::ElectronTools::EElectronEffectiveAreaTarget &eraEle)
366 – // {
367 – //   if( !ctrl.mc && ctrl.era == 2011 ) {
368 – //     eraMu  = mithep::MuonTools::kMuEAData2011;
369 – //     eraEle = mithep::ElectronTools::kEleEAData2011;
370 – //   } else if( !ctrl.mc && ctrl.era == 2012 ) {
371 – //     eraMu  = mithep::MuonTools::kMuEAData2012;
372 – //     eraEle = mithep::ElectronTools::kEleEAData2012;
373 – //   } else if( ctrl.mc && ctrl.era == 2011 ) {
374 – //     eraMu  = mithep::MuonTools::kMuEAFall11MC;
375 – //     eraEle = mithep::ElectronTools::kEleEAFall11MC;
376 – //   } else if( ctrl.mc && ctrl.era == 2012 ) {
377 – //     eraMu  = mithep::MuonTools::kMuEAData2012;
378 – //     eraEle = mithep::ElectronTools::kEleEAData2012;
379 – //   } else {
380 – //     cerr << "unknown era for effective areas ... quitting." << endl;
381 – //     exit(1);
382 – //   }
383 – // }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines