[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
    kTriggerScalersTime    = 32  ///< Trigger scalers acquisition time 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
    
    kESDnTotClustersPerCh       = 1000, ///< total number of associated clusters per chamber
    kESDnTotClustersPerDE       = 1001, ///< total number of associated clusters per DE
    kESDnTotFullClustersPerDE   = 1002, ///< total number of associated clusters containing pad info per DE
    kESDSumClusterChargePerDE   = 1003, ///< sum of cluster charge per DE
    kESDSumClusterSizePerDE     = 1004, ///< sum of cluster size per DE
    kESDSumResidualXPerDE       = 1005, ///< sum of cluster-track residual-X per DE
    kESDSumResidualYPerDE       = 1006, ///< sum of cluster-track residual-Y per DE
    kESDSumResidualX2PerDE      = 1007, ///< sum of cluster-track residual-X**2 per DE
    kESDSumResidualY2PerDE      = 1008, ///< sum of cluster-track residual-Y**2 per DE
    kESDSumLocalChi2XPerDE      = 1009, ///< sum of local chi2-X per DE
    kESDSumLocalChi2YPerDE      = 1010  ///< sum of local chi2-X per DE
  };
  
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
  
  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
c92a8c85	56	kTriggerScalersNBP = 23, ///< Trigger scalers on NBP histogram per chamber index
eb3886c8	57	kTriggerScalersDisplay = 24, ///< Trigger scalers display histogram per chamber index
eb3886c8	58	kTriggerScalersTime = 32 ///< Trigger scalers acquisition time index
b62881f2	59	};
	60
	61	/// Rec points histograms indices
	62	enum ERecPoints {
c92a8c85	63	kTriggerDigitsBendPlane = 0, ///< Trigger digits on BP histogram index
2346e1ad	64	kTriggerDigitsNonBendPlane = 1, ///< Trigger digits on NBP histogram index
	65	kTriggeredBoards = 2, ///< Triggered boards histogram index
	66	kTriggerDigitsDisplay = 3, ///< Trigger digits display histogram per plane index
c92a8c85	67	kTriggerBoardsDisplay = 11, ///< Triggered boards display histogram index
49e110ec	68	kTriggerRPCi = 12, ///< Trigger chamber currents index
	69	kTriggerRPChv = 13, ///< Trigger chamber HV index
	70	kTriggerIDisplay = 14, ///< Trigger chamber currents display histogram index
	71	kTriggerHVDisplay = 18, ///< Trigger chamber HV display histogram index
c92a8c85	72
	73	kTrackerNumberOfClustersPerChamber = 100, ///< Tracker: # of clusters per chamber
	74	kTrackerClusterMultiplicityPerChamber = 200, ///< Tracker: cluster multiplicity per chamber
	75	kTrackerClusterChargePerChamber = 300, ///< Tracker: cluster charge per chamber
f587a77d	76
c92a8c85	77	kTrackerNumberOfClustersPerDE = 1000, ///< Tracker : number of clusters per DE
	78	kTrackerClusterMultiplicityPerDE = 3000, ///< Tracker : cluster multiplicity per DE
	79	kTrackerClusterChargePerDE = 5000 ///< Tracker : cluster charge per DE
b62881f2	80	};
c92a8c85	81
b62881f2	82	/// ESD histograms indices
b62881f2	83	enum EESD {
5a240757	84	kESDnTracks = 0, ///< number of tracks
	85	kESDMatchTrig = 1, ///< number of tracks matched with trigger
	86	kESDMomentum = 2, ///< P distribution
	87	kESDPt = 3, ///< Pt distribution
	88	kESDRapidity = 4, ///< rapidity distribution
	89	kESDChi2 = 5, ///< normalized chi2 distribution
	90	kESDProbChi2 = 6, ///< distribution of probability of chi2
c92a8c85	91
5a240757	92	kESDClusterHitMap = 7, ///< cluster position distribution in chamber i
	93	kESDnClustersPerTrack = 17, ///< number of clusters per track
	94	kESDnClustersPerCh = 18, ///< number of clusters per chamber per track
	95	kESDnClustersPerDE = 19, ///< number of clusters per DE per track
	96	kESDClusterChargeInCh = 20, ///< cluster charge distribution in chamber i
	97	kESDClusterChargePerChMean = 30, ///< cluster charge per Ch: mean
	98	kESDClusterChargePerChSigma = 31, ///< cluster charge per Ch: dispersion
	99	kESDClusterChargePerDE = 32, ///< cluster charge per DE: mean
	100	kESDClusterSizeInCh = 33, ///< cluster size distribution in chamber i
	101	kESDClusterSizePerChMean = 43, ///< cluster size per Ch: mean
	102	kESDClusterSizePerChSigma = 44, ///< cluster size per Ch: dispersion
	103	kESDClusterSizePerDE = 45, ///< cluster size per DE: mean
c92a8c85	104
5a240757	105	kESDResidualXInCh = 46, ///< cluster-track residual-X distribution in chamber i
	106	kESDResidualYInCh = 56, ///< cluster-track residual-Y distribution in chamber i
	107	kESDResidualXPerChMean = 66, ///< cluster-track residual-X per Ch: mean
	108	kESDResidualYPerChMean = 67, ///< cluster-track residual-Y per Ch: mean
	109	kESDResidualXPerChSigma = 68, ///< cluster-track residual-X per Ch: dispersion
	110	kESDResidualYPerChSigma = 69, ///< cluster-track residual-Y per Ch: dispersion
	111	kESDResidualXPerDEMean = 70, ///< cluster-track residual-X per DE: mean
	112	kESDResidualYPerDEMean = 71, ///< cluster-track residual-Y per DE: mean
	113	kESDResidualXPerDESigma = 72, ///< cluster-track residual-X per DE: dispersion
	114	kESDResidualYPerDESigma = 73, ///< cluster-track residual-Y per DE: dispersion
	115	kESDLocalChi2XInCh = 74, ///< local chi2-X distribution in chamber i
	116	kESDLocalChi2YInCh = 84, ///< local chi2-Y distribution in chamber i
	117	kESDLocalChi2XPerChMean = 94, ///< local chi2-X per Ch: mean
	118	kESDLocalChi2YPerChMean = 95, ///< local chi2-Y per Ch: mean
	119	kESDLocalChi2XPerDEMean = 96, ///< local chi2-X per DE: mean
abcfd6d0	120	kESDLocalChi2YPerDEMean = 97, ///< local chi2-Y per DE: mean
	121
	122	kESDnTotClustersPerCh = 1000, ///< total number of associated clusters per chamber
	123	kESDnTotClustersPerDE = 1001, ///< total number of associated clusters per DE
	124	kESDnTotFullClustersPerDE = 1002, ///< total number of associated clusters containing pad info per DE
	125	kESDSumClusterChargePerDE = 1003, ///< sum of cluster charge per DE
	126	kESDSumClusterSizePerDE = 1004, ///< sum of cluster size per DE
	127	kESDSumResidualXPerDE = 1005, ///< sum of cluster-track residual-X per DE
	128	kESDSumResidualYPerDE = 1006, ///< sum of cluster-track residual-Y per DE
	129	kESDSumResidualX2PerDE = 1007, ///< sum of cluster-track residual-X**2 per DE
	130	kESDSumResidualY2PerDE = 1008, ///< sum of cluster-track residual-Y**2 per DE
	131	kESDSumLocalChi2XPerDE = 1009, ///< sum of local chi2-X per DE
	132	kESDSumLocalChi2YPerDE = 1010 ///< sum of local chi2-X per DE
b62881f2	133	};
abcfd6d0	134
f587a77d	135	private:
f587a77d	136
21973f26	137	void Ctor();
2346e1ad	138	void DisplayTriggerInfo(AliQA::TASKINDEX_t task);
49e110ec	139	Bool_t FillTriggerDCSHistos();
21973f26	140	void InitRecPointsTracker();
	141	void InitRecPointsTrigger();
	142	void MakeRawsTracker(AliRawReader* rawReader);
	143	void MakeRawsTrigger(AliRawReader* rawReader);
