[u/mrichter/AliRoot.git] / HLT / TPCLib / AliHLTTPCAgent.cxx

// @(#) $Id$

//**************************************************************************
//* This file is property of and copyright by the ALICE HLT Project        * 
//* ALICE Experiment at CERN, All rights reserved.                         *
//*                                                                        *
//* Primary Authors: Matthias Richter <Matthias.Richter@ift.uib.no>        *
//*                  for The ALICE HLT Project.                            *
//*                                                                        *
//* 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          *
//* provided "as is" without express or implied warranty.                  *
//**************************************************************************

/** @file   AliHLTTPCAgent.cxx
    @author Matthias Richter
    @date   
    @brief  Agent of the libAliHLTTPC library
*/

#include "AliHLTTPCAgent.h"
#include "AliHLTConfiguration.h"
#include "AliHLTTPCDefinitions.h"
#include "AliHLTOUT.h"
#include "AliHLTOUTHandlerChain.h"
#include "AliRunLoader.h"
#include "AliCDBManager.h"
#include "AliCDBEntry.h"
#include "AliTPCParam.h"

/** global instance for agent registration */
AliHLTTPCAgent gAliHLTTPCAgent;

// component headers
#include "AliHLTTPCRunStatisticsProducerComponent.h"
#include "AliHLTTPCEventStatisticsProducerComponent.h"
#include "AliHLTTPCCompModelInflaterComponent.h"
#include "AliHLTTPCCompModelDeflaterComponent.h"
#include "AliHLTTPCCompModelDeconverterComponent.h"
#include "AliHLTTPCCompModelConverterComponent.h"
#include "AliHLTTPCCompDumpComponent.h"
//#include "AliHLTTPCCalibCEComponent.h"
//#include "AliHLTTPCCalibPulserComponent.h"
//#include "AliHLTTPCCalibPedestalComponent.h"
#include "AliHLTTPCCAInputDataCompressorComponent.h"
#include "AliHLTTPCCATrackerComponent.h"
#include "AliHLTTPCTrackMCMarkerComponent.h"
#include "AliHLTTPCCAGlobalMergerComponent.h"
#include "AliHLTTPCdEdxComponent.h"
#include "AliHLTTPCGlobalMergerComponent.h"
#include "AliHLTTPCSliceTrackerComponent.h"
#include "AliHLTTPCVertexFinderComponent.h"
#include "AliHLTTPCClusterFinderComponent.h"
#include "AliHLTTPCRawDataUnpackerComponent.h"
#include "AliHLTTPCDigitPublisherComponent.h"
#include "AliHLTTPCZeroSuppressionComponent.h"
#include "AliHLTTPCDigitDumpComponent.h"
#include "AliHLTTPCClusterDumpComponent.h"
#include "AliHLTTPCEsdWriterComponent.h"
#include "AliHLTTPCOfflineClustererComponent.h"
#include "AliHLTTPCOfflineTrackerComponent.h"
#include "AliHLTTPCOfflineTrackerCalibComponent.h"
#include "AliHLTTPCOfflineCalibrationComponent.h" // to be added to the calibration library agent
#include "AliHLTTPCClusterHistoComponent.h"
#include "AliHLTTPCNoiseMapComponent.h"
#include "AliHLTTPCHistogramHandlerComponent.h"
//#include "AliHLTTPCCalibTracksComponent.h"
#include "AliHLTTPCTrackHistoComponent.h"
#include "AliHLTTPCTrackDumpComponent.h"
#include "AliHLTTPCHWCFDataReverterComponent.h"
#include "AliHLTTPCHWClusterTransformComponent.h"
// #include "AliHLTTPCCalibSeedMakerComponent.h"
// #include "AliHLTTPCCalibTimeComponent.h"
// #include "AliHLTTPCCalibTimeGainComponent.h"
// #include "AliHLTTPCCalibrationComponent.h"

/** ROOT macro for the implementation of ROOT specific class methods */
ClassImp(AliHLTTPCAgent)

AliHLTTPCAgent::AliHLTTPCAgent()
  :
  AliHLTModuleAgent("TPC"),
  fRawDataHandler(NULL),
  fTracksegsDataHandler(NULL)
{
  // see header file for class documentation
  // or
  // refer to README to build package
  // or
  // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
}

AliHLTTPCAgent::~AliHLTTPCAgent()
{
  // see header file for class documentation
}

