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

#ifndef ALIANALYSISTASKCOMPAREHF_H
#define ALIANALYSISTASKCOMPAREHF_H

/* Copyright(c) 1998-2008, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

//*************************************************************************
// Class AliAnalysisTaskSECompareHF
// AliAnalysisTaskSE for the comparison of heavy-flavour decay candidates
// to MC truth (kinematics stored in the AOD)
// Author: A.Dainese, andrea.dainese@lnl.infn.it
//*************************************************************************

#include <TNtuple.h>
#include <TH1F.h>

#include "AliAnalysisTaskSE.h"

class AliAnalysisTaskSECompareHF : public AliAnalysisTaskSE
{
 public:

  AliAnalysisTaskSECompareHF();
  AliAnalysisTaskSECompareHF(const char *name);
  virtual ~AliAnalysisTaskSECompareHF();


  // Implementation of interface methods
  virtual void UserCreateOutputObjects();
  virtual void Init();
  virtual void LocalInit() {Init();}
  virtual void UserExec(Option_t *option);
  virtual void Terminate(Option_t *option);

  void SetD0toKpiCuts(Double_t cut0=1000.,Double_t cut1=100000.,
		      Double_t cut2=1.1,Double_t cut3=0.,Double_t cut4=0.,
		      Double_t cut5=100000.,Double_t cut6=100000.,
		      Double_t cut7=100000000.,Double_t cut8=-1.1); 
  // cuts[0] = inv. mass half width [GeV]   
  // cuts[1] = dca [cm]
  // cuts[2] = cosThetaStar 
  // cuts[3] = pTK [GeV/c]
  // cuts[4] = pTPi [GeV/c]
  // cuts[5] = d0K [cm]   upper limit!
  // cuts[6] = d0Pi [cm]  upper limit!
  // cuts[7] = d0d0 [cm^2]
  // cuts[8] = cosThetaPoint
  void SetD0toKpiCuts(const Double_t cuts[9]); 
  
 private:

  AliAnalysisTaskSECompareHF(const AliAnalysisTaskSECompareHF &source);
  AliAnalysisTaskSECompareHF& operator=(const AliAnalysisTaskSECompareHF& source); 
  TList   *fOutput; //! list send on output slot 0
  TNtuple *fNtupleD0Cmp; // output ntuple
  TH1F    *fHistMass;    // output histogram
  Double_t fD0toKpiCuts[9]; // cuts for D0->Kpi selection
  
  ClassDef(AliAnalysisTaskSECompareHF,1); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates
};

#endif
Commit	Line	Data
3a219f60	1	#ifndef ALIANALYSISTASKCOMPAREHF_H
	2	#define ALIANALYSISTASKCOMPAREHF_H
	3
	4	/* Copyright(c) 1998-2008, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	//*************************************************************************
	8	// Class AliAnalysisTaskSECompareHF
	9	// AliAnalysisTaskSE for the comparison of heavy-flavour decay candidates
	10	// to MC truth (kinematics stored in the AOD)
	11	// Author: A.Dainese, andrea.dainese@lnl.infn.it
	12	//*************************************************************************
	13
	14	#include <TNtuple.h>
	15	#include <TH1F.h>
	16
	17	#include "AliAnalysisTaskSE.h"
	18
	19	class AliAnalysisTaskSECompareHF : public AliAnalysisTaskSE
	20	{
	21	public:
	22
	23	AliAnalysisTaskSECompareHF();
	24	AliAnalysisTaskSECompareHF(const char *name);
	25	virtual ~AliAnalysisTaskSECompareHF();
	26
	27
	28	// Implementation of interface methods
	29	virtual void UserCreateOutputObjects();
	30	virtual void Init();
	31	virtual void LocalInit() {Init();}
	32	virtual void UserExec(Option_t *option);
	33	virtual void Terminate(Option_t *option);
	34
	35	void SetD0toKpiCuts(Double_t cut0=1000.,Double_t cut1=100000.,
	36	Double_t cut2=1.1,Double_t cut3=0.,Double_t cut4=0.,
	37	Double_t cut5=100000.,Double_t cut6=100000.,
	38	Double_t cut7=100000000.,Double_t cut8=-1.1);
	39	// cuts[0] = inv. mass half width [GeV]
	40	// cuts[1] = dca [cm]
	41	// cuts[2] = cosThetaStar
	42	// cuts[3] = pTK [GeV/c]
	43	// cuts[4] = pTPi [GeV/c]
	44	// cuts[5] = d0K [cm] upper limit!
	45	// cuts[6] = d0Pi [cm] upper limit!
	46	// cuts[7] = d0d0 [cm^2]
	47	// cuts[8] = cosThetaPoint
	48	void SetD0toKpiCuts(const Double_t cuts[9]);
	49
	50	private:
	51
	52	AliAnalysisTaskSECompareHF(const AliAnalysisTaskSECompareHF &source);
	53	AliAnalysisTaskSECompareHF& operator=(const AliAnalysisTaskSECompareHF& source);
	54	TList *fOutput; //! list send on output slot 0
	55	TNtuple *fNtupleD0Cmp; // output ntuple
	56	TH1F *fHistMass; // output histogram
	57	Double_t fD0toKpiCuts[9]; // cuts for D0->Kpi selection
	58
	59	ClassDef(AliAnalysisTaskSECompareHF,1); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates
	60	};
	61
	62	#endif
	63