A new particle specie was created: a PHOTON HADRON

[u/mrichter/AliRoot.git] / PHOS / AliPHOSPIDv1.cxx

diff --git a/PHOS/AliPHOSPIDv1.cxx b/PHOS/AliPHOSPIDv1.cxx
index f4e934e96c920bcc4f000574c8077ede028907f1..ec95def00ab9f4ebe945aedce61d22540e30f66e 100644 (file)
--- a/PHOS/AliPHOSPIDv1.cxx
+++ b/PHOS/AliPHOSPIDv1.cxx
@@ -69,9 +69,17 @@ void  AliPHOSPIDv1::MakeParticles(TrackSegmentsList * trsl, RecParticlesList * r
       
       if( tracksegment->GetPpsdLow() == 0 )    // Neutral  
        type = kNEUTRAL ;   
-      else                           // Gamma
-       type = kGAMMA ;               
-    }
+      else {    // check the shower profile       
+       AliPHOSEmcRecPoint * recp = tracksegment->GetEmcRecPoint() ; 
+       Float_t * lambda = new Float_t[2]; 
+       recp->GetElipsAxis(lambda) ; 
+       if ( ( lambda[0] > fLambda1m && lambda[0] < fLambda1M ) && // shower profile cut
+            ( lambda[1] > fLambda2m && lambda[1] < fLambda2M ) )       
+         type = kGAMMA ;                      // a well identified photon 
+       else 
+         type = kGAMMAHADRON ;                // looks like a photon but is a hadron (most likely)  
+      }
+    } // Neutral
     else                            // Charged           
       type = kCHARGED ;