int AliHLTTPCAgent::CreateConfigurations(AliHLTConfigurationHandler* handler,
					 AliRawReader* rawReader,
					 AliRunLoader* runloader) const
{
  // see header file for class documentation
  if (handler) {
    const char* cdbEntry="TPC/Calib/Parameters";
    AliCDBManager* pMan=AliCDBManager::Instance();
    AliTPCParam* pTPCParam=NULL;
    if (pMan) {
      AliCDBEntry *pEntry = pMan->Get(cdbEntry);
      if (pEntry && 
	  pEntry->GetObject() &&
	  (pTPCParam=dynamic_cast<AliTPCParam*>(pEntry->GetObject()))) {
      } else {
	HLTWarning("can not load AliTPCParam from CDB entry %s", cdbEntry);
      }
    }

    // This the tracking configuration for the full TPC
    // - 216 clusterfinders (1 per partition)
    // - 36 slice trackers
    // - one global merger
    // - the esd converter
    // The ESD is shipped embedded into a TTree
    int iMinSlice=0; 
    int iMaxSlice=35;
    int iMinPart=0;
    int iMaxPart=5;
    TString mergerInput;
    TString sinkClusterInput;
    for (int slice=iMinSlice; slice<=iMaxSlice; slice++) {
      TString trackerInput;
      for (int part=iMinPart; part<=iMaxPart; part++) {
	TString arg, publisher, cf;

	// digit publisher components
	publisher.Form("TPC-DP_%02d_%d", slice, part);
	if (rawReader || !runloader) {
	  // AliSimulation: use the AliRawReaderPublisher if the raw reader is available
	  // Alireconstruction: indicated by runloader==NULL, run always on raw data
	  int ddlno=768;
	  if (part>1) ddlno+=72+4*slice+(part-2);
	  else ddlno+=2*slice+part;
	  arg.Form("-minid %d -datatype 'DDL_RAW ' 'TPC '  -dataspec 0x%02x%02x%02x%02x -silent", ddlno, slice, slice, part, part);
	  handler->CreateConfiguration(publisher.Data(), "AliRawReaderPublisher", NULL , arg.Data());
	} else {
	  arg.Form("-slice %d -partition %d", slice, part);
	  handler->CreateConfiguration(publisher.Data(), "TPCDigitPublisher", NULL , arg.Data());
	}

	// cluster finder components
	cf.Form("TPC-CF_%02d_%d", slice, part);
	arg="";
	if (pTPCParam) {
	  arg+=" -timebins "; arg+=pTPCParam->GetMaxTBin()+1;
	}
	if (!rawReader && runloader) {
	  arg+=" -do-mc";
	  handler->CreateConfiguration(cf.Data(), "TPCClusterFinderUnpacked", publisher.Data(), arg.Data());
	} else {
#ifndef HAVE_NOT_ALTRORAWSTREAMV3
	  handler->CreateConfiguration(cf.Data(), "TPCClusterFinder32Bit", publisher.Data(),"");
#else
	  // using the AltroDecoder if the official V3 decoder is not
	  // available in the offline code
	  handler->CreateConfiguration(cf.Data(), "TPCClusterFinderDecoder", publisher.Data(), arg.Data());
#endif 
	}
	if (trackerInput.Length()>0) trackerInput+=" ";
	trackerInput+=cf;
	if (sinkClusterInput.Length()>0) sinkClusterInput+=" ";
	sinkClusterInput+=cf;
      }
      TString tracker;
      // tracker finder components
      tracker.Form("TPC-TR_%02d", slice);
      handler->CreateConfiguration(tracker.Data(), "TPCCATracker", trackerInput.Data(), "");

      if (mergerInput.Length()>0) mergerInput+=" ";
      mergerInput+=tracker;

    }

    // GlobalMerger component
    handler->CreateConfiguration("TPC-globalmerger","TPCCAGlobalMerger",mergerInput.Data(),"");

    // the esd converter configuration
    TString converterInput="TPC-globalmerger";
    if (!rawReader && runloader) {
      // propagate cluster info to the esd converter in order to fill the MC information
      handler->CreateConfiguration("TPC-clustermc-info", "BlockFilter"   , sinkClusterInput.Data(), "-datatype 'CLMCINFO' 'TPC '");  
      handler->CreateConfiguration("TPC-mcTrackMarker","TPCTrackMCMarker","TPC-globalmerger TPC-clustermc-info","" );
      converterInput+=" ";
      converterInput+="TPC-mcTrackMarker";
    }
    handler->CreateConfiguration("TPC-esd-converter", "TPCEsdConverter"   , converterInput.Data(), "");

    // cluster dump collection
    handler->CreateConfiguration("TPC-clusters", "BlockFilter"   , sinkClusterInput.Data(), "-datatype 'CLUSTERS' 'TPC '");

    /////////////////////////////////////////////////////////////////////////////////////
    //
    // a kChain HLTOUT configuration for processing of {'TRAKSEGS':'TPC '} data blocks
    // collects the data blocks, merges the tracks and produces an ESD object

    // publisher component
    handler->CreateConfiguration("TPC-hltout-tracksegs-publisher", "AliHLTOUTPublisher"   , NULL, "");

    // GlobalMerger component
    handler->CreateConfiguration("TPC-hltout-tracksegs-merger", "TPCGlobalMerger", "TPC-hltout-tracksegs-publisher", "");

    // the esd converter configuration
    handler->CreateConfiguration("TPC-hltout-tracksegs-esd-converter", "TPCEsdConverter", "TPC-hltout-tracksegs-merger", "");

    /////////////////////////////////////////////////////////////////////////////////////
    //
    // a kChain HLTOUT configuration for processing of {'TRACKS  ':'TPC '} data blocks
    // produces an ESD object from the track structure

    // publisher component
    handler->CreateConfiguration("TPC-hltout-tracks-publisher", "AliHLTOUTPublisher"   , NULL, "");

    // the esd converter configuration
    handler->CreateConfiguration("TPC-hltout-tracks-esd-converter", "TPCEsdConverter", "TPC-hltout-tracks-publisher", "");
  }
  return 0;
}

const char* AliHLTTPCAgent::GetReconstructionChains(AliRawReader* /*rawReader*/,
						    AliRunLoader* runloader) const
{
  // see header file for class documentation
  if (runloader) {
    // reconstruction chains for AliRoot simulation
    // Note: run loader is only available while running embedded into
    // AliRoot simulation
    //if (runloader->GetLoader("TPCLoader") != NULL)
      //return "TPC-esd-converter TPC-clusters";
      return "TPC-clusters";
  }
  return NULL;
}

const char* AliHLTTPCAgent::GetRequiredComponentLibraries() const
{
  // see header file for class documentation

  // actually, the TPC library has dependencies to Util and RCU
  // so the two has to be loaded anyhow before we get here
  //return "libAliHLTUtil.so libAliHLTRCU.so";
  return "libAliHLTUtil.so";
}

