+ // *** Fill data array
+ // ** ADC header
+ UInt_t lADCHeaderGEO1 = 0;
+ UInt_t lADCHeaderGEO2 = 1;
+ UInt_t lADCHeaderGEO3 = 2;
+ UInt_t lADCHeaderGEO4 = 3;
+ UInt_t lADCHeaderCRATE = 0;
+ UInt_t lADCHeaderCNT1 = knADCData1;
+ UInt_t lADCHeaderCNT2 = knADCData2;
+ UInt_t lADCHeaderCNT3 = knADCData3;
+ UInt_t lADCHeaderCNT4 = knADCData4;
+
+ lADCHeader1 = lADCHeaderGEO1 << 27 | 0x1 << 25 | lADCHeaderCRATE << 16 |
+ lADCHeaderCNT1 << 8 ;
+ lADCHeader2 = lADCHeaderGEO2 << 27 | 0x1 << 25 | lADCHeaderCRATE << 16 |
+ lADCHeaderCNT2 << 8 ;
+ lADCHeader3 = lADCHeaderGEO3 << 27 | 0x1 << 25 | lADCHeaderCRATE << 16 |
+ lADCHeaderCNT3 << 8 ;
+ lADCHeader4 = lADCHeaderGEO4 << 27 | 0x1 << 25 | lADCHeaderCRATE << 16 |
+ lADCHeaderCNT4 << 8 ;
+
+ // ** ADC data word
+ UInt_t lADCDataGEO = 0;
+ //
+ UInt_t lADCDataValue1[2*knADCData1];
+ UInt_t lADCDataValue2[2*knADCData2];
+ UInt_t lADCDataValue3[2*knADCData3];
+ UInt_t lADCDataValue4[2*knADCData4];
+ //
+ UInt_t lADCDataOvFlwHG = 0;
+ UInt_t lADCDataOvFlwLG = 0;
+ //
+ for(Int_t i=0; i<2*knADCData1 ; i++) lADCDataValue1[i] = 0;
+ for(Int_t i=0; i<2*knADCData2 ; i++) lADCDataValue2[i] = 0;
+ for(Int_t i=0; i<2*knADCData3 ; i++) lADCDataValue3[i] = 0;
+ for(Int_t i=0; i<2*knADCData4 ; i++) lADCDataValue4[i] = 0;
+ //
+ UInt_t lADCDataChannel = 0;
+
+ Int_t indADC0=0, indADC1=0, indADC2=0, indADC3=0;
+
+ // loop over digits
+ for(Int_t iDigit=0; iDigit<(Int_t) (treeD->GetEntries()); iDigit++){
+ treeD->GetEntry(iDigit);
+ if(!pdigit) continue;
+ //digit.Print("");
+
+ // *** ADC data
+ // Scan of the map to assign the correct ADC module-channel
+ for(Int_t k=0; k<nCh; k++){
+ if(iDigit<knADCData1+knADCData2){
+ if(digit.GetSector(0)==mapADC[k][2] && digit.GetSector(1)==mapADC[k][3]){
+ lADCDataGEO = (UInt_t) mapADC[k][0];
+ lADCDataChannel = (UInt_t) mapADC[k][1];
+ break;
+ }
+ }
+ else{
+ if(digit.GetSector(0)==mapADC[k][2] && digit.GetSector(1)==mapADC[k][3]){
+ lADCDataGEO = (UInt_t) mapADC[k][0];
+ lADCDataChannel = (UInt_t) mapADC[k][1];
+ if(k>knADCData1+knADCData2) break;
+ }
+ }
+ }
+ // Ch. debug
+ //printf("iDigit %d det %d sec %d -> lADCDataGEO %d lADCDataChannel %d\n",
+ // iDigit,digit.GetSector(0),digit.GetSector(1),lADCDataGEO,lADCDataChannel);
+
+ if(lADCDataGEO==0){
+ if(indADC0>=knADCData1){
+ AliWarning(" Problem with digit index 4 ADC0\n");
+ return;
+ }
+ Int_t indLG = indADC0+knADCData1;
+ // High gain ADC ch.
+ if(digit.GetADCValue(0) > 2047) lADCDataOvFlwHG = 1;
+ lADCDataValue1[indADC0] = digit.GetADCValue(0);
+ lADCData1[indADC0] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlwHG << 12 | (lADCDataValue1[indADC0] & 0xfff);
+ // Low gain ADC ch.
+ if(digit.GetADCValue(1) > 2047) lADCDataOvFlwLG = 1;
+ lADCDataValue1[indLG] = digit.GetADCValue(1);
+ lADCData1[indLG] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlwLG << 12 | (lADCDataValue1[indLG] & 0xfff);
+ // Ch. debug
+ //printf(" lADCDataGEO %d ADCdataHG[%d] %d ADCdataLG[%d] %d\n",
+ // lADCDataGEO,indADC0,lADCDataValue1[indADC0],indLG,lADCDataValue1[indLG]);
+
+ indADC0++;
+ }
+ else if(lADCDataGEO==1){
+ if(indADC1>=knADCData2){
+ AliWarning(" Problem with digit index 4 ADC1\n");
+ return;
+ }
+ Int_t indLG = indADC1+knADCData2;
+ // High gain ADC ch.
+ if(digit.GetADCValue(0) > 2047) lADCDataOvFlwHG = 1;
+ lADCDataValue2[indADC1] = digit.GetADCValue(0);
+ lADCData2[indADC1] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlwHG << 12 | (lADCDataValue2[indADC1] & 0xfff);
+ // Low gain ADC ch.
+ if(digit.GetADCValue(1) > 2047) lADCDataOvFlwLG = 1;
+ lADCDataValue2[indLG] = digit.GetADCValue(1);
+ lADCData2[indLG] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlwLG << 12 | (lADCDataValue2[indLG] & 0xfff);
+ // Ch. debug
+ //printf(" lADCDataGEO %d ADCdataHG[%d] %d ADCdataLG[%d] %d\n",
+ // lADCDataGEO,indADC1,lADCDataValue2[indADC1],indLG,lADCDataValue2[indLG]);
+
+ indADC1++;
+ }
+ else if(lADCDataGEO==2){
+ if(indADC2>=knADCData3){
+ AliWarning(" Problem with digit index 4 ADC2\n");
+ return;
+ }
+ Int_t indLG = indADC2+knADCData3;
+ // High gain ADC ch.
+ if(digit.GetADCValue(0) > 2047) lADCDataOvFlwHG = 1;
+ lADCDataValue3[indADC1] = digit.GetADCValue(0);
+ lADCData3[indADC1] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlwHG << 12 | (lADCDataValue3[indADC2] & 0xfff);
+ // Low gain ADC ch.
+ if(digit.GetADCValue(1) > 2047) lADCDataOvFlwLG = 1;
+ lADCDataValue3[indLG] = digit.GetADCValue(1);
+ lADCData3[indLG] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlwLG << 12 | (lADCDataValue3[indLG] & 0xfff);
+ // Ch. debug
+ //printf(" lADCDataGEO %d ADCdataHG[%d] %d ADCdataLG[%d] %d\n",
+ // lADCDataGEO,indADC2,lADCDataValue3[indADC2],indLG,lADCDataValue3[indLG]);
+
+ indADC2++;
+ }
+ else if(lADCDataGEO==3){
+ if(indADC3>=knADCData4){
+ AliWarning(" Problem with digit index 4 ADC2\n");
+ return;
+ }
+ Int_t indLG = indADC3+knADCData4;
+ // High gain ADC ch.
+ if(digit.GetADCValue(0) > 2047) lADCDataOvFlwHG = 1;
+ lADCDataValue4[indADC3] = digit.GetADCValue(0);
+ lADCData4[indADC3] = lADCDataGEO << 27 | lADCDataChannel << 17 |
+ lADCDataOvFlwHG << 12 | (lADCDataValue4[indADC3] & 0xfff);
+ // Low gain ADC ch.
+ if(digit.GetADCValue(1) > 2047) lADCDataOvFlwLG = 1;
+ lADCDataValue4[indLG] = digit.GetADCValue(1);
+ lADCData4[indLG] = lADCDataGEO << 27 | lADCDataChannel << 17 | 0x1 << 16 |
+ lADCDataOvFlwLG << 12 | (lADCDataValue4[indLG] & 0xfff);
+ // Ch. debug
+ //printf(" lADCDataGEO %d ADCdataHG[%d] %d ADCdataLG[%d] %d\n",
+ // lADCDataGEO,indADC3,lADCDataValue4[indADC3],indLG,lADCDataValue4[indLG]);
+
+ indADC3++;
+ }