[u/mrichter/AliRoot.git] / MUON / AliMUONQADataMakerRec.h

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

// $Id$

/// \ingroup rec
/// \class AliMUONQADataMakerRec
/// \brief MUON Quality assurance data maker
///

// --- AliRoot header files ---
#include "AliQADataMakerRec.h"
#include "AliMUONRecoParam.h"

class AliMUONDigitMaker;
class AliMUONVDigitStore;
class AliMUONVTriggerStore;
class AliMUONVClusterStore;
class AliMUONVTrackerData;
class AliMUONVTrackerDataMaker;

class AliMUONQADataMakerRec: public AliQADataMakerRec {

public:
  AliMUONQADataMakerRec();         
  AliMUONQADataMakerRec(const AliMUONQADataMakerRec& qadm);   
  AliMUONQADataMakerRec& operator=(const AliMUONQADataMakerRec& qadm);
  virtual ~AliMUONQADataMakerRec();
  
  AliMUONVTrackerData* GetTrackerData() const;
  
protected:
	
  virtual void StartOfDetectorCycle(); 

  virtual void InitRaws(); 
  virtual void InitRecPoints(); 
  virtual void InitESDs(); 
  
  virtual void MakeRaws(AliRawReader* rawReader); 
  virtual void MakeRecPoints(TTree* recpo); 
  virtual void MakeESDs(AliESDEvent* esd) ;
  
  virtual void DefaultEndOfDetectorCycle(AliQA::TASKINDEX_t) {}

  virtual void EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray** list);