int AliHLTTPCAgent::RegisterComponents(AliHLTComponentHandler* pHandler) const
{
  // see header file for class documentation
  if (!pHandler) return -EINVAL;

  pHandler->AddComponent(new AliHLTTPCRunStatisticsProducerComponent);
  pHandler->AddComponent(new AliHLTTPCEventStatisticsProducerComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibCEComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibPulserComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibPedestalComponent);
  pHandler->AddComponent(new AliHLTTPCCompModelInflaterComponent);
  pHandler->AddComponent(new AliHLTTPCCompModelDeflaterComponent);
  pHandler->AddComponent(new AliHLTTPCCompModelDeconverterComponent);
  pHandler->AddComponent(new AliHLTTPCCompModelConverterComponent);
  pHandler->AddComponent(new AliHLTTPCCompDumpComponent);
  pHandler->AddComponent(new AliHLTTPCCAInputDataCompressorComponent);
  pHandler->AddComponent(new AliHLTTPCCATrackerComponent);
  pHandler->AddComponent(new AliHLTTPCCAGlobalMergerComponent);
  pHandler->AddComponent(new AliHLTTPCTrackMCMarkerComponent);
  pHandler->AddComponent(new AliHLTTPCGlobalMergerComponent);
  pHandler->AddComponent(new AliHLTTPCdEdxComponent);
  pHandler->AddComponent(new AliHLTTPCSliceTrackerComponent);
  pHandler->AddComponent(new AliHLTTPCVertexFinderComponent);
  pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinderPacked));
  pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinderUnpacked));
  pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinderDecoder));
  pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinder32Bit));
  pHandler->AddComponent(new AliHLTTPCRawDataUnpackerComponent);
  pHandler->AddComponent(new AliHLTTPCDigitPublisherComponent);
  pHandler->AddComponent(new AliHLTTPCZeroSuppressionComponent);
  pHandler->AddComponent(new AliHLTTPCDigitDumpComponent);
  pHandler->AddComponent(new AliHLTTPCClusterDumpComponent);
  pHandler->AddComponent(new AliHLTTPCEsdWriterComponent::AliWriter);
  pHandler->AddComponent(new AliHLTTPCEsdWriterComponent::AliConverter);
  pHandler->AddComponent(new AliHLTTPCOfflineClustererComponent);
  pHandler->AddComponent(new AliHLTTPCOfflineTrackerComponent);
  pHandler->AddComponent(new AliHLTTPCOfflineTrackerCalibComponent);
  pHandler->AddComponent(new AliHLTTPCOfflineCalibrationComponent);
  pHandler->AddComponent(new AliHLTTPCClusterHistoComponent);
  pHandler->AddComponent(new AliHLTTPCNoiseMapComponent);
  pHandler->AddComponent(new AliHLTTPCHistogramHandlerComponent);
  //pHandler->AddComponent(new AliHLTTPCCalibTracksComponent);
  pHandler->AddComponent(new AliHLTTPCTrackHistoComponent);
  pHandler->AddComponent(new AliHLTTPCTrackDumpComponent);
  pHandler->AddComponent(new AliHLTTPCHWCFDataReverterComponent);
  pHandler->AddComponent(new AliHLTTPCHWClusterTransformComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibSeedMakerComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibTimeComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibTimeGainComponent);
//   pHandler->AddComponent(new AliHLTTPCCalibrationComponent);

  return 0;
}

int AliHLTTPCAgent::GetHandlerDescription(AliHLTComponentDataType dt,
					  AliHLTUInt32_t spec,
					  AliHLTOUTHandlerDesc& desc) const
{
  // see header file for class documentation

  // raw data blocks to be fed into offline reconstruction
  if (dt==(kAliHLTDataTypeDDLRaw|kAliHLTDataOriginTPC)) {
    int slice=AliHLTTPCDefinitions::GetMinSliceNr(spec);
    int part=AliHLTTPCDefinitions::GetMinPatchNr(spec);
    if (slice==AliHLTTPCDefinitions::GetMaxSliceNr(spec) &&
	part==AliHLTTPCDefinitions::GetMaxPatchNr(spec)) {
      desc=AliHLTOUTHandlerDesc(kRawReader, dt, GetModuleId());
      return 1;
    } else {
      HLTWarning("handler can not process merged data from multiple ddls:"
		 " min slice %d, max slice %d, min part %d, max part %d",
		 slice, AliHLTTPCDefinitions::GetMaxSliceNr(spec),
		 part, AliHLTTPCDefinitions::GetMaxPatchNr(spec));
      return 0;
    }
  }

  // dump for {'CLUSTERS':'TPC '} currently not used any more
  if (dt==AliHLTTPCDefinitions::fgkClustersDataType) {
      desc=AliHLTOUTHandlerDesc(kProprietary, dt, GetModuleId());
      return 1;
  }

  // afterburner for {'TRAKSEGS':'TPC '} blocks to be converted to ESD format
  if (dt==AliHLTTPCDefinitions::fgkTrackSegmentsDataType) {
      desc=AliHLTOUTHandlerDesc(kChain, dt, GetModuleId());
      return 1;
  }

  // afterburner for {'TRACKS  ':'TPC '} block to be converted to ESD format
  // there is only one data block
  if (dt==AliHLTTPCDefinitions::fgkTracksDataType) {
      desc=AliHLTOUTHandlerDesc(kChain, dt, GetModuleId());
      return 1;
  }
  return 0;
}

