AliAnalysisTaskRhoBase.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKRHOBASE_H
#define ALIANALYSISTASKRHOBASE_H

// $Id$

class TString;
class TF1;
class TH1F;
class TH2F;
class TH3F;
class AliRhoParameter;

#include "AliAnalysisTaskEmcalJet.h"

class AliAnalysisTaskRhoBase : public AliAnalysisTaskEmcalJet {
 public:
  AliAnalysisTaskRhoBase();
  AliAnalysisTaskRhoBase(const char *name, Bool_t histo=kFALSE);
  virtual ~AliAnalysisTaskRhoBase() {}

  void                   UserCreateOutputObjects();

  void                   SetOutRhoName(const char *name)                       { fOutRhoName           = name    ;
                                                                                 fOutRhoScaledName     = Form("%s_Scaled",name);     }
  void                   SetCompareRhoName(const char *name)                   { fCompareRhoName       = name    ;                   }
  void                   SetCompareRhoScaledName(const char *name)             { fCompareRhoScaledName = name    ;                   }
  void                   SetScaleFunction(TF1* sf)                             { fScaleFunction        = sf      ;                   }
  void                   SetRhoFunction(TF1* rf)                               { fRhoFunction          = rf      ;                   }
  TF1*                   LoadRhoFunction(const char* path, const char* name);
  void                   SetInEventSigmaRho(Double_t s)                        { fInEventSigmaRho      = s       ;                   }
  void                   SetAttachToEvent(Bool_t a)                            { fAttachToEvent        = a       ;                   }
  void                   SetSmallSystem(Bool_t setter = kTRUE)                 { fIsPbPb               = !setter ;                   }

  const char*            GetOutRhoName() const                                 { return fOutRhoName.Data()       ;                   }
  const char*            GetOutRhoScaledName() const                           { return fOutRhoScaledName.Data() ;                   }

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

  virtual Double_t       GetRhoFactor(Double_t cent);
  virtual Double_t       GetScaleFactor(Double_t cent);

  TString                fOutRhoName;                    // name of output rho object
  TString                fOutRhoScaledName;              // name of output scaled rho object
  TString                fCompareRhoName;                // name of rho object to compare
  TString                fCompareRhoScaledName;          // name of scaled rho object to compare
  TF1                   *fRhoFunction;                   // pre-computed rho as a function of centrality
  TF1                   *fScaleFunction;                 // pre-computed scale factor as a function of centrality
  Double_t               fInEventSigmaRho;               // in-event sigma rho
  Bool_t                 fAttachToEvent;                 // whether or not attach rho to the event objects list
  Bool_t                 fIsPbPb;                        // different histogram ranges for pp/pPb and PbPb
  
  AliRhoParameter       *fOutRho;                        
  AliRhoParameter       *fOutRhoScaled;                  
  AliRhoParameter       *fCompareRho;                    
  AliRhoParameter       *fCompareRhoScaled;              

  TH2F                  *fHistJetPtvsCent;               
  TH2F                  *fHistJetAreavsCent;             
  TH2F                  *fHistJetRhovsCent;              
  TH2F                  *fHistNjetvsCent;                
  TH2F                  *fHistJetPtvsNtrack;             
  TH2F                  *fHistJetAreavsNtrack;           
  TH2F                  *fHistNjetvsNtrack;              
  TH2F                  *fHistNjUEoverNjVsNj[12];        
  TH2F                  *fHistJetNconstVsPt[4];          
  TH2F                  *fHistJetRhovsEta[4];               
  TH2F                  *fHistRhovsCent;                 
  TH2F                  *fHistRhoScaledvsCent;           
  TH2F                  *fHistDeltaRhovsCent;            
  TH2F                  *fHistDeltaRhoScalevsCent;       

  TH3F                  *fHistRhovsNtrackvsV0Mult;       
  TH3F                  *fHistRhoScaledvsNtrackvsV0Mult; 
  TH2F                  *fHistDeltaRhovsNtrack;          
  TH2F                  *fHistDeltaRhoScalevsNtrack;     
 
  TH2F                  *fHistRhovsNcluster;             
  TH2F                  *fHistRhoScaledvsNcluster;       

  AliAnalysisTaskRhoBase(const AliAnalysisTaskRhoBase&);             // not implemented
  AliAnalysisTaskRhoBase& operator=(const AliAnalysisTaskRhoBase&);  // not implemented
  
  ClassDef(AliAnalysisTaskRhoBase, 11); // Rho base task
};
#endif