+
+ } //if phys event
+ AliDebug(1,Form(" timeDiff %f #channel, meanTime %f #ps, TOFmean%f vertex %f cm meanVertex %f \n",timeDiff, meanTime,timeclock, vertex,meanVertex));
+ frecpoints.SetT0clock(timeclock);
+ frecpoints.SetVertex(vertex);
+ frecpoints.SetMeanTime(meanTime);
+ frecpoints.SetOnlineMean(Int_t(onlineMean));
+
+ // Set triggers
+ Bool_t tr[5];
+ Int_t trchan[5] = {50,51,52,55,56};
+ Float_t lowtr[5] = {meanTVDC-700, meanOrA-700, meanOrC-700, meanOrC-1000, meanOrC-1000 };
+ Float_t hightr[5] = {meanTVDC+700, meanOrA+700, meanOrC+700, meanOrC+1000, meanOrC+1000};
+
+ for (Int_t i=0; i<5; i++) tr[i] = false;
+ for (Int_t itr=0; itr<5; itr++) {
+ for (Int_t iHit=0; iHit<1; iHit++)
+ {
+ Int_t trr=trchan[itr];
+ if( allData[trr][iHit] > lowtr[itr] && allData[trr][iHit] < hightr[itr]) tr[itr]=true;
+
+ AliDebug(15,Form("Reconstruct ::: T0 triggers iHit %i tvdc %d orA %d orC %d centr %d semicentral %d",iHit, tr[0],tr[1],tr[2],tr[3],tr[4]));
+ }
+ }
+ frecpoints.SetT0Trig(tr);
+
+ // all times with amplitude correction
+ Float_t timecent;
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ timefull = timecent = -9999;
+ tvdc = ora = orc = -9999;
+ if(allData[50][iHit]>0)
+ tvdc = (Float_t(allData[50][iHit]) - meanTVDC) * channelWidth* 0.001;
+ if(allData[51][iHit]>0)
+ ora = (Float_t(allData[51][iHit]) - meanOrA) * channelWidth* 0.001;
+
+ if(allData[52][iHit]>0)
+ orc = (Float_t(allData[52][iHit]) - meanOrC) * channelWidth* 0.001;
+
+ frecpoints.SetOrC( iHit, orc);
+ frecpoints.SetOrA( iHit, ora);
+ frecpoints.SetTVDC( iHit, tvdc);
+ for (Int_t i0=0; i0<12; i0++) {
+ if (equalize ==0 )
+ timecent = fTime0vertex[i0] + timeDelayCFD[i0];
+ else
+ timecent = fTime0vertex[i0];
+ timefull = -9999;
+ if(allData[i0+1][iHit]>1)
+ timefull = (Float_t(allData[i0+1][iHit]) - timecent)* channelWidth* 0.001;
+ frecpoints.SetTimeFull(i0, iHit,timefull) ;
+ // if(allData[i0+1][iHit]>1) printf("i0 %d iHit %d data %d fTime0vertex %f timefull %f \n",i0, iHit, allData[i0+1][iHit], fTime0vertex[i0], timefull);
+
+ }
+
+ for (Int_t i0=12; i0<24; i0++) {
+ timefull = -9999;
+ if (equalize ==0 )
+ timecent = fTime0vertex[i0] + timeDelayCFD[i0];
+ else
+ timecent = fTime0vertex[i0];
+ if(allData[i0+45][iHit]>1) {
+ timefull = (Float_t(allData[i0+45][iHit]) - timecent)* channelWidth* 0.001;
+ }
+ // if(allData[i0+45][iHit]>1) printf("i0 %d iHit %d data %d fTime0vertex %f timefull %f \n",i0, iHit, allData[i0+45][iHit], fTime0vertex[i0], timefull);
+ frecpoints.SetTimeFull(i0, iHit, timefull) ;
+ }
+ }
+
+
+ //Set MPD
+ if(allData[53][0]>0 && allData[54][0])
+ frecpoints.SetMultA(allData[53][0]-allData[54][0]);
+ if(allData[105][0]>0 && allData[106][0])
+ frecpoints.SetMultC(allData[105][0]-allData[106][0]);
+
+