AliHLTOUTHandler* AliHLTTPCAgent::GetOutputHandler(AliHLTComponentDataType dt,
						   AliHLTUInt32_t /*spec*/)
{
  // see header file for class documentation

  // raw data blocks to be fed into offline reconstruction
  if (dt==(kAliHLTDataTypeDDLRaw|kAliHLTDataOriginTPC)) {
    if (!fRawDataHandler) {
      fRawDataHandler=new AliHLTTPCAgent::AliHLTTPCRawDataHandler;
    }
    return fRawDataHandler;
  }

  // afterburner for {'TRAKSEGS':'TPC '} blocks to be converted to ESD format
  // in a kChain HLTOUT handler
  if (dt==AliHLTTPCDefinitions::fgkTrackSegmentsDataType) {
    if (fTracksegsDataHandler==NULL)
      fTracksegsDataHandler=new AliHLTOUTHandlerChain("chains=TPC-hltout-tracksegs-esd-converter");
    return fTracksegsDataHandler;
  }

  // afterburner for {'TRACKS  ':'TPC '} block to be converted to ESD format
  // there is only one data block
  if (dt==AliHLTTPCDefinitions::fgkTracksDataType) {
    return new AliHLTOUTHandlerChain("chains=TPC-hltout-tracks-esd-converter");
  }

  return NULL;
}

int AliHLTTPCAgent::DeleteOutputHandler(AliHLTOUTHandler* pInstance)
{
  // see header file for class documentation
  if (pInstance==NULL) return -EINVAL;

  if (pInstance==fRawDataHandler) {
    delete fRawDataHandler;
    fRawDataHandler=NULL;
  }

  if (pInstance==fTracksegsDataHandler) {
    delete fTracksegsDataHandler;
    fTracksegsDataHandler=NULL;
  }
  return 0;
}

AliHLTTPCAgent::AliHLTTPCRawDataHandler::AliHLTTPCRawDataHandler()
{
  // see header file for class documentation
}

AliHLTTPCAgent::AliHLTTPCRawDataHandler::~AliHLTTPCRawDataHandler()
{
  // see header file for class documentation
}

int AliHLTTPCAgent::AliHLTTPCRawDataHandler::ProcessData(AliHLTOUT* pData)
{
  // see header file for class documentation
  if (!pData) return -EINVAL;
  AliHLTComponentDataType dt=kAliHLTVoidDataType;
  AliHLTUInt32_t spec=kAliHLTVoidDataSpec;
  int iResult=pData->GetDataBlockDescription(dt, spec);
  if (iResult>=0) {
    int slice=AliHLTTPCDefinitions::GetMinSliceNr(spec);
    int part=AliHLTTPCDefinitions::GetMinPatchNr(spec);
    if (slice==AliHLTTPCDefinitions::GetMaxSliceNr(spec) &&
	part==AliHLTTPCDefinitions::GetMaxPatchNr(spec)) {
      iResult=768;
      if (part>1) iResult+=72+4*slice+(part-2);
      else iResult+=2*slice+part;
    } else {
      HLTError("handler can not process merged data from multiple ddls:"
	       " min slice %d, max slice %d, min part %d, max part %d",
	       slice, AliHLTTPCDefinitions::GetMaxSliceNr(spec),
	       part, AliHLTTPCDefinitions::GetMaxPatchNr(spec));
      iResult=-EBADMSG;
    }
  }
  return iResult;
}
Commit	Line	Data
5578cf60	1	// @(#) $Id$
5578cf60	2
0efebbac	3	//**************************************************************************
	4	//* This file is property of and copyright by the ALICE HLT Project *
	5	//* ALICE Experiment at CERN, All rights reserved. *
	6	//* *
	7	//* Primary Authors: Matthias Richter <Matthias.Richter@ift.uib.no> *
	8	//* for The ALICE HLT Project. *
	9	//* *
	10	//* Permission to use, copy, modify and distribute this software and its *
	11	//* documentation strictly for non-commercial purposes is hereby granted *
	12	//* without fee, provided that the above copyright notice appears in all *
	13	//* copies and that both the copyright notice and this permission notice *
	14	//* appear in the supporting documentation. The authors make no claims *
	15	//* about the suitability of this software for any purpose. It is *
	16	//* provided "as is" without express or implied warranty. *
	17	//**************************************************************************
5578cf60	18
	19	/** @file AliHLTTPCAgent.cxx
	20	@author Matthias Richter
	21	@date
	22	@brief Agent of the libAliHLTTPC library
	23	*/
	24
	25	#include "AliHLTTPCAgent.h"
	26	#include "AliHLTConfiguration.h"
626bfcc1	27	#include "AliHLTTPCDefinitions.h"
626bfcc1	28	#include "AliHLTOUT.h"
73449074	29	#include "AliHLTOUTHandlerChain.h"
352670bf	30	#include "AliRunLoader.h"
286c9940	31	#include "AliCDBManager.h"
	32	#include "AliCDBEntry.h"
	33	#include "AliTPCParam.h"
5578cf60	34
	35	/** global instance for agent registration */
	36	AliHLTTPCAgent gAliHLTTPCAgent;
	37
e1440dab	38	// component headers
672f8b8c	39	#include "AliHLTTPCRunStatisticsProducerComponent.h"
	40	#include "AliHLTTPCEventStatisticsProducerComponent.h"
	41	#include "AliHLTTPCCompModelInflaterComponent.h"
	42	#include "AliHLTTPCCompModelDeflaterComponent.h"
	43	#include "AliHLTTPCCompModelDeconverterComponent.h"
	44	#include "AliHLTTPCCompModelConverterComponent.h"
	45	#include "AliHLTTPCCompDumpComponent.h"
32e0c022	46	//#include "AliHLTTPCCalibCEComponent.h"
	47	//#include "AliHLTTPCCalibPulserComponent.h"
	48	//#include "AliHLTTPCCalibPedestalComponent.h"
