Changeset 409 in svn

Timestamp:

May 20, 2009, 10:48:41 AM (16 years ago)

Author:

Xavier Rouby

Message:

new class TRootGenFwdParticle; modifs in ZDC and RomanPot classes. New leaves in the branches

Location:

trunk/Utilities/ExRootAnalysis/src

Files:

• BlockClasses.cc (modified) (4 diffs)
• BlockClassesLinkDef.h (modified) (1 diff)

trunk/Utilities/ExRootAnalysis/src/BlockClasses.cc

@@ -36 +36 @@
 
 TCompare *TRootGenParticle::fgCompare = 0;
+TCompare *TRootGenFwdParticle::fgCompare = 0;
 TCompare *TRootC::GenParticle::fgCompare = 0;
 TCompare *TRootElectron::fgCompare = TComparePT<TRootElectron>::Instance();
@@ -63 +64 @@
 }
 
+void TRootGenFwdParticle::Set(const TLorentzVector& momentum) {
+        genE   = momentum.E();
+        genPx  = momentum.Px();
+        genPy  = momentum.Py();
+        genPz  = momentum.Pz();
+        genPT  = momentum.Pt();
+        genEta = momentum.Eta();
+        genPhi = momentum.Phi();
+}
+
+
 void TRootTauJet::Set(const TLorentzVector& momentum) {
 
@@ -87 +99 @@
 }
 
+void TRootGenFwdParticle::Set(const float px, const float py, const float pz, const float e) {
+        TLorentzVector toFill;
+        toFill.SetPxPyPzE(px,py,pz,e);
+        genE   = e;
+        genPx  = px;
+        genPy  = py;
+        genPz  = pz;
+        genPT  = toFill.Pt();
+        genEta = toFill.Eta();
+        genPhi = toFill.Phi();
+}
+
+
 void TRootParticle::SetEtaPhiEET(const float eta, const float phi, const float e, const float et) {
    Eta = eta; Phi = phi; E = e; PT = et;
    //float theta = 2* atan(exp(-Eta);
-   Px = PT*cos(Phi);
-   Py = PT*sin(Phi);
+   Px = PT*cos(Phi);     Py = PT*sin(Phi);
+}
+
+void TRootGenFwdParticle::SetEtaPhiEET(const float eta, const float phi, const float e, const float et) {
+   genEta = eta; genPhi = phi; genE = e; genPT = et;
+   genPx = genPT*cos(genPhi);  genPy = genPT*sin(genPhi);
+}
+
+
+TRootGenFwdParticle::TRootGenFwdParticle(const TRootGenFwdParticle& p): 
+        genE(p.genE), genPx(p.genPx), genPy(p.genPy), genPz(p.genPz), genPT(p.genPT), genEta(p.genEta), genPhi(p.genPhi) {}
+
+TRootGenFwdParticle& TRootGenFwdParticle::operator=(const TRootGenFwdParticle& p){
+        genE=p.genE;    genPx=p.genPx;   genPy=p.genPy;   genPz=p.genPz;   
+        genPT=p.genPT;  genEta=p.genEta; genPhi=p.genPhi;
 }
 
@@ -198 +236 @@
 }
 
-/*TRootZdcHits::TRootZdcHits(const float e, const float t, const int side, const bool had) : E(e), T(t), side(side), hadronic_hit(had) {}
+TRootZdcHits::TRootZdcHits(const float e, const float t, const int side, const bool had) : E(e), T(t), side(side), hadronic_hit(had) {}
 TRootZdcHits::TRootZdcHits(const TRootZdcHits& zdc) : E(zdc.E), T(zdc.T), side(zdc.side),  hadronic_hit(zdc.hadronic_hit) {}
-TRootZdcHits::TRootZdcHits(const TRootZdcHits& zdc) : T(zdc.T), side(zdc.side) {}
 TRootZdcHits& TRootZdcHits::operator=(const TRootZdcHits& zdc) {
         if(this==&zdc) return *this;
-        //E = zdc.E; T = zdc.T;  side = zdc.side; hadronic_hit = zdc.hadronic_hit;
-        T = zdc.T;  side = zdc.side; 
+        E = zdc.E; T = zdc.T;  side = zdc.side; hadronic_hit = zdc.hadronic_hit;
   return *this;
 }
-*/
+
 
 void TRootGenParticle:: print(){

trunk/Utilities/ExRootAnalysis/src/BlockClassesLinkDef.h

@@ -42 +42 @@
 #pragma link C++ class TRootSelectorInfo;
 #pragma link C++ class TRootGenParticle;
+#pragma link C++ class TRootGenFwdParticle;
 #pragma link C++ class TRootC::GenParticle;
 #pragma link C++ class TRootElectron;