Adding eta info at hit level
authorcoppedis <oppedisa@alz112xl.(none)>
Fri, 28 Feb 2014 11:54:42 +0000 (12:54 +0100)
committercoppedis <oppedisa@alz112xl.(none)>
Fri, 28 Feb 2014 11:55:57 +0000 (12:55 +0100)
ZDC/AliZDC.cxx
ZDC/AliZDCHit.cxx
ZDC/AliZDCHit.h
ZDC/AliZDCv3.cxx
ZDC/AliZDCv4.cxx

index 73e2c4a..662ace9 100644 (file)
@@ -160,7 +160,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
   //
   //           Add a ZDC hit to the hit list.
   
-  static Float_t trackTime=0., primKinEn=0., xImpact=0., yImpact=0., sFlag=0.;
+  static Float_t trackTime=0., trackEta=0., primKinEn=0., xImpact=0., yImpact=0., sFlag=0.;
   static Int_t   pcPDGcode, motPDGcode;
 
   AliZDCHit *newquad, *curprimquad;
@@ -186,6 +186,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
       pcPDGcode         = newquad->GetPDGCode();
       motPDGcode = newquad->GetMotherPDGCode();
       trackTime  = newquad->GetTrackTOF();
+      trackEta   = newquad->GetTrackEta();
    }
    else{       
       newquad->SetPrimKinEn(primKinEn);
@@ -195,6 +196,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
       newquad->SetPDGCode(pcPDGcode);
       newquad->SetMotherPDGCode(motPDGcode);
       newquad->SetTrackTOF(trackTime);
+      newquad->SetTrackEta(trackEta);
    }
  
   Int_t j;
index 1978424..ae07219 100644 (file)
@@ -37,7 +37,8 @@ AliZDCHit::AliZDCHit() :
   fEnergy(0.), 
   fPDGCode(0),
   fMotherPDGCode(0),
-  fTrackTOF(0.)
+  fTrackTOF(0.),
+  fTrackEta(0.)
 
 {
   //
@@ -58,7 +59,8 @@ AliZDCHit::AliZDCHit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits) :
   fEnergy(hits[9]), 
   fPDGCode((Int_t) hits[10]),
   fMotherPDGCode((Int_t) hits[11]),
