+ if(ndig == 0) continue;
+ for(Int_t k=0; k<ndig; k++){
+ AliVZEROdigit* fVZERODigit = (AliVZEROdigit*) VZEROdigits->At(k);
+ // Convert aliroot channel k into FEE channel iChannel before writing data
+ Int_t iChannel = AliVZEROCalibData::GetBoardNumber(fVZERODigit->PMNumber()) * 8 +
+ AliVZEROCalibData::GetFEEChannelNumber(fVZERODigit->PMNumber());
+ for(Int_t iClock = 0; iClock < AliVZEROdigit::kNClocks; ++iClock) aADC[iChannel][iClock] = fVZERODigit->ChargeADC(iClock);
+ aTime[iChannel] = fVZERODigit->Time();
+ aWidth[iChannel] = fVZERODigit->Width();
+ aIntegrator[iChannel]= fVZERODigit->Integrator();
+ if(fVZERODigit->PMNumber() < 32) {
+ aBBflag[iChannel] = (aBBflagsV0C >> fVZERODigit->PMNumber()) & 0x1;
+ aBGflag[iChannel] = (aBGflagsV0C >> fVZERODigit->PMNumber()) & 0x1;
+ }
+ else {
+ aBBflag[iChannel] = (aBBflagsV0A >> (fVZERODigit->PMNumber()-32)) & 0x1;
+ aBGflag[iChannel] = (aBGflagsV0A >> (fVZERODigit->PMNumber()-32)) & 0x1;
+ }
+ AliDebug(1,Form("DDL: %s\tdigit number: %d\tPM number: %d\tADC: %d\tTime: %f",
+ fileName,k,fVZERODigit->PMNumber(),aADC[k][AliVZEROdigit::kNClocks/2],aTime[k]));
+ }
+ }
+
+ // Now fill raw data
+
+ for (Int_t iCIU = 0; iCIU < 8; iCIU++) {
+
+ // decoding of one Channel Interface Unit numbered iCIU - there are 8 channels per CIU (and 8 CIUs) :
+
+ for(Int_t iChannel_Offset = iCIU*8; iChannel_Offset < (iCIU*8)+8; iChannel_Offset=iChannel_Offset+4) {
+ for(Int_t iChannel = iChannel_Offset; iChannel < iChannel_Offset+4; iChannel++) {
+ buffer->WriteChannel(iChannel, aADC[iChannel], aIntegrator[iChannel]);
+ }
+ buffer->WriteBeamFlags(&aBBflag[iChannel_Offset],&aBGflag[iChannel_Offset]);
+ buffer->WriteMBInfo();
+ buffer->WriteMBFlags();
+ buffer->WriteBeamScalers();
+ }
+
+ for(Int_t iChannel = iCIU*8 + 7; iChannel >= iCIU*8; iChannel--) {
+ buffer->WriteTiming(aTime[iChannel], aWidth[iChannel]);
+ }
+
+ // End of decoding of one CIU card
+
+ } // end of decoding the eight CIUs
+
+ delete buffer;
+ fLoader->UnloadDigits();
+}
+
+//_____________________________________________________________________________
+Bool_t AliVZERO::Raw2SDigits(AliRawReader* rawReader){
+ // Converts the VZERO raw data into digits
+ // The method is used for merging simulated and
+ // real data events
+ TStopwatch timer;
+ timer.Start();
+
+ if(!fLoader) {
+ AliError("no VZERO loader found");
+ return kFALSE;
+ }
+ fLoader->LoadSDigits("UPDATE");
+
+ if (!fLoader->TreeS()) fLoader->MakeTree("S");
+ fLoader->MakeSDigitsContainer();
+ TTree* treeS = fLoader->TreeS();
+
+ TClonesArray *sdigits = new TClonesArray("AliVZEROSDigit", 64);
+ treeS->Branch("VZEROSDigit", &sdigits);
+
+ {
+ rawReader->Reset();
+ AliVZERORawStream rawStream(rawReader);
+
+ if (!rawStream.Next()) return kFALSE; // No VZERO data found
+
+ GetCalibData();
+
+ Int_t nSDigits = 0;
+ Float_t *charges = NULL;
+ Int_t nbins = 0;
+ for(Int_t iChannel=0; iChannel < 64; ++iChannel) {
+ Int_t offlineCh = rawStream.GetOfflineChannel(iChannel);
+ Short_t chargeADC[AliVZEROdigit::kNClocks];
+ for(Int_t iClock=0; iClock < AliVZEROdigit::kNClocks; ++iClock) {
+ chargeADC[iClock] = rawStream.GetPedestal(iChannel,iClock);
+ }
+ // Integrator flag
+ Bool_t integrator = rawStream.GetIntegratorFlag(iChannel,AliVZEROdigit::kNClocks/2);
+ // HPTDC data (leading time and width)
+ Int_t board = AliVZEROCalibData::GetBoardNumber(offlineCh);
+ Float_t time = rawStream.GetTime(iChannel)*fCalibData->GetTimeResolution(board);
+ // Float_t width = rawStream.GetWidth(iChannel)*fCalibData->GetWidthResolution(board);
+ Float_t adc = 0;
+
+ // Pedestal retrieval and suppression
+ Float_t maxadc = 0;
+ Int_t imax = -1;
+ Float_t adcPedSub[AliVZEROdigit::kNClocks];
+ Float_t integral = fSignalShape->Integral(0,200);
+ for(Int_t iClock=0; iClock < AliVZEROdigit::kNClocks; ++iClock) {
+ Bool_t iIntegrator = (iClock%2 == 0) ? integrator : !integrator;
+ Int_t k = offlineCh + 64*iIntegrator;
+ adcPedSub[iClock] = (Float_t)chargeADC[iClock] - fCalibData->GetPedestal(k);
+ if(adcPedSub[iClock] <= fRecoParam->GetNSigmaPed()*fCalibData->GetSigma(k)) {
+ adcPedSub[iClock] = 0;
+ continue;
+ }
+ if(iClock < fRecoParam->GetStartClock() || iClock > fRecoParam->GetEndClock()) continue;
+ if(adcPedSub[iClock] > maxadc) {
+ maxadc = adcPedSub[iClock];
+ imax = iClock;
+ }
+ }
+ if (imax != -1) {
+ Int_t start = imax - fRecoParam->GetNPreClocks();
+ if (start < 0) start = 0;
+ Int_t end = imax + fRecoParam->GetNPostClocks();
+ if (end > 20) end = 20;
+ for(Int_t iClock = start; iClock <= end; iClock++) {
+ adc += adcPedSub[iClock];