[u/mrichter/AliRoot.git] / PWGHF / hfe / AliAnalysisTaskHFECal.h

/**************************************************************************
 * Copyright(c) 1998-1999, 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.                  *               
 **************************************************************************/

#ifndef ALIANALYSISTASKHFECAL_H
#define ALIANALYSISTASKHFECAL_H

class THnSparse;
class TH2F;
class TLorentzVector;

class AliEMCALTrack;
class AliMagF;
class AliESDEvent;
class AliAODEvent;
class AliEMCALGeometry;
class AliEMCALRecoUtils;
class AliAnalysisFilter;
class AliESDtrackCuts;
class AliESDtrack;
class AliHFEcontainer;
class AliHFEcuts;
class AliHFEpid;
class AliHFEpidQAmanager;
class AliCFManager;

#include "AliAnalysisTaskSE.h"

class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
 public:
  AliAnalysisTaskHFECal();
  AliAnalysisTaskHFECal(const char *name);
  virtual ~AliAnalysisTaskHFECal();
  
  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
  virtual void   Terminate(Option_t *);
  
  void SetHFECuts(AliHFEcuts * const cuts) { fCuts = cuts; };
  void SetOpeningAngleCut (Double_t openingAngle) {fOpeningAngleCut = openingAngle;};
  void SetInvariantMassCut (Double_t invmass) {fInvmassCut = invmass;};
  void SetQAHist (int qahist) {fqahist = qahist;};
  AliHFEpid *GetPID() const { return fPID; }
  void SetRejectKinkMother(Bool_t rejectKinkMother = kFALSE) { fRejectKinkMother = rejectKinkMother; };
  void SelectPhotonicElectron(Int_t itrack, Double_t cent, AliESDtrack *track, Bool_t &fFlagPhotonicElec, Bool_t &fFlagConvinatElec, Double_t nSig, Double_t shower, Double_t ep, Double_t mce);
  void FindTriggerClusters();
 private:
  
  Bool_t ProcessCutStep(Int_t cutStep, AliVParticle *track);
  
  AliESDEvent 		*fESD;			//!ESD object
  AliMCEvent 		*fMC;			//!MC object
  AliEMCALGeometry  	*fGeom; 		// emcal geometry 
    
  TList       		*fOutputList;		//! output list
  Int_t                 fqahist;  

  AliESDtrackCuts	*fTrackCuts;		//! ESD track cuts
  AliHFEcuts 		*fCuts;                 //! Cut Collection
  Bool_t 		fIdentifiedAsOutInz;    //Out Of Range in z
  Bool_t 		fPassTheEventCut;       //Pass The Event Cut
  Bool_t 		fRejectKinkMother;      //Reject Kink Mother
  Bool_t                fmcData;
  Double_t 		fVz;                    //z position of the primary vertex
  AliCFManager 		*fCFM;                  //! Correction Framework Manager
  AliHFEpid 		*fPID;                  //! PID
  AliHFEpidQAmanager 	*fPIDqa;		//! PID QA manager
  Double_t 		fOpeningAngleCut;	//openingAngle cut value
  Double_t		fInvmassCut;		//invariant mass cut value
 
  int ftriggers[48][60];//!
  int ftriggersCut[48][60];//!
  int ftriggersTime[48][60];//!
 

  TH1F			*fNoEvents;		//! no of events
  THnSparseD		*fEMCAccE;		//! EMC acc
  TH1F			*fTrkpt;		//! track pt
  TH2F			*fTrkEovPBef;		//! track E/p before HFE pid
  TH2F			*fTrkEovPAft;		//! track E/p after HFE pid
  TH2F			*fdEdxBef;		//! track dEdx vs p before HFE pid
  TH2F			*fdEdxAft;		//! track dEdx vs p after HFE pid
  TH2F			*fIncpT;		//! HFE pid electron vs centrality
  TH2F			*fIncpTM20;		//! HFE pid electron vs centrality
  THnSparseD		*fInvmassLS;		//! Inv mass of LS (e,e)
  THnSparseD		*fInvmassULS;		//! Inv mass of ULS (e,e)
  TH1F			*fOpeningAngleLS;	//! opening angle for LS pairs
  TH1F			*fOpeningAngleULS;	//! opening angle for ULS pairs
  TH1F			*fPhotoElecPt;		//! photonic elec pt 
  TH2F			*fPhoElecPt;	        //! Pho inclusive ele pt
  TH2F			*fPhoElecPtM20;	        //! Pho inclusive ele pt
  TH2F			*fSameElecPt;	        //! Same inclusive ele pt
  TH2F			*fSameElecPtM20;	        //! Same inclusive ele pt

  TH1F			*fTrackPtBefTrkCuts;	//! Track pt before track cuts	
  TH1F			*fTrackPtAftTrkCuts;	//! Track pt after track cuts
  TH2F			*fTPCnsigma;		//! TPC n sigma vs p	
  
  TH1F			*fCent;			//! centrality
  THnSparseD		*fEleInfo;		//! EMC acc
  /*
  //<---- trigger info
  TH1F	      *fClsEBftTrigCut;	//Cluster E before trigger selection
  TH1F        *fClsEAftTrigCut;	//Cluster E after trigger selection
  TH1F	      *fClsEAftTrigCut1;	//Cluster E after trigger selection
  TH1F	      *fClsEAftTrigCut2;	//Cluster E after trigger selection
  TH1F	      *fClsEAftTrigCut3;	//Cluster E after trigger selection
  TH1F	      *fClsEAftTrigCut4;	//Cluster E after trigger selection
  TH2F        *fClsETime; //ClsE vs time distribution
  TH2F        *fClsETime1; //ClsE vs time distribution
  TH1F        *fTrigTimes;// trigger time
  TH2F        *fCellCheck;// trigger time
  */
  //<------ MC
  TH2F                  *fInputHFEMC;
  TH2F                  *fInputAlle;
  TH2F			*fIncpTMChfe;		//! MC HFE pid electron vs centrality
  TH2F			*fIncpTMChfeAll;		//! MC HFE pid electron vs centrality
  TH2F			*fIncpTMCM20hfe;	//! MC HFE pid electron vs centrality
  TH2F			*fIncpTMCM20hfeAll;	//! MC HFE pid electron vs centrality
 THnSparseD		*fIncpTMCpho;		//! MC HFE pid electron vs centrality
 THnSparseD 		*fIncpTMCM20pho;	//! MC HFE pid electron vs centrality
 THnSparseD 		*fPhoElecPtMC;	        //! Pho inclusive ele pt
 THnSparseD 		*fPhoElecPtMCM20;	        //! Pho inclusive ele pt
 THnSparseD 		*fSameElecPtMC;	        //! Same inclusive ele pt
 THnSparseD 		*fSameElecPtMCM20;	        //! Same inclusive ele pt
 TH1D                  *CheckNclust;  
 TH1D                  *CheckNits;  
  
  AliAnalysisTaskHFECal(const AliAnalysisTaskHFECal&); // not implemented
  AliAnalysisTaskHFECal& operator=(const AliAnalysisTaskHFECal&); // not implemented
  
  ClassDef(AliAnalysisTaskHFECal, 1); //!example of analysis
};

