AliAnaGeneratorKine.h

Go to the documentation of this file.

#ifndef ALIANAGENERATORKINE_H
#define ALIANAGENERATORKINE_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */

//___________________________________________________________________________
//_________________________________________________________________________

// --- ROOT ---
class TH2F ;
class TParticle ;
class AliStack ;
class TLorentzVector ;

// --- ANALYSIS ---
#include "AliAnaCaloTrackCorrBaseClass.h"

class AliAnaGeneratorKine : public AliAnaCaloTrackCorrBaseClass {
       
public:
  
  AliAnaGeneratorKine() ;

  virtual ~AliAnaGeneratorKine() { delete fFidCutTrigger ; }
  
  enum mcPrimTypes { kmcPrimPhoton = 0, kmcPrimPi0Decay = 1, kmcPrimEtaDecay  = 2, kmcPrimOtherDecay  = 3,
                     kmcPrimPi0    = 4, kmcPrimEta      = 5                                               } ;
  
  static const Int_t fgkNmcPrimTypes = 6;

  static const Int_t fgkNLead = 2;

  static const Int_t fgkNIso  = 4;

  Bool_t CorrelateWithPartonOrJet(Int_t   indexTrig,
                                  Int_t   pdgTrig,
                                  Bool_t  leading [fgkNIso],
                                  Bool_t  isolated[fgkNIso],
                                  Int_t & iparton) ; 
  
  TList * GetCreateOutputObjects() ;
  
  void    GetPartonsAndJets() ;
    
  void    GetXE(Int_t   indexTrig,
                Int_t   pdgTrig,
                Bool_t  leading [fgkNIso],
                Bool_t  isolated[fgkNIso],
                Int_t   iparton) ;
  
  void    InitParameters() ;
  
  void    IsLeadingAndIsolated(Int_t  indexTrig,
                               Int_t  pdgTrig,
                               Bool_t leading [fgkNIso],     
                               Bool_t isolated[fgkNIso]) ;
    
  void    MakeAnalysisFillHistograms() ;
  
  void    SetTriggerDetector( TString & det ) ;
  void    SetTriggerDetector( Int_t  det )    ;
  
  void    SetMinChargedPt   ( Float_t pt )   { fMinChargedPt    = pt   ; }
  void    SetMinNeutralPt   ( Float_t pt )   { fMinNeutralPt    = pt   ; }
  
  // Detector for trigger particles acceptance
  AliFiducialCut * GetFiducialCutForTrigger()
  { if(!fFidCutTrigger)  fFidCutTrigger  = new AliFiducialCut(); return  fFidCutTrigger  ; }
  virtual void     SetFiducialCut(AliFiducialCut * fc)
  { delete fFidCutTrigger;  fFidCutTrigger  = fc      ; }
  
private:

  Int_t            fTriggerDetector ;       
  TString          fTriggerDetectorString ; 

  AliFiducialCut * fFidCutTrigger;          
  
  Float_t          fMinChargedPt;           
  Float_t          fMinNeutralPt;           
  
  AliStack       * fStack;                  
  TClonesArray   * fAODMCparticles ;        

//  TParticle      * fParton2;              //!<! Initial state Parton
//  TParticle      * fParton3;              //!<! Initial state Parton
  
  TLorentzVector   fParton6;                
  TLorentzVector   fParton7;                
  
  Int_t            fParton6PDG;             
  Int_t            fParton7PDG;             
  
  TLorentzVector   fJet6;                   
  TLorentzVector   fJet7;                   

  TLorentzVector   fTrigger;                
  TLorentzVector   fLVTmp;                  
  
  Int_t            fNPrimaries;             
  Float_t          fPtHard;                 
  
  // Histograms
  
  TH1F      * fhPtHard;                     
  TH1F      * fhPtParton;                   
  TH1F      * fhPtJet;                      
  
  TH2F      * fhPtPartonPtHard;             
  TH2F      * fhPtJetPtHard;                
  TH2F      * fhPtJetPtParton;              

  TH1F      * fhPt      [fgkNmcPrimTypes];  
  TH2F      * fhPtOrigin[fgkNmcPrimTypes];  
  TH2F      * fhPtOriginNotFinal[fgkNmcPrimTypes];  
  TH2F      * fhPtOtherDecayMesonId;        
  
  TH1F      * fhPhi      [fgkNmcPrimTypes]; 
  TH1F      * fhEta      [fgkNmcPrimTypes]; 
  TH2F      * fhEtaPhi   [fgkNmcPrimTypes]; 

  TH2F      * fhPhiStatus[fgkNmcPrimTypes]; 
  TH2F      * fhEtaStatus[fgkNmcPrimTypes]; 
  
  TH2F      * fhPtPi0Status;                
  TH2F      * fhPtEtaStatus;                
  TH2F      * fhPtPi0DecayStatus;           
  TH2F      * fhPtEtaDecayStatus;           
  TH2F      * fhPtOtherDecayStatus;         

  TH1F      * fhPtPi0Not2Gamma;             
  TH1F      * fhPtEtaNot2Gamma;             
  TH1F      * fhPtGammaFromPi0Not2Gamma;    
  TH1F      * fhPtGammaFromEtaNot2Gamma;    
  
  // Histograms arrays for 4 isolation options and 2 options on leading or not leading particle

  TH2F      * fhPtAcceptedGammaJet                       [fgkNLead][fgkNIso]; 

  
  TH1F      * fhPtLeading               [fgkNmcPrimTypes]          [fgkNIso]; 

  TH2F      * fhPtLeadingSumPt          [fgkNmcPrimTypes]          [fgkNIso]; 
  
  TH1F      * fhPtLeadingIsolated       [fgkNmcPrimTypes]          [fgkNIso]; 

  TH2F      * fhPtPartonTypeNear        [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
  
  TH2F      * fhPtPartonTypeNearIsolated[fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhPtPartonTypeAway        [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhPtPartonTypeAwayIsolated[fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
  
  TH2F      * fhZHard                   [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
 
  TH2F      * fhZHardIsolated           [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
  
  TH2F      * fhZParton                 [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhZPartonIsolated         [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhZJet                    [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhZJetIsolated            [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
  
  TH2F      * fhXE                      [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhXEIsolated              [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
  
  TH2F      * fhXEUE                    [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 

  TH2F      * fhXEUEIsolated            [fgkNmcPrimTypes][fgkNLead][fgkNIso]; 
  
  AliAnaGeneratorKine              (const AliAnaGeneratorKine & gk) ;

  AliAnaGeneratorKine & operator = (const AliAnaGeneratorKine & gk) ; 
  
  ClassDef(AliAnaGeneratorKine,7) ;

} ;


#endif //ALIANAGENERATORKINE_H