751d16ac	49	#include "AliHLTTPCCAInputDataCompressorComponent.h"
400d857a	50	#include "AliHLTTPCCATrackerComponent.h"
a6fb14f5	51	#include "AliHLTTPCTrackMCMarkerComponent.h"
400d857a	52	#include "AliHLTTPCCAGlobalMergerComponent.h"
0973c527	53	#include "AliHLTTPCdEdxComponent.h"
672f8b8c	54	#include "AliHLTTPCGlobalMergerComponent.h"
	55	#include "AliHLTTPCSliceTrackerComponent.h"
	56	#include "AliHLTTPCVertexFinderComponent.h"
	57	#include "AliHLTTPCClusterFinderComponent.h"
	58	#include "AliHLTTPCRawDataUnpackerComponent.h"
	59	#include "AliHLTTPCDigitPublisherComponent.h"
	60	#include "AliHLTTPCZeroSuppressionComponent.h"
e1440dab	61	#include "AliHLTTPCDigitDumpComponent.h"
0efebbac	62	#include "AliHLTTPCClusterDumpComponent.h"
f32b83e1	63	#include "AliHLTTPCEsdWriterComponent.h"
1ac82ce6	64	#include "AliHLTTPCOfflineClustererComponent.h"
1ac82ce6	65	#include "AliHLTTPCOfflineTrackerComponent.h"
0e442a0a	66	#include "AliHLTTPCOfflineTrackerCalibComponent.h"
32e0c022	67	#include "AliHLTTPCOfflineCalibrationComponent.h" // to be added to the calibration library agent
0efebbac	68	#include "AliHLTTPCClusterHistoComponent.h"
3a3550fa	69	#include "AliHLTTPCNoiseMapComponent.h"
3a3550fa	70	#include "AliHLTTPCHistogramHandlerComponent.h"
32e0c022	71	//#include "AliHLTTPCCalibTracksComponent.h"
22f34176	72	#include "AliHLTTPCTrackHistoComponent.h"
22f34176	73	#include "AliHLTTPCTrackDumpComponent.h"
a912b63b	74	#include "AliHLTTPCHWCFDataReverterComponent.h"
6af4f584	75	#include "AliHLTTPCHWClusterTransformComponent.h"
32e0c022	76	// #include "AliHLTTPCCalibSeedMakerComponent.h"
	77	// #include "AliHLTTPCCalibTimeComponent.h"
	78	// #include "AliHLTTPCCalibTimeGainComponent.h"
	79	// #include "AliHLTTPCCalibrationComponent.h"
e1440dab	80
5578cf60	81	/** ROOT macro for the implementation of ROOT specific class methods */
	82	ClassImp(AliHLTTPCAgent)
	83
	84	AliHLTTPCAgent::AliHLTTPCAgent()
626bfcc1	85	:
	86	AliHLTModuleAgent("TPC"),
	87	fRawDataHandler(NULL),
