Changes needed to properly add and subtract the SDD low threshold with data compresse...

[u/mrichter/AliRoot.git] / PWG1 / AliComparisonRes.h

#ifndef ALICOMPARISONRES_H\r
#define ALICOMPARISONRES_H\r
\r
//------------------------------------------------------------------------------\r
// Class to keep information from comparison of \r
// reconstructed and MC particle tracks (TPC resolution).   \r
// \r
// Author: J.Otwinowski 04/02/2008 \r
//------------------------------------------------------------------------------\r
\r
class TFile;\r
class AliMCInfo;\r
class AliESDRecInfo;\r
class AliESDEvent; \r
class AliESD;\r
class AliESDfriend;\r
class AliMCInfoCuts;\r
class AliRecInfoCuts;\r
class TH1I;\r
class TH3F;\r
class TH3;\r
class TProfile;\r
class TProfile2D;\r
class AliESDVertex;\r
\r
#include "TNamed.h"\r
\r
class AliComparisonRes : public TNamed {\r
public :\r
  AliComparisonRes(); \r
  virtual ~AliComparisonRes();\r
  void      InitHisto();\r
  void      InitCuts();\r
  void      Exec(AliMCInfo* infoMC, AliESDRecInfo *infoRC);\r
  void      ProcessConstrained(AliMCInfo* infoMC, AliESDRecInfo *infoRC);\r
  void      Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC);\r
\r
  // Selection cuts\r
  void SetAliRecInfoCuts(AliRecInfoCuts* cuts=0) {fCutsRC = cuts;}   \r
  void SetAliMCInfoCuts(AliMCInfoCuts* cuts=0) {fCutsMC = cuts;}  \r
   \r
  AliRecInfoCuts*  GetAliRecInfoCuts() const {return fCutsRC;}  \r
  AliMCInfoCuts*   GetAliMCInfoCuts()  const {return fCutsMC;}  \r
\r
  // Merge output objects (needed by PROOF) \r
  virtual Long64_t Merge(TCollection* list);\r
\r
  // Analyse output histograms\r
  void Analyse();\r
  static TH1F*       MakeResol(TH2F * his, Int_t integ, Bool_t type); \r
\r
\r
private:\r
  //\r
  // Control histograms\r
  //\r
  TH2F* fPtResolLPT;        //-> pt resolution - low pt\r
  TH2F* fPtResolHPT;        //-> pt resolution - high pt \r
  TH2F* fPtPullLPT;         //-> pt resolution - low pt\r
  TH2F* fPtPullHPT;         //-> pt resolution - high pt \r
\r
  //\r
  // Resolution constrained param\r
  //\r
  TH2F   *fCPhiResolTan;   //-> angular resolution -  constrained\r
  TH2F   *fCTanResolTan;   //-> angular resolution -  constrained\r
  TH2F   *fCPtResolTan;    //-> pt resolution      -  constrained\r
  TH2F   *fCPhiPullTan;    //-> angular resolution -  constrained\r
  TH2F   *fCTanPullTan;    //-> angular resolution -  constrained\r
  TH2F   *fCPtPullTan;     //-> pt resolution      -  constrained\r
\r
  //\r
  // Histograms for track resolution parameterisation\r
  //\r
\r
  TH2F* f1Pt2Resol1PtTPC;      //-> (1/mcpt-1/pt)/(1+1/pt)^2 vs 1/pt (TPC)\r
  TH2F* f1Pt2Resol1PtTPCITS;   //-> (1/mcpt-1/pt)/(1+1/pt)^2 vs 1/pt (TPC+ITS)\r
  TH2F* fYResol1PtTPC;         //-> (mcy-y)/(0.2+1/pt) vs 1/pt (TPC) \r
  TH2F* fYResol1PtTPCITS;      //-> (mcy-y)/(0.2+1/pt) vs 1/pt (TPC + ITS) \r
  TH2F* fZResol1PtTPC;         //-> (mcz-z)/(0.2+1/pt) vs 1/pt (TPC)\r
  TH2F* fZResol1PtTPCITS;      //-> (mcz-z)/(0.2+1/pt) vs 1/pt (TPC+ITS)\r
  TH2F* fPhiResol1PtTPC;       //-> (mcphi-phi)/(0.1+1/pt) vs 1/pt (TPC)\r
  TH2F* fPhiResol1PtTPCITS;    //-> (mcphi-phi)/(0.1+1/pt) vs 1/pt (TPC+ITS)\r
  TH2F* fThetaResol1PtTPC;     //-> (mctheta-theta)/(0.1+1/pt) vs 1/pt (TPC)\r
  TH2F* fThetaResol1PtTPCITS;  //-> (mctheta-theta)/(0.1+1/pt) vs 1/pt (TPC+ITS)\r
  \r
  // constrained\r
  TH2F* fC1Pt2Resol1PtTPC;      //-> (1/mcpt-1/pt)/(1+1/pt)^2 vs 1/pt (TPC)\r
  TH2F* fC1Pt2Resol1PtTPCITS;   //-> (1/mcpt-1/pt)/(1+1/pt)^2 vs 1/pt (TPC+ITS)\r
  TH2F* fCYResol1PtTPC;         //-> (mcy-y)/(0.2+1/pt) vs 1/pt (TPC) \r
  TH2F* fCYResol1PtTPCITS;      //-> (mcy-y)/(0.2+1/pt) vs 1/pt (TPC + ITS) \r
  TH2F* fCZResol1PtTPC;         //-> (mcz-z)/(0.2+1/pt) vs 1/pt (TPC)\r
  TH2F* fCZResol1PtTPCITS;      //-> (mcz-z)/(0.2+1/pt) vs 1/pt (TPC+ITS)\r
  TH2F* fCPhiResol1PtTPC;       //-> (mcphi-phi)/(0.1+1/pt) vs 1/pt (TPC)\r
  TH2F* fCPhiResol1PtTPCITS;    //-> (mcphi-phi)/(0.1+1/pt) vs 1/pt (TPC+ITS)\r
  TH2F* fCThetaResol1PtTPC;     //-> (mctheta-theta)/(0.1+1/pt) vs 1/pt (TPC)\r
  TH2F* fCThetaResol1PtTPCITS;  //-> (mctheta-theta)/(0.1+1/pt) vs 1/pt (TPC+ITS)\r
\r
  AliESDVertex *fVertex;  //! \r
\r
  // Global cuts objects\r
  AliRecInfoCuts*  fCutsRC;      // selection cuts for reconstructed tracks\r
  AliMCInfoCuts*  fCutsMC;       // selection cuts for MC tracks\r
\r
  AliComparisonRes(const AliComparisonRes&); // not implemented\r
  AliComparisonRes& operator=(const AliComparisonRes&); // not implemented\r
\r
  ClassDef(AliComparisonRes,1);\r
};\r
\r
#endif\r