f587a77d	144	void MakeRecPointsTracker(TTree* treeR);
	145	void MakeRecPointsTrigger(TTree* treeR);
	146
21973f26	147	/// Return reco parameters
b8bd1ab8	148	const AliMUONRecoParam* GetRecoParam() const { return dynamic_cast<const AliMUONRecoParam *>(fRecoParam); }
a0dc65b4	149
b1341ab0	150	Bool_t fIsInitRaws; //!< info if InitRaws() went ok
f587a77d	151	Bool_t fIsInitRecPointsTracker; //!< info if InitRecPoints() went ok
f587a77d	152	Bool_t fIsInitRecPointsTrigger; //!< info if InitRecPoints() went ok
b1341ab0	153	Bool_t fIsInitESDs; //!< info if InitESDs() went ok
	154
	155	AliMUONVDigitStore* fDigitStore; //!< pointer to digits store
7035694f	156	AliMUONVTriggerStore* fTriggerStore; //!< pointer to trigger store
b1341ab0	157	AliMUONDigitMaker* fDigitMaker; //!< pointer to digit maker
f587a77d	158	AliMUONVClusterStore* fClusterStore; //!< pointer to cluster store
f587a77d	159
b8bd1ab8	160	AliMUONVTrackerDataMaker* fTrackerDataMaker; //!< tracker data accumulation
5a240757	161
5a240757	162	ClassDef(AliMUONQADataMakerRec,5) // MUON Quality assurance data maker
04236e67	163
	164	};
	165	#endif