From b2fb60b825b0828972620081bbedf1bc28643dea Mon Sep 17 00:00:00 2001 From: spyropo Date: Thu, 14 Feb 2013 14:32:39 +0000 Subject: [PATCH] setter for upper and lower limits --- PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx | 8 +++++--- PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.h | 3 +++ 2 files changed, 8 insertions(+), 3 deletions(-) diff --git a/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx b/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx index 0ae0bf88725..b5a226d76d2 100644 --- a/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx +++ b/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx @@ -61,7 +61,7 @@ AliAnalysisKinkESDat::AliAnalysisKinkESDat(const char *name) fRatioCrossedRows(0), fRatioCrossedRowsKink(0),fRadiusPt(0), fRadiusPtcln(0), fInvMassMuNuPt(0), fPtCut1(0), fPtCut2(0), fPtCut3(0), fAngMomKKinks(0), f1(0), f2(0), - fListOfHistos(0),fLowMulcut(-1),fUpMulcut(-1),fCutsMul(0), fMaxDCAtoVtxCut(0), fPIDResponse(0) + fListOfHistos(0),fLowMulcut(-1),fUpMulcut(-1), fKinkRadUp(200), fKinkRadLow(130), fCutsMul(0), fMaxDCAtoVtxCut(0), fPIDResponse(0) { // Constructor @@ -582,7 +582,8 @@ void AliAnalysisKinkESDat::UserExec(Option_t *) // fRadiusPt->Fill( kink->GetR(), trackPt); // // radius and Minv selection - if( ( kink->GetR()> 120 ) && ( kink->GetR() < 210 ) ) { + //if( ( kink->GetR()> 120 ) && ( kink->GetR() < 210 ) ) { + if( ( kink->GetR()> fKinkRadLow ) && ( kink->GetR() GetR()> 130 ) && ( kink->GetR() < 200 ) ) { if (qT>0.12) fAngMomKC->Fill(track->P(), kinkAngle); if ( qT>0.12) fM1kaon->Fill(invariantMassKmu); @@ -602,7 +603,8 @@ void AliAnalysisKinkESDat::UserExec(Option_t *) // fHistPtKPDG->Fill(track->Pt()); // ALL K-candidates until now // if((kinkAngle>maxDecAngpimu)&&(qT>0.12)&&(qT<0.30)&&((kink->GetR()>=120.)&&(kink->GetR()<=210.))&&(TMath::Abs(rapiditK)<0.7)&&(invariantMassKmu<0.6)){ - if((kinkAngle>maxDecAngpimu)&&(qT>0.12)&&(qT<0.30)&&((kink->GetR()>=120.)&&(kink->GetR()<=210.))&&(TMath::Abs(rapiditK)<0.7)&&(invariantMassKmu<0.8)){ + //if((kinkAngle>maxDecAngpimu)&&(qT>0.12)&&(qT<0.30)&&((kink->GetR()>=120.)&&(kink->GetR()<=210.))&&(TMath::Abs(rapiditK)<0.7)&&(invariantMassKmu<0.8)){ + if((kinkAngle>maxDecAngpimu)&&(qT>0.12)&&(qT<0.30)&&((kink->GetR()>= fKinkRadLow )&&(kink->GetR()<= fKinkRadUp ))&&(TMath::Abs(rapiditK)<0.7)&&(invariantMassKmu<0.8)){ // systematics if((kinkAngle>maxDecAngpimu)&&(qT>0.12)&&(qT<0.30)&&((kink->GetR()>=130.)&&(kink->GetR()<=200.))&&(TMath::Abs(rapiditK)<0.7)&&(invariantMassKmu<0.8)){ // fAngMomKKinks->Fill(track->P(), kinkAngle); diff --git a/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.h b/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.h index 4b9f90ced7b..addcfc683a3 100644 --- a/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.h +++ b/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.h @@ -38,6 +38,7 @@ class AliAnalysisKinkESDat : public AliAnalysisTaskSE { const AliESDVertex *GetEventVertex(const AliESDEvent* esd) const; // Marek multiplicity bins void SetMulCut(Int_t low, Int_t up){fLowMulcut=low;fUpMulcut=up;} + void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp) { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;} private: TH1F *fHistPtESD; //Pt spectrum of all ESD inside eta, Pt cuts @@ -123,6 +124,8 @@ class AliAnalysisKinkESDat : public AliAnalysisTaskSE { // Marek Multiplicity Int_t fLowMulcut; // Int_t fUpMulcut; +Int_t fKinkRadUp; +Int_t fKinkRadLow; AliESDtrackCuts* fCutsMul; AliESDtrackCuts* fMaxDCAtoVtxCut; AliPIDResponse *fPIDResponse; //! PID response object -- 2.43.0