3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
9 #include <TProfile2D.h>
10 #include "AliRecoParam.h"
15 class TTreeSRedirector;
18 class AliTPCRawStreamV3;
20 class AliTPCAltroMapping;
24 struct eventHeaderStruct;
26 class AliTPCdataQA : public TH1F {
30 AliTPCdataQA(const AliTPCdataQA &ped);
31 AliTPCdataQA(const TMap *config);
32 virtual ~AliTPCdataQA();
34 AliTPCdataQA& operator = (const AliTPCdataQA &source);
35 void DumpToFile(const Char_t *filename, const Char_t *dir="", Bool_t append=kFALSE);
36 void MakeTree(const char *fname="QApad.root") const;
39 Bool_t ProcessEvent(AliTPCRawStreamV3 *const rawStreamV3);
40 Bool_t ProcessEvent(AliRawReader *const rawReader);
41 Bool_t ProcessEvent(eventHeaderStruct *const event);
46 void SetPedestal(AliTPCCalPad *const pedestalCal){ fPedestal = pedestalCal;}
47 void SetNoise(AliTPCCalPad *const noiseCal){ fNoise = noiseCal;}
49 void SetMinQMax (Float_t minQmax ) { fMinQMax = minQmax; }
50 void SetRequireNeighbouringPad(Bool_t req=kTRUE) { fRequireNeighbouringPad = req; }
53 void FillOccupancyProfile();
57 AliTPCCalPad *GetNoThreshold() const { return fNoThreshold;}
58 AliTPCCalPad *GetMaxCharge() const { return fMaxCharge;}
59 AliTPCCalPad *GetMeanCharge() const { return fMeanCharge;}
60 AliTPCCalPad *GetNLocalMaxima() const { return fNLocalMaxima;}
61 AliTPCCalPad *GetOverThreshold10() const { return fOverThreshold10;}
62 AliTPCCalPad *GetOverThreshold20() const { return fOverThreshold20;}
63 AliTPCCalPad *GetOverThreshold30() const { return fOverThreshold30;}
64 AliTPCCalPad *GetNTimeBins() const { return fNTimeBins;}
65 AliTPCCalPad *GetNPads() const { return fNPads;}
66 AliTPCCalPad *GetTimePosition() const { return fTimePosition;}
67 TProfile* GetHistQVsTimeSideA() const {return fHistQVsTimeSideA;}
68 TProfile* GetHistQVsTimeSideC() const {return fHistQVsTimeSideC;}
69 TProfile* GetHistQMaxVsTimeSideA() const {return fHistQMaxVsTimeSideA;}
70 TProfile* GetHistQMaxVsTimeSideC() const {return fHistQMaxVsTimeSideC;}
71 TH1F* GetHistOccupancyVsEventConst() const {return fHistOccupancyVsEvent;}
72 TH1F* GetHistNclustersVsEventConst() const {return fHistNclustersVsEvent;}
73 TH1F* GetHistOccupancyVsEvent();
74 TH1F* GetHistNclustersVsEvent();
77 TProfile* GetHistOccVsSector() const { return fHistOccVsSector; }
78 TProfile2D* GetHistOcc2dVsSector() const { return fHistOcc2dVsSector; }
79 TProfile* GetHistQVsSector() const { return fHistQVsSector; }
80 TProfile* GetHistQmaxVsSector() const { return fHistQmaxVsSector; }
83 AliTPCAltroMapping **GetAltroMapping() const { return fMapping; };
84 void SetAltroMapping(AliTPCAltroMapping **mapp) { fMapping = mapp; };
87 Int_t GetFirstTimeBin() const { return fFirstTimeBin; }
88 Int_t GetLastTimeBin() const { return fLastTimeBin; }
89 Int_t GetAdcMin() const { return fAdcMin; }
90 Int_t GetAdcMax() const { return fAdcMax; }
91 Int_t GetEventCounter() const { return fEventCounter; }
92 Bool_t GetIsAnalysed() const { return fIsAnalysed; }
93 Int_t GetMaxEvents() const { return fMaxEvents; }
94 Int_t GetEventsPerBin() const { return fEventsPerBin; }
95 Int_t GetSignalCounter() const { return fSignalCounter; }
96 Int_t GetClusterCounter() const { return fClusterCounter;}
98 Float_t GetMinQMax() const { return fMinQMax; }
99 Bool_t GetRequireNeighbouringPad() const { return fRequireNeighbouringPad; }
102 Bool_t GetIsDQM() const { return fIsDQM; }
104 void SetRangeTime(Int_t tMin, Int_t tMax){ fFirstTimeBin=tMin; fLastTimeBin=tMax;} // Set time bin range that is used for the pedestal calibration
105 void SetRangeAdc (Int_t aMin, Int_t aMax){ fAdcMin=aMin; fAdcMax=aMax; } // Set adc range for the pedestal calibration
106 void SetMaxEvents (Int_t value) { fMaxEvents = value; }
107 void SetEventsPerBin(Int_t value) { fEventsPerBin = value; }
110 void SetIsDQM(Bool_t value) { fIsDQM = value; }
114 void SetChamberStatus(UInt_t roc, Bool_t status) { fActiveChambers.SetBitNumber(roc,status); }
115 Bool_t GetChamberStatus(UInt_t roc) {return fActiveChambers.TestBitNumber(roc);}
117 Int_t Update(const Int_t iSector, const Int_t iRow, const Int_t iPad,
118 const Int_t iTimeBin, Float_t signal,
119 const Int_t iPatch=-1, const Int_t iBranch=-1);
120 void FindLocalMaxima(const Int_t iSector);
122 void MakeArrays(); // Create arrays for random data acces
123 void CleanArrays(); // Clean arrays for random data acces
124 void SetExpandDigit(const Int_t iRow, Int_t iPad, Int_t iTimeBin,
125 const Float_t signal); // Fill arrays with signals
126 void GetPadAndTimeBin(Int_t bin, Int_t& iPad, Int_t& iTimeBin); // Get pad and time bin corresponding to the 1d bin
127 Float_t GetQ(const Float_t* adcArray, const Int_t time,
128 const Int_t pad, const Int_t maxTimeBins,
129 Int_t& timeMin,Int_t& timeMax,Int_t& padMin,Int_t& padMax) const;
130 void UpdateEventHistograms();
134 TObjArray *ConfigArrRocs(TObjArray *arr, const Text_t* name);
136 Int_t fFirstTimeBin; // First Time bin needed for analysis
137 Int_t fLastTimeBin; // Last Time bin needed for analysis
138 Int_t fAdcMin; // min adc channel of pedestal value
139 Int_t fAdcMax; // max adc channel of pedestal value
140 Float_t fMinQMax; // Minimun charge for Maximum ADC in cluster
141 Bool_t fRequireNeighbouringPad; // If clusterer should require a neighbouring pad to accept it
144 AliTPCAltroMapping **fMapping; //! Altro Mapping object
147 AliTPCCalPad * fPedestal; //! option to set pedestal cal object
148 AliTPCCalPad * fNoise; //! option to set noise cal object
149 AliTPCCalPad * fNLocalMaxima; // local maximas found
150 AliTPCCalPad * fMaxCharge; // max charge
151 AliTPCCalPad * fMeanCharge; // mean charge
152 AliTPCCalPad * fNoThreshold; // number of digits
153 AliTPCCalPad * fNTimeBins; // timebins width of cluster
154 AliTPCCalPad * fNPads; // pads with of cluster
155 AliTPCCalPad * fTimePosition; // Time position of local maximum
156 AliTPCCalPad * fOverThreshold10; //! local maxima with qMax over threshold
157 AliTPCCalPad * fOverThreshold20; //! local maxima with qMax over threshold
158 AliTPCCalPad * fOverThreshold30; //! local maxima with qMax over threshold
160 TProfile* fHistQVsTimeSideA; // Q vs time (side A)
161 TProfile* fHistQVsTimeSideC; // Q vs time (side C)
162 TProfile* fHistQMaxVsTimeSideA; // QMax vs time (side A)
163 TProfile* fHistQMaxVsTimeSideC; // QMax vs time (side C)
165 TH1F* fHistOccupancyVsEvent; // Occupancy vs event number (~time)
166 TH1F* fHistNclustersVsEvent; // Nclusters vs event number (~time)
168 Int_t fEventCounter; // event Counter
169 Bool_t fIsAnalysed; // Set to true after Analyse has been called
171 Int_t fMaxEvents; // Max events for event histograms
172 Int_t fEventsPerBin; // Events per bin for event histograms
173 Int_t fSignalCounter; // Signal counter
174 Int_t fClusterCounter; // Cluster counter
176 TBits fActiveChambers; // configured ROCs
181 Float_t** fAllBins; //! array for digit using random access
182 Int_t** fAllSigBins; //! array of pointers to the indexes over threshold
183 Int_t* fAllNSigBins; //!
184 Int_t fRowsMax; //! Maximum number of time bins
185 Int_t fPadsMax; //! Maximum number of time bins
186 Int_t fTimeBinsMax; //! Maximum number of time bins
189 Bool_t fIsDQM; //! Is DQM -> Simple output (no 2D!)
190 TProfile* fHistOccVsSector; //! Occ vs sector (for DQM only)
191 TProfile2D* fHistOcc2dVsSector; //! Occ vs sector 2D (for DQM only)
192 TProfile* fHistQVsSector; //! Q vs sector (for DQM only)
193 TProfile* fHistQmaxVsSector; //! QMax vs sector (for DQM only)
194 TArrayD* fOccVec; //! Occupancy help counter for DQM
195 TArrayD* fOccMaxVec; //! Occupancy help normlization for DQM
196 TArrayD* fOccVecFine; //! "2D" occupancy help counter for DQM
197 TArrayD* fOccMaxVecFine; //! "2D" occupancy help normlization for DQM
200 ClassDef(AliTPCdataQA, 6) // Implementation of the TPC Raw QA