- AliDebug(10," before read data ");
- AliT0RawReader myrawreader(rawReader);
- if (!myrawreader.Next())
- AliDebug(1,Form(" no raw data found!!"));
- else
- {
- for (Int_t i=0; i<105; i++) {
- for (Int_t iHit=0; iHit<5; iHit++)
- {
- allData[i][iHit] = myrawreader.GetData(i,iHit);
- }
- }
-
- Float_t channelWidth = fParam->GetChannelWidth() ;
-
- // Int_t meanT0 = fParam->GetMeanT0();
- if(option == "pdc"){
- for (Int_t in=0; in<24; in++)
- {
-
- timeLED[in] = allData[in+1][0] ;
- timeCFD[in] = allData[in+25][0] ;
- chargeQT1[in] = allData[in+57][0] ;
- chargeQT0[in] = allData[in+80][0] ;
- }
- }
-
- if(option == "cosmic") {
- for (Int_t in=0; in<12; in++)
- {
- timeCFD[in] = allData[in+1][0] ;
- timeCFD[in+12] = allData[in+56+1][0] ;
- timeLED[in] = allData[in+12+1][0] ;
- timeLED[in+12] = allData[in+68+1][0] ;
- }
-
- for (Int_t in=0; in<24; in=in+2)
- {
- Int_t cc=in/2;
- chargeQT1[cc]=allData[in+25][0];
- chargeQT0[cc]=allData[in+26][0];
- }
- for (Int_t in=24; in<48; in=in+2)
- {
- Int_t cc=in/2;
- chargeQT1[cc]=allData[in+57][0];
- chargeQT0[cc]=allData[in+58][0];
+ AliDebug(10," before read data ");
+ AliT0RawReader myrawreader(rawReader);
+
+ UInt_t type =rawReader->GetType();
+
+ if (!myrawreader.Next())
+ AliDebug(1,Form(" no raw data found!!"));
+ else
+ {
+ for (Int_t i=0; i<24; i++)
+ {
+ timeCFD[i]=0; timeLED[i]=0; chargeQT0[i]=0; chargeQT1[i]=0;
+ }
+ Int_t fBCID=Int_t (rawReader->GetBCID());
+ Int_t trmbunch= myrawreader.GetTRMBunchID();
+ AliDebug(10,Form(" CDH BC ID %i, TRM BC ID %i \n", fBCID, trmbunch ));
+
+ if(type == 7 ) { //only physics
+ for (Int_t i=0; i<105; i++) {
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ allData[i][iHit] = myrawreader.GetData(i,iHit);
+ }
+ }
+ Int_t ref=0;
+ if (refPoint>0)
+ ref = allData[refPoint][0]-5000;
+
+ Float_t channelWidth = fParam->GetChannelWidth() ;
+
+ // Int_t meanT0 = fParam->GetMeanT0();
+
+ for (Int_t in=0; in<12; in++)
+ {
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ if(allData[in+1][iHit] > low[in+1] &&
+ allData[in+1][iHit] < high[in+1])
+ {
+ timeCFD[in] = allData[in+1][iHit] ;
+ break;
+ }
+ }
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ if(allData[in+12+1][iHit] > low[in+12+1] &&
+ allData[in+1+12][iHit] < high[in+12+1])
+ {
+ timeLED[in] = allData[in+12+1][iHit] ;
+ break;
+ }
+ }
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ if(allData[in+1+56][iHit] > low[in+1+56] &&
+ allData[in+1+56][iHit] < high[in+1+56])
+ {
+ timeCFD[in+12] = allData[in+56+1][iHit] ;
+ break;
+ }
+ }
+
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ if(allData[in+1+68][iHit] > low[in+1+68] &&
+ allData[in+1+68][iHit] < high[in+1+68])
+ {
+ timeLED[in+12] = allData[in+68+1][iHit] ;
+ break;
+ }
+ }
+ AliDebug(10, Form(" readed i %i cfdC %i cfdA %i ledC %i ledA%i ",
+ in, timeCFD[in],timeCFD[in+12],timeLED[in],
+ timeLED[in+12]));
+
+ }
+
+
+ for (Int_t in=0; in<12; in++)
+ {
+ chargeQT0[in]=allData[2*in+25][0];
+ chargeQT1[in]=allData[2*in+26][0];
+ AliDebug(10, Form(" readed Raw %i %i %i",
+ in, chargeQT0[in],chargeQT1[in]));
+ }
+ for (Int_t in=12; in<24; in++)
+ {
+ chargeQT0[in]=allData[2*in+57][0];
+ chargeQT1[in]=allData[2*in+58][0];
+ AliDebug(10, Form(" readed Raw %i %i %i",
+ in, chargeQT0[in],chargeQT1[in]));
+
+ }
+
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ if(allData[49][iHit] > low[49] &&
+ allData[49][iHit] < high[49]){
+ onlineMean = allData[49][iHit];
+ break;
+ }
+ }
+ Double32_t time[24], adc[24], noncalibtime[24];
+ for (Int_t ipmt=0; ipmt<24; ipmt++) {
+ if(timeCFD[ipmt]>0 ){
+ //for simulated data
+ //for physics data
+ if(( chargeQT0[ipmt] - chargeQT1[ipmt])>0) {
+ adc[ipmt] = chargeQT0[ipmt] - chargeQT1[ipmt];