AltModuleSliceList.cxx

Go to the documentation of this file.

// $Id: AltModuleSliceList.cxx,v 1.16 2004/07/26 12:22:52 west Exp $
//
// AltModuleSliceList.cxx
//
// A JobControl Module for reconstructing CandSliceList 
//
// Costas Andreopoulos <C.V.Andreopoulos@rl.ac.uk>
// CCLRC, Rutherford Appleton Laboratory
// July 01, 2003

#include <cassert>

#include "AltModuleSliceList.h"

#include "RecoBase/CandSliceListHandle.h"
#include "RecoBase/CandSliceList.h"
#include "RecoBase/CandSliceHandle.h"

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "CandData/CandHeader.h"
#include "CandData/CandRecord.h"
#include "Candidate/CandContext.h"
#include "Conventions/Munits.h"
#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "RecoBase/CandStripListHandle.h"
#include "Validity/VldContext.h"

ClassImp(AltModuleSliceList)

//____________________________________________________________________________
CVSID("$Id: AltModuleSliceList.cxx,v 1.16 2004/07/26 12:22:52 west Exp $");
//____________________________________________________________________________
JOBMODULE(AltModuleSliceList, "AltModuleSliceList", "CandSliceList Reco.");
//____________________________________________________________________________
AltModuleSliceList::AltModuleSliceList() : 
fListIn("CandStripListHandle"), 
fListOut("CandSliceList")
{
  MSG("Alt", Msg::kVerbose) << "AltModuleSliceList::Constructor\n";

  //-- Get AlgFactory Instance.

  AlgFactory &af = AlgFactory::GetInstance();

  //-- Register Algorithm configuration 

  af.Register("AltAlgSlice",     "default", "libAltReco.so", "AlgConfig");
  af.Register("AltAlgSliceList", "default", "libAltReco.so", "AlgConfig");            
}

