ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitHzz4l/Angles/src/ZGAngles.cc
(Generate patch)

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.3 by dkralph, Mon Dec 17 12:35:22 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