[u/mrichter/AliRoot.git] / PWGCF / EBYE / BalanceFunctions / AliAnalysisTaskBFPsi.h

#ifndef ALIANALYSISTASKBFPSI_H\r
#define ALIANALYSISTASKBFPSI_H\r
\r
// Analysis task for the BF vs Psi code\r
// Authors: Panos Cristakoglou@cern.ch\r
\r
class TList;\r
class TH1F;\r
class TH2F;\r
class TF1;\r
\r
class AliBalancePsi;\r
class AliESDtrackCuts;\r
class AliEventPoolManager;\r
\r
\r
#include "AliAnalysisTaskSE.h"\r
#include "AliBalancePsi.h"\r
\r
#include "AliPID.h"  \r
#include "AliPIDResponse.h"\r
#include "AliPIDCombined.h"\r
 \r
\r
class AliAnalysisTaskBFPsi : public AliAnalysisTaskSE {\r
 public:\r
  AliAnalysisTaskBFPsi(const char *name = "AliAnalysisTaskBFPsi");\r
  virtual ~AliAnalysisTaskBFPsi(); \r
  \r
  \r
  virtual void   UserCreateOutputObjects();\r
  virtual void   UserExec(Option_t *option);\r
  virtual void   FinishTaskOutput();\r
  virtual void   Terminate(Option_t *);\r
\r
  void SetAnalysisObject(AliBalancePsi *const analysis) {\r
    fBalance         = analysis;\r
    }\r
  void SetShufflingObject(AliBalancePsi *const analysisShuffled) {\r
    fRunShuffling = kTRUE;\r
    fShuffledBalance = analysisShuffled;\r
  }\r
  void SetMixingObject(AliBalancePsi *const analysisMixed) {\r
    fRunMixing = kTRUE;\r
    fMixedBalance = analysisMixed;\r
  }\r
  void SetMixingTracks(Int_t tracks) { fMixingTracks = tracks; }\r
  void SetAnalysisCutObject(AliESDtrackCuts *const trackCuts) {\r
    fESDtrackCuts = trackCuts;}\r
  void SetVertexDiamond(Double_t vx, Double_t vy, Double_t vz) {\r
    fVxMax = vx;\r
    fVyMax = vy;\r
    fVzMax = vz;\r
  }\r
\r
  //==============AOD analysis==============//\r
  void SetAODtrackCutBit(Int_t bit){\r
    nAODtrackCutBit = bit;\r
  }\r
\r
  void SetKinematicsCutsAOD(Double_t ptmin, Double_t ptmax, Double_t etamin, Double_t etamax){\r
    fPtMin  = ptmin;  fPtMax  = ptmax;\r
    fEtaMin = etamin; fEtaMax = etamax;\r
  }\r
\r
  void SetExtraDCACutsAOD(Double_t DCAxy, Double_t DCAz){\r
    fDCAxyCut  = DCAxy;\r
    fDCAzCut = DCAz;\r
  }\r
\r
   void SetExtraTPCCutsAOD(Double_t maxTPCchi2, Int_t minNClustersTPC){\r
    fTPCchi2Cut      = maxTPCchi2;\r
    fNClustersTPCCut = minNClustersTPC;\r
  }\r
\r
  //==============MC analysis==============//\r
  void SetKinematicsCutsMC(Double_t ptmin, Double_t ptmax,\r
                           Double_t etamin, Double_t etamax){\r
    fPtMin  = ptmin; fPtMax  = ptmax;\r
    fEtaMin = etamin; fEtaMax = etamax;\r
  }\r
  void UseFlowAfterBurner(TF1 *gDifferentialV2) {\r
    fDifferentialV2 = gDifferentialV2;\r
    fUseFlowAfterBurner = kTRUE;\r
  }\r
  void ExcludeResonancesInMC() {fExcludeResonancesInMC = kTRUE;}\r
\r
  void SetPDGCode(Int_t gPdgCode) {\r
    fUseMCPdgCode = kTRUE;\r
    fPDGCodeToBeAnalyzed = gPdgCode;\r
  }\r
\r
  //Centrality\r
  void SetCentralityEstimator(const char* centralityEstimator) {fCentralityEstimator = centralityEstimator;}\r
  const char* GetCentralityEstimator(void)  const              {return fCentralityEstimator;}\r
  void SetCentralityPercentileRange(Double_t min, Double_t max) { \r
    fUseCentrality = kTRUE;\r
    fCentralityPercentileMin=min;\r
    fCentralityPercentileMax=max;\r
  }\r
  void SetImpactParameterRange(Double_t min, Double_t max) { \r
    fUseCentrality = kTRUE;\r
    fImpactParameterMin=min;\r
    fImpactParameterMax=max;\r
  }\r
\r
  //multiplicity\r
  void SetMultiplicityRange(Int_t min, Int_t max) {\r
    fUseMultiplicity = kTRUE;\r
    fNumberOfAcceptedTracksMin = min;\r
    fNumberOfAcceptedTracksMax = max;}\r
  \r
  void UseOfflineTrigger() {fUseOfflineTrigger = kTRUE;}\r
  \r
  //Acceptance filter\r
  void SetAcceptanceParameterization(TF1 *parameterization) {\r
    fAcceptanceParameterization = parameterization;}\r
\r
  //pid\r
  enum kDetectorUsedForPID { kTPCpid, kTOFpid, kTPCTOF }; // default TPC & TOF pid (via GetTPCpid & GetTOFpid)  \r
  enum kParticleOfInterest { kMuon, kElectron, kPion, kKaon, kProton };  \r
\r
  void SetUseBayesianPID(Double_t gMinProbabilityValue) {\r
    fUsePID = kTRUE; fUsePIDnSigma = kFALSE; fUsePIDPropabilities = kTRUE;\r
    fMinAcceptedPIDProbability = gMinProbabilityValue; }\r
\r
  void SetUseNSigmaPID(Double_t gMaxNSigma) {\r
    fUsePID = kTRUE; fUsePIDPropabilities = kFALSE; fUsePIDnSigma = kTRUE;\r
    fPIDNSigma = gMaxNSigma; }\r
\r
  void SetParticleOfInterest(kParticleOfInterest poi) {\r
    fParticleOfInterest = poi;}\r
  void SetDetectorUsedForPID(kDetectorUsedForPID detConfig) {\r
    fPidDetectorConfig = detConfig;}\r
\r
 private:\r
  Double_t    IsEventAccepted(AliVEvent* event);\r
  Double_t    GetEventPlane(AliVEvent* event);\r
  TObjArray* GetAcceptedTracks(AliVEvent* event, Double_t fCentrality, Double_t gReactionPlane);\r
  TObjArray* GetShuffledTracks(TObjArray* tracks);\r
\r
  AliBalancePsi *fBalance; //BF object\r
  Bool_t fRunShuffling;//run shuffling or not\r
  AliBalancePsi *fShuffledBalance; //BF object (shuffled)\r
  Bool_t fRunMixing;//run mixing or not\r
  Int_t  fMixingTracks;\r
  AliBalancePsi *fMixedBalance; //TriggeredBF object (mixed)\r
  AliEventPoolManager*     fPoolMgr;         //! event pool manager\r
\r
  TList *fList; //fList object\r
  TList *fListBF; //fList object\r
  TList *fListBFS; //fList object\r
  TList *fListBFM; //fList object\r
  TList *fHistListPIDQA;  //! list of histograms\r
\r
  TH2F *fHistEventStats; //event stats\r
  TH2F *fHistCentStats; //centrality stats\r
  TH1F *fHistTriggerStats; //trigger stats\r
  TH1F *fHistTrackStats; //Track filter bit stats\r
  TH1F *fHistVx; //x coordinate of the primary vertex\r
  TH1F *fHistVy; //y coordinate of the primary vertex\r
  TH2F *fHistVz; //z coordinate of the primary vertex\r
\r
  TH2F *fHistEventPlane; //event plane distribution\r
\r
  TH2F *fHistClus;//number of clusters (QA histogram)\r
  TH2F *fHistDCA;//DCA  (QA histogram)\r
  TH2F *fHistChi2;//track chi2 (QA histogram)\r
  TH2F *fHistPt;//transverse momentum (QA histogram)\r
  TH2F *fHistEta;//pseudorapidity (QA histogram)\r
  TH2F *fHistRapidity;//rapidity (QA histogram)\r
  TH2F *fHistPhi;//phi (QA histogram)\r
  TH2F *fHistPhiBefore;//phi before v2 afterburner (QA histogram)\r
  TH2F *fHistPhiAfter;//phi after v2 afterburner (QA histogram)\r
  TH2F *fHistPhiPos;//phi for positive particles (QA histogram)\r
  TH2F *fHistPhiNeg;//phi for negative particles (QA histogram)\r
  TH2F *fHistV0M;//V0 multiplicities (QA histogram)\r
  TH2F *fHistRefTracks;//reference track multiplicities (QA histogram)\r
\r
  //============PID============//\r
  TH2D *fHistdEdxVsPTPCbeforePID;//TPC dEdx vs momentum before PID cuts (QA histogram)\r
  TH2D *fHistBetavsPTOFbeforePID;//beta vs momentum before PID cuts (QA histogram)\r
  TH2D *fHistProbTPCvsPtbeforePID; //TPC probability vs pT before PID cuts (QA histogram)\r
  TH2D *fHistProbTOFvsPtbeforePID;//TOF probability vs pT before PID cuts (QA histogram)\r
  TH2D *fHistProbTPCTOFvsPtbeforePID;//TOF/TPC probability vs pT before PID cuts (QA histogram)\r
  TH2D *fHistNSigmaTPCvsPtbeforePID;//TPC nsigma vs pT before PID cuts (QA histogram)\r
  TH2D *fHistNSigmaTOFvsPtbeforePID;//TOF nsigma vs pT before PID cuts (QA histogram)\r
  TH2D *fHistdEdxVsPTPCafterPID;//TPC dEdx vs momentum after PID cuts (QA histogram)\r
  TH2D *fHistBetavsPTOFafterPID;//beta vs momentum after PID cuts (QA histogram)\r
  TH2D *fHistProbTPCvsPtafterPID; //TPC probability vs pT after PID cuts (QA histogram)\r
  TH2D *fHistProbTOFvsPtafterPID;//TOF probability vs pT after PID cuts (QA histogram)\r
  TH2D *fHistProbTPCTOFvsPtafterPID;//TOF/TPC probability vs pT after PID cuts (QA histogram)\r
  TH2D *fHistNSigmaTPCvsPtafterPID;//TPC nsigma vs pT after PID cuts (QA histogram)\r
  TH2D *fHistNSigmaTOFvsPtafterPID;//TOF nsigma vs pT after PID cuts (QA histogram)\r
\r
  AliPIDResponse *fPIDResponse;     //! PID response object\r
  AliPIDCombined       *fPIDCombined;     //! combined PID object\r
  \r
  kParticleOfInterest  fParticleOfInterest;//analyzed particle\r
  kDetectorUsedForPID   fPidDetectorConfig;//used detector for PID\r
\r
  Bool_t fUsePID; //flag to use PID \r
  Bool_t fUsePIDnSigma;//flag to use nsigma method for PID\r
  Bool_t fUsePIDPropabilities;//flag to use probability method for PID\r
  Double_t fPIDNSigma;//nsigma cut for PID\r
  Double_t fMinAcceptedPIDProbability;//probability cut for PID\r
  //============PID============//\r
\r
  AliESDtrackCuts *fESDtrackCuts; //ESD track cuts\r
\r
  TString fCentralityEstimator;      //"V0M","TRK","TKL","ZDC","FMD"\r
  Bool_t fUseCentrality;//use the centrality (PbPb) or not (pp)\r
  Double_t fCentralityPercentileMin;//centrality percentile min\r
  Double_t fCentralityPercentileMax;//centrality percentile max\r
  Double_t fImpactParameterMin;//impact parameter min (used for MC)\r
  Double_t fImpactParameterMax;//impact parameter max (used for MC)\r
\r
  Bool_t fUseMultiplicity;//use the multiplicity cuts\r
  Int_t fNumberOfAcceptedTracksMin;//min. number of number of accepted tracks (used for the multiplicity dependence study - pp)\r
  Int_t fNumberOfAcceptedTracksMax;//max. number of number of accepted tracks (used for the multiplicity dependence study - pp)\r
  TH2F *fHistNumberOfAcceptedTracks;//hisot to store the number of accepted tracks\r
\r
  Bool_t fUseOfflineTrigger;//Usage of the offline trigger selection\r
\r
  Double_t fVxMax;//vxmax\r
  Double_t fVyMax;//vymax\r
  Double_t fVzMax;//vzmax\r
\r
  Int_t nAODtrackCutBit;//track cut bit from track selection (only used for AODs)\r
\r
  Double_t fPtMin;//only used for AODs\r
  Double_t fPtMax;//only used for AODs\r
  Double_t fEtaMin;//only used for AODs\r
  Double_t fEtaMax;//only used for AODs\r
\r
  Double_t fDCAxyCut;//only used for AODs\r
  Double_t fDCAzCut;//only used for AODs\r
\r
  Double_t fTPCchi2Cut;//only used for AODs\r
  Int_t fNClustersTPCCut;//only used for AODs\r
\r
  TF1 *fAcceptanceParameterization;//acceptance filter used for MC\r
\r
  TF1 *fDifferentialV2;//pt-differential v2 (from real data)\r
  Bool_t fUseFlowAfterBurner;//Usage of a flow after burner\r
\r
  Bool_t fExcludeResonancesInMC;//flag to exclude the resonances' decay products from the MC analysis\r
  Bool_t fUseMCPdgCode; //Boolean to analyze a set of particles in MC\r
  Int_t fPDGCodeToBeAnalyzed; //Analyze a set of particles in MC\r
\r
  \r
\r
  AliAnalysisTaskBFPsi(const AliAnalysisTaskBFPsi&); // not implemented\r
  AliAnalysisTaskBFPsi& operator=(const AliAnalysisTaskBFPsi&); // not implemented\r
  \r
  ClassDef(AliAnalysisTaskBFPsi, 5); // example of analysis\r
};\r
\r
\r
\r
#endif\r
Commit	Line	Data
9fd4b54e	1	#ifndef ALIANALYSISTASKBFPSI_H\r
9fd4b54e	2	#define ALIANALYSISTASKBFPSI_H\r
0879e280	3	\r
	4	// Analysis task for the BF vs Psi code\r
	5	// Authors: Panos Cristakoglou@cern.ch\r
	6	\r
	7	class TList;\r
	8	class TH1F;\r
	9	class TH2F;\r
	10	class TF1;\r
	11	\r
	12	class AliBalancePsi;\r
	13	class AliESDtrackCuts;\r