//____________________________________________________________________________
AltModuleSliceList::~AltModuleSliceList() 
{
  MSG("Alt", Msg::kVerbose) << "AltModuleSliceList::Destructor\n";
}
//____________________________________________________________________________
void AltModuleSliceList::BeginJob() 
{
  MSG("Alt", Msg::kVerbose) << "AltModuleSliceList::BeginJob\n";
}
//____________________________________________________________________________
const Registry & AltModuleSliceList::DefaultConfig() const
{
  MSG("Alt", Msg::kDebug) << "AltModuleSliceList::DefaultConfig" << endl;
    
  static Registry r;
  
  std::string name = this->JobCModule::GetName();
  name += ".config.default";
  r.SetName(name.c_str());
  r.UnLockValues();
  r.LockValues();

  //-- Get AlgFactory Instance.
  AlgFactory &af = AlgFactory::GetInstance();

  //-- Get Algorithm handle and retrieve its config  
  AlgHandle alglisthandle = af.GetAlgHandle("AltAlgSliceList","default");
  AlgConfig &alglistconfig = alglisthandle.GetAlgConfig();

  //-- Set default configuration parameters
  alglistconfig.Set("Grfx_DebugGraphics",                  0);
  alglistconfig.Set("Grfx_TimeProfileLogView",             1);  
  alglistconfig.Set("Pkf_PeakThreshold",                   40);
  alglistconfig.Set("Pkf_NSuccessiveEmptyBins",            1);
  alglistconfig.Set("Pkf_MuSpecPeakThreshold",             20);
  alglistconfig.Set("Pkf_MuSpecNSuccessiveEmptyBins",      2);
  alglistconfig.Set("Pkf_LowPeakThreshold",                30);
  alglistconfig.Set("Pkf_LowNSuccessiveEmptyBins",         2);
  alglistconfig.Set("Pkf_NofMergedTimeBins",               1);
  alglistconfig.Set("Pkf_WeightProfileWithCharge",         1);
  alglistconfig.Set("Pkf_RecursivePeakSearch",             0);
  alglistconfig.Set("Pkf_TimeWindowBefPeak",               2);
  alglistconfig.Set("Pkf_TimeWindowAftPeak",               3);
  alglistconfig.Set("Refinement_Dissolving",               1);
  alglistconfig.Set("Refinement_Merging",                  0);
  alglistconfig.Set("Refinement_KMeansClustering",         1);
  alglistconfig.Set("Refinement_MSTClustering",            0);
  alglistconfig.Set("OrphanStrips_TimeWindow",             60*Munits::ns);
  alglistconfig.Set("OrphanStrips_PlaneWindow",            10);
  alglistconfig.Set("OrphanStrips_QWeight",                1);
  alglistconfig.Set("KMeans_TimeWindow",                   40.*Munits::ns);
  alglistconfig.Set("KMeans_PlaneWindow",                  3);
  alglistconfig.Set("KMeans_TPosWindow",                   30.*Munits::cm);
  alglistconfig.Set("MuSpec_TimeAftUpstrActivity",         80.*Munits::ns);
  alglistconfig.Set("MuSpec_TimeBefUpstrActivity",         40.*Munits::ns);
  alglistconfig.Set("MuSpec_NUpstrPlanes",                 30);
  alglistconfig.Set("MuSpec_HitStripsBefSpectr",           1);
  alglistconfig.Set("MuSpec_NPlnBefSpectr",                15);   
  alglistconfig.Set("MuSpec_SuppressUnmatchedSlices",      1);   
  alglistconfig.Set("MinHitStripsInSlice",                 4);
  alglistconfig.Set("MinChargeInSlice",                    5.);
  alglistconfig.Set("TimeResolution",                      19.*Munits::ns);
  alglistconfig.Set("TimeDiffBetweenPeaks",                70.*Munits::ns);
  alglistconfig.Set("ZDiffBetweenPeaks",                   200.*Munits::cm);
  alglistconfig.Set("ZDiffBetweenEnds",                    40.*Munits::cm);
  alglistconfig.Set("UVDiffBetweenPeaks",                  40.*Munits::cm);
        
  //alglistconfig.Print();

  return r;
}
//____________________________________________________________________________
void AltModuleSliceList::Config(const Registry &r)
{
  MSG("Alt", Msg::kDebug) << "AltModuleSliceList::Config" << endl;
    
  //-- Get AlgFactory Instance.
  AlgFactory &af = AlgFactory::GetInstance();

  //-- Get Algorithm handle and retrieve its config.    
  AlgHandle alglisthandle = af.GetAlgHandle("AltAlgSliceList","default");
  AlgConfig &alglistconfig = alglisthandle.GetAlgConfig();
    
  //-- Set configuration parameters    

  int          tmpi = 0;
  const char * tmpc = 0;
  double       tmpd = 0;

  if (r.Get("ListIn",tmpc))     fListIn  = tmpc;
  if (r.Get("ListOut",tmpc))    fListOut = tmpc;
  
  alglistconfig.UnLockValues();

  if (r.Get("Grfx_DebugGraphics",tmpi)) 
                                alglistconfig.Set("Grfx_DebugGraphics", tmpi);
  if (r.Get("Grfx_TimeProfileLogView",tmpi)) 
                           alglistconfig.Set("Grfx_TimeProfileLogView", tmpi);
  if (r.Get("Pkf_PeakThreshold",tmpi)) 
                                 alglistconfig.Set("Pkf_PeakThreshold", tmpi);
  if (r.Get("Pkf_NSuccessiveEmptyBins",tmpi)) 
                          alglistconfig.Set("Pkf_NSuccessiveEmptyBins", tmpi);
  if (r.Get("Pkf_MuSpecPeakThreshold",tmpi)) 
                           alglistconfig.Set("Pkf_MuSpecPeakThreshold", tmpi);
  if (r.Get("Pkf_MuSpecNSuccessiveEmptyBins",tmpi)) 
                    alglistconfig.Set("Pkf_MuSpecNSuccessiveEmptyBins", tmpi);
  if (r.Get("Pkf_LowPeakThreshold",tmpi)) 
                              alglistconfig.Set("Pkf_LowPeakThreshold", tmpi);
  if (r.Get("Pkf_LowNSuccessiveEmptyBins",tmpi)) 
                       alglistconfig.Set("Pkf_LowNSuccessiveEmptyBins", tmpi);
  if (r.Get("Pkf_NofMergedTimeBins",tmpi)) 
                            alglistconfig.Set("NumberOfMergedTimeBins", tmpi);
  if (r.Get("Pkf_WeightProfileWithCharge",tmpi)) 
                       alglistconfig.Set("Pkf_WeightProfileWithCharge", tmpi);
  if (r.Get("Pkf_RecursivePeakSearch",tmpi)) 
                           alglistconfig.Set("Pkf_RecursivePeakSearch", tmpi);
  if (r.Get("Pkf_TimeWindowBefPeak",tmpi)) 
                             alglistconfig.Set("Pkf_TimeWindowBefPeak", tmpi);
  if (r.Get("Pkf_TimeWindowAftPeak",tmpi)) 
                             alglistconfig.Set("Pkf_TimeWindowAftPeak", tmpi);
  if (r.Get("Refinement_Dissolving",tmpi)) 
                             alglistconfig.Set("Refinement_Dissolving", tmpi);
  if (r.Get("Refinement_Merging",tmpi)) 
                                alglistconfig.Set("Refinement_Merging", tmpi);
  if (r.Get("Refinement_KMeansClustering",tmpi)) 
                       alglistconfig.Set("Refinement_KMeansClustering", tmpi);
  if (r.Get("Refinement_MSTClustering",tmpi)) 
                          alglistconfig.Set("Refinement_MSTClustering", tmpi);
  if (r.Get("OrphanStrips_QWeight",tmpi)) 
                              alglistconfig.Set("OrphanStrips_QWeight", tmpi);
  if (r.Get("OrphanStrips_PlaneWindow",tmpi)) 
                          alglistconfig.Set("OrphanStrips_PlaneWindow", tmpi);
  if (r.Get("OrphanStrips_TimeWindow",tmpd)) 
                           alglistconfig.Set("OrphanStrips_TimeWindow", tmpd);
  if (r.Get("KMeans_TimeWindow",tmpd)) 
                                 alglistconfig.Set("KMeans_TimeWindow", tmpd);
  if (r.Get("KMeans_PlaneWindow",tmpi)) 
                                alglistconfig.Set("KMeans_PlaneWindow", tmpi);
  if (r.Get("KMeans_TPosWindow",tmpd)) 
                                 alglistconfig.Set("KMeans_TPosWindow", tmpd);
  if (r.Get("MuSpec_TimeAftUpstrActivity",tmpd)) 
                       alglistconfig.Set("MuSpec_TimeAftUpstrActivity", tmpd);
  if (r.Get("MuSpec_TimeBefUpstrActivity",tmpd)) 
                       alglistconfig.Set("MuSpec_TimeBefUpstrActivity", tmpd);
  if (r.Get("MuSpec_NUpstrPlanes",tmpi)) 
                               alglistconfig.Set("MuSpec_NUpstrPlanes", tmpi);
  if (r.Get("MuSpec_HitStripsBefSpectr",tmpi)) 
                         alglistconfig.Set("MuSpec_HitStripsBefSpectr", tmpi);
  if (r.Get("MuSpec_NPlnBefSpectr",tmpi)) 
                              alglistconfig.Set("MuSpec_NPlnBefSpectr", tmpi);
  if (r.Get("MuSpec_SuppressUnmatchedSlices",tmpi)) 
                    alglistconfig.Set("MuSpec_SuppressUnmatchedSlices", tmpi);
  if (r.Get("TimeResolution",tmpd)) 
                                    alglistconfig.Set("TimeResolution", tmpd);
  if (r.Get("MinHitStripsInSlice",tmpi)) 
                               alglistconfig.Set("MinHitStripsInSlice", tmpi);
  if (r.Get("TimeDiffBetweenPeaks",tmpd)) 
                              alglistconfig.Set("TimeDiffBetweenPeaks", tmpd);
  if (r.Get("ZDiffBetweenPeaks",tmpd)) 
                                 alglistconfig.Set("ZDiffBetweenPeaks", tmpd);
  if (r.Get("ZDiffBetweenEnds",tmpd))
                                  alglistconfig.Set("ZDiffBetweenEnds", tmpd);
  if (r.Get("UVDiffBetweenPeaks",tmpd))
                                alglistconfig.Set("UVDiffBetweenPeaks", tmpd);
  if (r.Get("MinChargeInSlice",tmpd)) 
                                  alglistconfig.Set("MinChargeInSlice", tmpd);

  alglistconfig.LockValues();
}
//____________________________________________________________________________
JobCResult AltModuleSliceList::Reco(MomNavigator *mom)
{
  MSG("Alt", Msg::kVerbose) << "AltModuleSliceList::Reco\n";
  
  JobCResult result(JobCResult::kPassed);

  //-- Find PrimaryCandidateRecord fragment in MOM & check it is not null.

  CandRecord *cdrec = dynamic_cast<CandRecord *>
                    (mom->GetFragment("CandRecord", "PrimaryCandidateRecord"));      
  if (cdrec == 0) {
     MSG("Alt", Msg::kWarning) 
                              << "NULL PrimaryCandidateRecord pointer" << endl;
     result.SetWarning().SetFailed();
     return result;
  }

  //-- Get a handle to CandStripList (previous step reconstruction objects).

  MSG("Alt", Msg::kDebug) << "Asking PrimaryCandidateRecord for fListIn: "  
                          << fListIn.Data() << endl; 

  CandStripListHandle *cslh = dynamic_cast<CandStripListHandle *>
                                       (cdrec->FindCandHandle(fListIn.Data()));
  
  //-- Require number of CandStrips to be non-zero.

  if (!cslh) {
     MSG("Alt", Msg::kDebug) << "NULL CandStripListHandle pointer" << endl;
     result.SetFailed();
     return result;         
      
  } else {
      if(cslh->GetNDaughters() < 1) {
          MSG("Alt", Msg::kDebug) << "Empty CandStrip list." << endl;
          result.SetFailed();
          return result;
      }
  }

  //-- Get AlgFactory instance.

  AlgFactory &af = AlgFactory::GetInstance();

  //-- Build a CandSliceList.

  AlgHandle adlh = af.GetAlgHandle("AltAlgSliceList", "default");

  MSG("Alt", Msg::kDebug) << "CandContext cx(this, mom);" << endl;
  CandContext cx(this, mom);

  MSG("Alt", Msg::kDebug) << "cx.SetDataIn(cslh);"        << endl;
  cx.SetDataIn(cslh);
   
  MSG("Alt", Msg::kDebug) << "cx.SetCandRecord(cdrec);"   << endl;   
  cx.SetCandRecord(cdrec);

  //-- Initiate Reco : CandSliceList:MakeCandidate
  MSG("Alt", Msg::kDebug)
            << "csllh = CandSliceList::MakeCandidate(adlh, cx);" << endl;
  CandSliceListHandle csllh = CandSliceList::MakeCandidate(adlh, cx);
  MSG("Alt", Msg::kDebug)
            << "CandSliceListHandle::SetName() & SetTitle()" << endl;
  csllh.SetName(fListOut.Data());
  csllh.SetTitle(
        TString("Brought to life by AltModuleSliceList from candstriplist").
        Append(cslh->GetName())
  );

  //-- Secure the CandHandle in the CandRecord
  MSG("Alt", Msg::kDebug) << "candrec->SecureCandHandle(csllh);" << endl;
  cdrec->SecureCandHandle(csllh);

  return result;
}
//____________________________________________________________________________
void AltModuleSliceList::HandleCommand(JobCommand *command) 
{
  // Process configuration commands
  MSG("Alt", Msg::kDebug) 
                         << "AltModuleSliceList::HandleCommand" << endl;

  TString cmd = command->PopCmd();
  if (cmd == "Set") {
    TString opt = command->PopOpt();
    if      (opt == "ListIn" ) fListIn  = command->PopOpt();
    else if (opt == "ListOut") fListOut = command->PopOpt();
    else {
      MSG("Alt", Msg::kWarning)
             << "AltModuleSliceList::HandleCommand: Unrecognized option "
                                                         << opt << endl;
    }
  } else {
      MSG("Alt", Msg::kWarning)
            << "AltModuleSliceList::HandleCommand: Unrecognized command "
                                                         << cmd << endl;
  }
}
//____________________________________________________________________________

---