#include <TH1F.h>
-#include <TObjArray.h>
+#include <TProfile.h>
+#include <TProfile2D.h>
+#include "AliRecoParam.h"
+
+#include <TArray.h>
-class TArrayF;
class TH2F;
class TTreeSRedirector;
class AliTPCROC;
class AliTPCCalROC;
class AliTPCRawStream;
-class AliTPCRawStreamFast;
+class AliTPCRawStreamV3;
class AliRawReader;
class AliTPCAltroMapping;
class AliTPCCalPad;
public:
AliTPCdataQA();
AliTPCdataQA(const AliTPCdataQA &ped);
- AliTPCdataQA(TMap *config);
+ AliTPCdataQA(const TMap *config);
virtual ~AliTPCdataQA();
AliTPCdataQA& operator = (const AliTPCdataQA &source);
- void DumpToFile(const Char_t *filename, const Char_t *dir="", const Bool_t append=kFALSE);
- void MakeTree(const char *fname="QApad.root");
+ void DumpToFile(const Char_t *filename, const Char_t *dir="", Bool_t append=kFALSE);
+ void MakeTree(const char *fname="QApad.root") const;
//
- Bool_t ProcessEventFast(AliTPCRawStreamFast *rawStreamFast);
- Bool_t ProcessEventFast(AliRawReader *rawReader);
- Bool_t ProcessEvent(AliTPCRawStream *rawStream);
- Bool_t ProcessEvent(AliRawReader *rawReader);
- Bool_t ProcessEvent(eventHeaderStruct *event);
+ Bool_t ProcessEvent(AliTPCRawStream *const rawStream);
+ Bool_t ProcessEvent(AliTPCRawStreamV3 *const rawStreamV3);
+ Bool_t ProcessEvent(AliRawReader *const rawReader);
+ Bool_t ProcessEventOld(AliRawReader *const rawReader);
+ Bool_t ProcessEvent(eventHeaderStruct *const event);
void Analyse();
//
//
- void SetPedestal(AliTPCCalPad *pedestalCal){ fPedestal = pedestalCal;}
- void SetNoise(AliTPCCalPad *noiseCal){ fNoise = noiseCal;}
+ void SetPedestal(AliTPCCalPad *const pedestalCal){ fPedestal = pedestalCal;}
+ void SetNoise(AliTPCCalPad *const noiseCal){ fNoise = noiseCal;}
+
+ // DQM methods
+ void FillOccupancyProfile();
+ void ResetProfiles();
+
AliTPCCalPad *GetNoThreshold() const { return fNoThreshold;}
AliTPCCalPad *GetMaxCharge() const { return fMaxCharge;}
AliTPCCalPad *GetNTimeBins() const { return fNTimeBins;}
AliTPCCalPad *GetNPads() const { return fNPads;}
AliTPCCalPad *GetTimePosition() const { return fTimePosition;}
+ TProfile* GetHistQVsTimeSideA() const {return fHistQVsTimeSideA;}
+ TProfile* GetHistQVsTimeSideC() const {return fHistQVsTimeSideC;}
+ TProfile* GetHistQMaxVsTimeSideA() const {return fHistQMaxVsTimeSideA;}
+ TProfile* GetHistQMaxVsTimeSideC() const {return fHistQMaxVsTimeSideC;}
+ TH1F* GetHistOccupancyVsEventConst() const {return fHistOccupancyVsEvent;}
+ TH1F* GetHistNclustersVsEventConst() const {return fHistNclustersVsEvent;}
+ TH1F* GetHistOccupancyVsEvent();
+ TH1F* GetHistNclustersVsEvent();
+
+ // DQM output
+ TProfile* GetHistOccVsSector() const { return fHistOccVsSector; }
+ TProfile2D* GetHistOcc2dVsSector() const { return fHistOcc2dVsSector; }
+ TProfile* GetHistQVsSector() const { return fHistQVsSector; }
+ TProfile* GetHistQmaxVsSector() const { return fHistQmaxVsSector; }
//
- AliTPCAltroMapping **GetAltroMapping() { return fMapping; };
+ AliTPCAltroMapping **GetAltroMapping() const { return fMapping; };
void SetAltroMapping(AliTPCAltroMapping **mapp) { fMapping = mapp; };
//
//
Int_t GetAdcMax() const { return fAdcMax; }
Int_t GetEventCounter() const { return fEventCounter; }
Bool_t GetIsAnalysed() const { return fIsAnalysed; }
- void SetRangeTime(Int_t tMin, Int_t tMax){ fFirstTimeBin=tMin; fLastTimeBin=tMax;} // Set time bin range that is used for the pedestal calibration
- void SetRangeAdc (Int_t aMin, Int_t aMax){ fAdcMin=aMin; fAdcMax=aMax; } // Set adc range for the pedestal calibration
+ Int_t GetMaxEvents() const { return fMaxEvents; }
+ Int_t GetEventsPerBin() const { return fEventsPerBin; }
+ Int_t GetSignalCounter() const { return fSignalCounter; }
+ Int_t GetClusterCounter() const { return fClusterCounter;}
+ // DQM getter
+ Bool_t GetIsDQM() const { return fIsDQM; }
+ void SetRangeTime(Int_t tMin, Int_t tMax){ fFirstTimeBin=tMin; fLastTimeBin=tMax;} // Set time bin range that is used for the pedestal calibration
+ void SetRangeAdc (Int_t aMin, Int_t aMax){ fAdcMin=aMin; fAdcMax=aMax; } // Set adc range for the pedestal calibration
+ void SetMaxEvents (Int_t value) { fMaxEvents = value; }
+ void SetEventsPerBin(Int_t value) { fEventsPerBin = value; }
+
+ // DQM setter
+ void SetIsDQM(Bool_t value) { fIsDQM = value; }
private:
Int_t Update(const Int_t iSector, const Int_t iRow, const Int_t iPad,
- const Int_t iTimeBin, Float_t signal);
+ const Int_t iTimeBin, Float_t signal,
+ const Int_t iPatch=-1, const Int_t iBranch=-1);
void FindLocalMaxima(const Int_t iSector);
void MakeArrays(); // Create arrays for random data acces
void GetPadAndTimeBin(Int_t bin, Int_t& iPad, Int_t& iTimeBin); // Get pad and time bin corresponding to the 1d bin
Float_t GetQ(const Float_t* adcArray, const Int_t time,
const Int_t pad, const Int_t maxTimeBins,
- Int_t& timeMin,Int_t& timeMax,Int_t& padMin,Int_t& padMax);
+ Int_t& timeMin,Int_t& timeMax,Int_t& padMin,Int_t& padMax) const;
+ void UpdateEventHistograms();
Int_t fFirstTimeBin; // First Time bin needed for analysis
Int_t fLastTimeBin; // Last Time bin needed for analysis
AliTPCCalPad * fOverThreshold20; //! local maxima with qMax over threshold
AliTPCCalPad * fOverThreshold30; //! local maxima with qMax over threshold
+ TProfile* fHistQVsTimeSideA; // Q vs time (side A)
+ TProfile* fHistQVsTimeSideC; // Q vs time (side C)
+ TProfile* fHistQMaxVsTimeSideA; // QMax vs time (side A)
+ TProfile* fHistQMaxVsTimeSideC; // QMax vs time (side C)
+
+ TH1F* fHistOccupancyVsEvent; // Occupancy vs event number (~time)
+ TH1F* fHistNclustersVsEvent; // Nclusters vs event number (~time)
+
Int_t fEventCounter; // event Counter
Bool_t fIsAnalysed; // Set to true after Analyse has been called
+
+ Int_t fMaxEvents; // Max events for event histograms
+ Int_t fEventsPerBin; // Events per bin for event histograms
+ Int_t fSignalCounter; // Signal counter
+ Int_t fClusterCounter; // Cluster counter
//
// Expand buffer
//
Int_t fPadsMax; //! Maximum number of time bins
Int_t fTimeBinsMax; //! Maximum number of time bins
-
-public:
- ClassDef(AliTPCdataQA, 3) // Implementation of the TPC pedestal and noise calibration
+ // DQM variables
+ Bool_t fIsDQM; //! Is DQM -> Simple output (no 2D!)
+ TProfile* fHistOccVsSector; //! Occ vs sector (for DQM only)
+ TProfile2D* fHistOcc2dVsSector; //! Occ vs sector 2D (for DQM only)
+ TProfile* fHistQVsSector; //! Q vs sector (for DQM only)
+ TProfile* fHistQmaxVsSector; //! QMax vs sector (for DQM only)
+ TArrayD* fOccVec; //! Occupancy help counter for DQM
+ TArrayD* fOccMaxVec; //! Occupancy help normlization for DQM
+ TArrayD* fOccVecFine; //! "2D" occupancy help counter for DQM
+ TArrayD* fOccMaxVecFine; //! "2D" occupancy help normlization for DQM
+
+
+ ClassDef(AliTPCdataQA, 5) // Implementation of the TPC Raw QA
};