Use event specie to identufy laser events
[u/mrichter/AliRoot.git] / TPC / AliTPCCalibRaw.h
1 #ifndef ALITPCCALIBRAW_H
2 #define ALITPCCALIBRAW_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 ////////////////////////////////////////////////////////////////////////////////////////
7 //                                                                                    //
8 //                        TPC ALTRO Header analysis                                   //
9 //                                                                                    //
10 ////////////////////////////////////////////////////////////////////////////////////////
11
12 #include <TVectorF.h>
13 #include <TObjArray.h>
14 #include <THnSparse.h>
15
16 #include "AliTPCCalibRawBase.h"
17 #include "AliTPCCalPad.h"
18 #include "AliTPCROC.h"
19
20 class TH2C;
21 class TH1F;
22 class TMap;
23 class TGraph;
24 class TCanvas;
25
26 class AliTPCCalibRaw : public AliTPCCalibRawBase {
27 public:
28   AliTPCCalibRaw();
29   AliTPCCalibRaw(const TMap *config);
30   
31   virtual ~AliTPCCalibRaw();
32
33   enum {kNRCU=216};
34   
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();
41
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;}
46   
47   const TVectorF* GetALTROL1PhaseEvents() const {return &fArrALTROL1Phase;}
48
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);}
51
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;}
59   
60   void  SetRangePeakDetection(Int_t minus, Int_t plus) { fPeakDetMinus=minus; fPeakDetPlus=plus;}
61   //Phase info
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();
67   //Occupancy info
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);
71
72   const THnSparseI *GetHnDrift() const {return fHnDrift;}
73 //   AliTPCCalPad *CreateCalPadL1Mean();
74 //   AliTPCCalPad *CreateCalPadL1RMS();
75   
76   void Merge(AliTPCCalibRaw * const sig);
77   virtual Long64_t Merge(TCollection * const list);
78   
79 private:
80   Int_t   fPeakDetMinus;             //  Consecutive timebins on rising edge to be regarded as a signal
81   Int_t   fPeakDetPlus;              //  Consecutive timebins on falling edge to be regarded as a signal
82   UInt_t  fNFailL1Phase;             //Number of failures in L1 phase
83   UInt_t  fNFailL1PhaseEvent;        //Number of events with L1 phase failures
84   //binning dv hist
85   UInt_t  fNSecTime;                 //Number of seconds per bin in time
86   UInt_t  fNBinsTime;                //Number of bin in time
87   //processing information
88   Bool_t    fPadProcessed;           //! if last pead has already been filled for the current pad
89   Int_t     fCurrentChannel;         //! current channel processed
90   Int_t     fCurrentSector;          //! current sector processed
91   Int_t     fLastSector;             //! current sector processed
92   Int_t     fCurrentRow;             //! current row processed
93   Int_t     fCurrentPad;             //! current pad processed
94   Int_t     fLastTimeBinProc;        //! last time bin processed
95   Int_t     fPeakTimeBin;            //! time bin with local maximum
96   Int_t     fLastSignal;             //! last signal processed
97   Int_t     fNOkPlus;                //! number of processed time bins fullfilling peak criteria
98   Int_t     fNOkMinus;               //! number of processed time bins fullfilling peak criteria
99   Int_t     fNanoSec;                //! current nano seconds stamp
100 //
101   //L1 phase stuff
102   TVectorF fArrCurrentPhaseDist;       //!Phase distribution of the current event
103   TVectorF fArrCurrentPhase;           //!Current phase of all RCUs
104   TVectorF fArrFailEventNumber;        //event numbers of failed events;
105   TVectorF fArrALTROL1Phase;           //Array of L1 phases on an event bases;
106   TObjArray fArrALTROL1PhaseEvent;     //L1 phase for each RCU and event
107   TObjArray fArrALTROL1PhaseFailEvent; //L1 failure for each RCU and event
108   //drift velocity stuff
109   enum {kHnBinsDV=3};
110   THnSparseI *fHnDrift;                //Histogram last time bin vs. ROC, Time
111   //occupancy
112   TVectorF fVOccupancyEvent;           //occupancy per event (number of samples above threshold)
113   TVectorF fVSignalSumEvent;           //occupancy per event (sum of all adc values)
114   TVectorF fVOccupancySenEvent;        //occupancy per event (number of samples abouve threshold) in sensitive regions
115   TVectorF fVSignalSumSenEvent;        //occupancy per event (sum of all adc values) in sensitive regions
116   TVectorF fVNfiredPadsSenEvent;       //number of pads with a signal above threshold in sensitive regions
117   TVectorF fVTimeStampEvent;           //timestamp for all events
118   
119   TVectorF *MakeArrL1PhaseRCU(Int_t rcu, Bool_t force=kFALSE);
120   TVectorF *MakeArrL1PhaseFailRCU(Int_t rcu, Bool_t force=kFALSE);
121   
122   Bool_t IsEdgePad(Int_t sector, Int_t row, Int_t pad) const;
123   void CreateDVhist();
124   
125   AliTPCCalibRaw(const AliTPCCalibRaw &calib);
126   AliTPCCalibRaw& operator = (const  AliTPCCalibRaw &source);
127
128   ClassDef(AliTPCCalibRaw,4) //  Analysis of the Altro header information
129 };
130
131 //----------------------
132 // Inline Functions
133 //----------------------
134 inline TVectorF *AliTPCCalibRaw::MakeArrL1PhaseRCU(Int_t rcu, Bool_t force)
135 {
136   TVectorF *arr=(TVectorF*)fArrALTROL1PhaseEvent.UncheckedAt(rcu);
137   if (!arr && force) {
138     arr=new TVectorF(100);
139     fArrALTROL1PhaseEvent.AddAt(arr,rcu);
140   }
141   return arr;
142 }
143 //
144 inline TVectorF *AliTPCCalibRaw::MakeArrL1PhaseFailRCU(Int_t rcu, Bool_t force)
145 {
146   TVectorF *arr=(TVectorF*)fArrALTROL1PhaseFailEvent.UncheckedAt(rcu);
147   if (!arr && force) {
148     arr=new TVectorF(100);
149     fArrALTROL1PhaseFailEvent.AddAt(arr,rcu);
150   }
151   return arr;
152 }
153 //_____________________________________________________________________
154 inline Bool_t AliTPCCalibRaw::IsEdgePad(Int_t sector, Int_t row, Int_t pad) const
155 {
156   //
157   // return true if pad is on the edge of a row
158   //
159   Int_t edge1   = 0;
160   if ( pad == edge1 ) return kTRUE;
161   Int_t edge2   = fROC->GetNPads(sector,row)-1;
162   if ( pad == edge2 ) return kTRUE;
163   
164   return kFALSE;
165 }
166
167
168 #endif