fCharge(0x0),
fPad(0x0),
fDDLNumber(-1),
+ fnDDLInStream(0x0),
+ fnDDLOutStream(0x0),
fLDCNumber( 0),
fTimeStamp( 0),
fRawReader(rawReader),
fNumOfErr[i] = new Int_t [kSumErr];
}
+ fnDDLInStream=new Int_t[kNDDL];
+ fnDDLOutStream=new Int_t[kNDDL];
+ for(Int_t iddl=0;iddl<kNDDL;iddl++) { fnDDLInStream[iddl]=-1;fnDDLOutStream[iddl]=-1;}
+
for(Int_t iddl=0;iddl<kNDDL;iddl++)
for(Int_t ierr=0; ierr < kSumErr; ierr++) fNumOfErr[iddl][ierr]=0; //reset errors
fRawReader->Reset();
fCharge(0x0),
fPad(0x0),
fDDLNumber(-1),
+ fnDDLInStream(0x0),
+ fnDDLOutStream(0x0),
fLDCNumber( 0),
fTimeStamp( 0),
fRawReader(0x0),
for(Int_t i=0;i<kNDDL;i++) {
fNumOfErr[i] = new Int_t [kSumErr];
}
+ fnDDLInStream=new Int_t[kNDDL];
+ fnDDLOutStream=new Int_t[kNDDL];
+ for(Int_t iddl=0;iddl<kNDDL;iddl++) { fnDDLInStream[iddl]=-1;fnDDLOutStream[iddl]=-1;}
+
+
for(Int_t iddl=0;iddl<kNDDL;iddl++)
for(Int_t ierr=0; ierr < kSumErr; ierr++) fNumOfErr[iddl][ierr]=0; //reset errors
Int_t rawDataSize=0;
if(fRawReader->GetType() == 7) { //New: Select Physics events, Old: Raw data size is not 0 and not 47148 (pedestal)
fDDLNumber = fRawReader->GetDDLID();
+ //if(fnDDLInStream[fDDLNumber]==-1) { fnDDLInStream[fDDLNumber]=0; fnDDLOutStream[fDDLNumber]=0; }//Printf("Setting DDL number: %d to fnDDLInStream[]= %d",fDDLNumber,fnDDLInStream[fDDLNumber]);}
+ fnDDLInStream[fDDLNumber]=1; fnDDLOutStream[fDDLNumber]=0;
+// fnDDLInStream[fDDLNumber]++;
+
fLDCNumber = fRawReader->GetLDCId();
fTimeStamp = fRawReader->GetTimestamp();
-
+
if(stDeb) Printf("DDL %i started to be decoded!",fDDLNumber);
rawDataSize=fRawReader->GetDataSize()/4;
InitVars(rawDataSize);
}
// stDeb=kFALSE;
}
-// return status;
- return kTRUE;
+ if(status==kTRUE) {fnDDLOutStream[fDDLNumber]++; }//Printf("fnDDLOutStream[%d]=%d",fDDLNumber,fnDDLOutStream[fDDLNumber]); } //Count the number of events when the DDL was succesfully decoded
+ return status;
+ // return kTRUE;
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPIDRawStream::InitVars(Int_t n)
for(Int_t i=0;i<kSumErr;i++) delete [] fNumOfErr[i];
delete [] fNumOfErr;
-
-
+ if(fnDDLInStream) { delete [] fnDDLInStream; fnDDLInStream = 0x0; }
+
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Bool_t AliHMPIDRawStream::ReadHMPIDRawData()
if(!CheckPad(pad)) continue;
Int_t charge = fWord & 0xfff;
fPad[fNPads] = GetPad(fDDLNumber,row,dilogic,pad);
- fCharge[fNPads] = charge;
+ fCharge[fNPads] = charge;
fNPads++;
- }
+
+ if(charge==0)
+ {
+ if(stDeb){AliWarning(Form("If PEDESTAL run -> WARNING: ZERO charge is read from DDL: %d row: %d dil: %d pad: %d",fDDLNumber,row,dilogic,pad));}
+ fNumOfErr[fDDLNumber][kPedQZero]++;
+ }
+
+ }//iDil
cntDilogic -= cnt;
return kTRUE;
static inline Int_t GetNDDL() { return kNDDL;} //return the number of max # of DDLs
static inline Int_t GetNErrors() { return kSumErr;} //return the number of max # of Error Types
- Int_t GetNPads() const { return fNPads;}
- Int_t* GetPadArray() const { return fPad;}
- Int_t* GetChargeArray() const { return fCharge;}
+ Int_t GetNPads() const { return fNPads;} //Get number of pads present in the stream
+ Int_t* GetPadArray() const { return fPad;} //Get pad array from stream decoded
+ Int_t* GetChargeArray() const { return fCharge;} //Get the charge of the pads from dedcoded stream
+ Int_t* GetnDDLInStream() const { return fnDDLInStream;} //Get the DDL input check array
+ Int_t* GetnDDLOutStream() const { return fnDDLOutStream;} //Get the DDL output check array
Int_t Pc ( Int_t ddl,Int_t row,Int_t dil,Int_t pad ) {return AliHMPIDParam::A2P(GetPad(ddl,row,dil,pad));} //PC position number
Int_t PadPcX ( Int_t ddl,Int_t row,Int_t dil,Int_t pad ) {return AliHMPIDParam::A2X(GetPad(ddl,row,dil,pad));} //pad pc x # 0..79
Int_t PadPcY ( Int_t ddl,Int_t row,Int_t dil,Int_t pad ) {return AliHMPIDParam::A2Y(GetPad(ddl,row,dil,pad));} //pad pc y # 0..47
enum EHMPIDRawStreamError { kRawDataSizeErr = 0, kRowMarkerErr = 1, kWrongRowErr = 2, kWrongDilogicErr = 3,
kWrongPadErr = 4, kEoEFlagErr = 5, kEoESizeErr = 6, kEoEDILOGICErr = 7,
kEoERowErr = 8, kBadSegWordErr = 9, kWrongSegErr = 10, kRowMarkerSizeErr = 11,
- kSumErr = 12 //This is always the last one, to retreive the number of errors
+ kPedQZero =12, kSumErr = 13 //This is always the last one, to retreive the number of errors
}; //Always check the updated list of names in the .cxx file for print-out!
enum {
Int_t *fCharge; // Array for charge values for all channels in one DDL
Int_t *fPad; // Array for abs pad values for all channels in one DDL
Int_t fDDLNumber; // index of current DDL number
+ Int_t *fnDDLInStream; // if the DDL is in the raw data
+ Int_t *fnDDLOutStream; // if the DDL is in the raw data
UInt_t fLDCNumber; // index of current LDC number
UInt_t fTimeStamp; // TimeStamp
AliRawReader *fRawReader; // object for reading the raw data
UInt_t fWord; // current position in fData
Bool_t fZeroSup; // set if zero suppression is applied
Int_t *fPos; // for debug purposes
- Int_t fiPos; // counter for debug
-
+ Int_t fiPos; // counter for debug
ClassDef(AliHMPIDRawStream, 2) // base class for reading HMPID raw digits
};
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Int_t a2y[6]={3,2,4,1,5,0}; //pady for a given padress (for single DILOGIC chip)
Int_t ch=ddl/2;
- Int_t tmp=(row-1)/8; Int_t pc=(ddl%2)? 5-2*tmp:2*tmp;
- Int_t px=(dil-1)*8+pad/6;
- tmp=(ddl%2)?(24-row):row-1;
- Int_t py=6*(tmp%8)+a2y[pad%6];
-
+ Int_t tmp=(24-row)/8;
+ Int_t pc=(ddl%2)?5-2*tmp:2*tmp;
+ Int_t px=(dil-1)*8+pad/6;
+
+ tmp=(ddl%2)?row-1:(24-row);
+ Int_t py=6*(tmp%8)+a2y[pad%6];
+
return AliHMPIDParam::Abs(ch,pc,px,py);
}//GetPad()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
for(Int_t iDig=0;iDig<pDigCh->GetEntriesFast();iDig++){//digits loop
AliHMPIDDigit *pDig1=(AliHMPIDDigit*)pDigCh->At(iDig);
- pDig1->Raw(w32,ddl,r,d,a);
+ pDig1->Raw(w32,ddl,r,d,a); //??????????
isDigThere[ddl][r][d][a]=iDig;
}
Char_t *eName[]={ "kRawDataSizeErr", "kRowMarkerErr" , "kWrongRowErr" , "kWrongDilogicErr",
"kWrongPadErr" , "kEoEFlagErr" , "kEoESizeErr" , "kEoEDILOGICErr",
"kEoERowErr" , "kBadSegWordErr", "kWrongSegErr" , "kRowMarkerSizeErr",
- "kSumErr" };
+ "kPedQZero" , "kSumErr" };
Char_t *eNoErr="NotDefinedErrorType";
if(eType<0 || eType>kSumErr) return eNoErr;
else return eName[eType];