AnaExample: Bug with handling calo cells corrected
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Jun 2009 09:32:24 +0000 (09:32 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Jun 2009 09:32:24 +0000 (09:32 +0000)
CaloPID: In EMCAL trach matching not used in PID until now, photon and electron weights summed to get the correct photon weight.

PWG4/PartCorrBase/AliCaloPID.cxx
PWG4/PartCorrDep/AliAnaExample.cxx

index 19b3648..d124a37 100755 (executable)
@@ -153,11 +153,11 @@ Int_t AliCaloPID::GetPdg(const TString calo, const Double_t * pid, const Float_t
     abort();
   }
   
-  Float_t wPh =  fPHOSPhotonWeight ;
+  Float_t wPh  =  fPHOSPhotonWeight ;
   Float_t wPi0 =  fPHOSPi0Weight ;
-  Float_t wE =  fPHOSElectronWeight ;
-  Float_t wCh =  fPHOSChargeWeight ;
-  Float_t wNe =  fPHOSNeutralWeight ;
+  Float_t wE   =  fPHOSElectronWeight ;
+  Float_t wCh  =  fPHOSChargeWeight ;
+  Float_t wNe  =  fPHOSNeutralWeight ;
   
   
   if(calo == "PHOS" && fPHOSWeightFormula){
@@ -202,9 +202,10 @@ Int_t AliCaloPID::GetPdg(const TString calo, const Double_t * pid, const Float_t
       pdg = kNeutralUnknown ;
   }
   else{//EMCAL
-    if(pid[AliAODCluster::kPhoton]  > wPh) pdg = kPhoton ;
+    if(pid[AliAODCluster::kPhoton]+pid[AliAODCluster::kElectron]  > wPh) pdg = kPhoton ; //temporal sollution until track matching for electrons is considered
+    //if(pid[AliAODCluster::kPhoton]  > wPh) pdg = kPhoton ;
     else if(pid[AliAODCluster::kPi0] > wPi0) pdg = kPi0 ; 
-    else if(pid[AliAODCluster::kElectron]  > wE) pdg = kElectron ;
+    //else if(pid[AliAODCluster::kElectron]  > wE) pdg = kElectron ;
     else if(chargedHadronWeight + neutralHadronWeight > wCh) pdg = kChargedHadron ;  
     else if(neutralHadronWeight + chargedHadronWeight > wNe) pdg = kNeutralHadron ; 
     else pdg =  kNeutralUnknown ;
index b509552..cc092a2 100755 (executable)
@@ -262,11 +262,11 @@ void  AliAnaExample::MakeAnalysisFillAOD()
       
       AliESDCaloCells * esdCell = new AliESDCaloCells ;
       if(fDetector == "PHOS") {
-       //ConnectAODPHOSCells(); //Do Only when filling AODCaloCells
+       ConnectAODPHOSCells(); //Do Only when filling AODCaloCells
        esdCell = (AliESDCaloCells *) GetPHOSCells();
       }
       else  {
-       //ConnectAODEMCALCells(); //Do Only when filling AODCaloCells
+       ConnectAODEMCALCells(); //Do Only when filling AODCaloCells
        esdCell = (AliESDCaloCells *) GetEMCALCells();
       }