- if(option == "pdc"){
- Double_t qt0 = Double_t(chargeQT0[ipmt]);
- Double_t qt1 = Double_t(chargeQT1[ipmt]);
- if((qt1-qt0)>0) adc[ipmt] = Int_t(TMath::Exp( Double_t (channelWidth*(qt1-qt0)/1000.)));
- time[ipmt] = fCalib-> WalkCorrection( ipmt,Int_t(qt1) , timeCFD[ipmt], "pdc" ) ;
- Double_t sl = (timeLED[ipmt] - time[ipmt])*channelWidth;
- if(fAmpLEDrec.At(ipmt))
- qt=((TGraph*)fAmpLEDrec.At(ipmt))->Eval(sl/1000.);
- frecpoints->SetTime(ipmt,time[ipmt]);
- frecpoints->SetAmp(ipmt,adc[ipmt]/Float_t(mv2MIP));
- frecpoints->SetAmpLED(ipmt,qt/Float_t(mv2MIP));
- AliDebug(10,Form(" QTC %f mv, time in chann %f ampLED %f",adc[ipmt] ,time[ipmt], qt));
- AliDebug(10,Form(" Amlitude in MIPS LED %f , QTC %f \n ", adc[ipmt]/Float_t(mv2MIP),qt/Float_t(mv2MIP)));
- }
- if(option == "cosmic") {
- // if(ipmt == 15) continue; //skip crashed PMT
- if(( chargeQT1[ipmt] - chargeQT0[ipmt])>0)
- adc[ipmt] = chargeQT1[ipmt] - chargeQT0[ipmt];
- else
- adc[ipmt] = 0;
- // time[ipmt] = fCalib-> WalkCorrection( ipmt, adc[ipmt], timeCFD[ipmt],"cosmic" ) ;
-
- Double_t sl = timeLED[ipmt] - timeCFD[ipmt];
- time[ipmt] = fCalib-> WalkCorrection( ipmt, Int_t(sl), timeCFD[ipmt],"cosmic" ) ;
- // if(fAmpLEDrec.At(ipmt))
- // qt=((TGraph*)fAmpLEDrec.At(ipmt))->Eval(sl);
- time[ipmt] = time[ipmt] - allData[0][0] + 5000;
- AliDebug(10,Form(" ipmt %i QTC %i , time in chann %i (led-cfd) %i ",
- ipmt, Int_t(adc[ipmt]) ,Int_t(time[ipmt]),Int_t( sl)));
- Double_t ampMip =( (TGraph*)fAmpLED.At(ipmt))->Eval(sl);
- Double_t qtMip = ((TGraph*)fQTC.At(ipmt))->Eval(adc[ipmt]);
- AliDebug(10,Form(" Amlitude in MIPS LED %f , QTC %f \n ",ampMip,qtMip));
-
- frecpoints->SetTime(ipmt, Float_t(time[ipmt]) );
- frecpoints->SetAmp(ipmt, Float_t( ampMip)); //for cosmic &pp beam
- frecpoints->SetAmpLED(ipmt, Float_t(qtMip));
-