[u/mrichter/AliRoot.git] / PWG3 / muon / AliAnalysisTaskMuonQA.h

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

/* $Id$ */ 

/// \ingroup muondep
/// \class AliAnalysisTaskMuonQA
/// \brief Quality assurance of MUON ESDs
//Author: Philippe Pillot - SUBATECH Nantes

class TMap;
class TList;
class TObjArray;
class AliCounterCollection;

class AliAnalysisTaskMuonQA : public AliAnalysisTaskSE {
 public:
  
  AliAnalysisTaskMuonQA();
  AliAnalysisTaskMuonQA(const char *name);
  virtual ~AliAnalysisTaskMuonQA();
  
  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *);
  virtual void   Terminate(Option_t *);
  
  /// Select negative (<0), positive (>0) or all (==0) tracks to fill histograms
  void SelectCharge(Short_t charge = 0) {fSelectCharge = charge;}
  
  /// Select events passing the physics selection to fill histograms
  void SelectPhysics(Bool_t flag = kTRUE) {fSelectPhysics = flag;}
  
  /// Select events belonging to at least one of the trigger classes selected by the mask to fill histograms:
  /// - if the physics selection is used, apply the mask to the trigger word returned by the physics selection
  /// - if not, apply the mask to the trigger word built by looking for triggers listed in "fSelectTriggerClass"
  void SelectTrigger(Bool_t flag = kTRUE, UInt_t mask = AliVEvent::kMUON) {fSelectTrigger = flag; fTriggerMask = mask;}
	
  /// Select track matching the trigger to fill histograms
  void SelectMatched(Bool_t flag = kTRUE) {fSelectMatched = flag;}
  
  /// Use only tracks passing the acceptance cuts (Rabs, eta)
  void ApplyAccCut(Bool_t flag = kTRUE) { fApplyAccCut = flag; }
  
 private:
  
  /// Not implemented
  AliAnalysisTaskMuonQA(const AliAnalysisTaskMuonQA& rhs);
  /// Not implemented
  AliAnalysisTaskMuonQA& operator = (const AliAnalysisTaskMuonQA& rhs);
  
  Double_t ChangeThetaRange(Double_t theta);
  
  UInt_t BuildTriggerWord(TString& FiredTriggerClasses);
  
  TList* BuildListOfTriggerCases(TString& FiredTriggerClasses);
  TList* BuildListOfAllTriggerCases(TString& FiredTriggerClasses);
  TList* BuildListOfSelectedTriggerCases(TString& FiredTriggerClasses);
	
 private:
  
  enum eList {
    kNTracks                 = 0,  ///< number of tracks
    kMatchTrig               = 1,  ///< number of tracks matched with trigger
    kSign                    = 2,  ///< track sign
    kDCA                     = 3,  ///< DCA distribution
    kP                       = 4,  ///< P distribution
    kPMuPlus                 = 5,  ///< P distribution of mu+
    kPMuMinus                = 6,  ///< P distribution of mu-
    kPt                      = 7,  ///< Pt distribution
    kPtMuPlus                = 8,  ///< Pt distribution of mu+
    kPtMuMinus               = 9,  ///< Pt distribution of mu-
    kRapidity                = 10, ///< rapidity distribution
    kThetaX                  = 11, ///< thetaX distribution
    kThetaY                  = 12, ///< thetaY distribution
    kChi2                    = 13, ///< normalized chi2 distribution
    kProbChi2                = 14, ///< distribution of probability of chi2
    kNClustersPerTrack       = 15, ///< number of clusters per track
    kNChamberHitPerTrack     = 16,  ///< number of chamber hit per track
    kPtMatchLpt              = 17, ///< Pt distribution match Lpt
    kPtMatchHpt              = 18, ///< Pt distribution match Hpt
    kPtMuPlusMatchLpt        = 19,  ///< Pt distribution of mu+ match Lpt
    kPtMuPlusMatchHpt        = 20,  ///< Pt distribution of mu+ match Hpt
    kPtMuMinusMatchLpt       = 21,  ///< Pt distribution of mu- match Lpt
    kPtMuMinusMatchHpt       = 22   ///< Pt distribution of mu- match Hpt
  };
  
  enum eListExpert {
    kNClustersPerCh          = 0,  ///< number of clusters per chamber
    kNClustersPerDE          = 1,  ///< number of clusters per DE
    kClusterHitMapInCh       = 2,  ///< cluster position distribution in chamber i
    kClusterChargeInCh       = 12, ///< cluster charge distribution in chamber i
    kClusterChargePerDE      = 22, ///< cluster charge distribution per DE
    kClusterSizeInCh         = 23, ///< cluster size distribution in chamber i
    kClusterSizePerDE        = 33  ///< cluster size distribution per DE
  };
  
  enum eListNorm {
    kClusterChargePerChMean  = 0,  ///< cluster charge per Ch: mean
    kClusterChargePerChSigma = 1,  ///< cluster charge per Ch: dispersion
    kClusterChargePerDEMean  = 2,  ///< cluster charge per DE: mean
    kClusterChargePerDESigma = 3,  ///< cluster charge per DE: dispersion
    kClusterSizePerChMean    = 4,  ///< cluster size per Ch: mean
    kClusterSizePerChSigma   = 5,  ///< cluster size per Ch: dispersion
    kClusterSizePerDEMean    = 6,  ///< cluster size per DE: mean
    kClusterSizePerDESigma   = 7,  ///< cluster size per DE: dispersion
    kNClustersPerChPerTrack  = 8,  ///< number of clusters per chamber per track
    kNClustersPerDEPerTrack  = 9   ///< number of clusters per DE per track
  };
  
  TObjArray*  fList;       //!< List of output object for everybody
  TObjArray*  fListExpert; //!< List of output object for experts
  TObjArray*  fListNorm;   //!< Normalized histograms
  
  AliCounterCollection* fTrackCounters; //!< track statistics
  AliCounterCollection* fEventCounters; //!< event statistics
  
  Short_t fSelectCharge;  ///< Fill histograms only with negative/position tracks (0=all)
  Bool_t  fSelectPhysics; ///< Fill histograms only with events passing the physics selection
  Bool_t  fSelectTrigger; ///< Fill histograms only with events passing the trigger selection
  UInt_t  fTriggerMask;   ///< Trigger mask to be used when selecting events
  Bool_t  fSelectMatched; ///< Fill histograms only with tracks matching the trigger
  Bool_t  fApplyAccCut;   ///< use only tracks passing the acceptance cuts (Rabs, eta)
  
  TMap*  fTriggerClass;       //!< map of trigger class name associated to short name
  TList* fSelectTriggerClass; //!< list of trigger class that can be selected to fill histograms
  
  static const Int_t nCh;       ///< number of tracking chambers
  static const Int_t nDE;       ///< number of DE
  static const Float_t dMax[5]; ///< maximum diameter of each station
  
  ClassDef(AliAnalysisTaskMuonQA, 6);
};

