PWG3/vertexingHF/AliHFSystErr.h

   1 #ifndef ALIHFSYSTERR_H
   2 #define ALIHFSYSTERR_H
   3 /* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //***********************************************************
   9 // Class AliRDHFSystErr
  10 // to handle systematic errors for charm hadrons
  11 // Author: A.Dainese, andrea.dainese@pd.infn.it
  12 //***********************************************************
  13
  14 #include <TNamed.h>
  15 #include <TH1F.h>
  16 #include "TGraphAsymmErrors.h"
  17
  18
  19 class AliHFSystErr : public TNamed
  20 {
  21  public:
  22
  23   AliHFSystErr(const Char_t* name="HFSystErr", const Char_t* title="");
  24   AliHFSystErr(Int_t decay,const Char_t* name="HFSystErr", const Char_t* title="");
  25
  26   virtual ~AliHFSystErr();
  27
  28   void DrawErrors(TGraphAsymmErrors *grErrFeeddown=0) const;
  29
  30   Double_t GetNormErr() const {return (fNorm ? fNorm->GetBinContent(0) : 0.);}
  31   Double_t GetBRErr() const {return (fBR ? fBR->GetBinContent(0) : 0.);}
  32   Double_t GetCutsEffErr(Double_t pt) const;
  33   Double_t GetMCPtShapeErr(Double_t pt) const;
  34   Double_t GetSeleEffErr(Double_t pt) const;
  35   Double_t GetPartAntipartErr(Double_t pt) const;
  36   Double_t GetPIDEffErr(Double_t pt) const;
  37   Double_t GetRawYieldErr(Double_t pt) const;
  38   Double_t GetTrackingEffErr(Double_t pt) const;
  39   Double_t GetTotalSystErr(Double_t pt,Double_t feeddownErr=0) const;
  40
  41
  42  private:
  43
  44   AliHFSystErr(const AliHFSystErr& source);
  45   AliHFSystErr& operator=(const AliHFSystErr& source);
  46
  47
  48   void InitD0toKpi();
  49   void InitDplustoKpipi();
  50   void InitDstartoD0pi();
  51
  52   TH1F* ReflectHisto(TH1F *hin) const;
  53
  54   TH1F *fNorm;            // normalization
  55   TH1F *fRawYield;        // raw yield
  56   TH1F *fTrackingEff;     // tracking efficiency
  57   TH1F *fBR;              // branching ratio
  58   TH1F *fCutsEff;         // cuts efficiency
  59   TH1F *fPIDEff;          // PID efficiency
  60   TH1F *fMCPtShape;       // MC dNdpt
  61   TH1F *fPartAntipart;    // particle=antiparticle
  62
  63   ClassDef(AliHFSystErr,1);  // class for systematic errors of charm hadrons
  64 };
  65
  66 #endif
  67