X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PMD%2FAliPMDDDLRawData.cxx;h=faca049fd0debbf11ca96ceed6b59c381c5121fe;hb=ed5512107390ee1d6b033403dcd839093b8069d0;hp=32aaae91f1e417c9e7606338c4f0107917df74e9;hpb=08f92f1423b62a716d143535edd275e8b020d76a;p=u%2Fmrichter%2FAliRoot.git diff --git a/PMD/AliPMDDDLRawData.cxx b/PMD/AliPMDDDLRawData.cxx index 32aaae91f1e..faca049fd0d 100644 --- a/PMD/AliPMDDDLRawData.cxx +++ b/PMD/AliPMDDDLRawData.cxx @@ -22,20 +22,27 @@ #include #include "AliLog.h" -#include "AliRawDataHeader.h" +#include "AliRawDataHeaderSim.h" #include "AliBitPacking.h" #include "AliPMDdigit.h" #include "AliPMDBlockHeader.h" #include "AliPMDDspHeader.h" #include "AliPMDPatchBusHeader.h" #include "AliPMDRawStream.h" +#include "AliPMDddlinfoData.h" +#include "AliPMDMappingData.h" #include "AliPMDDDLRawData.h" #include "AliDAQ.h" #include "AliFstream.h" +#include "AliCDBManager.h" +#include "AliCDBStorage.h" +#include "AliCDBEntry.h" ClassImp(AliPMDDDLRawData) AliPMDDDLRawData::AliPMDDDLRawData(): + fDdlinfo(GetDdlinfoData()), + fMapData(GetMappingData()), fDigits(new TClonesArray("AliPMDdigit", 1000)) { // Default Constructor @@ -45,6 +52,8 @@ AliPMDDDLRawData::AliPMDDDLRawData(): //____________________________________________________________________________ AliPMDDDLRawData::AliPMDDDLRawData(const AliPMDDDLRawData& ddlraw): TObject(ddlraw), + fDdlinfo(ddlraw.fDdlinfo), + fMapData(ddlraw.fMapData), fDigits(ddlraw.fDigits) { //Copy Constructor @@ -55,7 +64,9 @@ AliPMDDDLRawData & AliPMDDDLRawData::operator=(const AliPMDDDLRawData& ddlraw) //Assignment operator if(this != &ddlraw) { - fDigits = ddlraw.fDigits; + fDdlinfo = ddlraw.fDdlinfo; + fMapData = ddlraw.fMapData; + fDigits = ddlraw.fDigits; } return *this; } @@ -87,74 +98,64 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) AliDebug(1,Form("Number of modules inside treeD = %d",nmodules)); const Int_t kDDL = AliDAQ::NumberOfDdls("PMD"); - Int_t modulePerDDL = 0; - AliRawDataHeader header; + AliRawDataHeaderSim header; UInt_t sizeRawData = 0; + const Int_t kbusSize = 51; const Int_t kSize = 1536; UInt_t buffer[kSize]; - UInt_t busPatch[50][1536]; + UInt_t busPatch[kbusSize][1536]; - Int_t contentsBus[50]; + Int_t contentsBus[kbusSize]; Char_t filename[80]; - + Int_t modulePerDDL = 0; Int_t mmodule = 0; + Int_t ddlno; + Int_t modulenoddl[12] = {-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1}; + for(Int_t iddl = 0; iddl < kDDL; iddl++) { - strcpy(filename,AliDAQ::DdlFileName("PMD",iddl)); - - outfile = new AliFstream(filename); - - if (iddl < 4) - { - modulePerDDL = 6; - mmodule = 6*iddl; - } - else if (iddl == 4) - { - modulePerDDL = 12; - mmodule = 24; - } - else if (iddl == 5) + ddlno = iddl; + modulePerDDL = fDdlinfo->GetNoOfModulePerDdl(iddl); + if (modulePerDDL == 0) continue; + for (Int_t im = 0; im < 12; im++) { - modulePerDDL = 12; - mmodule = 30; + modulenoddl[im] = fDdlinfo->GetModulesPerDdl(iddl,im);; } - - + strncpy(filename,AliDAQ::DdlFileName("PMD",iddl),79); + + 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 ibus = 0; ibus < 50; ibus++) + for (Int_t ibus = 0; ibus < kbusSize; ibus++) { contentsBus[ibus] = 0; - for (Int_t ich = 0; ich < 1536; ich++) + for (Int_t ich = 0; ich < kSize; ich++) { busPatch[ibus][ich] = 0; } } - for(Int_t ium = 0; ium < modulePerDDL; ium++) + for(Int_t ium = 0; ium < 12; 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 + mmodule = modulenoddl[ium]; + if(mmodule == -1) continue; GetUMDigitsData(treeD, mmodule, iddl, contentsBus, busPatch); - mmodule++; } - - Int_t ij = 0; Int_t dsp[10]; Int_t dspBus[10]; @@ -164,12 +165,12 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) dspBus[i] = 0; for (Int_t ibus=0; ibus < 5; ibus++) { + ij++; if (contentsBus[ij] > 0) { dsp[i] += contentsBus[ij]; - dspBus[i]++; } - ij++; + dspBus[i]++; } // Add the patch Bus header to the DSP contents dsp[i] += 4*dspBus[i]; @@ -177,6 +178,8 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) Int_t dspBlockARDL = 0; Int_t dspBlockBRDL = 0; + Int_t remainder = 0; + for (Int_t i = 0; i < 5; i++) { @@ -185,7 +188,7 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) if (dsp[ieven] > 0) { dspBlockARDL += dsp[ieven]; - Int_t remainder = dsp[ieven]%2; + remainder = dsp[ieven]%2; if (remainder == 1) { dspBlockARDL++; @@ -194,14 +197,17 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) if (dsp[iodd] > 0) { dspBlockBRDL += dsp[iodd]; - Int_t remainder = dsp[ieven]%2; + remainder = dsp[iodd]%2; if (remainder == 1) { - dspBlockARDL++; + dspBlockBRDL++; } } } - + + dspBlockARDL += 50; + dspBlockBRDL += 50; + // Start writing the DDL file AliPMDBlockHeader blHeader; @@ -216,8 +222,8 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) UInt_t dspBlockHeaderWord[8]; UInt_t dspHeaderWord[10]; UInt_t patchBusHeaderWord[4]; - Int_t iskip[5]; - + Int_t iskip[5]; + UInt_t ddlEndWord[2] = {0xDEADFACE, 0xDEADFACE}; for (Int_t iblock = 0; iblock < 2; iblock++) { @@ -276,19 +282,22 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) { dspHeaderWord[i] = 0; } - Int_t remainder = dspRDL%2; + remainder = dspRDL%2; if (remainder == 1) dspRDL++; dspHeaderWord[1] = dspRDL + kdspHLen; dspHeaderWord[2] = dspRDL; dspHeaderWord[3] = dspno; if (remainder == 1) dspHeaderWord[8] = 1; // setting the padding word + + outfile->WriteBuffer((char*)dspHeaderWord,kdspHLen*sizeof(UInt_t)); for (Int_t ibus = 0; ibus < 5; ibus++) { // Patch Bus Header - Int_t busno = iskip[idsp] + ibus; + + Int_t busno = iskip[idsp] + ibus + 1; Int_t patchbusRDL = contentsBus[busno]; if (patchbusRDL > 0) @@ -309,7 +318,6 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) outfile->WriteBuffer((char*)patchBusHeaderWord,4*sizeof(UInt_t)); - for (Int_t iword = 0; iword < patchbusRDL; iword++) { buffer[iword] = busPatch[busno][iword]; @@ -329,6 +337,9 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) } } + // 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 @@ -345,25 +356,23 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) delete outfile; } // DDL Loop over - } //____________________________________________________________________________ + void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, - Int_t ddlno, Int_t *contentsBus, + Int_t ddlno, Int_t *contentsBus, UInt_t busPatch[][1536]) { - // Retrives digits data UnitModule by UnitModule + // Retrieves digits data UnitModule by UnitModule - UInt_t baseword; - UInt_t mcmno, chno; - UInt_t adc; - Int_t det, smn, irow, icol; - Int_t parity; + const Int_t kMaxBus = 51; + + UInt_t baseword = 0; + UInt_t mcmno = 0, chno = 0; + UInt_t adc = 0; + Int_t det = 0, smn = 0, irow = 0, icol = 0; + UInt_t parity = 0; - const Int_t kMaxBus = 50; - 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]; @@ -380,83 +389,40 @@ void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, endRowBus[i] = -1; endColBus[i] = -1; } - Int_t modulePerDDL = 0; - if (ddlno < 4) - { - modulePerDDL = 6; - } - else if (ddlno == 4 || ddlno == 5) - { - modulePerDDL = 12; - } - TString fileName(gSystem->Getenv("ALICE_ROOT")); + // Fetch the DDL mapping info from the mapping database - 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) - { - fileName += "/PMD/PMD_Mapping_ddl4.dat"; - } - else if(ddlno == 5) - { - fileName += "/PMD/PMD_Mapping_ddl5.dat"; - } + DdlMapping(ddlno, imodule, beginPatchBus, endPatchBus, + patchBusNo, mcmperBus, startRowBus, endRowBus, + startColBus, endColBus); - 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)); + // Read if some chains are off from the ddlinfo database - for (Int_t im = 0; im < modulePerDDL; im++) - { - infile >> moduleno; - infile >> totPatchBus >> bPatchBus >> ePatchBus; + 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]; - if (moduleno == imodule) - { - beginPatchBus = bPatchBus; - endPatchBus = ePatchBus; - } - - 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; - } + srowoff1[idet][im] = fDdlinfo->GetStartRowA(idet,im); + erowoff1[idet][im] = fDdlinfo->GetEndRowA(idet,im); + scoloff1[idet][im] = fDdlinfo->GetStartColA(idet,im); + ecoloff1[idet][im] = fDdlinfo->GetEndColA(idet,im); + srowoff2[idet][im] = fDdlinfo->GetStartRowB(idet,im); + erowoff2[idet][im] = fDdlinfo->GetEndRowB(idet,im); + scoloff2[idet][im] = fDdlinfo->GetStartColB(idet,im); + ecoloff2[idet][im] = fDdlinfo->GetEndColB(idet,im); } - } - infile.close(); - treeD->GetEntry(imodule); Int_t nentries = fDigits->GetLast(); Int_t totword = nentries+1; - AliPMDdigit *pmddigit; + AliPMDdigit *pmddigit = 0x0; for (Int_t ient = 0; ient < totword; ient++) { @@ -466,11 +432,35 @@ void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, smn = pmddigit->GetSMNumber(); irow = pmddigit->GetRow(); icol = pmddigit->GetColumn(); - adc = (UInt_t) pmddigit->GetADC(); + Float_t aadc = pmddigit->GetADC(); + if (aadc < 0.) aadc = 0.; + adc = (UInt_t) aadc; TransformS2H(smn,irow,icol); + + // 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(ddlno, irow, icol, beginPatchBus, endPatchBus, + + GetMCMCh(imodule, irow, icol, beginPatchBus, endPatchBus, mcmperBus, startRowBus, startColBus, endRowBus, endColBus, busno, mcmno, chno); @@ -522,7 +512,7 @@ void AliPMDDDLRawData::TransformS2H(Int_t smn, Int_t &irow, Int_t &icol) //____________________________________________________________________________ -void AliPMDDDLRawData::GetMCMCh(Int_t ddlno, Int_t row, Int_t col, +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, @@ -532,27 +522,140 @@ void AliPMDDDLRawData::GetMCMCh(Int_t ddlno, Int_t row, Int_t col, // This converts row col to hardware channel number // This is the final version of mapping supplied by Mriganka - static const UInt_t kCh[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} }; - + 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 = kCh[irownew][icolnew]; + chno = iCh[irownew][icolnew]; for (Int_t ibus = beginPatchBus; ibus <= endPatchBus; ibus++) @@ -562,120 +665,193 @@ void AliPMDDDLRawData::GetMCMCh(Int_t ddlno, Int_t row, Int_t col, Int_t scol = startColBus[ibus]; Int_t ecol = endColBus[ibus]; Int_t tmcm = mcmperBus[ibus]; + if ((row >= srow && row <= erow) && (col >= scol && col <= ecol)) { busno = ibus; // Find out the MCM Number // - - if (ddlno == 0 || ddlno == 1) + + if (imodule < 6) mcmno = (col-scol)/4 + 1; + if (imodule >= 6 && imodule < 12) mcmno = (col-scol)/4 + 1; + + if (imodule >= 12 && imodule < 18) { - // PRE plane, SU Mod = 0, 1 - mcmno = (col-scol)/4; - - } - else if (ddlno == 2 || ddlno == 3) + icolnew = (col - scol)/4; + mcmno = tmcm - icolnew; + } + if (imodule >= 18 && imodule < 24) { - // PRE plane, SU Mod = 2, 3 - Int_t icolnew = (col - scol)/4; - mcmno = tmcm - 1 - icolnew; - } - else if (ddlno == 4 ) + icolnew = (col - scol)/4; + mcmno = tmcm - icolnew; + } + + // DDL = 4 + if (imodule >= 24 && imodule < 30) { - // CPV plane, SU Mod = 0, 3 : ddl = 4 - - if(ibus <= 17) + + //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; + mcmno = 12 + (col-scol)/4 + 1; } else if(row >= midrow && row <= erow) - + { - mcmno = (col-scol)/4; + mcmno = (col-scol)/4 + 1; } } - - else if (ibus > 17) + else if (rowdiff < 16) { - Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; - 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) { - Int_t midrow = srow + 16; - if (row >= midrow && row <= erow) - { - mcmno = 12 + (ecol -col)/4; - } - else if (row >= srow && row < midrow) - { - mcmno = (ecol - col)/4; - } + mcmno = 12 + (ecol -col)/4 + 1; } - else if (rowdiff < 16) + else if (row >= srow && row < midrow) { - mcmno = (ecol - col)/4; - } + mcmno = (ecol - col)/4 + 1; + } } - } - else if ( ddlno == 5) + else if (rowdiff < 16) + { + mcmno = (ecol - col)/4 + 1; + } + } + + // DDL = 5 + if (imodule >= 30 && imodule < 36) { - // CPV plane, SU Mod = 0, 3 : ddl = 4 + // CPV plane, SU Mod = 1, 2 : ddl = 5 - if(ibus <= 17) + //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; + mcmno = 12 + (col-scol)/4 + 1; } else if(row >= midrow && row <= erow) { - mcmno = (col-scol)/4; + mcmno = (col-scol)/4 + 1; } } + else if(rowdiff < 16) + { + mcmno = (col-scol)/4 + 1; + } - else if (ibus > 17) + } + if (imodule >= 36 && imodule < 42) + { + Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; + if(rowdiff > 16) { - Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; - if(rowdiff > 16) + Int_t midrow = srow + 16; + if (row >= midrow && row <= erow) { - Int_t midrow = srow + 16; - if (row >= midrow && row <= erow) - { - mcmno = 12 + (ecol -col)/4; - } - else if (row >= srow && row < midrow) - { - mcmno = (ecol - col)/4; - } + mcmno = 12 + (ecol - col)/4 + 1; } - else if (rowdiff < 16) + else if (row >= srow && row < midrow) { - mcmno = (ecol - col)/4; - } + mcmno = (ecol - col)/4 + 1; + } + } + else if (rowdiff < 16) + { + mcmno = (ecol - col)/4 + 1; } } + } } } //____________________________________________________________________________ -Int_t AliPMDDDLRawData::ComputeParity(UInt_t baseword) +UInt_t AliPMDDDLRawData::ComputeParity(UInt_t baseword) { // Generate the parity bit - Int_t count = 0; + UInt_t count = 0; for(Int_t j=0; j<29; j++) { if (baseword & 0x01 ) count++; baseword >>= 1; } - Int_t parity = count%2; + UInt_t parity = count%2; return parity; } +//____________________________________________________________________________ +void AliPMDDDLRawData::DdlMapping(Int_t iddl, Int_t imodule, + Int_t &beginPatchBus, Int_t &endPatchBus, + Int_t patchBusNo[], Int_t mcmperBus[], + Int_t startRowBus[], Int_t endRowBus[], + Int_t startColBus[], Int_t endColBus[]) +{ + // DDL Mapping fetching from mapping database + + beginPatchBus = fMapData->GetBeginPatchBus(iddl,imodule); + endPatchBus = fMapData->GetEndPatchBus(iddl,imodule); + + for(Int_t ibus = beginPatchBus; ibus < endPatchBus+1; ibus++) + { + patchBusNo[ibus] = ibus; + mcmperBus[ibus] = fMapData->GetMcmperBus(iddl,ibus); + startRowBus[ibus] = fMapData->GetStartRowBus(iddl,ibus); + startColBus[ibus] = fMapData->GetStartColBus(iddl,ibus); + endRowBus[ibus] = fMapData->GetEndRowBus(iddl,ibus); + endColBus[ibus] = fMapData->GetEndColBus(iddl,ibus); + } + +} +//____________________________________________________________________________ +AliPMDddlinfoData* AliPMDDDLRawData::GetDdlinfoData() const +{ + // The run number will be centralized in AliCDBManager, + // you don't need to set it here! + AliCDBEntry *entry = AliCDBManager::Instance()->Get("PMD/Calib/Ddlinfo"); + + if(!entry) AliFatal("ddlinfo object retrieval failed!"); + + AliPMDddlinfoData *ddlinfo = 0; + if (entry) ddlinfo = (AliPMDddlinfoData*) entry->GetObject(); + + if (!ddlinfo) AliFatal("No ddl info data from database !"); + + return ddlinfo; +} +//____________________________________________________________________________ +AliPMDMappingData* AliPMDDDLRawData::GetMappingData() const +{ + // The run number will be centralized in AliCDBManager, + // you don't need to set it here! + AliCDBEntry *entry = AliCDBManager::Instance()->Get("PMD/Calib/Mapping"); + + if(!entry) AliFatal("Mapping object retrieval failed!"); + + AliPMDMappingData *mapda = 0; + if (entry) mapda = (AliPMDMappingData*) entry->GetObject(); + + if (!mapda) AliFatal("No mapping data from database !"); + + return mapda; +} //____________________________________________________________________________