From ba1eeb1f8f18d03c2a62a4a38d2e428b3fcb7be7 Mon Sep 17 00:00:00 2001 From: gconesab Date: Sun, 10 Apr 2011 14:18:42 +0000 Subject: [PATCH] Change pt cut by E cut for clusters --- PWG4/PartCorrBase/AliCaloTrackReader.cxx | 6 ++++-- PWG4/PartCorrDep/AliAnaPhoton.cxx | 8 ++++---- PWG4/PartCorrDep/AliAnaPhoton.h | 1 - PWG4/PartCorrDep/AliAnaPi0.cxx | 2 +- 4 files changed, 9 insertions(+), 8 deletions(-) diff --git a/PWG4/PartCorrBase/AliCaloTrackReader.cxx b/PWG4/PartCorrBase/AliCaloTrackReader.cxx index d841ef15d1b..139e641f3a5 100755 --- 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 67a745ab68e..39a7982c757 100755 --- 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 72d90a9a99d..ecece04ce74 100755 --- 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 6b7e00fac6f..359b5cc1d39 100755 --- a/PWG4/PartCorrDep/AliAnaPi0.cxx +++ b/PWG4/PartCorrDep/AliAnaPi0.cxx @@ -1928,7 +1928,7 @@ void AliAnaPi0::MakeAnalysisFillHistograms() for(Int_t icell=0; icellPt() > 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) ; -- 2.43.0