f06d59b3	14	class AliEventPoolManager;\r
f06d59b3	15	\r
0879e280	16	\r
	17	#include "AliAnalysisTaskSE.h"\r
	18	#include "AliBalancePsi.h"\r
	19	\r
	20	#include "AliPID.h" \r
	21	#include "AliPIDResponse.h"\r
	22	#include "AliPIDCombined.h"\r
	23	\r
	24	\r
	25	class AliAnalysisTaskBFPsi : public AliAnalysisTaskSE {\r
	26	public:\r
	27	AliAnalysisTaskBFPsi(const char *name = "AliAnalysisTaskBFPsi");\r
	28	virtual ~AliAnalysisTaskBFPsi(); \r
	29	\r
	30	\r
	31	virtual void UserCreateOutputObjects();\r
	32	virtual void UserExec(Option_t *option);\r
	33	virtual void FinishTaskOutput();\r
	34	virtual void Terminate(Option_t *);\r
	35	\r
	36	void SetAnalysisObject(AliBalancePsi *const analysis) {\r
	37	fBalance = analysis;\r
	38	}\r
	39	void SetShufflingObject(AliBalancePsi *const analysisShuffled) {\r
	40	fRunShuffling = kTRUE;\r
	41	fShuffledBalance = analysisShuffled;\r
	42	}\r
f06d59b3	43	void SetMixingObject(AliBalancePsi *const analysisMixed) {\r
	44	fRunMixing = kTRUE;\r
	45	fMixedBalance = analysisMixed;\r
	46	}\r
	47	void SetMixingTracks(Int_t tracks) { fMixingTracks = tracks; }\r
0879e280	48	void SetAnalysisCutObject(AliESDtrackCuts *const trackCuts) {\r
	49	fESDtrackCuts = trackCuts;}\r
	50	void SetVertexDiamond(Double_t vx, Double_t vy, Double_t vz) {\r
	51	fVxMax = vx;\r
	52	fVyMax = vy;\r
	53	fVzMax = vz;\r
	54	}\r
	55	\r
	56	//==============AOD analysis==============//\r
	57	void SetAODtrackCutBit(Int_t bit){\r
	58	nAODtrackCutBit = bit;\r
	59	}\r
	60	\r
	61	void SetKinematicsCutsAOD(Double_t ptmin, Double_t ptmax, Double_t etamin, Double_t etamax){\r
9fd4b54e	62	fPtMin = ptmin; fPtMax = ptmax;\r
9fd4b54e	63	fEtaMin = etamin; fEtaMax = etamax;\r
0879e280	64	}\r
	65	\r
	66	void SetExtraDCACutsAOD(Double_t DCAxy, Double_t DCAz){\r
	67	fDCAxyCut = DCAxy;\r
	68	fDCAzCut = DCAz;\r
	69	}\r
	70	\r
	71	void SetExtraTPCCutsAOD(Double_t maxTPCchi2, Int_t minNClustersTPC){\r
	72	fTPCchi2Cut = maxTPCchi2;\r
	73	fNClustersTPCCut = minNClustersTPC;\r
	74	}\r
	75	\r
	76	//==============MC analysis==============//\r
	77	void SetKinematicsCutsMC(Double_t ptmin, Double_t ptmax,\r
	78	Double_t etamin, Double_t etamax){\r
	79	fPtMin = ptmin; fPtMax = ptmax;\r
	80	fEtaMin = etamin; fEtaMax = etamax;\r
	81	}\r
	82	void UseFlowAfterBurner(TF1 *gDifferentialV2) {\r
	83	fDifferentialV2 = gDifferentialV2;\r
	84	fUseFlowAfterBurner = kTRUE;\r
	85	}\r
	86	void ExcludeResonancesInMC() {fExcludeResonancesInMC = kTRUE;}\r
	87	\r
	88	void SetPDGCode(Int_t gPdgCode) {\r
	89	fUseMCPdgCode = kTRUE;\r
	90	fPDGCodeToBeAnalyzed = gPdgCode;\r
	91	}\r
	92	\r
	93	//Centrality\r
	94	void SetCentralityEstimator(const char* centralityEstimator) {fCentralityEstimator = centralityEstimator;}\r
9fd4b54e	95	const char* GetCentralityEstimator(void) const {return fCentralityEstimator;}\r
0879e280	96	void SetCentralityPercentileRange(Double_t min, Double_t max) { \r
	97	fUseCentrality = kTRUE;\r
	98	fCentralityPercentileMin=min;\r
	99	fCentralityPercentileMax=max;\r
	100	}\r
	101	void SetImpactParameterRange(Double_t min, Double_t max) { \r
	102	fUseCentrality = kTRUE;\r
	103	fImpactParameterMin=min;\r
	104	fImpactParameterMax=max;\r
	105	}\r
	106	\r
	107	//multiplicity\r
	108	void SetMultiplicityRange(Int_t min, Int_t max) {\r
	109	fUseMultiplicity = kTRUE;\r
	110	fNumberOfAcceptedTracksMin = min;\r
	111	fNumberOfAcceptedTracksMax = max;}\r
	112	\r
	113	void UseOfflineTrigger() {fUseOfflineTrigger = kTRUE;}\r
	114	\r
	115	//Acceptance filter\r
	116	void SetAcceptanceParameterization(TF1 *parameterization) {\r
	117	fAcceptanceParameterization = parameterization;}\r
	118	\r
	119	//pid\r
	120	enum kDetectorUsedForPID { kTPCpid, kTOFpid, kTPCTOF }; // default TPC & TOF pid (via GetTPCpid & GetTOFpid) \r
	121	enum kParticleOfInterest { kMuon, kElectron, kPion, kKaon, kProton }; \r
	122	\r
	123	void SetUseBayesianPID(Double_t gMinProbabilityValue) {\r
	124	fUsePID = kTRUE; fUsePIDnSigma = kFALSE; fUsePIDPropabilities = kTRUE;\r
	125	fMinAcceptedPIDProbability = gMinProbabilityValue; }\r
	126	\r
	127	void SetUseNSigmaPID(Double_t gMaxNSigma) {\r
	128	fUsePID = kTRUE; fUsePIDPropabilities = kFALSE; fUsePIDnSigma = kTRUE;\r
	129	fPIDNSigma = gMaxNSigma; }\r
	130	\r
	131	void SetParticleOfInterest(kParticleOfInterest poi) {\r
	132	fParticleOfInterest = poi;}\r
	133	void SetDetectorUsedForPID(kDetectorUsedForPID detConfig) {\r
	134	fPidDetectorConfig = detConfig;}\r
	135	\r
	136	private:\r
f06d59b3	137	Double_t IsEventAccepted(AliVEvent* event);\r
	138	Double_t GetEventPlane(AliVEvent* event);\r
	139	TObjArray* GetAcceptedTracks(AliVEvent* event, Double_t fCentrality, Double_t gReactionPlane);\r
	140	TObjArray* GetShuffledTracks(TObjArray* tracks);\r
	141	\r
0879e280	142	AliBalancePsi *fBalance; //BF object\r
	143	Bool_t fRunShuffling;//run shuffling or not\r
	144	AliBalancePsi *fShuffledBalance; //BF object (shuffled)\r
f06d59b3	145	Bool_t fRunMixing;//run mixing or not\r
	146	Int_t fMixingTracks;\r
	147	AliBalancePsi *fMixedBalance; //TriggeredBF object (mixed)\r
	148	AliEventPoolManager* fPoolMgr; //! event pool manager\r
	149	\r
0879e280	150	TList *fList; //fList object\r
	151	TList *fListBF; //fList object\r
	152	TList *fListBFS; //fList object\r
f06d59b3	153	TList *fListBFM; //fList object\r
0879e280	154	TList *fHistListPIDQA; //! list of histograms\r
0879e280	155	\r
6b046844	156	TH2F *fHistEventStats; //event stats\r
0879e280	157	TH2F *fHistCentStats; //centrality stats\r
	158	TH1F *fHistTriggerStats; //trigger stats\r
	159	TH1F *fHistTrackStats; //Track filter bit stats\r
	160	TH1F *fHistVx; //x coordinate of the primary vertex\r
	161	TH1F *fHistVy; //y coordinate of the primary vertex\r
6b046844	162	TH2F *fHistVz; //z coordinate of the primary vertex\r
0879e280	163	\r
6b046844	164	TH2F *fHistEventPlane; //event plane distribution\r
0879e280	165	\r
9fd4b54e	166	TH2F *fHistClus;//number of clusters (QA histogram)\r
	167	TH2F *fHistDCA;//DCA (QA histogram)\r
	168	TH2F *fHistChi2;//track chi2 (QA histogram)\r
	169	TH2F *fHistPt;//transverse momentum (QA histogram)\r
	170	TH2F *fHistEta;//pseudorapidity (QA histogram)\r
	171	TH2F *fHistRapidity;//rapidity (QA histogram)\r
	172	TH2F *fHistPhi;//phi (QA histogram)\r
	173	TH2F *fHistPhiBefore;//phi before v2 afterburner (QA histogram)\r
	174	TH2F *fHistPhiAfter;//phi after v2 afterburner (QA histogram)\r
	175	TH2F *fHistPhiPos;//phi for positive particles (QA histogram)\r
	176	TH2F *fHistPhiNeg;//phi for negative particles (QA histogram)\r
	177	TH2F *fHistV0M;//V0 multiplicities (QA histogram)\r
	178	TH2F *fHistRefTracks;//reference track multiplicities (QA histogram)\r
0879e280	179	\r
0879e280	180	//============PID============//\r
9fd4b54e	181	TH2D *fHistdEdxVsPTPCbeforePID;//TPC dEdx vs momentum before PID cuts (QA histogram)\r
	182	TH2D *fHistBetavsPTOFbeforePID;//beta vs momentum before PID cuts (QA histogram)\r
	183	TH2D *fHistProbTPCvsPtbeforePID; //TPC probability vs pT before PID cuts (QA histogram)\r
	184	TH2D *fHistProbTOFvsPtbeforePID;//TOF probability vs pT before PID cuts (QA histogram)\r
	185	TH2D *fHistProbTPCTOFvsPtbeforePID;//TOF/TPC probability vs pT before PID cuts (QA histogram)\r
	186	TH2D *fHistNSigmaTPCvsPtbeforePID;//TPC nsigma vs pT before PID cuts (QA histogram)\r
	187	TH2D *fHistNSigmaTOFvsPtbeforePID;//TOF nsigma vs pT before PID cuts (QA histogram)\r
	188	TH2D *fHistdEdxVsPTPCafterPID;//TPC dEdx vs momentum after PID cuts (QA histogram)\r
	189	TH2D *fHistBetavsPTOFafterPID;//beta vs momentum after PID cuts (QA histogram)\r
	190	TH2D *fHistProbTPCvsPtafterPID; //TPC probability vs pT after PID cuts (QA histogram)\r
	191	TH2D *fHistProbTOFvsPtafterPID;//TOF probability vs pT after PID cuts (QA histogram)\r
	192	TH2D *fHistProbTPCTOFvsPtafterPID;//TOF/TPC probability vs pT after PID cuts (QA histogram)\r
	193	TH2D *fHistNSigmaTPCvsPtafterPID;//TPC nsigma vs pT after PID cuts (QA histogram)\r
	194	TH2D *fHistNSigmaTOFvsPtafterPID;//TOF nsigma vs pT after PID cuts (QA histogram)\r
0879e280	195	\r
	196	AliPIDResponse *fPIDResponse; //! PID response object\r
	197	AliPIDCombined *fPIDCombined; //! combined PID object\r
	198	\r
9fd4b54e	199	kParticleOfInterest fParticleOfInterest;//analyzed particle\r
	200	kDetectorUsedForPID fPidDetectorConfig;//used detector for PID\r
	201	\r
	202	Bool_t fUsePID; //flag to use PID \r
	203	Bool_t fUsePIDnSigma;//flag to use nsigma method for PID\r
	204	Bool_t fUsePIDPropabilities;//flag to use probability method for PID\r
	205	Double_t fPIDNSigma;//nsigma cut for PID\r
	206	Double_t fMinAcceptedPIDProbability;//probability cut for PID\r
0879e280	207	//============PID============//\r
	208	\r
	209	AliESDtrackCuts *fESDtrackCuts; //ESD track cuts\r
	210	\r
	211	TString fCentralityEstimator; //"V0M","TRK","TKL","ZDC","FMD"\r
	212	Bool_t fUseCentrality;//use the centrality (PbPb) or not (pp)\r
	213	Double_t fCentralityPercentileMin;//centrality percentile min\r
	214	Double_t fCentralityPercentileMax;//centrality percentile max\r
	215	Double_t fImpactParameterMin;//impact parameter min (used for MC)\r
	216	Double_t fImpactParameterMax;//impact parameter max (used for MC)\r
	217	\r
	218	Bool_t fUseMultiplicity;//use the multiplicity cuts\r
	219	Int_t fNumberOfAcceptedTracksMin;//min. number of number of accepted tracks (used for the multiplicity dependence study - pp)\r
	220	Int_t fNumberOfAcceptedTracksMax;//max. number of number of accepted tracks (used for the multiplicity dependence study - pp)\r
6b046844	221	TH2F *fHistNumberOfAcceptedTracks;//hisot to store the number of accepted tracks\r
0879e280	222	\r
	223	Bool_t fUseOfflineTrigger;//Usage of the offline trigger selection\r
	224	\r
	225	Double_t fVxMax;//vxmax\r
	226	Double_t fVyMax;//vymax\r
	227	Double_t fVzMax;//vzmax\r
	228	\r
	229	Int_t nAODtrackCutBit;//track cut bit from track selection (only used for AODs)\r
	230	\r
	231	Double_t fPtMin;//only used for AODs\r
	232	Double_t fPtMax;//only used for AODs\r
	233	Double_t fEtaMin;//only used for AODs\r
	234	Double_t fEtaMax;//only used for AODs\r
	235	\r
	236	Double_t fDCAxyCut;//only used for AODs\r
	237	Double_t fDCAzCut;//only used for AODs\r
	238	\r
	239	Double_t fTPCchi2Cut;//only used for AODs\r
	240	Int_t fNClustersTPCCut;//only used for AODs\r
	241	\r
	242	TF1 *fAcceptanceParameterization;//acceptance filter used for MC\r
	243	\r
	244	TF1 *fDifferentialV2;//pt-differential v2 (from real data)\r
	245	Bool_t fUseFlowAfterBurner;//Usage of a flow after burner\r
	246	\r
	247	Bool_t fExcludeResonancesInMC;//flag to exclude the resonances' decay products from the MC analysis\r
	248	Bool_t fUseMCPdgCode; //Boolean to analyze a set of particles in MC\r
	249	Int_t fPDGCodeToBeAnalyzed; //Analyze a set of particles in MC\r
	250	\r
	251	\r
	252	\r
	253	AliAnalysisTaskBFPsi(const AliAnalysisTaskBFPsi&); // not implemented\r
	254	AliAnalysisTaskBFPsi& operator=(const AliAnalysisTaskBFPsi&); // not implemented\r
	255	\r
	256	ClassDef(AliAnalysisTaskBFPsi, 5); // example of analysis\r
	257	};\r
	258	\r
f06d59b3	259	\r
f06d59b3	260	\r
0879e280	261	#endif\r