1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
20 // This class provides access to ZDC digits in raw data. //
22 // It loops over all ZDC digits in the raw data given by the AliRawReader. //
23 // The Next method goes to the next digit. If there are no digits left //
24 // it returns kFALSE. //
25 // Getters provide information about the current digit. //
27 ///////////////////////////////////////////////////////////////////////////////
29 #include "AliZDCRawStream.h"
30 #include "AliRawReader.h"
33 ClassImp(AliZDCRawStream)
36 //_____________________________________________________________________________
37 AliZDCRawStream::AliZDCRawStream(AliRawReader* rawReader) :
38 fRawReader(rawReader),
44 fIsADCDataWord(kFALSE)
46 // Create an object to read ZDC raw digits
48 fRawReader->Select("ZDC");
51 //_____________________________________________________________________________
52 AliZDCRawStream::AliZDCRawStream(const AliZDCRawStream& stream) :
54 fRawADC(stream.GetADCRaw()),
55 fADCModule(stream.GetADCModule()),
56 fADCChannel(stream.GetADCChannel()),
57 fADCValue(stream.GetADCValue()),
58 fADCGain(stream.GetADCGain()),
59 fIsADCDataWord(stream.IsADCDataWord())
62 for(Int_t j=0; j<2; j++) fSector[j] = stream.GetSector(j);
66 //_____________________________________________________________________________
67 AliZDCRawStream& AliZDCRawStream::operator = (const AliZDCRawStream&
70 // Assignment operator
71 Fatal("operator =", "assignment operator not implemented");
75 //_____________________________________________________________________________
76 AliZDCRawStream::~AliZDCRawStream()
83 //_____________________________________________________________________________
84 Bool_t AliZDCRawStream::Next()
86 // Read the next raw digit
87 // Returns kFALSE if there is no digit left
89 if(!fRawReader->ReadNextInt((UInt_t&) fRawADC)) return kFALSE;
90 fIsADCDataWord = kFALSE;
93 if((fRawADC == 0xe52b6300) || (fRawADC == 0xe52c0300) || (fRawADC == 0xffffffff) ||
94 (fRawADC == 0xdeadface) || (fRawADC == 0xdeadbeef)){
95 //printf(" This is a DARC header!!!\n");
98 else if(fRawADC == 0xcafefade){
99 //printf(" End of data!\n");
103 Bool_t firstADCHeader = kTRUE;
105 if((fRawADC & 0x6000000) == 0x6000000){ // Not valid datum BEFORE valid data!
106 firstADCHeader = kFALSE;
107 //printf(" AliZDCRawStream -> Not valid datum in ADC module!!! %d\n",fADCModule);
109 else if((fRawADC & 0x4000000) == 0x4000000){ // ADC EOB
110 //printf(" AliZDCRawStream -> ADC %d End Of Block - Event no. %d\n\n",fADCModule, (fRawADC & 0xffffff));
111 fSector[0] = fSector[1] = 99;
113 else if((fRawADC & 0x2000000) == 0x2000000){ // ADC Header
114 //printf(" fRawADC %x\n",fRawADC);
115 // Reading the GEO address to determine ADC module
116 fADCModule = (fRawADC & 0xf8000000) >> 27;
118 // If the not valid datum isn't followed by the 1st ADC header
119 // the event is corrupted (i.e., 2 gates arrived before trigger)
120 if(fADCModule == 1) firstADCHeader = kTRUE;
121 //if(fADCModule == 1) printf(" fADCModule %d, firstADCHeader %d\n",fADCModule,firstADCHeader);
122 //printf(" AliZDCRawStream -> HEADER: ADC mod. %d contains %d data words \n",fADCModule,((fRawADC & 0x3f00) >> 8));
124 else{ // ADC data word
125 fIsADCDataWord = kTRUE;
126 //printf(" fRawADC = %x firstADCHeader %d\n",fRawADC,firstADCHeader);
128 if(!(fRawADC & 0x1000) && !(fRawADC & 0x2000) && firstADCHeader){ // Valid ADC data
129 fADCGain = (fRawADC & 0x10000) >> 16;
130 fADCValue = (fRawADC & 0xfff);
132 fADCChannel = (fRawADC & 0x1e0000) >> 17;
134 // If raw data corresponds to an ADC ch. without physical signal
135 // fsector[0] = fSector[1] = -1
136 for(Int_t i=0; i<2; i++)fSector[i] = -1;
137 //if(fSector[0]!=-1) printf(" \t \t ADC Data Word: ADC ch. %d",fADCChannel);
139 if(fADCModule==1 || fADCModule==3){ //1st & 3rd ADC modules
140 if(fADCChannel >= 0 && fADCChannel <= 4){
141 fSector[0] = 1; // ZN1
142 fSector[1] = fADCChannel;
144 else if(fADCChannel >= 8 && fADCChannel <= 12){
145 fSector[0] = 2; // ZP1
146 fSector[1] = fADCChannel-8;
148 else if(fADCChannel == 5 || fADCChannel == 13){
149 fSector[0] = 3; // ZEM 1,2
150 fSector[1] = ((fADCChannel-5)/8)+1;
153 else if(fADCModule==2 || fADCModule==4){ //2nd & 4rth ADC modules
154 if(fADCChannel >= 0 && fADCChannel <= 4){
155 fSector[0] = 4; // ZN2
156 fSector[1] = fADCChannel;
158 else if(fADCChannel >= 8 && fADCChannel <= 12){
159 fSector[0] = 5; // ZP2
160 fSector[1] = fADCChannel-8;
162 else if(fADCChannel == 5 || fADCChannel == 13){
163 fSector[0] = (fADCChannel-5)*3/8+1; // PM Ref 1,2
168 AliWarning(" AliZDCRawStream -> No valid ADC module\n");
169 fRawReader->AddMajorErrorLog(kInvalidADCModule);
171 //printf(" ADC %d ch %d range %d -> det %d quad %d value %x\n",
172 // fADCModule, fADCChannel, fADCGain, fSector[0], fSector[1], fADCValue);//Chiara debugging
175 else if(fRawADC & 0x1000){
176 //printf(" ADC %d ch %d range %d overflow\n",fADCModule, (fRawADC & 0x1e0000) >> 17, fADCGain); // Overflow
178 else if(fRawADC & 0x2000){
179 //printf(" ADC %d ch %d range %d underflow\n",fADCModule, (fRawADC & 0x1e0000) >> 17, fADCGain); // Underflow