X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PMD%2FAliPMDDDLRawData.cxx;h=4140dbacf7e9e98e8dca83e0305d52409950295d;hb=bcd135ec6f4696011285128efcdbb12f6008e7d9;hp=b752b394c27833d704b27c02bc36bdd62735f2bb;hpb=28328eebdab01e366d35993630bc574389417151;p=u%2Fmrichter%2FAliRoot.git diff --git a/PMD/AliPMDDDLRawData.cxx b/PMD/AliPMDDDLRawData.cxx index b752b394c27..4140dbacf7e 100644 --- a/PMD/AliPMDDDLRawData.cxx +++ b/PMD/AliPMDDDLRawData.cxx @@ -18,10 +18,20 @@ #include #include #include +#include +#include +#include "AliLog.h" +#include "AliRawDataHeaderSim.h" +#include "AliBitPacking.h" #include "AliPMDdigit.h" +#include "AliPMDBlockHeader.h" +#include "AliPMDDspHeader.h" +#include "AliPMDPatchBusHeader.h" +#include "AliPMDRawStream.h" #include "AliPMDDDLRawData.h" - +#include "AliDAQ.h" +#include "AliFstream.h" ClassImp(AliPMDDDLRawData) @@ -30,6 +40,24 @@ AliPMDDDLRawData::AliPMDDDLRawData(): { // Default Constructor // + +} +//____________________________________________________________________________ +AliPMDDDLRawData::AliPMDDDLRawData(const AliPMDDDLRawData& ddlraw): + TObject(ddlraw), + fDigits(ddlraw.fDigits) +{ + //Copy Constructor +} +//____________________________________________________________________________ +AliPMDDDLRawData & AliPMDDDLRawData::operator=(const AliPMDDDLRawData& ddlraw) +{ + //Assignment operator + if(this != &ddlraw) + { + fDigits = ddlraw.fDigits; + } + return *this; } //____________________________________________________________________________ @@ -41,437 +69,816 @@ AliPMDDDLRawData::~AliPMDDDLRawData() } //____________________________________________________________________________ -void AliPMDDDLRawData::WritePMDRawData(TTree *treeD, Int_t evtno) +void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) { // write digits into raw data format - ofstream outfile; + AliFstream *outfile; TBranch *branch = treeD->GetBranch("PMDDigit"); - if (!branch) return; + if (!branch) + { + AliError("PMD Digit branch not found"); + return; + } branch->SetAddress(&fDigits); Int_t nmodules = (Int_t) treeD->GetEntries(); - // cout << " nmodules = " << nmodules << endl; + AliDebug(1,Form("Number of modules inside treeD = %d",nmodules)); - const Int_t kSize = 4608; - const Int_t kDDL = 6; + const Int_t kDDL = AliDAQ::NumberOfDdls("PMD"); Int_t modulePerDDL = 0; + AliRawDataHeaderSim header; UInt_t sizeRawData = 0; - UInt_t magicWord = 0x123456; - UInt_t detectorID = 0; - UInt_t ddlID = 0; - Int_t flag = 0; - Int_t version = 1; - UInt_t mHSize = 3*sizeof(UInt_t); - Int_t totword = 0; - - UInt_t miniHeader[3]; + + const Int_t kbusSize = 51; + const Int_t kSize = 1536; UInt_t buffer[kSize]; + UInt_t busPatch[kbusSize][1536]; + + Int_t contentsBus[kbusSize]; + Char_t filename[80]; + // open the ddl file info to know the modules per DDL + TString ddlinfofileName(gSystem->Getenv("ALICE_ROOT")); + ddlinfofileName += "/PMD/PMD_ddl_info.dat"; + + ifstream infileddl; + infileddl.open(ddlinfofileName.Data(), ios::in); // ascii file + if(!infileddl) AliError("Could not read the ddl info file"); + + + Int_t mmodule = 0; + Int_t ddlno; + Int_t modno; + Int_t *modulenoddl = 0x0; + for(Int_t iddl = 0; iddl < kDDL; iddl++) { - sprintf(filename,"Ev%dPMDddl%d.dat",evtno,iddl); -#ifndef __DECCXX - outfile.open(filename,ios::binary); -#else - outfile.open(filename); -#endif - - if (iddl < 4) - { - modulePerDDL = 6; - mmodule = 6*iddl; - detectorID = 0; - } - else if (iddl == 4) - { - modulePerDDL = 12; - mmodule = 24; - detectorID = 1; - } - else if (iddl == 5) - { - modulePerDDL = 12; - mmodule = 36; - detectorID = 1; - } + if (infileddl.eof()) break; + infileddl >> ddlno >> modulePerDDL; - miniHeader[0] = sizeRawData; - PackWord(0,23,magicWord,miniHeader[1]); - PackWord(24,31,detectorID,miniHeader[1]); - ddlID = iddl; - PackWord(0,15,ddlID,miniHeader[2]); - PackWord(16,23,flag,miniHeader[2]); - PackWord(24,31,version,miniHeader[2]); + if (modulePerDDL == 0) continue; + modulenoddl = new Int_t [modulePerDDL]; + for (Int_t im = 0; im < modulePerDDL; im++) + { + infileddl >> modno; + modulenoddl[im] = modno; + } - // Write the Dummy Mini Header into the file - Int_t bHPosition = outfile.tellp(); - outfile.write((char*)(miniHeader),mHSize); - + strcpy(filename,AliDAQ::DdlFileName("PMD",iddl)); + + outfile = new AliFstream(filename); + + // Write the Dummy Data Header into the file + Int_t bHPosition = outfile->Tellp(); + outfile->WriteBuffer((char*)(&header),sizeof(header)); - for(Int_t ium = 0; ium < modulePerDDL; ium++) + for (Int_t ibus = 0; ibus < kbusSize; ibus++) { - - for (Int_t i = 0; i < kSize; i++) + contentsBus[ibus] = 0; + for (Int_t ich = 0; ich < kSize; ich++) { - buffer[i] = 0; + busPatch[ibus][ich] = 0; } + } + + for(Int_t ium = 0; ium < modulePerDDL; ium++) + { + //if (iddl == 4 && ium == 6) mmodule = 42; + // Extract energy deposition per cell and pack it // in a 32-bit word and returns all the total words // per one unit-module - GetUMDigitsData(treeD, mmodule, ium, iddl, totword, buffer); - - outfile.write((char*)buffer,totword*sizeof(UInt_t)); + mmodule = modulenoddl[ium]; + GetUMDigitsData(treeD, mmodule, iddl, modulePerDDL, contentsBus, busPatch); + } + + Int_t ij = 0; + Int_t dsp[10]; + Int_t dspBus[10]; + for (Int_t i = 0; i < 10; i++) + { + dsp[i] = 0; + dspBus[i] = 0; + for (Int_t ibus=0; ibus < 5; ibus++) + { + ij++; + if (contentsBus[ij] > 0) + { + dsp[i] += contentsBus[ij]; + } + dspBus[i]++; + } + // Add the patch Bus header to the DSP contents + dsp[i] += 4*dspBus[i]; + } + + Int_t dspBlockARDL = 0; + Int_t dspBlockBRDL = 0; + Int_t remainder = 0; - mmodule++; + for (Int_t i = 0; i < 5; i++) + { + Int_t ieven = 2*i; + Int_t iodd = 2*i + 1; + if (dsp[ieven] > 0) + { + dspBlockARDL += dsp[ieven]; + remainder = dsp[ieven]%2; + if (remainder == 1) + { + dspBlockARDL++; + } + } + if (dsp[iodd] > 0) + { + dspBlockBRDL += dsp[iodd]; + remainder = dsp[iodd]%2; + if (remainder == 1) + { + dspBlockBRDL++; + } + } } - // Write real mini header - // take the pointer to the beginning of the mini header - // write the total number of words per ddl and bring the - // pointer to the current file position and close it - UInt_t cFPosition = outfile.tellp(); - sizeRawData = cFPosition - bHPosition - mHSize; - outfile.seekp(bHPosition); - outfile.write((char*)(&sizeRawData),sizeof(UInt_t)); - outfile.seekp(cFPosition); + dspBlockARDL += 50; + dspBlockBRDL += 50; - outfile.close(); - } // DDL Loop over + // Start writing the DDL file + AliPMDBlockHeader blHeader; + AliPMDDspHeader dspHeader; + AliPMDPatchBusHeader pbusHeader; -} -//____________________________________________________________________________ + const Int_t kblHLen = blHeader.GetHeaderLength(); + const Int_t kdspHLen = dspHeader.GetHeaderLength(); + const Int_t kpbusHLen = pbusHeader.GetHeaderLength(); -void AliPMDDDLRawData::ReadPMDRawData(Int_t evtno) -{ - // reads the raw data - ifstream infile; + UInt_t dspRDL = 0; + UInt_t dspBlockHeaderWord[8]; + UInt_t dspHeaderWord[10]; + UInt_t patchBusHeaderWord[4]; + Int_t iskip[5]; + UInt_t ddlEndWord[2] = {0xDEADFACE, 0xDEADFACE}; - Char_t filename[80]; + Int_t bknJunk = 0; - const Int_t kDDL = 6; - Int_t imodule = 0; - for(Int_t iddl = 0; iddl < kDDL; iddl++) - { - sprintf(filename,"Ev%dPMDddl%d.dat",evtno,iddl); -#ifndef __DECCXX - infile.open(filename,ios::binary); -#else - infile.open(filename); -#endif - - Int_t ium; - Int_t irow; - Int_t icol; - UInt_t baseword; - UInt_t miniHeader[3]; - UInt_t sizeRawData; - UInt_t mHSize = 3*sizeof(UInt_t); - - infile.read((char*)(miniHeader), mHSize); - - sizeRawData = miniHeader[0]; - Int_t totword = sizeRawData/4; - - UInt_t adc; - UInt_t chno; - UInt_t mcmno; - - for (Int_t iword = 0; iword < totword; iword++) + for (Int_t iblock = 0; iblock < 2; iblock++) { - infile.read((char*)(&baseword),sizeof(UInt_t)); + // DSP Block Header - UnpackWord(0,11,adc,baseword); - UnpackWord(12,17,chno,baseword); - UnpackWord(18,28,mcmno,baseword); - - GetRowCol(iddl, mcmno, chno, ium, irow, icol); - if (iddl < 4) + for (Int_t i=0; iWriteBuffer((char*)dspBlockHeaderWord,kblHLen*sizeof(UInt_t)); + + if (iblock == 0) + { + iskip[0] = 0; + iskip[1] = 10; + iskip[2] = 20; + iskip[3] = 30; + iskip[4] = 40; + } + else if (iblock == 1) + { + iskip[0] = 5; + iskip[1] = 15; + iskip[2] = 25; + iskip[3] = 35; + iskip[4] = 45; + } + + for (Int_t idsp = 0; idsp < 5; idsp++) + { + // DSP Header + Int_t dspno = 0; + if (iblock == 0) + { + dspno = 2*idsp; + dspRDL = (UInt_t) dsp[dspno]; + } + else if (iblock == 1) + { + dspno = 2*idsp + 1; + dspRDL = (UInt_t) dsp[dspno]; + } + + for (Int_t i=0; iWriteBuffer((char*)dspHeaderWord,kdspHLen*sizeof(UInt_t)); + + for (Int_t ibus = 0; ibus < 5; ibus++) + { + // Patch Bus Header + + Int_t busno = iskip[idsp] + ibus + 1; + Int_t patchbusRDL = contentsBus[busno]; + + if (patchbusRDL > 0) + { + patchBusHeaderWord[0] = 0; + patchBusHeaderWord[1] = (UInt_t) (patchbusRDL + kpbusHLen); + patchBusHeaderWord[2] = (UInt_t) patchbusRDL; + patchBusHeaderWord[3] = (UInt_t) busno; + } + else if (patchbusRDL == 0) + { + patchBusHeaderWord[0] = 0; + patchBusHeaderWord[1] = (UInt_t) kpbusHLen; + patchBusHeaderWord[2] = (UInt_t) 0; + patchBusHeaderWord[3] = (UInt_t) busno; + } + + + outfile->WriteBuffer((char*)patchBusHeaderWord,4*sizeof(UInt_t)); + + bknJunk += patchbusRDL; + + for (Int_t iword = 0; iword < patchbusRDL; iword++) + { + buffer[iword] = busPatch[busno][iword]; + } + + outfile->WriteBuffer((char*)buffer,patchbusRDL*sizeof(UInt_t)); + + } // End of patch bus loop + + + // Adding a padding word if the total words odd + if (remainder == 1) + { + UInt_t paddingWord = dspHeader.GetDefaultPaddingWord(); + outfile->WriteBuffer((char*)(&paddingWord),sizeof(UInt_t)); + } + } } - } // ddl loop + // Write two extra word at the end of each DDL file + outfile->WriteBuffer((char*)ddlEndWord,2*sizeof(UInt_t)); + + // Write real data header + // take the pointer to the beginning of the data header + // write the total number of words per ddl and bring the + // pointer to the current file position and close it + UInt_t cFPosition = outfile->Tellp(); + sizeRawData = cFPosition - bHPosition - sizeof(header); + + header.fSize = cFPosition - bHPosition; + header.SetAttribute(0); // valid data + outfile->Seekp(bHPosition); + outfile->WriteBuffer((char*)(&header),sizeof(header)); + outfile->Seekp(cFPosition); + + delete outfile; + } // DDL Loop over + + delete [] modulenoddl; + infileddl.close(); } //____________________________________________________________________________ -void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, Int_t ium, - Int_t ddlno, Int_t & totword, - UInt_t *buffer) +void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, + Int_t ddlno, Int_t modulePerDDL, + Int_t *contentsBus, + UInt_t busPatch[][1536]) { + // Retrieves digits data UnitModule by UnitModule - UInt_t dataword, baseword; - + UInt_t baseword; UInt_t mcmno, chno; UInt_t adc; - Int_t irownew = 0; - Int_t icolnew = 0; Int_t det, smn, irow, icol; + Int_t parity; + + const Int_t kMaxBus = 51; // BKN + Int_t totPatchBus, bPatchBus, ePatchBus; + Int_t ibus, totmcm, rows, cols, rowe, cole; + Int_t moduleno; + Int_t busno = 0; + Int_t patchBusNo[kMaxBus], mcmperBus[kMaxBus]; + Int_t startRowBus[kMaxBus], startColBus[kMaxBus]; + Int_t endRowBus[kMaxBus], endColBus[kMaxBus]; + + Int_t beginPatchBus = -1; + Int_t endPatchBus = -1; + for(Int_t i = 0; i < kMaxBus; i++) + { + patchBusNo[i] = -1; + mcmperBus[i] = -1; + startRowBus[i] = -1; + startColBus[i] = -1; + endRowBus[i] = -1; + endColBus[i] = -1; + } - treeD->GetEntry(imodule); - Int_t nentries = fDigits->GetLast(); - totword = nentries+1; + TString fileName(gSystem->Getenv("ALICE_ROOT")); - for (Int_t ient = 0; ient < totword; ient++) + if(ddlno == 0) + { + fileName += "/PMD/PMD_Mapping_ddl0.dat"; + } + else if(ddlno == 1) + { + fileName += "/PMD/PMD_Mapping_ddl1.dat"; + } + else if(ddlno == 2) + { + fileName += "/PMD/PMD_Mapping_ddl2.dat"; + } + else if(ddlno == 3) + { + fileName += "/PMD/PMD_Mapping_ddl3.dat"; + } + else if(ddlno == 4) { - fPMDdigit = (AliPMDdigit*)fDigits->UncheckedAt(ient); + fileName += "/PMD/PMD_Mapping_ddl4.dat"; + } + else if(ddlno == 5) + { + fileName += "/PMD/PMD_Mapping_ddl5.dat"; + } + + ifstream infile; + infile.open(fileName.Data(), ios::in); // ascii file + if(!infile) + AliError(Form("Could not read the mapping file for DDL No = %d",ddlno)); + + for (Int_t im = 0; im < modulePerDDL; im++) + { + infile >> moduleno; + infile >> totPatchBus >> bPatchBus >> ePatchBus; - det = fPMDdigit->GetDetector(); - smn = fPMDdigit->GetSMNumber(); - irow = fPMDdigit->GetRow(); - icol = fPMDdigit->GetColumn(); - adc = (UInt_t) fPMDdigit->GetADC(); + if (totPatchBus == 0) continue; // BKN - if(smn < 12) + if (moduleno == imodule) { - irownew = icol; - icolnew = irow; + beginPatchBus = bPatchBus; + endPatchBus = ePatchBus; } - else if( smn >= 12 && smn < 24) + + for(Int_t i=0; i> ibus >> totmcm >> rows >> rowe >> cols >> cole; + + if (moduleno == imodule) + { + patchBusNo[ibus] = ibus; + mcmperBus[ibus] = totmcm; + startRowBus[ibus] = rows; + startColBus[ibus] = cols; + endRowBus[ibus] = rowe; + endColBus[ibus] = cole; + } } + } + infile.close(); - GetMCMCh(ddlno, ium, irownew, icolnew, mcmno, chno); + // Read if some chains are off + TString rchainName(gSystem->Getenv("ALICE_ROOT")); + rchainName += "/PMD/PMD_removed_chains.dat"; - fprintf(fpw2,"%d %d %d %d %d %d\n",imodule,irownew,icolnew,mcmno,chno,adc); + ifstream rchainfile; + rchainfile.open(rchainName.Data(), ios::in); // ascii file + if(!rchainfile)AliError("Could not read the removed cahins file"); - baseword = 0; - dataword = adc; - PackWord(0, 11,dataword,baseword); - dataword = chno; - PackWord(12,17,dataword,baseword); - dataword = mcmno; - PackWord(18,28,dataword,baseword); - dataword = 0; - PackWord(29,29,dataword,baseword); - dataword = 0; - PackWord(30,30,dataword,baseword); - dataword = 1; - PackWord(31,31,dataword,baseword); - - buffer[ient] = baseword; - + Int_t srowoff1[2][24], erowoff1[2][24]; + Int_t scoloff1[2][24], ecoloff1[2][24]; + Int_t srowoff2[2][24], erowoff2[2][24]; + Int_t scoloff2[2][24], ecoloff2[2][24]; + Int_t rows1, rowe1, cols1, cole1; + Int_t rows2, rowe2, cols2, cole2; - // fprintf(fpw1,"%d %d %d %d %d\n",det,smn,irow,icol,adc); + for (Int_t im = 0; im < 48; im++) + { + rchainfile >> det >> smn >> rows1 >> rowe1 >> cols1 >> cole1 + >> rows2 >> rowe2 >> cols2 >> cole2; + srowoff1[det][smn] = rows1; + erowoff1[det][smn] = rowe1; + scoloff1[det][smn] = cols1; + ecoloff1[det][smn] = cole1; + srowoff2[det][smn] = rows2; + erowoff2[det][smn] = rowe2; + scoloff2[det][smn] = cols2; + ecoloff2[det][smn] = cole2; } + rchainfile.close(); -} -//____________________________________________________________________________ -void AliPMDDDLRawData::GetMCMCh(Int_t ddlno, Int_t um, - Int_t row, Int_t col, - UInt_t &mcmno, UInt_t &chno) -{ - // This part will be modified once the final track layout on the PCB is - // designed. This will only change the coordinate of the individual cell - UInt_t ch[16][4] = { {3, 2, 29, 28}, - {5, 1, 31, 27}, - {6, 0, 30, 26}, - {7, 4, 25, 24}, - {8, 9, 20, 23}, - {10, 14, 16, 22}, - {11, 15, 17, 21}, - {12, 13, 18, 19}, - {35, 34, 61, 60}, - {37, 33, 63, 59}, - {38, 32, 62, 58}, - {39, 36, 57, 56}, - {40, 41, 52, 55}, - {42, 46, 48, 54}, - {43, 47, 49, 53}, - {44, 45, 50, 51} }; - - if (ddlno == 0 || ddlno == 1) - { - // PRE plane, SU Mod = 1, 2 - Int_t irownew = row%16; - Int_t icolnew = col%4; - Int_t irowdiv = row/16; - Int_t icoldiv = col/4; - - mcmno = 72*um + 6*icoldiv + irowdiv; - chno = ch[irownew][icolnew]; - } - else if (ddlno == 2 || ddlno == 3) + treeD->GetEntry(imodule); + Int_t nentries = fDigits->GetLast(); + Int_t totword = nentries+1; + + AliPMDdigit *pmddigit; + + for (Int_t ient = 0; ient < totword; ient++) { - // PRE plane, SU Mod = 3, 4 - Int_t irownew = row%16; - Int_t icolnew = col%4; - Int_t irowdiv = row/16; - Int_t icoldiv = col/4; + pmddigit = (AliPMDdigit*)fDigits->UncheckedAt(ient); - mcmno = 72*um + 3*icoldiv + irowdiv; - chno = ch[irownew][icolnew]; - } - else if (ddlno == 4) - { - // CPV plane, SU Mod = 1, 2 - Int_t irownew = row%16; - Int_t icolnew = col%4; - Int_t irowdiv = row/16; - Int_t icoldiv = col/4; - - mcmno = 72*um + 6*icoldiv + irowdiv; - chno = ch[irownew][icolnew]; - } - else if (ddlno == 5) - { - // CPV plane, SU Mod = 3, 4 - Int_t irownew = row%16; - Int_t icolnew = col%4; - Int_t irowdiv = row/16; - Int_t icoldiv = col/4; + det = pmddigit->GetDetector(); + smn = pmddigit->GetSMNumber(); + irow = pmddigit->GetRow(); + icol = pmddigit->GetColumn(); + adc = (UInt_t) pmddigit->GetADC(); + + TransformS2H(smn,irow,icol); - mcmno = 72*um + 3*icoldiv + irowdiv; - chno = ch[irownew][icolnew]; + // remove the non-existence channels + + //printf("%d %d %d %d\n",det,smn,irow,icol); + //printf("--- %d %d %d %d\n",srowoff[det][smn],erowoff[det][smn], + // scoloff[det][smn],ecoloff[det][smn]); + + if (irow >= srowoff1[det][smn] && irow <= erowoff1[det][smn]) + { + if (icol >= scoloff1[det][smn] && icol <= ecoloff1[det][smn]) + { + continue; + } + } + if (irow >= srowoff2[det][smn] && irow <= erowoff2[det][smn]) + { + if (icol >= scoloff2[det][smn] && icol <= ecoloff2[det][smn]) + { + continue; + } + } + + + GetMCMCh(imodule, irow, icol, beginPatchBus, endPatchBus, + mcmperBus, startRowBus, startColBus, + endRowBus, endColBus, busno, mcmno, chno); + + baseword = 0; + AliBitPacking::PackWord(adc,baseword,0,11); + AliBitPacking::PackWord(chno,baseword,12,17); + AliBitPacking::PackWord(mcmno,baseword,18,28); + AliBitPacking::PackWord(0,baseword,29,30); + parity = ComputeParity(baseword); // generate the parity bit + AliBitPacking::PackWord(parity,baseword,31,31); + + Int_t jj = contentsBus[busno]; + busPatch[busno][jj] = baseword; + + contentsBus[busno]++; } } -//____________________________________________________________________________ -void AliPMDDDLRawData::GetRowCol(Int_t ddlno, UInt_t mcmno, UInt_t chno, - Int_t &um, Int_t &row, Int_t &col) +//____________________________________________________________________________ +void AliPMDDDLRawData::TransformS2H(Int_t smn, Int_t &irow, Int_t &icol) { - UInt_t ch[64] = { 9, 5, 1, 0, 13, 4, 8, 12, - 16, 17, 20, 24, 28, 29, 21, 25, - 22, 26, 30, 31, 18, 27, 23, 19, - 15, 14, 11, 7, 3, 2, 10, 6, - 41, 37, 33, 32, 45, 36, 40, 44, - 48, 49, 52, 56, 60, 61, 53, 57, - 54, 58, 62, 63, 50, 59, 55, 51, - 47, 46, 43, 39, 35, 34, 42, 38 }; - - if (ddlno == 1 || ddlno == 2) - { - um = mcmno/72; - Int_t mcmnonew = mcmno - 72*um; - Int_t rowcol = ch[chno]; - Int_t irownew = rowcol/4; - Int_t icolnew = rowcol%4; - - Int_t remmcm = mcmnonew%6; - Int_t divmcm = mcmnonew/6; - - row = 16*remmcm + irownew; - col = 4*divmcm + icolnew; - } - else if (ddlno == 3 || ddlno == 4) + // Does the Software to Hardware coordinate transformation + // + + Int_t irownew = 0; + Int_t icolnew = 0; + + // First in digits we have all dimension 48x96 + // Transform into the realistic one, i.e, For SM 0&1 96(row)x48(col) + // and for SM 2&3 48(row)x96(col) + // + if(smn < 12) { - um = mcmno/72; - Int_t mcmnonew = mcmno - 72*um; - Int_t rowcol = ch[chno]; - Int_t irownew = rowcol/4; - Int_t icolnew = rowcol%4; - - Int_t remmcm = mcmnonew%3; - Int_t divmcm = mcmnonew/3; - - row = 16*remmcm + irownew; - col = 4*divmcm + icolnew; + irownew = icol; + icolnew = irow; } - else if (ddlno == 4) + else if( smn >= 12 && smn < 24) { - um = mcmno/144; - Int_t mcmnonew = mcmno - 72*um; - Int_t rowcol = ch[chno]; - Int_t irownew = rowcol/4; - Int_t icolnew = rowcol%4; - - Int_t remmcm = mcmnonew%6; - Int_t divmcm = mcmnonew/6; - - row = 16*remmcm + irownew; - col = 4*divmcm + icolnew; - } - else if (ddlno == 5) - { - um = mcmno/144; - Int_t mcmnonew = mcmno - 72*um; - Int_t rowcol = ch[chno]; - Int_t irownew = rowcol/4; - Int_t icolnew = rowcol%4; - - Int_t remmcm = mcmnonew%3; - Int_t divmcm = mcmnonew/3; - - row = 16*remmcm + irownew; - col = 4*divmcm + icolnew; + irownew = irow; + icolnew = icol; } + irow = irownew; + icol = icolnew; } + + //____________________________________________________________________________ -void AliPMDDDLRawData::PackWord(UInt_t startbit, UInt_t stopbit, - UInt_t dataword, UInt_t &packedword) +void AliPMDDDLRawData::GetMCMCh(Int_t imodule, Int_t row, Int_t col, + Int_t beginPatchBus, Int_t endPatchBus, + Int_t *mcmperBus, + Int_t *startRowBus, Int_t *startColBus, + Int_t *endRowBus, Int_t *endColBus, + Int_t & busno, UInt_t &mcmno, UInt_t &chno) { - UInt_t bitLength = stopbit - startbit + 1; - UInt_t bitContent = (UInt_t) (TMath::Power(2,bitLength) - 1); - if(bitContent < dataword) + // This converts row col to hardware channel number + // This is the final version of mapping supplied by Mriganka + + UInt_t iCh[16][4]; + + static const UInt_t kChDdl01[16][4] = { {6, 4, 5, 7}, + {10, 2, 1, 9}, + {12, 0, 3, 11}, + {14, 8, 13, 15}, + {16, 18, 23, 17}, + {20, 28, 31, 19}, + {22, 30, 29, 21}, + {24, 26, 27, 25}, + {38, 36, 37, 39}, + {42, 34, 33, 41}, + {44, 32, 35, 43}, + {46, 40, 45, 47}, + {48, 50, 55, 49}, + {52, 60, 63, 51}, + {54, 62, 61, 53}, + {56, 58, 59, 57} }; + + + static const UInt_t kChDdl23[16][4] = { {57, 59, 58, 56}, + {53, 61, 62, 54}, + {51, 63, 60, 52}, + {49, 55, 50, 48}, + {47, 45, 40, 46}, + {43, 35, 32, 44}, + {41, 33, 34, 42}, + {39, 37, 36, 38}, + {25, 27, 26, 24}, + {21, 29, 30, 22}, + {19, 31, 28, 20}, + {17, 23, 18, 16}, + {15, 13, 8, 14}, + {11, 3, 0, 12}, + {9, 1, 2, 10}, + {7, 5, 4, 6} }; + + + static const UInt_t kChDdl41[16][4] = { {56, 58, 59, 57}, + {54, 62, 61, 53}, + {52, 60, 63, 51}, + {48, 50, 55, 49}, + {46, 40, 45, 47}, + {44, 32, 35, 43}, + {42, 34, 33, 41}, + {38, 36, 37, 39}, + {24, 26, 27, 25}, + {22, 30, 29, 21}, + {20, 28, 31, 19}, + {16, 18, 23, 17}, + {14, 8, 13, 15}, + {12, 0, 3, 11}, + {10, 2, 1, 9}, + {6, 4, 5, 7} }; + + + static const UInt_t kChDdl42[16][4] = { {7, 5, 4, 6}, + {9, 1, 2, 10}, + {11, 3, 0, 12}, + {15, 13, 8, 14}, + {17, 23, 18, 16}, + {19, 31, 28, 20}, + {21, 29, 30, 22}, + {25, 27, 26, 24}, + {39, 37, 36, 38}, + {41, 33, 34, 42}, + {43, 35, 32, 44}, + {47, 45, 40, 46}, + {49, 55, 50, 48}, + {51, 63, 60, 52}, + {53, 61, 62, 54}, + {57, 59, 58, 56} }; + + + static const UInt_t kChDdl51[16][4] = { {7, 5, 4, 6}, + {9, 1, 2, 10}, + {11, 3, 0, 12}, + {15, 13, 8, 14}, + {17, 23, 18, 16}, + {19, 31, 28, 20}, + {21, 29, 30, 22}, + {25, 27, 26, 24}, + {39, 37, 36, 38}, + {41, 33, 34, 42}, + {43, 35, 32, 44}, + {47, 45, 40, 46}, + {49, 55, 50, 48}, + {51, 63, 60, 52}, + {53, 61, 62, 54}, + {57, 59, 58, 56} }; + + + + static const UInt_t kChDdl52[16][4] = { {56, 58, 59, 57}, + {54, 62, 61, 53}, + {52, 60, 63, 51}, + {48, 50, 55, 49}, + {46, 40, 45, 47}, + {44, 32, 35, 43}, + {42, 34, 33, 41}, + {38, 36, 37, 39}, + {24, 26, 27, 25}, + {22, 30, 29, 21}, + {20, 28, 31, 19}, + {16, 18, 23, 17}, + {14, 8, 13, 15}, + {12, 0, 3, 11}, + {10, 2, 1, 9}, + {6, 4, 5, 7} }; + + + for (Int_t i = 0; i < 16; i++) + { + for (Int_t j = 0; j < 4; j++) + { + + if(imodule < 6) iCh[i][j] = kChDdl01[i][j]; + if(imodule >= 6 && imodule <= 11) iCh[i][j] = kChDdl01[i][j]; + if(imodule >= 12 && imodule <= 17) iCh[i][j] = kChDdl23[i][j]; + if(imodule >= 18 && imodule <= 23) iCh[i][j] = kChDdl23[i][j]; + + if(imodule >= 24 && imodule <= 29) iCh[i][j] = kChDdl41[i][j]; + if(imodule >= 42 && imodule <= 47) iCh[i][j] = kChDdl42[i][j]; + if(imodule >= 36 && imodule <= 41) iCh[i][j] = kChDdl51[i][j]; + if(imodule >= 30 && imodule <= 35) iCh[i][j] = kChDdl52[i][j]; + + } + } + + + Int_t irownew = row%16; + Int_t icolnew = col%4; + + chno = iCh[irownew][icolnew]; + + + for (Int_t ibus = beginPatchBus; ibus <= endPatchBus; ibus++) { - cout << " *** ERROR *** bitContent is less than the dataword" << endl; - return; - } - UInt_t packedBits = 0; - if (packedword != 0) - packedBits = (UInt_t) (TMath::Log(packedword)/TMath::Log(2)); + Int_t srow = startRowBus[ibus]; + Int_t erow = endRowBus[ibus]; + Int_t scol = startColBus[ibus]; + Int_t ecol = endColBus[ibus]; + Int_t tmcm = mcmperBus[ibus]; - UInt_t counter; - if (packedBits <= stopbit) - { - counter = 31 - stopbit; - } - else - { - counter = 31 - packedBits; - } - UInt_t dummyword = 0xFFFFFFFF; - dummyword >>= counter; - UInt_t lword = dataword << startbit; - UInt_t nword = lword | packedword; - packedword = dummyword & nword; + if ((row >= srow && row <= erow) && (col >= scol && col <= ecol)) + { + busno = ibus; + + // Find out the MCM Number + // + if (imodule < 6) mcmno = (col-scol)/4 + 1; + if (imodule >= 6 && imodule < 12) mcmno = (col-scol)/4 + 1; + + if (imodule >= 12 && imodule < 18) + { + icolnew = (col - scol)/4; + mcmno = tmcm - icolnew; + } + if (imodule >= 18 && imodule < 24) + { + icolnew = (col - scol)/4; + mcmno = tmcm - icolnew; + } + + // DDL = 4 + if (imodule >= 24 && imodule < 30) + { + + //if (tmcm == 24) + Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; + if(rowdiff > 16) + { + Int_t midrow = srow + 16; + if(row >= srow && row < midrow) + { + mcmno = 12 + (col-scol)/4 + 1; + } + else if(row >= midrow && row <= erow) + + { + mcmno = (col-scol)/4 + 1; + } + } + else if (rowdiff < 16) + { + mcmno = (col-scol)/4 + 1; + } + + } + if (imodule >= 42 && imodule < 48) + { + Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; + if(rowdiff > 16) + { + Int_t midrow = srow + 16; + if (row >= midrow && row <= erow) + { + mcmno = 12 + (ecol -col)/4 + 1; + } + else if (row >= srow && row < midrow) + { + mcmno = (ecol - col)/4 + 1; + } + } + else if (rowdiff < 16) + { + mcmno = (ecol - col)/4 + 1; + } + } + + // DDL = 5 + if (imodule >= 30 && imodule < 36) + { + // CPV plane, SU Mod = 1, 2 : ddl = 5 + + //if(tmcm == 24) + Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; + if(rowdiff > 16) + { + Int_t midrow = srow + 16; + if(row >= srow && row < midrow) + { + mcmno = 12 + (col-scol)/4 + 1; + } + else if(row >= midrow && row <= erow) + { + mcmno = (col-scol)/4 + 1; + } + } + else if(rowdiff < 16) + { + mcmno = (col-scol)/4 + 1; + } + + } + if (imodule >= 36 && imodule < 42) + { + Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; + if(rowdiff > 16) + { + Int_t midrow = srow + 16; + if (row >= midrow && row <= erow) + { + mcmno = 12 + (ecol - col)/4 + 1; + } + else if (row >= srow && row < midrow) + { + mcmno = (ecol - col)/4 + 1; + } + } + else if (rowdiff < 16) + { + mcmno = (ecol - col)/4 + 1; + } + } + + } + } +} -} //____________________________________________________________________________ -void AliPMDDDLRawData::UnpackWord(UInt_t startbit, UInt_t stopbit, - UInt_t &dataword, UInt_t packedword) + +Int_t AliPMDDDLRawData::ComputeParity(UInt_t baseword) { - UInt_t bitLength = stopbit - startbit + 1; - UInt_t bitContent = (UInt_t) (TMath::Power(2,bitLength) - 1); - bitContent <<= startbit; - dataword = packedword & bitContent; - dataword >>= startbit; + // Generate the parity bit + Int_t count = 0; + for(Int_t j=0; j<29; j++) + { + if (baseword & 0x01 ) count++; + baseword >>= 1; + } + Int_t parity = count%2; + return parity; } -//____________________________________________________________________________ +//____________________________________________________________________________