1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
17 ///////////////////////////////////////////////////////////////////////////////
19 // Class with T0 reconstruction parameters //
21 // for low flux time-amplitude correction equalize time to amplitude 1 MIP;
22 // for high flux - to 15MIP
23 // To have nice time spectra after reconstruction we need to know
24 // reference point to write t(i) - RefPoint.
25 // It can be apparatus RefPoint or one of PMT //
26 // fRefPoint - number of channel with RF
28 // Alla.Maevskaya@cern.ch
29 /////////////////////////////////////////////////////////////////////////
33 #include "AliT0RecoParam.h"
34 #include "Riostream.h"
38 ClassImp(AliT0RecoParam)
43 //_____________________________________________________________________________
44 AliT0RecoParam::AliT0RecoParam():
45 AliDetectorRecoParam(),
60 fSatelliteThresholds[0] = -15;
61 fSatelliteThresholds[1] = -1.5;
62 for (Int_t i=0; i<500; i++) {
65 if( i<24) fBadChannels[i]=-1;
70 //_____________________________________________________________________________
71 AliT0RecoParam::~AliT0RecoParam()
78 //_____________________________________________________________________________
80 AliT0RecoParam::AliT0RecoParam(const AliT0RecoParam &p):
81 AliDetectorRecoParam(p),
83 fRefPoint(p.fRefPoint),
84 fLatencyL1(p.fLatencyL1),
85 fLatencyL1A(p.fLatencyL1A),
86 fLatencyL1C(p.fLatencyL1C),
87 fLatencyHPTDC(p.fLatencyHPTDC),
88 fVertexShift(p.fVertexShift),
89 fEqualised( p.fEqualised)
92 fSatelliteThresholds[0] = (p.fSatelliteThresholds[0]);
93 fSatelliteThresholds[1] = (p.fSatelliteThresholds[1]);
94 for (Int_t i=0; i<500; i++) {
96 fHigh[i] = p.fHigh[i];
97 if( i<24) fBadChannels[i] = p.fBadChannels[i];
102 //_____________________________________________________________________________
104 AliT0RecoParam& AliT0RecoParam:: operator=(const AliT0RecoParam &p)
113 AliDetectorRecoParam::operator=(p);
115 fRefPoint = p.fRefPoint;
116 fLatencyL1 = p.fLatencyL1;
117 fLatencyL1A = p.fLatencyL1A;
118 fLatencyL1C = p.fLatencyL1C;
120 fLatencyHPTDC = p.fLatencyHPTDC;
121 fVertexShift = p.fVertexShift;
123 fSatelliteThresholds[0] = (p.fSatelliteThresholds[0]);
124 fSatelliteThresholds[1] = (p.fSatelliteThresholds[1]);
125 fEqualised = p.fEqualised;
126 for (Int_t i=0; i<500; i++) {
128 fHigh[i] = p.fHigh[i];
129 if( i<24) fBadChannels[i] = p.fBadChannels[i];
135 //_____________________________________________________________________________
137 AliT0RecoParam *AliT0RecoParam::GetLowFluxParam()
140 // make default reconstruction parameters for low flux env.
142 AliT0RecoParam *param = new AliT0RecoParam();
144 param->fRefPoint = 0;
145 param->fLatencyL1 = 7782.05;
146 param->fLatencyL1A = 7781.90;
147 param->fLatencyL1C = 7782.19;
148 param->fLatencyHPTDC = 22000;
149 param->fVertexShift = 0;
150 for (Int_t i=0; i<500; i++)
153 param-> fHigh[i]=10000.;
155 param->SetName("Low Flux");
156 param->SetTitle("Low Flux");
157 param->SetSatelliteThresholds(-15, -1.5);
162 //_____________________________________________________________________________
164 AliT0RecoParam *AliT0RecoParam::GetHighFluxParam()
167 // make reco parameters for high flux env.
170 AliT0RecoParam *param = new AliT0RecoParam();
172 param->fRefPoint = 0;
173 param->fLatencyL1 = 7782.05;
174 param->fLatencyL1A = 7781.90;
175 param->fLatencyL1C = 7782.19;
176 param->fVertexShift = 0;
177 param->fLatencyHPTDC = 22000;
178 for (Int_t i=0; i<500; i++)
181 param-> fHigh[i]=20000.;
184 param->SetSatelliteThresholds(-15, -1.5);
187 param->SetName("High Flux");
188 param->SetTitle("High Flux");
193 //_____________________________________________________________________________
195 AliT0RecoParam *AliT0RecoParam::GetLaserTestParam()
198 // special setting for laser
200 AliT0RecoParam *param = new AliT0RecoParam();
202 param->fRefPoint = 1;
203 param->fLatencyL1 = 7782.05;
204 param->fLatencyL1A = 7781.90;
205 param->fLatencyL1C = 7782.19;
206 param->fLatencyHPTDC = 22000;
207 param->fVertexShift = 0;
208 param->SetSatelliteThresholds(-15, -1.5);
211 for (Int_t i=0; i<500; i++)
214 param-> fHigh[i]=12000.;
217 param->SetName("Laser Flux");
218 param->SetTitle("Laser Flux");
221 //_____________________________________________________________________________
223 void AliT0RecoParam::PrintParameters() const
226 // Printing of the used T0 reconstruction parameters
228 AliInfo(Form(" Reference amplitude for walk corerection : %f", fRefAmp));
229 AliInfo(Form(" Bad channel in channel : %i", fRefPoint));
230 cout<<" AliT0RecoParam::PrintParameters() "<<endl;
231 for (Int_t i=0; i<105; i++) cout<<i<<" "<<fLow[i]<<" "<<fHigh[i]<<endl;