AliAnalysisTaskSoftDrop.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSOFTDROP_H
#define ALIANALYSISTASKSOFTDROP_H

// $Id$

class TH1;
class TH2;
class TH3;
class AliJetContainer;
class AliParticleContainer;
class AliClusterContainer;

#include "AliAnalysisTaskEmcalJet.h"
#include "FJ_includes.h"

class AliAnalysisTaskSoftDrop : public AliAnalysisTaskEmcalJet {
 public:

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

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

  static Float_t              SoftDropDeclustering(fastjet::PseudoJet jet, const Float_t zcut, const Float_t beta);

  static AliAnalysisTaskSoftDrop* AddTaskSoftDrop(
    const char *ntracks            = "usedefault",
    const char *nclusters          = "usedefault",
    const char *njets              = "Jets",
    const char *nrho               = "Rho",
    Int_t       nCentBins          = 1,
    Double_t    jetradius          = 0.2,
    Double_t    jetptcut           = 1,
    Double_t    jetareacut         = 0.6,
    const char *type               = "EMCAL",
    Int_t       leadhadtype        = 0,
    const char *taskname           = "AliAnalysisTaskSoftDrop"
  );

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

  void                        SoftDropDeepDeclustering(fastjet::PseudoJet jet, const Float_t inpt); 

  // General histograms
  TH1                       **fHistTracksPt;            
  TH1                       **fHistClustersPt;          
  TH1                       **fHistLeadingJetPt;        
  TH2                       **fHistJetsPhiEta;          
  TH2                       **fHistJetsPtArea;          
  TH2                       **fHistJetsPtLeadHad;       
  TH2                       **fHistJetsCorrPtArea;      
  TH3                        *fHistPtDEtaDPhiTrackClus; 
  TH3                        *fHistPtDEtaDPhiClusTrack; 
  TH1                       **fHistNTracks;             

  TH1                        *fHistClustDx;             
  TH1                        *fHistClustDz;             
  TH1                        *fNAccJets;                
  TH1                        *fhZg;                     
  TH2                        *fhCorrPtZg;               
  TH2                        *fhCorrPtRg;               
  TH2                        *fhCorrPtPtfrac;           
  TH2                        *fhCorrPtDropCount;        

  TH2                        *fhCorrPtZg2;               
  TH2                        *fhCorrPtZgD;               

  TH3                        *fhCorrPtZgRg;              
  TH3                        *fhCorrPtZgSDstep;          
  TH3                        *fhCorrPtRgSDstep;          

  TH2                        *fhCorrPtRgD;               

  AliJetContainer            *fJetsCont;                   
  AliParticleContainer       *fTracksCont;                 
  AliClusterContainer        *fCaloClustersCont;           

 private:
  AliAnalysisTaskSoftDrop(const AliAnalysisTaskSoftDrop&);            // not implemented
  AliAnalysisTaskSoftDrop &operator=(const AliAnalysisTaskSoftDrop&); // not implemented

  Int_t                       fSDM;                     

  ClassDef(AliAnalysisTaskSoftDrop, 1) // jet sample analysis task
};
#endif