From f224d35b35fc6ef9e058626bdf366af5b8c60a5e Mon Sep 17 00:00:00 2001 From: loizides Date: Wed, 23 Mar 2011 10:17:02 +0000 Subject: [PATCH] Add some more cuts --- PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx | 16 +++++++++++++--- PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h | 8 ++++++-- 2 files changed, 19 insertions(+), 5 deletions(-) diff --git a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx index d4d282ccab5..33eae5ba03d 100644 --- a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx +++ b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx @@ -25,7 +25,6 @@ ClassImp(AliAnalysisTaskEMCALPi0PbPb) //________________________________________________________________________ AliAnalysisTaskEMCALPi0PbPb::AliAnalysisTaskEMCALPi0PbPb(const char *name) : AliAnalysisTaskSE(name), - fAsymMax(1), fCentVar("V0M"), fCentFrom(0), fCentTo(100), @@ -35,7 +34,10 @@ AliAnalysisTaskEMCALPi0PbPb::AliAnalysisTaskEMCALPi0PbPb(const char *name) fClusName(), fDoNtuple(0), fDoAfterburner(0), + fAsymMax(1), fNminCells(1), + fMinErat(0), + fMinEcc(0), fNEvs(0), fGeom(0), fOutput(0), @@ -393,7 +395,7 @@ void AliAnalysisTaskEMCALPi0PbPb::Terminate(Option_t *) fNtuple->Write(); } - AliInfo(Form("\nAccepted %lld events", fNEvs)); + AliInfo(Form("\n%s: Accepted %lld events", GetName(), fNEvs)); } //________________________________________________________________________ @@ -469,7 +471,7 @@ void AliAnalysisTaskEMCALPi0PbPb::FillClusHists() GetSigma(clus,maxAxis,minAxis); if (maxAxis > 0) clusterEcc = TMath::Sqrt(1.0 - minAxis*minAxis/(maxAxis*maxAxis)); - + clus->SetChi2(clusterEcc); // store ecc in chi2 fHClustEccentricity->Fill(clusterEcc); fHClustEtaPhi->Fill(clusterVec.Eta(),clusterVec.Phi()); fHClustEnergyPt->Fill(clusterVec.E(),clusterVec.Pt()); @@ -531,6 +533,10 @@ void AliAnalysisTaskEMCALPi0PbPb::FillPionHists() continue; if (clus1->GetNCells()E()Chi2()GetMomentum(clusterVec1,vertex); for (Int_t j = i+1; j(clusters->At(j)); @@ -542,6 +548,10 @@ void AliAnalysisTaskEMCALPi0PbPb::FillPionHists() continue; if (clus2->GetNCells()E()Chi2()GetMomentum(clusterVec2,vertex); pionVec = clusterVec1 + clusterVec2; Double_t pionZgg = TMath::Abs(clusterVec1.E()-clusterVec2.E())/pionVec.E(); diff --git a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h index e6f86c25814..69cf27277ee 100644 --- a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h +++ b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h @@ -34,6 +34,8 @@ class AliAnalysisTaskEMCALPi0PbPb : public AliAnalysisTaskSE { void SetClusName(const char *name) { fClusName = name; } void SetDoAfterburner(Bool_t b) { fDoAfterburner = b; } void SetFillNtuple(Bool_t b) { fDoNtuple = b; } + void SetMinEcc(Double_t ecc) { fMinEcc = ecc; } + void SetMinErat(Double_t erat) { fMinErat = erat; } void SetNminCells(Int_t n) { fNminCells = n; } void SetUseQualFlag(Bool_t b) { fUseQualFlag = b; } void SetVertexRange(Double_t z1, Double_t z2) { fVtxZMin=z1; fVtxZMax=z2; } @@ -48,7 +50,6 @@ class AliAnalysisTaskEMCALPi0PbPb : public AliAnalysisTaskSE { void GetSigma(AliVCluster *c, Double_t &sigmaMax, Double_t &sigmaMin); // input members - Double_t fAsymMax; // energy asymmetry max (def=1) TString fCentVar; // variable for centrality determination Double_t fCentFrom; // min centrality (def=0) Double_t fCentTo; // max centrality (def=100) @@ -58,7 +59,10 @@ class AliAnalysisTaskEMCALPi0PbPb : public AliAnalysisTaskSE { TString fClusName; // cluster branch name (def="") Bool_t fDoNtuple; // if true write out ntuple Bool_t fDoAfterburner; // if true run after burner - Int_t fNminCells; // minimum number of cells attached to cluster + Double_t fAsymMax; // maximum energy asymmetry (def=1) + Int_t fNminCells; // minimum number of cells attached to cluster (def=1) + Double_t fMinErat; // minimum emax/ec ratio (def=0) + Double_t fMinEcc; // minimum eccentricity (def=0) // derived members (ie with ! after //) ULong64_t fNEvs; //!accepted events AliEMCALGeoUtils *fGeom; //!geometry utils -- 2.43.5