private:
  /// Raw histograms indices
  enum ERaw { 
    kTrackerData           = 3,  ///< Accumulated data
    kTrackerBusPatchOccupancy = 4, ///< Bus patch occupancies
    kTriggerScalersBP      = 22, ///< Trigger scalers on BP histogram per chamber index
    kTriggerScalersNBP     = 23, ///< Trigger scalers on NBP histogram per chamber index
    kTriggerScalersDisplay = 24  ///< Trigger scalers display histogram per chamber index
  };
         
  /// Rec points histograms indices
  enum ERecPoints { 
    kTriggerDigitsBendPlane    = 0,  ///< Trigger digits on BP histogram index
    kTriggerDigitsNonBendPlane = 1,  ///< Trigger digits on NBP histogram index
    kTriggeredBoards           = 2,  ///< Triggered boards histogram index
    kTriggerDigitsDisplay      = 3,  ///< Trigger digits display histogram per plane index
    kTriggerBoardsDisplay      = 11, ///< Triggered boards display histogram index
    kTriggerRPCi               = 12, ///< Trigger chamber currents index
    kTriggerRPChv              = 13, ///< Trigger chamber HV index
    kTriggerIDisplay           = 14, ///< Trigger chamber currents display histogram index
    kTriggerHVDisplay          = 18, ///< Trigger chamber HV display histogram index
    
    kTrackerNumberOfClustersPerChamber    = 100, ///< Tracker: # of clusters per chamber
    kTrackerClusterMultiplicityPerChamber = 200, ///< Tracker: cluster multiplicity per chamber
    kTrackerClusterChargePerChamber       = 300, ///< Tracker: cluster charge per chamber
				
    kTrackerNumberOfClustersPerDE    = 1000, ///< Tracker : number of clusters per DE		
    kTrackerClusterMultiplicityPerDE = 3000, ///< Tracker : cluster multiplicity per DE		
    kTrackerClusterChargePerDE       = 5000  ///< Tracker : cluster charge per DE
  };
  
  /// ESD histograms indices
  enum EESD { 
    kESDnTracks                 = 0,  ///< number of tracks
    kESDMatchTrig               = 1,  ///< number of tracks matched with trigger
    kESDMomentum                = 2,  ///< P distribution
    kESDPt                      = 3,  ///< Pt distribution
    kESDRapidity                = 4,  ///< rapidity distribution
    kESDChi2                    = 5,  ///< normalized chi2 distribution
    kESDProbChi2                = 6,  ///< distribution of probability of chi2
    
    kESDClusterHitMap           = 7,  ///< cluster position distribution in chamber i
    kESDnClustersPerTrack       = 17, ///< number of clusters per track
    kESDnClustersPerCh          = 18, ///< number of clusters per chamber per track
    kESDnClustersPerDE          = 19, ///< number of clusters per DE per track
    kESDClusterChargeInCh       = 20, ///< cluster charge distribution in chamber i
    kESDClusterChargePerChMean  = 30, ///< cluster charge per Ch: mean
    kESDClusterChargePerChSigma = 31, ///< cluster charge per Ch: dispersion
    kESDClusterChargePerDE      = 32, ///< cluster charge per DE: mean
    kESDClusterSizeInCh         = 33, ///< cluster size distribution in chamber i
    kESDClusterSizePerChMean    = 43, ///< cluster size per Ch: mean
    kESDClusterSizePerChSigma   = 44, ///< cluster size per Ch: dispersion
    kESDClusterSizePerDE        = 45, ///< cluster size per DE: mean
    
    kESDResidualXInCh           = 46, ///< cluster-track residual-X distribution in chamber i
    kESDResidualYInCh           = 56, ///< cluster-track residual-Y distribution in chamber i
    kESDResidualXPerChMean      = 66, ///< cluster-track residual-X per Ch: mean
    kESDResidualYPerChMean      = 67, ///< cluster-track residual-Y per Ch: mean
    kESDResidualXPerChSigma     = 68, ///< cluster-track residual-X per Ch: dispersion
    kESDResidualYPerChSigma     = 69, ///< cluster-track residual-Y per Ch: dispersion
    kESDResidualXPerDEMean      = 70, ///< cluster-track residual-X per DE: mean
    kESDResidualYPerDEMean      = 71, ///< cluster-track residual-Y per DE: mean
    kESDResidualXPerDESigma     = 72, ///< cluster-track residual-X per DE: dispersion
    kESDResidualYPerDESigma     = 73, ///< cluster-track residual-Y per DE: dispersion
    kESDLocalChi2XInCh          = 74, ///< local chi2-X distribution in chamber i
    kESDLocalChi2YInCh          = 84, ///< local chi2-Y distribution in chamber i
    kESDLocalChi2XPerChMean     = 94, ///< local chi2-X per Ch: mean
    kESDLocalChi2YPerChMean     = 95, ///< local chi2-Y per Ch: mean
    kESDLocalChi2XPerDEMean     = 96, ///< local chi2-X per DE: mean
    kESDLocalChi2YPerDEMean     = 97  ///< local chi2-Y per DE: mean
  };

