From: basanta <basanta@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 3 Jul 2008 15:13:40 +0000 (+0000)
Subject: photon trackno corrected to assign the PMD cell
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=dbb41941d1a3622b755e2144ffbfa8859e0f303c;p=u%2Fmrichter%2FAliRoot.git

photon trackno corrected to assign the PMD cell
---

diff --git a/PMD/AliPMDDigitizer.cxx b/PMD/AliPMDDigitizer.cxx
index 0e2a9fd839d..f8891e67130 100644
--- a/PMD/AliPMDDigitizer.cxx
+++ b/PMD/AliPMDDigitizer.cxx
@@ -325,6 +325,7 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
 		  statusOld   = -1;
 		}
 	      Int_t igstatus = 0;
+	      Int_t trnotemp = trackno;
 	      //------------------modified by Mriganka ----------------------
 	      if(ks==1||(imo = mparticle->GetFirstMother())<0 ){
 		vx = mparticle->Vx();
@@ -338,24 +339,31 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
 	      
 	      while(((imo = mparticle->GetFirstMother()) >= 0)&& 
 		    (ks = mparticle->GetStatusCode() <1) )
-		{
+	      {
 		  mparticle =  gAlice->GetMCApp()->Particle(imo);
 		  trackpid = mparticle->GetPdgCode();
 		  ks = mparticle->GetStatusCode();
 		  vx = mparticle->Vx();
 		  vy = mparticle->Vy();
 		  vz = mparticle->Vz();
+		  trnotemp = trackno;
+                  if(trackpid == 111)
+		  {
+		      trackno = trnotemp;
+		  }
+                  if(trackpid != 111)
+		  {
+		      trackno = imo;
+		  }
 		  
-		  trackno=imo;
-		  		
-		    }
+	      }
 	 
-	      if(trackpid==kGamma||trackpid==11||trackpid==-11||
-		 trackpid==kPi0)igstatus=1;
-	      mtrackpid=trackpid;
-	      mtrackno=trackno;
-	      trackpid=trackpidOld;
-	      trackno=tracknoOld;
+	      if(trackpid == kGamma || trackpid == 11 ||
+		 trackpid == -11 || trackpid == kPi0) igstatus = 1;
+	      mtrackpid = trackpid;
+	      mtrackno  = trackno;
+	      trackpid  = trackpidOld;
+	      trackno   = tracknoOld;
 	      
 	      //-----------------end of modification----------------
 	      xPos = fPMDHit->X();
@@ -393,7 +401,7 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
 		}
 
 	      AliDebug(2,Form("Zposition = %f Edeposition = %f",zPos,edep));
-	      //Float_t zposition = TMath::Abs(zPos);
+
 	      if (zPos < fZPos)
 		{
 		  // CPV