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

Comparing UserCode/MitHzz4l/LeptonSelection/src/MuonSelection.cc (file contents):
Revision 1.27 by anlevin, Wed Jun 6 18:17:21 2012 UTC vs.
Revision 1.30 by dkralph, Thu Jul 5 13:43:13 2012 UTC

# Line 19 | Line 19 | extern Float_t computePFMuonIso(const mi
19                          const mithep::Array<mithep::PFCandidate> * pfCandidates,
20                          const Double_t dRMax);
21  
22 < bool muon2012CutBasedIDTight(const mithep::Muon *mu, const mithep::Vertex * vtx, const mithep::Array<mithep::PFCandidate> * pfCandidates, const mithep::Array<mithep::PileupEnergyDensity> * puEnergyDensity)
22 > bool muon2012CutBasedIDTight(ControlFlags &ctrl, const mithep::Muon *mu, const mithep::Vertex * vtx, const mithep::Array<mithep::PFCandidate> * pfCandidates,
23 >                             const mithep::Array<mithep::PileupEnergyDensity> * puEnergyDensity,   mithep::MuonTools::EMuonEffectiveAreaTarget EffectiveAreaVersion
24 >                             )
25   {
26  
27    mithep::MuonTools   muT;
# Line 52 | Line 54 | bool muon2012CutBasedIDTight(const mithe
54        neutral_iso += pf->Pt();
55      }
56    }
57 <  
58 <  mithep::MuonTools::EMuonEffectiveAreaTarget EffectiveAreaVersion  = mithep::MuonTools::kMuEAData2012;
59 <  float rho = puEnergyDensity->At(0)->RhoKt6PFJetsCentralNeutral();
60 <  
61 <  double iso_sum = charged_iso + fmax(0.0,(gamma_iso + neutral_iso
62 <                                        -rho*muT.MuonEffectiveArea(muT.kMuGammaAndNeutralHadronIso04,
63 <                                                                   mu->Eta(),EffectiveAreaVersion)));
57 >
58 >  double iso_sum=999;
59 >  double rho;
60 >  if(ctrl.era == 2012) {
61 >    EffectiveAreaVersion  = mithep::MuonTools::kMuEAData2012;
62 >    rho = puEnergyDensity->At(0)->RhoKt6PFJetsCentralNeutral();
63 >    
64 >    iso_sum = charged_iso + fmax(0.0,(gamma_iso + neutral_iso
65 >                                      -rho*muT.MuonEffectiveArea(muT.kMuGammaAndNeutralHadronIso04,
66 >                                                                 mu->Eta(),EffectiveAreaVersion)));
67 >  } else if(ctrl.era == 2011) {
68 >    setMuonRhoAndEaVersion(rho,EffectiveAreaVersion,puEnergyDensity);
69 >    iso_sum = charged_iso + fmax(0.0,(gamma_iso + neutral_iso
70 >                                      -rho*muT.MuonEffectiveArea(muT.kMuGammaAndNeutralHadronIso04,
71 >                                                                 mu->Eta(),EffectiveAreaVersion)));
72 >  } else assert(0);
73  
74    bool pass = kTRUE;
75    if(!mu->IsGlobalMuon()) pass = kFALSE;
# Line 66 | Line 77 | bool muon2012CutBasedIDTight(const mithe
77    //there is not necessarily a global track or a tracker track if is global muon is false
78    if(!pass) return pass;
79  
80 <  if(!mu->IsPFMuon()) pass = kFALSE;
81 <  if(mu->GlobalTrk()->RChi2() > 10) pass= kFALSE;
82 <  if(mu->NValidHits() == 0 ) pass= kFALSE;
83 <  if(mu->NMatches() <= 1) pass = kFALSE;
84 <  if(fabs(mu->TrackerTrk()->D0Corrected(*vtx)) > 0.2) pass = kFALSE;
85 <  if(fabs(mu->TrackerTrk()->DzCorrected(*vtx)) > 0.5) pass = kFALSE;
86 <  if(mu->BestTrk()->NPixelHits() == 0) pass = kFALSE;
87 <  if(mu->NTrkLayersHit() <= 5) pass = kFALSE;
88 <  if(iso_sum/mu->Pt() > 0.12 ) pass = kFALSE;
80 >  if(!mu->IsPFMuon() && ctrl.era==2012)                 pass = kFALSE;
81 >  if(mu->GlobalTrk()->RChi2() > 10)                     pass = kFALSE;
82 >  if(mu->NValidHits() == 0 )                            pass = kFALSE;
83 >  if(mu->NMatches() <= 1)                               pass = kFALSE;
84 >  if(fabs(mu->TrackerTrk()->D0Corrected(*vtx)) > 0.2)   pass = kFALSE;
85 >  if(fabs(mu->TrackerTrk()->DzCorrected(*vtx)) > 0.5)   pass = kFALSE;
86 >  if(mu->BestTrk()->NPixelHits() == 0)                  pass = kFALSE;
87 >  if(mu->NTrkLayersHit() <= 5)                          pass = kFALSE;
88 >  if(iso_sum/mu->Pt() > 0.12 )                          pass = kFALSE;
89  
90    return pass;
91  
# Line 197 | Line 208 | SelectionStatus muonReferencePreSelectio
208    if(ctrl.debug) cout << "muon presel returning status : " << status.getStatus() << endl;  
209    return status;
210   }
200 –
211   //--------------------------------------------------------------------------------------------------
212   SelectionStatus muonPreSelectionNoD0DzIP( ControlFlags &ctrl,  
213                                            const mithep::Muon * mu,
# Line 330 | Line 340 | SelectionStatus muonIDPFSelection(Contro
340   //      if(mu->IsPFMuon()) pass = true;
341   //      else pass = false;
342  
333 –  //
334 –  // NB : ipsig here for reference synch ...
335 –  //
336 – #ifdef SYNC
337 –  if( pass && fabs(mu->Ip3dPVSignificance()) < 4 ) {
338 –    status.orStatus(SelectionStatus::LOOSEID);
339 –    status.orStatus(SelectionStatus::TIGHTID);
340 –  }
341 – #else
343    if( pass ) {
344      status.orStatus(SelectionStatus::LOOSEID);
345      status.orStatus(SelectionStatus::TIGHTID);
346    }
346 – #endif
347  
348    return status;
349   }
350 + //----------------------------------------------------------------------------------------
351 + SelectionStatus PassWwMuonSel(const mithep::Muon *mu, const mithep::Vertex *vtx,
352 +                              const mithep::Array<mithep::PFCandidate> *pfArr
353 + )
354 + {
355 +  bool passID=true;
356 +  bool passIso=true;
357 +  
358 +  if (mu->Pt() < 5)
359 +    passID = false;
360 +  if (fabs(mu->Eta()) > 2.4)
361 +    passID = false;
362 +  if (mu->BestTrk()->PtErr()/mu->Pt() > 0.1)
363 +    passID = false;
364 +  if (fabs(mu->BestTrk()->DzCorrected(*vtx)) > 0.1)
365 +    passID = false;
366 +  Bool_t isGlobal  = (mu->IsGlobalMuon()) && (mu->BestTrk()->RChi2() < 10) &&
367 +                      (mu->NMatches() > 1) && (mu->NValidHits() > 0);
368 +  Bool_t isTracker = (mu->IsTrackerMuon()) &&
369 +                      (mu->Quality().Quality(mithep::MuonQuality::TMLastStationTight));
370 +  if (!isGlobal && !isTracker)
371 +    passID = false;
372 +  int ntrkhits = (mu->HasTrackerTrk()) ? mu->TrackerTrk()->NHits() : 0;
373 +  if (ntrkhits < 10)
374 +    passID = false;
375 +  if (mu->BestTrk()->NPixelHits() < 1)
376 +    passID = false;
377 +  if (fabs(mu->BestTrk()->D0Corrected(*vtx)) > 0.02)
378 +    passID = false;
379 +  
380 +  // note: this isn't really ww muon isolation
381 +  float pfiso = computePFMuonIso(mu,vtx,pfArr,0.3);
382 +
383 +  if (pfiso > 0.2*mu->Pt())
384 +    passIso = false;
385 +
386 +  SelectionStatus status;
387 +  if(passID)  status |= SelectionStatus::TIGHTIDANDPRE;
388 +  if(passIso) status |= SelectionStatus::TIGHTISOANDPRE;
389 +  
390 +  return status;
391 +
392 + }

Diff Legend

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