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

#ifndef ALICFTASKVERTEXINGHF_H
#define ALICFTASKVERTEXINGHF_H
/**************************************************************************
 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

//-----------------------------------------------------------------------
// Class for HF corrections as a function of many variables and step 
// Author : C. Zampolli, CERN
//			D. Caffarri, Univ & INFN Padova caffarri@pd.infn.it
// Base class for HF Unfolding - agrelli@uu.nl
//-----------------------------------------------------------------------


#include "AliAnalysisTaskSE.h"
#include "AliCFVertexingHF2Prong.h"
#include "AliCFVertexingHF.h"

class TH1I;
class TParticle ;
class TFile ;
class TClonesArray ;
class AliCFManager;
class AliAODRecoDecay;
class AliAODRecoDecayHF2Prong;
class AliAODMCParticle;
class THnSparse;
class AliRDHFCuts;
class AliCFVertexingHF2Prong;
class AliCFVertexingHF3Prong;

class AliCFTaskVertexingHF: public AliAnalysisTaskSE {
public:
	
	enum {
	        kStepGeneratedLimAcc = 0,
		kStepGenerated       = 1,
		kStepAcceptance      = 2,
		kStepVertex          = 3,
		kStepRefit           = 4,
		kStepReconstructed   = 5,
		kStepRecoAcceptance  = 6,
		kStepRecoITSClusters = 7,
		kStepRecoPPR         = 8,
		kStepRecoPID         = 9
	};
	
	AliCFTaskVertexingHF();
	AliCFTaskVertexingHF(const Char_t* name, AliRDHFCuts* cuts);
	AliCFTaskVertexingHF& operator= (const AliCFTaskVertexingHF& c);
	AliCFTaskVertexingHF(const AliCFTaskVertexingHF& c);
 	virtual ~AliCFTaskVertexingHF();
	
	// ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects
	void     UserCreateOutputObjects();
	void     UserExec(Option_t *option);
	void     Init();
	void     LocalInit() {Init();}
	void     Terminate(Option_t *);
	
	// UNFOLDING
	void     SetCorrelationMatrix(THnSparse* h) {fCorrelation=h;}
	void     SetAcceptanceUnf(Bool_t AcceptanceUnf) {fAcceptanceUnf = AcceptanceUnf;}
	Bool_t   GetAcceptanceUnf() const {return fAcceptanceUnf;}
	
	
	// CORRECTION FRAMEWORK RELATED FUNCTIONS
	void           SetCFManager(AliCFManager* io) {fCFManager = io;}   // global correction manager
	AliCFManager * GetCFManager()                 {return fCFManager;} // get corr manager
	
	// Setters (and getters) for the config macro
	void    SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;}
	Bool_t  GetFillFromGenerated() const {return fFillFromGenerated;}
	void    SetDecayChannel (Int_t decayChannel) {fDecayChannel = decayChannel;}
	Int_t   GetDecayChannel () {return fDecayChannel;}
	void     SetUseWeight(Bool_t useWeight){fUseWeight=useWeight;}
	Bool_t   GetUseWeight() const {return fUseWeight;}
	Double_t GetWeight(Float_t pt);
	Double_t dNdptFit(Float_t pt, Double_t* par);
 	
	void   SetDselection(UShort_t originDselection) {fOriginDselection=originDselection;}
	UShort_t GetDselection (){return fOriginDselection;}
	void SetSign(Char_t isSign) {fSign = isSign;}
	Char_t GetSign() {return fSign;}
	 
	void SetCentralitySelection(Bool_t centSelec = kTRUE) {fCentralitySelection = centSelec;}        Bool_t GetCentralitySelection() {return fCentralitySelection;} 

protected:
	AliCFManager   *fCFManager;   //  pointer to the CF manager
	TH1I *fHistEventsProcessed;   //! simple histo for monitoring the number of events processed
	THnSparse* fCorrelation;      //  response matrix for unfolding
	Int_t fCountMC;               //  MC particle found
	Int_t fCountAcc;              //  MC particle found that satisfy acceptance cuts
	Int_t fCountVertex;       //  Reco particle found that satisfy vertex constrained
	Int_t fCountRefit;        //  Reco particle found that satisfy kTPCrefit and kITSrefit
	Int_t fCountReco;             //  Reco particle found that satisfy cuts
	Int_t fCountRecoAcc;          //  Reco particle found that satisfy cuts in requested acceptance
	Int_t fCountRecoITSClusters;  //  Reco particle found that satisfy cuts in n. of ITS clusters
	Int_t fCountRecoPPR;          //  Reco particle found that satisfy cuts in PPR
	Int_t fCountRecoPID;          //Reco PID step 
	Int_t fEvents;                //  n. of events
	Int_t fDecayChannel;          // decay channel to configure the task
	Bool_t fFillFromGenerated;    //  flag to indicate whether data container should be filled with generated values also for reconstructed particles
	UShort_t fOriginDselection;      // flag to select D0 origins. 0 Only from charm 1 only from beauty 2 both from charm and beauty
	Bool_t fAcceptanceUnf;        //  flag for unfolding before or after cuts.
	AliRDHFCuts* fCuts;            // cuts
	Bool_t fUseWeight;             //flag to decide whether to use weights != 1 when filling the container or not
	Double_t fWeight;              //weight used to fill the container
	Int_t fNvar;                   // number of variables for the container
	TString fPartName;    // D meson name
	TString fDauNames;    // daughter in fin state
	Char_t fSign;                 // flag to decide wheter to keep D0 only (0), D0bar only (1), or both D0 and D0bar (2)
        Bool_t fCentralitySelection;  //flag to switch off the centrality selection


	ClassDef(AliCFTaskVertexingHF,5); // class for HF corrections as a function of many variables
};

#endif
Commit	Line	Data
379592af	1	#ifndef ALICFTASKVERTEXINGHF_H
	2	#define ALICFTASKVERTEXINGHF_H
	3	/**************************************************************************
	4	* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
	5	* *
	6	* Author: The ALICE Off-line Project. *
	7	* Contributors are mentioned in the code where appropriate. *
	8	* *
	9	* Permission to use, copy, modify and distribute this software and its *
	10	* documentation strictly for non-commercial purposes is hereby granted *
	11	* without fee, provided that the above copyright notice appears in all *
	12	* copies and that both the copyright notice and this permission notice *
	13	* appear in the supporting documentation. The authors make no claims *
	14	* about the suitability of this software for any purpose. It is *
	15	* provided "as is" without express or implied warranty. *
	16	**************************************************************************/
	17
	18	//-----------------------------------------------------------------------
	19	// Class for HF corrections as a function of many variables and step
	20	// Author : C. Zampolli, CERN
	21	// D. Caffarri, Univ & INFN Padova caffarri@pd.infn.it
	22	// Base class for HF Unfolding - agrelli@uu.nl
	23	//-----------------------------------------------------------------------
	24
	25
	26	#include "AliAnalysisTaskSE.h"
	27	#include "AliCFVertexingHF2Prong.h"
	28	#include "AliCFVertexingHF.h"
	29
	30	class TH1I;
	31	class TParticle ;
	32	class TFile ;
	33	class TClonesArray ;
	34	class AliCFManager;
	35	class AliAODRecoDecay;
	36	class AliAODRecoDecayHF2Prong;
	37	class AliAODMCParticle;
	38	class THnSparse;
