]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - HLT/MUON/OnlineAnalysis/AliHLTMUONTriggerReconstructorComponent.cxx
Cleanup of documentation and code.
[u/mrichter/AliRoot.git] / HLT / MUON / OnlineAnalysis / AliHLTMUONTriggerReconstructorComponent.cxx
index 9009826c315ec5ec92b0aed1fe3037450bbf8086..2a70a3ef444e5abbe3332a2a48c385c01f3cf752 100644 (file)
 /**************************************************************************
 /**************************************************************************
- * This file is property of and copyright by the ALICE HLT Project        * 
+ * This file is property of and copyright by the ALICE HLT Project        *
  * All rights reserved.                                                   *
  *                                                                        *
  * Primary Authors:                                                       *
  *   Indranil Das <indra.das@saha.ac.in>                                  *
  * All rights reserved.                                                   *
  *                                                                        *
  * Primary Authors:                                                       *
  *   Indranil Das <indra.das@saha.ac.in>                                  *
+ *   Artur Szostak <artursz@iafrica.com>                                  *
  *                                                                        *
  * Permission to use, copy, modify and distribute this software and its   *
  * documentation strictly for non-commercial purposes is hereby granted   *
  * without fee, provided that the above copyright notice appears in all   *
  * copies and that both the copyright notice and this permission notice   *
  * appear in the supporting documentation. The authors make no claims     *
  *                                                                        *
  * Permission to use, copy, modify and distribute this software and its   *
  * documentation strictly for non-commercial purposes is hereby granted   *
  * without fee, provided that the above copyright notice appears in all   *
  * copies and that both the copyright notice and this permission notice   *
  * appear in the supporting documentation. The authors make no claims     *
- * about the suitability of this software for any purpose. It is          * 
+ * about the suitability of this software for any purpose. It is          *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
 /* $Id$ */
 
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
 /* $Id$ */
 
-/** @file   AliHLTMUONTriggerReconstructorComponent.cxx
-    @author Indranil Das
-    @date   
-    @brief  A processing component for the dHLT TrigRec. */
+///
+/// @file   AliHLTMUONTriggerReconstructorComponent.cxx
+/// @author Indranil Das <indra.das@saha.ac.in>, Artur Szostak <artursz@iafrica.com>
+/// @date
+/// @brief  Implementation of the trigger DDL reconstructor component.
+///
 
 
-#if __GNUC__ >= 3
-using namespace std;
-#endif
-
-#include "AliHLTSystem.h"
 #include "AliHLTMUONTriggerReconstructorComponent.h"
 #include "AliHLTMUONTriggerReconstructorComponent.h"
-#include "AliHLTDefinitions.h"
-#include <stdlib.h>
-#include <errno.h>
+#include "AliHLTMUONTriggerReconstructor.h"
+#include "AliHLTMUONHitReconstructor.h"
+#include "AliHLTMUONConstants.h"
+#include "AliHLTMUONUtils.h"
+#include "AliHLTMUONDataBlockWriter.h"
+#include "AliCDBManager.h"
+#include "AliCDBStorage.h"
+#include "AliGeomManager.h"
+#include "AliMUONGeometryTransformer.h"
+#include "AliMUONGeometryDetElement.h"
+#include "AliMpCDB.h"
+#include "AliMpDDLStore.h"
+#include "AliMpPad.h"
+#include "AliMpSegmentation.h"
+#include "AliMpDEIterator.h"
+#include "AliMpVSegmentation.h"
+#include "AliMpDEManager.h"
+#include "AliMpLocalBoard.h"
+#include "AliMpTriggerCrate.h"
+#include <cstdlib>
+#include <cerrno>
+#include <cassert>
+#include <fstream>
 
 
-// this is a global object used for automatic component registration, do not use this
-AliHLTMUONTriggerReconstructorComponent gAliHLTMUONTriggerReconstructorComponent;
 
 ClassImp(AliHLTMUONTriggerReconstructorComponent)
 
 ClassImp(AliHLTMUONTriggerReconstructorComponent)
-    
-AliHLTMUONTriggerReconstructorComponent::AliHLTMUONTriggerReconstructorComponent()
-  :
-  fOutputPercentage(100), // By default we copy to the output exactly what we got as input
-  fDDLDir(""),
-  fDDL(0),
-  fTrigRec(NULL)
-    {
-    }
+
+
+AliHLTMUONTriggerReconstructorComponent::AliHLTMUONTriggerReconstructorComponent() :
+       AliHLTProcessor(),
+       fTrigRec(NULL),
+       fDDL(-1),
+       fWarnForUnexpecedBlock(false),
+       fSuppressPartialTrigs(false)
+{
+       ///
+       /// Default constructor.
+       ///
+}
+
 
 AliHLTMUONTriggerReconstructorComponent::~AliHLTMUONTriggerReconstructorComponent()
 
 AliHLTMUONTriggerReconstructorComponent::~AliHLTMUONTriggerReconstructorComponent()
