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.16 by khahn, Wed May 9 20:03:49 2012 UTC vs.
Revision 1.30 by dkralph, Thu Jul 5 13:43:13 2012 UTC

# Line 13 | Line 13
13   mithep::MuonIDMVA * muIDMVA;
14   mithep::MuonTools muTools;
15  
16 + extern vector<bool> PFnoPUflag;
17   extern Float_t computePFMuonIso(const mithep::Muon *mu,
18 <                        const mithep::Vertex & vtx,
18 >                        const mithep::Vertex * vtx,
19                          const mithep::Array<mithep::PFCandidate> * pfCandidates,
20                          const Double_t dRMax);
21  
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;
28 +
29 +  Float_t charged_iso = 0;
30 +  Float_t neutral_iso = 0;
31 +  Float_t gamma_iso = 0;
32 +  
33 +  for(int k=0; k<pfCandidates->GetEntries(); ++k) {
34 +    
35 +    if( !(PFnoPUflag[k]) ) continue;
36 +    const mithep::PFCandidate *pf = (mithep::PFCandidate*)((*pfCandidates)[k]);
37 +    
38 +    Double_t dr = mithep::MathUtils::DeltaR(mu->Phi(),mu->Eta(), pf->Phi(), pf->Eta());
39 +    if (dr > 0.4) continue;
40 +    
41 +    if (pf->Charge() != 0 && (pf->HasTrackerTrk()||pf->HasGsfTrk()) ) {
42 +      
43 +      if (abs(pf->PFType()) == mithep::PFCandidate::eElectron || abs(pf->PFType()) == mithep::PFCandidate::eMuon) continue;
44 +      
45 +      charged_iso += pf->Pt();
46 +    }
47 +    
48 +    else if (abs(pf->PFType()) == mithep::PFCandidate::eGamma) {
49 +      
50 +      gamma_iso += pf->Pt();
51 +    }
52 +    
53 +    else {
54 +      neutral_iso += pf->Pt();
55 +    }
56 +  }
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;
76 +
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() && 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 +
92 +
93 + }
94 +
95 + bool trackToPFMuonSelection( const mithep::Track *track, const mithep::MuonCol * muoncol )
96 + {
97 +        bool pass = false;
98 +  for(UInt_t i=0; i<muoncol->GetEntries(); ++i) {
99 +    const mithep::Muon *mu = muoncol->At(i);
100 +                if(mu->TrackerTrk() == track && (mu->IsGlobalMuon() || mu->IsTrackerMuon()) && mu->IsPFMuon())
101 +                        pass = true;
102 +        }              
103 +                
104 +        return pass;            
105 +                
106 + }
107 +
108   //--------------------------------------------------------------------------------------------------
109   SelectionStatus muonDummyVeto(ControlFlags &ctrl,
110                                const mithep::Muon *muon,
111 <                              const mithep::Vertex &vtx)
111 >                              const mithep::Vertex * vtx)
112   //--------------------------------------------------------------------------------------------------
113   {
114    SelectionStatus status;
# Line 32 | Line 119 | SelectionStatus muonDummyVeto(ControlFla
119   //--------------------------------------------------------------------------------------------------
120   SelectionStatus muonCutBasedVeto(ControlFlags &ctrl,
121                                   const mithep::Muon *muon,
122 <                                 const mithep::Vertex &vtx)
122 >                                 const mithep::Vertex *vtx)
123   //--------------------------------------------------------------------------------------------------
124   {
125    //
# Line 64 | Line 151 | SelectionStatus noPreSelection( ControlF
151   //--------------------------------------------------------------------------------------------------
152   SelectionStatus muonPreSelection( ControlFlags &ctrl,  
153                                    const mithep::Muon * mu,
154 <                                  const mithep::Vertex & vtx,
154 >                                  const mithep::Vertex * vtx,
155                                    const mithep::Array<mithep::PFCandidate> * pfCandidates )
156   //--------------------------------------------------------------------------------------------------
157   {
158    bool ret = true;
159 <  if(ctrl.debug) cout << "inside muonpresel ... " << endl;
160 <  //  bool ret = isMuFO(mu,vtx,pfCandidates);
74 <  if(ctrl.debug) cout << "iS fo? ... " << ret << endl;  
75 <  ret &= ( fabs(mu->Ip3dPVSignificance()) < 4 );
76 <  if(ctrl.debug) cout << "and pass IP (" << mu->Ip3dPVSignificance() << ") ? ... " << ret << endl;  
159 >
160 >  ret &= ( fabs(mu->Ip3dPVSignificance()) < 100 );
161    ret &= ( mu->Pt() >= 5 );
78 –  if(ctrl.debug) cout << "and >5 GeV (" << mu->Pt() << ") ? ... " << ret << endl;  
162    ret &= ( fabs(mu->Eta()) <= 2.4 );
80 –  if(ctrl.debug) cout << "and < 2.4 eta (" << mu->Eta() << ")? ... " << ret << endl;  
163    ret &= ( mu->IsTrackerMuon() || mu->IsGlobalMuon() );
164 <  if(ctrl.debug) cout << "is Tracker or Global? ... " << ret << endl;  
165 <  if( mu->IsTrackerMuon() && mu->HasTrackerTrk() && !mu->IsGlobalMuon() )  
166 <    ret &=  ( mu->Quality().QualityMask().Mask() & mithep::MuonQuality::AllArbitrated);
85 <  if(ctrl.debug) cout << "if isTrackerMuon, arbitrated ? ... " << ret << endl;  
86 <  ret &= ( fabs(mu->BestTrk()->DzCorrected(vtx)) < 0.1 );  
87 <  if( ctrl.debug ) cout << "elepresel : ret after dZcorr : " << ret <<  endl;
88 <
164 > //   if( mu->IsTrackerMuon() && mu->HasTrackerTrk() && !mu->IsGlobalMuon() )  
165 > //     ret &=  ( mu->Quality().QualityMask().Mask() & mithep::MuonQuality::AllArbitrated);
166 >  ret &= (  mu->IsoR03SumPt()/mu->Pt() < 0.7 );
167  
168    SelectionStatus status;
169    if( ret ) status.setStatus(SelectionStatus::PRESELECTION);
# Line 98 | Line 176 | SelectionStatus muonPreSelection( Contro
176   //--------------------------------------------------------------------------------------------------
177   SelectionStatus muonReferencePreSelection( ControlFlags &ctrl,  
178                                             const mithep::Muon * mu,
179 <                                           const mithep::Vertex & vtx,
179 >                                           const mithep::Vertex * vtx,
180                                             const mithep::Array<mithep::PFCandidate> * pfCandidates )
181   //--------------------------------------------------------------------------------------------------
182   {
183    bool ret = true;
184 +
185    if(ctrl.debug) cout << "inside muonpresel ... " << endl;
186    ret &= ( fabs(mu->Ip3dPVSignificance()) < 100 );
187    if(ctrl.debug) cout << "and pass IP (" << mu->Ip3dPVSignificance() << ") ? ... " << ret << endl;  
# Line 111 | Line 190 | SelectionStatus muonReferencePreSelectio
190    ret &= ( fabs(mu->Eta()) <= 2.4 );
191    if(ctrl.debug) cout << "and < 2.4 eta (" << mu->Eta() << ")? ... " << ret << endl;  
192    ret &= ( mu->IsTrackerMuon() || mu->IsGlobalMuon() );
193 <  if(ctrl.debug) cout << "is Tracker or Global? ... " << ret << endl;  
193 >  if(ctrl.debug) cout << "is Tracker or Global? ... " << ret << endl;
194 >  ret &= (mu->HasTrackerTrk() &&
195 >          fabs(mu->TrackerTrk()->D0Corrected(*vtx)) < 0.5 &&
196 >          fabs(mu->TrackerTrk()->DzCorrected(*vtx)) < 1.0);
197 >  if(ctrl.debug) cout << "d0 and dz? ... " << ret << endl;  
198 >          
199 >  //  if( mu->IsTrackerMuon() && mu->HasTrackerTrk() && !mu->IsGlobalMuon() )  
200 >  //    ret &=  ( mu->Quality().QualityMask().Mask() & mithep::MuonQuality::AllArbitrated);
201 >  //  if(ctrl.debug) cout << "if isTrackerMuon, arbitrated ? ... " << ret << endl;  
202 >  //  ret &= ( mu->BestTrk()->DzCorrected(vtx) < 0.1);  
203 >  //  if( ctrl.debug ) cout << "elepresel : ret after dZcorr : " << ret <<  endl;
204 >
205 >
206 >  SelectionStatus status;
207 >  if( ret ) status.setStatus(SelectionStatus::PRESELECTION);
208 >  if(ctrl.debug) cout << "muon presel returning status : " << status.getStatus() << endl;  
209 >  return status;
210 > }
211 > //--------------------------------------------------------------------------------------------------
212 > SelectionStatus muonPreSelectionNoD0DzIP( ControlFlags &ctrl,  
213 >                                          const mithep::Muon * mu,
214 >                                          const mithep::Vertex * vtx,
215 >                                          const mithep::Array<mithep::PFCandidate> * pfCandidates )
216 > //--------------------------------------------------------------------------------------------------
217 > {
218 >  bool ret = true;
219 >
220 >  if(ctrl.debug) cout << "inside muonpresel ... " << endl;
221 >  ret &= ( mu->Pt() >= 5 );
222 >  if(ctrl.debug) cout << "and >5 GeV (" << mu->Pt() << ") ? ... " << ret << endl;  
223 >  ret &= ( fabs(mu->Eta()) <= 2.4 );
224 >  if(ctrl.debug) cout << "and < 2.4 eta (" << mu->Eta() << ")? ... " << ret << endl;  
225 >  ret &= ( mu->IsTrackerMuon() || mu->IsGlobalMuon() );
226 >  if(ctrl.debug) cout << "is Tracker or Global? ... " << ret << endl;
227 >          
228    //  if( mu->IsTrackerMuon() && mu->HasTrackerTrk() && !mu->IsGlobalMuon() )  
229    //    ret &=  ( mu->Quality().QualityMask().Mask() & mithep::MuonQuality::AllArbitrated);
230    //  if(ctrl.debug) cout << "if isTrackerMuon, arbitrated ? ... " << ret << endl;  
# Line 131 | Line 244 | SelectionStatus muonReferencePreSelectio
244   //--------------------------------------------------------------------------------------------------
245   SelectionStatus muonIDMVASelection(ControlFlags &ctrl,
246                                   const mithep::Muon *mu,
247 <                                 const mithep::Vertex & vtx   )
247 >                                 const mithep::Vertex * vtx   )
248   //--------------------------------------------------------------------------------------------------
249   {
250    double global_rchi2 = (mu->IsGlobalMuon()) ? mu->GlobalTrk()->RChi2() : 0.;
# Line 205 | Line 318 | void initMuonIDMVA()
318   //--------------------------------------------------------------------------------------------------
319   SelectionStatus muonIDPFSelection(ControlFlags &ctrl,
320                                    const mithep::Muon *mu,
321 <                                  const mithep::Vertex & vtx,
321 >                                  const mithep::Vertex * vtx,
322                                    const mithep::Array<mithep::PFCandidate> * pfCandidates )
323   //--------------------------------------------------------------------------------------------------
324   {
# Line 214 | Line 327 | SelectionStatus muonIDPFSelection(Contro
327  
328    // check that it matches to a PF muon
329    for( int i=0; i<pfCandidates->GetEntries(); i++ ) {
330 +    // tmp
331 +    if( !(PFnoPUflag[i]) ) continue; // my PF no PU hack
332 +
333      const mithep::PFCandidate *pf = (mithep::PFCandidate*)((*pfCandidates)[i]);
334 <    if( pf->HasTrackerTrk() ) {
335 <      if( (pf->TrackerTrk() == mu->TrackerTrk()) && abs(pf->PFType()) == mithep::PFCandidate::eMuon ) {
336 <        pass = true;
221 <        break;
222 <      }
334 >    if( (pf->TrackerTrk() == mu->TrackerTrk()) && abs(pf->PFType()) == mithep::PFCandidate::eMuon ) {
335 >      pass = true;
336 >      break;
337      }
338    }
339 + //      
340 + //      if(mu->IsPFMuon()) pass = true;
341 + //      else pass = false;
342  
343 <  //
227 <  // NB : ipsig here for reference synch ...
228 <  //
229 <  if( pass && fabs(mu->Ip3dPVSignificance()) < 4 ) {
343 >  if( pass ) {
344      status.orStatus(SelectionStatus::LOOSEID);
345      status.orStatus(SelectionStatus::TIGHTID);
346    }
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