#endif
Commit	Line	Data
bfc7c23b	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	#ifndef ALIANALYSISTASKHFECAL_H
	17	#define ALIANALYSISTASKHFECAL_H
	18
	19	class THnSparse;
	20	class TH2F;
	21	class TLorentzVector;
	22
	23	class AliEMCALTrack;
	24	class AliMagF;
	25	class AliESDEvent;
	26	class AliAODEvent;
	27	class AliEMCALGeometry;
	28	class AliEMCALRecoUtils;
	29	class AliAnalysisFilter;
	30	class AliESDtrackCuts;
	31	class AliESDtrack;
	32	class AliHFEcontainer;
	33	class AliHFEcuts;
	34	class AliHFEpid;
	35	class AliHFEpidQAmanager;
	36	class AliCFManager;
	37
	38	#include "AliAnalysisTaskSE.h"
	39
	40	class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
	41	public:
	42	AliAnalysisTaskHFECal();
	43	AliAnalysisTaskHFECal(const char *name);
	44	virtual ~AliAnalysisTaskHFECal();
	45
	46	virtual void UserCreateOutputObjects();
	47	virtual void UserExec(Option_t *option);
	48	virtual void Terminate(Option_t *);
	49
	50	void SetHFECuts(AliHFEcuts * const cuts) { fCuts = cuts; };
	51	void SetOpeningAngleCut (Double_t openingAngle) {fOpeningAngleCut = openingAngle;};
	52	void SetInvariantMassCut (Double_t invmass) {fInvmassCut = invmass;};
85985bb0	53	void SetQAHist (int qahist) {fqahist = qahist;};
bfc7c23b	54	AliHFEpid *GetPID() const { return fPID; }
bfc7c23b	55	void SetRejectKinkMother(Bool_t rejectKinkMother = kFALSE) { fRejectKinkMother = rejectKinkMother; };
c2801a73	56	void SelectPhotonicElectron(Int_t itrack, Double_t cent, AliESDtrack *track, Bool_t &fFlagPhotonicElec, Bool_t &fFlagConvinatElec, Double_t nSig, Double_t shower, Double_t ep, Double_t mce);
fb87d707	57	void FindTriggerClusters();
bfc7c23b	58	private:
	59
	60	Bool_t ProcessCutStep(Int_t cutStep, AliVParticle *track);
	61
	62	AliESDEvent *fESD; //!ESD object
e7d87aef	63	AliMCEvent *fMC; //!MC object
fb87d707	64	AliEMCALGeometry *fGeom; // emcal geometry
bfc7c23b	65
bfc7c23b	66	TList *fOutputList; //! output list
85985bb0	67	Int_t fqahist;
85985bb0	68
bfc7c23b	69	AliESDtrackCuts *fTrackCuts; //! ESD track cuts
	70	AliHFEcuts *fCuts; //! Cut Collection
	71	Bool_t fIdentifiedAsOutInz; //Out Of Range in z
	72	Bool_t fPassTheEventCut; //Pass The Event Cut
	73	Bool_t fRejectKinkMother; //Reject Kink Mother