73449074	88	fTracksegsDataHandler(NULL)
5578cf60	89	{
	90	// see header file for class documentation
	91	// or
	92	// refer to README to build package
	93	// or
	94	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
	95	}
	96
	97	AliHLTTPCAgent::~AliHLTTPCAgent()
	98	{
	99	// see header file for class documentation
	100	}
	101
	102	int AliHLTTPCAgent::CreateConfigurations(AliHLTConfigurationHandler* handler,
4fc08fff	103	AliRawReader* rawReader,
7c6a3cd3	104	AliRunLoader* runloader) const
5578cf60	105	{
	106	// see header file for class documentation
	107	if (handler) {
286c9940	108	const char* cdbEntry="TPC/Calib/Parameters";
	109	AliCDBManager* pMan=AliCDBManager::Instance();
	110	AliTPCParam* pTPCParam=NULL;
	111	if (pMan) {
	112	AliCDBEntry *pEntry = pMan->Get(cdbEntry);
	113	if (pEntry &&
	114	pEntry->GetObject() &&
	115	(pTPCParam=dynamic_cast<AliTPCParam*>(pEntry->GetObject()))) {
	116	} else {
	117	HLTWarning("can not load AliTPCParam from CDB entry %s", cdbEntry);
	118	}
	119	}
672f8b8c	120
	121	// This the tracking configuration for the full TPC
	122	// - 216 clusterfinders (1 per partition)
	123	// - 36 slice trackers
	124	// - one global merger
	125	// - the esd converter
	126	// The ESD is shipped embedded into a TTree
dfa1592e	127	int iMinSlice=0;
672f8b8c	128	int iMaxSlice=35;
dfa1592e	129	int iMinPart=0;
672f8b8c	130	int iMaxPart=5;
672f8b8c	131	TString mergerInput;
84bb5f35	132	TString sinkClusterInput;
dfa1592e	133	for (int slice=iMinSlice; slice<=iMaxSlice; slice++) {
	134	TString trackerInput;
	135	for (int part=iMinPart; part<=iMaxPart; part++) {
	136	TString arg, publisher, cf;
	137
	138	// digit publisher components
0ba35c53	139	publisher.Form("TPC-DP_%02d_%d", slice, part);
7c6a3cd3	140	if (rawReader \|\| !runloader) {
	141	// AliSimulation: use the AliRawReaderPublisher if the raw reader is available
	142	// Alireconstruction: indicated by runloader==NULL, run always on raw data
4fc08fff	143	int ddlno=768;
	144	if (part>1) ddlno+=72+4*slice+(part-2);
	145	else ddlno+=2*slice+part;
	146	arg.Form("-minid %d -datatype 'DDL_RAW ' 'TPC ' -dataspec 0x%02x%02x%02x%02x -silent", ddlno, slice, slice, part, part);
	147	handler->CreateConfiguration(publisher.Data(), "AliRawReaderPublisher", NULL , arg.Data());
7c6a3cd3	148	} else {
	149	arg.Form("-slice %d -partition %d", slice, part);
	150	handler->CreateConfiguration(publisher.Data(), "TPCDigitPublisher", NULL , arg.Data());
4fc08fff	151	}
dfa1592e	152
dfa1592e	153	// cluster finder components
0ba35c53	154	cf.Form("TPC-CF_%02d_%d", slice, part);
286c9940	155	arg="";
	156	if (pTPCParam) {
	157	arg+=" -timebins "; arg+=pTPCParam->GetMaxTBin()+1;
	158	}
9c790772	159	if (!rawReader && runloader) {
fc7c0369	160	arg+=" -do-mc";
286c9940	161	handler->CreateConfiguration(cf.Data(), "TPCClusterFinderUnpacked", publisher.Data(), arg.Data());
4fc08fff	162	} else {
9c790772	163	#ifndef HAVE_NOT_ALTRORAWSTREAMV3
c1e9d798	164	handler->CreateConfiguration(cf.Data(), "TPCClusterFinder32Bit", publisher.Data(),"");
9c790772	165	#else
	166	// using the AltroDecoder if the official V3 decoder is not
	167	// available in the offline code
286c9940	168	handler->CreateConfiguration(cf.Data(), "TPCClusterFinderDecoder", publisher.Data(), arg.Data());
9c790772	169	#endif
4fc08fff	170	}
dfa1592e	171	if (trackerInput.Length()>0) trackerInput+=" ";
dfa1592e	172	trackerInput+=cf;
84bb5f35	173	if (sinkClusterInput.Length()>0) sinkClusterInput+=" ";
84bb5f35	174	sinkClusterInput+=cf;
dfa1592e	175	}
	176	TString tracker;
	177	// tracker finder components
0ba35c53	178	tracker.Form("TPC-TR_%02d", slice);
fc7c0369	179	handler->CreateConfiguration(tracker.Data(), "TPCCATracker", trackerInput.Data(), "");
dfa1592e	180
672f8b8c	181	if (mergerInput.Length()>0) mergerInput+=" ";
672f8b8c	182	mergerInput+=tracker;
dfa1592e	183
dfa1592e	184	}
5578cf60	185
672f8b8c	186	// GlobalMerger component
fc7c0369	187	handler->CreateConfiguration("TPC-globalmerger","TPCCAGlobalMerger",mergerInput.Data(),"");
672f8b8c	188
672f8b8c	189	// the esd converter configuration
791fbeb7	190	TString converterInput="TPC-globalmerger";
9c790772	191	if (!rawReader && runloader) {
791fbeb7	192	// propagate cluster info to the esd converter in order to fill the MC information
6595fdc6	193	handler->CreateConfiguration("TPC-clustermc-info", "BlockFilter" , sinkClusterInput.Data(), "-datatype 'CLMCINFO' 'TPC '");
6595fdc6	194	handler->CreateConfiguration("TPC-mcTrackMarker","TPCTrackMCMarker","TPC-globalmerger TPC-clustermc-info","" );
791fbeb7	195	converterInput+=" ";
6595fdc6	196	converterInput+="TPC-mcTrackMarker";
791fbeb7	197	}
791fbeb7	198	handler->CreateConfiguration("TPC-esd-converter", "TPCEsdConverter" , converterInput.Data(), "");
73449074	199
84bb5f35	200	// cluster dump collection
	201	handler->CreateConfiguration("TPC-clusters", "BlockFilter" , sinkClusterInput.Data(), "-datatype 'CLUSTERS' 'TPC '");
	202
73449074	203	/////////////////////////////////////////////////////////////////////////////////////
	204	//
	205	// a kChain HLTOUT configuration for processing of {'TRAKSEGS':'TPC '} data blocks
	206	// collects the data blocks, merges the tracks and produces an ESD object
	207
	208	// publisher component
	209	handler->CreateConfiguration("TPC-hltout-tracksegs-publisher", "AliHLTOUTPublisher" , NULL, "");
	210
	211	// GlobalMerger component
	212	handler->CreateConfiguration("TPC-hltout-tracksegs-merger", "TPCGlobalMerger", "TPC-hltout-tracksegs-publisher", "");
	213
	214	// the esd converter configuration
	215	handler->CreateConfiguration("TPC-hltout-tracksegs-esd-converter", "TPCEsdConverter", "TPC-hltout-tracksegs-merger", "");
	216
	217	/////////////////////////////////////////////////////////////////////////////////////
	218	//
	219	// a kChain HLTOUT configuration for processing of {'TRACKS ':'TPC '} data blocks
	220	// produces an ESD object from the track structure
	221
	222	// publisher component
	223	handler->CreateConfiguration("TPC-hltout-tracks-publisher", "AliHLTOUTPublisher" , NULL, "");
	224
	225	// the esd converter configuration
	226	handler->CreateConfiguration("TPC-hltout-tracks-esd-converter", "TPCEsdConverter", "TPC-hltout-tracks-publisher", "");
5578cf60	227	}
	228	return 0;
	229	}
	230
