[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskCLQA.h

#ifndef ALIANALYSISTASKCLQA_H
#define ALIANALYSISTASKCLQA_H

// $Id: AliAnalysisTaskCLQA.h 58847 2012-09-30 18:11:49Z loizides $

class TClonesArray;
class TString;
class TH1F;
class TH2F;
class TH3F;
class TNtuple;
class TNtupleD;
class TTree;

#include "AliAnalysisTaskEmcalJet.h"

class AliNtupCumInfo;

class AliAnalysisTaskCLQA : public AliAnalysisTaskEmcalJet {
 public:
  AliAnalysisTaskCLQA();
  AliAnalysisTaskCLQA(const char *name);
  virtual ~AliAnalysisTaskCLQA();

  void                        UserCreateOutputObjects();
  void                        SetDoCumulants(Bool_t b)          { fDoCumulants = b; }
  void                        SetCumParams(Double_t Mmin, Double_t ptmin, Double_t ptmax, Double_t etamin, Double_t etamax);

 protected:
  Bool_t                      FillHistograms();
  Bool_t                      RetrieveEventObjects();
  Bool_t                      Run();
  void                        RunCumulants(Double_t Mmin, Double_t ptmin, Double_t ptmax, Double_t etamin, Double_t etamax);

  Bool_t                      fDoCumulants;      // if true run cumulant analysis
  Double_t                    fCumPtMin;         // minimum pt for cumulants
  Double_t                    fCumPtMax;         // maximum pt for cumulants
  Double_t                    fCumEtaMin;        // minimum eta for cumulants
  Double_t                    fCumEtaMax;        // maximum eta for cumulants
  Double_t                    fCumMmin;          // minimum number of tracks for cumulants 
  TTree                      *fNtupCum;          //!ntuple for cumulant analysis
  AliNtupCumInfo             *fNtupCumInfo;      //!object holding cumulant results

 private:
  AliAnalysisTaskCLQA(const AliAnalysisTaskCLQA&);            // not implemented
  AliAnalysisTaskCLQA &operator=(const AliAnalysisTaskCLQA&); // not implemented

  ClassDef(AliAnalysisTaskCLQA, 1) // Constantin's Task
};

class AliNtupCumInfo {
 public:
    AliNtupCumInfo() : fTrig(0), fRun(0), fVz(0), fIsFEC(0), fIsVSel(0), fIsP(0),
                       fMall(0), fMall2(0), fPtMaxall(0), fMPtall(0), 
                       fMPt2all(0), fMPtall2(0), fTSall(0),
                       fM(0), fQ2abs(0), fQ4abs(0), fQ42re(0),
                       fPtMax(0), fMPt(0), fMPt2(0), fTS(0), fMV0M(0), 
                       fCl1(0), fV0M(0), fV0MEq(0), fV0A(0), fV0AEq(0), fZNA(0) {;}
  virtual ~AliNtupCumInfo() {;}

 public:
  UInt_t        fTrig;         // trigger selection
  Int_t         fRun;          // run number
  Double_t      fVz;           // vertex z
  Bool_t        fIsFEC;        // is first event in chunk
  Bool_t        fIsVSel;       // is vertex selected
  Bool_t        fIsP;          // is SPD pileup
  Int_t         fMall;         // multiplicity (tracks in eta range)
  Int_t         fMall2;        // multiplicity (tracks above 2 GeV/c in eta range)
  Double32_t    fPtMaxall;     //[0,0,16] maximum pT
  Double32_t    fMPtall;       //[0,0,16] mean pT
  Double32_t    fMPt2all;      //[0,0,16] mean pT2
  Double32_t    fMPtall2;      //[0,0,16] mean pT truncated above 2 GeV/c
  Double32_t    fTSall;        //[0,0,16] transverse sphericity
  Int_t         fM;            // multiplicity (tracks in pT range)
  Double_t      fQ2abs;        // Q2 absolute
  Double_t      fQ4abs;        // Q4 absolute
  Double_t      fQ42re;        // Re(Q2Q*Q*)
  Double32_t    fPtMax;        //[0,0,16] maximum pT
  Double32_t    fMPt;          //[0,0,16] mean pT
  Double32_t    fMPt2;         //[0,0,16] mean pT2
  Double32_t    fTS;           //[0,0,16] transverse sphericity
  Double32_t    fMV0M;         // V0M amplitude
  Double32_t    fCl1;          //[0,0,16] class CL1
  Double32_t    fV0M;          //[0,0,16] class V0M
  Double32_t    fV0MEq;        //[0,0,16] class V0M Eq
  Double32_t    fV0A;          //[0,0,16] class V0A
  Double32_t    fV0AEq;        //[0,0,16] class V0A Eq
  Double32_t    fZNA;          //[0,0,16] class ZNA
 
  ClassDef(AliNtupCumInfo,2) // Cumulant storage class
};

#endif
Commit	Line	Data
d87c9c78	1	#ifndef ALIANALYSISTASKCLQA_H
	2	#define ALIANALYSISTASKCLQA_H
	3
	4	// $Id: AliAnalysisTaskCLQA.h 58847 2012-09-30 18:11:49Z loizides $
	5
	6	class TClonesArray;
	7	class TString;
	8	class TH1F;
	9	class TH2F;
	10	class TH3F;
96a87a98	11	class TNtuple;
	12	class TNtupleD;
	13	class TTree;
d87c9c78	14
	15	#include "AliAnalysisTaskEmcalJet.h"
	16
96a87a98	17	class AliNtupCumInfo;
96a87a98	18
d87c9c78	19	class AliAnalysisTaskCLQA : public AliAnalysisTaskEmcalJet {
	20	public:
	21	AliAnalysisTaskCLQA();
	22	AliAnalysisTaskCLQA(const char *name);
	23	virtual ~AliAnalysisTaskCLQA();
	24
	25	void UserCreateOutputObjects();
96a87a98	26	void SetDoCumulants(Bool_t b) { fDoCumulants = b; }
96a87a98	27	void SetCumParams(Double_t Mmin, Double_t ptmin, Double_t ptmax, Double_t etamin, Double_t etamax);
d87c9c78	28
d87c9c78	29	protected:
96a87a98	30	Bool_t FillHistograms();
	31	Bool_t RetrieveEventObjects();
	32	Bool_t Run();
	33	void RunCumulants(Double_t Mmin, Double_t ptmin, Double_t ptmax, Double_t etamin, Double_t etamax);
d87c9c78	34
96a87a98	35	Bool_t fDoCumulants; // if true run cumulant analysis
	36	Double_t fCumPtMin; // minimum pt for cumulants
	37	Double_t fCumPtMax; // maximum pt for cumulants
	38	Double_t fCumEtaMin; // minimum eta for cumulants
	39	Double_t fCumEtaMax; // maximum eta for cumulants
	40	Double_t fCumMmin; // minimum number of tracks for cumulants
	41	TTree *fNtupCum; //!ntuple for cumulant analysis
	42	AliNtupCumInfo *fNtupCumInfo; //!object holding cumulant results
d87c9c78	43
	44	private:
	45	AliAnalysisTaskCLQA(const AliAnalysisTaskCLQA&); // not implemented
	46	AliAnalysisTaskCLQA &operator=(const AliAnalysisTaskCLQA&); // not implemented
	47
	48	ClassDef(AliAnalysisTaskCLQA, 1) // Constantin's Task
	49	};
96a87a98	50
	51	class AliNtupCumInfo {
	52	public:
	53	AliNtupCumInfo() : fTrig(0), fRun(0), fVz(0), fIsFEC(0), fIsVSel(0), fIsP(0),
e29fb8cd	54	fMall(0), fMall2(0), fPtMaxall(0), fMPtall(0),
e29fb8cd	55	fMPt2all(0), fMPtall2(0), fTSall(0),
96a87a98	56	fM(0), fQ2abs(0), fQ4abs(0), fQ42re(0),
e29fb8cd	57	fPtMax(0), fMPt(0), fMPt2(0), fTS(0), fMV0M(0),
96a87a98	58	fCl1(0), fV0M(0), fV0MEq(0), fV0A(0), fV0AEq(0), fZNA(0) {;}
	59	virtual ~AliNtupCumInfo() {;}
	60
	61	public:
	62	UInt_t fTrig; // trigger selection
	63	Int_t fRun; // run number
	64	Double_t fVz; // vertex z
	65	Bool_t fIsFEC; // is first event in chunk
	66	Bool_t fIsVSel; // is vertex selected
	67	Bool_t fIsP; // is SPD pileup
	68	Int_t fMall; // multiplicity (tracks in eta range)
e29fb8cd	69	Int_t fMall2; // multiplicity (tracks above 2 GeV/c in eta range)
96a87a98	70	Double32_t fPtMaxall; //[0,0,16] maximum pT
	71	Double32_t fMPtall; //[0,0,16] mean pT
	72	Double32_t fMPt2all; //[0,0,16] mean pT2
e29fb8cd	73	Double32_t fMPtall2; //[0,0,16] mean pT truncated above 2 GeV/c
96a87a98	74	Double32_t fTSall; //[0,0,16] transverse sphericity
	75	Int_t fM; // multiplicity (tracks in pT range)
	76	Double_t fQ2abs; // Q2 absolute
	77	Double_t fQ4abs; // Q4 absolute
	78	Double_t fQ42re; // Re(Q2QQ)
	79	Double32_t fPtMax; //[0,0,16] maximum pT
	80	Double32_t fMPt; //[0,0,16] mean pT
	81	Double32_t fMPt2; //[0,0,16] mean pT2
	82	Double32_t fTS; //[0,0,16] transverse sphericity
	83	Double32_t fMV0M; // V0M amplitude
	84	Double32_t fCl1; //[0,0,16] class CL1
	85	Double32_t fV0M; //[0,0,16] class V0M
	86	Double32_t fV0MEq; //[0,0,16] class V0M Eq
	87	Double32_t fV0A; //[0,0,16] class V0A
	88	Double32_t fV0AEq; //[0,0,16] class V0A Eq
	89	Double32_t fZNA; //[0,0,16] class ZNA
	90
e29fb8cd	91	ClassDef(AliNtupCumInfo,2) // Cumulant storage class
96a87a98	92	};
96a87a98	93
d87c9c78	94	#endif