setter for upper and lower limits

diff --git a/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx b/PWGLF/SPECTRA/Kinks/AliAnalysisKinkESDat.cxx
index 0ae0bf887257f4d1d11e9f62207146ff239f888a..b5a226d76d2ae8a48068ec3558ab4375a756d4fc 100644 (file)
--- 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() <fKinkRadUp   )  )  {
     //  for systematics   if( ( kink->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 4b9f90ced7b3b8a388e04af5118ff5c97464591f..addcfc683a38fce0fe8b86124cea482cabaecd65 100644 (file)
--- 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