[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliHFSystErr.h

#ifndef ALIHFSYSTERR_H
#define ALIHFSYSTERR_H
/* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

//***********************************************************
// Class AliRDHFSystErr
// to handle systematic errors for charm hadrons
// Author: A.Dainese, andrea.dainese@pd.infn.it
//***********************************************************

#include <TNamed.h>
#include <TH1F.h>
#include "TGraphAsymmErrors.h"


class AliHFSystErr : public TNamed 
{
 public:

  AliHFSystErr(const Char_t* name="HFSystErr", const Char_t* title="");
  AliHFSystErr(Int_t decay,const Char_t* name="HFSystErr", const Char_t* title="");
  
  virtual ~AliHFSystErr();
  
  void DrawErrors(TGraphAsymmErrors *grErrFeeddown=0) const; 

  Double_t GetNormErr() const {return (fNorm ? fNorm->GetBinContent(0) : 0.);}
  Double_t GetBRErr() const {return (fBR ? fBR->GetBinContent(0) : 0.);}
  Double_t GetCutsEffErr(Double_t pt) const;
  Double_t GetMCPtShapeErr(Double_t pt) const;
  Double_t GetSeleEffErr(Double_t pt) const;
  Double_t GetPartAntipartErr(Double_t pt) const;
  Double_t GetPIDEffErr(Double_t pt) const;
  Double_t GetRawYieldErr(Double_t pt) const;
  Double_t GetTrackingEffErr(Double_t pt) const;
  Double_t GetTotalSystErr(Double_t pt,Double_t feeddownErr=0) const;


 private:

  AliHFSystErr(const AliHFSystErr& source);
  AliHFSystErr& operator=(const AliHFSystErr& source); 
 

  void InitD0toKpi();
  void InitDplustoKpipi();
  void InitDstartoD0pi();

  TH1F* ReflectHisto(TH1F *hin) const;

  TH1F *fNorm;            // normalization
  TH1F *fRawYield;        // raw yield 
  TH1F *fTrackingEff;     // tracking efficiency
  TH1F *fBR;              // branching ratio
  TH1F *fCutsEff;         // cuts efficiency
  TH1F *fPIDEff;          // PID efficiency
  TH1F *fMCPtShape;       // MC dNdpt
  TH1F *fPartAntipart;    // particle=antiparticle

  ClassDef(AliHFSystErr,1);  // class for systematic errors of charm hadrons
};

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