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

#ifndef ALICFHEAVYFLAVOURTASKMULTIVARMULTISTEP_H
#define ALICFHEAVYFLAVOURTASKMULTIVARMULTISTEP_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
// Base class for HF Unfolding - agrelli@uu.nl
//-----------------------------------------------------------------------


#include "AliAnalysisTaskSE.h"

class TH1I;
class TParticle ;
class TFile ;
class TClonesArray ;
class AliCFManager;
class AliAODRecoDecay;
class AliAODRecoDecayHF2Prong;
class AliAODMCParticle;
class THnSparse;

class AliCFHeavyFlavourTaskMultiVarMultiStep : public AliAnalysisTaskSE {
  public:

  enum {
    kStepGenerated       = 0,
    kStepAcceptance      = 1,
    kStepVertex          = 2,
    kStepRefit           = 3,
    kStepReconstructed   = 4,
    kStepRecoAcceptance  = 5,
    kStepRecoITSClusters = 6,
    kStepRecoPPR         = 7
  };

  AliCFHeavyFlavourTaskMultiVarMultiStep();
  AliCFHeavyFlavourTaskMultiVarMultiStep(const Char_t* name);
  AliCFHeavyFlavourTaskMultiVarMultiStep& operator= (const AliCFHeavyFlavourTaskMultiVarMultiStep& c);
  AliCFHeavyFlavourTaskMultiVarMultiStep(const AliCFHeavyFlavourTaskMultiVarMultiStep& c);
  virtual ~AliCFHeavyFlavourTaskMultiVarMultiStep();

  // ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects
  void     UserCreateOutputObjects();
  void     UserExec(Option_t *option);
  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

