#include "AliTPCcalibDB.h"
#include "AliTPCCalPad.h"
#include "AliTPCParam.h"
+//#include "AliHLTTPCCAInputDataCompressorComponent.h"
+//#include "AliHLTTPCCADef.h"
#include <cstdlib>
#include <cerrno>
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
#include "AliCDBStorage.h"
+#include "TGeoGlobalMagField.h"
+#include "AliGeomManager.h"
#include <sys/time.h>
/** ROOT macro for the implementation of ROOT specific class methods */
ClassImp(AliHLTTPCClusterFinderComponent)
+const char* AliHLTTPCClusterFinderComponent::fgkOCDBEntryPacked="HLT/ConfigTPC/TPCClusterFinderPacked";
+const char* AliHLTTPCClusterFinderComponent::fgkOCDBEntryUnpacked="HLT/ConfigTPC/TPCClusterFinderUnpacked";
+const char* AliHLTTPCClusterFinderComponent::fgkOCDBEntryDecoder="HLT/ConfigTPC/TPCClusterFinderDecoder";
+const char* AliHLTTPCClusterFinderComponent::fgkOCDBEntry32Bit="HLT/ConfigTPC/TPCClusterFinder32Bit";
+
AliHLTTPCClusterFinderComponent::AliHLTTPCClusterFinderComponent(int mode)
:
fClusterFinder(NULL),
fGetActivePads(0),
fFirstTimeBin(-1),
fLastTimeBin(-1),
- fDoMC(kFALSE)
+ fDoMC(kFALSE),
+ fReleaseMemory( kFALSE ),
+ fBenchmark("TPCClusterFinder")
{
// see header file for class documentation
// or
tgtList.clear();
tgtList.push_back(AliHLTTPCDefinitions::fgkClustersDataType);
tgtList.push_back(kAliHLTDataTypeHwAddr16);
+ //tgtList.push_back(AliHLTTPCCADefinitions::fgkCompressedInputDataType);
return tgtList.size();
}
{
// see header file for class documentation
if ( fClusterFinder )
- return EINPROGRESS;
+ return -EINPROGRESS;
//Test if the OCDB entries used by AliTPCTransform is availible
AliTPCcalibDB* calib=AliTPCcalibDB::Instance();
//
+ if(!calib){
+ HLTError("AliTPCcalibDB does not exist");
+ return -ENOENT;
+ }
+ calib->SetRun(GetRunNo());
+ calib->UpdateRunInformations(GetRunNo());
AliTPCCalPad * time0TPC = calib->GetPadTime0();
if(!time0TPC){
HLTError("OCDB entry TPC/Calib/PadTime0 (AliTPCcalibDB::GetPadTime0()) is not available.");
return -ENOENT;
}
+ // Check field
+ if (!TGeoGlobalMagField::Instance()) {
+ HLTError("magnetic field not initialized, please set up TGeoGlobalMagField and AliMagF");
+ return -ENODEV;
+ }
+ calib->SetExBField(GetBz());
+
+ if(!AliGeomManager::GetGeometry()){
+ AliGeomManager::LoadGeometry();
+ }
fClusterFinder = new AliHLTTPCClusterFinder();
+ // first configure the default
+ int iResult = 0;
+ switch(fModeSwitch){
+ case kClusterFinderPacked:
+ iResult=ConfigureFromCDBTObjString(fgkOCDBEntryPacked);
+ break;
+ case kClusterFinderUnpacked:
+ iResult=ConfigureFromCDBTObjString(fgkOCDBEntryUnpacked);
+ break;
+ case kClusterFinderDecoder:
+ iResult=ConfigureFromCDBTObjString(fgkOCDBEntryDecoder);
+ break;
+ case kClusterFinder32Bit:
+ iResult=ConfigureFromCDBTObjString(fgkOCDBEntry32Bit);
+ break;
+ }
+
+ // configure from the command line parameters if specified
+ if (iResult>=0 && argc>0)
+ iResult=ConfigureFromArgumentString(argc, argv);
+ // return iResult;
+
+ /*
Int_t iResult=0;
TString configuration="";
TString argument="";
} else {
iResult=Reconfigure(NULL, NULL);
}
+ */
//Checking for conflicting arguments
if(fClusterDeconv){
if(fLastTimeBin>0 && fLastTimeBin>fFirstTimeBin && fLastTimeBin<=AliHLTTPCTransform::GetNTimeBins()){
fClusterFinder->SetLastTimeBin(fLastTimeBin);
}
+ fBenchmark.Reset();
+ fBenchmark.SetTimer(0,"total");
+ fBenchmark.SetTimer(1,"reco");
- return 0;
+ return iResult;
}
int AliHLTTPCClusterFinderComponent::DoDeinit()
{
// see header file for class documentation
int iResult=0;
+ //SG!!!
+ //fDeconvTime = kTRUE;
+ //fClusterFinder->SetDeconvTime(fDeconvTime);
+ // fDeconvPad = kTRUE;
+ //fClusterFinder->SetDeconvPad(fDeconvPad);
+
if(fReader == NULL){
HLTFatal("Digit reader not initialized, skipping HLT TPC cluster reconstruction.");
size=0;
- return 0;
+ return -ENODEV;
}
- if(GetFirstInputBlock( kAliHLTDataTypeSOR ) || GetFirstInputBlock( kAliHLTDataTypeEOR )){
+ if(!IsDataEvent()){
size=0;
return 0;
}
+ fBenchmark.StartNewEvent();
+ fBenchmark.Start(0);
+
// == init iter (pointer to datablock)
const AliHLTComponentBlockData* iter = NULL;
unsigned long ndx;
}
+ fBenchmark.AddInput(iter->fSize);
+
slice = AliHLTTPCDefinitions::GetMinSliceNr( *iter );
patch = AliHLTTPCDefinitions::GetMinPatchNr( *iter );
fClusterFinder->InitSlice( slice, patch, maxPoints );
fClusterFinder->SetOutputArray( (AliHLTTPCSpacePointData*)outPtr->fSpacePoints );
+ fBenchmark.Start(1);
+
if(fUnsorted){
if(fGetActivePads){
fClusterFinder->SetDoPadSelection(kTRUE);
fClusterFinder->Read(iter->fPtr, iter->fSize );
fClusterFinder->ProcessDigits();
}
+ fBenchmark.Stop(1);
+
fReader->Reset();
realPoints = fClusterFinder->GetNumberOfClusters();
bd.fSpecification = iter->fSpecification;
bd.fDataType = AliHLTTPCDefinitions::fgkClustersDataType;
outputBlocks.push_back( bd );
-
+
+ fBenchmark.AddOutput(bd.fSize);
+
tSize += mysize;
outBPtr += mysize;
outPtr = (AliHLTTPCClusterData*)outBPtr;
bdHW.fDataType = kAliHLTDataTypeHwAddr16;
outputBlocks.push_back( bdHW );
+ fBenchmark.AddOutput(bdHW.fSize);
tSize+=nHWAdd*sizeof(AliHLTUInt16_t);
}
+ /*
+ { // compressed output for the CA tracker
+
+ AliHLTUInt32_t dSize = 0;
+
+ int ret = AliHLTTPCCAInputDataCompressorComponent::Compress( (AliHLTTPCClusterData*)( outputPtr + bd.fOffset ),
+ size - tSize,
+ outputPtr+tSize,
+ dSize );
+
+ if ( ret!=0 || tSize + dSize > size )
+ {
+ Logging( kHLTLogFatal, "HLT::TPCClusterFinder::DoEvent", "Too much data",
+ "Data written over allowed buffer. Amount written: %lu, allowed amount: %lu.",
+ tSize + dSize, size );
+ iResult=-ENOSPC;
+ break;
+ }
+
+ AliHLTComponentBlockData bdCompressed;
+ FillBlockData( bdCompressed );
+ bdCompressed.fOffset = tSize ;
+ bdCompressed.fSize = dSize;
+ bdCompressed.fSpecification = iter->fSpecification;
+ bdCompressed.fDataType = AliHLTTPCCADefinitions::fgkCompressedInputDataType;
+ outputBlocks.push_back( bdCompressed );
+
+ tSize += dSize;
+ outBPtr += dSize;
+ outPtr = (AliHLTTPCClusterData*)outBPtr;
+ }
+ */
+
if(fDoMC){
Int_t maxNumberOfClusterMCInfo = (Int_t)((size-tSize)/sizeof(AliHLTTPCClusterFinder::ClusterMCInfo)-1);
AliHLTTPCClusterFinder::ClusterMCInfo* outputMCInfo= (AliHLTTPCClusterFinder::ClusterMCInfo*)(outputPtr+tSize);
bdMCInfo.fSpecification = iter->fSpecification;
bdMCInfo.fDataType = AliHLTTPCDefinitions::fgkAliHLTDataTypeClusterMCInfo;
outputBlocks.push_back( bdMCInfo );
-
+ fBenchmark.AddOutput(bdMCInfo.fSize);
+
tSize+=nMCInfo*sizeof(AliHLTTPCClusterFinder::ClusterMCInfo);
}
if (iResult>=0)
size = tSize;
+ fBenchmark.Stop(0);
+ HLTInfo(fBenchmark.GetStatistics());
return iResult;
}
+int AliHLTTPCClusterFinderComponent::ScanConfigurationArgument(int argc, const char** argv){
+
+ // see header file for class documentation
+
+ if (argc<=0) return 0;
+ int i=0;
+ TString argument=argv[i];
+
+ if (argument.CompareTo("-solenoidBz")==0){
+ if (++i>=argc) return -EPROTO;
+ HLTWarning("argument -solenoidBz is deprecated, magnetic field set up globally (%f)", GetBz());
+ return 2;
+ }
+
+ if (argument.CompareTo("-update-calibdb")==0 || argument.CompareTo("-update-transform")==0 ){
+ if(fClusterFinder->UpdateCalibDB()){
+ HLTDebug("CalibDB and offline transform successfully updated.");
+ }
+ else{
+ HLTError("CalibDB could not be updated.");
+ }
+ return 1;
+ }
+
+ if (argument.CompareTo("-deconvolute-time")==0){
+ HLTDebug("Switching on deconvolution in time direction.");
+ fDeconvTime = kTRUE;
+ fClusterFinder->SetDeconvTime(fDeconvTime);
+ return 1;
+ }
+
+ if (argument.CompareTo("-deconvolute-pad")==0){
+ HLTDebug("Switching on deconvolution in pad direction.");
+ fDeconvPad = kTRUE;
+ fClusterFinder->SetDeconvPad(fDeconvPad);
+ return 1;
+ }
+
+ if (argument.CompareTo("-timebins")==0 || argument.CompareTo("timebins" )==0){
+ HLTWarning("Argument %s is depreciated after moving to the offline AliTPCTransform class for xyz calculations.",argument.Data());
+ /*
+ if (++i>=argc) return -EPROTO;
+ argument=argv[i];
+ AliHLTTPCTransform::SetNTimeBins(argument.Atoi());
+ fClusterFinder->UpdateLastTimeBin();
+ HLTInfo("number of timebins set to %d, zbin=%f", AliHLTTPCTransform::GetNTimeBins(), AliHLTTPCTransform::GetZWidth());
+ return 2;
+ */
+ if(argument.CompareTo("timebins")==0){
+ HLTWarning("Argument 'timebins' is old, please switch to new argument naming convention (-timebins). The timebins argument will still work, but please change anyway.");
+ }
+ return 2;
+ }
+
+ if (argument.CompareTo("-first-timebin")==0){
+ if (++i>=argc) return -EPROTO;
+ argument=argv[i];
+ fFirstTimeBin = argument.Atoi();
+ if(fFirstTimeBin>=0){
+ HLTDebug("fFirstTimeBin set to %d",fFirstTimeBin);
+ fClusterFinder->SetFirstTimeBin(fFirstTimeBin);
+ }
+ else{
+ HLTError("-first-timebin specifier is negative: %d",fFirstTimeBin);
+ }
+ return 2;
+ }
+
+ if (argument.CompareTo("-last-timebin")==0){
+ if (++i>=argc) return -EPROTO;
+ argument=argv[i];
+ fLastTimeBin = argument.Atoi();
+ if(fLastTimeBin<AliHLTTPCTransform::GetNTimeBins()){
+ HLTDebug("fLastTimeBin set to %d",fLastTimeBin);
+ }
+ else{
+ HLTError("fLastTimeBins is too big: %d. Maximum: %d",fLastTimeBin,AliHLTTPCTransform::GetNTimeBins());
+ }
+ return 2;
+ }
+
+ if (argument.CompareTo("-sorted")==0) {
+ fUnsorted=0;
+ HLTDebug("Swithching unsorted off.");
+ fClusterFinder->SetUnsorted(0);
+ return 1;
+ }
+
+ if (argument.CompareTo("-do-mc")==0) {
+ fDoMC=kTRUE;
+ fClusterFinder->SetDoMC(fDoMC);
+ HLTDebug("Setting fDoMC to true.");
+ return 1;
+ }
+ if (argument.CompareTo("-release-memory")==0) {
+ fReleaseMemory=kTRUE;
+ fClusterFinder->SetReleaseMemory( fReleaseMemory );
+ HLTDebug("Setting fReleaseMemory to true.");
+ return 1;
+ }
+
+ if (argument.CompareTo("-active-pads")==0 || argument.CompareTo("activepads")==0){
+ if(argument.CompareTo("activepads" )==0){
+ HLTWarning("Please change to new component argument naming scheme and use '-active-pads' instead of 'activepads'");
+ }
+ HLTDebug("Switching on ActivePads");
+ fGetActivePads = 1;
+ fClusterFinder->SetDoPadSelection(kTRUE);
+ return 1;
+ }
+
+ if (argument.CompareTo("-occupancy-limit")==0 || argument.CompareTo("occupancy-limit")==0){
+ if(argument.CompareTo("occupancy-limit" )==0){
+ HLTWarning("Please switch to new component argument naming convention, use '-occupancy-limit' instead of 'occupancy-limit'");
+ }
+ if (++i>=argc) return -EPROTO;
+ argument=argv[i];
+ fClusterFinder->SetOccupancyLimit(argument.Atof());
+ HLTDebug("Occupancy limit set to occulimit %f", argument.Atof());
+ return 2;
+ }
+
+ if (argument.CompareTo("rawreadermode")==0){
+ if (++i>=argc) return -EPROTO;
+ HLTWarning("Argument 'rawreadermode' is deprecated");
+ return 2;
+ }
+
+ if (argument.CompareTo("pp-run")==0){
+ HLTWarning("Argument 'pp-run' is obsolete, deconvolution is swiched off in both time and pad directions by default.");
+ fClusterDeconv = false;
+ return 1;
+ }
+
+ if (argument.CompareTo("adc-threshold" )==0){
+ if (++i>=argc) return -EPROTO;
+ HLTWarning("'adc-threshold' is no longer a valid argument, please use TPCZeroSuppression component if you want to zerosuppress data.");
+ return 2;
+ }
+
+ if (argument.CompareTo("oldrcuformat" )==0){
+ if (++i>=argc) return -EPROTO;
+ HLTWarning("Argument 'oldrcuformat' is deprecated.");
+ return 2;
+ }
+
+ if (argument.CompareTo("unsorted" )==0 || argument.CompareTo("-unsorted" )==0){
+ HLTWarning("Argument is obsolete, unsorted reading is default.");
+ // fClusterFinder->SetUnsorted(1);
+ return 1;
+ }
+ if (argument.CompareTo("nsigma-threshold")==0){
+ if (++i>=argc) return -EPROTO;
+ HLTWarning("Argument 'nsigma-threshold' argument is obsolete.");
+ return 2;
+ }
+
+ // unknown argument
+ return -EINVAL;
+}
+
+int AliHLTTPCClusterFinderComponent::Reconfigure(const char* cdbEntry, const char* /*chainId*/)
+{
+ // see header file for class documentation
+
+ const char* entry=cdbEntry;
+ if (!entry || entry[0]==0){
+ switch(fModeSwitch){
+ case kClusterFinderPacked:
+ entry=fgkOCDBEntryPacked;
+ break;
+ case kClusterFinderUnpacked:
+ entry=fgkOCDBEntryUnpacked;
+ break;
+ case kClusterFinderDecoder:
+ entry=fgkOCDBEntryDecoder;
+ break;
+ case kClusterFinder32Bit:
+ entry=fgkOCDBEntry32Bit;
+ break;
+ }
+ }
+
+ return ConfigureFromCDBTObjString(entry);
+
+ /*
+ int iResult=0;
+
+ const char* path="HLT/ConfigTPC/ClusterFinderComponent";
+ if (cdbEntry) path=cdbEntry;
+ if (path) {
+ HLTInfo("reconfigure from entry %s, chain id %s", path, (chainId!=NULL && chainId[0]!=0)?chainId:"<none>");
+ AliCDBEntry *pEntry = AliCDBManager::Instance()->Get(path);//,GetRunNo());
+ if (pEntry) {
+ TObjString* pString=dynamic_cast<TObjString*>(pEntry->GetObject());
+ if (pString) {
+ HLTInfo("received configuration object: %s", pString->GetString().Data());
+ iResult = Configure(pString->GetString().Data());
+ } else {
+ HLTError("configuration object \"%s\" has wrong type, required TObjString", path);
+ }
+ } else {
+ HLTError("can not fetch object \"%s\" from CDB", path);
+ }
+ }
+ return iResult;
+ */
+}
+
int AliHLTTPCClusterFinderComponent::Configure(const char* arguments){
// see header file for class documentation
int iResult=0;
fClusterFinder->SetDeconvPad(fDeconvPad);
}
else if (argument.CompareTo("-timebins")==0 || argument.CompareTo("timebins" )==0){
+ HLTWarning("Argument %s is depreciated after moving to the offline AliTPCTransform class for xyz calculations.",argument.Data());
+ /*
if ((bMissingParam=(++i>=pTokens->GetEntries()))) break;
AliHLTTPCTransform::SetNTimeBins(((TObjString*)pTokens->At(i))->GetString().Atoi());
fClusterFinder->UpdateLastTimeBin();
HLTInfo("number of timebins set to %d, zbin=%f", AliHLTTPCTransform::GetNTimeBins(), AliHLTTPCTransform::GetZWidth());
+ */
if(argument.CompareTo("timebins")==0){
HLTWarning("Argument 'timebins' is old, please switch to new argument naming convention (-timebins). The timebins argument will still work, but please change anyway.");
}
+
}
else if (argument.CompareTo("-first-timebin")==0){
if ((bMissingParam=(++i>=pTokens->GetEntries()))) break;
fClusterFinder->SetDoMC(fDoMC);
HLTInfo("Setting fDoMC to true.");
}
+ else if (argument.CompareTo("-release-memory")==0) {
+ fReleaseMemory = kTRUE;
+ fClusterFinder->SetReleaseMemory( kTRUE );
+ HLTInfo("Setting fReleaseMemory to true.");
+ }
else if (argument.CompareTo("-active-pads")==0 || argument.CompareTo("activepads")==0){
if(argument.CompareTo("activepads" )==0){
HLTWarning("Please change to new component argument naming scheme and use '-active-pads' instead of 'activepads'");
if ((bMissingParam=(++i>=pTokens->GetEntries()))) break;
HLTWarning("Argument 'nsigma-threshold' argument is obsolete.");
}
+ else if (argument.CompareTo("-update-calibdb")==0){
+ fClusterFinder->UpdateCalibDB();
+ }
else {
HLTError("unknown argument %s", argument.Data());
iResult=-EINVAL;
return iResult;
}
-int AliHLTTPCClusterFinderComponent::Reconfigure(const char* cdbEntry, const char* chainId)
-{
+void AliHLTTPCClusterFinderComponent::GetOCDBObjectDescription( TMap* const targetMap){
+// Get a list of OCDB object description needed for the particular component
- int iResult=0;
- // see header file for class documentation
- const char* path="HLT/ConfigTPC/ClusterFinderComponent";
- if (cdbEntry) path=cdbEntry;
- if (path) {
- HLTInfo("reconfigure from entry %s, chain id %s", path, (chainId!=NULL && chainId[0]!=0)?chainId:"<none>");
- AliCDBEntry *pEntry = AliCDBManager::Instance()->Get(path/*,GetRunNo()*/);
- if (pEntry) {
- TObjString* pString=dynamic_cast<TObjString*>(pEntry->GetObject());
- if (pString) {
- HLTInfo("received configuration object: %s", pString->GetString().Data());
- iResult = Configure(pString->GetString().Data());
- } else {
- HLTError("configuration object \"%s\" has wrong type, required TObjString", path);
- }
- } else {
- HLTError("can not fetch object \"%s\" from CDB", path);
- }
- }
- return iResult;
+ if (!targetMap) return;
+
+ // OCDB entries for component arguments
+ targetMap->Add(new TObjString("HLT/ConfigTPC/TPCClusterFinder32Bit"), new TObjString("component arguments, empty at the moment"));
+ targetMap->Add(new TObjString("HLT/ConfigTPC/TPCClusterFinderPacked"), new TObjString("component arguments, empty at the moment"));
+ targetMap->Add(new TObjString("HLT/ConfigTPC/TPCClusterFinderUnpacked"), new TObjString("component arguments, empty at the moment"));
+ targetMap->Add(new TObjString("HLT/ConfigTPC/TPCClusterFinderDecoder"), new TObjString("component arguments, empty at the moment"));
+
+ // OCDB entries to be fetched by the TAXI (access via the AliTPCcalibDB class)
+ targetMap->Add(new TObjString("TPC/Calib/Parameters"), new TObjString("unknown content"));
+ targetMap->Add(new TObjString("TPC/Calib/TimeDrift"), new TObjString("drift velocity calibration"));
+ targetMap->Add(new TObjString("TPC/Calib/Temperature"), new TObjString("temperature map"));
+ targetMap->Add(new TObjString("TPC/Calib/PadGainFactor"), new TObjString("gain factor pad by pad"));
+ targetMap->Add(new TObjString("TPC/Calib/ClusterParam"), new TObjString("cluster parameters"));
+
+ // OCDB entries needed to be fetched by the Pendolino
+ targetMap->Add(new TObjString("TPC/Calib/AltroConfig"), new TObjString("contains the altro config, e.g. info about the L0 trigger timing"));
+ targetMap->Add(new TObjString("GRP/CTP/CTPtiming"), new TObjString("content used in the cluster coordinate transformation in relation to the L0 trigger timing"));
+
+ // OCDB entries necessary for replaying data on the HLT cluster
+ targetMap->Add(new TObjString("GRP/GRP/Data"), new TObjString("contains magnetic field info"));
+
+ // OCDB entries needed to suppress fatals/errors/warnings during reconstruction
+ targetMap->Add(new TObjString("TPC/Calib/PadTime0"), new TObjString("time0 offset pad by pad"));
+ targetMap->Add(new TObjString("TPC/Calib/PadNoise"), new TObjString("pad noise values"));
+ targetMap->Add(new TObjString("TPC/Calib/Pedestals"), new TObjString("pedestal info"));
+ targetMap->Add(new TObjString("TPC/Calib/Pulser"), new TObjString("pulser info"));
+ targetMap->Add(new TObjString("TPC/Calib/CE"), new TObjString("CE laser calibration result"));
+ targetMap->Add(new TObjString("TPC/Calib/Raw"), new TObjString("unknown content"));
+ targetMap->Add(new TObjString("TPC/Calib/QA"), new TObjString("not important"));
+ targetMap->Add(new TObjString("TPC/Calib/Mapping"), new TObjString("unknown content"));
+ targetMap->Add(new TObjString("TPC/Calib/Goofie"), new TObjString("Goofie values, not used at the moment (05.03.2010)"));
+ targetMap->Add(new TObjString("TPC/Calib/HighVoltage"), new TObjString("high voltage values, not used"));
+ targetMap->Add(new TObjString("TPC/Calib/Ref"), new TObjString("unknown content"));
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff