ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/MitPhysics/Utils/src/JetTools.cc (file contents):
Revision 1.8 by ceballos, Wed Jan 19 10:18:19 2011 UTC vs.
Revision 1.9 by ceballos, Thu Feb 17 14:09:07 2011 UTC

#	Line 300 | Line 300 | Double_t JetTools::MtHiggs(const Particl
300		return mtHiggs;
301		}
302
303	<	void JetTools::Alpha(Double_t AlphaVar[2], const TrackCol *tracks, Jet *jet, const VertexCol *vertices, Double_t  delta_z, Double_t delta_cone){
304	<	AlphaVar[0] = -1.0;
305	<	AlphaVar[1] = -1.0;
306	<	if(tracks->GetEntries() <= 0) return;
303	>	Double_t JetTools::Beta(const TrackCol *tracks, Jet *jet, const Vertex *vertex, Double_t  delta_z, Double_t delta_cone){
304	>
305	>	if(tracks->GetEntries() <= 0) return 1.0;
306
307		double Pt_jets_X = 0. ;
308		double Pt_jets_Y = 0. ;
#	Line 314 | Line 313 | void JetTools::Alpha(Double_t AlphaVar[2
313		if(MathUtils::DeltaR(tracks->At(i)->Mom(),jet->Mom()) < delta_cone){
314		Pt_jets_X_tot += tracks->At(i)->Px();
315		Pt_jets_Y_tot += tracks->At(i)->Py();
316	<	double pDz = TMath::Abs(tracks->At(i)->DzCorrected(*vertices->At(0)));
316	>	double pDz = TMath::Abs(tracks->At(i)->DzCorrected(*vertex));
317		if(pDz < delta_z){
318		Pt_jets_X += tracks->At(i)->Px();
319		Pt_jets_Y += tracks->At(i)->Py();
#	Line 322 | Line 321 | void JetTools::Alpha(Double_t AlphaVar[2
321		}
322		}
323
325	–	if(jet->Pt() > 0)
326	–	AlphaVar[0] = sqrt(Pt_jets_X*Pt_jets_X + Pt_jets_Y*Pt_jets_Y) / jet->Pt();
324		if(sqrt(Pt_jets_X_tot*Pt_jets_X_tot + Pt_jets_Y_tot*Pt_jets_Y_tot) > 0)
325	<	AlphaVar[1] = sqrt(Pt_jets_X*Pt_jets_X + Pt_jets_Y*Pt_jets_Y) / sqrt(Pt_jets_X_tot*Pt_jets_X_tot + Pt_jets_Y_tot*Pt_jets_Y_tot);
329	<	}
325	>	return sqrt(Pt_jets_X*Pt_jets_X + Pt_jets_Y*Pt_jets_Y) / sqrt(Pt_jets_X_tot*Pt_jets_X_tot + Pt_jets_Y_tot*Pt_jets_Y_tot);
326
327	+	return 1.0;
328	+	}
329
332	–	void JetTools::Alpha(Double_t AlphaVar[2], const PFJet *jet, const VertexCol *vertices, Double_t  delta_z){
333	–	AlphaVar[0] = -1.0;
334	–	AlphaVar[1] = -1.0;
330
331	+	Double_t JetTools::Beta(const PFJet *jet, const Vertex *vertex, Double_t  delta_z){
332		double Pt_jets_X = 0. ;
333		double Pt_jets_Y = 0. ;
334		double Pt_jets_X_tot = 0. ;
#	Line 342 | Line 338 | void JetTools::Alpha(Double_t AlphaVar[2
338		if(jet->PFCand(i)->BestTrk()){
339		Pt_jets_X_tot += jet->PFCand(i)->BestTrk()->Px();
340		Pt_jets_Y_tot += jet->PFCand(i)->BestTrk()->Py();
341	<	double pDz = TMath::Abs(jet->PFCand(i)->BestTrk()->DzCorrected(*vertices->At(0)));
341	>	double pDz = TMath::Abs(jet->PFCand(i)->BestTrk()->DzCorrected(*vertex));
342		if(pDz < delta_z){
343		Pt_jets_X += jet->PFCand(i)->BestTrk()->Px();
344	<	Pt_jets_Y += jet->PFCand(i)->BestTrk()->Py();
344	>	Pt_jets_Y += jet->PFCand(i)->BestTrk()->Py();
345		}
346		}
347		}
348
353	–	if(jet->Pt() > 0)
354	–	AlphaVar[0] = sqrt(Pt_jets_X*Pt_jets_X + Pt_jets_Y*Pt_jets_Y) / jet->Pt();
349		if(sqrt(Pt_jets_X_tot*Pt_jets_X_tot + Pt_jets_Y_tot*Pt_jets_Y_tot) > 0)
350	<	AlphaVar[1] = sqrt(Pt_jets_X*Pt_jets_X + Pt_jets_Y*Pt_jets_Y) / sqrt(Pt_jets_X_tot*Pt_jets_X_tot + Pt_jets_Y_tot*Pt_jets_Y_tot);
350	>	return sqrt(Pt_jets_X*Pt_jets_X + Pt_jets_Y*Pt_jets_Y) / sqrt(Pt_jets_X_tot*Pt_jets_X_tot + Pt_jets_Y_tot*Pt_jets_Y_tot);
351	>
352	>	return 1.0;
353		}

Diff Legend

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