added ntuple

[u/mrichter/AliRoot.git] / PWG4 / GammaConv / AliAnalysisTaskGammaJet.h

#ifndef AliAnalysisTaskGammaJet_cxx\r
#define AliAnalysisTaskGammaJet_cxx\r
\r
// example of an analysis task creating a p_t spectrum\r
// Authors: Panos Cristakoglou, Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing\r
\r
class TH1F;\r
class AliESDEvent;\r
class AliGammaConversionAODObject;\r
class AliAODConversionParticle;\r
class AliAODPWG4ParticleCorrelation;\r
class AliAODPWG4Particle;\r
class TClonesArray;\r
class TString;\r
\r
#include "AliAnalysisTaskSE.h"\r
\r
class AliAnalysisTaskGammaJet : public AliAnalysisTaskSE {\r
\r
public:\r
  \r
  AliAnalysisTaskGammaJet(); \r
  AliAnalysisTaskGammaJet(const char *name);\r
  virtual ~AliAnalysisTaskGammaJet();\r
  \r
  virtual void   UserCreateOutputObjects();\r
  virtual void   UserExec(Option_t *option);\r
  virtual void   Terminate(Option_t *);\r
\r
  void SetDeltaAODFileName(TString string) { fDeltaAODFileName = string;}\r
 \r
  inline Float_t GetMinPt() const { return fMinPt;} \r
  void SetMinPt( const Float_t pt ) { fMinPt = pt; }\r
  inline Float_t GetConeSize () const { return fConeSize; }\r
  void SetConeSize ( const Float_t cs ) { fConeSize = cs; }\r
  inline Float_t GetPtThreshold () const { return fPtThreshold; }\r
  void SetPtThreshold ( Float_t ptt ) { fPtThreshold = ptt; }\r
 \r
  \r
 private:\r
\r
  //Clean up\r
  void CleanUp();\r
\r
  //Get the AOD event from whereever it might be accessible\r
  AliAODEvent * GetAODEvent();\r
\r
  //Get Conversion gammas branch\r
  TClonesArray * GetConversionGammas(const AliAODEvent * aodEvent);\r
\r
  //Create PWG4 particles from the aod objects\r
  AliAODPWG4ParticleCorrelation * PWG4PartFromGammaConvAODObject(AliGammaConversionAODObject * gcObject, TString detector);\r
\r
  //Fill AOD tree with PWG4 particles\r
  AliAODPWG4ParticleCorrelation * AddToAOD(AliGammaConversionAODObject * aodO, TClonesArray * branch, TString detector);\r
  AliAODPWG4ParticleCorrelation * AddToAOD(AliAODConversionParticle * aodO, TClonesArray * branch, TString detector);\r
  \r
  //Is particle isolated\r
  Bool_t IsIsolated( AliAODPWG4Particle * particle, TClonesArray * tracks, Float_t coneSize, Float_t ptThreshold);\r
\r
  //Process conv gamma\r
  void ProcessConvGamma( const AliAODEvent * const aodEvent );\r
\r
  //Process calorimeters\r
  void ProcessCalorimeters( const AliAODEvent * const aodEvent );\r
  \r
  //Correlate particle with jets\r
  void CorrelateWithJets(AliAODPWG4ParticleCorrelation * photon, const TClonesArray * const jets);\r
  void CorrelateWithJets(AliAODPWG4Particle * photon, const TClonesArray * const jets, Bool_t const isolated);\r
  void CorrelateWithHadrons(AliAODPWG4Particle * photon, const TClonesArray * tracks, Bool_t const isolated);\r
\r
  //Is eta - phi distance smaller than conesize ?\r
  inline Bool_t IsInCone(Float_t dEta, Float_t dPhi, Float_t coneSize) {   \r
    return ( (dEta*dEta + dPhi*dPhi) < coneSize*coneSize);\r
  }\r
\r
  TList       *fOutputList; //! Output list\r
  TH1F        *fHistPt; //! Pt spectrum\r
  TH1F        *fHistPtPhos; //! Pt spectrum\r
  TH1F        *fHistPtEmcal; //! Pt spectrum\r
\r
  TH1F        *fHistPhotPhi;\r
  TH1F        *fHistHadPhi;\r
  TH1F        *fHistJetPhi;\r
\r
\r
  TH1F        *fHistPtJets; //! Pt spectrum\r
  TH1F        *fHistGammaJets; //!Phi correlations\r
  TH1F        *fHistGammaJetsIso; //!Phi correlations\r
  TH1F        *fHistMaxdPhi; //!Phi correlations\r
  TH1F        *fHistMaxdPhiIso; //!Phi correlations\r
  TH1F        *fHistMaxdPhiIsoPt; //!Phi correlations\r
\r
  TH1F        *fHadHistPt; //! Pt spectrum\r
  TH1F        *fHadHistdPhi; //!Phi correlations\r
  TH1F        *fHadHistdPhiIso; //!Phi correlations\r
  TH1F        *fHadHistMaxdPhi; //!Phi correlations\r
  TH1F        *fHadHistMaxdPhiIso; //!Phi correlations\r
  TH1F        *fHadHistMaxdPhiIsoPt; //!Phi correlations\r
  \r
  \r
  Float_t fMinPt; //Minimum pt for correlation\r
  Float_t fConeSize; //cone size for isolation\r
  Float_t fPtThreshold; //Threshold pt for isolation\r
\r
  TString     fDeltaAODFileName;//! File where Gamma Conv AOD is located, if not in default AOD\r
\r
  TClonesArray * fPhotons;\r
\r
  AliAnalysisTaskGammaJet(const AliAnalysisTaskGammaJet&); // not implemented\r
  AliAnalysisTaskGammaJet& operator=(const AliAnalysisTaskGammaJet&); // not implemented\r
  \r
  ClassDef(AliAnalysisTaskGammaJet, 2); // example of analysis\r
};\r
\r
#endif\r