From 50878ee186f78534d2e9e626b425969166f9b6f2 Mon Sep 17 00:00:00 2001 From: lmilano Date: Tue, 3 Dec 2013 11:35:35 +0000 Subject: [PATCH] rapidity cut according to the mass and not using AliSpectraNameSpace --- .../SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskSpectraAOD.cxx | 8 +++++--- PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.cxx | 6 ++---- PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.h | 5 +---- 3 files changed, 8 insertions(+), 11 deletions(-) diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskSpectraAOD.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskSpectraAOD.cxx index f5af53758e8..43b548e1451 100644 --- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskSpectraAOD.cxx +++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskSpectraAOD.cxx @@ -135,6 +135,7 @@ void AliAnalysisTaskSpectraAOD::UserExec(Option_t *) } } + Double_t mass[3]={1.39570000000000000e-01,4.93676999999999977e-01,9.38271999999999995e-01};//FIXME masses to be taken from AliHelperPID //main loop on tracks for (Int_t iTracks = 0; iTracks < fAOD->GetNumberOfTracks(); iTracks++) { AliAODTrack* track = fAOD->GetTrack(iTracks); @@ -160,8 +161,9 @@ void AliAnalysisTaskSpectraAOD::UserExec(Option_t *) Int_t charge = track->Charge() > 0 ? kChPos : kChNeg; // Fill histograms, only if inside y and nsigma acceptance - if(idRec != kSpUndefined && fTrackCuts->CheckYCut ((AODParticleSpecies_t)idRec))fHistMan->GetHistogram2D(kHistPtRecSigma,idRec,charge)->Fill(track->Pt(),dca); - //can't put a continue because we still have to fill allcharged primaries, done later + if(idRec != kSpUndefined){ + if(fTrackCuts->CheckYCut (mass[idRec]))fHistMan->GetHistogram2D(kHistPtRecSigma,idRec,charge)->Fill(track->Pt(),dca); + }//can't put a continue because we still have to fill allcharged primaries, done later /* MC Part */ if (arrayMC) { @@ -196,7 +198,7 @@ void AliAnalysisTaskSpectraAOD::UserExec(Option_t *) if(idRec == kSpUndefined) continue; // rapidity cut (reconstructed pt and identity) - if(!fTrackCuts->CheckYCut ((AODParticleSpecies_t)idRec)) continue; + if(!fTrackCuts->CheckYCut (mass[idRec])) continue; // Get true ID Int_t idGen = fPID->GetParticleSpecie(partMC); diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.cxx index b99328405ee..aef00914795 100644 --- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.cxx +++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.cxx @@ -157,13 +157,11 @@ Bool_t AliSpectraAODTrackCuts::CheckEtaCut() return kFALSE; } -Bool_t AliSpectraAODTrackCuts::CheckYCut(AODParticleSpecies_t species) +Bool_t AliSpectraAODTrackCuts::CheckYCut(Double_t mass) { // check if the rapidity is within the set range Double_t y=-1000; - if (species == kSpProton) { y = fTrack->Y(9.38271999999999995e-01); } - if ( species == kSpKaon ) { y = fTrack->Y(4.93676999999999977e-01); } - if ( species == kSpPion) { y = fTrack->Y(1.39570000000000000e-01); } + if (mass > 0.) { y = fTrack->Y(mass); }//negative mass for unidentified particles if (TMath::Abs(y) > fYCut || y < -998.) return kFALSE; return kTRUE; } diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.h b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.h index bdfa8954eda..3c2bf175361 100644 --- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.h +++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.h @@ -18,10 +18,7 @@ class AliPIDResponse; class AliAODMCParticle; class AliAODTrack; -#include "AliSpectraAODHistoManager.h" #include "TNamed.h" - -using namespace AliSpectraNameSpace; class AliSpectraAODTrackCuts : public TNamed { @@ -65,7 +62,7 @@ class AliSpectraAODTrackCuts : public TNamed Bool_t CheckTrackType(); Bool_t CheckTrackCuts(); Bool_t CheckEtaCut(); - Bool_t CheckYCut(AODParticleSpecies_t specie); // not included in standard cuts + Bool_t CheckYCut(Double_t mass); // not included in standard cuts Bool_t CheckDCACut(); Bool_t CheckPCut(); Bool_t CheckPtCut(); -- 2.43.0