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 fRawReader(stream.fRawReader),
55 fRawADC(stream.GetADCRaw()),
56 fADCModule(stream.GetADCModule()),
57 fADCChannel(stream.GetADCChannel()),
58 fADCValue(stream.GetADCValue()),
59 fADCGain(stream.GetADCGain()),
60 fIsADCDataWord(stream.IsADCDataWord())
63 for(Int_t j=0; j<2; j++) fSector[j] = stream.GetSector(j);
67 //_____________________________________________________________________________
68 AliZDCRawStream& AliZDCRawStream::operator = (const AliZDCRawStream&
71 // Assignment operator
72 Fatal("operator =", "assignment operator not implemented");
76 //_____________________________________________________________________________
77 AliZDCRawStream::~AliZDCRawStream()
84 //_____________________________________________________________________________
85 Bool_t AliZDCRawStream::Next()
87 // Read the next raw digit
88 // Returns kFALSE if there is no digit left
90 if(!fRawReader->ReadNextInt((UInt_t&) fRawADC)) return kFALSE;
91 fIsADCDataWord = kFALSE;
94 if((fRawADC == 0xe52b6300) || (fRawADC == 0xe52c0300) || (fRawADC == 0xffffffff) ||
95 (fRawADC == 0xdeadface) || (fRawADC == 0xdeadbeef)){
96 //printf(" This is a DARC header!!!\n");
99 else if(fRawADC == 0xcafefade){
100 //printf(" End of data!\n");
104 Bool_t firstADCHeader = kTRUE;
106 if((fRawADC & 0x6000000) == 0x6000000){ // Not valid datum BEFORE valid data!
107 firstADCHeader = kFALSE;
108 //printf(" AliZDCRawStream -> Not valid datum in ADC module!!! %d\n",fADCModule);
110 else if((fRawADC & 0x4000000) == 0x4000000){ // ADC EOB
111 //printf(" AliZDCRawStream -> ADC %d End Of Block - Event no. %d\n\n",fADCModule, (fRawADC & 0xffffff));
112 fSector[0] = fSector[1] = 99;
114 else if((fRawADC & 0x2000000) == 0x2000000){ // ADC Header
115 //printf(" fRawADC %x\n",fRawADC);
116 // Reading the GEO address to determine ADC module
117 fADCModule = (fRawADC & 0xf8000000) >> 27;
119 // If the not valid datum isn't followed by the 1st ADC header
120 // the event is corrupted (i.e., 2 gates arrived before trigger)
121 if(fADCModule == 1) firstADCHeader = kTRUE;
122 //if(fADCModule == 1) printf(" fADCModule %d, firstADCHeader %d\n",fADCModule,firstADCHeader);
123 //printf(" AliZDCRawStream -> HEADER: ADC mod. %d contains %d data words \n",fADCModule,((fRawADC & 0x3f00) >> 8));
125 else{ // ADC data word
126 fIsADCDataWord = kTRUE;
127 //printf(" fRawADC = %x firstADCHeader %d\n",fRawADC,firstADCHeader);
129 if(!(fRawADC & 0x1000) && !(fRawADC & 0x2000) && firstADCHeader){ // Valid ADC data
130 fADCGain = (fRawADC & 0x10000) >> 16;
131 fADCValue = (fRawADC & 0xfff);
133 fADCChannel = (fRawADC & 0x1e0000) >> 17;
135 // If raw data corresponds to an ADC ch. without physical signal
136 // fsector[0] = fSector[1] = -1
137 for(Int_t i=0; i<2; i++)fSector[i] = -1;
138 //if(fSector[0]!=-1) printf(" \t \t ADC Data Word: ADC ch. %d",fADCChannel);
140 if(fADCModule==1 || fADCModule==3){ //1st & 3rd ADC modules
141 if(fADCChannel >= 0 && fADCChannel <= 4){
142 fSector[0] = 1; // ZN1
143 fSector[1] = fADCChannel;
145 else if(fADCChannel >= 8 && fADCChannel <= 12){
146 fSector[0] = 2; // ZP1
147 fSector[1] = fADCChannel-8;
149 else if(fADCChannel == 5 || fADCChannel == 13){
150 fSector[0] = 3; // ZEM 1,2
151 fSector[1] = ((fADCChannel-5)/8)+1;
154 else if(fADCModule==2 || fADCModule==4){ //2nd & 4rth ADC modules
155 if(fADCChannel >= 0 && fADCChannel <= 4){
156 fSector[0] = 4; // ZN2
157 fSector[1] = fADCChannel;
159 else if(fADCChannel >= 8 && fADCChannel <= 12){
160 fSector[0] = 5; // ZP2
161 fSector[1] = fADCChannel-8;
163 else if(fADCChannel == 5 || fADCChannel == 13){
164 fSector[0] = (fADCChannel-5)*3/8+1; // PM Ref 1,2
169 AliWarning(" AliZDCRawStream -> No valid ADC module\n");
170 fRawReader->AddMajorErrorLog(kInvalidADCModule);
172 //printf(" ADC %d ch %d range %d -> det %d quad %d value %x\n",
173 // fADCModule, fADCChannel, fADCGain, fSector[0], fSector[1], fADCValue);//Chiara debugging
176 else if(fRawADC & 0x1000){
177 //printf(" ADC %d ch %d range %d overflow\n",fADCModule, (fRawADC & 0x1e0000) >> 17, fADCGain); // Overflow
179 else if(fRawADC & 0x2000){
180 //printf(" ADC %d ch %d range %d underflow\n",fADCModule, (fRawADC & 0x1e0000) >> 17, fADCGain); // Underflow