Changes to fast embedding and jetresponse: QA plot for delta AOD, Jet eta phi selecti...

[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliAnalysisTaskMinijet.h

#ifndef ALIANALYSISTASKMINIJET_H\r
#define ALIANALYSISTASKMINIJET_H\r
\r
// analysis task performing mini jet analysis (use ESD or AOD as input)\r
// Author: eva.sicking@cern.ch\r
\r
class TList;\r
class TH1F;\r
class TH2F;\r
class TProfile;\r
\r
class AliESDtrackCuts;\r
\r
\r
#include "AliAnalysisTaskSE.h"\r
\r
class AliAnalysisTaskMinijet : public AliAnalysisTaskSE {\r
 public:\r
  AliAnalysisTaskMinijet(const char *name="<default name>");\r
  virtual ~AliAnalysisTaskMinijet();\r
  \r
  virtual void   UserCreateOutputObjects();\r
  virtual void   UserExec(Option_t* option);\r
  virtual void   Terminate(Option_t *);\r
  \r
  Int_t LoopESD  (Float_t **pt, Float_t **eta, Float_t **phi, Int_t **nTracksTracklets);\r
  Int_t LoopESDMC(Float_t **pt, Float_t **eta, Float_t **phi, Int_t **nTracksTracklets);\r
  Int_t LoopAOD  (Float_t **pt, Float_t **eta, Float_t **phi, Int_t **nTracksTracklets);\r
  Int_t LoopAODMC(Float_t **pt, Float_t **eta, Float_t **phi, Int_t **nTracksTracklets);\r
  void  Analyse  (Float_t* pt, Float_t* eta, Float_t* phi, Int_t ntacks, Int_t ntacklets=0,Int_t nAll=0, Int_t mode=0);\r
  void  CleanArrays(Float_t* pt, Float_t* eta, Float_t* phi, Int_t* nTracksTracklets=0);\r
  Bool_t SelectParticlePlusCharged(Short_t charge, Int_t pdg, Bool_t prim);\r
  Bool_t SelectParticle(Short_t charge, Int_t pdg, Bool_t prim);\r
\r
\r
  void UseMC(Bool_t useMC=kTRUE, Bool_t mcOnly=kFALSE)    {fUseMC = useMC; fMcOnly=mcOnly;}\r
\r
  virtual void   SetCuts(AliESDtrackCuts* cuts)           {fCuts = cuts;}\r
\r
  void   SetRadiusCut(Float_t radiusCut)           {fRadiusCut = radiusCut;}  \r
  void   SetTriggerPtCut(Float_t triggerPtCut)     {fTriggerPtCut = triggerPtCut;}  \r
  void   SetAssociatePtCut(Float_t associatePtCut) {fAssociatePtCut = associatePtCut;} \r
  void   SetEventAxis(Int_t leadingOrRandom)       {fLeadingOrRandom = leadingOrRandom;}  \r
  void   SetMode(Int_t mode)                       {fMode = mode;}\r
  void   SetMaxVertexZ(Float_t vertexZCut)         {fVertexZCut = vertexZCut;}\r
  void   SetMaxEta(Float_t etaCut)                 {fEtaCut = etaCut;}\r
  void   SetMaxEtaSeed(Float_t etaCutSeed)         {fEtaCutSeed = etaCutSeed;}\r
\r
  void   SelectParticles(Int_t selectParticles)    {fSelectParticles = selectParticles;}\r
\r
\r
 private:\r
\r
  Bool_t       fUseMC;\r
  Bool_t       fMcOnly;\r
  AliESDtrackCuts* fCuts;                   // List of cuts for ESDs\r
  Float_t      fRadiusCut;                  // radius cut \r
  Float_t      fTriggerPtCut;               // cut on particle pt used as event axis\r
  Float_t      fAssociatePtCut;             // cut on particle pt used for correlations\r
  Int_t        fLeadingOrRandom;            // event axis:leading track or random track\r
  Int_t        fMode;                       // ESD(=0) of AOD(=1) reading \r
  Float_t      fVertexZCut;                 // vertex cut\r
  Float_t      fEtaCut;                     // eta acceptance cut\r
  Float_t      fEtaCutSeed;                 // eta acceptance cut for seed\r
  Int_t        fSelectParticles;            // only in cas of MC: use also neutral particles or not \r
\r
  AliESDEvent *fESDEvent;                   //! esd event\r
  AliAODEvent *fAODEvent;                   //! aod event\r
  \r
  TList	      *fHists;                      // output list\r
  TH1F        *fHistPt;                     // Pt spectrum ESD\r
  TH1F        *fHistPtMC;                   // Pt spectrum MC\r
  TH2F        *fChargedPi0;                 // charged versus charged+Pi0\r
  TH1F       * fVertexZ[4];                 // z of vertex\r
\r
  TH1F       * fPt[4];                      // pt\r
  TH1F       * fEta[4];                     // et\r
  TH1F       * fPhi[4];                     // phi\r
\r
  TH1F       * fPtSeed[4];                  // pt of seed (event axis)\r
  TH1F       * fEtaSeed[4];                 // eta of seed \r
  TH1F       * fPhiSeed[4];                 // phi of seed\r
\r
  TH1F       * fPtOthers[4];                // pt of all other particels used in dEtadPhi\r
  TH1F       * fEtaOthers[4];               // eta of all other particels used in dEtadPhi\r
  TH1F       * fPhiOthers[4];               // phi of all other particels used in dEtadPhi\r
  TH2F       * fPtEtaOthers[4];             // pt-eta of all other particels used in dEtadPhi\r
\r
\r
  TH2F       * fPhiEta[4];                  // eta - phi\r
  TH2F       * fDPhiDEtaEventAxis[4];       // correlation dEta-dPhi towards event axis\r
  TH1F       * fTriggerNch[4];              // number of triggers with accepted-track number\r
  TH1F       * fTriggerTracklet[4];         // number of triggers with accepted-tracklet number\r
  TH2F       * fNch07Nch[4];                // nCharged with pT>fTriggerPtCut vs nCharged\r
  TProfile   * pNch07Nch[4];                // nCharged with pT>fTriggerPtCut vs nCharged\r
  TH2F       * fNch07Tracklet[4];           // nCharged with pT>fTriggerPtCut vs nTracklet\r
  TH2F       * fNchTracklet[4];             // nCharged vs nTracklet\r
  TProfile   * pNch07Tracklet[4];           // nCharged with pT>fTriggerPtCut vs nTracklet\r
\r
  TH1F       * fDPhiEventAxisNchBin[4][150];// delta phi of associate tracks to event axis per Nch bin\r
  TH1F       * fDPhiEventAxisNchBinTrig[4][150];// "" for all possoble trigger particles\r
\r
  TH1F       * fDPhiEventAxisTrackletBin[4][150]; // delta phi of associate tracks to event axis per Nch bin\r
  TH1F       * fDPhiEventAxisTrackletBinTrig[4][150]; // "" for all possible trigger particles\r
\r
  AliAnalysisTaskMinijet(const AliAnalysisTaskMinijet&); // not implemented\r
  AliAnalysisTaskMinijet& operator=(const AliAnalysisTaskMinijet&); // not implemented\r
  \r
  ClassDef(AliAnalysisTaskMinijet, 1); // example of analysis\r
};\r
\r
#endif\r