private:
	
  void Ctor();
  void DisplayTriggerInfo(AliQA::TASKINDEX_t task);
  Bool_t FillTriggerDCSHistos();
  void InitRecPointsTracker();
  void InitRecPointsTrigger();
  void MakeRawsTracker(AliRawReader* rawReader);
  void MakeRawsTrigger(AliRawReader* rawReader);
  void MakeRecPointsTracker(TTree* treeR);
  void MakeRecPointsTrigger(TTree* treeR);
	
  /// Return reco parameters
  const AliMUONRecoParam* GetRecoParam() const { return dynamic_cast<const AliMUONRecoParam *>(fRecoParam); }
  
  Bool_t  fIsInitRaws;       //!<  info if InitRaws() went ok
  Bool_t  fIsInitRecPointsTracker;  //!<  info if InitRecPoints() went ok
  Bool_t  fIsInitRecPointsTrigger;  //!<  info if InitRecPoints() went ok
  Bool_t  fIsInitESDs;       //!<  info if InitESDs() went ok
  
  AliMUONVDigitStore*   fDigitStore; //!< pointer to digits store
  AliMUONVTriggerStore* fTriggerStore; //!< pointer to trigger store
  AliMUONDigitMaker*    fDigitMaker;  //!< pointer to digit maker
  AliMUONVClusterStore* fClusterStore; //!< pointer to cluster store
	
  AliMUONVTrackerDataMaker* fTrackerDataMaker; //!< tracker data accumulation
  
  TH1F* fhESDnTotClustersPerCh;      //!< total number of associated clusters per chamber
  TH1F* fhESDnTotClustersPerDE;      //!< total number of associated clusters per DE
  TH1F* fhESDnTotFullClustersPerDE;  //!< total number of associated clusters containing pad info per DE
  TH1F* fhESDSumClusterChargePerDE;  //!< sum of cluster charge per DE
  TH1F* fhESDSumClusterSizePerDE;    //!< sum of cluster size per DE
  TH1F* fhESDSumResidualXPerDE;      //!< sum of cluster-track residual-X per DE
  TH1F* fhESDSumResidualYPerDE;      //!< sum of cluster-track residual-Y per DE
  TH1F* fhESDSumResidualX2PerDE;     //!< sum of cluster-track residual-X**2 per DE
  TH1F* fhESDSumResidualY2PerDE;     //!< sum of cluster-track residual-Y**2 per DE
  TH1F* fhESDSumLocalChi2XPerDE;     //!< sum of local chi2-X per DE
  TH1F* fhESDSumLocalChi2YPerDE;     //!< sum of local chi2-X per DE
  
  ClassDef(AliMUONQADataMakerRec,5)  // MUON Quality assurance data maker

};
#endif
Commit	Line	Data
7d5d0cc5	1	#ifndef ALIMUONQADATAMAKERREC_H
7d5d0cc5	2	#define ALIMUONQADATAMAKERREC_H
04236e67	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
bf4d93eb	6	// $Id$
bf4d93eb	7
04236e67	8	/// \ingroup rec
	9	/// \class AliMUONQADataMakerRec
	10	/// \brief MUON Quality assurance data maker
	11	///
04236e67	12
04236e67	13	// --- AliRoot header files ---
c92a8c85	14	#include "AliQADataMakerRec.h"
	15	#include "AliMUONRecoParam.h"
	16
7035694f	17	class AliMUONDigitMaker;
	18	class AliMUONVDigitStore;
	19	class AliMUONVTriggerStore;
