[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliAnalysisTaskJetServices.h

#ifndef ALIANALYSISTASKJETSERVICES_H
#define ALIANALYSISTASKJETSERVICES_H
 
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

// **************************************
// Task used for the correction of determiantion of reconstructed jet spectra
// Compares input (gen) and output (rec) jets   
// *******************************************

#include  "AliAnalysisTaskSE.h"
#include  "THnSparse.h" // cannot forward declare ThnSparseF

////////////////
class AliJetHeader;
class AliESDEvent;
class AliESDVertex;
class AliAODEvent;
class AliAODVertex;
class AliAODJet;
class AliGenPythiaEventHeader;
class AliCFManager;
class AliTriggerAnalysis;

class TList;
class TClonesArray;
class TChain;
class TH1F;
class TH2F;
class TH3F;
class TProfile;


class AliAnalysisTaskJetServices : public AliAnalysisTaskSE
{
 public:
    AliAnalysisTaskJetServices();
    AliAnalysisTaskJetServices(const char* name);
    virtual ~AliAnalysisTaskJetServices();
    // 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);
    virtual void SetZVertexCut(Float_t f){fVtxZCut = f;}
    virtual void SetPtMinCosmic(Float_t ptMin) {fPtMinCosmic = ptMin;}
    virtual void SetRMinCosmic(Float_t rMin) {fRIsolMinCosmic = rMin;}
    virtual void SetMaxCosmicAngle(Float_t angle)   {fMaxCosmicAngle = angle;}
    virtual Bool_t Notify();

    virtual void SetAODInput(Bool_t b){fUseAODInput = b;}
    virtual void SetRunRange(Float_t fLo,Float_t fUp){fRunRange[0] = fLo;fRunRange[1] = fUp;}
    virtual void SetMCData(Bool_t b){fMC = b;}
    virtual void SetUsePhysicsSelection(Bool_t b){fUsePhysicsSelection = b;}
    virtual void SetPhysicsSelectionFlag(Int_t i){fPhysicsSelectionFlag = i;}
    virtual void SetFilterAODCollisions(Bool_t b){fFilterAODCollisions = b;}

    virtual void SetNonStdFile(char *c){fNonStdFile = c;}
    Bool_t IsEventSelected(const AliESDEvent* esd);
    Bool_t IsEventSelected(const AliAODEvent* aod) const;

    Bool_t IsEventPileUp(const AliESDEvent* esd) const;

    Bool_t IsEventCosmic(const AliESDEvent* esd) const;

    Bool_t IsVertexValid(const AliESDVertex *vtx);
    Bool_t IsVertexValid(const AliAODVertex *vtx) const;

    Bool_t IsVertexIn(const AliESDVertex *vtx);
    Bool_t IsVertexIn(const AliAODVertex *vtx) const;
    Int_t GetEventClass(AliESDEvent *esd);
    Int_t GetEventClass(AliAODEvent *aod);


    virtual void SetFilterMask(UInt_t i){fFilterMask = i;}
    virtual void SetMinTrackPt(Float_t f){fMinTrackPt = f;}
    virtual void SetTrackEtaWindow(Float_t f){fTrackRecEtaWindow = f;}

    virtual void SetPhiWeights(TH3F *phiw){fh3PhiWeights = phiw;}
    virtual void SetFlatteningCoeff(Float_t *fA,Float_t *fB){
      fFlatA[0] = fA[0];fFlatA[1] = fA[1];
      fFlatA[0] = fB[0];fFlatB[1] = fB[1];
    }
    virtual void SetDeltaQxy(Float_t *fD){
      fDeltaQxy[0] = fD[0];
      fDeltaQxy[1] = fD[1];
    }

    Bool_t   CalculateReactionPlaneAngle(const TList *trackList);
    Double_t GetPhiWeight(Double_t phi,Double_t signedpt);
    Int_t   GetListOfTracks(TList *list);

    enum { kAllTriggered = 0,kTriggeredVertex,kTriggeredVertexIn,kSelectedALICE,kSelectedALICEVertexValid,kSelectedALICEVertexIn,kSelected,kConstraints};

    enum { kNoEventCut=1<<0,
	   kPhysicsSelectionCut=1<<1,
	   kContributorsCut1=1<<2,
	   kContributorsCut2=1<<3,
	   kContributorsCut3=1<<4,
	   kVertexTPC=1<<5,
	   kVertexSPD=1<<6,
	   kVertexGlobal=1<<7,
	   kSPDDispersionCut=1<<8,
	   kVertexZCut=1<<9,
	   kVertexRCut=1<<10,
	   kTotalEventCuts=(1<<11)-1};

 private:

    AliAnalysisTaskJetServices(const AliAnalysisTaskJetServices&);
    AliAnalysisTaskJetServices& operator=(const AliAnalysisTaskJetServices&);

    Bool_t        fUseAODInput;        // take jet from input AOD not from ouptu AOD
    Bool_t        fUsePhysicsSelection;// decide wether we take into account physicsselction task
    Bool_t        fMC;                 // true for MC data to allow correct trigger slection
    Bool_t        fFilterAODCollisions; // filter out collision canditates to the AOD
    UInt_t        fPhysicsSelectionFlag; // defines the glag for acceptance of events from physics selection
    UInt_t        fSelectionInfoESD;   // slection info bit mask
    UInt_t        fEventCutInfoESD;   // event selection info of what is cutted after physics selection
    UInt_t        fFilterMask;         // filter bit for slecected tracks
    Int_t         fRPSubeventMethod;   // method for subevent calculation
    Float_t       fAvgTrials;          // Average number of trials
    Float_t       fVtxXMean;           // mean x for cuts
    Float_t       fVtxYMean;           // mean y for cuts
    Float_t       fVtxZMean;           // mean z for cuts
    Float_t       fVtxRCut;            // vtx cut in R
    Float_t       fVtxZCut;            // vtzx cut in Z
    Float_t       fPtMinCosmic;        // Minimum pT to be considered as cosmic candidate
    Float_t       fRIsolMinCosmic;     // Minimum R = sqrt{deltaPhi^2 + deltaEta^2} to be considered as cosmic candidate
    Float_t       fMaxCosmicAngle;     // Max deviation from pi (angle between two tracks) in case of cosmic candidate
    Float_t       fRunRange[2];        // only important for real data for 
    Float_t       fCentrality;         // ! centrality
    Float_t       fTrackRecEtaWindow;     // eta window for rec tracks
    Float_t       fMinTrackPt;            // limits the track p_T 
    Float_t       fRPAngle;               // ! RP angle of the reaction plane
    Float_t       fFlatA[2];              // flattening for RP
    Float_t       fFlatB[2];              // flattening for RP
    Float_t       fDeltaQxy[2];           // centering of QX QY

    TRandom3      *fRandomizer;        // ! randomizer

    TString       fNonStdFile;         // outputName for replication
    TProfile*     fh1Xsec;             //! pythia cross section and trials
    TH1F*         fh1Trials;           //! trials are added
    TH1F*         fh1PtHard;           //! Pt har of the event...       
    TH1F*         fh1PtHardTrials;     //! Number of trials 
    TH1F*         fh1SelectionInfoESD; //! Masks that satisfy fSelectionInfo
    TH1F*         fh1EventCutInfoESD;  //! Masks that satisfy fSelectionInfo
    TH1F*         fh1CentralityESD;    //! centrality 
    TH1F*         fh1Centrality;       //! centrality 
    TH1F*         fh1RP;               //! RP distribution
    TH2F*         fh2TriggerCount;     //! number of fire triggers in each case
    TH2F*         fh2ESDTriggerCount;  //! number of fire triggers in each case
    TH2F*         fh2TriggerVtx;       //! vtx. position vs. trigger decision
    TH2F*         fh2ESDTriggerVtx;    //! vtx. position vs. trigger decision 
    TH2F*         fh2ESDTriggerRun;    //! fired triggers vs. run number
    TH2F*         fh2VtxXY;            //! XY position of VTX were available
    TH1F*         fh1NCosmicsPerEvent;  //! Number of coscmic candidates found in event
    TH2F*         fh2RPSubevents;   //! subevent RP 
    TH2F*         fh2RPCentrality;   //! RP vs centrality
    TH2F*         fh2RPDeltaRP;     //! centrality vs. RP dela  
    TH2F*         fh2RPQxQy;        //! QX QY moments
    TH2F*         fh2RPCosDeltaRP;  //! RP resolution

    TH3F*         fh3PhiWeights;  // RP phi weights, need to be set externally
    TH3F*         fh3RPPhiTracks; //! RP angle

    AliTriggerAnalysis *fTriggerAnalysis; //! Trigger Analysis to get the background rates etc.
    TList *fHistList; //! Output list

        // Provisions for replication
    static AliAODHeader*    fgAODHeader;        //! Header for replication
    static TClonesArray*  fgAODVertices;        //! primary vertex for replication
    ClassDef(AliAnalysisTaskJetServices,12)
};
 
#endif
Commit	Line	Data
bf7b8731	1	#ifndef ALIANALYSISTASKJETSERVICES_H
	2	#define ALIANALYSISTASKJETSERVICES_H
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	// **************************************
	8	// Task used for the correction of determiantion of reconstructed jet spectra
	9	// Compares input (gen) and output (rec) jets
	10	// *******************************************
	11
	12	#include "AliAnalysisTaskSE.h"
	13	#include "THnSparse.h" // cannot forward declare ThnSparseF
	14
	15	////////////////
	16	class AliJetHeader;
	17	class AliESDEvent;
6d597ca2	18	class AliESDVertex;
bf7b8731	19	class AliAODEvent;
6d597ca2	20	class AliAODVertex;
bf7b8731	21	class AliAODJet;
	22	class AliGenPythiaEventHeader;
	23	class AliCFManager;
3493c3a9	24	class AliTriggerAnalysis;
bf7b8731	25
bf7b8731	26	class TList;
80a790fd	27	class TClonesArray;
bf7b8731	28	class TChain;
c0997643	29	class TH1F;
bf7b8731	30	class TH2F;
	31	class TH3F;
	32	class TProfile;
	33
	34
	35
	36	class AliAnalysisTaskJetServices : public AliAnalysisTaskSE
	37	{
	38	public:
	39	AliAnalysisTaskJetServices();
	40	AliAnalysisTaskJetServices(const char* name);
80a790fd	41	virtual ~AliAnalysisTaskJetServices();
bf7b8731	42	// Implementation of interface methods
	43	virtual void UserCreateOutputObjects();
	44	virtual void Init();
	45	virtual void LocalInit() { Init(); }
	46	virtual void UserExec(Option_t *option);
	47	virtual void Terminate(Option_t *option);
c0997643	48	virtual void SetZVertexCut(Float_t f){fVtxZCut = f;}
d49d814a	49	virtual void SetPtMinCosmic(Float_t ptMin) {fPtMinCosmic = ptMin;}
	50	virtual void SetRMinCosmic(Float_t rMin) {fRIsolMinCosmic = rMin;}
	51	virtual void SetMaxCosmicAngle(Float_t angle) {fMaxCosmicAngle = angle;}
bf7b8731	52	virtual Bool_t Notify();
	53
	54	virtual void SetAODInput(Bool_t b){fUseAODInput = b;}
b5a3f310	55	virtual void SetRunRange(Float_t fLo,Float_t fUp){fRunRange[0] = fLo;fRunRange[1] = fUp;}
3493c3a9	56	virtual void SetMCData(Bool_t b){fMC = b;}
cc0649e4	57	virtual void SetUsePhysicsSelection(Bool_t b){fUsePhysicsSelection = b;}
39e7e8ab	58	virtual void SetPhysicsSelectionFlag(Int_t i){fPhysicsSelectionFlag = i;}
f4132e7d	59	virtual void SetFilterAODCollisions(Bool_t b){fFilterAODCollisions = b;}
	60
	61	virtual void SetNonStdFile(char *c){fNonStdFile = c;}
6d597ca2	62	Bool_t IsEventSelected(const AliESDEvent* esd);
	63	Bool_t IsEventSelected(const AliAODEvent* aod) const;
	64
	65	Bool_t IsEventPileUp(const AliESDEvent* esd) const;
	66
	67	Bool_t IsEventCosmic(const AliESDEvent* esd) const;
bf7b8731	68
6d597ca2	69	Bool_t IsVertexValid(const AliESDVertex *vtx);
	70	Bool_t IsVertexValid(const AliAODVertex *vtx) const;
	71
	72	Bool_t IsVertexIn(const AliESDVertex *vtx);
	73	Bool_t IsVertexIn(const AliAODVertex *vtx) const;
f4132e7d	74	Int_t GetEventClass(AliESDEvent *esd);
ffab794c	75	Int_t GetEventClass(AliAODEvent *aod);
6d597ca2	76
40445651	77
	78	virtual void SetFilterMask(UInt_t i){fFilterMask = i;}
	79	virtual void SetMinTrackPt(Float_t f){fMinTrackPt = f;}
	80	virtual void SetTrackEtaWindow(Float_t f){fTrackRecEtaWindow = f;}
	81
	82	virtual void SetPhiWeights(TH3F *phiw){fh3PhiWeights = phiw;}
	83	virtual void SetFlatteningCoeff(Float_t fA,Float_t fB){
	84	fFlatA[0] = fA[0];fFlatA[1] = fA[1];
	85	fFlatA[0] = fB[0];fFlatB[1] = fB[1];
	86	}
	87	virtual void SetDeltaQxy(Float_t *fD){
	88	fDeltaQxy[0] = fD[0];
	89	fDeltaQxy[1] = fD[1];
	90	}
	91
	92	Bool_t CalculateReactionPlaneAngle(const TList *trackList);
	93	Double_t GetPhiWeight(Double_t phi,Double_t signedpt);
	94	Int_t GetListOfTracks(TList *list);
	95
6d597ca2	96	enum { kAllTriggered = 0,kTriggeredVertex,kTriggeredVertexIn,kSelectedALICE,kSelectedALICEVertexValid,kSelectedALICEVertexIn,kSelected,kConstraints};
f2dd0695	97
6d597ca2	98	enum { kNoEventCut=1<<0,
	99	kPhysicsSelectionCut=1<<1,
	100	kContributorsCut1=1<<2,
	101	kContributorsCut2=1<<3,
	102	kContributorsCut3=1<<4,
	103	kVertexTPC=1<<5,
	104	kVertexSPD=1<<6,
	105	kVertexGlobal=1<<7,
	106	kSPDDispersionCut=1<<8,
	107	kVertexZCut=1<<9,
	108	kVertexRCut=1<<10,
	109	kTotalEventCuts=(1<<11)-1};
f2dd0695	110
bf7b8731	111	private:
	112
	113	AliAnalysisTaskJetServices(const AliAnalysisTaskJetServices&);
	114	AliAnalysisTaskJetServices& operator=(const AliAnalysisTaskJetServices&);
	115
	116	Bool_t fUseAODInput; // take jet from input AOD not from ouptu AOD
cc0649e4	117	Bool_t fUsePhysicsSelection;// decide wether we take into account physicsselction task
3493c3a9	118	Bool_t fMC; // true for MC data to allow correct trigger slection
f4132e7d	119	Bool_t fFilterAODCollisions; // filter out collision canditates to the AOD
39e7e8ab	120	UInt_t fPhysicsSelectionFlag; // defines the glag for acceptance of events from physics selection
d49d814a	121	UInt_t fSelectionInfoESD; // slection info bit mask
6d597ca2	122	UInt_t fEventCutInfoESD; // event selection info of what is cutted after physics selection
40445651	123	UInt_t fFilterMask; // filter bit for slecected tracks
40445651	124	Int_t fRPSubeventMethod; // method for subevent calculation
bf7b8731	125	Float_t fAvgTrials; // Average number of trials
d49d814a	126	Float_t fVtxXMean; // mean x for cuts
	127	Float_t fVtxYMean; // mean y for cuts
	128	Float_t fVtxZMean; // mean z for cuts
	129	Float_t fVtxRCut; // vtx cut in R
	130	Float_t fVtxZCut; // vtzx cut in Z
	131	Float_t fPtMinCosmic; // Minimum pT to be considered as cosmic candidate
	132	Float_t fRIsolMinCosmic; // Minimum R = sqrt{deltaPhi^2 + deltaEta^2} to be considered as cosmic candidate
	133	Float_t fMaxCosmicAngle; // Max deviation from pi (angle between two tracks) in case of cosmic candidate
b5a3f310	134	Float_t fRunRange[2]; // only important for real data for
40445651	135	Float_t fCentrality; // ! centrality
	136	Float_t fTrackRecEtaWindow; // eta window for rec tracks
	137	Float_t fMinTrackPt; // limits the track p_T
	138	Float_t fRPAngle; // ! RP angle of the reaction plane
	139	Float_t fFlatA[2]; // flattening for RP
	140	Float_t fFlatB[2]; // flattening for RP
	141	Float_t fDeltaQxy[2]; // centering of QX QY
	142
	143	TRandom3 *fRandomizer; // ! randomizer
	144
f4132e7d	145	TString fNonStdFile; // outputName for replication
3493c3a9	146	TProfile* fh1Xsec; //! pythia cross section and trials
	147	TH1F* fh1Trials; //! trials are added
	148	TH1F* fh1PtHard; //! Pt har of the event...
	149	TH1F* fh1PtHardTrials; //! Number of trials
	150	TH1F* fh1SelectionInfoESD; //! Masks that satisfy fSelectionInfo
	151	TH1F* fh1EventCutInfoESD; //! Masks that satisfy fSelectionInfo
742ee86c	152	TH1F* fh1CentralityESD; //! centrality
742ee86c	153	TH1F* fh1Centrality; //! centrality
40445651	154	TH1F* fh1RP; //! RP distribution
3493c3a9	155	TH2F* fh2TriggerCount; //! number of fire triggers in each case
	156	TH2F* fh2ESDTriggerCount; //! number of fire triggers in each case
	157	TH2F* fh2TriggerVtx; //! vtx. position vs. trigger decision
	158	TH2F* fh2ESDTriggerVtx; //! vtx. position vs. trigger decision
	159	TH2F* fh2ESDTriggerRun; //! fired triggers vs. run number
	160	TH2F* fh2VtxXY; //! XY position of VTX were available
	161	TH1F* fh1NCosmicsPerEvent; //! Number of coscmic candidates found in event
40445651	162	TH2F* fh2RPSubevents; //! subevent RP
	163	TH2F* fh2RPCentrality; //! RP vs centrality
	164	TH2F* fh2RPDeltaRP; //! centrality vs. RP dela
	165	TH2F* fh2RPQxQy; //! QX QY moments
	166	TH2F* fh2RPCosDeltaRP; //! RP resolution
	167
	168	TH3F* fh3PhiWeights; // RP phi weights, need to be set externally
	169	TH3F* fh3RPPhiTracks; //! RP angle
	170
3493c3a9	171	AliTriggerAnalysis *fTriggerAnalysis; //! Trigger Analysis to get the background rates etc.
f4132e7d	172	TList *fHistList; //! Output list
	173
	174	// Provisions for replication
	175	static AliAODHeader* fgAODHeader; //! Header for replication
80a790fd	176	static TClonesArray* fgAODVertices; //! primary vertex for replication
80a790fd	177	ClassDef(AliAnalysisTaskJetServices,12)
bf7b8731	178	};
	179
	180	#endif