Change pt cut by E cut for clusters
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 10 Apr 2011 14:18:42 +0000 (14:18 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 10 Apr 2011 14:18:42 +0000 (14:18 +0000)
PWG4/PartCorrBase/AliCaloTrackReader.cxx
PWG4/PartCorrDep/AliAnaPhoton.cxx
PWG4/PartCorrDep/AliAnaPhoton.h
PWG4/PartCorrDep/AliAnaPi0.cxx

index d841ef1..139e641 100755 (executable)
@@ -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;
index 67a745a..39a7982 100755 (executable)
@@ -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
index 72d90a9..ecece04 100755 (executable)
@@ -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) {
index 6b7e00f..359b5cc 100755 (executable)
@@ -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) ;