-  fTrackTOF(hits[12])
+  fTrackTOF(hits[12]),
+  fTrackEta(hits[13])
 
 {
   //
@@ -83,7 +85,8 @@ AliZDCHit::AliZDCHit(const AliZDCHit &oldhit) :
   fEnergy(oldhit.GetEnergy()),
   fPDGCode(oldhit.GetPDGCode()),
   fMotherPDGCode(oldhit.GetMotherPDGCode()),
-  fTrackTOF(oldhit.GetTrackTOF())
+  fTrackTOF(oldhit.GetTrackTOF()),
+  fTrackEta(oldhit.GetTrackEta())
 {
   // Copy constructor
   fX = oldhit.X();
@@ -108,6 +111,7 @@ AliZDCHit &AliZDCHit::operator= (const AliZDCHit &hit)
   fPDGCode = hit.GetPDGCode();
   fMotherPDGCode = hit.GetMotherPDGCode();
   fTrackTOF = hit.GetTrackTOF();
+  fTrackEta = hit.GetTrackEta();
 
   fX = hit.X();
   fY = hit.Y();
@@ -122,10 +126,10 @@ AliZDCHit &AliZDCHit::operator= (const AliZDCHit &hit)
 void AliZDCHit::Print(Option_t *) const 
 {
    // Print method
-   printf("\t AliZDCHit: track %d PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV SFlag = %1.0f\n" 
+   printf("\t AliZDCHit: track %d eta %f PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV \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,
+          fTrack,fTrackEta,fPDGCode,fTrackTOF,fPrimKinEn,
          fVolume[0],fVolume[1],fXImpact,fYImpact,
           fLightPMQ,fLightPMC,fEnergy);
 }
index b4cd99e..74eeca3 100644 (file)
@@ -30,6 +30,7 @@ public:
   virtual Float_t GetLightPMC() const      {return fLightPMC;}
   virtual Float_t GetEnergy() const        {return fEnergy;}
   virtual Float_t GetTrackTOF() const      {return fTrackTOF;}
+  virtual Float_t GetTrackEta() const      {return fTrackEta;}
 
   // Setters 
   virtual void SetVolume(Int_t i, Int_t val) {fVolume[i]=val;} 
@@ -42,6 +43,7 @@ public:
   virtual void SetXImpact(Float_t value)     {fXImpact=value;}
   virtual void SetYImpact(Float_t value)     {fYImpact=value;}
   virtual void SetTrackTOF(Float_t value)    {fTrackTOF=value;}
+  virtual void SetTrackEta(Float_t value)    {fTrackEta=value;}
 
   // Operators
   Int_t operator == (AliZDCHit &quad){
@@ -74,8 +76,9 @@ protected:
   Int_t   fPDGCode;      //PDG code of particle in the ZDC
   Int_t   fMotherPDGCode; //PDG code of particle in the ZDC
   Float_t fTrackTOF;     //Track time in ns
+  Float_t fTrackEta;     //Track eta
 
-  ClassDef(AliZDCHit,4)  // Hits for the Zero Degree Calorimeters
+  ClassDef(AliZDCHit,5)  // Hits for the Zero Degree Calorimeters
 };
  
 #endif
index 6f226ab..949d02b 100644 (file)
@@ -2252,12 +2252,12 @@ void AliZDCv3::StepManager()
   // Routine called at every step in the Zero Degree Calorimeters
   //
   Int_t   j, vol[2]={0,0}, ibeta=0, ialfa=0, ibe=0, nphe=0;
-  Float_t hits[13], x[3], xdet[3]={999.,999.,999.}, um[3], ud[3];
+  Float_t hits[14], x[3], xdet[3]={999.,999.,999.}, um[3], ud[3];
   Float_t destep=0., be=0., out=0.;
   Double_t s[3], p[4];
   const char *knamed;
   //
-  for(j=0;j<13;j++) hits[j]=-999.;
+  for(j=0;j<14;j++) hits[j]=-999.;
   //
   // --- This part is for no shower developement in beam pipe, TDI, VColl
   // If particle interacts with beam pipe, TDI, VColl -> return
@@ -2464,6 +2464,7 @@ void AliZDCv3::StepManager()
       //
       hits[12] = 1.0e09*gMC->TrackTime(); // in ns!
       //printf("\t TrackTime = %f\n", hits[12]);
+      hits[13] = part->Eta();
 
       AddHit(curTrackN, vol, hits);
 
index d6c1296..e5407fa 100644 (file)
@@ -2254,12 +2254,12 @@ void AliZDCv4::StepManager()
   // Routine called at every step in the Zero Degree Calorimeters
   //
   Int_t   j, vol[2]={0,0}, ibeta=0, ialfa=0, ibe=0, nphe=0;
-  Float_t hits[13], x[3], xdet[3]={999.,999.,999.}, um[3], ud[3];
+  Float_t hits[14], x[3], xdet[3]={999.,999.,999.}, um[3], ud[3];
   Float_t destep=0., be=0., out=0.;
   Double_t s[3], p[4];
   const char *knamed;
   //
-  for(j=0;j<13;j++) hits[j]=-999.;
+  for(j=0;j<14;j++) hits[j]=-999.;
   //
   // --- This part is for no shower developement in beam pipe, TDI, VColl
   // If particle interacts with beam pipe, TDI, VColl -> return
@@ -2466,6 +2466,7 @@ void AliZDCv4::StepManager()
       //
       hits[12] = 1.0e09*gMC->TrackTime(); // in ns!
       //printf("\t TrackTime = %f\n", hits[12]);
+      hits[13] = part->Eta();
 
       AddHit(curTrackN, vol, hits);