X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONRawWriter.cxx;h=ba24638909ffa071e351789d98966b5d7eef621c;hb=73cfbcd5e63f1980fe1edd720deb9a3f68a878a1;hp=01476bd1cdf18e2a8a426a1d54bf14739ec7930d;hpb=a19e254322c6f837c1c4da1fd4f7a733952ec499;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/AliMUONRawWriter.cxx b/MUON/AliMUONRawWriter.cxx index 01476bd1cdf..ba24638909f 100644 --- a/MUON/AliMUONRawWriter.cxx +++ b/MUON/AliMUONRawWriter.cxx @@ -13,611 +13,708 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ -//////////////////////////////////// -// -// MUON Raw Data generator and reader in ALICE-MUON -// This class version 3 (further details could be found in Alice-note) -// -// Implemented non-constant buspatch numbers for tracking -// with correct DDL id (first guess) -// (Ch. Finck, dec 2005) -// -// Digits2Raw: -// Generates raw data for MUON tracker and finally for trigger -// Using real mapping (inverse) for tracker -// For trigger there is no mapping (mapping could be found in AliMUONTriggerCircuit) -// Ch. Finck july 04 -// -// Raw2Digits: -// Using real mapping for tracker -// Indranil Das (Adapted for runloader: Ch. Finck) july 05 -// -//////////////////////////////////// - -#include -#include - -#include - -#include "AliLoader.h" -#include "AliBitPacking.h" -#include "AliRawReader.h" -#include "AliLog.h" -#include "AliRun.h" +/* $Id$ */ + +//----------------------------------------------------------------------------- +/// \class AliMUONRawWriter +/// MUON Raw Data generaton in ALICE-MUON +/// Raw data structure could be found in Alice-note. +/// +/// Implemented non-constant buspatch numbers for tracking +/// with correct DDL id (first guess) +/// (Ch. Finck, dec 2005) +/// +/// Digits2Raw: +/// Generates raw data for MUON tracker and finally for trigger +/// Using real mapping (inverse) for tracker +/// For trigger there is no mapping (mapping could be found in AliMUONTriggerCircuit) +/// Ch. Finck, July 04 +/// Use memcpy instead of assignment elt by elt +/// Introducing variable DSP numbers, real manu numbers per buspatch for st12 +/// Implemented scaler event for Trigger +/// Ch. Finck, Jan. 06 +/// Using bus itr in DDL instead of simple incrementation +/// treat correctly the DDL & buspatch for station 3. +/// Using informations from AliMUONTriggerCrateStore for +/// empty slots and non-notified cards in trigger crates. +/// Ch. Finck, August 06. +/// Using AliMpDDLStore::GetBusPatchId. +/// +/// \author Ch. Finck, Feb. 07. +//----------------------------------------------------------------------------- + -#include "AliMUON.h" #include "AliMUONRawWriter.h" -#include "AliMUONDigit.h" +#include "AliMUONBlockHeader.h" +#include "AliMUONBusStruct.h" #include "AliMUONConstants.h" -#include "AliMUONData.h" - -#include "AliMUONSubEventTrigger.h" -#include "AliMUONDDLTracker.h" -#include "AliMUONDDLTrigger.h" - -#include "AliMUONLocalTrigger.h" +#include "AliMUONDarcHeader.h" +#include "AliMUONVDigit.h" +#include "AliMUONVDigitStore.h" +#include "AliMUONDspHeader.h" #include "AliMUONGlobalTrigger.h" - -#include "AliMUONGeometrySegmentation.h" -#include "AliMUONGeometryModule.h" -#include "AliMUONGeometryStore.h" -#include "AliMpSegFactory.h" +#include "AliMUONLocalStruct.h" +#include "AliMUONLocalTrigger.h" +#include "AliMUONLocalTriggerBoard.h" +#include "AliMUONRegionalTrigger.h" +#include "AliMUONRegHeader.h" + +#include "AliMUONVTriggerStore.h" +#include "AliCodeTimer.h" + +#include "AliMpCDB.h" +#include "AliMpDDLStore.h" +#include "AliMpDDL.h" +#include "AliMpRegionalTrigger.h" +#include "AliMpTriggerCrate.h" +#include "AliMpLocalBoard.h" +#include "AliMpDetElement.h" +#include "AliMpDEManager.h" +#include "AliMpExMap.h" +#include "AliMpConstants.h" #include "AliMpPlaneType.h" +#include "AliMpSegmentation.h" +#include "AliMpStationType.h" #include "AliMpVSegmentation.h" -#include "AliMpHelper.h" -#include "AliMpPad.h" +#include "AliRawReader.h" +#include "AliRawDataHeaderSim.h" +#include "AliBitPacking.h" +#include "AliDAQ.h" +#include "AliLog.h" -ClassImp(AliMUONRawWriter) // Class implementation in ROOT context -//__________________________________________________________________________ -AliMUONRawWriter::AliMUONRawWriter(AliLoader* loader, AliMUONData* data) - : TObject() -{ - // Standard Constructor - - // initialize loader's - fLoader = loader; - - // initialize segmentation factory - fSegFactory = new AliMpSegFactory(); - - // initialize container -// fMUONData = new AliMUONData(fLoader,"MUON","MUON"); - fMUONData = data; - - // initialize array - fSubEventArray = new TClonesArray("AliMUONSubEventTracker",1000); - - - // ddl pointer - fDDLTracker = new AliMUONDDLTracker(); - fDDLTrigger = new AliMUONDDLTrigger(); - - fBusPatchManager = new AliMpBusPatch(); - fBusPatchManager->ReadBusPatchFile(); +#include "TObjArray.h" +#include "TStopwatch.h" +#include -} +/// \cond CLASSIMP +ClassImp(AliMUONRawWriter) // Class implementation in ROOT context +/// \endcond //__________________________________________________________________________ AliMUONRawWriter::AliMUONRawWriter() : TObject(), - fMUONData(0), - fLoader(0), - fSegFactory(0), - fDDLTracker(0), - fDDLTrigger(0) + fBlockHeader(new AliMUONBlockHeader()), + fDspHeader(new AliMUONDspHeader()), + fDarcHeader(new AliMUONDarcHeader()), + fRegHeader(new AliMUONRegHeader()), + fLocalStruct(new AliMUONLocalStruct()), + fDDLStore(AliMpDDLStore::Instance()), + fScalerEvent(kFALSE), + fHeader(0x0), + fBufferSize((((43*AliMpConstants::ManuNofChannels() + 4)*5 + 10)*5 + 8)*2), + fBuffer(new Int_t [fBufferSize]) { - // Default Constructor - fFile[0] = fFile[1] = 0x0; - -} + /// Standard Constructor -//_______________________________________________________________________ -AliMUONRawWriter::AliMUONRawWriter (const AliMUONRawWriter& rhs) - : TObject(rhs) -{ -// Protected copy constructor + AliDebug(1,"Standard ctor"); - AliFatal("Not implemented."); -} - -//_______________________________________________________________________ -AliMUONRawWriter & -AliMUONRawWriter::operator=(const AliMUONRawWriter& rhs) -{ -// Protected assignement operator - - if (this == &rhs) return *this; + // setting data key to default value (only for writting) + fBlockHeader->SetDataKey(fBlockHeader->GetDefaultDataKey()); + fDspHeader->SetDataKey(fDspHeader->GetDefaultDataKey()); - AliFatal("Not implemented."); - - return *this; + // Load mapping + if ( ! fDDLStore ) { + if ( ! AliMpCDB::LoadDDLStore() ) { + AliFatal("Could not access mapping from OCDB !"); + } + fDDLStore = AliMpDDLStore::Instance(); + } } //__________________________________________________________________________ AliMUONRawWriter::~AliMUONRawWriter(void) { - if (fSegFactory) - fSegFactory->DeleteSegmentations(); - delete fSegFactory; + /// Destructor - if (fSubEventArray) - fSubEventArray->Delete(); //using delete cos allocating memory in copy ctor. + AliDebug(1,"dtor"); + + delete fBlockHeader; + delete fDspHeader; + delete fDarcHeader; + delete fRegHeader; + delete fLocalStruct; + delete[] fBuffer; +} - if (fDDLTracker) - delete fDDLTracker; - if (fDDLTrigger) - delete fDDLTrigger; +//____________________________________________________________________ +void AliMUONRawWriter::LocalWordPacking(UInt_t& word, UInt_t locId, UInt_t locDec, + UInt_t trigY, UInt_t posY, UInt_t posX, + UInt_t sdevX, UInt_t devX) +{ +/// pack local trigger word - fBusPatchManager->Delete(); + AliBitPacking::PackWord(locId,word,19,22); //card id number in crate + AliBitPacking::PackWord(locDec,word,15,18); + AliBitPacking::PackWord(trigY,word,14,14); + AliBitPacking::PackWord(posY,word,10,13); + AliBitPacking::PackWord(sdevX,word,9,9); + AliBitPacking::PackWord(devX,word,5,8); + AliBitPacking::PackWord(posX,word,0,4); - return; } + //____________________________________________________________________ -Int_t AliMUONRawWriter::Digits2Raw() +Int_t AliMUONRawWriter::Digits2Raw(const AliMUONVDigitStore* digitStore, + const AliMUONVTriggerStore* triggerStore) { - // convert digits of the current event to raw data + /// convert digits of the current event to raw data + AliCodeTimerAuto("",0) + Int_t idDDL; - Char_t name[20]; - - fLoader->LoadDigits("READ"); - - fMUONData->SetTreeAddress("D,GLT"); - // tracking chambers + + if ( digitStore ) + { + AliCodeTimerAuto("for Tracker",1) - for (Int_t ich = 0; ich < AliMUONConstants::NTrackingCh(); ich++) { - - // open files - idDDL = ich * 2 + 0x900; // official number for MUON - sprintf(name, "MUON_%d.ddl",idDDL); - fFile[0] = fopen(name,"w"); + AliMpExMap busPatchMap; - idDDL = (ich * 2) + 1 + 0x900; - sprintf(name, "MUON_%d.ddl",idDDL); - fFile[1] = fopen(name,"w"); + Int_t nDDLs = AliDAQ::NumberOfDdls("MUONTRK"); + + Int_t nofBusPatches(0); + + for (Int_t iDDL = 0; iDDL < nDDLs; ++iDDL ) + { + AliMpDDL* ddl = fDDLStore->GetDDL(iDDL); + nofBusPatches += ddl->GetNofBusPatches(); + } + + busPatchMap.SetSize(nofBusPatches); + + Digits2BusPatchMap(*digitStore,busPatchMap); - WriteTrackerDDL(ich); - - // reset and close - fclose(fFile[0]); - fclose(fFile[1]); - fMUONData->ResetDigits(); + for (Int_t iDDL = 0; iDDL < nDDLs; ++iDDL ) + { + WriteTrackerDDL(busPatchMap,iDDL); + } + AliDebug(1,"Tracker written"); } - // trigger chambers - - // open files - idDDL = 0xA00;// official number for MUTR - sprintf(name, "MUTR_%d.ddl",idDDL); - fFile[0] = fopen(name,"w"); - - idDDL = 0xA00 + 1; - sprintf(name, "MUTR_%d.ddl",idDDL); - fFile[1] = fopen(name,"w"); + if ( triggerStore ) + { + AliCodeTimerAuto("for Trigger",1) - WriteTriggerDDL(); - - // reset and close - fclose(fFile[0]); - fclose(fFile[1]); - fMUONData->ResetTrigger(); + // trigger chambers + + AliFstream* file[2]; + + // open files + idDDL = 0;// MUTR + file[0] = new AliFstream(AliDAQ::DdlFileName("MUONTRG",idDDL)); + + idDDL = 1;// MUTR + file[1] = new AliFstream(AliDAQ::DdlFileName("MUONTRG",idDDL)); + + WriteTriggerDDL(*triggerStore,file); + + // reset and close + delete file[0]; + delete file[1]; + + AliDebug(1,"Trigger written"); + } - fLoader->UnloadDigits(); - return kTRUE; } -//____________________________________________________________________ -Int_t AliMUONRawWriter::WriteTrackerDDL(Int_t iCh) -{ - // writing DDL for tracker - // used inverse mapping - - // resets - TClonesArray* muonDigits = 0; - fSubEventArray->Clear(); - - // - TArrayI nbInBus; - - nbInBus.Set(5000); - - nbInBus.Reset(); - - // DDL header - AliRawDataHeader header = fDDLTracker->GetHeader(); - Int_t headerSize = fDDLTracker->GetHeaderSize(); +//______________________________________________________________________________ +void +AliMUONRawWriter::Digits2BusPatchMap(const AliMUONVDigitStore& digitStore, + AliMpExMap& busPatchMap) +{ + /// Create bus patch structures corresponding to digits in the store + + AliCodeTimerAuto("",0) + + static const Int_t kMAXADC = (1<<12)-1; // We code the charge on a 12 bits ADC. + // DDL event one per half chamber - AliMUONSubEventTracker* subEvent; - - // data format + + // raw data Char_t parity = 0x4; UShort_t manuId = 0; UChar_t channelId = 0; UShort_t charge = 0; Int_t busPatchId = 0; - + Int_t currentBusPatchId = -1; UInt_t word; - Int_t nEntries = 0; - Int_t* buffer = 0; - Int_t index; - Int_t indexDsp; - Int_t indexBlk; - Int_t padX; - Int_t padY; - Int_t cathode = 0; - Int_t detElemId; - Int_t nDigits; - - const AliMUONDigit* digit; - - AliDebug(3, Form("WriteDDL chamber %d\n", iCh+1)); - - // getting digits - fMUONData->ResetDigits(); - fMUONData->GetDigits(); - muonDigits = fMUONData->Digits(iCh); - - nDigits = muonDigits->GetEntriesFast(); - AliDebug(3,Form("ndigits = %d\n",nDigits)); - - // loop over digit - for (Int_t idig = 0; idig < nDigits; idig++) { - - digit = (AliMUONDigit*) muonDigits->UncheckedAt(idig); - - padX = digit->PadX(); - padY = digit->PadY(); - charge = digit->Signal(); - charge &= 0xFFF; - cathode = digit->Cathode(); - detElemId = digit->DetElemId(); - + + AliMUONBusStruct* busStruct(0x0); + + TIter next(digitStore.CreateTrackerIterator()); + AliMUONVDigit* digit; + + while ( ( digit = static_cast(next()) ) ) + { + charge = digit->ADC(); + if ( charge > kMAXADC ) + { + // This is most probably an error in the digitizer (which should insure + // the adc is below kMAXADC), so make it a (non-fatal) error indeed. + AliError(Form("adc value %d above 0x%x for DE %d . Setting to 0x%x. Digit is:", + charge,kMAXADC,digit->DetElemId(),kMAXADC)); + StdoutToAliError(digit->Print()); + charge = kMAXADC; + } + // inverse mapping - Int_t error = GetInvMapping(digit, busPatchId, manuId, channelId); - if (error) continue; - - AliDebug(3,Form("input IdDE %d busPatchId %d PadX %d PadY %d iCath %d \n", - detElemId, busPatchId, padX, padY, cathode)); - - AliDebug(3,Form("busPatchId %d, manuId %d channelId %d\n", busPatchId, manuId, channelId )); + busPatchId = GetBusPatch(*digit); + if (busPatchId<0) continue; + + if ( digit->ManuId() > 0x7FF || + digit->ManuChannel() > 0x3F ) + { + StdoutToAliError(digit->Print();); + AliFatal("ManuId,ManuChannel are invalid for the digit above."); + } + + manuId = ( digit->ManuId() & 0x7FF ); // 11 bits + channelId = ( digit->ManuChannel() & 0x3F ); // 6 bits + //packing word - AliBitPacking::PackWord((UInt_t)parity,word,29,31); + word = 0; AliBitPacking::PackWord((UInt_t)manuId,word,18,28); AliBitPacking::PackWord((UInt_t)channelId,word,12,17); AliBitPacking::PackWord((UInt_t)charge,word,0,11); + + // parity word + parity = word & 0x1; + for (Int_t i = 1; i <= 30; ++i) + { + parity ^= ((word >> i) & 0x1); + } + AliBitPacking::PackWord((UInt_t)parity,word,31,31); - // set sub Event - subEvent = new AliMUONSubEventTracker(); - subEvent->AddData(word); - subEvent->SetBusPatchId(busPatchId); - - // storing the number of identical buspatches - nbInBus[busPatchId]++; - AddData(subEvent); - - delete subEvent; - } - - // sorting by buspatch - fSubEventArray->Sort(); - - // gather datas from same bus patch - nEntries = fSubEventArray->GetEntriesFast(); - - for (Int_t i = 0; i < nEntries; i++) { - AliMUONSubEventTracker* temp = (AliMUONSubEventTracker*)fSubEventArray->At(i); - busPatchId = temp->GetBusPatchId(); - - // add bus patch header, length and total length managed by subevent class - temp->SetTriggerWord(0xdeadbeef); - for (Int_t j = 0; j < nbInBus[busPatchId]-1; j++) { - AliMUONSubEventTracker* temp1 = (AliMUONSubEventTracker*)fSubEventArray->At(++i); - temp->AddData(temp1->GetData(0)); - fSubEventArray->RemoveAt(i) ; + if ( currentBusPatchId != busPatchId ) + { + busStruct = + static_cast(busPatchMap.GetValue(busPatchId)); + currentBusPatchId = busPatchId; } - } - fSubEventArray->Compress(); - - if (AliLog::GetGlobalDebugLevel() == 3) { - nEntries = fSubEventArray->GetEntriesFast(); - for (Int_t i = 0; i < nEntries; i++) { - AliMUONSubEventTracker* temp = (AliMUONSubEventTracker*)fSubEventArray->At(i); - printf("busPatchid back %d\n",temp->GetBusPatchId()); - for (Int_t j = 0; j < temp->GetLength(); j++) { - printf("manuId back %d, ",temp->GetManuId(j)); - printf("channelId back %d, ",temp->GetChannelId(j)); - printf("charge back %d\n",temp->GetCharge(j)); - } + + if (!busStruct) + { + busStruct = new AliMUONBusStruct; + busStruct->SetDataKey(busStruct->GetDefaultDataKey()); + busStruct->SetBusPatchId(busPatchId); + busStruct->SetLength(0); + busPatchMap.Add(busPatchId,busStruct); } - printf("\n"); + + // set sub Event + busStruct->AddData(word); } +} + +//______________________________________________________________________________ +void +AliMUONRawWriter::WriteTrackerDDL(AliMpExMap& busPatchMap, Int_t iDDL) +{ + /// Write DDL file for one tracker DDL - // getting info for the number of buspatches - Int_t iBusPatch; - Int_t length; - Int_t iBusPerDSP[5];//number of bus patches per DSP - Int_t iDspMax; //number max of DSP per block + // buffer size (max'ed out) + // (((43 manus max per bus patch *64 channels + 4 bus patch words) * 5 bus patch + // + 10 dsp words)*5 dsps + 8 block words)*2 blocks - Int_t iFile = 0; - fBusPatchManager->GetDspInfo(iCh, iDspMax, iBusPerDSP); - - TArrayI* vec = fBusPatchManager->GetBusfromDE((iCh+1)*100); - - Int_t iBus0AtCh = vec->At(0); //get first bus patch id for a given ich - - AliDebug(3,Form("iBus0AtCh %d", iBus0AtCh)); - - iBusPatch = iBus0AtCh - 1; // starting point for each chamber - - // nEntries = fSubEventArray->GetEntriesFast(); - AliMUONSubEventTracker* temp = 0x0; + AliCodeTimerAuto("",0) - // open DDL file, on per 1/2 chamber - for (Int_t iDDL = 0; iDDL < 2; iDDL++) { + if (fHeader == 0x0) { + AliError("Raw data header must be set"); + return; + } + memset(fBuffer,0,fBufferSize*sizeof(Int_t)); + + AliMpDDL* ddl = fDDLStore->GetDDL(iDDL); + Int_t iDspMax = ddl->GetMaxDsp(); + Int_t iBusPerDSP[5]; //number of bus patches per DSP + ddl->GetBusPerDsp(iBusPerDSP); + Int_t busIter = 0; + + Int_t totalDDLLength = 0; + + Int_t index = 0; + + // two blocks A and B per DDL + for (Int_t iBlock = 0; iBlock < 2; ++iBlock) + { + // block header + Int_t length = fBlockHeader->GetHeaderLength(); + memcpy(&fBuffer[index],fBlockHeader->GetHeader(),length*4); + Int_t indexBlk = index; + index += length; - - // filling buffer - buffer = new Int_t [(2048+24)*50]; // 24 words in average for one buspatch and 2048 manu info at most - - indexBlk = 0; - indexDsp = 0; - index = 0; - - // two blocks A and B per DDL - for (Int_t iBlock = 0; iBlock < 2; iBlock++) { - - // block header - length = fDDLTracker->GetBlkHeaderLength(); - memcpy(&buffer[index],fDDLTracker->GetBlkHeader(),length*4); - indexBlk = index; - index += length; - - // 5 DSP's max per block - for (Int_t iDsp = 0; iDsp < iDspMax; iDsp++) { - - // DSP header - length = fDDLTracker->GetDspHeaderLength(); - memcpy(&buffer[index],fDDLTracker->GetDspHeader(),length*4); - indexDsp = index; - index += length; - - // 5 buspatches max per DSP - for (Int_t i = 0; i < iBusPerDSP[iDsp]; i++) { - - iBusPatch ++; - if ((fBusPatchManager->GetDDLfromBus(iBusPatch) % 2) == 1) // comparing to DDL file - iFile = 0; - else - iFile = 1; - - AliDebug(3,Form("iCh %d iDDL %d iBlock %d iDsp %d busPatchId %d", iCh, iDDL, iBlock, iDsp, iBusPatch)); - - nEntries = fSubEventArray->GetEntriesFast(); - - for (Int_t iEntries = 0; iEntries < nEntries; iEntries++) { // method "bourrique"... - temp = (AliMUONSubEventTracker*)fSubEventArray->At(iEntries); - busPatchId = temp->GetBusPatchId(); - if (busPatchId == iBusPatch) break; - busPatchId = -1; - AliDebug(3,Form("busPatchId %d", temp->GetBusPatchId())); - } - - // check if buspatchid has digit - if (busPatchId != -1) { - // add bus patch structure - length = temp->GetHeaderLength(); - memcpy(&buffer[index],temp->GetAddress(),length*4); - index += length; - for (Int_t j = 0; j < temp->GetLength(); j++) { - buffer[index++] = temp->GetData(j); - AliDebug(3,Form("busPatchId %d, manuId %d channelId %d\n", temp->GetBusPatchId(), - temp->GetManuId(j), temp->GetChannelId(j) )); - } - // fSubEventArray->RemoveAt(iEntries); - // fSubEventArray->Compress(); - } else { - // writting anyhow buspatch structure (empty ones) - buffer[index++] = 4; // total length - buffer[index++] = 0; // raw data length - buffer[index++] = iBusPatch; // bus patch - buffer[index++] = 0xdeadbeef; // trigger word - } - } // bus patch - buffer[indexDsp] = index - indexDsp; // dsp length - buffer[indexDsp+1] = index - indexDsp - fDDLTracker->GetDspHeaderLength(); - if ((index - indexDsp) % 2 == 0) - buffer[indexDsp+7] = 0; - else - buffer[indexDsp+7] = 1; - } // dsp - buffer[indexBlk] = index - indexBlk; // block length - buffer[indexBlk+1] = index - indexBlk - fDDLTracker->GetBlkHeaderLength(); - } + // 5 DSP's max per block + for (Int_t iDsp = 0; iDsp < iDspMax; ++iDsp) + { + // DSP header + Int_t dspHeaderLength = fDspHeader->GetHeaderLength(); + memcpy(&fBuffer[index],fDspHeader->GetHeader(),dspHeaderLength*4); + Int_t indexDsp = index; + index += dspHeaderLength; + + // 5 buspatches max per DSP + for (Int_t i = 0; i < iBusPerDSP[iDsp]; ++i) + { + Int_t iBusPatch = ddl->GetBusPatchId(busIter++); + + // iteration over bus patch in DDL + if (iBusPatch == -1) + { + AliWarning(Form("Error in bus itr in DDL %d\n", iDDL)); + continue; + } + + AliMUONBusStruct* busStructPtr = static_cast(busPatchMap.GetValue(iBusPatch)); + + // check if buspatchid has digit + if (busStructPtr) + { + // add bus patch structure header + Int_t busHeaderLength = busStructPtr->GetHeaderLength(); + memcpy(&fBuffer[index],busStructPtr->GetHeader(),busHeaderLength*4); + index += busHeaderLength; + + // add bus patch data + Int_t busLength = busStructPtr->GetLength(); + memcpy(&fBuffer[index],busStructPtr->GetData(),busLength*4); + index += busLength; + } + else + { + // writting anyhow buspatch structure (empty ones) + fBuffer[index++] = busStructPtr->GetDefaultDataKey(); // fill it also for empty data size + fBuffer[index++] = busStructPtr->GetHeaderLength(); // header length + fBuffer[index++] = 0; // raw data length + fBuffer[index++] = iBusPatch; // bus patch + } + } // bus patch + + // check if totalLength even + // set padding word in case + // Add one word 0xBEEFFACE at the end of DSP structure + Int_t totalDspLength = index - indexDsp; + if ((totalDspLength % 2) == 1) + { + fBuffer[indexDsp + fDspHeader->GetHeaderLength() - 2] = 1; + fBuffer[index++] = fDspHeader->GetDefaultPaddingWord(); + totalDspLength++; + } + + Int_t dspLength = totalDspLength - fDspHeader->GetHeaderLength(); + + fBuffer[indexDsp+1] = totalDspLength; // dsp total length + fBuffer[indexDsp+2] = dspLength; // data length + + } // dsp - //writting onto disk - // write DDL 1 & 2 - header.fSize = (index + headerSize) * 4;// total length in bytes - fwrite((char*)(&header),headerSize*4,1,fFile[iFile]); - fwrite(buffer,sizeof(int),index,fFile[iFile]); - - delete[] buffer; - } + Int_t totalBlkLength = index - indexBlk; + Int_t blkLength = totalBlkLength - fBlockHeader->GetHeaderLength(); + totalDDLLength += totalBlkLength; + + fBuffer[indexBlk+1] = totalBlkLength; // total block length + fBuffer[indexBlk+2] = blkLength; + + } // block + + // add twice the end of CRT structure data key + // hope it's good placed (ChF) + fBuffer[index++] = fBlockHeader->GetDdlDataKey(); + fBuffer[index++] = fBlockHeader->GetDdlDataKey(); + totalDDLLength += 2; + + // writting onto disk + // total length in bytes + // DDL header - return kTRUE; + Int_t headerSize = sizeof(AliRawDataHeader)/4; + + fHeader->fSize = (totalDDLLength + headerSize) * 4; + + AliFstream* file = new AliFstream(AliDAQ::DdlFileName("MUONTRK",iDDL)); + + file->WriteBuffer((char*)fHeader,headerSize*4); + file->WriteBuffer((char*)fBuffer,sizeof(int)*index); + delete file; } -//____________________________________________________________________ -Int_t AliMUONRawWriter::WriteTriggerDDL() -{ - // DDL event one per half chamber - AliMUONSubEventTrigger* subEvent = 0x0; +//______________________________________________________________________________ +Int_t AliMUONRawWriter::GetBusPatch(const AliMUONVDigit& digit) const +{ + /// Determine the BusPatch this digit belongs to. + return fDDLStore->GetBusPatchId(digit.DetElemId(),digit.ManuId()); +} - // stored local id number - TArrayI isFired(256); - isFired.Reset(); +//______________________________________________________________________________ +Int_t AliMUONRawWriter::WriteTriggerDDL(const AliMUONVTriggerStore& triggerStore, AliFstream* file[2]) +{ + /// Write trigger DDL + + AliCodeTimerAuto("",0) + if (fHeader == 0x0) { + AliError("Raw data header must be set"); + return 0; + } - // DDL header - AliRawDataHeader header = fDDLTrigger->GetHeader(); - Int_t headerSize = fDDLTrigger->GetHeaderSize(); - Int_t length; - TClonesArray* localTrigger; - TClonesArray* globalTrigger; - AliMUONGlobalTrigger* gloTrg; - AliMUONLocalTrigger* locTrg = 0x0; + // DDL event one per half chamber - // getting information from trigger - fMUONData->GetTriggerD(); + // DDL header size + Int_t headerSize = sizeof(AliRawDataHeader)/4; // global trigger for trigger pattern - globalTrigger = fMUONData->GlobalTrigger(); - gloTrg = (AliMUONGlobalTrigger*)globalTrigger->UncheckedAt(0); - Int_t gloTrigPat = GetGlobalTriggerPattern(gloTrg); - - // local trigger - localTrigger = fMUONData->LocalTrigger(); + AliMUONGlobalTrigger* gloTrg = triggerStore.Global(); + if (!gloTrg) + { + return 0; + } + + Int_t gloTrigResp = gloTrg->GetGlobalResponse(); + UInt_t *gloTrigInput = gloTrg->GetGlobalInput(); UInt_t word; Int_t* buffer = 0; Int_t index; - Int_t iEntries = 0; - Int_t iLocCard, locCard; - Char_t locDec, trigY, posY, posX,regOut; - Int_t devX; - Int_t version = 1; // software version - Int_t eventType = 1; // trigger type: 1 for physics ? - Int_t serialNb = 0xF; // serial nb of card: all bits on for the moment - Int_t globalFlag = 1; // set to 2 if global info present in DDL else set to 1 - - Int_t nEntries = (Int_t) (localTrigger->GetEntries());// 234 local cards - // stored the local card id that's fired - for (Int_t i = 0; i < nEntries; i++) { - locTrg = (AliMUONLocalTrigger*)localTrigger->At(i); - isFired[locTrg->LoCircuit()] = 1; // storing local boards with informations + Int_t locCard; + UChar_t locDec, trigY, posY, posX, regOut; + UInt_t regInpLpt; + UInt_t regInpHpt; + + UInt_t devX; + UChar_t sdevX; + UInt_t version = 1; // software version + UInt_t eventPhys = 1; // trigger type: 1 for physics, 0 for software + UInt_t serialNb = 0xF; // serial nb of card: all bits on for the moment + Int_t globalFlag = 0; // set to 1 if global info present in DDL else set to 0 + + // size of headers + static const Int_t kDarcHeaderLength = fDarcHeader->GetDarcHeaderLength(); + static const Int_t kGlobalHeaderLength = fDarcHeader->GetGlobalHeaderLength(); + static const Int_t kDarcScalerLength = fDarcHeader->GetDarcScalerLength(); + static const Int_t kGlobalScalerLength = fDarcHeader->GetGlobalScalerLength(); + static const Int_t kRegHeaderLength = fRegHeader->GetHeaderLength(); + static const Int_t kRegScalerLength = fRegHeader->GetScalerLength(); + static const Int_t kLocHeaderLength = fLocalStruct->GetLength(); + static const Int_t kLocScalerLength = fLocalStruct->GetScalerLength(); + + // [16(local)*6 words + 6 words]*8(reg) + 8 words = 824 + static const Int_t kBufferSize = (16 * (kLocHeaderLength+1) + (kRegHeaderLength+1))* 8 + + kDarcHeaderLength + kGlobalHeaderLength + 2; + + // [16(local)*51 words + 16 words]*8(reg) + 8 + 10 + 8 words scaler event 6682 words + static const Int_t kScalerBufferSize = (16 * (kLocHeaderLength + kLocScalerLength +1) + + (kRegHeaderLength + kRegScalerLength +1))* 8 + + (kDarcHeaderLength + kDarcScalerLength + + kGlobalHeaderLength + kGlobalScalerLength + 2); + if(fScalerEvent) { + eventPhys = 0; //set to generate scaler events + fHeader->fWord2 |= (0x1 << 14); // set L1SwC bit on } + if(fScalerEvent) + buffer = new Int_t [kScalerBufferSize]; + else + buffer = new Int_t [kBufferSize]; - if (!nEntries) - AliError("No Trigger information available"); - - buffer = new Int_t [672]; // [16(local)*5 words + 3 words]*8(reg) + 8 words = 672 + // reset crate // open DDL file, on per 1/2 chamber - for (Int_t iDDL = 0; iDDL < 2; iDDL++) { - + for ( Int_t iDDL = 0; iDDL < 2; ++iDDL ) + { index = 0; - // DDL enhanced header - word = 0; - AliBitPacking::PackWord((UInt_t)iDDL+1,word,28,31); //see AliMUONDDLTrigger.h for details - AliBitPacking::PackWord((UInt_t)serialNb,word,24,27); - AliBitPacking::PackWord((UInt_t)version,word,16,23); - AliBitPacking::PackWord((UInt_t)eventType,word,12,15); - if (iDDL == 0) // suppose global info in DDL one - globalFlag = 2; - else globalFlag = 1; - AliBitPacking::PackWord((UInt_t)globalFlag,word,8,11); - fDDLTrigger->SetDDLWord(word); - - if (iDDL == 0) - fDDLTrigger->SetGlobalOutput(gloTrigPat);// no global input for the moment.... else - fDDLTrigger->SetGlobalOutput(0); - length = fDDLTrigger->GetHeaderLength(); - memcpy(&buffer[index],fDDLTrigger->GetEnhancedHeader(),length*4); - index += length; + globalFlag = 0; - // 8 regional cards per DDL - for (Int_t iReg = 0; iReg < 8; iReg++) { + word = 0; + // set darc status word + // see AliMUONDarcHeader.h for details + AliBitPacking::PackWord((UInt_t)eventPhys,word,30,30); + AliBitPacking::PackWord((UInt_t)serialNb,word,20,23); + AliBitPacking::PackWord((UInt_t)globalFlag,word,10,10); + AliBitPacking::PackWord((UInt_t)version,word,12,19); + fDarcHeader->SetWord(word); + + memcpy(&buffer[index], fDarcHeader->GetHeader(), (kDarcHeaderLength)*4); + index += kDarcHeaderLength; + + // no global input for the moment.... + if (iDDL == 0) { + fDarcHeader->SetGlobalOutput(gloTrigResp); + for (Int_t ii = 0; ii < 4; ii++) { + fDarcHeader->SetGlobalInput(gloTrigInput[ii],ii); + } + } else { + fDarcHeader->SetGlobalOutput(0); + } + + if (fScalerEvent) { + // 6 DARC scaler words + memcpy(&buffer[index], fDarcHeader->GetDarcScalers(),kDarcScalerLength*4); + index += kDarcScalerLength; + } + // end of darc word + buffer[index++] = fDarcHeader->GetEndOfDarc(); + + // 4 words of global board input + Global board output + memcpy(&buffer[index], fDarcHeader->GetGlobalInput(), (kGlobalHeaderLength)*4); + index += kGlobalHeaderLength; - subEvent = new AliMUONSubEventTrigger(); + if (fScalerEvent) { + // 10 Global scaler words + memcpy(&buffer[index], fDarcHeader->GetGlobalScalers(), kGlobalScalerLength*4); + index += kGlobalScalerLength; + } + + // end of global word + buffer[index++] = fDarcHeader->GetEndOfGlobal(); + const AliMpRegionalTrigger* reg = AliMpDDLStore::Instance()->GetRegionalTrigger(); + Int_t nCrate = reg->GetNofTriggerCrates()/2; + // 8 regional cards per DDL + for (Int_t iReg = 0; iReg < nCrate; ++iReg) { + + // crate info + AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->GetTriggerCrate(iDDL, iReg); + + if (!crate) { + AliError(Form("Missing crate number %d in DDL %d\n", iReg, iDDL)); + continue; + } + + // regional info tree, make sure that no reg card missing + AliMUONRegionalTrigger* regTrg = triggerStore.FindRegional(crate->GetId()); + if (!regTrg) { + AliError(Form("Missing regional board %d in trigger Store\n", crate->GetId())); + continue; + } + // Regional card header word = 0; - regOut = 0; - AliBitPacking::PackWord((UInt_t)serialNb,word,24,28); //see AliMUONSubEventTrigger.h for details - AliBitPacking::PackWord((UInt_t)version,word,16,23); - AliBitPacking::PackWord((UInt_t)iReg,word,12,15); - AliBitPacking::PackWord((UInt_t)regOut,word,0,7); // whenever regional output will be implemented - subEvent->SetRegWord(word); - memcpy(&buffer[index++],subEvent->GetAddress(),4); + // set darc status word + fRegHeader->SetDarcWord(word); - buffer[index++] = 0;// 2 words of regional input - buffer[index++] = 0; + regOut = regTrg->GetOutput(); + regInpLpt = regTrg->GetLocalOutput(0); + regInpHpt = regTrg->GetLocalOutput(1); - // 16 local card per regional board - for (Int_t iLoc = 0; iLoc < 16; iLoc++) { - - iLocCard = iLoc + iReg*16 + iDDL*128; - - if (isFired[iLocCard]) { - locTrg = (AliMUONLocalTrigger*)localTrigger->At(iEntries); - locCard = locTrg->LoCircuit(); - locDec = locTrg->GetLoDecision(); - trigY = 0; - posY = locTrg->LoStripY(); - posX = locTrg->LoStripX(); - devX = locTrg->LoDev(); - AliDebug(4,Form("loctrg %d, posX %d, posY %d, devX %d\n", - locTrg->LoCircuit(),locTrg->LoStripX(),locTrg->LoStripY(),locTrg->LoDev())); - } else { //no trigger (see PRR chpt 3.4) - locCard = -1; - locDec = 0; - trigY = 1; - posY = 15; - posX = 0; - devX = 0x8; - } + // fill darc word, not darc status for the moment (empty) + //see AliMUONRegHeader.h for details + AliBitPacking::PackWord((UInt_t)eventPhys,word,31,31); + AliBitPacking::PackWord((UInt_t)serialNb,word,20,25); + AliBitPacking::PackWord((UInt_t)version,word,8,15); + AliBitPacking::PackWord((UInt_t)crate->GetId(),word,16,19); + AliBitPacking::PackWord((UInt_t)regOut,word,0,7); + fRegHeader->SetWord(word); - //packing word - word = 0; - AliBitPacking::PackWord((UInt_t)(iLocCard % 16),word,19,22); //card id number in crate - AliBitPacking::PackWord((UInt_t)locDec,word,15,18); - AliBitPacking::PackWord((UInt_t)trigY,word,14,14); - AliBitPacking::PackWord((UInt_t)posY,word,10,13); - AliBitPacking::PackWord((UInt_t)devX,word,5,9); - AliBitPacking::PackWord((UInt_t)posX,word,0,4); - - if (locCard == iLocCard) { - // add local cards structure - buffer[index++] = (locTrg->GetX1Pattern() | (locTrg->GetX2Pattern() << 16)); - buffer[index++] = (locTrg->GetX3Pattern() | (locTrg->GetX4Pattern() << 16)); - buffer[index++] = (locTrg->GetY1Pattern() | (locTrg->GetY2Pattern() << 16)); - buffer[index++] = (locTrg->GetY3Pattern() | (locTrg->GetY4Pattern() << 16)); - buffer[index++] = (Int_t)word; // data word - if (iEntries < nEntries-1) - iEntries++; - } else { - buffer[index++] = 0; // 4 words for x1, x2, y1, y2 - buffer[index++] = 0; - buffer[index++] = 0; - buffer[index++] = 0; - buffer[index++] = (Int_t)word; // data word - } - } // local card + // fill header later, need local response + Int_t indexReg = index; + index += kRegHeaderLength; - delete subEvent; + // 11 regional scaler word + if (fScalerEvent) { + memcpy(&buffer[index], fRegHeader->GetScalers(), kRegScalerLength*4); + index += kRegScalerLength; + } + // end of regional word + buffer[index++] = fRegHeader->GetEndOfReg(); + + // 16 local card per regional board + // UShort_t localMask = 0x0; + + Int_t nLocalBoard = AliMpConstants::LocalBoardNofChannels(); + + for (Int_t iLoc = 0; iLoc < nLocalBoard; iLoc++) { + + // slot zero for Regional card + Int_t localBoardId = crate->GetLocalBoardId(iLoc); + + if (localBoardId) { // if not empty slot + AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(localBoardId); + + if (localBoard->IsNotified()) {// if notified board + AliMUONLocalTrigger* locTrg = triggerStore.FindLocal(localBoardId); + + locCard = locTrg->LoCircuit(); + locDec = locTrg->GetLoDecision(); + trigY = locTrg->LoTrigY(); + posY = locTrg->LoStripY(); + posX = locTrg->LoStripX(); + devX = locTrg->LoDev(); + sdevX = locTrg->LoSdev(); + + AliDebug(4,Form("loctrg %d, posX %d, posY %d, devX %d\n", + locTrg->LoCircuit(),locTrg->LoStripX(),locTrg->LoStripY(),locTrg->LoDev())); + //packing word + word = 0; + LocalWordPacking(word, (UInt_t)iLoc, (UInt_t)locDec, (UInt_t)trigY, (UInt_t)posY, + (UInt_t)posX, (UInt_t)sdevX, (UInt_t)devX); + + buffer[index++] = (locTrg->GetX1Pattern() | (locTrg->GetX2Pattern() << 16)); + buffer[index++] = (locTrg->GetX3Pattern() | (locTrg->GetX4Pattern() << 16)); + buffer[index++] = (locTrg->GetY1Pattern() | (locTrg->GetY2Pattern() << 16)); + buffer[index++] = (locTrg->GetY3Pattern() | (locTrg->GetY4Pattern() << 16)); + buffer[index++] = (Int_t)word; // data word + + + } + // fill copy card X-Y inputs from the notified cards + if (localBoard->GetInputXfrom() && localBoard->GetInputYfrom()) + { + // not triggered + locDec = 0; trigY = 1; posY = 15; + posX = 0; devX = 0; sdevX = 1; + LocalWordPacking(word, (UInt_t)iLoc, (UInt_t)locDec, (UInt_t)trigY, (UInt_t)posY, + (UInt_t)posX, (UInt_t)sdevX, (UInt_t)devX); + + Int_t localFromId = localBoard->GetInputXfrom(); + AliMUONLocalTrigger* locTrgfrom = triggerStore.FindLocal(localFromId); + + buffer[index++] = 0; // copy only X3-4 & Y1-4 + buffer[index++] = (locTrgfrom->GetX3Pattern() | (locTrgfrom->GetX4Pattern() << 16)); + buffer[index++] = (locTrgfrom->GetY1Pattern() | (locTrgfrom->GetY2Pattern() << 16)); + buffer[index++] = (locTrgfrom->GetY3Pattern() | (locTrgfrom->GetY4Pattern() << 16)); + buffer[index++] = word; + } + + } else { + // fill with 10CDEAD word for empty slots + for (Int_t i = 0; i < fLocalStruct->GetLength(); i++) + buffer[index++] = fLocalStruct->GetDisableWord(); + }// condition localBoard + + // 45 regional scaler word + if (fScalerEvent) { + memcpy(&buffer[index], fLocalStruct->GetScalers(), kLocScalerLength*4); + index += kLocScalerLength; + } + + // end of local structure words + buffer[index++] = fLocalStruct->GetEndOfLocal(); + + } // local card + // fill regional header with local output + fRegHeader->SetInput(regInpLpt, 0); + fRegHeader->SetInput(regInpHpt, 1); + memcpy(&buffer[indexReg],fRegHeader->GetHeader(),kRegHeaderLength*4); + } // Regional card - buffer[index++] = fDDLTrigger->GetEoD(); // End of DDL word - buffer[index++] = fDDLTrigger->GetEoD(); // End of DDL word for 64 bits transfer purpose // writting onto disk - // write DDL 1 - header.fSize = (index + headerSize) * 4;// total length in bytes - fwrite((char*)(&header),headerSize*4,1,fFile[iDDL]); - fwrite(buffer,sizeof(int),index,fFile[iDDL]); + // write DDL's + fHeader->fSize = (index + headerSize) * 4;// total length in bytes + file[iDDL]->WriteBuffer((char*)fHeader,headerSize*4); + file[iDDL]->WriteBuffer((char*)buffer,sizeof(int)*index); } delete[] buffer; @@ -626,109 +723,16 @@ Int_t AliMUONRawWriter::WriteTriggerDDL() } //____________________________________________________________________ -Int_t AliMUONRawWriter::GetInvMapping(const AliMUONDigit* digit, - Int_t &busPatchId, UShort_t &manuId, UChar_t &channelId) +void AliMUONRawWriter::SetScalersNumbers() { + /// set numbers for scaler events for trigger headers + /// since this is provided by the experiment + /// put dummy numbers to check the monitoring - // Inverse mapping for tracker - - // information from digits - Int_t iCath = digit->Cathode(); - Int_t idDE = digit->DetElemId(); - Int_t padX = digit->PadX(); - Int_t padY = digit->PadY(); - - if (idDE >= 500) { // Since in AliMpSlat pads begin at (0,0) - padX--; // while in AliMUONSt345Seg. they begin at (1,1) - padY--; - } - - // segmentation - AliMpPlaneType plane; - AliMpPlaneType plane1 = kBendingPlane; - AliMpPlaneType plane2 = kNonBendingPlane; - - if (idDE < 500) { // should use GetDirection somehow (ChF) - if ( ((idDE % 100) % 2) != 0 ) { - plane1 = kNonBendingPlane; - plane2 = kBendingPlane; - } - } - // station 345 bending == cath0 for the moment - plane = (iCath == 0) ? plane1 : plane2; - - //AliMpVSegmentation* seg = AliMUONSegmentationManager::Segmentation(idDE, plane); - AliMpVSegmentation* seg = fSegFactory->CreateMpSegmentation(idDE, iCath); - AliMpPad pad = seg->PadByIndices(AliMpIntPair(padX,padY),kTRUE); - - if(!pad.IsValid()) { - AliWarning(Form("No elec. for idDE: %d, padx: %d pady %d, charge: %d\n", - idDE, digit->PadX(), digit->PadY(), digit->Signal())); - return kTRUE; - } - - // Getting Manu id - manuId = pad.GetLocation().GetFirst(); - manuId &= 0x7FF; // 11 bits - - // Getting channel id - channelId = pad.GetLocation().GetSecond(); - channelId &= 0x3F; // 6 bits - - // Getting buspatch id - TArrayI* vec = fBusPatchManager->GetBusfromDE(idDE); - Int_t pos; - - if (idDE < 500) { // station 1 & 2 - // set 32 manus for one bus patch ? (ChF) - pos = manuId/32; - } else { - // offset of 100 in manuId for following bus patch - pos = manuId/100; - } - - // if (pos >(int_t) vec.size()) -// AliWarning("pos greater than size\n"); - busPatchId = vec->At(pos); - - if (plane == kNonBendingPlane) // for Non-Bending manuid+= 1000; - manuId += 1000; // tmp solution til one finds something better (ChF) - - AliDebug(3,Form("idDE: %d, busPatchId %d, manuId: %d, channelId:%d\n", - idDE, busPatchId, manuId, channelId)); - - AliDebug(3,Form("idDE: %d, busPatchId %d, manuId: %d, channelId: %d, padx: %d pady: %d, charge: %d\n", - idDE, busPatchId, manuId, channelId, digit->PadX(), digit->PadY(), digit->Signal())); - - return kFALSE; // no error -} - -//____________________________________________________________________ -Int_t AliMUONRawWriter::GetGlobalTriggerPattern(const AliMUONGlobalTrigger* gloTrg) const -{ - // global trigger pattern calculation - - Int_t gloTrigPat = 0; - - if (gloTrg->SinglePlusLpt()) gloTrigPat|= 0x1; - if (gloTrg->SinglePlusHpt()) gloTrigPat|= 0x2; - if (gloTrg->SinglePlusApt()) gloTrigPat|= 0x4; - - if (gloTrg->SingleMinusLpt()) gloTrigPat|= 0x8; - if (gloTrg->SingleMinusHpt()) gloTrigPat|= 0x10; - if (gloTrg->SingleMinusApt()) gloTrigPat|= 0x20; - - if (gloTrg->SingleUndefLpt()) gloTrigPat|= 0x40; - if (gloTrg->SingleUndefHpt()) gloTrigPat|= 0x80; - if (gloTrg->SingleUndefApt()) gloTrigPat|= 0x100; + fDarcHeader->SetScalersNumbers(); + fRegHeader->SetScalersNumbers(); + fLocalStruct->SetScalersNumbers(); - if (gloTrg->PairUnlikeLpt()) gloTrigPat|= 0x200; - if (gloTrg->PairUnlikeHpt()) gloTrigPat|= 0x400; - if (gloTrg->PairUnlikeApt()) gloTrigPat|= 0x800; - - if (gloTrg->PairLikeLpt()) gloTrigPat|= 0x1000; - if (gloTrg->PairLikeHpt()) gloTrigPat|= 0x2000; - if (gloTrg->PairLikeApt()) gloTrigPat|= 0x4000; - - return gloTrigPat; + fScalerEvent = kTRUE; } +