From 783b974ca6273ce6efb4eeab398207aa748e1559 Mon Sep 17 00:00:00 2001 From: gconesab Date: Sun, 10 Aug 2014 11:25:42 +0200 Subject: [PATCH] remove hardcoded minimum energy for charged and neutrals used for isolation --- .../AliAnaGeneratorKine.cxx | 24 ++++++++++--------- .../AliAnaGeneratorKine.h | 5 ++++ 2 files changed, 18 insertions(+), 11 deletions(-) diff --git a/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.cxx b/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.cxx index 2355034cd3e..d1d283af89a 100755 --- a/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.cxx +++ b/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.cxx @@ -37,8 +37,8 @@ ClassImp(AliAnaGeneratorKine) //__________________________________________ AliAnaGeneratorKine::AliAnaGeneratorKine() : AliAnaCaloTrackCorrBaseClass(), -fTriggerDetector("EMCAL"), -fCalorimeter("EMCAL"), +fTriggerDetector(""),fCalorimeter(""), +fMinChargedPt(0), fMinNeutralPt(0), fStack(0), fParton2(0), fParton3(0), fParton6(0), fParton7(0), @@ -646,8 +646,6 @@ void AliAnaGeneratorKine::GetXE(TLorentzVector trigger, // Calculate the real XE and the UE XE - Float_t ptThresTrack = 0.2; - Float_t ptTrig = trigger.Pt(); Float_t etaTrig = trigger.Eta(); Float_t phiTrig = trigger.Phi(); @@ -676,7 +674,7 @@ void AliAnaGeneratorKine::GetXE(TLorentzVector trigger, Float_t phi = particle->Phi(); if(phi < 0 ) phi += TMath::TwoPi(); - if( pt < ptThresTrack) continue ; + if( pt < fMinChargedPt) continue ; if(TMath::Abs(eta) > 0.8) continue ; // TPC acceptance cut @@ -772,6 +770,12 @@ void AliAnaGeneratorKine::InitParameters() //Initialize the parameters of the analysis. AddToHistogramsName("AnaGenKine_"); + fCalorimeter = "EMCAL"; + fTriggerDetector = "EMCAL"; + + fMinChargedPt = 0.2; + fMinNeutralPt = 0.3; + } //_____________________________________________________________________ @@ -806,12 +810,10 @@ void AliAnaGeneratorKine::IsLeadingAndIsolated(TLorentzVector trigger, if(phiTrig < 0 ) phiTrig += TMath::TwoPi(); // Minimum track or cluster energy - Float_t ptThresTrack = 0.2; - Float_t ptThresCalo = 0.3; //Isolation cuts - Float_t ptThresIC = 0.5; - Float_t rThresIC = 0.4; + Float_t ptThresIC = GetIsolationCut()->GetPtTreshold(); + Float_t rThresIC = GetIsolationCut()->GetConeSize(); Int_t nICTrack = 0; Int_t nICNeutral = 0; Int_t nICNeutEMCAL = 0; @@ -849,7 +851,7 @@ void AliAnaGeneratorKine::IsLeadingAndIsolated(TLorentzVector trigger, if(charge==0) { - if(pt < ptThresCalo) continue ; + if(pt < fMinNeutralPt) continue ; if( ptMaxNeutral < pt ) ptMaxNeutral = pt; if( pt > ptThresIC && radius < rThresIC ) nICNeutral++ ; @@ -881,7 +883,7 @@ void AliAnaGeneratorKine::IsLeadingAndIsolated(TLorentzVector trigger, } else { - if( pt < ptThresTrack) continue ; + if( pt < fMinChargedPt) continue ; if( ptMaxCharged < pt ) ptMaxCharged = pt; diff --git a/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.h b/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.h index 986a497387b..0c59e3be2b9 100755 --- a/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.h +++ b/PWGGA/CaloTrackCorrelations/AliAnaGeneratorKine.h @@ -57,6 +57,8 @@ public: void SetTriggerDetector( TString name ) { fTriggerDetector = name ; } void SetCalorimeter ( TString name ) { fCalorimeter = name ; } + void SetMinChargedPt ( Float_t pt ) { fMinChargedPt = pt ; } + void SetMinNeutralPt ( Float_t pt ) { fMinNeutralPt = pt ; } private: @@ -64,6 +66,9 @@ private: TString fTriggerDetector; //! trigger detector, for fiducial region TString fCalorimeter; //! detector neutral particles, for fiducial region + Float_t fMinChargedPt; //! Minimum energy for charged particles in correlation + Float_t fMinNeutralPt; //! Minimum energy for neutral particles in correlation + AliStack * fStack; //! access stack TParticle * fParton2; //! Initial state Parton -- 2.43.0