1 #ifndef ALITPCCALIBALTROHEADERS_H
2 #define ALITPCCALIBALTROHEADERS_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 ////////////////////////////////////////////////////////////////////////////////////////
8 // TPC ALTRO Header analysis //
10 ////////////////////////////////////////////////////////////////////////////////////////
13 #include <TObjArray.h>
14 #include <THnSparse.h>
16 #include "AliTPCCalibRawBase.h"
17 #include "AliTPCCalPad.h"
18 #include "AliTPCROC.h"
26 class AliTPCCalibRaw : public AliTPCCalibRawBase {
29 AliTPCCalibRaw(const TMap *config);
31 virtual ~AliTPCCalibRaw();
35 virtual Int_t Update(const Int_t isector, const Int_t iRow, const Int_t iPad,
36 const Int_t iTimeBin, const Float_t signal);
37 virtual void UpdateDDL();
38 virtual void EndEvent();
39 virtual void ResetEvent();
40 virtual void Analyse();
42 UInt_t GetNFailL1Phase() const {return fNFailL1Phase;}
43 UInt_t GetNFailL1PhaseEvents() const {return fNFailL1PhaseEvent;}
44 Int_t GetPeakDetectionMinus() const {return fPeakDetMinus;}
45 Int_t GetPeakDetectionPlus() const {return fPeakDetPlus;}
47 const TVectorF* GetALTROL1PhaseEvents() const {return &fArrALTROL1Phase;}
49 const TVectorF *GetALTROL1PhaseEventsRCU(Int_t rcu) const {return (TVectorF*)fArrALTROL1PhaseEvent.At(rcu);}
50 const TVectorF *GetALTROL1PhaseFailEventsRCU(Int_t rcu) const {return (TVectorF*)fArrALTROL1PhaseFailEvent.At(rcu);}
52 const TVectorF *GetOccupancyEvent() const {return &fVOccupancyEvent;}
53 const TVectorF *GetOccupancyEventSensitive() const {return &fVOccupancySenEvent;}
54 const TVectorF *GetSignalSumEvent() const {return &fVSignalSumEvent;}
55 const TVectorF *GetSignalSumEventSensitive() const {return &fVSignalSumSenEvent;}
56 const TVectorF *GetFiredPadsSensitive() const {return &fVNfiredPadsSenEvent;}
57 const TVectorF *GetEventTimeStamps() const {return &fVTimeStampEvent;}
58 UInt_t GetFirstTimeStamp() const {return fFirstTimeStamp;}
60 void SetRangePeakDetection(Int_t minus, Int_t plus) { fPeakDetMinus=minus; fPeakDetPlus=plus;}
62 TH2C *MakeHistL1RCUEvents(Int_t type=0);
63 TH2C *MakeHistL1RCUEventsIROC(Int_t type=0);
64 TH2C *MakeHistL1RCUEventsOROC(Int_t type=0);
65 TH1F *MakeHistL1PhaseDist();
66 TVectorF *MakeVectL1PhaseDist();
68 TGraph* MakeGraphOccupancy(const Int_t type=0, const Int_t xType=0);
69 TGraph* MakeGraphNoiseEvents();
70 TCanvas* MakeCanvasOccupancy(const Int_t xType=1, Bool_t sen=kFALSE);
72 const THnSparseI *GetHnDrift() const {return fHnDrift;}
73 // AliTPCCalPad *CreateCalPadL1Mean();
74 // AliTPCCalPad *CreateCalPadL1RMS();
77 Int_t fPeakDetMinus; // Consecutive timebins on rising edge to be regarded as a signal
78 Int_t fPeakDetPlus; // Consecutive timebins on falling edge to be regarded as a signal
79 UInt_t fNFailL1Phase; //Number of failures in L1 phase
80 UInt_t fNFailL1PhaseEvent; //Number of events with L1 phase failures
81 UInt_t fFirstTimeStamp; //Time Stamp from first event
83 UInt_t fNSecTime; //Number of seconds per bin in time
84 UInt_t fNBinsTime; //Number of bin in time
85 //processing information
86 Bool_t fPadProcessed; //! if last pead has already been filled for the current pad
87 Int_t fCurrentChannel; //! current channel processed
88 Int_t fCurrentSector; //! current sector processed
89 Int_t fLastSector; //! current sector processed
90 Int_t fCurrentRow; //! current row processed
91 Int_t fCurrentPad; //! current pad processed
92 Int_t fLastTimeBinProc; //! last time bin processed
93 Int_t fPeakTimeBin; //! time bin with local maximum
94 Int_t fLastSignal; //! last signal processed
95 Int_t fNOkPlus; //! number of processed time bins fullfilling peak criteria
96 Int_t fNOkMinus; //! number of processed time bins fullfilling peak criteria
97 Int_t fNanoSec; //! current nano seconds stamp
100 TVectorF fArrCurrentPhaseDist; //!Phase distribution of the current event
101 TVectorF fArrCurrentPhase; //!Current phase of all RCUs
102 TVectorF fArrFailEventNumber; //event numbers of failed events;
103 TVectorF fArrALTROL1Phase; //Array of L1 phases on an event bases;
104 TObjArray fArrALTROL1PhaseEvent; //L1 phase for each RCU and event
105 TObjArray fArrALTROL1PhaseFailEvent; //L1 failure for each RCU and event
106 //drift velocity stuff
108 THnSparseI *fHnDrift; //Histogram last time bin vs. ROC, Time
110 TVectorF fVOccupancyEvent; //occupancy per event (number of samples above threshold)
111 TVectorF fVSignalSumEvent; //occupancy per event (sum of all adc values)
112 TVectorF fVOccupancySenEvent; //occupancy per event (number of samples abouve threshold) in sensitive regions
113 TVectorF fVSignalSumSenEvent; //occupancy per event (sum of all adc values) in sensitive regions
114 TVectorF fVNfiredPadsSenEvent; //number of pads with a signal above threshold in sensitive regions
115 TVectorF fVTimeStampEvent; //timestamp for all events
117 TVectorF *MakeArrL1PhaseRCU(Int_t rcu, Bool_t force=kFALSE);
118 TVectorF *MakeArrL1PhaseFailRCU(Int_t rcu, Bool_t force=kFALSE);
120 Bool_t IsEdgePad(Int_t sector, Int_t row, Int_t pad) const;
123 AliTPCCalibRaw(AliTPCCalibRaw &calib);
124 AliTPCCalibRaw& operator = (const AliTPCCalibRaw &source);
126 ClassDef(AliTPCCalibRaw,3) // Analysis of the Altro header information
129 //----------------------
131 //----------------------
132 inline TVectorF *AliTPCCalibRaw::MakeArrL1PhaseRCU(Int_t rcu, Bool_t force)
134 TVectorF *arr=(TVectorF*)fArrALTROL1PhaseEvent.UncheckedAt(rcu);
136 arr=new TVectorF(100);
137 fArrALTROL1PhaseEvent.AddAt(arr,rcu);
142 inline TVectorF *AliTPCCalibRaw::MakeArrL1PhaseFailRCU(Int_t rcu, Bool_t force)
144 TVectorF *arr=(TVectorF*)fArrALTROL1PhaseFailEvent.UncheckedAt(rcu);
146 arr=new TVectorF(100);
147 fArrALTROL1PhaseFailEvent.AddAt(arr,rcu);
151 //_____________________________________________________________________
152 inline Bool_t AliTPCCalibRaw::IsEdgePad(Int_t sector, Int_t row, Int_t pad) const
155 // return true if pad is on the edge of a row
158 if ( pad == edge1 ) return kTRUE;
159 Int_t edge2 = fROC->GetNPads(sector,row)-1;
160 if ( pad == edge2 ) return kTRUE;