ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitHzz4l/Angles/src/ZGAngles.cc

Comparing UserCode/MitHzz4l/Angles/src/ZGAngles.cc (file contents):
Revision 1.1 by khahn, Wed Aug 8 08:50:03 2012 UTC vs.
Revision 1.2 by khahn, Mon Aug 13 19:13:34 2012 UTC

#	Line 19 | Line 19 | using namespace std;
19		// fill angle struct from ZZ-type GenInfo objects
20		//
21		void
22	<	fillZGAngles( GenInfoStruct &geninfo, ZGAngles &a )
22	>	fillZGAngles( GenInfoStruct &geninfo, ZGAngles &a, bool debug )
23		// --------------------------------------------------------------------------------------------------
24		{
25		ZGLabVectors l;
#	Line 29 | Line 29 | fillZGAngles( GenInfoStruct &geninfo, ZG
29
30		if( geninfo.id_1_a < 0 ) {
31		l.veclp.SetPtEtaPhiM( geninfo.pt_1_a, geninfo.eta_1_a, geninfo.phi_1_a,
32	<	(geninfo.id_1_a == 11 ? ELECTRON_MASS : MUON_MASS)  );
32	>	(abs(geninfo.id_1_a) == 11 ? ELECTRON_MASS : MUON_MASS)  );
33		l.veclm.SetPtEtaPhiM( geninfo.pt_2_a, geninfo.eta_2_a, geninfo.phi_2_a,
34	<	(geninfo.id_2_a == 11 ? ELECTRON_MASS : MUON_MASS)  );
34	>	(abs(geninfo.id_2_a) == 11 ? ELECTRON_MASS : MUON_MASS)  );
35		} else {
36		l.veclp.SetPtEtaPhiM( geninfo.pt_2_a, geninfo.eta_2_a, geninfo.phi_2_a,
37	<	(geninfo.id_2_a == 11 ? ELECTRON_MASS : MUON_MASS)  );
37	>	(abs(geninfo.id_2_a) == 11 ? ELECTRON_MASS : MUON_MASS)  );
38		l.veclm.SetPtEtaPhiM( geninfo.pt_1_a, geninfo.eta_1_a, geninfo.phi_1_a,
39	<	(geninfo.id_1_a == 11 ? ELECTRON_MASS : MUON_MASS)  );
39	>	(abs(geninfo.id_1_a) == 11 ? ELECTRON_MASS : MUON_MASS)  );
40		}
41
42
#	Line 45 | Line 45 | fillZGAngles( GenInfoStruct &geninfo, ZG
45		l.vecz.SetPtEtaPhiM( geninfo.vpt_a, geninfo.veta_a, geninfo.vphi_a, geninfo.vmass_a);
46		l.veczg = (l.vecz+l.vecg);
47
48	<	getZGAngles( l, a );
48	>	getZGAngles( l, a, debug );
49		}
50
51		// --------------------------------------------------------------------------------------------------
#	Line 261 | Line 261 | getZGAngles( ZGLabVectors &l, ZGAngles &
261		a.ptl2  = l.veclm.Pt();
262		a.etal2 = l.veclm.Eta();
263		a.ptz   = (l.veclp+l.veclm).Pt();
264	<	a.etaz = (l.veclp+l.veclm).Eta();
264	>	a.etaz  = (l.veclp+l.veclm).Eta();
265	>	a.mzg   = l.veczg.M();
266	>	a.mz    = l.vecz.M();
267		};

Diff Legend

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