make dQdl plot in p bins, reduce PID resolution for track IN

[u/mrichter/AliRoot.git] / PWGPP / TRD / AliTRDrecoTask.h

#ifndef ALITRDRECOTASK_H\r
#define ALITRDRECOTASK_H\r
\r
///////////////////////////////////////////////////////\r
//\r
// Basic class for Performance/Calibration TRD tasks\r
//\r
// Author: Alexandru Bercuci, 10/09/2008\r
//\r
//////////////////////////////////////////////////\r
\r
#ifndef ALIANALYSISTASKSE_H\r
#include "AliAnalysisTaskSE.h"\r
#endif\r
\r
#ifndef ALITRDTRACKINFO_H\r
#include "info/AliTRDtrackInfo.h"\r
#endif\r
\r
#ifndef ALITRDEVENTINFO_H\r
#include "info/AliTRDeventInfo.h"\r
#endif\r
\r
class TAxis;\r
class TH1;\r
class TH2;\r
class TH3;\r
class TF1;\r
class TList;\r
class TObjArray;\r
class TTreeSRedirector;\r
class AliTRDtrackV1;\r
template <typename Value> class TVectorT;\r
typedef class TVectorT<Float_t> TVector;\r
class AliTRDrecoTask : public AliAnalysisTaskSE \r
{\r
friend class AliEveTRDTrackList;\r
public:\r
  enum AliTRDrecoSteeringBits{\r
     kMCdata      = BIT(18)\r
    ,kFriends     = BIT(19)\r
    ,kPostProcess = BIT(20)\r
    ,kHeavyIon    = BIT(21)\r
  };\r
  \r
  class AliTRDrecoProjection : public TNamed\r
  {\r
  public:\r
    AliTRDrecoProjection();\r
    virtual ~AliTRDrecoProjection();\r
    AliTRDrecoProjection& operator+=(const AliTRDrecoProjection& other);\r
    AliTRDrecoProjection& operator=(const AliTRDrecoProjection& other);\r
    void  Build(const Char_t *n, const Char_t *t, Int_t ix, Int_t iy, Int_t iz, TAxis *aa[]);\r
    void  Increment(Int_t bin[], Double_t v);\r
    TH3*  H() const { return fH;}\r
    TH2*  Projection2D(const Int_t nstat, const Int_t ncol, const Int_t mid=0, Bool_t del=kTRUE);\r
    void  SetRebinStrategy(Int_t n, Int_t rebx[], Int_t reby[]);\r
    void  SetShowRange(Float_t zm, Float_t zM, Float_t em=0., Float_t eM=0.) {fRange[0] = zm; fRange[1] = zM; fRange[2] = em; fRange[3] = eM;}\r
  private:\r
    AliTRDrecoProjection(const AliTRDrecoProjection&);\r
  protected:\r
    TH3  *fH;          // data container\r
    Int_t fAx[3];      // projection axes\r
    Int_t fNrebin;     // no. of rebinning steps\r
    Int_t *fRebinX;    //[fNrebin] rebinning of the X axis\r
    Int_t *fRebinY;    //[fNrebin] rebinning of the Y axis\r
    Float_t fRange[4]; //show range of the z processed\r
\r
    ClassDef(AliTRDrecoProjection, 2)  // wrapper for a projection container THnSparse -> TH3\r
  };\r
\r
  AliTRDrecoTask();\r
  AliTRDrecoTask(const char *name, const char *title);\r
  virtual ~AliTRDrecoTask();\r
  \r
  \r
  virtual void   UserCreateOutputObjects();\r
  virtual void   UserExec(Option_t *opt);\r
  virtual void   SetDebugLevel(Int_t level);\r
  \r
    \r
  static Float_t GetMeanStat(TH1 *h, Float_t cut=0., Option_t *opt="");\r
  Int_t          GetNRefFigures() const; \r
  const Char_t*  GetNameId() const       { return fNameId;}\r
  TList*         GetPlotFunctors() const { return fPlotFuncList;}\r
  static Int_t   GetPtBinSignificant(Float_t pt);\r
  virtual Bool_t GetRefFigure(Int_t ifig);\r
  virtual void   MakeSummary();\r
  void           MakeDetectorPlot(Int_t ly=0, const Option_t *opt="eta");\r
  void           MakeDetectorPlotNEW(Int_t ly=0, const Option_t *opt="eta");\r
  Bool_t         IsHeavyIon() const      { return TestBit(kHeavyIon);};\r
  Bool_t         IsPP() const            { return !TestBit(kHeavyIon);};\r
  Bool_t         HasFriends() const      { return TestBit(kFriends);};\r
  Bool_t         HasMCdata() const       { return TestBit(kMCdata);};\r
  Bool_t         HasPostProcess() const  { return TestBit(kPostProcess);};\r
  Bool_t         HasRunTerminate() const { return fRunTerminate; }\r
  virtual TObjArray* Histos()            { return fContainer;}\r
\r
  virtual Bool_t Load(const Char_t *file = "AnalysisResults.root", const Char_t *dir = "TRD_Performance");\r
  virtual Bool_t LoadDetectorMap(const Char_t *file = "AnalysisResults.root", const Char_t *dir = "TRD_Performance");\r
  virtual Bool_t Save(TObjArray * const res);\r
  virtual Bool_t PostProcess();\r
  virtual Bool_t PutTrendValue(const Char_t *name, Double_t val, Double_t err=0);\r
  virtual void   SetFriends(Bool_t fr = kTRUE) {SetBit(kFriends, fr);}\r
  virtual void   SetMCdata(Bool_t mc = kTRUE) {SetBit(kMCdata, mc);}\r
  virtual void   SetNameId(const Char_t *nid) {snprintf(fNameId, 10, "%s", nid);}\r
  virtual void   SetPostProcess(Bool_t pp = kTRUE) {SetBit(kPostProcess, pp);}\r
  static Float_t SetNormZ(TH2 *h2, Int_t bxmin=1, Int_t bxmax=-1, Int_t bymin=1, Int_t bymax=-1, Float_t thr=0.);\r
  static void    SetRangeZ(TH2 *h2, Float_t m, Float_t M, Float_t thr=0.);\r
  void SetRunTerminate(Bool_t runTerminate = kTRUE) { fRunTerminate = runTerminate; }\r
  virtual void   Terminate(Option_t *);\r
\r
protected:\r
  static TTreeSRedirector* DebugStream() { return fgDebugStream;}\r
  virtual void   InitFunctorList();\r
  Bool_t         HasFunctorList() const { return fPlotFuncList != NULL; }\r
\r
  Char_t                fNameId[10];       // unique identifier of task particularity\r
  UChar_t               fNRefFigures;      // no of reference figures reported by task\r
  TObjArray             *fDets;            //! OLD container to store detector position and status support should be discontinued \r
  TVector               *fDetsV;           //! NEW container to store detector position and status\r
  TObjArray             *fContainer;       //! container to store results\r
  AliTRDeventInfo       *fEvent;           //! Event Info\r
  TObjArray             *fTracks;          //! Array of tracks\r
  TObjArray             *fClusters;        //! Array of clusters\r
  const TObjArray       *fkClusters;       //! current detector clusters array\r
  const AliTRDtrackV1   *fkTrack;          //! current track\r
  const AliTRDtrackInfo::AliMCinfo  *fkMC; //! MC info\r
  const AliTRDtrackInfo::AliESDinfo *fkESD;//! ESD info\r
  Char_t                 fSpecies;         //! species index +1 with charge sign\r
  Float_t                fPt;              //! p_t of the track being analyzed\r
  Float_t                fPhi;             //! phi of the track being analyzed\r
  Float_t                fEta;             //! eta of the track being analyzed\r
\r
private:\r
  AliTRDrecoTask(const AliTRDrecoTask&);\r
  AliTRDrecoTask& operator=(const AliTRDrecoTask&);\r
\r
  TList             *fPlotFuncList;        //! track functors list\r
  TList             *fDetFuncList;         //! detector functors list\r
  Bool_t            fRunTerminate;         // Switch for Terminate Function\r
  static TTreeSRedirector *fgDebugStream;  //! Debug stream\r
  static const Int_t fgNPt0 = 4;           // No of significant pt bins \r
  static Float_t fgPt0[fgNPt0];            // Array with limits for significant pt bins \r
\r
  ClassDef(AliTRDrecoTask, 5) // base TRD reconstruction task\r
};\r
\r
#endif\r
\r