//__________________________________________
AliAnaGeneratorKine::AliAnaGeneratorKine() :
AliAnaCaloTrackCorrBaseClass(),
-fTriggerDetector("EMCAL"),
-fCalorimeter("EMCAL"),
+fTriggerDetector(""),fCalorimeter(""),
+fMinChargedPt(0), fMinNeutralPt(0),
fStack(0),
fParton2(0), fParton3(0),
fParton6(0), fParton7(0),
// 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();
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
//Initialize the parameters of the analysis.
AddToHistogramsName("AnaGenKine_");
+ fCalorimeter = "EMCAL";
+ fTriggerDetector = "EMCAL";
+
+ fMinChargedPt = 0.2;
+ fMinNeutralPt = 0.3;
+
}
//_____________________________________________________________________
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;
if(charge==0)
{
- if(pt < ptThresCalo) continue ;
+ if(pt < fMinNeutralPt) continue ;
if( ptMaxNeutral < pt ) ptMaxNeutral = pt;
if( pt > ptThresIC && radius < rThresIC ) nICNeutral++ ;
}
else
{
- if( pt < ptThresTrack) continue ;
+ if( pt < fMinChargedPt) continue ;
if( ptMaxCharged < pt ) ptMaxCharged = pt;
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:
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