e7d87aef	74	Bool_t fmcData;
bfc7c23b	75	Double_t fVz; //z position of the primary vertex
	76	AliCFManager *fCFM; //! Correction Framework Manager
	77	AliHFEpid *fPID; //! PID
	78	AliHFEpidQAmanager *fPIDqa; //! PID QA manager
	79	Double_t fOpeningAngleCut; //openingAngle cut value
	80	Double_t fInvmassCut; //invariant mass cut value
fb87d707	81
	82	int ftriggers[48][60];//!
	83	int ftriggersCut[48][60];//!
	84	int ftriggersTime[48][60];//!
	85
e7d87aef	86
bfc7c23b	87	TH1F *fNoEvents; //! no of events
	88	THnSparseD *fEMCAccE; //! EMC acc
	89	TH1F *fTrkpt; //! track pt
	90	TH2F *fTrkEovPBef; //! track E/p before HFE pid
	91	TH2F *fTrkEovPAft; //! track E/p after HFE pid
	92	TH2F *fdEdxBef; //! track dEdx vs p before HFE pid
	93	TH2F *fdEdxAft; //! track dEdx vs p after HFE pid
	94	TH2F *fIncpT; //! HFE pid electron vs centrality
fb87d707	95	TH2F *fIncpTM20; //! HFE pid electron vs centrality
f09766dd	96	THnSparseD *fInvmassLS; //! Inv mass of LS (e,e)
f09766dd	97	THnSparseD *fInvmassULS; //! Inv mass of ULS (e,e)
bfc7c23b	98	TH1F *fOpeningAngleLS; //! opening angle for LS pairs
	99	TH1F *fOpeningAngleULS; //! opening angle for ULS pairs
	100	TH1F *fPhotoElecPt; //! photonic elec pt
f09766dd	101	TH2F *fPhoElecPt; //! Pho inclusive ele pt
fb87d707	102	TH2F *fPhoElecPtM20; //! Pho inclusive ele pt
f09766dd	103	TH2F *fSameElecPt; //! Same inclusive ele pt
fb87d707	104	TH2F *fSameElecPtM20; //! Same inclusive ele pt
a6df418c	105
bfc7c23b	106	TH1F *fTrackPtBefTrkCuts; //! Track pt before track cuts
	107	TH1F *fTrackPtAftTrkCuts; //! Track pt after track cuts
	108	TH2F *fTPCnsigma; //! TPC n sigma vs p
	109
	110	TH1F *fCent; //! centrality
	111	THnSparseD *fEleInfo; //! EMC acc
feffe705	112	/*
fb87d707	113	//<---- trigger info
	114	TH1F *fClsEBftTrigCut; //Cluster E before trigger selection
	115	TH1F *fClsEAftTrigCut; //Cluster E after trigger selection
	116	TH1F *fClsEAftTrigCut1; //Cluster E after trigger selection
	117	TH1F *fClsEAftTrigCut2; //Cluster E after trigger selection
	118	TH1F *fClsEAftTrigCut3; //Cluster E after trigger selection
	119	TH1F *fClsEAftTrigCut4; //Cluster E after trigger selection
	120	TH2F *fClsETime; //ClsE vs time distribution
	121	TH2F *fClsETime1; //ClsE vs time distribution
	122	TH1F *fTrigTimes;// trigger time
42c75692	123	TH2F *fCellCheck;// trigger time
feffe705	124	*/
e7d87aef	125	//<------ MC
e7d87aef	126	TH2F *fInputHFEMC;
feffe705	127	TH2F *fInputAlle;
e7d87aef	128	TH2F *fIncpTMChfe; //! MC HFE pid electron vs centrality
3db00c72	129	TH2F *fIncpTMChfeAll; //! MC HFE pid electron vs centrality
e7d87aef	130	TH2F *fIncpTMCM20hfe; //! MC HFE pid electron vs centrality
3db00c72	131	TH2F *fIncpTMCM20hfeAll; //! MC HFE pid electron vs centrality
44be9e1d	132	THnSparseD *fIncpTMCpho; //! MC HFE pid electron vs centrality
	133	THnSparseD *fIncpTMCM20pho; //! MC HFE pid electron vs centrality
	134	THnSparseD *fPhoElecPtMC; //! Pho inclusive ele pt
	135	THnSparseD *fPhoElecPtMCM20; //! Pho inclusive ele pt
	136	THnSparseD *fSameElecPtMC; //! Same inclusive ele pt
	137	THnSparseD *fSameElecPtMCM20; //! Same inclusive ele pt
a6df418c	138	TH1D *CheckNclust;
a6df418c	139	TH1D *CheckNits;
feffe705	140
bfc7c23b	141	AliAnalysisTaskHFECal(const AliAnalysisTaskHFECal&); // not implemented
	142	AliAnalysisTaskHFECal& operator=(const AliAnalysisTaskHFECal&); // not implemented
	143
	144	ClassDef(AliAnalysisTaskHFECal, 1); //!example of analysis
	145	};
	146
	147	#endif
	148
	149