dee38f1b	231	const char* AliHLTTPCAgent::GetReconstructionChains(AliRawReader* /rawReader/,
ff6648de	232	AliRunLoader* runloader) const
5578cf60	233	{
5578cf60	234	// see header file for class documentation
ff6648de	235	if (runloader) {
	236	// reconstruction chains for AliRoot simulation
	237	// Note: run loader is only available while running embedded into
	238	// AliRoot simulation
6595fdc6	239	//if (runloader->GetLoader("TPCLoader") != NULL)
	240	//return "TPC-esd-converter TPC-clusters";
	241	return "TPC-clusters";
ff6648de	242	}
ff6648de	243	return NULL;
5578cf60	244	}
	245
	246	const char* AliHLTTPCAgent::GetRequiredComponentLibraries() const
	247	{
	248	// see header file for class documentation
5ad1e184	249
	250	// actually, the TPC library has dependencies to Util and RCU
	251	// so the two has to be loaded anyhow before we get here
	252	//return "libAliHLTUtil.so libAliHLTRCU.so";
9c790772	253	return "libAliHLTUtil.so";
5578cf60	254	}
f3506ea2	255
e1440dab	256	int AliHLTTPCAgent::RegisterComponents(AliHLTComponentHandler* pHandler) const
f3506ea2	257	{
e1440dab	258	// see header file for class documentation
e1440dab	259	if (!pHandler) return -EINVAL;
672f8b8c	260
	261	pHandler->AddComponent(new AliHLTTPCRunStatisticsProducerComponent);
	262	pHandler->AddComponent(new AliHLTTPCEventStatisticsProducerComponent);
32e0c022	263	// pHandler->AddComponent(new AliHLTTPCCalibCEComponent);
	264	// pHandler->AddComponent(new AliHLTTPCCalibPulserComponent);
	265	// pHandler->AddComponent(new AliHLTTPCCalibPedestalComponent);
672f8b8c	266	pHandler->AddComponent(new AliHLTTPCCompModelInflaterComponent);
	267	pHandler->AddComponent(new AliHLTTPCCompModelDeflaterComponent);
	268	pHandler->AddComponent(new AliHLTTPCCompModelDeconverterComponent);
	269	pHandler->AddComponent(new AliHLTTPCCompModelConverterComponent);
	270	pHandler->AddComponent(new AliHLTTPCCompDumpComponent);
751d16ac	271	pHandler->AddComponent(new AliHLTTPCCAInputDataCompressorComponent);
400d857a	272	pHandler->AddComponent(new AliHLTTPCCATrackerComponent);
400d857a	273	pHandler->AddComponent(new AliHLTTPCCAGlobalMergerComponent);
a6fb14f5	274	pHandler->AddComponent(new AliHLTTPCTrackMCMarkerComponent);
672f8b8c	275	pHandler->AddComponent(new AliHLTTPCGlobalMergerComponent);
0973c527	276	pHandler->AddComponent(new AliHLTTPCdEdxComponent);
672f8b8c	277	pHandler->AddComponent(new AliHLTTPCSliceTrackerComponent);
	278	pHandler->AddComponent(new AliHLTTPCVertexFinderComponent);
	279	pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinderPacked));
	280	pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinderUnpacked));
	281	pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinderDecoder));
6f23f6cc	282	pHandler->AddComponent(new AliHLTTPCClusterFinderComponent(AliHLTTPCClusterFinderComponent::kClusterFinder32Bit));
672f8b8c	283	pHandler->AddComponent(new AliHLTTPCRawDataUnpackerComponent);
	284	pHandler->AddComponent(new AliHLTTPCDigitPublisherComponent);
	285	pHandler->AddComponent(new AliHLTTPCZeroSuppressionComponent);
e1440dab	286	pHandler->AddComponent(new AliHLTTPCDigitDumpComponent);
0efebbac	287	pHandler->AddComponent(new AliHLTTPCClusterDumpComponent);
f32b83e1	288	pHandler->AddComponent(new AliHLTTPCEsdWriterComponent::AliWriter);
f32b83e1	289	pHandler->AddComponent(new AliHLTTPCEsdWriterComponent::AliConverter);
1ac82ce6	290	pHandler->AddComponent(new AliHLTTPCOfflineClustererComponent);
1ac82ce6	291	pHandler->AddComponent(new AliHLTTPCOfflineTrackerComponent);
0e442a0a	292	pHandler->AddComponent(new AliHLTTPCOfflineTrackerCalibComponent);
e642ae99	293	pHandler->AddComponent(new AliHLTTPCOfflineCalibrationComponent);
0efebbac	294	pHandler->AddComponent(new AliHLTTPCClusterHistoComponent);
3a3550fa	295	pHandler->AddComponent(new AliHLTTPCNoiseMapComponent);
3a3550fa	296	pHandler->AddComponent(new AliHLTTPCHistogramHandlerComponent);
32e0c022	297	//pHandler->AddComponent(new AliHLTTPCCalibTracksComponent);
22f34176	298	pHandler->AddComponent(new AliHLTTPCTrackHistoComponent);
22f34176	299	pHandler->AddComponent(new AliHLTTPCTrackDumpComponent);
a912b63b	300	pHandler->AddComponent(new AliHLTTPCHWCFDataReverterComponent);
6af4f584	301	pHandler->AddComponent(new AliHLTTPCHWClusterTransformComponent);
32e0c022	302	// pHandler->AddComponent(new AliHLTTPCCalibSeedMakerComponent);
	303	// pHandler->AddComponent(new AliHLTTPCCalibTimeComponent);
	304	// pHandler->AddComponent(new AliHLTTPCCalibTimeGainComponent);
	305	// pHandler->AddComponent(new AliHLTTPCCalibrationComponent);
