Fork me on GitHub

Ignore:
Timestamp:
Mar 3, 2009, 12:48:15 AM (16 years ago)
Author:
Xavier Rouby
Message:

tracks and Taujets have charge

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Utilities/ExRootAnalysis/src/BlockClasses.cc

    r278 r290  
    126126}
    127127
    128     float Eta, Phi; // (eta,phi) at the beginning of the track
     128/*    float Eta, Phi; // (eta,phi) at the beginning of the track
    129129    float EtaOuter, PhiOuter; // (eta,phi) at the end of the track
    130     float PT, E, Px, Py, Pz;
     130    float PT, E, Px, Py, Pz; */
    131131
    132132TRootTracks::TRootTracks() :
    133         TRootParticle(),
    134         EtaOuter(UNDEFINED), PhiOuter(UNDEFINED) {}
    135 /*        Eta = UNDEFINED;
    136         Phi = UNDEFINED;
    137         PT = UNDEFINED;
    138         E = UNDEFINED;
    139         Px = UNDEFINED;
    140         Py = UNDEFINED;
    141         Pz = UNDEFINED;
    142 
    143 }*/
     133        Eta(UNDEFINED), Phi(UNDEFINED), EtaOuter(UNDEFINED), PhiOuter(UNDEFINED),
     134        PT(UNDEFINED), E(UNDEFINED), Px(UNDEFINED), Py(UNDEFINED), Pz(UNDEFINED),
     135        Charge(UNDEFINED) {}
    144136
    145137TRootTracks::TRootTracks(const float inEta, const float inPhi, const float outEta, const float outPhi, const float pt) :
    146         EtaOuter(outEta), PhiOuter(outPhi) {
    147         Eta = inEta;
    148         Phi = inPhi;
    149         PT = pt;
     138        Eta(inEta), Phi(inPhi), EtaOuter(outEta), PhiOuter(outPhi),
     139        PT(pt) {
    150140        TLorentzVector p; p.SetPtEtaPhiE(PT,Eta,Phi,PT*cosh(Eta));
    151141        Px = p.Px(); Py = p.Py(); Pz = p.Pz();  E=p.E();
     142        Charge = UNDEFINED;
    152143}
    153144TRootTracks::TRootTracks(const TRootTracks& track) :
    154         EtaOuter(track.EtaOuter), PhiOuter(track.PhiOuter)
    155         {
    156           Eta = track.Eta;
    157           Phi = track.Phi;
    158           PT = track.PT;
    159           E = track.E;
    160           Px = track.Px;
    161           Py = track.Py;
    162           Pz = track.Pz;
    163         }
     145        Eta(track.Eta), Phi(track.Phi),
     146        EtaOuter(track.EtaOuter), PhiOuter(track.PhiOuter),
     147        PT(track.PT), E(track.E), Px(track.Px), Py(track.Py), Pz(track.Pz), Charge(track.Charge) {}
    164148
    165149TRootTracks& TRootTracks::operator=(const TRootTracks& track) {
     
    168152        EtaOuter = track.EtaOuter;   PhiOuter = track.PhiOuter;
    169153        PT = track.PT; E = track.E;  Px = track.Px;  Py= track.Py; Pz =track.Pz;
     154        Charge = track.Charge;
    170155  return *this;
    171156}
    172157
    173 void TRootTracks::Set(const float inEta, const float inPhi, const float outEta, const float outPhi, const float pt) {
     158void TRootTracks::Set(const float inEta, const float inPhi, const float outEta, const float outPhi, const float pt, const float charge) {
    174159        Eta= inEta; Phi=inPhi;  EtaOuter=outEta; PhiOuter=outPhi; PT=pt;
    175160        TLorentzVector p; p.SetPtEtaPhiE(PT,Eta,Phi,PT*cosh(Eta));
    176161        Px = p.Px(); Py = p.Py(); Pz = p.Pz();  E=p.E();
     162        Charge = charge;
    177163}
    178164
Note: See TracChangeset for help on using the changeset viewer.