fRefAmp(7),
fRefPoint(0),
fLatencyL1(0),
+ fLatencyL1A(0),
+ fLatencyL1C(0),
fLatencyHPTDC(0),
fVertexShift(0)
{
fRefAmp(p.fRefAmp),
fRefPoint(p.fRefPoint),
fLatencyL1(p.fLatencyL1),
- fLatencyHPTDC(p.fLatencyHPTDC),
+ fLatencyL1A(p.fLatencyL1A),
+ fLatencyL1C(p.fLatencyL1C),
+ fLatencyHPTDC(p.fLatencyHPTDC),
fVertexShift(p.fVertexShift)
{
fRefAmp = p.fRefAmp;
fRefPoint = p.fRefPoint;
fLatencyL1 = p.fLatencyL1;
+ fLatencyL1A = p.fLatencyL1A;
+ fLatencyL1C = p.fLatencyL1C;
+
fLatencyHPTDC = p.fLatencyHPTDC;
fVertexShift = p.fVertexShift;
return *this;
AliT0RecoParam *param = new AliT0RecoParam();
param->fRefAmp = 1;
param->fRefPoint = 0;
- param->fLatencyL1 = 7782.01;
+ param->fLatencyL1 = 7782.05;
+ param->fLatencyL1A = 7781.90;
+ param->fLatencyL1C = 7782.19;
param->fLatencyHPTDC = 22000;
param->fVertexShift = 4.6;
for (Int_t i=0; i<500; i++)
//
AliT0RecoParam *param = new AliT0RecoParam();
- param->fRefAmp = 5;
+ param->fRefAmp = 10;
param->fRefPoint = 0;
- param->fLatencyL1 = 7782.01;
+ param->fLatencyL1 = 7782.05;
+ param->fLatencyL1A = 7781.90;
+ param->fLatencyL1C = 7782.19;
param->fLatencyHPTDC = 22000;
+
param->fVertexShift = 4.6;
for (Int_t i=0; i<500; i++)
{
AliT0RecoParam *param = new AliT0RecoParam();
param->fRefAmp = 1;
param->fRefPoint = 1;
- param->fLatencyL1 = 7782.01;
+ param->fLatencyL1 = 7782.05;
+ param->fLatencyL1A = 7781.90;
+ param->fLatencyL1C = 7782.19;
param->fLatencyHPTDC = 22000;
+
param->fVertexShift = 0;
for (Int_t i=0; i<500; i++)
{
AliInfo(Form(" Reference point in channel : %i", fRefPoint));
AliInfo(Form(" Current latency : %i ns", fLatencyL1));
AliInfo(Form(" HPTDC latency : %i ns", fLatencyHPTDC));
- cout<<" AliT0RecoParam::PrintParameters() "<<endl;
- for (Int_t i=0; i<500; i++) cout<<i<<" "<<fLow[i]<<" "<<fHigh[i]<<endl;
+ // cout<<" AliT0RecoParam::PrintParameters() "<<endl;
+ // for (Int_t i=0; i<500; i++) cout<<i<<" "<<fLow[i]<<" "<<fHigh[i]<<endl;
}
Float_t GetLatencyL1() const {return fLatencyL1;}
void SetLatencyL1(Float_t lat) {fLatencyL1 = lat;}
+ Float_t GetLatencyL1A() const {return fLatencyL1A;}
+ void SetLatencyL1A(Float_t lat) {fLatencyL1A = lat;}
+ Float_t GetLatencyL1C() const {return fLatencyL1C;}
+ void SetLatencyL1C(Float_t lat) {fLatencyL1C = lat;}
Float_t GetLatencyHPTDC() const {return fLatencyHPTDC;}
void SetLatencyHPTDC(Float_t lat) {fLatencyHPTDC = lat;}
Float_t GetVertexShift() const {return fVertexShift;}
Float_t fLow[500]; //low limit of monitoring histograms
Float_t fHigh[500]; //high limit of monitoring histograms
Float_t fLatencyL1; //Latency L1
+ Float_t fLatencyL1A; //Latency L1 for OrA
+ Float_t fLatencyL1C; //Latency L1 for orC
Float_t fLatencyHPTDC; //Latency HPTDC
Float_t fVertexShift; // for slewing correcton
- ClassDef(AliT0RecoParam, 3);
+ ClassDef(AliT0RecoParam, 4);
};
#endif
Int_t refAmp = GetRecoParam()->GetRefAmp();
Int_t refPoint = GetRecoParam()->GetRefPoint();
Float_t latencyL1 = GetRecoParam()->GetLatencyL1();
+ Float_t latencyL1A = GetRecoParam()->GetLatencyL1A();
+ Float_t latencyL1C = GetRecoParam()->GetLatencyL1C();
Float_t latencyHPTDC = GetRecoParam()->GetLatencyHPTDC();
Float_t vertexshift = GetRecoParam()->GetVertexShift();
-
+ AliDebug(10,Form(" LatencyL1 %f latencyL1A %f latencyL1C %f latencyHPTDC %f \n",latencyL1,latencyL1A,latencyL1C, latencyHPTDC));
Int_t allData[110][5];
Int_t timeCFD[24], timeLED[24], chargeQT0[24], chargeQT1[24];
Double32_t besttimeC=9999999;
Int_t pmtBestA=99999;
Int_t pmtBestC=99999;
- Float_t meanVertex = fParam->GetMeanVertex();
+ // Float_t meanVertex = fParam->GetMeanVertex();
+ Float_t meanVertex = 0;
AliT0RecPoint* frecpoints= new AliT0RecPoint ();
}
}
if(besttimeA < 999999)
- frecpoints->SetTimeBestA( 0.001 * besttimeA * channelWidth - latencyHPTDC + latencyL1);
+ frecpoints->SetTimeBestA(0.001* besttimeA * channelWidth - latencyHPTDC + latencyL1A);
if( besttimeC < 999999 )
- frecpoints->SetTimeBestC(0.001 * besttimeC * channelWidth - latencyHPTDC + latencyL1);
+ frecpoints->SetTimeBestC( 0.001 *besttimeC * channelWidth - latencyHPTDC + latencyL1C);
AliDebug(10,Form(" pmtA %i besttimeA %f ps, pmtC %i besttimeC %f ps",
pmtBestA,besttimeA, pmtBestC, besttimeC));
if(besttimeA <999999 && besttimeC < 999999 ){
timeDiff = ( besttimeA - besttimeC) *0.001 * channelWidth;
- timeclock = 0.001 * channelWidth * Float_t(besttimeA+besttimeC)/2. - latencyHPTDC + latencyL1;
+ timeclock = 0.001*channelWidth * Float_t( besttimeA+besttimeC)/2. - latencyHPTDC + latencyL1;
meanTime = (besttimeA+besttimeC-2.*Float_t(ref))/2.;
vertex = meanVertex - c*(timeDiff)/2. + vertexshift; //+ (fdZonA - fdZonC)/2;
}