ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitHzz4l/LeptonSelection/src/MuonSelection.cc

Comparing UserCode/MitHzz4l/LeptonSelection/src/MuonSelection.cc (file contents):
Revision 1.22 by anlevin, Mon May 28 12:37:10 2012 UTC vs.
Revision 1.27 by anlevin, Wed Jun 6 18:17:21 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)
23	+	{
24	+
25	+	mithep::MuonTools   muT;
26	+
27	+	Float_t charged_iso = 0;
28	+	Float_t neutral_iso = 0;
29	+	Float_t gamma_iso = 0;
30	+
31	+	for(int k=0; k<pfCandidates->GetEntries(); ++k) {
32	+
33	+	if( !(PFnoPUflag[k]) ) continue;
34	+	const mithep::PFCandidate *pf = (mithep::PFCandidate*)((*pfCandidates)[k]);
35	+
36	+	Double_t dr = mithep::MathUtils::DeltaR(mu->Phi(),mu->Eta(), pf->Phi(), pf->Eta());
37	+	if (dr > 0.4) continue;
38	+
39	+	if (pf->Charge() != 0 && (pf->HasTrackerTrk()||pf->HasGsfTrk()) ) {
40	+
41	+	if (abs(pf->PFType()) == mithep::PFCandidate::eElectron || abs(pf->PFType()) == mithep::PFCandidate::eMuon) continue;
42	+
43	+	charged_iso += pf->Pt();
44	+	}
45	+
46	+	else if (abs(pf->PFType()) == mithep::PFCandidate::eGamma) {
47	+
48	+	gamma_iso += pf->Pt();
49	+	}
50	+
51	+	else {
52	+	neutral_iso += pf->Pt();
53	+	}
54	+	}
55	+
56	+	mithep::MuonTools::EMuonEffectiveAreaTarget EffectiveAreaVersion  = mithep::MuonTools::kMuEAData2012;
57	+	float rho = puEnergyDensity->At(0)->RhoKt6PFJetsCentralNeutral();
58	+
59	+	double iso_sum = charged_iso + fmax(0.0,(gamma_iso + neutral_iso
60	+	-rho*muT.MuonEffectiveArea(muT.kMuGammaAndNeutralHadronIso04,
61	+	mu->Eta(),EffectiveAreaVersion)));
62	+
63	+	bool pass = kTRUE;
64	+	if(!mu->IsGlobalMuon()) pass = kFALSE;
65	+
66	+	//there is not necessarily a global track or a tracker track if is global muon is false
67	+	if(!pass) return pass;
68	+
69	+	if(!mu->IsPFMuon()) pass = kFALSE;
70	+	if(mu->GlobalTrk()->RChi2() > 10) pass= kFALSE;
71	+	if(mu->NValidHits() == 0 ) pass= kFALSE;
72	+	if(mu->NMatches() <= 1) pass = kFALSE;
73	+	if(fabs(mu->TrackerTrk()->D0Corrected(*vtx)) > 0.2) pass = kFALSE;
74	+	if(fabs(mu->TrackerTrk()->DzCorrected(*vtx)) > 0.5) pass = kFALSE;
75	+	if(mu->BestTrk()->NPixelHits() == 0) pass = kFALSE;
76	+	if(mu->NTrkLayersHit() <= 5) pass = kFALSE;
77	+	if(iso_sum/mu->Pt() > 0.12 ) pass = kFALSE;
78	+
79	+	return pass;
80	+
81	+
82	+	}
83	+
84	+	bool trackToPFMuonSelection( const mithep::Track *track, const mithep::MuonCol * muoncol )
85	+	{
86	+	bool pass = false;
87	+	for(UInt_t i=0; i<muoncol->GetEntries(); ++i) {
88	+	const mithep::Muon *mu = muoncol->At(i);
89	+	if(mu->TrackerTrk() == track && (mu->IsGlobalMuon() || mu->IsTrackerMuon()) && mu->IsPFMuon())
90	+	pass = true;
91	+	}
92	+
93	+	return pass;
94	+
95	+	}
96	+
97		//--------------------------------------------------------------------------------------------------
98		SelectionStatus muonDummyVeto(ControlFlags &ctrl,
99		const mithep::Muon *muon,
#	Line 105 | Line 180 | SelectionStatus muonReferencePreSelectio
180		if(ctrl.debug) cout << "and < 2.4 eta (" << mu->Eta() << ")? ... " << ret << endl;
181		ret &= ( mu->IsTrackerMuon() || mu->IsGlobalMuon() );
182		if(ctrl.debug) cout << "is Tracker or Global? ... " << ret << endl;
183	<	ret &= (mu->HasTrackerTrk() && mu->TrackerTrk()->D0Corrected(*vtx) < 0.5 && mu->TrackerTrk()->DzCorrected(*vtx) < 1.0);
183	>	ret &= (mu->HasTrackerTrk() &&
184	>	fabs(mu->TrackerTrk()->D0Corrected(*vtx)) < 0.5 &&
185	>	fabs(mu->TrackerTrk()->DzCorrected(*vtx)) < 1.0);
186		if(ctrl.debug) cout << "d0 and dz? ... " << ret << endl;
187
188		//  if( mu->IsTrackerMuon() && mu->HasTrackerTrk() && !mu->IsGlobalMuon() )
189		//    ret &=  ( mu->Quality().QualityMask().Mask() & mithep::MuonQuality::AllArbitrated);
190	+	//  if(ctrl.debug) cout << "if isTrackerMuon, arbitrated ? ... " << ret << endl;
191	+	//  ret &= ( mu->BestTrk()->DzCorrected(vtx) < 0.1);
192	+	//  if( ctrl.debug ) cout << "elepresel : ret after dZcorr : " << ret <<  endl;
193	+
194	+
195	+	SelectionStatus status;
196	+	if( ret ) status.setStatus(SelectionStatus::PRESELECTION);
197	+	if(ctrl.debug) cout << "muon presel returning status : " << status.getStatus() << endl;
198	+	return status;
199	+	}
200	+
201	+	//--------------------------------------------------------------------------------------------------
202	+	SelectionStatus muonPreSelectionNoD0DzIP( ControlFlags &ctrl,
203	+	const mithep::Muon * mu,
204	+	const mithep::Vertex * vtx,
205	+	const mithep::Array<mithep::PFCandidate> * pfCandidates )
206	+	//--------------------------------------------------------------------------------------------------
207	+	{
208	+	bool ret = true;
209	+
210	+	if(ctrl.debug) cout << "inside muonpresel ... " << endl;
211	+	ret &= ( mu->Pt() >= 5 );
212	+	if(ctrl.debug) cout << "and >5 GeV (" << mu->Pt() << ") ? ... " << ret << endl;
213	+	ret &= ( fabs(mu->Eta()) <= 2.4 );
214	+	if(ctrl.debug) cout << "and < 2.4 eta (" << mu->Eta() << ")? ... " << ret << endl;
215	+	ret &= ( mu->IsTrackerMuon() || mu->IsGlobalMuon() );
216	+	if(ctrl.debug) cout << "is Tracker or Global? ... " << ret << endl;
217	+
218	+	//  if( mu->IsTrackerMuon() && mu->HasTrackerTrk() && !mu->IsGlobalMuon() )
219	+	//    ret &=  ( mu->Quality().QualityMask().Mask() & mithep::MuonQuality::AllArbitrated);
220		//  if(ctrl.debug) cout << "if isTrackerMuon, arbitrated ? ... " << ret << endl;
221		//  ret &= ( mu->BestTrk()->DzCorrected(vtx) < 0.1);
222		//  if( ctrl.debug ) cout << "elepresel : ret after dZcorr : " << ret <<  endl;

Diff Legend

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