-    {
-    }
+{
+       ///
+       /// Default destructor.
+       ///
+       
+       if (fTrigRec != NULL) delete fTrigRec;
+}
+
 
 const char* AliHLTMUONTriggerReconstructorComponent::GetComponentID()
 
 const char* AliHLTMUONTriggerReconstructorComponent::GetComponentID()
-    {
-    return "MUONTrigRec"; // The ID of this component
-    }
-
-void AliHLTMUONTriggerReconstructorComponent::GetInputDataTypes( vector<AliHLTComponentDataType>& list)
-    {
-      /* in order to be backward compatible we have to keep the old code, at
-       * least for a while. Remember to use the new const kAliHLTVoidDataType
-       * if you are using a more recent AliRoot version (from Jan 07)
-       list.push_back(kAliHLTAnyDataType); // We do not have any requirements for our input data type(s).
-      */
-
-      list.clear();
-      list.push_back( AliHLTMUONConstants::TriggerDDLRawDataType() );
-    }
+{
+       ///
+       /// Inherited from AliHLTComponent. Returns the component ID.
+       ///
+       
+       return AliHLTMUONConstants::TriggerReconstructorId();
+}
+
+
+void AliHLTMUONTriggerReconstructorComponent::GetInputDataTypes( std::vector<AliHLTComponentDataType>& list)
+{
+       ///
+       /// Inherited from AliHLTProcessor. Returns the list of expected input data types.
+       ///
+       
+       list.clear();
+       list.push_back( AliHLTMUONConstants::DDLRawDataType() );
+}
+
 
 AliHLTComponentDataType AliHLTMUONTriggerReconstructorComponent::GetOutputDataType()
 
 AliHLTComponentDataType AliHLTMUONTriggerReconstructorComponent::GetOutputDataType()
-    {
-      /* in order to be backward compatible we have to keep the old code, at
-       * least for a while. Remember to use the new const kAliHLTVoidDataType
-       * if you are using a more recent AliRoot version (from Jan 07)
-      return kAliHLTVoidDataType;
-      */
-      return AliHLTMUONConstants::TriggerRecordsBlockDataType();
-    }
+{
+       ///
+       /// Inherited from AliHLTComponent. Returns the output data type.
+       ///
+       
+       return AliHLTMUONConstants::TriggerRecordsBlockDataType();
+}
 
 
-void AliHLTMUONTriggerReconstructorComponent::GetOutputDataSize( unsigned long& constBase, double& inputMultiplier )
-    {
-    constBase = 0;
-    inputMultiplier = ((double)fOutputPercentage)/100.0;
-    }
 
 
+void AliHLTMUONTriggerReconstructorComponent::GetOutputDataSize(
+               unsigned long& constBase, double& inputMultiplier
+       )
+{
+       ///
+       /// Inherited from AliHLTComponent. Returns an estimate of the expected output data size.
+       ///
+       
+       constBase = sizeof(AliHLTMUONTriggerRecordsBlockWriter::HeaderType);
+       inputMultiplier = 4;
+}
 
 
 
 
-// Spawn function, return new instance of this class
 AliHLTComponent* AliHLTMUONTriggerReconstructorComponent::Spawn()
 AliHLTComponent* AliHLTMUONTriggerReconstructorComponent::Spawn()
