Fixing a bug
authorcoppedis <Chiara.Oppedisano@cern.ch>
Tue, 3 Feb 2015 16:08:11 +0000 (17:08 +0100)
committercoppedis <Chiara.Oppedisano@cern.ch>
Tue, 3 Feb 2015 16:08:11 +0000 (17:08 +0100)
ZDC/ZDCsim/AliZDC.cxx
ZDC/ZDCsim/AliZDCHit.cxx
ZDC/ZDCsim/AliZDCHit.h

index 4a497b4..c30bb01 100644 (file)
@@ -178,7 +178,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
   
   if(fNhits==0){
       // First hit -> setting flag for primary or secondary particle
-      /*TParticle * p = gAlice->GetMCApp()->Particle(track);
+      TParticle * p = gAlice->GetMCApp()->Particle(track);
       Int_t imo = p->GetFirstMother();
       //
       if(track != imo){
@@ -186,7 +186,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
       }
       else if(track == imo){
         newquad->SetSFlag(0);  // PRIMARY particle entering the ZDC
-      }*/
+      }
       //  
       sFlag     = newquad->GetSFlag();
       primKinEn  = newquad->GetPrimKinEn();
@@ -196,6 +196,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
       motPDGcode = newquad->GetMotherPDGCode();
       trackTime  = newquad->GetTrackTOF();
       trackEta   = newquad->GetTrackEta();
+
    }
    else{       
       newquad->SetPrimKinEn(primKinEn);
@@ -207,7 +208,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
       newquad->SetTrackTOF(trackTime);
       newquad->SetTrackEta(trackEta);
    }
+  
   Int_t j;
   for(j=0; j<fNhits; j++){
     // If hits are equal (same track, same volume), sum them.
@@ -222,6 +223,7 @@ void AliZDC::AddHit(Int_t track, Int_t *vol, Float_t *hits)
        return;
      } 
   }
+  //printf( "PDG from hits[10] = %f\n", hits[10]);
 
     //Otherwise create a new hit
     new(lhits[fNhits]) AliZDCHit(*newquad);
index e641968..d51ce58 100644 (file)
@@ -125,7 +125,7 @@ AliZDCHit &AliZDCHit::operator= (const AliZDCHit &hit)
 void AliZDCHit::Print(Option_t *) const 
 {
    // Print method
-   printf("\t **** AliZDCHit: track %d eta %f PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV \n" 
+   printf(" **** AliZDCHit: track %d eta %f PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV \n" 
          "\t DETECTOR (%d, %d)  (X, Y)_impact (%1.2f, %1.2f) cm  Secflag %1.0f\n"
           "\t PMQLight %1.0f, PMCLight %1.0f,  E_dep %1.2f\n ", 
           fTrack,fTrackEta,fPDGCode,fTrackTOF,fPrimKinEn,
index 74eeca3..8118e3f 100644 (file)
@@ -35,7 +35,7 @@ public:
   // Setters 
   virtual void SetVolume(Int_t i, Int_t val) {fVolume[i]=val;} 
   virtual void SetPDGCode(Int_t code)        {fPDGCode=code;}
-  virtual void SetMotherPDGCode(Int_t code)  {fPDGCode=code;}
+  virtual void SetMotherPDGCode(Int_t code)  {fMotherPDGCode=code;}
   virtual void SetLightPMQ(Float_t value)    {fLightPMQ=value;}
   virtual void SetLightPMC(Float_t value)    {fLightPMC=value;}
   virtual void SetSFlag(Float_t value)       {fSFlag=value;}