a912b63b	306
f3506ea2	307	return 0;
f3506ea2	308	}
626bfcc1	309
	310	int AliHLTTPCAgent::GetHandlerDescription(AliHLTComponentDataType dt,
	311	AliHLTUInt32_t spec,
	312	AliHLTOUTHandlerDesc& desc) const
	313	{
	314	// see header file for class documentation
73449074	315
73449074	316	// raw data blocks to be fed into offline reconstruction
626bfcc1	317	if (dt==(kAliHLTDataTypeDDLRaw\|kAliHLTDataOriginTPC)) {
	318	int slice=AliHLTTPCDefinitions::GetMinSliceNr(spec);
	319	int part=AliHLTTPCDefinitions::GetMinPatchNr(spec);
	320	if (slice==AliHLTTPCDefinitions::GetMaxSliceNr(spec) &&
	321	part==AliHLTTPCDefinitions::GetMaxPatchNr(spec)) {
	322	desc=AliHLTOUTHandlerDesc(kRawReader, dt, GetModuleId());
	323	return 1;
	324	} else {
	325	HLTWarning("handler can not process merged data from multiple ddls:"
	326	" min slice %d, max slice %d, min part %d, max part %d",
	327	slice, AliHLTTPCDefinitions::GetMaxSliceNr(spec),
	328	part, AliHLTTPCDefinitions::GetMaxPatchNr(spec));
	329	return 0;
	330	}
	331	}
73449074	332
84bb5f35	333	// dump for {'CLUSTERS':'TPC '} currently not used any more
	334	if (dt==AliHLTTPCDefinitions::fgkClustersDataType) {
	335	desc=AliHLTOUTHandlerDesc(kProprietary, dt, GetModuleId());
	336	return 1;
	337	}
	338
73449074	339	// afterburner for {'TRAKSEGS':'TPC '} blocks to be converted to ESD format
	340	if (dt==AliHLTTPCDefinitions::fgkTrackSegmentsDataType) {
	341	desc=AliHLTOUTHandlerDesc(kChain, dt, GetModuleId());
	342	return 1;
	343	}
	344
	345	// afterburner for {'TRACKS ':'TPC '} block to be converted to ESD format
	346	// there is only one data block
	347	if (dt==AliHLTTPCDefinitions::fgkTracksDataType) {
	348	desc=AliHLTOUTHandlerDesc(kChain, dt, GetModuleId());
	349	return 1;
	350	}
626bfcc1	351	return 0;
	352	}
	353
	354	AliHLTOUTHandler* AliHLTTPCAgent::GetOutputHandler(AliHLTComponentDataType dt,
b0914d2e	355	AliHLTUInt32_t /spec/)
626bfcc1	356	{
626bfcc1	357	// see header file for class documentation
73449074	358
73449074	359	// raw data blocks to be fed into offline reconstruction
626bfcc1	360	if (dt==(kAliHLTDataTypeDDLRaw\|kAliHLTDataOriginTPC)) {
	361	if (!fRawDataHandler) {
	362	fRawDataHandler=new AliHLTTPCAgent::AliHLTTPCRawDataHandler;
	363	}
626bfcc1	364	return fRawDataHandler;
626bfcc1	365	}
73449074	366
	367	// afterburner for {'TRAKSEGS':'TPC '} blocks to be converted to ESD format
	368	// in a kChain HLTOUT handler
	369	if (dt==AliHLTTPCDefinitions::fgkTrackSegmentsDataType) {
	370	if (fTracksegsDataHandler==NULL)
	371	fTracksegsDataHandler=new AliHLTOUTHandlerChain("chains=TPC-hltout-tracksegs-esd-converter");
	372	return fTracksegsDataHandler;
	373	}
	374
	375	// afterburner for {'TRACKS ':'TPC '} block to be converted to ESD format
	376	// there is only one data block
	377	if (dt==AliHLTTPCDefinitions::fgkTracksDataType) {
	378	return new AliHLTOUTHandlerChain("chains=TPC-hltout-tracks-esd-converter");
	379	}
	380
626bfcc1	381	return NULL;
	382	}
	383
	384	int AliHLTTPCAgent::DeleteOutputHandler(AliHLTOUTHandler* pInstance)
	385	{
	386	// see header file for class documentation
	387	if (pInstance==NULL) return -EINVAL;
	388
	389	if (pInstance==fRawDataHandler) {
73449074	390	delete fRawDataHandler;
	391	fRawDataHandler=NULL;
	392	}
	393
	394	if (pInstance==fTracksegsDataHandler) {
	395	delete fTracksegsDataHandler;
	396	fTracksegsDataHandler=NULL;
626bfcc1	397	}
	398	return 0;
	399	}
	400
	401	AliHLTTPCAgent::AliHLTTPCRawDataHandler::AliHLTTPCRawDataHandler()
	402	{
	403	// see header file for class documentation
	404	}
	405
	406	AliHLTTPCAgent::AliHLTTPCRawDataHandler::~AliHLTTPCRawDataHandler()
	407	{
	408	// see header file for class documentation
	409	}
	410
	411	int AliHLTTPCAgent::AliHLTTPCRawDataHandler::ProcessData(AliHLTOUT* pData)
	412	{
	413	// see header file for class documentation
	414	if (!pData) return -EINVAL;
	415	AliHLTComponentDataType dt=kAliHLTVoidDataType;
	416	AliHLTUInt32_t spec=kAliHLTVoidDataSpec;
	417	int iResult=pData->GetDataBlockDescription(dt, spec);
	418	if (iResult>=0) {
	419	int slice=AliHLTTPCDefinitions::GetMinSliceNr(spec);
	420	int part=AliHLTTPCDefinitions::GetMinPatchNr(spec);
	421	if (slice==AliHLTTPCDefinitions::GetMaxSliceNr(spec) &&
	422	part==AliHLTTPCDefinitions::GetMaxPatchNr(spec)) {
	423	iResult=768;
	424	if (part>1) iResult+=72+4*slice+(part-2);
	425	else iResult+=2*slice+part;
	426	} else {
	427	HLTError("handler can not process merged data from multiple ddls:"
	428	" min slice %d, max slice %d, min part %d, max part %d",
	429	slice, AliHLTTPCDefinitions::GetMaxSliceNr(spec),
	430	part, AliHLTTPCDefinitions::GetMaxPatchNr(spec));
	431	iResult=-EBADMSG;
	432	}
	433	}
	434	return iResult;
	435	}