Fill ntuple only on request; more mass histos (Francesco, Renu)

[u/mrichter/AliRoot.git] / ZDC / AliZDCHit.cxx

diff --git a/ZDC/AliZDCHit.cxx b/ZDC/AliZDCHit.cxx
index 204709f61f9a275e955c9cbb17e8a920bbe31e67..551c0875f00abf93c84a0ebd1814ecd4588ca705 100644 (file)
--- a/ZDC/AliZDCHit.cxx
+++ b/ZDC/AliZDCHit.cxx
@@ -34,7 +34,10 @@ AliZDCHit::AliZDCHit() :
   fSFlag(0),
   fLightPMQ(0.),
   fLightPMC(0.),
-  fEnergy(0.) 
+  fEnergy(0.), 
+  fPDGCode(0),
+  fMotherPDGCode(0),
+  fTrackTOF(0.)
 
 {
   //
@@ -52,7 +55,10 @@ AliZDCHit::AliZDCHit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits) :
   fSFlag(hits[6]),
   fLightPMQ(hits[7]),
   fLightPMC(hits[8]),
-  fEnergy(hits[9]) 
+  fEnergy(hits[9]), 
+  fPDGCode((Int_t) hits[10]),
+  fMotherPDGCode((Int_t) hits[11]),
+  fTrackTOF(hits[12])
 
 {
   //
@@ -67,20 +73,23 @@ AliZDCHit::AliZDCHit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits) :
   
 //_____________________________________________________________________________
 AliZDCHit::AliZDCHit(const AliZDCHit &oldhit) :
-  AliHit(0,oldhit.GetTrack())
+  AliHit(0,oldhit.GetTrack()),
+  fPrimKinEn(oldhit.GetPrimKinEn()),
+  fXImpact(oldhit.GetXImpact()),  
+  fYImpact(oldhit.GetYImpact()),  
+  fSFlag(oldhit.GetSFlag()),
+  fLightPMQ(oldhit.GetLightPMQ()), 
+  fLightPMC(oldhit.GetLightPMC()),
+  fEnergy(oldhit.GetEnergy()),
+  fPDGCode(oldhit.GetPDGCode()),
+  fMotherPDGCode(oldhit.GetMotherPDGCode()),
+  fTrackTOF(oldhit.GetTrackTOF())
 {
   // Copy constructor
   fX = oldhit.X();
   fY = oldhit.Y();
   fZ = oldhit.Z();
   for(Int_t i=0; i<2; i++) fVolume[i] = oldhit.GetVolume(i);
-  fPrimKinEn = oldhit.GetPrimKinEn();
-  fXImpact = oldhit.GetXImpact();  
-  fYImpact = oldhit.GetYImpact();  
-  fSFlag = oldhit.GetSFlag();    
-  fLightPMQ = oldhit.GetLightPMQ(); 
-  fLightPMC = oldhit.GetLightPMC(); 
-  fEnergy = oldhit.GetEnergy();   
 }
   
   
@@ -88,9 +97,10 @@ AliZDCHit::AliZDCHit(const AliZDCHit &oldhit) :
 void AliZDCHit::Print(Option_t *) const 
 {
    // Print method
-   printf(" -> HIT: vol[0] =  %d vol[1] =  %d Track: %d \n" 
-         "  Primary E = %f, Ximpact = %f, Yimpact = %f, SFlag = %f\n"
-          "  PMQLight = %f, PMCLight = %f,  Deposited E = %f\n ", 
-          fVolume[0],fVolume[1],fTrack,fPrimKinEn,fXImpact,fYImpact,
-          fSFlag,fLightPMQ,fLightPMC,fEnergy);
+   printf("\t AliZDCHit: track %d PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV SFlag = %1.0f\n" 
+         "\t DETECTOR (%d, %d)  (X, Y)_impact (%f, %f) cm\n"
+          "\t PMQLight %1.0f, PMCLight %1.0f,  E_dep %1.2f\n ", 
+          fTrack,fPDGCode,fTrackTOF,fPrimKinEn,fSFlag,
+         fVolume[0],fVolume[1],fXImpact,fYImpact,
+          fLightPMQ,fLightPMC,fEnergy);
 }