+ const int knADCData1=24, knADCData2=24; // In principle the 2 numbers can be different!
+ UInt_t lADCHeader1;
+ UInt_t lADCHeader2;
+ UInt_t lADCData1[knADCData1];
+ UInt_t lADCData2[knADCData2];
+ UInt_t lADCData3[knADCData1];
+ UInt_t lADCData4[knADCData2];
+ //
+ UInt_t lADCEndBlock;
+
+ // load the digits
+ fLoader->LoadDigits("read");
+ AliZDCDigit digit;
+ AliZDCDigit* pdigit = &digit;
+ TTree* treeD = fLoader->TreeD();
+ if(!treeD) return;
+ treeD->SetBranchAddress("ZDC", &pdigit);
+ //printf("\t AliZDC::Digits2Raw -> TreeD has %d entries\n",(Int_t) treeD->GetEntries());
+
+ // Fill data array
+ // ADC header
+ UInt_t lADCHeaderGEO = 0;
+ UInt_t lADCHeaderCRATE = 0;
+ UInt_t lADCHeaderCNT1 = knADCData1;
+ UInt_t lADCHeaderCNT2 = knADCData2;
+
+ lADCHeader1 = lADCHeaderGEO << 27 | 0x1 << 25 | lADCHeaderCRATE << 16 |
+ lADCHeaderCNT1 << 8 ;
+ lADCHeader2 = lADCHeaderGEO << 27 | 0x1 << 25 | lADCHeaderCRATE << 16 |
+ lADCHeaderCNT2 << 8 ;
+
+ // ADC data word
+ UInt_t lADCDataGEO = lADCHeaderGEO;
+ //
+ UInt_t lADCDataValue1[knADCData1];
+ UInt_t lADCDataValue2[knADCData2];
+ UInt_t lADCDataValue3[knADCData1];
+ UInt_t lADCDataValue4[knADCData2];
+ //
+ UInt_t lADCDataOvFlw1[knADCData1];
+ UInt_t lADCDataOvFlw2[knADCData2];
+ UInt_t lADCDataOvFlw3[knADCData1];
+ UInt_t lADCDataOvFlw4[knADCData2];
+ //
+ for(Int_t i=0; i<knADCData1 ; i++){
+ lADCDataValue1[i] = 0;
+ lADCDataOvFlw1[i] = 0;
+ lADCDataValue3[i] = 0;
+ lADCDataOvFlw3[i] = 0;
+ }
+ for(Int_t i=0; i<knADCData2 ; i++){
+ lADCDataValue2[i] = 0;
+ lADCDataOvFlw2[i] = 0;
+ lADCDataValue4[i] = 0;
+ lADCDataOvFlw4[i] = 0;
+ }
+ //
+ UInt_t lADCDataChannel = 0;
+
+ // loop over digits
+ for(Int_t iDigit=0; iDigit<treeD->GetEntries(); iDigit++){
+ treeD->GetEntry(iDigit);
+ if(!pdigit) continue;
+ //digit.Print("");
+
+ // *** ADC data
+ Int_t index=0;
+ if(digit.GetSector(1)!=5){ // ZDC signal channels
+ // *** ADC1 (ZN1, ZP1, ZEM1,2) or ADC3 (ZN1, ZP1, ZEM1,2 o.o.t.)
+ if(digit.GetSector(0)==1 || digit.GetSector(0)==2 || digit.GetSector(0)==3){
+ if(digit.GetSector(0)==1 || digit.GetSector(0)==2){
+ index = (digit.GetSector(0)-1) + 4*digit.GetSector(1); // ZN1 or ZP1
+ lADCDataChannel = 8*(digit.GetSector(0)-1) + digit.GetSector(1);
+ }
+ else if(digit.GetSector(0)==3){ // ZEM 1,2
+ index = 20 + (digit.GetSector(1)-1);
+ lADCDataChannel = 5 + 8*(digit.GetSector(1)-1);
+ }
+ //
+ /*printf("\t AliZDC::Digits2Raw -> idig%d det %d quad %d index %d, ADCch %d ADCVal[%d, %d]\n",
+ iDigit,digit.GetSector(0),digit.GetSector(1),index,lADCDataChannel,
+ digit.GetADCValue(0),digit.GetADCValue(1));// Ch. debug
+ */
+ //
+ if(iDigit<knADCData1){ // *** In-time signals
+ lADCDataValue1[index] = digit.GetADCValue(0); // High gain ADC ch.
+ if(lADCDataValue1[index] > 2047) lADCDataOvFlw1[index] = 1;
+ lADCDataValue1[index+2] = digit.GetADCValue(1); // Low gain ADC ch.
+ if(lADCDataValue1[index+2] > 2047) lADCDataOvFlw1[index+2] = 1;
+
+ lADCData1[index] = lADCDataGEO << 27 | 0x1 << 24 | lADCDataChannel << 17 |
+ lADCDataOvFlw1[index] << 12 | (lADCDataValue1[index] & 0xfff);
+ lADCData1[index+2] = lADCDataGEO << 27 | 0x1 << 24 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlw1[index+2] << 12 | (lADCDataValue1[index+2] & 0xfff);
+ }
+ else{ // *** Out-of-time signals
+ lADCDataValue3[index] = digit.GetADCValue(0); // High gain ADC ch.
+ if(lADCDataValue3[index] > 2047) lADCDataOvFlw3[index] = 1;
+ lADCDataValue3[index+2] = digit.GetADCValue(1); // Low gain ADC ch.
+ if(lADCDataValue3[index+2] > 2047) lADCDataOvFlw3[index+2] = 1;
+
+ lADCData3[index] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlw3[index] << 12 | (lADCDataValue3[index] & 0xfff);
+ lADCData3[index+2] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlw3[index+2] << 12 | (lADCDataValue3[index+2] & 0xfff);
+ }
+ }
+ // *** ADC2 (ZN2, ZP2) or ADC4 (ZN2, ZP2 o.o.t.)
+ else if(digit.GetSector(0)==4 || digit.GetSector(0)==5){
+ index = (digit.GetSector(0)-4) + 4*digit.GetSector(1); // ZN2 or ZP2
+ lADCDataChannel = 8*(digit.GetSector(0)-4) + digit.GetSector(1);
+ //
+ /*printf("\t AliZDC::Digits2Raw -> idig%d det %d quad %d index %d, ADCch %d ADCVal[%d, %d]\n",
+ iDigit,digit.GetSector(0),digit.GetSector(1),index,lADCDataChannel,
+ digit.GetADCValue(0),digit.GetADCValue(1));// Ch. debug
+ */
+ //
+ if(iDigit<knADCData2){ // *** In-time signals
+ lADCDataValue2[index] = digit.GetADCValue(0);
+ if(lADCDataValue2[index] > 2047) lADCDataOvFlw2[index] = 1;
+ lADCDataValue2[index+2] = digit.GetADCValue(1);
+ if(lADCDataValue2[index+2] > 2047) lADCDataOvFlw2[index+2] = 1;
+ //
+ lADCData2[index] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlw2[index] << 12 | (lADCDataValue2[index] & 0xfff);
+ lADCData2[index+2] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlw2[index+2] << 12 | (lADCDataValue2[index+2] & 0xfff);
+ }
+ else{ // *** Out-of-time signals
+ lADCDataValue4[index] = digit.GetADCValue(0);
+ if(lADCDataValue4[index] > 2047) lADCDataOvFlw4[index] = 1;
+ lADCDataValue4[index+2] = digit.GetADCValue(1);
+ if(lADCDataValue4[index+2] > 2047) lADCDataOvFlw4[index+2] = 1;
+ //
+ lADCData4[index] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlw4[index] << 12 | (lADCDataValue4[index] & 0xfff);
+ lADCData4[index+2] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlw4[index+2] << 12 | (lADCDataValue4[index+2] & 0xfff);
+ }
+ }
+ }
+ // *** ADC2 (Reference PTMs) or ADC4 (Reference PTMs o.o.t.)
+ else if(digit.GetSector(1)==5){
+ index = 20 + (digit.GetSector(0)-1)/3;
+ lADCDataChannel = 5 + 8*(digit.GetSector(0)-1)/3;
+ //
+ /*printf("\t AliZDC::Digits2Raw -> idig%d det %d quad %d index %d, ADCch %d ADCVal[%d, %d]\n",
+ iDigit,digit.GetSector(0),digit.GetSector(1),index,lADCDataChannel,
+ digit.GetADCValue(0),digit.GetADCValue(1));// Ch. debug
+ */
+ //
+ if(iDigit<knADCData2){ // *** In-time signals
+ lADCDataValue2[index] = digit.GetADCValue(0);
+ if(lADCDataValue2[index] > 2047) lADCDataOvFlw2[index] = 1;
+ lADCDataValue2[index+2] = digit.GetADCValue(1);
+ if(lADCDataValue2[index+2] > 2047) lADCDataOvFlw2[index+2] = 1;
+ //
+ lADCData2[index] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlw2[index] << 12 | (lADCDataValue2[index] & 0xfff);
+ lADCData2[index+2] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlw2[index+2] << 12 | (lADCDataValue2[index+2] & 0xfff);
+ }
+ else{ // *** Out-of-time signals
+ lADCDataValue4[index] = digit.GetADCValue(0);
+ if(lADCDataValue4[index] > 2047) lADCDataOvFlw4[index] = 1;
+ lADCDataValue4[index+2] = digit.GetADCValue(1);
+ if(lADCDataValue4[index+2] > 2047) lADCDataOvFlw4[index+2] = 1;
+ //
+ lADCData4[index] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlw4[index] << 12 | (lADCDataValue4[index] & 0xfff);
+ lADCData4[index+2] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlw4[index+2] << 12 | (lADCDataValue4[index+2] & 0xfff);
+ }
+
+ }
+ if((index<0) || (index>23)) {
+ Error("Digits2Raw", "sector[0] = %d, sector[1] = %d",
+ digit.GetSector(0), digit.GetSector(1));
+ continue;
+ }
+
+
+ }
+ //
+ /*
+ for(Int_t i=0;i<knADCData1;i++) printf("\t ADCData1[%d] = %x\n",i,lADCData1[i]);
+ for(Int_t i=0;i<knADCData2;i++) printf("\t ADCData2[%d] = %x\n",i,lADCData2[i]);
+ for(Int_t i=0;i<knADCData1;i++) printf("\t ADCData3[%d] = %x\n",i,lADCData3[i]);
+ for(Int_t i=0;i<knADCData2;i++) printf("\t ADCData4[%d] = %x\n",i,lADCData4[i]);
+ */
+
+ // End of Block
+ UInt_t lADCEndBlockGEO = lADCHeaderGEO;
+ UInt_t lADCEndBlockEvCount = gAlice->GetEventNrInRun();
+ //
+ lADCEndBlock = lADCEndBlockGEO << 27 | 0x1 << 26 | lADCEndBlockEvCount;
+ //printf("\t AliZDC::Digits2Raw -> ADCEndBlock = %d\n",lADCEndBlock);
+
+
+ // open the output file
+ char fileName[30];
+ strcpy(fileName,AliDAQ::DdlFileName("ZDC",0));
+
+ AliFstream* file = new AliFstream(fileName);
+
+ // write the DDL data header
+ AliRawDataHeaderSim header;
+ header.fSize = sizeof(header) +
+ sizeof(lADCHeader1) + sizeof(lADCData1) + sizeof(lADCEndBlock) +
+ sizeof(lADCHeader2) + sizeof(lADCData2) + sizeof(lADCEndBlock) +
+ sizeof(lADCHeader1) + sizeof(lADCData3) + sizeof(lADCEndBlock) +
+ sizeof(lADCHeader2) + sizeof(lADCData4) + sizeof(lADCEndBlock);