f2703bd2	39	class AliRDHFCuts;
379592af	40	class AliCFVertexingHF2Prong;
6e2e4f50	41	class AliCFVertexingHF3Prong;
379592af	42
	43	class AliCFTaskVertexingHF: public AliAnalysisTaskSE {
	44	public:
	45
	46	enum {
	47	kStepGeneratedLimAcc = 0,
	48	kStepGenerated = 1,
	49	kStepAcceptance = 2,
	50	kStepVertex = 3,
	51	kStepRefit = 4,
	52	kStepReconstructed = 5,
	53	kStepRecoAcceptance = 6,
	54	kStepRecoITSClusters = 7,
	55	kStepRecoPPR = 8,
	56	kStepRecoPID = 9
	57	};
	58
	59	AliCFTaskVertexingHF();
f2703bd2	60	AliCFTaskVertexingHF(const Char_t* name, AliRDHFCuts* cuts);
379592af	61	AliCFTaskVertexingHF& operator= (const AliCFTaskVertexingHF& c);
	62	AliCFTaskVertexingHF(const AliCFTaskVertexingHF& c);
	63	virtual ~AliCFTaskVertexingHF();
	64
	65	// ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects
	66	void UserCreateOutputObjects();
	67	void UserExec(Option_t *option);
f2703bd2	68	void Init();
f2703bd2	69	void LocalInit() {Init();}
379592af	70	void Terminate(Option_t *);
	71
	72	// UNFOLDING
	73	void SetCorrelationMatrix(THnSparse* h) {fCorrelation=h;}
	74	void SetAcceptanceUnf(Bool_t AcceptanceUnf) {fAcceptanceUnf = AcceptanceUnf;}
	75	Bool_t GetAcceptanceUnf() const {return fAcceptanceUnf;}
	76
	77
	78	// CORRECTION FRAMEWORK RELATED FUNCTIONS
	79	void SetCFManager(AliCFManager* io) {fCFManager = io;} // global correction manager
	80	AliCFManager * GetCFManager() {return fCFManager;} // get corr manager
	81
f2703bd2	82	// Setters (and getters) for the config macro
	83	void SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;}
	84	Bool_t GetFillFromGenerated() const {return fFillFromGenerated;}
	85	void SetDecayChannel (Int_t decayChannel) {fDecayChannel = decayChannel;}