f587a77d	20	class AliMUONVClusterStore;
de1d4a53	21	class AliMUONVTrackerData;
f587a77d	22	class AliMUONVTrackerDataMaker;
f587a77d	23
04236e67	24	class AliMUONQADataMakerRec: public AliQADataMakerRec {
	25
	26	public:
b8bd1ab8	27	AliMUONQADataMakerRec();
04236e67	28	AliMUONQADataMakerRec(const AliMUONQADataMakerRec& qadm);
	29	AliMUONQADataMakerRec& operator=(const AliMUONQADataMakerRec& qadm);
	30	virtual ~AliMUONQADataMakerRec();
	31
de1d4a53	32	AliMUONVTrackerData* GetTrackerData() const;
de1d4a53	33
7d5d0cc5	34	protected:
	35
	36	virtual void StartOfDetectorCycle();
	37
	38	virtual void InitRaws();
	39	virtual void InitRecPoints();
	40	virtual void InitESDs();
	41
	42	virtual void MakeRaws(AliRawReader* rawReader);
	43	virtual void MakeRecPoints(TTree* recpo);
	44	virtual void MakeESDs(AliESDEvent* esd) ;
	45
	46	virtual void DefaultEndOfDetectorCycle(AliQA::TASKINDEX_t) {}
	47
21973f26	48	virtual void EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray** list);
7d5d0cc5	49
04236e67	50	private:
b62881f2	51	/// Raw histograms indices
b62881f2	52	enum ERaw {
c92a8c85	53	kTrackerData = 3, ///< Accumulated data
5bb54f70	54	kTrackerBusPatchOccupancy = 4, ///< Bus patch occupancies
c92a8c85	55	kTriggerScalersBP = 22, ///< Trigger scalers on BP histogram per chamber index
	56	kTriggerScalersNBP = 23, ///< Trigger scalers on NBP histogram per chamber index
	57	kTriggerScalersDisplay = 24 ///< Trigger scalers display histogram per chamber index
b62881f2	58	};
	59
	60	/// Rec points histograms indices
	61	enum ERecPoints {
c92a8c85	62	kTriggerDigitsBendPlane = 0, ///< Trigger digits on BP histogram index
2346e1ad	63	kTriggerDigitsNonBendPlane = 1, ///< Trigger digits on NBP histogram index
	64	kTriggeredBoards = 2, ///< Triggered boards histogram index
	65	kTriggerDigitsDisplay = 3, ///< Trigger digits display histogram per plane index
c92a8c85	66	kTriggerBoardsDisplay = 11, ///< Triggered boards display histogram index
49e110ec	67	kTriggerRPCi = 12, ///< Trigger chamber currents index
	68	kTriggerRPChv = 13, ///< Trigger chamber HV index
	69	kTriggerIDisplay = 14, ///< Trigger chamber currents display histogram index
	70	kTriggerHVDisplay = 18, ///< Trigger chamber HV display histogram index
c92a8c85	71
	72	kTrackerNumberOfClustersPerChamber = 100, ///< Tracker: # of clusters per chamber
	73	kTrackerClusterMultiplicityPerChamber = 200, ///< Tracker: cluster multiplicity per chamber
	74	kTrackerClusterChargePerChamber = 300, ///< Tracker: cluster charge per chamber
f587a77d	75
c92a8c85	76	kTrackerNumberOfClustersPerDE = 1000, ///< Tracker : number of clusters per DE
	77	kTrackerClusterMultiplicityPerDE = 3000, ///< Tracker : cluster multiplicity per DE
	78	kTrackerClusterChargePerDE = 5000 ///< Tracker : cluster charge per DE
b62881f2	79	};
c92a8c85	80
b62881f2	81	/// ESD histograms indices
b62881f2	82	enum EESD {
5a240757	83	kESDnTracks = 0, ///< number of tracks
	84	kESDMatchTrig = 1, ///< number of tracks matched with trigger
	85	kESDMomentum = 2, ///< P distribution
	86	kESDPt = 3, ///< Pt distribution
	87	kESDRapidity = 4, ///< rapidity distribution
	88	kESDChi2 = 5, ///< normalized chi2 distribution
	89	kESDProbChi2 = 6, ///< distribution of probability of chi2
c92a8c85	90
5a240757	91	kESDClusterHitMap = 7, ///< cluster position distribution in chamber i
	92	kESDnClustersPerTrack = 17, ///< number of clusters per track
	93	kESDnClustersPerCh = 18, ///< number of clusters per chamber per track
	94	kESDnClustersPerDE = 19, ///< number of clusters per DE per track
	95	kESDClusterChargeInCh = 20, ///< cluster charge distribution in chamber i
	96	kESDClusterChargePerChMean = 30, ///< cluster charge per Ch: mean
	97	kESDClusterChargePerChSigma = 31, ///< cluster charge per Ch: dispersion
	98	kESDClusterChargePerDE = 32, ///< cluster charge per DE: mean
	99	kESDClusterSizeInCh = 33, ///< cluster size distribution in chamber i
	100	kESDClusterSizePerChMean = 43, ///< cluster size per Ch: mean
	101	kESDClusterSizePerChSigma = 44, ///< cluster size per Ch: dispersion
	102	kESDClusterSizePerDE = 45, ///< cluster size per DE: mean
c92a8c85	103
5a240757	104	kESDResidualXInCh = 46, ///< cluster-track residual-X distribution in chamber i
	105	kESDResidualYInCh = 56, ///< cluster-track residual-Y distribution in chamber i
	106	kESDResidualXPerChMean = 66, ///< cluster-track residual-X per Ch: mean
	107	kESDResidualYPerChMean = 67, ///< cluster-track residual-Y per Ch: mean
	108	kESDResidualXPerChSigma = 68, ///< cluster-track residual-X per Ch: dispersion
	109	kESDResidualYPerChSigma = 69, ///< cluster-track residual-Y per Ch: dispersion
	110	kESDResidualXPerDEMean = 70, ///< cluster-track residual-X per DE: mean
	111	kESDResidualYPerDEMean = 71, ///< cluster-track residual-Y per DE: mean
	112	kESDResidualXPerDESigma = 72, ///< cluster-track residual-X per DE: dispersion
	113	kESDResidualYPerDESigma = 73, ///< cluster-track residual-Y per DE: dispersion
	114	kESDLocalChi2XInCh = 74, ///< local chi2-X distribution in chamber i
	115	kESDLocalChi2YInCh = 84, ///< local chi2-Y distribution in chamber i
	116	kESDLocalChi2XPerChMean = 94, ///< local chi2-X per Ch: mean
	117	kESDLocalChi2YPerChMean = 95, ///< local chi2-Y per Ch: mean
	118	kESDLocalChi2XPerDEMean = 96, ///< local chi2-X per DE: mean
	119	kESDLocalChi2YPerDEMean = 97 ///< local chi2-Y per DE: mean
b62881f2	120	};
b62881f2	121
f587a77d	122	private:
f587a77d	123
21973f26	124	void Ctor();
2346e1ad	125	void DisplayTriggerInfo(AliQA::TASKINDEX_t task);
49e110ec	126	Bool_t FillTriggerDCSHistos();
21973f26	127	void InitRecPointsTracker();
	128	void InitRecPointsTrigger();
	129	void MakeRawsTracker(AliRawReader* rawReader);
	130	void MakeRawsTrigger(AliRawReader* rawReader);
