ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitPhysics/Mods/src/MuonIDMod.cc

Comparing UserCode/MitPhysics/Mods/src/MuonIDMod.cc (file contents):
Revision 1.8 by loizides, Fri Nov 28 09:56:42 2008 UTC vs.
Revision 1.41 by ceballos, Fri Mar 11 15:13:09 2011 UTC

#	Line 2 | Line 2
2
3		#include "MitPhysics/Mods/interface/MuonIDMod.h"
4		#include "MitCommon/MathTools/interface/MathUtils.h"
5	+	#include "MitAna/DataTree/interface/MuonFwd.h"
6	+	#include "MitAna/DataTree/interface/ElectronFwd.h"
7	+	#include "MitAna/DataTree/interface/VertexCol.h"
8		#include "MitPhysics/Init/interface/ModNames.h"
9
10		using namespace mithep;
#	Line 13 | Line 16 | ClassImp(mithep::MuonIDMod)
16		BaseMod(name,title),
17		fMuonBranchName(Names::gkMuonBrn),
18		fCleanMuonsName(ModNames::gkCleanMuonsName),
19	<	fMuonIDType("Loose"),
20	<	fMuonIsoType("TrackCaloSliding"),
19	>	fNonIsolatedMuonsName("random"),
20	>	fNonIsolatedElectronsName("random"),
21	>	fVertexName(ModNames::gkGoodVertexesName),
22	>	fTrackName(Names::gkTrackBrn),
23	>	fPFCandidatesName(Names::gkPFCandidatesBrn),
24	>	fMuonIDType("WWMuId"),
25	>	fMuonIsoType("TrackCaloSliding"),
26		fMuonClassType("Global"),
27		fTrackIsolationCut(3.0),
28		fCaloIsolationCut(3.0),
29	<	fCombIsolationCut(5.0),
29	>	fCombIsolationCut(0.15),
30		fMuonPtMin(10),
31	<	fMuons(0)
31	>	fApplyD0Cut(kTRUE),
32	>	fD0Cut(0.020),
33	>	fEtaCut(2.4),
34	>	fReverseIsoCut(kFALSE),
35	>	fReverseD0Cut(kFALSE),
36	>	fMuIDType(kIdUndef),
37	>	fMuIsoType(kIsoUndef),
38	>	fMuClassType(kClassUndef),
39	>	fMuons(0),
40	>	fVertices(0),
41	>	fTracks(0),
42	>	fPFCandidates(0),
43	>	fNonIsolatedMuons(0),
44	>	fNonIsolatedElectrons(0),
45	>	fPileupEnergyDensityName(Names::gkPileupEnergyDensityBrn),
46	>	fPileupEnergyDensity(0)
47		{
48		// Constructor.
49		}
50
28	–
51		//--------------------------------------------------------------------------------------------------
52		void MuonIDMod::Process()
53		{
54		// Process entries of the tree.
55
56	<	LoadBranch(fMuonBranchName);
57	<
56	>	if(fMuIsoType != kPFIsoNoL) {
57	>	LoadEventObject(fMuonBranchName, fMuons);
58	>	}
59	>	else {
60	>	fMuons = GetObjThisEvt<MuonOArr>(fMuonBranchName);
61	>	}
62	>	LoadEventObject(fTrackName, fTracks);
63	>	LoadEventObject(fPFCandidatesName, fPFCandidates);
64	>	if(fMuIsoType == kTrackCaloSliding) {
65	>	LoadEventObject(fPileupEnergyDensityName, fPileupEnergyDensity);
66	>	}
67		MuonOArr *CleanMuons = new MuonOArr;
68		CleanMuons->SetName(fCleanMuonsName);
69
70	+	fVertices = GetObjThisEvt<VertexOArr>(fVertexName);
71	+
72		for (UInt_t i=0; i<fMuons->GetEntries(); ++i) {
73		const Muon *mu = fMuons->At(i);
74
75		Bool_t pass = kFALSE;
76	<	Double_t pt = -1; // make sure pt is taken from the correct track!
76	>	Double_t pt = 0;  // make sure pt is taken from the correct track!
77	>	Double_t eta = 0; // make sure eta is taken from the correct track!
78		switch (fMuClassType) {
79		case kAll:
80		pass = kTRUE;
81	<	pt = mu->Pt();
81	>	if (mu->HasTrk()) {
82	>	pt  = mu->Pt();
83	>	eta = TMath::Abs(mu->Eta());
84	>	}
85		break;
86		case kGlobal:
87	<	pass = (mu->GlobalTrk() != 0);
88	<	if (pass)
89	<	pt = mu->GlobalTrk()->Pt();
90	<	break;
87	>	pass = mu->HasGlobalTrk() && mu->IsTrackerMuon();
88	>	if (pass && mu->TrackerTrk()) {
89	>	pt  = mu->TrackerTrk()->Pt();
90	>	eta = TMath::Abs(mu->TrackerTrk()->Eta());
91	>	}
92	>	else {
93	>	pt  = mu->Pt();
94	>	eta = TMath::Abs(mu->Eta());
95	>	}
96	>	break;
97		case kSta:
98	<	pass = (mu->StandaloneTrk() != 0);
99	<	if (pass)
100	<	pt = mu->StandaloneTrk()->Pt();
101	<	break;
102	<	case kTrackerOnly:
103	<	pass = (mu->TrackerTrk() != 0);
104	<	if (pass)
105	<	pt = mu->TrackerTrk()->Pt();
98	>	pass = mu->HasStandaloneTrk();
99	>	if (pass) {
100	>	pt  = mu->StandaloneTrk()->Pt();
101	>	eta = TMath::Abs(mu->StandaloneTrk()->Eta());
102	>	}
103	>	break;
104	>	case kTrackerMuon:
105	>	pass = mu->HasTrackerTrk() && mu->IsTrackerMuon() &&
106	>	mu->Quality().Quality(MuonQuality::TrackerMuonArbitrated);
107	>	if (pass) {
108	>	pt  = mu->TrackerTrk()->Pt();
109	>	eta = TMath::Abs(mu->TrackerTrk()->Eta());
110	>	}
111	>	break;
112	>	case kCaloMuon:
113	>	pass = mu->HasTrackerTrk() && mu->IsCaloMuon();
114	>	if (pass) {
115	>	pt  = mu->TrackerTrk()->Pt();
116	>	eta = TMath::Abs(mu->TrackerTrk()->Eta());
117	>	}
118	>	break;
119	>	case kTrackerBased:
120	>	pass = mu->HasTrackerTrk();
121	>	if (pass) {
122	>	pt  = mu->TrackerTrk()->Pt();
123	>	eta = TMath::Abs(mu->TrackerTrk()->Eta());
124	>	}
125		break;
126		default:
127		break;
#	Line 71 | Line 133 | void MuonIDMod::Process()
133		if (pt <= fMuonPtMin)
134		continue;
135
136	+	if (eta >= fEtaCut)
137	+	continue;
138	+
139	+	Double_t RChi2 = 0.0;
140	+	if     (mu->HasGlobalTrk()) {
141	+	RChi2 = mu->GlobalTrk()->Chi2()/mu->GlobalTrk()->Ndof();
142	+	}
143	+	else if(mu->BestTrk() != 0){
144	+	RChi2 = mu->BestTrk()->Chi2()/mu->BestTrk()->Ndof();
145	+	}
146		Bool_t idpass = kFALSE;
147		switch (fMuIDType) {
148	+	case kWMuId:
149	+	idpass = mu->BestTrk() != 0 &&
150	+	mu->BestTrk()->NHits() > 10 &&
151	+	RChi2 < 10.0 &&
152	+	(mu->NSegments() > 1 || mu->NMatches() > 1) &&
153	+	mu->BestTrk()->NPixelHits() > 0 &&
154	+	mu->Quality().Quality(MuonQuality::GlobalMuonPromptTight);
155	+	break;
156	+	case kZMuId:
157	+	idpass = mu->BestTrk() != 0 &&
158	+	mu->BestTrk()->NHits() > 10 &&
159	+	(mu->NSegments() > 1 || mu->NMatches() > 1) &&
160	+	mu->BestTrk()->NPixelHits() > 0 &&
161	+	mu->Quality().Quality(MuonQuality::GlobalMuonPromptTight);
162	+	break;
163		case kLoose:
164	<	idpass = fMuonTools->IsGood(mu, MuonTools::kTMOneStationLoose) &&
165	<	fMuonTools->IsGood(mu, MuonTools::kTM2DCompatibilityLoose);
164	>	idpass = mu->BestTrk() != 0 &&
165	>	mu->Quality().Quality(MuonQuality::TMOneStationLoose) &&
166	>	mu->Quality().Quality(MuonQuality::TM2DCompatibilityLoose) &&
167	>	mu->BestTrk()->NHits() > 10 &&
168	>	RChi2 < 10.0 &&
169	>	mu->Quality().Quality(MuonQuality::GlobalMuonPromptTight);
170		break;
171		case kTight:
172	<	idpass = fMuonTools->IsGood(mu, MuonTools::kTMOneStationTight) &&
173	<	fMuonTools->IsGood(mu, MuonTools::kTM2DCompatibilityTight);
172	>	idpass = mu->BestTrk() !=  0 &&
173	>	mu->Quality().Quality(MuonQuality::TMOneStationTight) &&
174	>	mu->Quality().Quality(MuonQuality::TM2DCompatibilityTight) &&
175	>	mu->BestTrk()->NHits() > 10 &&
176	>	RChi2 < 10.0 &&
177	>	mu->Quality().Quality(MuonQuality::GlobalMuonPromptTight);
178	>	break;
179	>	case kWWMuId:
180	>	idpass = mu->BestTrk() != 0 &&
181	>	mu->BestTrk()->NHits() > 10 &&
182	>	RChi2 < 10.0 &&
183	>	(mu->NSegments() > 1 || mu->NMatches() > 1) &&
184	>	mu->BestTrk()->NPixelHits() > 0 &&
185	>	mu->Quality().Quality(MuonQuality::GlobalMuonPromptTight) &&
186	>	mu->BestTrk()->PtErr()/mu->BestTrk()->Pt() < 0.1;
187	>	break;
188	>	case kNoId:
189	>	idpass = kTRUE;
190		break;
191		default:
192		break;
#	Line 88 | Line 195 | void MuonIDMod::Process()
195		if (!idpass)
196		continue;
197
198	<	Bool_t isopass = kFALSE;
198	>	Bool_t isocut = kFALSE;
199		switch (fMuIsoType) {
200		case kTrackCalo:
201	<	isopass = (mu->IsoR03SumPt() < fTrackIsolationCut) &&
201	>	isocut = (mu->IsoR03SumPt() < fTrackIsolationCut) &&
202		(mu->IsoR03EmEt() + mu->IsoR03HadEt() < fCaloIsolationCut);
203		break;
204		case kTrackCaloCombined:
205	<	isopass = (1.0 * mu->IsoR03SumPt() + 1.0 * mu->IsoR03EmEt() +
206	<	1.0 * mu->IsoR03HadEt() < fCombIsolationCut);
205	>	isocut = (1.0 * mu->IsoR03SumPt() +
206	>	1.0 * mu->IsoR03EmEt()  +
207	>	1.0 * mu->IsoR03HadEt() < fCombIsolationCut);
208		break;
209		case kTrackCaloSliding:
210		{
211	<	Double_t totalIso = 1.0 * mu->IsoR03SumPt() +
212	<	1.0 * mu->IsoR03EmEt() +
213	<	1.0 * mu->IsoR03HadEt();
214	<	if ((totalIso < (mu->Pt()-10.0)*5.0/15.0) ||
215	<	(totalIso < 5.0 && mu->Pt() > 25))
216	<	isopass = kTRUE;
217	<	}
211	>	//Double_t beta = IsolationTools::BetaM(fTracks, mu, fVertices->At(0), 0.0, 0.2, 0.3, 0.02);
212	>	//if(beta == 0) beta = 1.0;
213	>	const PileupEnergyDensity *rho =  fPileupEnergyDensity->At(0);
214	>	Double_t totalIso =  mu->IsoR03SumPt() + TMath::Max(mu->IsoR03EmEt() + mu->IsoR03HadEt() - rho->Rho() * TMath::Pi() * 0.3 * 0.3, 0.0);
215	>	if (totalIso < (mu->Pt()*fCombIsolationCut) )
216	>	isocut = kTRUE;
217	>
218	>	if     (fReverseIsoCut == kTRUE &&
219	>	isocut == kFALSE && totalIso < 10)
220	>	isocut = kTRUE;
221	>	else if(fReverseIsoCut == kTRUE)
222	>	isocut = kFALSE;
223	>	}
224	>	break;
225	>	case kTrackCaloSlidingNoCorrection:
226	>	{
227	>	Double_t totalIso =  1.0 * mu->IsoR03SumPt() +
228	>	1.0 * mu->IsoR03EmEt()  +
229	>	1.0 * mu->IsoR03HadEt();
230	>	if (totalIso < (mu->Pt()*fCombIsolationCut) )
231	>	isocut = kTRUE;
232	>	}
233	>	break;
234	>	case kPFIso:
235	>	{
236	>	Double_t beta = IsolationTools::BetaM(fTracks, mu, fVertices->At(0), 0.0, 0.2, 0.3, 0.02);
237	>	if(beta == 0) beta = 1.0;
238	>	Double_t totalIso =  IsolationTools::PFMuonIsolation(mu, fPFCandidates, fVertices->At(0), 0.2, 0.5, 0.3, 0.02, 0, beta, fNonIsolatedMuons, fNonIsolatedElectrons);
239	>	if (totalIso < (mu->Pt()*fCombIsolationCut) )
240	>	isocut = kTRUE;
241	>	}
242	>	break;
243	>	case kPFIsoNoL:
244	>	{
245	>	fNonIsolatedMuons     = GetObjThisEvt<MuonCol>(fNonIsolatedMuonsName);
246	>	fNonIsolatedElectrons = GetObjThisEvt<ElectronCol>(fNonIsolatedElectronsName);
247	>
248	>	Double_t beta = IsolationTools::BetaM(fTracks, mu, fVertices->At(0), 0.0, 0.2, 0.3, 0.02);
249	>	if(beta == 0) beta = 1.0;
250	>	Double_t totalIso =  IsolationTools::PFMuonIsolation(mu, fPFCandidates, fVertices->At(0), 0.2, 0.5, 0.3, 0.02, 3, beta, fNonIsolatedMuons, fNonIsolatedElectrons);
251	>	if (totalIso < (mu->Pt()*fCombIsolationCut) )
252	>	isocut = kTRUE;
253	>	}
254	>	break;
255	>	case kNoIso:
256	>	isocut = kTRUE;
257		break;
111	–	case kNoIso:
112	–	isopass = kTRUE;
113	–	break;
258		case kCustomIso:
259		default:
260		break;
261		}
262
263	<	if (!isopass)
263	>	if (isocut == kFALSE)
264		continue;
265
266	+	if (fApplyD0Cut) {
267	+	Bool_t passD0cut = MuonTools::PassD0Cut(mu, fVertices, fD0Cut);
268	+	if (!passD0cut)
269	+	continue;
270	+	}
271	+
272		// add good muon
273		CleanMuons->Add(mu);
274		}
275
276	+	// sort according to pt
277	+	CleanMuons->Sort();
278	+
279		// add objects for other modules to use
280		AddObjThisEvt(CleanMuons);
281		}
282
130	–
283		//--------------------------------------------------------------------------------------------------
284		void MuonIDMod::SlaveBegin()
285		{
286		// Run startup code on the computer (slave) doing the actual analysis. Here,
287		// we just request the muon collection branch.
288
289	<	ReqBranch(fMuonBranchName, fMuons);
290	<
291	<	fMuonTools = new MuonTools;
289	>	// In this case we cannot have a branch
290	>	if (fMuonIsoType.CompareTo("PFIsoNoL") != 0) {
291	>	ReqEventObject(fMuonBranchName, fMuons, kTRUE);
292	>	}
293	>	ReqEventObject(fTrackName, fTracks, kTRUE);
294	>	ReqEventObject(fPFCandidatesName, fPFCandidates, kTRUE);
295	>	if (fMuonIsoType.CompareTo("TrackCaloSliding") == 0) {
296	>	ReqEventObject(fPileupEnergyDensityName, fPileupEnergyDensity, kTRUE);
297	>	}
298
299	<	if (fMuonIDType.CompareTo("Tight") == 0)
299	>	if (fMuonIDType.CompareTo("WMuId") == 0)
300	>	fMuIDType = kWMuId;
301	>	else if (fMuonIDType.CompareTo("ZMuId") == 0)
302	>	fMuIDType = kZMuId;
303	>	else if (fMuonIDType.CompareTo("Tight") == 0)
304		fMuIDType = kTight;
305		else if (fMuonIDType.CompareTo("Loose") == 0)
306		fMuIDType = kLoose;
307	+	else if (fMuonIDType.CompareTo("WWMuId") == 0)
308	+	fMuIDType = kWWMuId;
309	+	else if (fMuonIDType.CompareTo("NoId") == 0)
310	+	fMuIDType = kNoId;
311		else if (fMuonIDType.CompareTo("Custom") == 0) {
312		fMuIDType = kCustomId;
313		SendError(kWarning, "SlaveBegin",
#	Line 159 | Line 325 | void MuonIDMod::SlaveBegin()
325		fMuIsoType = kTrackCaloCombined;
326		else if (fMuonIsoType.CompareTo("TrackCaloSliding") == 0)
327		fMuIsoType = kTrackCaloSliding;
328	+	else if (fMuonIsoType.CompareTo("TrackCaloSlidingNoCorrection") == 0)
329	+	fMuIsoType = kTrackCaloSlidingNoCorrection;
330	+	else if (fMuonIsoType.CompareTo("PFIso") == 0)
331	+	fMuIsoType = kPFIso;
332	+	else if (fMuonIsoType.CompareTo("PFIsoNoL") == 0)
333	+	fMuIsoType = kPFIsoNoL;
334		else if (fMuonIsoType.CompareTo("NoIso") == 0)
335		fMuIsoType = kNoIso;
336		else if (fMuonIsoType.CompareTo("Custom") == 0) {
#	Line 178 | Line 350 | void MuonIDMod::SlaveBegin()
350		fMuClassType = kGlobal;
351		else if (fMuonClassType.CompareTo("Standalone") == 0)
352		fMuClassType = kSta;
353	<	else if (fMuonClassType.CompareTo("TrackerOnly") == 0)
354	<	fMuClassType = kTrackerOnly;
353	>	else if (fMuonClassType.CompareTo("TrackerMuon") == 0)
354	>	fMuClassType = kTrackerMuon;
355	>	else if (fMuonClassType.CompareTo("CaloMuon") == 0)
356	>	fMuClassType = kCaloMuon;
357	>	else if (fMuonClassType.CompareTo("TrackerBased") == 0)
358	>	fMuClassType = kTrackerBased;
359		else {
360		SendError(kAbortAnalysis, "SlaveBegin",
361		"The specified muon class %s is not defined.",

Diff Legend

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