AliAnalysisTaskEmcalJetCDF.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALJETCDF_H
#define ALIANALYSISTASKEMCALJETCDF_H

#include <cstdio>
#include <iostream>
#include <vector>
#include <utility>
#include <algorithm>

#include <Rtypes.h>
#include <TH1.h>
#include <TH1D.h>
#include <TProfile.h>
#include <TMath.h>

#include <AliEmcalJet.h>
#include <AliJetContainer.h>
#include <AliParticleContainer.h>
#include <AliClusterContainer.h>
#include <THistManager.h>

#include "AliAnalysisTaskEmcalJet.h"

using std::cout;
using std::endl;
using std::vector;
using std::pair;

class TH1;
class TProfile;
class AliAnalysisUtils;
class AliJetContainer;
class AliParticleContainer;
class AliClusterContainer;

class AliAnalysisTaskEmcalJetCDF : public AliAnalysisTaskEmcalJet
  {
  public:

    AliAnalysisTaskEmcalJetCDF();
    AliAnalysisTaskEmcalJetCDF ( const char *name );
    virtual ~AliAnalysisTaskEmcalJetCDF();

    void                        UserCreateOutputObjects();
    void                        Terminate ( Option_t *option );

    void                        SetProcessJetCont (Int_t i) {idx_jetcont = i;}

    Int_t                       GetProcessJetCont () {return idx_jetcont;}

    Double_t                    Z_ptot ( const AliEmcalJet* jet, const AliVParticle* trk )  const; // Get Z of constituent trk ; p total

    Double_t                    Z_pt   ( const AliEmcalJet* jet, const AliVParticle* trk )  const; // Get Z of constituent trk ; pt

    Double_t                    Xi ( Double_t z )  const { return TMath::Log ( 1/z ); } // Get Xi of value z

  protected:
    void                        ExecOnce();
    Bool_t                      Run() ;

    Bool_t                      FillHistograms()   ;

    Bool_t                      ProcessJetContainer(Int_t idx_jet_container = 0);

    std::vector<Int_t>          SortTracksPt ( AliVEvent *event ) const;

    std::vector<Int_t>          SortTracksPt ( AliParticleContainer *track_container ) const;

    Double_t                    DeltaR ( const AliVParticle *part1, const AliVParticle *part2 );

    Bool_t                      IdxInArray ( Int_t index, TArrayI &array );

    // Histogram list
    TH1D      *fH1;                
    TH1D      *fH2;                
    TH1D      *fH3;                

//________________________________________________________________________
    TH1D      *fH4;                
    TH1D      *fH4c;               

//________________________________________________________________________
    TH1D      *fH5;                

//________________________________________________________________________
    TH1D      *fH6;                
    TH1D      *fH6c;               

//________________________________________________________________________
    TProfile  *fH7;                
    TProfile  *fH7all;             

//________________________________________________________________________
    TH1D      *fH8;                
    TH1D      *fH8xi;              

//________________________________________________________________________
    TH1D      *fH8_p;                
    TH1D      *fH8xi_p;              

//________________________________________________________________________
    TH1D      *fH8_pt;                
    TH1D      *fH8xi_pt;              

//________________________________________________________________________
    TH1D      *fH8_all;            
    TH1D      *fH8xi_all;          

//________________________________________________________________________
    TH1D      *fH8_all_p;            
    TH1D      *fH8xi_all_p;          

//________________________________________________________________________
    TH1D      *fH8_all_pt;            
    TH1D      *fH8xi_all_pt;          

//________________________________________________________________________
    TProfile  *fH15;               

    TProfile  *fH15_n90;           
    TProfile  *fH15_n85;           
    TProfile  *fH15_n80;           
    TProfile  *fH15_n75;           
    TProfile  *fH15_n70;           

    TProfile  *fH15_pt90;          
    TProfile  *fH15_pt85;          
    TProfile  *fH15_pt80;          
    TProfile  *fH15_pt75;          
    TProfile  *fH15_pt70;          

//________________________________________________________________________
    TProfile  *fH15all;            

    TProfile  *fH15all_n90;        
    TProfile  *fH15all_n85;        
    TProfile  *fH15all_n80;        
    TProfile  *fH15all_n75;        
    TProfile  *fH15all_n70;        

    TProfile  *fH15all_pt90;       
    TProfile  *fH15all_pt85;       
    TProfile  *fH15all_pt80;       
    TProfile  *fH15all_pt75;       
    TProfile  *fH15all_pt70;       

//________________________________________________________________________
    TH1D      *fH15_bin;           

    TH1D      *fH15_bin_n90;       
    TH1D      *fH15_bin_n85;       
    TH1D      *fH15_bin_n80;       
    TH1D      *fH15_bin_n75;       
    TH1D      *fH15_bin_n70;       

    TH1D      *fH15_bin_pt90;      
    TH1D      *fH15_bin_pt85;      
    TH1D      *fH15_bin_pt80;      
    TH1D      *fH15_bin_pt75;      
    TH1D      *fH15_bin_pt70;      

//________________________________________________________________________
    TH1D      *fH15all_bin;       

    TH1D      *fH15all_bin_n90;   
    TH1D      *fH15all_bin_n85;   
    TH1D      *fH15all_bin_n80;   
    TH1D      *fH15all_bin_n75;   
    TH1D      *fH15all_bin_n70;   

    TH1D      *fH15all_bin_pt90;  
    TH1D      *fH15all_bin_pt85;  
    TH1D      *fH15all_bin_pt80;  
    TH1D      *fH15all_bin_pt75;  
    TH1D      *fH15all_bin_pt70;  

//________________________________________________________________________
    TH1D      *fH20;               

    TH1D      *fH20_n90;           
    TH1D      *fH20_n85;           
    TH1D      *fH20_n80;           
    TH1D      *fH20_n75;           
    TH1D      *fH20_n70;           

    TH1D      *fH20_pt90;          
    TH1D      *fH20_pt85;          
    TH1D      *fH20_pt80;          
    TH1D      *fH20_pt75;          
    TH1D      *fH20_pt70;          


//________________________________________________________________________
    TH1D      *fH20all;           

    TH1D      *fH20all_n90;       
    TH1D      *fH20all_n85;       
    TH1D      *fH20all_n80;       
    TH1D      *fH20all_n75;       
    TH1D      *fH20all_n70;       

    TH1D      *fH20all_pt90;      
    TH1D      *fH20all_pt85;      
    TH1D      *fH20all_pt80;      
    TH1D      *fH20all_pt75;      
    TH1D      *fH20all_pt70;      

//________________________________________________________________________
    TH1D      *fHg;                

    TH1D      *fHg_n90;            
    TH1D      *fHg_n85;            
    TH1D      *fHg_n80;            
    TH1D      *fHg_n75;            
    TH1D      *fHg_n70;            

    TH1D      *fHg_pt90;           
    TH1D      *fHg_pt85;           
    TH1D      *fHg_pt80;           
    TH1D      *fHg_pt75;           
    TH1D      *fHg_pt70;           

//________________________________________________________________________
    TH1D      *fHptd;                

    TH1D      *fHptd_n90;            
    TH1D      *fHptd_n85;            
    TH1D      *fHptd_n80;            
    TH1D      *fHptd_n75;            
    TH1D      *fHptd_n70;            

    TH1D      *fHptd_pt90;           
    TH1D      *fHptd_pt85;           
    TH1D      *fHptd_pt80;           
    TH1D      *fHptd_pt75;           
    TH1D      *fHptd_pt70;           

//________________________________________________________________________
    TProfile  *fH23;               
    TProfile  *fH24;               
    TProfile  *fH25;               

//________________________________________________________________________
    TProfile  *fH23all;            
    TProfile  *fH24all;            
    TProfile  *fH25all;            


    //_________________________________________________________________________
    AliJetContainer           *fJetsCont;                   
    AliParticleContainer      *fTracksCont;                 
    AliClusterContainer       *fCaloClustersCont;           
    TClonesArray              *fTracksContArray;            
    TClonesArray              *fCaloClustContArray;         

    Int_t                      idx_jetcont;                 
    Int_t                      fNJets_accepted;             
    Int_t                      fNaccPart;                   
    Int_t                      fNaccClus;                   
    THistManager               fHistManager;                

  private:
    typedef std::pair<Double_t, Int_t> ptidx_pair;

    struct sort_descend
      {
      bool operator () ( const ptidx_pair &p1, const ptidx_pair &p2 )  { return p1.first > p2.first ; }
      };

    AliAnalysisTaskEmcalJetCDF ( const AliAnalysisTaskEmcalJetCDF& );           // not implemented
    AliAnalysisTaskEmcalJetCDF &operator= ( const AliAnalysisTaskEmcalJetCDF& ); // not implemented

    ClassDef ( AliAnalysisTaskEmcalJetCDF, 4 );

  };
#endif

// kate: indent-mode none; indent-width 2; replace-tabs on;