Int_t ph2Mip = fParam->GetPh2Mip();
- Int_t channelWidth = fParam->GetChannelWidth() ;
+ Float_t channelWidth = fParam->GetChannelWidth() ;
Float_t delayVertex = fParam->GetTimeDelayTVD();
zdetC = TMath::Abs(fParam->GetZPosition("T0/C/PMT1"));
sl = gr->Eval(qt);
trLED = Int_t(( timeGaus[i] + 1000*sl )/channelWidth);
- qtCh=Int_t (1000.*TMath::Log(qt)) / channelWidth;
+ qtCh=Int_t (1000.*TMath::Log(qt) / channelWidth);
fADC0->AddAt(0,i);
fADC->AddAt(qtCh,i);
ftimeLED->AddAt(trLED,i);
// Set various `Fixed' parameters
void SetPh2Mip(Int_t r=300) { fPh2Mip = r; }
void SetmV2Mip(Int_t r=50) { fmV2Mip = r; }
- void SetChannelWidth(Int_t s=24.4) { fChannelWidth = s;}
+ void SetChannelWidth(Float_t s=24.4) { fChannelWidth = s;}
void SetmV2channel(Int_t size=320) { fmV2Channel = size; }
void SetQTmin(Int_t qt=13) {fQTmin = qt;}
void SetQTmax(Int_t qt=125) {fQTmax = qt;}
// Get `Fixed' various parameters
Int_t GetPh2Mip() const { return fPh2Mip; }
Int_t GetmV2Mip() const { return fmV2Mip; }
- Int_t GetChannelWidth() const { return fChannelWidth; }
+ Float_t GetChannelWidth() const { return fChannelWidth; }
Int_t GetmV2channel() const { return fmV2Channel; }
Int_t GetQTmin() const {return fQTmin;}
Int_t GetQTmax() const {return fQTmax;}
Float_t fT0zPosition[2] ; // z-position of the two T0s
Int_t fPh2Mip; // # photoelectrons per MIP in radiator
Int_t fmV2Mip; // # mV per MIP in radiator
- Int_t fChannelWidth; // channel width in ns
+ Float_t fChannelWidth; // channel width in ns
Int_t fmV2Channel; // ADC mv 2 channel # (200000ps/(25*25).
Int_t fQTmin; //min time for QTC
Int_t fQTmax; //max time fro QTC
// Int_t mV2Mip = param->GetmV2Mip();
//mV2Mip = param->GetmV2Mip();
- Int_t channelWidth = fParam->GetChannelWidth() ;
+ Float_t channelWidth = fParam->GetChannelWidth() ;
Int_t meanT0 = fParam->GetMeanT0();
AliDebug(1,Form("Start DIGITS reconstruction "));
// Int_t mV2Mip = param->GetmV2Mip();
//mV2Mip = param->GetmV2Mip();
- Int_t channelWidth = fParam->GetChannelWidth() ;
+ Float_t channelWidth = fParam->GetChannelWidth() ;
+
Int_t meanT0 = fParam->GetMeanT0();
for (Int_t in=0; in<24; in++)
frecpoints->SetTime(ipmt,time[ipmt]);
frecpoints->SetAmp(ipmt,adc[ipmt]);
frecpoints->SetAmpLED(ipmt,qt);
+ AliDebug(1,Form(" QTC %f mv, QTC %f MIPS time in chann %f time %f ",adc[ipmt], adc[ipmt]/50.,time[ipmt], time[ipmt]*channelWidth));
+
}
else {
time[ipmt] = 0;
Float_t vertex = 0;
if(besttimeA !=999999 && besttimeC != 999999 ){
timeDiff = (besttimeC - besttimeA)*channelWidth;
- // meanTime = (besttimeA + besttimeC)/2.;
- meanTime = (meanT0 - (besttimeA + besttimeC)/2) * channelWidth;
+ meanTime = (meanT0 - (besttimeA + besttimeC)/2) * channelWidth;
vertex = c*(timeDiff)/2. + (fdZ_A - fdZ_C)/2; //-(lenr-lenl))/2;
AliDebug(1,Form(" timeDiff %f ps, meanTime %f ps, vertex %f cm",timeDiff, meanTime,vertex ));
frecpoints->SetVertex(vertex);
}
//time in each channel as time[ipmt]-MeanTimeinThisChannel(with vertex=0)
+ /*
for (Int_t ipmt=0; ipmt<24; ipmt++) {
if(time[ipmt]>1) {
time[ipmt] = (time[ipmt] - fTime0vertex[ipmt])*channelWidth;
frecpoints->SetTime(ipmt,time[ipmt]);
}
- }
+ }
+ */
recTree->Fill();