Full PID delegated to AliPHOSPID
authorschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Feb 2000 16:44:58 +0000 (16:44 +0000)
committerschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Feb 2000 16:44:58 +0000 (16:44 +0000)
PHOS/AliPHOSAnalyze.cxx
PHOS/AliPHOSPID.h
PHOS/AliPHOSPIDv1.cxx
PHOS/AliPHOSPIDv1.h
PHOS/AliPHOSRecParticle.cxx
PHOS/AliPHOSReconstructioner.cxx

index ca704bb..3f2d8ff 100644 (file)
@@ -249,7 +249,7 @@ void AliPHOSAnalyze::AnalyzeOneEvent(Int_t evt)
                      //fhNeutralEMPositionY->Fill(recpart. ) ; 
                      //cout << "NEUTRAL EM" << endl;
                      break ;
-                   case kCHARGEDHADRON :
+                   case kCHARGEDHADRON:
                      fhChargedHadronEnergy->Fill(recparticle->Energy() ) ; 
                      //fhChargedHadronPositionX->Fill(recpart. ) ;
                      //fhChargedHadronPositionY->Fill(recpart. ) ; 
index 2f99a52..f538eb0 100644 (file)
@@ -36,7 +36,7 @@ public:
   AliPHOSPID() ;          // ctor            
   virtual ~AliPHOSPID() ; // dtor
 
-  virtual void GetParticleType(TrackSegmentsList * trsl, RecParticlesList * rpl) {} ; 
+  virtual void MakeParticles(TrackSegmentsList * trsl, RecParticlesList * rpl) {} ; 
   virtual void SetShowerProfileCuts(Float_t, Float_t, Float_t, Float_t) {} ; 
   virtual void SetDispersionCutOff(Float_t ) {}    
 
index 4454fce..f4e934e 100644 (file)
@@ -49,7 +49,7 @@ ClassImp( AliPHOSPIDv1)
 
 
 //____________________________________________________________________________
-void  AliPHOSPIDv1::GetParticleType(TrackSegmentsList * trsl, RecParticlesList * rpl)
+void  AliPHOSPIDv1::MakeParticles(TrackSegmentsList * trsl, RecParticlesList * rpl)
 {
   // main function, does the job
 
@@ -57,11 +57,26 @@ void  AliPHOSPIDv1::GetParticleType(TrackSegmentsList * trsl, RecParticlesList *
   AliPHOSTrackSegment * tracksegment ; 
   Int_t index = 0 ; 
   AliPHOSRecParticle * rp ; 
+  Int_t type ; 
+
   while ( (tracksegment = (AliPHOSTrackSegment *)next()) ) {
     new( (*rpl)[index] ) AliPHOSRecParticle(tracksegment) ;
     rp = (AliPHOSRecParticle *)(*rpl)[index] ; 
-    Int_t type =  rp->GetType() ;  
-    if ( type == kNEUTRAL ) { // resolve neutral baryon from photon
+
+    // try to figure out the type of particle:
+    //    1. just looking at the PPSD information 
+    if( tracksegment->GetPpsdUp() == 0 ) {     // Neutral
+      
+      if( tracksegment->GetPpsdLow() == 0 )    // Neutral  
+       type = kNEUTRAL ;   
+      else                           // Gamma
+       type = kGAMMA ;               
+    }
+    else                            // Charged           
+      type = kCHARGED ;   
+
+    //   2. from the shower profile analysis
+    if ( type == kNEUTRAL ) { 
       AliPHOSEmcRecPoint * recp = tracksegment->GetEmcRecPoint() ; 
       Float_t * lambda = new Float_t[2]; 
       recp->GetElipsAxis(lambda) ; 
@@ -72,6 +87,8 @@ void  AliPHOSPIDv1::GetParticleType(TrackSegmentsList * trsl, RecParticlesList *
        type = kNEUTRALHADRON ; 
       delete lambda ; 
    }
+
+    //   3. from the shower dispersion 
     if (type == kCHARGED) { 
       if( tracksegment->GetEmcRecPoint()->GetDispersion() > fCutOnDispersion)  // shower dispersion cut
        type = kCHARGEDHADRON ;
index 4ddb24d..e983074 100644 (file)
@@ -28,7 +28,7 @@ public:
   AliPHOSPIDv1() ;                     
   virtual ~ AliPHOSPIDv1() ; // dtor
 
-  virtual void GetParticleType(TrackSegmentsList * trsl, RecParticlesList * rpl ) ; // does the job
+  virtual void MakeParticles(TrackSegmentsList * trsl, RecParticlesList * rpl ) ; // does the job
   void Print() ; 
   virtual void SetDispersionCutOff(Float_t Dcut) {fCutOnDispersion = Dcut ; }    
   virtual void SetShowerProfileCuts(Float_t l1m, Float_t l1M, Float_t l2m, Float_t l2M) ; 
index 8add627..7efe03b 100644 (file)
@@ -40,19 +40,8 @@ ClassImp(AliPHOSRecParticle)
   fPx               = fE * momdir.X() ; 
   fPy               = fE * momdir.Y() ; 
   fPz               = fE * momdir.Z() ; 
-
-  fType = kUNDEFINED ; // undefined
-                            
-  if( ts->GetPpsdUp() == 0 ) {     // Neutral
-
-    if( ts->GetPpsdLow() == 0 )    // Neutral  
-      fType = kNEUTRAL ;   
-    else                           // Gamma
-      fType = kGAMMA ;               
-  }
-  else                            // Charged           
-    fType = kCHARGED ;   
-  
+  fType             = kUNDEFINED ;  
+                           
 }
 
 //____________________________________________________________________________
index eea1e4a..173879e 100644 (file)
@@ -75,5 +75,5 @@ AliPHOSReconstructioner::~AliPHOSReconstructioner()
 
   fTrackSegmentMaker->MakeTrackSegments(dl, emccl, ppsdl, trsl) ;
 
-  fPID->GetParticleType(trsl, rpl) ; 
+  fPID->MakeParticles(trsl, rpl) ; 
 }