3 /* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //***********************************************************
9 // Class AliRDHFSystErr
10 // to handle systematic errors for charm hadrons
11 // Author: A.Dainese, andrea.dainese@pd.infn.it
12 //***********************************************************
18 #include "TGraphAsymmErrors.h"
21 class AliHFSystErr : public TNamed
25 AliHFSystErr(const Char_t* name="HFSystErr", const Char_t* title="");
27 virtual ~AliHFSystErr();
29 void DrawErrors(TGraphAsymmErrors *grErrFeeddown=0) const;
31 Double_t GetNormErr() const {return (fNorm ? fNorm->GetBinContent(0) : 0.);}
32 Double_t GetBRErr() const {return (fBR ? fBR->GetBinContent(0) : 0.);}
33 Double_t GetCutsEffErr(Double_t pt) const;
34 Double_t GetMCPtShapeErr(Double_t pt) const;
35 Double_t GetSeleEffErr(Double_t pt) const;
36 Double_t GetPartAntipartErr(Double_t pt) const;
37 Double_t GetPIDEffErr(Double_t pt) const;
38 Double_t GetRawYieldErr(Double_t pt) const;
39 Double_t GetTrackingEffErr(Double_t pt) const;
40 Double_t GetTotalSystErr(Double_t pt,Double_t feeddownErr=0) const;
42 // Setting the run number
43 // set the two last numbers of the year (is 10 for 2010)
44 void SetRunNumber(Int_t number) {
46 AliInfo(Form(" Settings for run year 20%2d",fRunNumber));
48 // Setting the collision type
50 void SetCollisionType(Int_t type) {
51 fCollisionType = type;
52 if (fCollisionType==0) { AliInfo(" Settings for p-p collisions"); }
53 else if(fCollisionType==1) { AliInfo(" Settings for Pb-Pb collisions"); }
55 // Setting for the centrality class
56 // 0100 for MB, 020 (4080) for 0-20 (40-80) CC and so on
57 void SetCentrality(TString centrality) {
58 fCentralityClass = centrality;
59 AliInfo(Form(" Settings for centrality class %s",fCentralityClass.Data()));
61 void SetIsLowEnergy(Bool_t flag) {
63 AliInfo(" Settings for the low energy run");
66 // Function to initialize the variables/histograms
67 void Init(Int_t decay);
71 AliHFSystErr(const AliHFSystErr& source);
72 AliHFSystErr& operator=(const AliHFSystErr& source);
74 void InitD0toKpi2010pp();
75 void InitDplustoKpipi2010pp();
76 void InitDstartoD0pi2010pp();
78 void InitD0toKpi2010PbPb020();
79 void InitDplustoKpipi2010PbPb020();
80 void InitDstartoD0pi2010PbPb020();
82 void InitD0toKpi2010PbPb4080();
83 void InitDplustoKpipi2010PbPb4080();
84 void InitDstartoD0pi2010PbPb4080();
86 void InitD0toKpi2010ppLowEn();
87 void InitDplustoKpipi2010ppLowEn();
88 void InitDstartoD0pi2010ppLowEn();
90 void InitDstartoD0pi2010PbPb2040();
92 TH1F* ReflectHisto(TH1F *hin) const;
94 TH1F *fNorm; // normalization
95 TH1F *fRawYield; // raw yield
96 TH1F *fTrackingEff; // tracking efficiency
97 TH1F *fBR; // branching ratio
98 TH1F *fCutsEff; // cuts efficiency
99 TH1F *fPIDEff; // PID efficiency
100 TH1F *fMCPtShape; // MC dNdpt
101 TH1F *fPartAntipart; // particle=antiparticle
103 Int_t fRunNumber; // Run Number (year)
104 Int_t fCollisionType; // Collision type: pp=0, PbPb=1
105 TString fCentralityClass; // Centrality class
106 // MB:0100, 0-10:010, 0-20:020 ...40-80:4080...
107 Bool_t fIsLowEnergy; // flag for the low energy (2.76TeV) run
109 ClassDef(AliHFSystErr,3); // class for systematic errors of charm hadrons