3ee5eb83	86	Int_t GetDecayChannel () {return fDecayChannel;}
f2703bd2	87	void SetUseWeight(Bool_t useWeight){fUseWeight=useWeight;}
	88	Bool_t GetUseWeight() const {return fUseWeight;}
	89	Double_t GetWeight(Float_t pt);
	90	Double_t dNdptFit(Float_t pt, Double_t* par);
	91
6e2e4f50	92	void SetDselection(UShort_t originDselection) {fOriginDselection=originDselection;}
3ee5eb83	93	UShort_t GetDselection (){return fOriginDselection;}
	94	void SetSign(Char_t isSign) {fSign = isSign;}
	95	Char_t GetSign() {return fSign;}
379592af	96
b7af2639	97	void SetCentralitySelection(Bool_t centSelec = kTRUE) {fCentralitySelection = centSelec;} Bool_t GetCentralitySelection() {return fCentralitySelection;}
b7af2639	98
379592af	99	protected:
379592af	100	AliCFManager *fCFManager; // pointer to the CF manager
	101	TH1I *fHistEventsProcessed; //! simple histo for monitoring the number of events processed
	102	THnSparse* fCorrelation; // response matrix for unfolding
	103	Int_t fCountMC; // MC particle found
	104	Int_t fCountAcc; // MC particle found that satisfy acceptance cuts
	105	Int_t fCountVertex; // Reco particle found that satisfy vertex constrained
	106	Int_t fCountRefit; // Reco particle found that satisfy kTPCrefit and kITSrefit
	107	Int_t fCountReco; // Reco particle found that satisfy cuts
	108	Int_t fCountRecoAcc; // Reco particle found that satisfy cuts in requested acceptance
	109	Int_t fCountRecoITSClusters; // Reco particle found that satisfy cuts in n. of ITS clusters
	110	Int_t fCountRecoPPR; // Reco particle found that satisfy cuts in PPR
	111	Int_t fCountRecoPID; //Reco PID step
	112	Int_t fEvents; // n. of events
f2703bd2	113	Int_t fDecayChannel; // decay channel to configure the task
379592af	114	Bool_t fFillFromGenerated; // flag to indicate whether data container should be filled with generated values also for reconstructed particles
	115	UShort_t fOriginDselection; // flag to select D0 origins. 0 Only from charm 1 only from beauty 2 both from charm and beauty
	116	Bool_t fAcceptanceUnf; // flag for unfolding before or after cuts.
f2703bd2	117	AliRDHFCuts* fCuts; // cuts
	118	Bool_t fUseWeight; //flag to decide whether to use weights != 1 when filling the container or not
	119	Double_t fWeight; //weight used to fill the container
	120	Int_t fNvar; // number of variables for the container
31da6b05	121	TString fPartName; // D meson name
31da6b05	122	TString fDauNames; // daughter in fin state
3ee5eb83	123	Char_t fSign; // flag to decide wheter to keep D0 only (0), D0bar only (1), or both D0 and D0bar (2)
b7af2639	124	Bool_t fCentralitySelection; //flag to switch off the centrality selection
b7af2639	125
3ee5eb83	126
b7af2639	127	ClassDef(AliCFTaskVertexingHF,5); // class for HF corrections as a function of many variables
379592af	128	};
	129
	130	#endif