-    {
-    return new AliHLTMUONTriggerReconstructorComponent;
-    }
+{
+       ///
+       /// Inherited from AliHLTComponent. Creates a new object instance.
+       ///
+       
+       return new AliHLTMUONTriggerReconstructorComponent;
+}
 
 
-int AliHLTMUONTriggerReconstructorComponent::DoInit( int argc, const char** argv )
+
+int AliHLTMUONTriggerReconstructorComponent::DoInit(int argc, const char** argv)
 {
 {
-    // perform initialization. We check whether our relative output size is specified in the arguments.
-      
-      fTrigRec = new AliHLTMUONTriggerReconstructor();
-      
-      HLTInfo("dHLT trigrec");
-      if (argc==0 && argv==NULL) {
-       Logging( kHLTLogError, "AliHLTMUONTriggerReconstructorComponent::DoInit", "Arguments missing", " no arguments" );
-       // this is just to get rid of the warning "unused parameter"
-      }
-
-      //Int_t i = 0;
-      char lutFileName[500], ddlDir[500];
-
-      fOutputPercentage = 100;
-      int i = 0;
-      char* cpErr;
-      while ( i < argc )
+       ///
+       /// Inherited from AliHLTComponent.
+       /// Parses the command line parameters and initialises the component.
+       ///
+       
+       // perform initialization.
+       
+       HLTInfo("Initialising dHLT trigger reconstructor component.");
+       
+       // Make sure to cleanup fTrigRec if it is still there for some reason.
+       if (fTrigRec != NULL)
+       {
+               delete fTrigRec;
+               fTrigRec = NULL;
+       }
+       
+       try
        {
        {
-         Logging( kHLTLogDebug, "HLT::MUONTrigRec::DoInit", "Arguments", "argv[%d] == %s", i, argv[i] );
-         if ( !strcmp( argv[i], "output_percentage" ) )
-           {
-             if ( i+1>=argc )
+               fTrigRec = new AliHLTMUONTriggerReconstructor();
+       }
+       catch (const std::bad_alloc&)
+       {
+               HLTError("Could not allocate more memory for the trigger reconstructor component.");
+               return -ENOMEM;
+       }
+       
+       fDDL = -1;
+       fWarnForUnexpecedBlock = false;
+       fSuppressPartialTrigs = false;
+       
+       const char* lutFileName = NULL;
+       const char* cdbPath = NULL;
+       Int_t run = -1;
+       bool useCDB = false;
+       
+       for (int i = 0; i < argc; i++)
+       {
+               if (strcmp( argv[i], "-lut" ) == 0)
+               {
+                       if ( argc <= i+1 )
+                       {
+                               HLTError("LookupTable filename not specified." );
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+                       
+                       lutFileName = argv[i+1];
+                       
+                       i++;
+                       continue;
+               }
+               
+               if (strcmp( argv[i], "-ddl" ) == 0)
                {
                {
-                 Logging(kHLTLogError, "HLT::MUONTrigRec::DoInit", "Missing Argument", "Missing output_percentage parameter");
-                 return ENOTSUP;
+                       if ( argc <= i+1 )
+                       {
+                               HLTError("DDL number not specified. It must be in the range [21..22]" );
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+               
+                       char* cpErr = NULL;
+                       unsigned long num = strtoul(argv[i+1], &cpErr, 0);
+                       if (cpErr == NULL or *cpErr != '\0')
+                       {
+                               HLTError("Cannot convert '%s' to a DDL Number.", argv[i+1] );\
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+                       if (num < 21 or 22 < num)
+                       {
+                               HLTError("The DDL ID number must be in the range [21..22].");
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+                       fDDL = num - 1; // Convert to DDL number in the range 0..21
+                       
+                       i++;
+                       continue;
                }
                }
-             Logging( kHLTLogDebug, "HLT::MUONTrigRec::DoInit", "Arguments", "argv[%d+1] == %s", i, argv[i+1] );
-             fOutputPercentage = strtoul( argv[i+1], &cpErr, 0 );
-             if ( *cpErr )
+               
+               if (strcmp( argv[i], "-cdb" ) == 0)
+               {
+                       useCDB = true;
+                       continue;
+               }
+               
+               if (strcmp( argv[i], "-cdbpath" ) == 0)
+               {
+                       if ( argc <= i+1 )
+                       {
+                               HLTError("The CDB path was not specified." );
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+                       cdbPath = argv[i+1];
+                       useCDB = true;
+                       i++;
+                       continue;
+               }
+       
+               if (strcmp( argv[i], "-run" ) == 0)
                {
                {
-                 Logging(kHLTLogError, "HLT::MUONTrigRec::DoInit", "Wrong Argument", "Cannot convert output_percentage parameter '%s'", argv[i+1] );
-                 return EINVAL;
+                       if ( argc <= i+1 )
+                       {
+                               HLTError("The RUN number was not specified." );
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+                       
+                       char* cpErr = NULL;
+                       run = Int_t( strtoul(argv[i+1], &cpErr, 0) );
+                       if (cpErr == NULL or *cpErr != '\0')
+                       {
+                               HLTError("Cannot convert '%s' to a valid run number."
+                                       " Expected an integer value.", argv[i+1]
+                               );
+                               // Make sure to delete fTrigRec to avoid partial initialisation.
+                               delete fTrigRec;
+                               fTrigRec = NULL;
+                               return -EINVAL;
+                       }
+                       
+                       i++;
+                       continue;
                }
                }
-             Logging( kHLTLogInfo, "HLT::MUONTrigRec::DoInit", "Output percentage set", "Output percentage set to %lu %%", fOutputPercentage );
-             i += 2;
-           continue;
-           }
-       
-
-         if ( !strcmp( argv[i], "lut" ) ) {
-           if ( argc <= i+1 ) {
-             Logging( kHLTLogError, "AliHLTMUONTriggerReconstructorComponent::DoInit", "Missing LookupTable filename", "LookupTable filename not specified" );
-             return EINVAL; /* Invalid argument */ 
-           }
-           
-           sprintf(lutFileName,"%s",argv[i+1]);
-           
-           i += 2;
-           continue;
-         }// lut argument
-         
-         
-         if ( !strcmp( argv[i], "ddl" ) ) {
-           if ( argc <= i+1 ) {
-             Logging( kHLTLogError, "AliHLTMUONTriggerReconstructorComponent::DoInit", "Missing DDL argument", "DDL number not specified" );
-             HLTError("AliHLTMUONTriggerReconstructorComponent::DoInit : DDL number is not specified ");
-             return EINVAL;  /* Invalid argument */
-           }
-           
-           fDDL = atoi(argv[i+1]);
-           
-           i += 2;
-           continue;
-         }// ddl argument
-         
-
-         if ( !strcmp( argv[i], "rawdir" ) ) {
-           if ( argc <= i+1 ) {
-             Logging( kHLTLogError, "AliHLTMUONTriggerReconstructorComponent::DoInit", "Missing DDL directory", "DDL directory not specified" );
-             HLTError("AliHLTMUONTriggerReconstructorComponent::DoInit : DDL directory is not specified ");
-             return EINVAL;  /* Invalid argument */
-           }
-           
-           fDDLDir = argv[i+1] ;
-           
-           i += 2;
-           continue;
-         }// ddl directory argument
-
-         Logging(kHLTLogError, "HLT::MUONTrigRec::DoInit", "Unknown Option", "Unknown option '%s'", argv[i] );
-         return EINVAL;
-         
-       }//while loop
-
-    int lutline = fTrigRec->GetLutLine();
-    AliHLTMUONHitReconstructor::DHLTLut* lookupTable = new AliHLTMUONHitReconstructor::DHLTLut[lutline];
-    if(!ReadLookUpTable(lookupTable,lutFileName)){
-      Logging(kHLTLogInfo, "AliHLTMUONTriggerReconstructorComponent::DoInit", "Failed to read lut", "lut cannot be read, DoInit");
-      return ENOENT ; /* No such file or directory */
-    }else{
-      
-      fTrigRec->LoadLookUpTable(lookupTable,fDDL);
-      
-    }// reading lut
-
-    delete []lookupTable;
-
-    return 0;
+               
+               if (strcmp( argv[i], "-warn_on_unexpected_block" ) == 0)
+               {
+                       fWarnForUnexpecedBlock = true;
+                       continue;
+               }
+               
+               if (strcmp( argv[i], "-suppress_partial_triggers" ) == 0)
+               {
+                       fSuppressPartialTrigs = true;
+                       continue;
+               }
+               
+               HLTError("Unknown option '%s'.", argv[i] );
+               // Make sure to delete fTrigRec to avoid partial initialisation.
+               delete fTrigRec;
+               fTrigRec = NULL;
+               return -EINVAL;
+               
+       } // for loop
+       
+       if (lutFileName == NULL) useCDB = true;
+       
+       if (fDDL == -1)
+       {
+               HLTWarning("DDL number not specified. Cannot check if incomming data is valid.");
+       }
+       
+       int result = 0;
+       if (useCDB)
+       {
+               HLTInfo("Loading lookup table information from CDB for DDL %d.", fDDL+1);
+               if (fDDL == -1)
+                       HLTWarning("DDL number not specified. The lookup table loaded from CDB will be empty!");
+               result = ReadCDB(cdbPath, run);
+       }
+       else
+       {
+               HLTInfo("Loading lookup table information from file %s.", lutFileName);
+               result = ReadLookUpTable(lutFileName);
+       }
+       if (result != 0)
+       {
+               // Error messages already generated in ReadCDB or ReadLookUpTable.
+               
+               // Make sure to delete fTrigRec to avoid partial initialisation.
+               delete fTrigRec;
+               fTrigRec = NULL;
+               return result;
+       }
+       
+       return 0;
 }
 
 }
 
+
 int AliHLTMUONTriggerReconstructorComponent::DoDeinit()
 int AliHLTMUONTriggerReconstructorComponent::DoDeinit()
-    {
-      if(fTrigRec)
-       delete fTrigRec;
-      HLTInfo("dHLT trigrec");
-      return 0;
-  
-    return 0;
-    }
-
-int AliHLTMUONTriggerReconstructorComponent::DoEvent( const AliHLTComponentEventData& evtData, const AliHLTComponentBlockData* blocks, 
-                                     AliHLTComponentTriggerData& trigData, AliHLTUInt8_t* outputPtr, 
-                                     AliHLTUInt32_t& size, vector<AliHLTComponentBlockData>& outputBlocks )
-    {
-      // Process an event
-    unsigned long totalSize = 0;
-    Logging( kHLTLogInfo, "HLT::MUONTrigRec::DoEvent", "Output percentage set", "Output percentage set to %lu %% and totalSize %lu", fOutputPercentage,totalSize );
-
-    // Loop over all input blocks in the event
-    for ( unsigned long n = 0; n < evtData.fBlockCnt; n++ )
-      {
-
-       if ( totalSize > size )
-           break;
-
-       int totalDDLSize = blocks[n].fSize/4;
-       int  ddlRawDataSize = totalDDLSize - AliHLTMUONTriggerReconstructor::GetkDDLHeaderSize();
-       //cout<<"ddlRawDataSize :"<<ddlRawDataSize<<endl;
-       int *buffer = (int *)((int *)blocks[n].fPtr + AliHLTMUONTriggerReconstructor::GetkDDLHeaderSize()) ;
-       
-       AliHLTMUONTriggerRecordStruct trigRecord[300];
-       int nofTrigRec = 300;
-       
-       if(! (fTrigRec->Run(buffer,&ddlRawDataSize,trigRecord,&nofTrigRec))){
-         HLTError("ERROR In Processing of TrigRec Algo ");
-         return EIO;
+{
+       ///
+       /// Inherited from AliHLTComponent. Performs a cleanup of the component.
+       ///
+       
+       HLTInfo("Deinitialising dHLT trigger reconstructor component.");
+
+       if (fTrigRec != NULL)
+       {
+               delete fTrigRec;
+               fTrigRec = NULL;
+       }
+       return 0;
+}
+
+
+int AliHLTMUONTriggerReconstructorComponent::DoEvent(
+               const AliHLTComponentEventData& evtData,
+               const AliHLTComponentBlockData* blocks,
+               AliHLTComponentTriggerData& /*trigData*/,
+               AliHLTUInt8_t* outputPtr,
+               AliHLTUInt32_t& size,
+               std::vector<AliHLTComponentBlockData>& outputBlocks
+       )
+{
+       ///
+       /// Inherited from AliHLTProcessor. Processes the new event data.
+       ///
+       
+       // Process an event
+       unsigned long totalSize = 0; // Amount of memory currently consumed in bytes.
+
+       HLTDebug("Processing event %llu with %u input data blocks.",
+               evtData.fEventID, evtData.fBlockCnt
+       );
+       
+       // Loop over all input blocks in the event and run the trigger DDL
+       // reconstruction algorithm on the raw data.
+       for (AliHLTUInt32_t n = 0; n < evtData.fBlockCnt; n++)
+       {
+               HLTDebug("Handling block: %u, with fDataType = '%s', fPtr = %p and fSize = %u bytes.",
+                       n, DataType2Text(blocks[n].fDataType).c_str(), blocks[n].fPtr, blocks[n].fSize
+               );
+
+               if (blocks[n].fDataType != AliHLTMUONConstants::DDLRawDataType()
+                   or not AliHLTMUONUtils::IsTriggerDDL(blocks[n].fSpecification)
+                  )
+               {
+                       // Log a message indicating that we got a data block that we
+                       // do not know how to handle.
+                       if (fWarnForUnexpecedBlock)
+                               HLTWarning("Received a data block of a type we cannot handle: '%s', spec: 0x%X",
+                                       DataType2Text(blocks[n].fDataType).c_str(), blocks[n].fSpecification
+                               );
+                       else
+                               HLTDebug("Received a data block of a type we cannot handle: '%s', spec: 0x%X",
+                                       DataType2Text(blocks[n].fDataType).c_str(), blocks[n].fSpecification
+                               );
+                       
+                       continue;
+               }
+               
+               if (fDDL != -1)
+               {
+                       bool ddl[22];
+                       AliHLTMUONUtils::UnpackSpecBits(blocks[n].fSpecification, ddl);
+                       if (not ddl[fDDL])
+                       {
+                               HLTWarning("Received raw data from an unexpected DDL.");
+                       }
+               }
+               
+               // Create a new output data block and initialise the header.
+               AliHLTMUONTriggerRecordsBlockWriter block(outputPtr+totalSize, size-totalSize);
+               if (not block.InitCommonHeader())
+               {
+                       HLTError("There is not enough space in the output buffer for the new data block.",
+                                " We require at least %ufTrigRec->GetkDDLHeaderSize() bytes, but have %u bytes left.",
+                               sizeof(AliHLTMUONTriggerRecordsBlockWriter::HeaderType),
+                               block.BufferSize()
+                       );
+                       break;
+               }
+
+               AliHLTUInt32_t totalDDLSize = blocks[n].fSize / sizeof(AliHLTUInt32_t);
+               AliHLTUInt32_t ddlRawDataSize = totalDDLSize - 8;
+               AliHLTUInt32_t* buffer = reinterpret_cast<AliHLTUInt32_t*>(blocks[n].fPtr) + 8;
+               AliHLTUInt32_t nofTrigRec = block.MaxNumberOfEntries();
+
+               bool runOk = fTrigRec->Run(
+                               buffer, ddlRawDataSize,
+                               block.GetArray(), nofTrigRec,
+                               fSuppressPartialTrigs
+                       );
+               if (not runOk)
+               {
+                       HLTError("Error while processing of trigger DDL reconstruction algorithm.");
+                       size = totalSize; // Must tell the framework how much buffer space was used.
+                       return -EIO;
+               }
+               
+               // nofTrigRec should now contain the number of triggers actually found
+               // and filled into the output data block, so we can set this number.
+               assert( nofTrigRec <= block.MaxNumberOfEntries() );
+               block.SetNumberOfEntries(nofTrigRec);
+               
+               HLTDebug("Number of trigger records found is %d", nofTrigRec);
+               
+               // Fill a block data structure for our output block.
+               AliHLTComponentBlockData bd;
+               FillBlockData(bd);
+               bd.fPtr = outputPtr;
+               // This block's start (offset) is after all other blocks written so far.
+               bd.fOffset = totalSize;
+               bd.fSize = block.BytesUsed();
+               bd.fDataType = AliHLTMUONConstants::TriggerRecordsBlockDataType();
+               bd.fSpecification = blocks[n].fSpecification;
+               outputBlocks.push_back(bd);
+               
+               HLTDebug("Created a new output data block at fPtr = %p,"
+                         " with fOffset = %u (0x%.X) and fSize = %u bytes.", 
+                       bd.fPtr, bd.fOffset, bd.fOffset, bd.fSize
+               );
+               
+               // Increase the total amount of data written so far to our output memory.
+               totalSize += block.BytesUsed();
        }
        }
+       
+       // Finally we set the total size of output memory we consumed.
+       size = totalSize;
+       return 0;
+}
+
 
 
-//     if(! (fTrigRec->Run((int)evtData.fEventID,fDDL,trigRecord,&nofTrigRec))){
-//       HLTError("ERROR In Processing of TrigRec Algo ");
-//       return EIO;
-//     }
-       
-       unsigned long mySize = sizeof(AliHLTMUONTriggerRecordStruct)*nofTrigRec;
-    
-//     cout<<"nofHit "<<nofHit<<endl;
-//     for(int i=0;i<nofHit;i++)
-//       cout<<"\t 0 : recHit["<<i<<"].fX :"<<recHit[i].fX
-//           <<"  recHit["<<i<<"].fY :"<<recHit[i].fY
-//           <<"  recHit["<<i<<"].fZ :"<<recHit[i].fZ
-//           <<"  recHit["<<i<<"].fDetElemId :"<<recHit[i].fDetElemId
-//           <<endl;
-
-       //unsigned long mySize = (blocks[n].fSize * fOutputPercentage) / 100;
-
-       Logging( kHLTLogInfo, "HLT::MUONTrigRec::DoEvent", "mySize set (1)", "mySize == %lu B - blocks[%lu].fSize == %lu - fOutputPercentage == %lu", 
-                mySize, n, blocks[n].fSize, fOutputPercentage );
-
-       // Check how much space we have left and adapt this output block's size accordingly.
-       if ( totalSize + mySize > size )
-           mySize = size-totalSize;
-
-       Logging( kHLTLogInfo, "HLT::MUONTrigRec::DoEvent", "mySize set (2)", "mySize == %lu B - totalSize == %lu - size == %lu", 
-                mySize, totalSize, size );
-
-       if ( mySize<=0 )
-           continue; // No room left to write a further block.
-
-       // Now copy the input block
-       unsigned long copied = 0;
-       // First copy all full multiples of the input block
-
-       // And the copy the remaining fragment of the block
-       Logging( kHLTLogInfo, "1 : HLT::MUONTrigRec::DoEvent", "Copying", "Copying %lu B - Copied: %lu B - totalSize: %lu B", 
-                mySize-copied, copied, totalSize );
-       //memcpy( outputPtr+totalSize+copied, blocks[n].fPtr, mySize-copied );
-       memcpy( outputPtr+totalSize+copied, &trigRecord[0], mySize);
-       Logging( kHLTLogInfo, "HLT::MUONTrigRec::DoEvent", "Copied", "Copied: %lu B - totalSize: %lu B", 
-                copied, totalSize );
-       // Fill a block data structure for our output block.
-       AliHLTComponentBlockData ob;
-       // Let the structure be filled with the default values.
-       // This takes care of setting the shared memory and data type values to default values,
-       // so that they can be filled in by the calling code.
-       FillBlockData( ob );
-       // This block's start (offset) is after all other blocks written so far
-       ob.fOffset = totalSize;
-       // the size of this block's data.
-       ob.fSize = mySize;
-       // The specification of the data is copied from the input block.
-       ob.fSpecification = blocks[n].fSpecification;
-       // The data type is set automatically to the component's specified output data type.
-       // Place this block into the list of output blocks
-       outputBlocks.push_back( ob );
-       // Increase the total amount of data written so far to our output memory
-       totalSize += mySize;
+int AliHLTMUONTriggerReconstructorComponent::ReadLookUpTable(const char* lutpath)
+{
+       ///
+       /// Read in the lookup table from file.
+       ///
+       
+       assert(fTrigRec != NULL);
+
+       fstream file;
+       file.open(lutpath, fstream::binary | fstream::in);
+       if (not file)
+       {
+               HLTError("Could not open file: %s", lutpath);
+               return -ENOENT;
        }
        }
-    // Finally we set the total size of output memory we consumed.
-    size = totalSize;
+       
+       file.read(reinterpret_cast<char*>(fTrigRec->LookupTableBuffer()), fTrigRec->LookupTableSize());
+       if (file.eof())
+       {
+               HLTError("The file %s was too short to contain a valid lookup table for this component.", lutpath);
+               file.close();
+               return -EIO;
+       }
+       if (file.fail())
+       {
+               HLTError("Could not read from file: %s", lutpath);
+               file.close();
+               return -EIO;
+       }
+       
+       file.close();
+       return 0;
+}
 
 
-    return 0;
-    }
 
 
-bool AliHLTMUONTriggerReconstructorComponent::ReadLookUpTable(AliHLTMUONHitReconstructor::DHLTLut* lookupTable, const char* lutpath)
+int AliHLTMUONTriggerReconstructorComponent::ReadCDB(const char* cdbPath, Int_t run)
 {
 {
-  if (fDDL < AliHLTMUONTriggerReconstructor::GetkDDLOffSet() ||
-      fDDL >= AliHLTMUONTriggerReconstructor::GetkDDLOffSet() + AliHLTMUONTriggerReconstructor::GetkNofDDL()){
-    HLTError("DDL number is out of range");
-    return false;
-  }
-  
-  int lutLine = fTrigRec->GetLutLine();
-  
-  FILE* fin = fopen(lutpath, "r");
-  if (fin == NULL){
-    HLTError("Failed to open file: %s",lutpath);
-    return false;
-  }
-  
-  for(int i=0;i<lutLine;i++){
-    fscanf(
-          fin,
-          "%d\t%d\t%d\t%f\t%f\t%f\t%d\t%d\n",
-          &lookupTable[i].fIdManuChannel,
-          &lookupTable[i].fIX,
-          &lookupTable[i].fIY,
-          &lookupTable[i].fRealX,
-          &lookupTable[i].fRealY,
-          &lookupTable[i].fRealZ,
-          &lookupTable[i].fPlane,
-          &lookupTable[i].fPcbZone
-          );
-  }
-  
-  fclose(fin);
-  return true;
+       /// Loads the lookup table containing channel and geometrical position
+       /// information about trigger strips from CDB.
+       /// \param cdbPath  This specifies the CDB path to use to load from.
+       ///                 Can be set to NULL meaning the default storage is used.
+       /// \param run  Specifies the run number to use. If set to -1 then the
+       ///             default / current run number set for the CDB is used.
+       /// \return 0 on success and non zero codes for errors.
+
+       if (fDDL == -1)
+       {
+               HLTError("No DDL number specified for which to load LUT data from CDB.");
+               return -EINVAL;
+       }
+
+       Bool_t warn = kFALSE;
+       
+       AliCDBManager* cdbManager = AliCDBManager::Instance();
+       if (cdbManager == NULL)
+       {
+               HLTError("CDB manager instance does not exist.");
+               return -EIO;
+       }
+       
+       const char* cdbPathUsed = "unknown (not set)";
+       if (cdbPath != NULL)
+       {
+               cdbManager->SetDefaultStorage(cdbPath);
+               cdbPathUsed = cdbPath;
+       }
+       else
+       {
+               AliCDBStorage* store = cdbManager->GetDefaultStorage();
+               if (store != NULL) cdbPathUsed = store->GetURI().Data();
+       }
+       
+       if (run != -1) cdbManager->SetRun(run);
+       Int_t runUsed = cdbManager->GetRun();
+       
+       if (not AliMpCDB::LoadDDLStore(warn))
+       {
+               HLTError("Failed to load DDL store specified for CDB path '%s' and run no. %d",
+                       cdbPathUsed, runUsed
+               );
+               return -ENOENT;
+       }
+       AliMpDDLStore* ddlStore = AliMpDDLStore::Instance();
+       if (ddlStore == NULL)
+       {
+               HLTError("Could not find DDL store instance.");
+               return -EIO;
+       }
+       AliMpSegmentation* segmentation = AliMpSegmentation::Instance();
+       if (segmentation == NULL)
+       {
+               HLTError("Could not find segmentation mapping (AliMpSegmentation) instance.");
+               return -EIO;
+       }
+       
+       AliGeomManager::LoadGeometry();
+       AliMUONGeometryTransformer transformer;
+       if (not transformer.LoadGeometryData())
+       {
+               HLTError("Could not load geometry into transformer.");
+               return -ENOENT;
+       }
+       
+       AliHLTMUONTriggerRecoLookupTable* lookupTable = fTrigRec->LookupTableBuffer();
+       
+       for (Int_t i = 0; i < 8; i++)
+       for (Int_t j = 0; j < 16; j++)
+       for (Int_t k = 0; k < 4; k++)
+       for (Int_t n = 0; n < 2; n++)
+       for (Int_t m = 0; m < 16; m++)
+       {
+               lookupTable->fRow[i][j][k][n][m].fX = 0;
+               lookupTable->fRow[i][j][k][n][m].fY = 0;
+               lookupTable->fRow[i][j][k][n][m].fZ = 0;
+       }
+       
+       AliMpDEIterator detElemIter;
+       for (Int_t iReg = 0; iReg < 8; iReg++)
+       {
+               AliMpTriggerCrate* crate = ddlStore->GetTriggerCrate(fDDL, iReg);
+               if (crate == NULL)
+               {
+                       cerr << "ERROR: Could not get crate for regional header = " << iReg
+                               << ", and DDL ID = " << fDDL << endl;
+                       continue;
+               }
+               
+               for (Int_t iLocBoard = 0; iLocBoard < 16; iLocBoard++)
+               {
+                       Int_t boardId = crate->GetLocalBoardId(iLocBoard);
+                       if (boardId == 0) continue;
+                       
+                       AliMpLocalBoard* localBoard = ddlStore->GetLocalBoard(boardId);
+                       if (localBoard == NULL)
+                       {
+                               cerr << "ERROR: Could not get loacl board: " << boardId << endl;
+                               continue;
+                       }
+
+                       // skip copy cards
+                       if (! localBoard->IsNotified()) continue;
+               
+                       for (Int_t iChamber = 0; iChamber < 4; iChamber++)
+                       {
+                               Int_t detElemId = ddlStore->GetDEfromLocalBoard(boardId, iChamber);
+                               
+                               const AliMUONGeometryDetElement* detElemTransform = transformer.GetDetElement(detElemId);
+                               if (detElemTransform == NULL)
+                               {
+                                       cerr << "ERROR: Got NULL pointer for geometry transformer for detection element ID = "
+                                               << detElemId << endl;
+                                       continue;
+                               }
+                               
+                               for (Int_t iCathode = 0; iCathode <= 1; iCathode++)
+                               {
+                                       const AliMpVSegmentation* seg = segmentation->GetMpSegmentation(
+                                                       detElemId, AliMp::GetCathodType(iCathode)
+                                               );
+                                       
+                                       for (Int_t bitxy = 0; bitxy < 16; bitxy++)
+                                       {
+                                               Int_t offset = 0;
+                                               if (iCathode && localBoard->GetSwitch(6)) offset = -8;
+                                               
+                                               AliMpPad pad = seg->PadByLocation(AliMpIntPair(boardId, bitxy+offset), kFALSE);
+                                       
+                                               if (! pad.IsValid())
+                                               {
+                                                       // There is no pad associated with the given local board and bit pattern.
+                                                       continue;
+                                               }
+                                               
+                                               // Get the global coodinates of the pad.
+                                               Float_t lx = pad.Position().X();
+                                               Float_t ly = pad.Position().Y();
+                                               Float_t gx, gy, gz;
+                                               detElemTransform->Local2Global(lx, ly, 0, gx, gy, gz);
+                                               
+                                               // Fill the LUT
+                                               lookupTable->fRow[iReg][iLocBoard][iChamber][iCathode][bitxy].fX = gx;
+                                               lookupTable->fRow[iReg][iLocBoard][iChamber][iCathode][bitxy].fY = gy;
+                                               lookupTable->fRow[iReg][iLocBoard][iChamber][iCathode][bitxy].fZ = gz;
+                                       }
+                               }
+                       }
+               }
+       }
+
+       return 0;
 }
 
 }
 
-// implement this as well.
 
 
+bool AliHLTMUONTriggerReconstructorComponent::GenerateLookupTable(
+               AliHLTInt32_t ddl, const char* filename,
+               const char* cdbPath, Int_t run
+       )
+{
+       /// Generates a binary file containing the lookup table (LUT) from the
+       /// CDB, which can be used for the trigger reconstructor component later.
+       /// @param ddl  Must be the DDL for which to generate the DDL,
+       ///             in the range [20..21].
+       /// @param filename  The name of the LUT file to generate.
+       /// @param cdbPath  The CDB path to use.
+       /// @param run  The run number to use for the CDB.
+       /// @return  True if the generation of the LUT file succeeded.
+       
+       AliHLTMUONTriggerReconstructorComponent comp;
+       
+       if (ddl < 20 or 21 < ddl)
+       {
+               std::cerr << "ERROR: the DDL number must be in the range [20..21]." << std::endl;
+               return false;
+       }
+       
+       char ddlNum[32];
+       char runNum[32];
+       sprintf(ddlNum, "%d", ddl+1);
+       sprintf(runNum, "%d", run);
+       const char* argv[7] = {"-ddl", ddlNum, "-cdbpath", cdbPath, "-run", runNum, NULL};
+       int result = comp.DoInit(6, argv);
+       if (result != 0)
+       {
+               // Error message already generated in DoInit.
+               return false;
+       }
+       
+       std::fstream file(filename, std::ios::out);
+       if (not file)
+       {
+               std::cerr << "ERROR: could not open file: " << filename << std::endl;
+               return false;
+       }
+       
+       file.write(
+                       reinterpret_cast<char*>(comp.fTrigRec->LookupTableBuffer()),
+                       comp.fTrigRec->LookupTableSize()
+               );
+       if (not file)
+       {
+               cerr << "ERROR: There was a problem writing to the file: " << filename << endl;
+               return false;
+       }
+       file.close();
+       
+       comp.DoDeinit();
+       
+       return true;
+}