AliAnalysisTaskEmcalJetCDFUE.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALJETCDFUE_H
#define ALIANALYSISTASKEMCALJETCDFUE_H

#include <cstdio>
#include <Rtypes.h>
#include <TROOT.h>
#include <TSystem.h>
#include <TProfile.h>
#include <TH1.h>
#include <TH1F.h>
#include <TH1D.h>
#include <TProfile.h>
#include <TMath.h>

#include "AliLog.h"
#include "AliEmcalJet.h"
#include "AliJetContainer.h"
#include "AliParticleContainer.h"
#include "AliClusterContainer.h"

#include "AliAnalysisTaskEmcalJet.h"

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

class AliAnalysisTaskEmcalJetCDFUE : public AliAnalysisTaskEmcalJet
  {
  public:

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

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

  protected:

    Bool_t                      FillHistograms()   ;

    void                        ExecOnce();
    Bool_t                      Run() ;

    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 );

    TProfile  *fH21;               
    TProfile  *fH21Toward;         
    TProfile  *fH21Transverse_min; 
    TProfile  *fH21Transverse_max; 
    TProfile  *fH21Away;           

    TH1D      *fH21_bin;               
    TH1D      *fH21Toward_bin;         
    TH1D      *fH21Transverse_min_bin; 
    TH1D      *fH21Transverse_max_bin; 
    TH1D      *fH21Away_bin;           

    TH1D      *fH21_bin_wojet1;               
    TH1D      *fH21Toward_bin_wojet1;         
    TH1D      *fH21Transverse_min_bin_wojet1; 
    TH1D      *fH21Transverse_max_bin_wojet1; 
    TH1D      *fH21Away_bin_wojet1;           

    TProfile  *fH22;               
    TProfile  *fH22Toward;         
    TProfile  *fH22Transverse_min; 
    TProfile  *fH22Transverse_max; 
    TProfile  *fH22Away;           

    TH1D  *fH22_bin;               
    TH1D  *fH22Toward_bin;         
    TH1D  *fH22Transverse_min_bin; 
    TH1D  *fH22Transverse_max_bin; 
    TH1D  *fH22Away_bin;           

    TH1D  *fH22_bin_wojet1;               
    TH1D  *fH22Toward_bin_wojet1;         
    TH1D  *fH22Transverse_min_bin_wojet1; 
    TH1D  *fH22Transverse_max_bin_wojet1; 
    TH1D  *fH22Away_bin_wojet1;           

    TH1D      *fH23;               
    TH1D      *fH23Toward;         
    TH1D      *fH23Transverse_min; 
    TH1D      *fH23Transverse_max; 
    TH1D      *fH23Away;           

    TProfile  *fH40;           
    TProfile  *fH40toward;     
    TProfile  *fH40away;       
    TProfile  *fH40transmin;   
    TProfile  *fH40transmax;   

    TH1D  *fH40_bin;           
    TH1D  *fH40toward_bin;     
    TH1D  *fH40away_bin;       
    TH1D  *fH40transmin_bin;   
    TH1D  *fH40transmax_bin;   

    TProfile  *fH41;           
    TProfile  *fH41toward;     
    TProfile  *fH41away;       
    TProfile  *fH41transmin;   
    TProfile  *fH41transmax;   

    TH1D  *fH41_bin;           
    TH1D  *fH41toward_bin;     
    TH1D  *fH41away_bin;       
    TH1D  *fH41transmin_bin;   
    TH1D  *fH41transmax_bin;   

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

    AliEmcalJet               *fJet1;                       
    UInt_t                     fNJets_accepted;             
    UInt_t                     fNaccPart;                   
    Int_t                      fNaccClus;                   
    Double_t                   fEvPt;                       

    std::vector<Int_t>         fJet1_sorted_idxvec;         
    std::vector<Int_t>         fEvent_sorted_idxvec;        

  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 ; }
      };

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

    ClassDef ( AliAnalysisTaskEmcalJetCDFUE, 0 );

  };
#endif

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