Change pt cut by E cut for clusters

diff --git a/PWG4/PartCorrBase/AliCaloTrackReader.cxx b/PWG4/PartCorrBase/AliCaloTrackReader.cxx
index d841ef15d1b8e8e79e0923313a892fec0dd1c49c..139e641f3a5dd24ad4855182f884de1c2e85039c 100755 (executable)
--- a/PWG4/PartCorrBase/AliCaloTrackReader.cxx
+++ b/PWG4/PartCorrBase/AliCaloTrackReader.cxx
@@ -479,9 +479,11 @@ Bool_t AliCaloTrackReader::FillInputEvent(const Int_t iEntry, const char * /*cur
           bV0AND = fTriggerAnalysis->IsOfflineTriggerFired(esd, AliTriggerAnalysis::kV0AND);
         //else bV0AND = //FIXME FOR AODs
         if(!bV0AND) return kFALSE;
+        
       }
       
       if(!CheckForPrimaryVertex()) return kFALSE;
+      
     }//CaloFilter patch
     else{ 
       if(fInputEvent->GetNumberOfCaloClusters() > 0) {
@@ -809,7 +811,7 @@ void AliCaloTrackReader::FillInputEMCALAlgorithm(AliVCluster * clus, const Int_t
   
   clus->GetMomentum(momentum, fVertex[vindex]);      
   
-  if(fEMCALPtMin < momentum.Pt()){
+  if(fEMCALPtMin < momentum.E()){
     
     if(fCheckFidCut && !fFiducialCut->IsInFiducialCut(momentum,"EMCAL")) 
       return;
@@ -963,7 +965,7 @@ void AliCaloTrackReader::FillInputPHOS() {
         
         clus->GetMomentum(momentum, fVertex[vindex]);      
         
-        if(fPHOSPtMin < momentum.Pt()){
+        if(fPHOSPtMin < momentum.E()){
           
           if(fCheckFidCut && !fFiducialCut->IsInFiducialCut(momentum,"PHOS")) 
             continue;

diff --git a/PWG4/PartCorrDep/AliAnaPhoton.cxx b/PWG4/PartCorrDep/AliAnaPhoton.cxx
index 67a745ab68e14b65850c5f907af8aa7da509b604..39a7982c757384149cbcac1cd1e2cd102aab45b8 100755 (executable)
--- a/PWG4/PartCorrDep/AliAnaPhoton.cxx
+++ b/PWG4/PartCorrDep/AliAnaPhoton.cxx
@@ -158,8 +158,8 @@ TList *  AliAnaPhoton::GetCreateOutputObjects()
   fhNtraNclu->SetYTitle("# of clusters");
   outputContainer->Add(fhNtraNclu);
   
-  fhNCellsPt  = new TH2F ("hNCellsPt","# of tracks vs # of clusters", nptbins,ptmin, ptmax, 10,0,10); 
-  fhNCellsPt->SetXTitle("p_{T}");
+  fhNCellsPt  = new TH2F ("hNCellsPt","# of cells in cluster vs E of clusters", nptbins,ptmin, ptmax, 10,0,10); 
+  fhNCellsPt->SetXTitle("p_{T} (GeV/c)");
   fhNCellsPt->SetYTitle("# of cells in cluster");
   outputContainer->Add(fhNCellsPt);  
   
@@ -767,8 +767,8 @@ void  AliAnaPhoton::MakeAnalysisFillAOD()
  
     //.......................................
     //If too small or big pt, skip it
-    if(mom.Pt() < GetMinPt() || mom.Pt() > GetMaxPt() ) continue ; 
-    if(GetDebug() > 2) printf("\t Cluster %d Pass Pt Cut \n",icalo);
+    if(mom.E() < GetMinPt() || mom.E() > GetMaxPt() ) continue ; 
+    if(GetDebug() > 2) printf("\t Cluster %d Pass E Cut \n",icalo);
     
     //.......................................
     // TOF cut, BE CAREFUL WITH THIS CUT

diff --git a/PWG4/PartCorrDep/AliAnaPhoton.h b/PWG4/PartCorrDep/AliAnaPhoton.h
index 72d90a9a99db615df05535a0986041c2d7e5e855..ecece04ce742692f0cfebe7a9ec56f07b63ae236 100755 (executable)
--- a/PWG4/PartCorrDep/AliAnaPhoton.h
+++ b/PWG4/PartCorrDep/AliAnaPhoton.h
@@ -65,7 +65,6 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   TString GetCalorimeter()                  const {return fCalorimeter ; }
   void    SetCalorimeter(TString  & det)          {fCalorimeter = det  ; }
   
-  
   // ** Cluster selection methods **
   
   void    SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {

diff --git a/PWG4/PartCorrDep/AliAnaPi0.cxx b/PWG4/PartCorrDep/AliAnaPi0.cxx
index 6b7e00fac6f5be9ad466f59f34c8d5cae2b0240c..359b5cc1d397bdfe3d5e0e5cf279f8fbf80b8419 100755 (executable)
--- a/PWG4/PartCorrDep/AliAnaPi0.cxx
+++ b/PWG4/PartCorrDep/AliAnaPi0.cxx
@@ -1928,7 +1928,7 @@ void AliAnaPi0::MakeAnalysisFillHistograms()
             for(Int_t icell=0; icell<fNCellNCuts; icell++){
               for(Int_t iasym=0; iasym<fNAsymCuts; iasym++){
                 Int_t index = ((ipt*fNCellNCuts)+icell)*fNAsymCuts + iasym;
-                if(p1->Pt() >   fPtCuts[ipt]      && p2->Pt() > fPtCuts[ipt]        && 
+                if(p1->E() >   fPtCuts[ipt]      && p2->E() > fPtCuts[ipt]        && 
                    a        <   fAsymCuts[iasym]                                    && 
                    ncell1   >=  fCellNCuts[icell] && ncell2   >= fCellNCuts[icell]){
                     fhRePtNCellAsymCuts[index]->Fill(pt,m) ;