**************************************************************************/
/*
+$Log$
+Revision 1.19 2007/06/22 11:37:47 cvetan
+Fixes in order to write correct raw-data on big-endian platforms (Marco)
+
+Revision 1.18 2007/05/21 13:26:19 decaro
+Correction on matching_window control and bug fixing (R.Preghenella)
+
+Revision 1.17 2007/05/10 09:29:34 hristov
+Last moment fixes and changes from v4-05-Release (Silvia)
+
+Revision 1.16 2007/05/03 09:07:22 decaro
+Double digit in the same TDC channel. Wrong sequence during the raw data writing in unpacked mode: solved
+
+Revision 1.15 2007/04/23 16:51:39 decaro
+Digits-to-raw_data conversion: correction for a more real description (A.De Caro, R.Preghenella)
+
+Revision 1.14 2007/03/28 10:50:33 decaro
+Rounding off problem in rawData coding/decoding: solved
+
+Revision 1.13 2007/02/20 15:57:00 decaro
+Raw data update: to read the TOF raw data defined in UNPACKED mode
+
+Revision 1.12 2006/08/22 13:29:42 arcelli
+removal of effective c++ warnings (C.Zampolli)
+
+Revision 1.11 2006/08/10 14:46:54 decaro
+TOF raw data format: updated version
+
+Revision 1.10.1 2006/06/28 A.De Caro
+ Update TOF raw data format
+ according to the final version
+ (see ALICE internal note in preparation
+ 'ALICE TOF raw data format')
+
Revision 0.02 2005/7/25 A.De Caro
Update number of bits allocated for time-of-flight
and 'charge' measurements
GetDigits
*/
-//////////////////////////////////////////////////////////////////
-//
-// This class contains the methods to create the Raw Data files
-// for the TOF detector starting from the Digits.
-// In this preliminary implementation, we defined the structure
-// of the ALICE-TOF raw data starting from the current format
-// for the TOF digits and the TOF raw data.
-//
-//////////////////////////////////////////////////////////////////
+////////////////////////////////////////////////////////////////////
+// //
+// This class contains the methods to create the Raw Data files //
+// for the TOF detector starting from the Digits. //
+// In this implementation, we defined the structure //
+// of the ALICE-TOF raw data (according to the //
+// ALICE technical note, in preparation) //
+// starting from the TOF digit format. //
+// //
+////////////////////////////////////////////////////////////////////
#include "Riostream.h"
#include "TBranch.h"
#include "TClonesArray.h"
#include "TMath.h"
+#include "TRandom.h"
#include "AliBitPacking.h"
+#include "AliDAQ.h"
#include "AliLog.h"
-#include "AliRawDataHeader.h"
+//#include "AliRawDataHeader.h"
+#include "AliRawDataHeaderSim.h"
+#include "AliFstream.h"
#include "AliTOFDDLRawData.h"
+#include "AliTOFDigitMap.h"
#include "AliTOFdigit.h"
#include "AliTOFGeometry.h"
#include "AliTOFRawStream.h"
+//#include "AliTOFCableLengthMap.h"
+
+extern TRandom *gRandom;
ClassImp(AliTOFDDLRawData)
-//----------------------------------------------------------------------------------------
-AliTOFDDLRawData::AliTOFDDLRawData()
+//---------------------------------------------------------------------------
+AliTOFDDLRawData::AliTOFDDLRawData():
+ TObject(),
+ fVerbose(0),
+ fIndex(-1),
+ fPackedAcquisition(kFALSE),
+ fFakeOrphaneProduction(kFALSE),
+ fMatchingWindow(8192),
+ fTOFdigitMap(new AliTOFDigitMap()),
+ fTOFdigitArray(0x0),
+ fWordsPerDRM(0),
+ fWordsPerTRM(0),
+ fWordsPerChain(0)
{
//Default constructor
- fIndex=-1;
- fVerbose=0;
- fTOFgeometry = 0;
}
+//----------------------------------------------------------------------------
+AliTOFDDLRawData::AliTOFDDLRawData(const AliTOFDDLRawData &source) :
+ TObject(source),
+ fVerbose(source.fVerbose),
+ fIndex(source.fIndex),
+ fPackedAcquisition(source.fPackedAcquisition),
+ fFakeOrphaneProduction(source.fFakeOrphaneProduction),
+ fMatchingWindow(source.fMatchingWindow),
+ fTOFdigitMap(source.fTOFdigitMap),
+ fTOFdigitArray(source.fTOFdigitArray),
+ fWordsPerDRM(source.fWordsPerDRM),
+ fWordsPerTRM(source.fWordsPerTRM),
+ fWordsPerChain(source.fWordsPerChain)
+ {
+ //Copy Constructor
+ return;
+}
+
+//----------------------------------------------------------------------------
+AliTOFDDLRawData& AliTOFDDLRawData::operator=(const AliTOFDDLRawData &source) {
+ //Assigment operator
-//----------------------------------------------------------------------------------------
-AliTOFDDLRawData::AliTOFDDLRawData(AliTOFGeometry *tofGeom)
+ if (this == &source)
+ return *this;
+
+ TObject::operator=(source);
+ fIndex=source.fIndex;
+ fVerbose=source.fVerbose;
+ fPackedAcquisition=source.fPackedAcquisition;
+ fFakeOrphaneProduction=source.fFakeOrphaneProduction;
+ fMatchingWindow=source.fMatchingWindow;
+ fTOFdigitMap=source.fTOFdigitMap;
+ fTOFdigitArray=source.fTOFdigitArray;
+ fWordsPerDRM=source.fWordsPerDRM;
+ fWordsPerTRM=source.fWordsPerTRM;
+ fWordsPerChain=source.fWordsPerChain;
+ return *this;
+}
+
+//----------------------------------------------------------------------------
+AliTOFDDLRawData::~AliTOFDDLRawData()
{
- //Constructor
- fIndex=-1;
- fVerbose=0;
- fTOFgeometry = tofGeom;
+ // dtr
+
+ delete fTOFdigitMap;
}
+//----------------------------------------------------------------------------
+Int_t AliTOFDDLRawData::RawDataTOF(TBranch* branch)
+{
+ //
+ // This method creates the Raw data files for TOF detector
+ //
-//----------------------------------------------------------------------------------------
+ const Int_t kSize = 5000; // max number of digits per DDL file times 2
+
+ UInt_t buf[kSize];
+
+ // To clear the digit indices map for each event
+ fTOFdigitMap->Clear();
+
+ fIndex = -1;
+
+ fTOFdigitArray = * (TClonesArray**) branch->GetAddress();
+
+ AliFstream* outfile; // logical name of the output file
+
+ AliRawDataHeaderSim header;
+
+ UInt_t sizeRawData = 0;
+
+ branch->GetEvent();
+
+ GetDigits();
+
+ //if (!fPackedAcquisition) fMatchingWindow = 2097152;//AdC
+
+ Int_t jj = -1;
+ Int_t iDDL = -1;
+ Int_t nDDL = -1;
+ Int_t nTRM = 0;
+ Int_t iChain = -1;
+
+ //loop over TOF DDL files
+ for (nDDL=0; nDDL<AliDAQ::NumberOfDdls("TOF"); nDDL++) {
+
+ char fileName[256]="";
+ strncpy(fileName,AliDAQ::DdlFileName("TOF",nDDL),255); //The name of the output file
+
+ outfile = new AliFstream(fileName);
+ iDDL = AliTOFRawStream::GetDDLnumberPerSector(nDDL);
+
+ // write Dummy DATA HEADER
+ UInt_t dataHeaderPosition = outfile->Tellp();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
+
+ // DRM section: trailer
+ MakeDRMtrailer(buf);
+
+ // loop on TRM number per DRM
+ for (nTRM=AliTOFGeometry::NTRM(); nTRM>=3; nTRM--) {
+
+ fWordsPerTRM = 0;
+
+ // the slot number 3 of the even DRM (i.e. right) doesn't contain TDC digit data
+ if (iDDL%2==0 && nTRM==3) continue;
+
+ // loop on TRM chain number per TRM
+ for (iChain=AliTOFGeometry::NChain()-1; iChain>=0; iChain--) {
+
+ // TRM chain trailer
+ MakeTRMchainTrailer(iChain, buf); fWordsPerTRM++;
+
+ // TRM TDC digits
+ MakeTDCdigits(nDDL, nTRM, iChain, buf);
+
+ // TRM chain header
+ MakeTRMchainHeader(nTRM, iChain, buf); fWordsPerTRM++;
+
+ } // end loop on iChain
+
+ // TRM global header
+ MakeTRMheader(nTRM, buf); fWordsPerTRM++;
+
+ // TRM filler in case where TRM data number is odd
+ if ((fWordsPerTRM+1)%2!=0) {
+ MakeTRMfiller(buf); fWordsPerTRM++;
+ }
+
+ MakeTRMtrailer(buf); fWordsPerDRM++;
+
+ fWordsPerDRM += fWordsPerTRM;
+
+ } // end loop on nTRM
+
+
+ // LTM section
+ //fIndex++;
+ //buf[fIndex] = MakeFiller(); fWordsPerDRM++; // valid till when LTM word number was 33
+ MakeLTMtrailer(buf); fWordsPerDRM++;
+ MakeLTMdata(buf); fWordsPerDRM+=48;
+ MakeLTMheader(buf); fWordsPerDRM++;
+
+ // DRM section: in
+ MakeDRMheader(nDDL, buf);
+
+ ReverseArray(buf, fIndex+1);
+
+ outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
+
+ for (jj=0; jj<(fIndex+1); jj++) buf[jj]=0;
+ fIndex = -1;
+
+ //Write REAL DATA HEADER
+ UInt_t currentFilePosition = outfile->Tellp();
+ sizeRawData = currentFilePosition - dataHeaderPosition - sizeof(header);
+ header.fSize = currentFilePosition - dataHeaderPosition;
+ header.SetAttribute(0); // valid data
+ outfile->Seekp(dataHeaderPosition);
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
+ outfile->Seekp(currentFilePosition);
+
+ delete outfile;
+
+ } //end loop on DDL file number
+
+ return 0;
-AliTOFDDLRawData::AliTOFDDLRawData(const AliTOFDDLRawData &source) :
- TObject(source){
- //Copy Constructor
- this->fIndex=source.fIndex;
- this->fVerbose=source.fVerbose;
- this->fTOFgeometry=source.fTOFgeometry;
- return;
}
-//---------------------------------------------------------------------------------------
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::GetDigits()
+{
+ //
+ // Fill the TOF volumes' map with the TOF digit indices
+ //
-AliTOFDDLRawData& AliTOFDDLRawData::operator=(const AliTOFDDLRawData &source){
- //Assigment operator
- this->fIndex=source.fIndex;
- this->fVerbose=source.fVerbose;
- this->fTOFgeometry=source.fTOFgeometry;
- return *this;
+ Int_t vol[5] = {-1,-1,-1,-1,-1};
+
+ Int_t digit = -1;
+ Int_t ndigits = fTOFdigitArray->GetEntries();
+ AliDebug(2, Form(" Number of read digits = %d",ndigits));
+ AliTOFdigit *digs;
+
+ // loop on TOF digits
+ for (digit=0; digit<ndigits; digit++) {
+ digs = (AliTOFdigit*)fTOFdigitArray->UncheckedAt(digit);
+
+ vol[0] = digs->GetSector(); // Sector Number (0-17)
+ vol[1] = digs->GetPlate(); // Plate Number (0-4)
+ vol[2] = digs->GetStrip(); // Strip Number (0-14/18)
+ vol[3] = digs->GetPadx(); // Pad Number in x direction (0-47)
+ vol[4] = digs->GetPadz(); // Pad Number in z direction (0-1)
+
+ fTOFdigitMap->AddDigit(vol, digit);
+
+ } // close loop on digit del TOF
+ AliDebug(2,Form(" Number of mapped digits = %d",fTOFdigitMap->GetFilledCellNumber()));
}
-//---------------------------------------------------------------------------------------
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeDRMheader(Int_t nDDL, UInt_t *buf)
+{
+ //
+ // DRM global header
+ //
+
+ //Int_t iDDL = fTOFrawStream->GetDDLnumberPerSector(nDDL);
+ Int_t iDDL = AliTOFRawStream::GetDDLnumberPerSector(nDDL);
+
+ //Int_t iSector = fTOFrawStream->GetSectorNumber(nDDL);
+ Int_t iSector = AliTOFRawStream::GetSectorNumber(nDDL);
+
+ UInt_t baseWord=0;
+ UInt_t word;
+
+ // DRM event CRC
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // event CRC --> CHANGED
+ AliBitPacking::PackWord(word,baseWord, 4,19);
+ word = 0; // reserved for future use
+ AliBitPacking::PackWord(word,baseWord,20,27);
+ word = 4; // 0100 -> DRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // DRM status header 4
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // temperature
+ AliBitPacking::PackWord(word,baseWord, 4,13);
+ word = 0; // zero
+ AliBitPacking::PackWord(word,baseWord, 14,14);
+ word = 0; // ACK
+ AliBitPacking::PackWord(word,baseWord, 15,15);
+ word = 0; // Sens AD
+ AliBitPacking::PackWord(word,baseWord, 16,18);
+ word = 0; // zero
+ AliBitPacking::PackWord(word,baseWord, 19,19);
+ word = 0; // reserved for future use
+ AliBitPacking::PackWord(word,baseWord, 20,27);
+ word = 4; // 0100 -> DRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // DRM status header 3
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // L0 BCID
+ AliBitPacking::PackWord(word,baseWord, 4,15);
+ word = 0; // Run Time info
+ AliBitPacking::PackWord(word,baseWord, 16,27);
+ word = 4; // 0100 -> DRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // DRM status header 2
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+
+ if (iDDL%2==1) {
+ word = 2047; // enable ID: [00000000000;11111111111] for odd
+ // (i.e. right) crates
+ AliBitPacking::PackWord(word,baseWord, 4,14);
+ } else {
+ word = 2045; // enable ID: [00000000000;11111111101] for even
+ // (i.e. left) crates
+ AliBitPacking::PackWord(word,baseWord, 4,14);
+ }
+
+ word = 0; //
+ AliBitPacking::PackWord(word,baseWord,15,15);
+ word = 0; // fault ID
+ AliBitPacking::PackWord(word,baseWord,16,26);
+ word = 0; // RTO
+ AliBitPacking::PackWord(word,baseWord,27,27);
+ word = 4; // 0100 -> DRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // DRM status header 1
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
-void AliTOFDDLRawData::GetDigits(TClonesArray *TOFdigits,Int_t nDDL,UInt_t *buf)
+ if (iDDL%2==1) {
+ word = 2047; // slot ID: [00000000000;11111111111] for odd
+ // (i.e. right) crates
+ AliBitPacking::PackWord(word,baseWord, 4,14);
+ } else {
+ word = 2045; // slot ID: [00000000000;11111111101] for even
+ // (i.e. left) crates
+ AliBitPacking::PackWord(word,baseWord, 4,14);
+ }
+
+ word = 1; // LHC clock status: 1/0
+ AliBitPacking::PackWord(word,baseWord,15,15);
+ word = 0; // Vers ID
+ AliBitPacking::PackWord(word,baseWord,16,20);
+ word = 0; // DRMH size
+ AliBitPacking::PackWord(word,baseWord,21,24);
+ word = 0; // reserved for future use
+ AliBitPacking::PackWord(word,baseWord,25,27);
+ word = 4; // 0100 -> DRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // DRM global header
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = fIndex+1 + 1; // event words
+ AliBitPacking::PackWord(word,baseWord, 4,20);
+ word = iDDL; // crate ID [0;3]
+ AliBitPacking::PackWord(word,baseWord,21,22);
+ word = iSector; // sector ID [0;17]
+ AliBitPacking::PackWord(word,baseWord,23,27);
+ word = 4; // 0100 -> DRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeDRMtrailer(UInt_t *buf)
{
+ //
+ // DRM global trailer
+ //
+
+ UInt_t baseWord;
+ UInt_t word;
+
+ baseWord=0;
+ word = 1; // 0001 -> DRM data are coming from the VME slot number 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // local event counter --> TO BE CHANGED IN fWordsPerDRM+5
+ AliBitPacking::PackWord(word,baseWord, 4,15);
+ word = 0; // reserved for future use
+ AliBitPacking::PackWord(word,baseWord,16,27);
+ word = 5; // 0101 -> DRM trailer ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
- //This method packs the TOF digits in a proper 32 bits structure
-
- Int_t iDDL = nDDL%AliTOFGeometry::NDDL();
- Int_t iSector = (Int_t)((Float_t)nDDL/AliTOFGeometry::NDDL());
- Int_t iTRM = 0;
- Int_t iTDC = 0;
- Int_t iCH =-1;
- Int_t sector;
- Int_t plate;
- Int_t strip;
- Int_t padx;
- Int_t padz;
- Int_t totCharge;
- Int_t timeOfFlight;
- Int_t error=0;
- Int_t eureka;
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeLTMheader(UInt_t *buf)
+{
+ //
+ // LTM header
+ //
+
+ UInt_t baseWord;
UInt_t word;
+
+ baseWord=0;
+ word = 2; // 0010 -> LTM data are coming from the VME slot number 2
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 35; // event words
+ AliBitPacking::PackWord(word,baseWord, 4,16);
+ word = 0; // crc error
+ AliBitPacking::PackWord(word,baseWord,17,17);
+ word = 0; // fault
+ AliBitPacking::PackWord(word,baseWord,18,23);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 4; // 0100 -> LTM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeLTMdata(UInt_t *buf)
+{
+ //
+ // LTM data
+ //
+
UInt_t baseWord;
- Int_t ndigits = TOFdigits->GetEntries();
- AliTOFdigit *digs;
- ofstream ftxt;
- if(!ndigits)
- {
- AliError("No found TOF digits");
- return;
- }
+ UInt_t word;
- if (fVerbose==2) ftxt.open("TOFdigits.txt",ios::app);
+ baseWord=0;
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,0,9);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,10,19);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,20,29);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,30,30);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+
+ // OR45, OR46, OR47
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR42, OR43, OR44
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR39, OR40, OR41
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR36, OR37, OR38
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR33, OR34, OR35
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR30, OR31, OR32
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR27, OR28, OR29
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR24, OR25, OR26
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR21, OR22, OR23
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR18, OR19, OR20
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR15, OR16, OR17
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR12, OR12, OR24
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR9, OR10, OR11
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR6, OR7, OR8
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR3, OR4, OR5
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // OR0, OR1, OR2
+ fIndex++;
+ buf[fIndex]=baseWord;
- for (Int_t digit=0;digit<ndigits;digit++) {
- digs = (AliTOFdigit*)TOFdigits->UncheckedAt(digit);
- sector = digs->GetSector(); // Sector Number (0-17)
- plate = digs->GetPlate(); // Plate Number (0-4)
- strip = digs->GetStrip(); // Strip Number (0-14/18/19)
- padx = digs->GetPadx(); // Pad Number in x direction (0-47)
- padz = digs->GetPadz(); // Pad Number in z direction (0-1)
- eureka = digs->GetTotPad(fTOFgeometry); // Global Pad Number inside a Sector
- // totCharge = (Int_t)digs->GetAdc(); //Use realistic ToT, for Standard production with no miscalibration/Slewing it == fAdC in digit (see AliTOFDigitizer)
- totCharge = (Int_t)digs->GetToT();
- timeOfFlight = (Int_t)digs->GetTdc();
-
- if (sector!=iSector || (Int_t)((Float_t)eureka/AliTOFGeometry::NPadXTRM()/AliTOFGeometry::NTRM())!=iDDL) continue;
-
- if (fVerbose==2) ftxt << " Sector: " << sector << " Plate: " << plate << " Strip: " << strip << " PadZ: " << padz << " PadX: " << padx << " totPadNumber " << eureka << endl;
-
- iTRM = (Int_t)((Float_t)eureka/AliTOFGeometry::NPadXTRM()) - iDDL*AliTOFGeometry::NTRM();
-
- iTDC = (Int_t)((Float_t)eureka/AliTOFGeometry::NCh()) - (iDDL*AliTOFGeometry::NTRM() + iTRM) * AliTOFGeometry::NTdc();
- /*
- iTDC = (Int_t)(AliTOFGeometry::NTdc()*
- (
- (Float_t)eureka/AliTOFGeometry::NPadXTRM() -
- (Int_t)((Float_t)eureka/AliTOFGeometry::NPadXTRM())
- )
- );
- */
-
- iCH = eureka - ((iDDL*AliTOFGeometry::NTRM() + iTRM) * AliTOFGeometry::NTdc() + iTDC) * AliTOFGeometry::NCh();
- /*
- iCH = (Int_t)(AliTOFGeometry::NCh() *
- (
- (Float_t)eureka/AliTOFGeometry::NPadXTRM()*AliTOFGeometry::NTdc() - (Int_t)((Float_t)eureka/AliTOFGeometry::NPadXTRM()*AliTOFGeometry::NTdc()) -
- (Int_t)((Float_t)eureka/AliTOFGeometry::NPadXTRM()*AliTOFGeometry::NTdc() - (Int_t)((Float_t)eureka/AliTOFGeometry::NPadXTRM()*AliTOFGeometry::NTdc()))
- )
- );
- */
-
- if (fVerbose==2) ftxt << "DDL: "<<nDDL<<" TRM: "<<iTRM<<" TDC: "<<iTDC<<" Channel: "<<iCH<<" totCharge: "<<totCharge<<" tof: "<<timeOfFlight<<endl;
-
- AliDebug(2,Form("%2i %2i %2i %2i %2i %2i %2i %2i %2i %7i %8i",nDDL,iTRM,iTDC,iCH,sector,plate,strip,padz,padx,totCharge,timeOfFlight));
-
- baseWord=0;
-
- word=iTRM;
- AliBitPacking::PackWord(word,baseWord, 0, 3); // TRM ID
- word=iTDC;
- AliBitPacking::PackWord(word,baseWord, 4, 8); // TDC ID
- word=iCH;
- AliBitPacking::PackWord(word,baseWord, 9,11); // CH ID
-
- // temporary control
- if (totCharge<0) word=TMath::Abs(totCharge);
- else word=totCharge;
- AliBitPacking::PackWord(word,baseWord,12,31); // Charge (TOT) // v0.01
- //AliBitPacking::PackWord(word,baseWord,12,19); // Charge (TOT) // v0.02
- //AliBitPacking::PackWord(0,baseWord,20,31); // v0.02
-
- fIndex++;
- buf[fIndex]=baseWord;
-
- baseWord=0;
-
- word=error;
- AliBitPacking::PackWord(word,baseWord,0, 7); // Error flag
- word=timeOfFlight;
- AliBitPacking::PackWord(word,baseWord,8,31); // time-of-flight // v0.01
- //AliBitPacking::PackWord(word,baseWord,8,19); // time-of-flight // v0.02
- //AliBitPacking::PackWord(0,baseWord,20,30); // v0.02
- //AliBitPacking::PackWord(1,baseWord,31,31); // v0.02
-
- fIndex++;
- buf[fIndex]=baseWord;
- word=0;
- baseWord=0;
-
- }//end for
- if (fVerbose==2) ftxt.close();
+ baseWord=0;
+ word = 100; // Local temperature in LTM11 -> 4 X 25 degree (environment temperature)
+ AliBitPacking::PackWord(word,baseWord, 0, 9);
+ word = 100; // Local temperature in LTM10 -> 4 X 25 degree (environment temperature)
+ AliBitPacking::PackWord(word,baseWord,10,19);
+ word = 100; // Local temperature in LTM9 -> 4 X 25 degree (environment temperature)
+ AliBitPacking::PackWord(word,baseWord,20,29);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,30,30);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,31,31);
- return;
+ fIndex++;
+ buf[fIndex]=baseWord;
-}//end GetDigits
+ // Local temperature in LTM8, LMT7, LTM6 -> 4 X 25 degree (environment temperature)
+ fIndex++;
+ buf[fIndex]=baseWord;
-//---------------------------------------------------------------------------------------
+ // Local temperature in LTM5, LMT4, LTM3 -> 4 X 25 degree (environment temperature)
+ fIndex++;
+ buf[fIndex]=baseWord;
-Int_t AliTOFDDLRawData::RawDataTOF(TBranch* branch){
+ // Local temperature in LTM2, LMT1, LTM0 -> 4 X 25 degree (environment temperature)
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+
+
+ // Local temperature in T7, T6, T5 -> 4 X 25 degree (environment temperature)
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Local temperature in T4, T3, T2 -> 4 X 25 degree (environment temperature)
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Local temperature in T1, T0 -> 4 X 25 degree (environment temperature)
+ // Local temperature in VTH15 -> Thereshould voltage for FEAC15
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ baseWord=0;
+ word = 0; // VTH13 -> Thereshould voltage for FEAC16
+ AliBitPacking::PackWord(word,baseWord, 0, 9);
+ word = 0; // VTH14 -> Thereshould voltage for FEAC14
+ AliBitPacking::PackWord(word,baseWord,10,19);
+ word = 0; // GND-FEAC7 -> Voltage drop between GND and FEAC7
+ AliBitPacking::PackWord(word,baseWord,20,29);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,30,30);
+ word = 0;
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // VTH11 -> Thereshould voltage for FEAC11
+ // VTH12 -> Thereshould voltage for FEAC12
+ // GND-FEAC6 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // VTH9 -> Thereshould voltage for FEAC11
+ // VTH10 -> Thereshould voltage for FEAC12
+ // GND-FEAC5 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // VTH7 -> Thereshould voltage for FEAC11
+ // VTH8 -> Thereshould voltage for FEAC12
+ // GND-FEAC4 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // VTH5 -> Thereshould voltage for FEAC11
+ // VTH6 -> Thereshould voltage for FEAC12
+ // GND-FEAC3 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // VTH3 -> Thereshould voltage for FEAC11
+ // VTH4 -> Thereshould voltage for FEAC12
+ // GND-FEAC2 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // VTH1 -> Thereshould voltage for FEAC11
+ // VTH2 -> Thereshould voltage for FEAC12
+ // GND-FEAC1 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV15
+ // VTH0 -> Thereshould voltage for FEAC12
+ // GND-FEAC0 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV15
+ // VTH0 -> Thereshould voltage for FEAC12
+ // GND-FEAC0 -> Voltage drop between GND and FEAC6
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV12
+ // LV13
+ // LV14
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV9
+ // LV10
+ // LV11
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV6
+ // LV7
+ // LV8
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV3
+ // LV4
+ // LV5
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // LV0
+ // LV1
+ // LV2
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+
+
+ baseWord=0;
+ word = 0; // PDL45 -> Delay Line setting for PDL41
+ AliBitPacking::PackWord(word,baseWord, 0, 7);
+ word = 0; // PDL46 -> Delay Line setting for PDL42
+ AliBitPacking::PackWord(word,baseWord, 8,15);
+ word = 0; // PDL47 -> Delay Line setting for PDL43
+ AliBitPacking::PackWord(word,baseWord,16,23);
+ word = 0; // PDL48 -> Delay Line setting for PDL44
+ AliBitPacking::PackWord(word,baseWord,24,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL37, PDL38, PDL39, PDL40
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL33, PDL34, PDL35, PDL36
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL29, PDL30, PDL31, PDL32
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL25, PDL26, PDL27, PDL28
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL21, PDL22, PDL23, PDL24
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL17, PDL18, PDL19, PDL20
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL13, PDL14, PDL15, PDL16
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL9, PDL10, PDL11, PDL12
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL5, PDL6, PDL7, PDL8
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+ // Delay Line setting for PDL1, PDL2, PDL3, PDL4
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeLTMtrailer(UInt_t *buf)
+{
//
- // This method creates the Raw data files for TOF detector
+ // LTM trailer
//
+
+ UInt_t baseWord;
+ UInt_t word;
+
+ baseWord=0;
+ word = 2; // 0010 -> LTM data are coming from the VME slot number 2
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // event crc
+ AliBitPacking::PackWord(word,baseWord, 4,15);
+ word = 0; // event number
+ AliBitPacking::PackWord(word,baseWord,16,27);
+ word = 5; // 0101 -> LTM trailer ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
- const Int_t kSize = 5000; //max number of digits per DDL file times 2
+}
- UInt_t buf[kSize];
- //UInt_t baseWord; // v0.01
- //UInt_t word; // v0.01
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeTRMheader(Int_t nTRM, UInt_t *buf)
+{
+ //
+ // TRM header for the TRM number nTRM [ 3;12]
+ //
- fIndex=-1;
+ if (nTRM<3 || nTRM>12) {
+ AliWarning(Form(" TRM number is out of the right range [3;12] (nTRM = %3i",nTRM));
+ return;
+ }
- TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
- char fileName[15];
- ofstream outfile; // logical name of the output file
- AliRawDataHeader header;
- UInt_t sizeRawData = 0;
+ UInt_t baseWord;
+ UInt_t word;
- //loop over TOF DDL files
- for(Int_t i = 0; i<AliTOFGeometry::NDDL()*AliTOFGeometry::NSectors(); i++){
-
- sprintf(fileName,"TOF_%d.ddl",i+AliTOFRawStream::kDDLOffset); //The name of the output file
-#ifndef __DECCXX
- outfile.open(fileName,ios::binary);
-#else
- outfile.open(fileName);
-#endif
-
- //write Dummy DATA HEADER
- UInt_t dataHeaderPosition=outfile.tellp();
- outfile.write((char*)(&header),sizeof(header));
-
- /*
- // v0.01
- baseWord=0;
- word=i;
- //AliBitPacking::PackWord(word,baseWord,0, 31); // Number of DDL file
- AliBitPacking::PackWord(word,baseWord,0, 6); // Number of DDL file
- AliBitPacking::PackWord(0,baseWord,7,31);
-
- fIndex++;
- buf[fIndex]=baseWord;
- */
-
- branch->GetEvent();
-
- //For each DDL file, buf contains the array of data words in Binary format
- //fIndex gives the number of 32 bits words in the buffer for each DDL file
- GetDigits(digits,i,buf);
- outfile.write((char *)buf,((fIndex+1)*sizeof(UInt_t)));
-
- for(Int_t ii=0;ii<(fIndex+1);ii++) buf[ii]=0;
- fIndex=-1;
-
- //Write REAL DATA HEADER
- UInt_t currentFilePosition=outfile.tellp();
- sizeRawData = currentFilePosition - dataHeaderPosition - sizeof(header);
- header.fSize=currentFilePosition-dataHeaderPosition;
- header.SetAttribute(0); // valid data
- outfile.seekp(dataHeaderPosition);
- outfile.write((char*)(&header),sizeof(header));
- outfile.seekp(currentFilePosition);
+ baseWord = 0;
+ word = nTRM; // TRM data coming from the VME slot number nTRM
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // event words
+ AliBitPacking::PackWord(word,baseWord, 4,16);
+
+ if (fPackedAcquisition)
+ word = 0; // ACQuisition mode: [0;3] see document
+ else
+ word = 3; // ACQuisition mode: [0;3] see document
+ AliBitPacking::PackWord(word,baseWord,17,18);
+ word = 0; // description of a SEU inside LUT tables for INL compensation;
+ // the data are unaffected
+ AliBitPacking::PackWord(word,baseWord,19,19);
+ word = 0; // Must Be Zero (MBZ)
+ AliBitPacking::PackWord(word,baseWord,20,27);
+ word = 4; // 0100 -> TRM header ID
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeTRMtrailer(UInt_t *buf)
+{
+ //
+ // Set TRM Global Trailer
+ // with the calculated CRC
+ //
+
+ UInt_t baseWord;
+ UInt_t word;
+
+ baseWord=0;
+ word = 15; // 1111 -> TRM trailer ID 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+
+ UInt_t trmCRC=0x0;
+ for (Int_t ii=fIndex-(fWordsPerTRM-1); ii<fIndex; ii++)
+ trmCRC ^= buf[ii];
+ //PH printf(" A trmCRC=%d\n",trmCRC);
+
+ word = 0x0;
+ word ^= ( (trmCRC & 0x00000fff) >> 0);
+ word ^= ( (trmCRC & 0x00fff000) >> 12);
+ word ^= ( (trmCRC & 0xff000000) >> 24);
+
+ //PH printf(" B trmCRC=%d\n",word);
+
+ AliBitPacking::PackWord(word,baseWord, 4,15); // event CRC --> CHANGED
+
+ word = 0; // local event counter == DRM local event counter --> TO BE CHANGED
+ AliBitPacking::PackWord(word,baseWord,16,27);
+ word = 5; // 0101 -> TRM trailer ID 2
+ AliBitPacking::PackWord(word,baseWord,28,31);
+
+ fIndex++;
+ for (Int_t ii=fIndex; ii>fIndex-fWordsPerTRM; ii--)
+ buf[ii]=buf[ii-1];
+
+ buf[fIndex-fWordsPerTRM] = baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeTRMchainHeader(Int_t nTRM, Int_t iChain,
+ UInt_t *buf)
+{
+ //
+ // TRM chain header
+ //
+
+ UInt_t baseWord;
+ UInt_t word;
+
+ if (nTRM<3 || nTRM>12) {
+ AliWarning(Form(" TRM number is out of the right range [3;12] (nTRM = %3i", nTRM));
+ return;
+ }
+
+ if (iChain<0 || iChain>1) {
+ AliWarning(Form(" Chain number is out of the right range [0;1] (iChain = %3i", iChain));
+ return;
+ }
+
+ baseWord=0;
+ word = nTRM; // TRM data coming from the VME slot ID nTRM
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // bunch ID
+ AliBitPacking::PackWord(word,baseWord, 4,15);
+ word = 0;//100; // PB24 temperature -> 4 X 25 degree (environment temperature)
+ AliBitPacking::PackWord(word,baseWord,16,23);
+ word = 0;//(Int_t)(5 * gRandom->Rndm()); // PB24 ID [0;4]
+ AliBitPacking::PackWord(word,baseWord,24,26);
+ word = 0; // TS
+ AliBitPacking::PackWord(word,baseWord,27,27);
+ switch (iChain) {
+ case 0:
+ word = 0; // 0000 -> TRM chain 0 ID
+ break;
+ case 1:
+ word = 2; // 0010 -> TRM chain 1 ID
+ break;
+ }
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeTRMfiller(UInt_t *buf)
+{
+ //
+ // TRM filler
+ //
- outfile.close();
+ Int_t jj = -1;
- }//end for
+ fIndex++;
+ for (jj=fIndex; jj>fIndex-fWordsPerTRM; jj--)
+ buf[jj] = buf[jj-1];
+
+ buf[fIndex-fWordsPerTRM] = MakeFiller();
+
+}
- return 0;
+//----------------------------------------------------------------------------
+UInt_t AliTOFDDLRawData::MakeFiller() const
+{
+ //
+ // Filler word definition: to make even the number of words per TRM/LTM
+ //
+
+ UInt_t baseWord;
+ UInt_t word;
+
+ baseWord=0;
+ word = 0; // 0000 -> filler ID 1
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // MBZ
+ AliBitPacking::PackWord(word,baseWord, 4,27);
+ word = 7; // 0111 -> filler ID 2
+ AliBitPacking::PackWord(word,baseWord, 28,31);
+
+ return baseWord;
+
}
-//---------------------------------------------------------------------------
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeTRMchainTrailer(Int_t iChain, UInt_t *buf)
+{
+ //
+ // TRM chain trailer
+ //
+
+ if (iChain<0 || iChain>1) {
+ AliWarning(Form(" Chain number is out of the right range [0;1] (iChain = %3i", iChain));
+ return;
+ }
+
+ UInt_t baseWord;
+ UInt_t word;
+
+ baseWord=0;
+ word = 0; // status
+ AliBitPacking::PackWord(word,baseWord, 0, 3);
+ word = 0; // MBZ
+ AliBitPacking::PackWord(word,baseWord, 4,15);
+ word = 0; // event counter --> TO BE CHANGED
+ AliBitPacking::PackWord(word,baseWord,16,27);
+ switch (iChain) {
+ case 0:
+ word = 1; // 0001 -> TRM chain 0 trailer ID
+ break;
+ case 1:
+ word = 3; // 0011 -> TRM chain 1 trailer ID
+ break;
+ }
+ AliBitPacking::PackWord(word,baseWord,28,31);
+ fIndex++;
+ buf[fIndex]=baseWord;
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::MakeTDCdigits(Int_t nDDL, Int_t nTRM, Int_t iChain, UInt_t *buf)
+{
+ //
+ // TRM TDC digit
+ //
+
+ const Double_t kOneMoreFilledCell = 1./(AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors());
+ Double_t percentFilledCells = Double_t(fTOFdigitMap->GetFilledCellNumber())/(AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors());
+
+ if (nDDL<0 || nDDL>71) {
+ AliWarning(Form(" DDL number is out of the right range [0;71] (nDDL = %3i", nDDL));
+ return;
+ }
+
+ if (nTRM<3 || nTRM>12) {
+ AliWarning(Form(" TRM number is out of the right range [3;12] (nTRM = %3i", nTRM));
+ return;
+ }
+
+ if (iChain<0 || iChain>1) {
+ AliWarning(Form(" Chain number is out of the right range [0;1] (iChain = %3i", iChain));
+ return;
+ }
+
+ Int_t psArray[1000];
+ UInt_t localBuffer[1000];
+ Int_t localIndex = -1;
+
+ Int_t iDDL = nDDL%AliTOFGeometry::NDDL();
+
+ Int_t volume[5] = {-1, -1, -1, -1, -1};
+ Int_t indexDigit[3] = {-1, -1, -1};
+
+ Int_t totCharge = -1;
+ Int_t timeOfFlight = -1;
+
+ Int_t trailingSpurious = -1;
+ Int_t leadingSpurious = -1;
+
+ AliTOFdigit *digs;
+
+ UInt_t baseWord=0;
+ UInt_t word=0;
+
+ Int_t jj = -1;
+ Int_t nTDC = -1;
+ Int_t iCH = -1;
+
+ //Int_t numberOfMeasuresPerChannel = 0;
+ //Int_t maxMeasuresPerChannelInTDC = 0;
+
+ Bool_t outOut = HeadOrTail();
+
+ ofstream ftxt;
+
+ if (fVerbose==2) ftxt.open("TOFdigits.txt",ios::app);
+
+ for (jj=0; jj<5; jj++) volume[jj] = -1;
+
+ // loop on TDC number
+ for (nTDC=AliTOFGeometry::NTdc()-1; nTDC>=0; nTDC--) {
+
+ // the DRM odd (i.e. left) slot number 3 doesn't contain TDC digit data
+ // for TDC numbers 3-14
+ if (iDDL%2==1 && nTRM==3 && (Int_t)(nTDC/3)!=0) continue;
+
+ // loop on TDC channel number
+ for (iCH=AliTOFGeometry::NCh()-1; iCH>=0; iCH--) {
+
+ //numberOfMeasuresPerChannel = 0;
+
+ for (Int_t aa=0; aa<5; aa++) volume[aa]=-1;
+ AliTOFRawStream::EquipmentId2VolumeId(nDDL, nTRM, iChain, nTDC, iCH, volume);
+
+ AliDebug(3,Form(" volume -> %2d %1d %2d %2d %1d",volume[0],volume[1],volume[2],volume[3],volume[4]));
+
+ if (volume[0]==-1 || volume[1]==-1 || volume[2]==-1 ||
+ volume[3]==-1 || volume[4]==-1) continue;
+
+ AliDebug(3,Form(" ====== %2d %1d %2d %2d %1d",volume[0],volume[1],volume[2],volume[3],volume[4]));
+
+ for (jj=0; jj<3; jj++) indexDigit[jj] = -1;
+
+ fTOFdigitMap->GetDigitIndex(volume, indexDigit);
+
+ if (indexDigit[0]<0) {
+
+ trailingSpurious = Int_t(2097152*gRandom->Rndm());
+ leadingSpurious = Int_t(2097152*gRandom->Rndm());
+
+ if ( fFakeOrphaneProduction &&
+ ( ( fPackedAcquisition && percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) ) ||
+ (!fPackedAcquisition && percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) ) ) ) {
+
+ percentFilledCells+=kOneMoreFilledCell;
+
+ Int_t dummyPS = 0;
+
+ if (outOut) {
+ word = trailingSpurious; // trailing edge measurement
+ dummyPS = 2;
+ }
+ else {
+ word = leadingSpurious; // leading edge measurement
+ dummyPS = 1;
+ }
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ ftxt << " " << -1;
+ if (word<10) ftxt << " " << word;
+ else if (word>=10 && word<100) ftxt << " " << word;
+ else if (word>=100 && word<1000) ftxt << " " << word;
+ else if (word>=1000 && word<10000) ftxt << " " << word;
+ else if (word>=10000 && word<100000) ftxt << " " << word;
+ else if (word>=100000 && word<1000000) ftxt << " " << word;
+ else ftxt << " " << word;
+ ftxt << " " << dummyPS << endl;
+ }
+
+ AliBitPacking::PackWord(word,baseWord, 0,20);
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = dummyPS; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+ psArray[localIndex]=dummyPS;
+
+ baseWord=0;
+
+ } // if ( fFakeOrphaneProduction && ( ( fPackedAcquisition && percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) ) or ... ) )
+ } // if (indexDigit[0]<0)
+
+ for (jj=0; jj<3;jj++) {
+
+ if (indexDigit[jj]<0) continue;
+
+ AliDebug(3,Form(" ====== %2d %1d %2d %2d %1d -> %1d %d",volume[0],volume[1],volume[2],volume[3],volume[4],jj,indexDigit[jj]));
+
+ digs = (AliTOFdigit*)fTOFdigitArray->UncheckedAt(indexDigit[jj]);
+
+ if (digs->GetSector()!=volume[0] ||
+ digs->GetPlate() !=volume[1] ||
+ digs->GetStrip() !=volume[2] ||
+ digs->GetPadx() !=volume[3] ||
+ digs->GetPadz() !=volume[4]) AliWarning(Form(" --- ERROR --- %2i (%2i) %1i (%1i) %2i (%2i) %2i (%2i) %1i (%1i)",
+ digs->GetSector(), volume[0],
+ digs->GetPlate(), volume[1],
+ digs->GetStrip(), volume[2],
+ digs->GetPadx(), volume[3],
+ digs->GetPadz(), volume[4])
+ );
+
+ timeOfFlight = (Int_t)(digs->GetTdc());
+ /*+
+ fTOFCableLengthMap->GetCableTimeShiftBin(nDDL, nTRM, iChain, nTDC)*/;
+
+ //if (timeOfFlight>=fMatchingWindow+1000) continue; // adc
+ //if (timeOfFlight>=fMatchingWindow) continue; // adc
+ if (digs->GetTdcND()>=fMatchingWindow) {
+ AliDebug(2,"Out of matching window.");
+ continue; // adc
+ }
+ else AliDebug(2,"In matching window");
+
+ //numberOfMeasuresPerChannel++;
+
+ // totCharge = (Int_t)digs->GetAdc(); //Use realistic ToT, for Standard production with no miscalibration/Slewing it == fAdC in digit (see AliTOFDigitizer)
+ totCharge = (Int_t)(digs->GetToT());
+ // temporary control
+ if (totCharge<0) totCharge = 0;//TMath::Abs(totCharge);
+
+ if (fPackedAcquisition) {
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ if (totCharge<10) ftxt << " " << totCharge;
+ else if (totCharge>=10 && totCharge<100) ftxt << " " << totCharge;
+ else ftxt << " " << totCharge;
+ if (timeOfFlight<10) ftxt << " " << timeOfFlight << endl;
+ else if (timeOfFlight>=10 && timeOfFlight<100) ftxt << " " << timeOfFlight << endl;
+ else if (timeOfFlight>=100 && timeOfFlight<1000) ftxt << " " << timeOfFlight << endl;
+ else ftxt << " " << timeOfFlight << endl;
+ }
+
+ word = timeOfFlight%8192; // time-of-fligth measurement
+ AliBitPacking::PackWord(word,baseWord, 0,12);
+
+ if (totCharge>=256) totCharge = 255;
+ word = totCharge; // time-over-threshould measurement
+ AliBitPacking::PackWord(word,baseWord,13,20);
+
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = 0; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+
+ baseWord=0;
+
+ if ( fFakeOrphaneProduction &&
+ percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) ) {
+
+ percentFilledCells+=kOneMoreFilledCell;
+
+ trailingSpurious = Int_t(2097152*gRandom->Rndm());
+ leadingSpurious = Int_t(2097152*gRandom->Rndm());
+
+ Int_t dummyPS = 0;
+
+ if (outOut) {
+ word = trailingSpurious; // trailing edge measurement
+ dummyPS = 2;
+ }
+ else {
+ word = leadingSpurious; // leading edge measurement
+ dummyPS = 1;
+ }
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ ftxt << " " << -1;
+ if (word<10) ftxt << " " << word;
+ else if (word>=10 && word<100) ftxt << " " << word;
+ else if (word>=100 && word<1000) ftxt << " " << word;
+ else if (word>=1000 && word<10000) ftxt << " " << word;
+ else if (word>=10000 && word<100000) ftxt << " " << word;
+ else if (word>=100000 && word<1000000) ftxt << " " << word;
+ else ftxt << " " << word;
+ ftxt << " " << dummyPS << endl;
+ }
+
+ AliBitPacking::PackWord(word,baseWord, 0,20);
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = dummyPS; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+ psArray[localIndex]=dummyPS;
+
+ baseWord=0;
+
+ } // if ( fFakeOrphaneProduction && percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) )
+
+
+ } // if (fPackedAcquisition)
+ else { // if (!fPackedAcquisition)
+
+ if ( fFakeOrphaneProduction &&
+ percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && outOut ) {
+
+ percentFilledCells+=kOneMoreFilledCell;
+
+ trailingSpurious = Int_t(2097152*gRandom->Rndm());
+ word = trailingSpurious;
+ Int_t dummyPS = 2;
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ ftxt << " " << -1;
+ if (word<10) ftxt << " " << word;
+ else if (word>=10 && word<100) ftxt << " " << word;
+ else if (word>=100 && word<1000) ftxt << " " << word;
+ else if (word>=1000 && word<10000) ftxt << " " << word;
+ else if (word>=10000 && word<100000) ftxt << " " << word;
+ else if (word>=100000 && word<1000000) ftxt << " " << word;
+ else ftxt << " " << word;
+ ftxt << " " << dummyPS << endl;
+ }
+
+ AliBitPacking::PackWord(word,baseWord, 0,20);
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = dummyPS; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+ psArray[localIndex]=dummyPS;
+
+ baseWord=0;
+
+ } // if ( fFakeOrphaneProduction && percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && outOut )
+
+
+ word = (timeOfFlight + Int_t(totCharge*AliTOFGeometry::ToTBinWidth()/AliTOFGeometry::TdcBinWidth()))%2097152; // trailing edge measurement
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ ftxt << " " << -1;
+ if (word<10) ftxt << " " << word;
+ else if (word>=10 && word<100) ftxt << " " << word;
+ else if (word>=100 && word<1000) ftxt << " " << word;
+ else if (word>=1000 && word<10000) ftxt << " " << word;
+ else if (word>=10000 && word<100000) ftxt << " " << word;
+ else if (word>=100000 && word<1000000) ftxt << " " << word;
+ else ftxt << " " << word;
+ ftxt << " " << 2 << endl;
+ }
+
+ AliBitPacking::PackWord(word,baseWord, 0,20);
+
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = 2; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+ psArray[localIndex]=2;
+
+ baseWord=0;
+
+ word = timeOfFlight%2097152; // leading edge measurement
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ ftxt << " " << -1;
+ if (word<10) ftxt << " " << word;
+ else if (word>=10 && word<100) ftxt << " " << word;
+ else if (word>=100 && word<1000) ftxt << " " << word;
+ else if (word>=1000 && word<10000) ftxt << " " << word;
+ else if (word>=10000 && word<100000) ftxt << " " << word;
+ else if (word>=100000 && word<1000000) ftxt << " " << word;
+ else ftxt << " " << word;
+ ftxt << " " << 1 << endl;
+ }
+
+ AliBitPacking::PackWord(word,baseWord, 0,20);
+
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = 1; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+ psArray[localIndex]=1;
+
+ baseWord=0;
+
+
+ if ( fFakeOrphaneProduction &&
+ percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && !outOut ) {
+
+ percentFilledCells+=kOneMoreFilledCell;
+
+ leadingSpurious = Int_t(2097152*gRandom->Rndm());
+ word = leadingSpurious;
+ Int_t dummyPS = 1;
+
+ if (fVerbose==2) {
+ if (nDDL<10) ftxt << " " << nDDL;
+ else ftxt << " " << nDDL;
+ if (nTRM<10) ftxt << " " << nTRM;
+ else ftxt << " " << nTRM;
+ ftxt << " " << iChain;
+ if (nTDC<10) ftxt << " " << nTDC;
+ else ftxt << " " << nTDC;
+ ftxt << " " << iCH;
+ if (volume[0]<10) ftxt << " -> " << volume[0];
+ else ftxt << " -> " << volume[0];
+ ftxt << " " << volume[1];
+ if (volume[2]<10) ftxt << " " << volume[2];
+ else ftxt << " " << volume[2];
+ ftxt << " " << volume[4];
+ if (volume[3]<10) ftxt << " " << volume[3];
+ else ftxt << " " << volume[3];
+ ftxt << " " << -1;
+ if (word<10) ftxt << " " << word;
+ else if (word>=10 && word<100) ftxt << " " << word;
+ else if (word>=100 && word<1000) ftxt << " " << word;
+ else if (word>=1000 && word<10000) ftxt << " " << word;
+ else if (word>=10000 && word<100000) ftxt << " " << word;
+ else if (word>=100000 && word<1000000) ftxt << " " << word;
+ else ftxt << " " << word;
+ ftxt << " " << dummyPS << endl;
+ }
+
+ AliBitPacking::PackWord(word,baseWord, 0,20);
+ word = iCH; // TDC channel ID [0;7]
+ AliBitPacking::PackWord(word,baseWord,21,23);
+ word = nTDC; // TDC ID [0;14]
+ AliBitPacking::PackWord(word,baseWord,24,27);
+ word = 0; // error flag
+ AliBitPacking::PackWord(word,baseWord,28,28);
+ word = dummyPS; // Packing Status [0;3]
+ AliBitPacking::PackWord(word,baseWord,29,30);
+ word = 1; // TRM TDC digit ID
+ AliBitPacking::PackWord(word,baseWord,31,31);
+
+ localIndex++; fWordsPerTRM++;
+ localBuffer[localIndex]=baseWord;
+ psArray[localIndex]=dummyPS;
+
+ baseWord=0;
+
+ } // if ( fFakeOrphaneProduction && percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && !outOut )
+
+
+ } // if (!fPackedAcquisition)
+
+ } //end loop on digits in the same volume
+
+ //if (numberOfMeasuresPerChannel>maxMeasuresPerChannelInTDC)
+ //maxMeasuresPerChannelInTDC = numberOfMeasuresPerChannel;
+
+ } // end loop on TDC channel number
+
+ //AliDebug(3,Form(" TDC number %2i: numberOfMeasuresPerChannel = %2i --- maxMeasuresPerChannelInTDC = %2i ", nTDC, numberOfMeasuresPerChannel, maxMeasuresPerChannelInTDC));
+
+ if (localIndex==-1) continue;
+
+ if (fPackedAcquisition) {
+
+ for (jj=0; jj<=localIndex; jj++) {
+ fIndex++;
+ buf[fIndex] = localBuffer[jj];
+ localBuffer[jj] = 0;
+ psArray[jj] = -1;
+ }
+
+ }
+ else {
+ /*
+ if (maxMeasuresPerChannelInTDC = 1) {
+
+ for (Int_t jj=0; jj<=localIndex; jj++) {
+ if (psArray[jj]==2) {
+ fIndex++;
+ buf[fIndex] = localBuffer[jj];
+ localBuffer[jj] = 0;
+ psArray[jj] = -1;
+ }
+ }
+ for (Int_t jj=0; jj<=localIndex; jj++) {
+ if (psArray[jj]==1) {
+ fIndex++;
+ buf[fIndex] = localBuffer[jj];
+ localBuffer[jj] = 0;
+ psArray[jj] = -1;
+ }
+ }
+
+ } // if (maxMeasuresPerChannelInTDC = 1)
+ else if (maxMeasuresPerChannelInTDC>1) {
+
+ AliDebug(3,Form(" In the TOF DDL %2i, TRM %2i, TDC %2i, chain %1i, the maximum number of t.o.f. good measurements per channel is %2i",
+ nDDL, nTRM, iChain, nTDC, iCH, maxMeasuresPerChannelInTDC));
+ */
+ for (jj=0; jj<=localIndex; jj++) {
+ fIndex++;
+ buf[fIndex] = localBuffer[jj];
+ localBuffer[jj] = 0;
+ psArray[jj] = -1;
+ }
+
+ //} // else if (maxMeasuresPerChannelInTDC>1)
+
+ } // else (!fPackedAcquisition)
+
+ localIndex = -1;
+
+ //maxMeasuresPerChannelInTDC = 0;
+
+ } // end loop on TDC number
+
+
+ if (fVerbose==2) ftxt.close();
+
+}
+
+//----------------------------------------------------------------------------
+void AliTOFDDLRawData::ReverseArray(UInt_t a[], Int_t n) const
+{
+ //
+ // Reverses the n elements of array a
+ //
+
+ Int_t ii, temp;
+
+ for (ii=0; ii<n/2; ii++) {
+ temp = a[ii];
+ a[ii] = a[n-ii-1];
+ a[n-ii-1] = temp;
+ }
+
+ return;
+
+}
+
+//----------------------------------------------------------------------------
+Bool_t AliTOFDDLRawData::HeadOrTail() const
+{
+ //
+ // Returns the result of a 'pitch and toss'
+ //
+
+ Double_t dummy = gRandom->Rndm();
+
+ if (dummy<0.5) return kFALSE;
+ else return kTRUE;
+
+}