#endif
Commit	Line	Data
94ef1a28	1	#ifndef ALIANALYSISTASKMUONQA_H
	2	#define ALIANALYSISTASKMUONQA_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
27de2dfb	6	/* $Id$ */
27de2dfb	7
94ef1a28	8	/// \ingroup muondep
	9	/// \class AliAnalysisTaskMuonQA
	10	/// \brief Quality assurance of MUON ESDs
	11	//Author: Philippe Pillot - SUBATECH Nantes
	12
03ac5989	13	class TMap;
03ac5989	14	class TList;
94ef1a28	15	class TObjArray;
	16	class AliCounterCollection;
	17
	18	class AliAnalysisTaskMuonQA : public AliAnalysisTaskSE {
16dc5a7c	19	public:
94ef1a28	20
40be9f05	21	AliAnalysisTaskMuonQA();
40be9f05	22	AliAnalysisTaskMuonQA(const char *name);
94ef1a28	23	virtual ~AliAnalysisTaskMuonQA();
	24
	25	virtual void UserCreateOutputObjects();
	26	virtual void UserExec(Option_t *);
	27	virtual void Terminate(Option_t *);
	28
b4f7418a	29	/// Select negative (<0), positive (>0) or all (==0) tracks to fill histograms
94ef1a28	30	void SelectCharge(Short_t charge = 0) {fSelectCharge = charge;}
b4f7418a	31
b4f7418a	32	/// Select events passing the physics selection to fill histograms
94ef1a28	33	void SelectPhysics(Bool_t flag = kTRUE) {fSelectPhysics = flag;}
94ef1a28	34
03ac5989	35	/// Select events belonging to at least one of the trigger classes selected by the mask to fill histograms:
	36	/// - if the physics selection is used, apply the mask to the trigger word returned by the physics selection
	37	/// - if not, apply the mask to the trigger word built by looking for triggers listed in "fSelectTriggerClass"
	38	void SelectTrigger(Bool_t flag = kTRUE, UInt_t mask = AliVEvent::kMUON) {fSelectTrigger = flag; fTriggerMask = mask;}
16dc5a7c	39
76e3d5d4	40	/// Select track matching the trigger to fill histograms
	41	void SelectMatched(Bool_t flag = kTRUE) {fSelectMatched = flag;}
	42
417d6f63	43	/// Use only tracks passing the acceptance cuts (Rabs, eta)
	44	void ApplyAccCut(Bool_t flag = kTRUE) { fApplyAccCut = flag; }
	45
16dc5a7c	46	private:
94ef1a28	47
	48	/// Not implemented
	49	AliAnalysisTaskMuonQA(const AliAnalysisTaskMuonQA& rhs);
	50	/// Not implemented
	51	AliAnalysisTaskMuonQA& operator = (const AliAnalysisTaskMuonQA& rhs);
	52
	53	Double_t ChangeThetaRange(Double_t theta);
	54
03ac5989	55	UInt_t BuildTriggerWord(TString& FiredTriggerClasses);
	56
	57	TList* BuildListOfTriggerCases(TString& FiredTriggerClasses);
16dc5a7c	58	TList* BuildListOfAllTriggerCases(TString& FiredTriggerClasses);
	59	TList* BuildListOfSelectedTriggerCases(TString& FiredTriggerClasses);
	60
	61	private:
94ef1a28	62
03ac5989	63	enum eList {
94ef1a28	64	kNTracks = 0, ///< number of tracks
	65	kMatchTrig = 1, ///< number of tracks matched with trigger
	66	kSign = 2, ///< track sign
	67	kDCA = 3, ///< DCA distribution
	68	kP = 4, ///< P distribution
03ac5989	69	kPMuPlus = 5, ///< P distribution of mu+
	70	kPMuMinus = 6, ///< P distribution of mu-
	71	kPt = 7, ///< Pt distribution
	72	kPtMuPlus = 8, ///< Pt distribution of mu+
	73	kPtMuMinus = 9, ///< Pt distribution of mu-
	74	kRapidity = 10, ///< rapidity distribution
	75	kThetaX = 11, ///< thetaX distribution
	76	kThetaY = 12, ///< thetaY distribution
	77	kChi2 = 13, ///< normalized chi2 distribution
	78	kProbChi2 = 14, ///< distribution of probability of chi2
	79	kNClustersPerTrack = 15, ///< number of clusters per track
16dc5a7c	80	kNChamberHitPerTrack = 16, ///< number of chamber hit per track
b8e17107	81	kPtMatchLpt = 17, ///< Pt distribution match Lpt
	82	kPtMatchHpt = 18, ///< Pt distribution match Hpt
	83	kPtMuPlusMatchLpt = 19, ///< Pt distribution of mu+ match Lpt
	84	kPtMuPlusMatchHpt = 20, ///< Pt distribution of mu+ match Hpt
	85	kPtMuMinusMatchLpt = 21, ///< Pt distribution of mu- match Lpt
	86	kPtMuMinusMatchHpt = 22 ///< Pt distribution of mu- match Hpt
03ac5989	87	};
	88
	89	enum eListExpert {
	90	kNClustersPerCh = 0, ///< number of clusters per chamber
	91	kNClustersPerDE = 1, ///< number of clusters per DE
	92	kClusterHitMapInCh = 2, ///< cluster position distribution in chamber i
	93	kClusterChargeInCh = 12, ///< cluster charge distribution in chamber i
	94	kClusterChargePerDE = 22, ///< cluster charge distribution per DE
	95	kClusterSizeInCh = 23, ///< cluster size distribution in chamber i
	96	kClusterSizePerDE = 33 ///< cluster size distribution per DE
	97	};
	98
	99	enum eListNorm {
	100	kClusterChargePerChMean = 0, ///< cluster charge per Ch: mean
	101	kClusterChargePerChSigma = 1, ///< cluster charge per Ch: dispersion
	102	kClusterChargePerDEMean = 2, ///< cluster charge per DE: mean
	103	kClusterChargePerDESigma = 3, ///< cluster charge per DE: dispersion
	104	kClusterSizePerChMean = 4, ///< cluster size per Ch: mean
	105	kClusterSizePerChSigma = 5, ///< cluster size per Ch: dispersion
	106	kClusterSizePerDEMean = 6, ///< cluster size per DE: mean
	107	kClusterSizePerDESigma = 7, ///< cluster size per DE: dispersion
	108	kNClustersPerChPerTrack = 8, ///< number of clusters per chamber per track
	109	kNClustersPerDEPerTrack = 9 ///< number of clusters per DE per track
94ef1a28	110	};
	111
	112	TObjArray* fList; //!< List of output object for everybody
	113	TObjArray* fListExpert; //!< List of output object for experts
40be9f05	114	TObjArray* fListNorm; //!< Normalized histograms
94ef1a28	115
	116	AliCounterCollection* fTrackCounters; //!< track statistics
	117	AliCounterCollection* fEventCounters; //!< event statistics
	118
	119	Short_t fSelectCharge; ///< Fill histograms only with negative/position tracks (0=all)
76e3d5d4	120	Bool_t fSelectPhysics; ///< Fill histograms only with events passing the physics selection
	121	Bool_t fSelectTrigger; ///< Fill histograms only with events passing the trigger selection
	122	UInt_t fTriggerMask; ///< Trigger mask to be used when selecting events
	123	Bool_t fSelectMatched; ///< Fill histograms only with tracks matching the trigger
417d6f63	124	Bool_t fApplyAccCut; ///< use only tracks passing the acceptance cuts (Rabs, eta)
03ac5989	125
	126	TMap* fTriggerClass; //!< map of trigger class name associated to short name
	127	TList* fSelectTriggerClass; //!< list of trigger class that can be selected to fill histograms
94ef1a28	128
	129	static const Int_t nCh; ///< number of tracking chambers
	130	static const Int_t nDE; ///< number of DE
	131	static const Float_t dMax[5]; ///< maximum diameter of each station
94ef1a28	132
417d6f63	133	ClassDef(AliAnalysisTaskMuonQA, 6);
94ef1a28	134	};
	135
	136	#endif
	137