f587a77d	131	void MakeRecPointsTracker(TTree* treeR);
	132	void MakeRecPointsTrigger(TTree* treeR);
	133
21973f26	134	/// Return reco parameters
b8bd1ab8	135	const AliMUONRecoParam* GetRecoParam() const { return dynamic_cast<const AliMUONRecoParam *>(fRecoParam); }
a0dc65b4	136
b1341ab0	137	Bool_t fIsInitRaws; //!< info if InitRaws() went ok
f587a77d	138	Bool_t fIsInitRecPointsTracker; //!< info if InitRecPoints() went ok
f587a77d	139	Bool_t fIsInitRecPointsTrigger; //!< info if InitRecPoints() went ok
b1341ab0	140	Bool_t fIsInitESDs; //!< info if InitESDs() went ok
	141
	142	AliMUONVDigitStore* fDigitStore; //!< pointer to digits store
7035694f	143	AliMUONVTriggerStore* fTriggerStore; //!< pointer to trigger store
b1341ab0	144	AliMUONDigitMaker* fDigitMaker; //!< pointer to digit maker
f587a77d	145	AliMUONVClusterStore* fClusterStore; //!< pointer to cluster store
f587a77d	146
b8bd1ab8	147	AliMUONVTrackerDataMaker* fTrackerDataMaker; //!< tracker data accumulation
5a240757	148
	149	TH1F* fhESDnTotClustersPerCh; //!< total number of associated clusters per chamber
	150	TH1F* fhESDnTotClustersPerDE; //!< total number of associated clusters per DE
	151	TH1F* fhESDnTotFullClustersPerDE; //!< total number of associated clusters containing pad info per DE
	152	TH1F* fhESDSumClusterChargePerDE; //!< sum of cluster charge per DE
	153	TH1F* fhESDSumClusterSizePerDE; //!< sum of cluster size per DE
	154	TH1F* fhESDSumResidualXPerDE; //!< sum of cluster-track residual-X per DE
	155	TH1F* fhESDSumResidualYPerDE; //!< sum of cluster-track residual-Y per DE
	156	TH1F* fhESDSumResidualX2PerDE; //!< sum of cluster-track residual-X**2 per DE
	157	TH1F* fhESDSumResidualY2PerDE; //!< sum of cluster-track residual-Y**2 per DE
	158	TH1F* fhESDSumLocalChi2XPerDE; //!< sum of local chi2-X per DE
	159	TH1F* fhESDSumLocalChi2YPerDE; //!< sum of local chi2-X per DE
	160
	161	ClassDef(AliMUONQADataMakerRec,5) // MUON Quality assurance data maker
04236e67	162
	163	};
	164	#endif