ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitPhysics/Utils/src/ElectronTools.cc

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.8 by ceballos, Fri May 28 16:31:43 2010 UTC vs.
Revision 1.11 by ceballos, Sun Jun 6 17:27:21 2010 UTC

#	Line 75 | Line 75 | Bool_t ElectronTools::PassCustomID(const
75		{0.025, 0.025, 0.025, 0.025, 0.012,  0.012,  0.012,  0.012}, //hovere
76		{0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03 }, //sigmaetaeta
77		{0.03,  0.03,  0.03,  0.03,  0.02,   0.02,   0.02,   0.02 }, //deltaphiin
78	<	{0.003, 0.003, 0.003, 0.003, 0.005,  0.005,  0.005,  0.005}, //deltaetain
78	>	{0.004, 0.004, 0.004, 0.004, 0.005,  0.005,  0.005,  0.005}, //deltaetain
79		{0.0,   0.0,   0.0,   0.0,   0.0,    0.0,    0.0,    0.0  }, //eoverp
80		{0.0,   0.0,   0,     0,     0,      0,      0,      0    }  //extra cuts fbrem and E_Over_P
81		};
#	Line 187 | Line 187 | Bool_t ElectronTools::PassCustomIso(cons
187		};
188
189		Double_t VBTFWorkingPoint70[4][2] = {
190	<	{0.06 , 0.03   },   //TrkIso
191	<	{0.06 , 0.03   },   //ECALIso
192	<	{0.05 , 0.015  },   //HCALIso
193	<	{0.05,  0.04   }   //Combined
190	>	{0.05 , 0.025  },   //TrkIso
191	>	{0.06 , 0.025  },   //ECALIso
192	>	{0.03 , 0.020  },   //HCALIso
193	>	{0.04,  0.030  }   //Combined
194		};
195
196		switch (isoType) {
#	Line 214 | Line 214 | Bool_t ElectronTools::PassCustomIso(cons
214		break;
215		}
216
217	<	Double_t trkIso  = ele->TrackIsolationDr03();
218	<	Double_t ecalIso = ele->EcalRecHitIsoDr03();
219	<	Double_t hcalIso = ele->HcalTowerSumEtDr03();
220	<	Double_t combinedIso = trkIso + ecalIso + hcalIso;
221	<	if(ele->IsEB()) combinedIso = trkIso + TMath::Max(ecalIso - 1.0, 0.0) + hcalIso;
217	>	Double_t trkIso  = ele->TrackIsolationDr03() / ele->Pt();
218	>	Double_t ecalIso = ele->EcalRecHitIsoDr03() / ele->Pt();
219	>	Double_t hcalIso = ele->HcalTowerSumEtDr03() / ele->Pt();
220	>	Double_t combinedIso = ele->TrackIsolationDr03() + ele->EcalRecHitIsoDr03() + ele->HcalTowerSumEtDr03();
221	>	if(ele->IsEB()) combinedIso = ele->TrackIsolationDr03() + TMath::Max(ele->EcalRecHitIsoDr03() - 1.0, 0.0) + ele->HcalTowerSumEtDr03();
222		combinedIso = combinedIso / ele->Pt();
223
224		Int_t eb = 1;
#	Line 366 | Line 366 | Bool_t ElectronTools::PassSpikeRemovalFi
366		return passSpikeRemovalFilter;
367		}
368
369	+	Bool_t ElectronTools::PassTriggerMatching(const Electron *ele, const TriggerObjectCol *trigobjs)
370	+	{
371	+
372	+	for (UInt_t i=0; i<trigobjs->GetEntries(); ++i) {
373	+	const TriggerObject *trigobj = trigobjs->At(i);
374	+	if (trigobj->TriggerType()==TriggerObject::TriggerCluster || trigobj->TriggerType()==TriggerObject::TriggerElectron) {
375	+	if (MathUtils::DeltaR(ele,trigobj)<0.3) {
376	+	return kTRUE;
377	+	}
378	+	}
379	+	}
380	+
381	+	return kFALSE;
382	+
383	+
384	+	}
385	+
386		//--------------------------------------------------------------------------------------------------
387		Int_t ElectronTools::Classify(const Electron *ele) {
388

Diff Legend

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