  void     SetPDG(Int_t code) {fPDG = code; }     // defines the PDG code of searched HF
  Double_t CosThetaStar(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const;  // returns cos(ThetaStar) of the D0 decay
  Double_t CT(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const;            // returns cT of the D0 decay
  void     SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;}
  Bool_t   GetFillFromGenerated() const {return fFillFromGenerated;}
  Bool_t   GetGeneratedValuesFromMCParticle(AliAODMCParticle* mcPart, TClonesArray* mcArray, Double_t* vectorMC) const;
  void     SetMinITSClusters(Int_t minITSClusters) {fMinITSClusters = minITSClusters;}
  Int_t    GetMinITSClusters() const {return fMinITSClusters;}
  Int_t    CheckOrigin(AliAODMCParticle* mcPart, TClonesArray* mcArray) const;

 protected:
  Int_t           fPDG;         //  PDG code of searched V0's
  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 fEvents;                //  n. of events
  Bool_t fFillFromGenerated;    //  flag to indicate whether data container should be filled 
                                //  with generated values also for reconstructed particles
  Int_t fMinITSClusters;        //  min n. of ITS clusters for RecoDecay
  Bool_t fAcceptanceUnf;        //  flag for unfolding before or after cuts.
  
  ClassDef(AliCFHeavyFlavourTaskMultiVarMultiStep,3); // class for HF corrections as a function of many variables
};

#endif
Commit	Line	Data
f5ec6c30	1	#ifndef ALICFHEAVYFLAVOURTASKMULTIVARMULTISTEP_H
	2	#define ALICFHEAVYFLAVOURTASKMULTIVARMULTISTEP_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
a5d92cb4	21	// Base class for HF Unfolding - agrelli@uu.nl
f5ec6c30	22	//-----------------------------------------------------------------------
	23
	24
	25	#include "AliAnalysisTaskSE.h"
	26
	27	class TH1I;
	28	class TParticle ;
	29	class TFile ;
	30	class TClonesArray ;
	31	class AliCFManager;
	32	class AliAODRecoDecay;
	33	class AliAODRecoDecayHF2Prong;
	34	class AliAODMCParticle;
a5d92cb4	35	class THnSparse;
f5ec6c30	36
	37	class AliCFHeavyFlavourTaskMultiVarMultiStep : public AliAnalysisTaskSE {
	38	public:
	39
	40	enum {
	41	kStepGenerated = 0,
	42	kStepAcceptance = 1,
ae3e319c	43	kStepVertex = 2,
	44	kStepRefit = 3,
	45	kStepReconstructed = 4,
	46	kStepRecoAcceptance = 5,
	47	kStepRecoITSClusters = 6,
	48	kStepRecoPPR = 7
f5ec6c30	49	};
	50
	51	AliCFHeavyFlavourTaskMultiVarMultiStep();
	52	AliCFHeavyFlavourTaskMultiVarMultiStep(const Char_t* name);
	53	AliCFHeavyFlavourTaskMultiVarMultiStep& operator= (const AliCFHeavyFlavourTaskMultiVarMultiStep& c);
	54	AliCFHeavyFlavourTaskMultiVarMultiStep(const AliCFHeavyFlavourTaskMultiVarMultiStep& c);
	55	virtual ~AliCFHeavyFlavourTaskMultiVarMultiStep();
	56
	57	// ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects
	58	void UserCreateOutputObjects();
	59	void UserExec(Option_t *option);
	60	void Terminate(Option_t *);
a5d92cb4	61
	62	// UNFOLDING
	63	void SetCorrelationMatrix(THnSparse* h) {fCorrelation=h;}
	64	void SetAcceptanceUnf(Bool_t AcceptanceUnf) {fAcceptanceUnf = AcceptanceUnf;}
	65	Bool_t GetAcceptanceUnf() const {return fAcceptanceUnf;}
	66
f5ec6c30	67
	68	// CORRECTION FRAMEWORK RELATED FUNCTIONS
	69	void SetCFManager(AliCFManager* io) {fCFManager = io;} // global correction manager
	70	AliCFManager * GetCFManager() {return fCFManager;} // get corr manager
	71
	72	void SetPDG(Int_t code) {fPDG = code; } // defines the PDG code of searched HF
	73	Double_t CosThetaStar(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const; // returns cos(ThetaStar) of the D0 decay
	74	Double_t CT(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const; // returns cT of the D0 decay
	75	void SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;}
	76	Bool_t GetFillFromGenerated() const {return fFillFromGenerated;}
	77	Bool_t GetGeneratedValuesFromMCParticle(AliAODMCParticle* mcPart, TClonesArray* mcArray, Double_t* vectorMC) const;
	78	void SetMinITSClusters(Int_t minITSClusters) {fMinITSClusters = minITSClusters;}
	79	Int_t GetMinITSClusters() const {return fMinITSClusters;}
e8158e56	80	Int_t CheckOrigin(AliAODMCParticle* mcPart, TClonesArray* mcArray) const;
f5ec6c30	81
	82	protected:
	83	Int_t fPDG; // PDG code of searched V0's
4df1e470	84	AliCFManager *fCFManager; // pointer to the CF manager
f5ec6c30	85	TH1I *fHistEventsProcessed; //! simple histo for monitoring the number of events processed
a5d92cb4	86	THnSparse* fCorrelation; // response matrix for unfolding
f5ec6c30	87	Int_t fCountMC; // MC particle found
f5ec6c30	88	Int_t fCountAcc; // MC particle found that satisfy acceptance cuts
ae3e319c	89	Int_t fCountVertex; // Reco particle found that satisfy vertex constrained
ae3e319c	90	Int_t fCountRefit; // Reco particle found that satisfy kTPCrefit and kITSrefit
f5ec6c30	91	Int_t fCountReco; // Reco particle found that satisfy cuts
	92	Int_t fCountRecoAcc; // Reco particle found that satisfy cuts in requested acceptance
	93	Int_t fCountRecoITSClusters; // Reco particle found that satisfy cuts in n. of ITS clusters
	94	Int_t fCountRecoPPR; // Reco particle found that satisfy cuts in PPR
	95	Int_t fEvents; // n. of events
	96	Bool_t fFillFromGenerated; // flag to indicate whether data container should be filled
	97	// with generated values also for reconstructed particles
	98	Int_t fMinITSClusters; // min n. of ITS clusters for RecoDecay
a5d92cb4	99	Bool_t fAcceptanceUnf; // flag for unfolding before or after cuts.
f5ec6c30	100
e5450e63	101	ClassDef(AliCFHeavyFlavourTaskMultiVarMultiStep,3); // class for HF corrections as a function of many variables
f5ec6c30	102	};
	103
	104	#endif