FilterLIModule.cxx

Go to the documentation of this file.

#include <cassert>

#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "CandData/CandHeader.h"
#include "CandData/CandRecord.h"
#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigitList.h"
#include "CandDigit/CandDigitListHandle.h"
#include "Candidate/CandContext.h"
#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "RawData/RawRecord.h"
#include "RawData/RawDaqSnarlHeader.h"

#include "FilterLIModule.h"
#include "LIPatternFinder.h"
#include "LIFinderFactory.h"

ClassImp(FilterLIModule)

//______________________________________________________________________
CVSID("$Id: FilterLIModule.cxx,v 1.6 2004/07/09 16:08:20 costas Exp $");
//______________________________________________________________________
JOBMODULE(FilterLIModule, "FilterLIModule",
             "Filters LI-like snarls from physics reconstruction jobs");
//______________________________________________________________________
FilterLIModule::FilterLIModule()
{
  MSG("LIFilter", Msg::kVerbose) << "FilterLIModule::Constructor\n";
}
//______________________________________________________________________
FilterLIModule::~FilterLIModule()
{
  MSG("LIFilter", Msg::kVerbose) << "FilterLIModule::Destructor\n";
}
//______________________________________________________________________
void FilterLIModule::BeginJob()
{
  MSG("LIFilter", Msg::kVerbose) << "FilterLIModule::BeginJob\n";
}
//______________________________________________________________________
void FilterLIModule::Config(const Registry& /* r */)
{
  MSG("LIFilter", Msg::kDebug) << "FilterLIModule::Config" << endl;
}
//______________________________________________________________________
const Registry& FilterLIModule::DefaultConfig() const
{
  MSG("LIFilter", Msg::kDebug)
                             << "FilterLIModule::DefaultConfig" << endl;
  static Registry r;

  std::string name = this->JobCModule::GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  r.UnLockValues();
  
  r.Set("NameListIn",                    "canddigitlist");
  r.Set("NameListOut",                   "canddigitlist");
  r.Set("LIFinderName",                  "SimpleLIPatternFinder");

  //-- inputs relevant to [SimpleLIPatternFinder] only
  r.Set("MaxEWAsymmetry",                 4.e+5);
  r.Set("MaxNDigits",                     1400);
  r.Set("MaxCharge",                      5.e+5);
  r.Set("PulseHeightThreshold",           200);
  r.Set("HighActivityThreshold",          0.94);
  r.Set("LowActivityThreshold",           0.07);

  //-- inputs relevant to [FancyLIPatternFinder] only
  r.Set("FractionalChargeThreshold",      0.8);
   
  r.LockValues();

  return r;
}
//______________________________________________________________________
JobCResult FilterLIModule::Ana(const MomNavigator * mom)
{
  JobCResult result(JobCResult::kPassed);

  MSG("LIFilter", Msg::kVerbose) << "FilterLIModule::Reco" << endl;

// Cache JobModule Registry values for Reco method
  const char * tmps = 0;

  const char * namelistin  = 0;
  const char * namelistout = 0;
  const char * algname     = 0;
  
  Registry &r = GetConfig();     // Get this JobModule's Registry object
  if (r.Get("NameListIn",       tmps)) namelistin       = tmps;
  if (r.Get("NameListOut",      tmps)) namelistout      = tmps;
  if (r.Get("LIFinderName",     tmps)) algname          = tmps;

  //-- Get PrimaryCandidateRecord fragment from MOM.
  CandRecord *candrec = dynamic_cast<CandRecord *>
             (mom->GetFragment("CandRecord", "PrimaryCandidateRecord"));
  if (candrec == 0) {
    MSG("LIFilter", Msg::kWarning)
                          << "No PrimaryCandidateRecord in MOM" << endl;
    result.SetWarning().SetFailed();
    return result;
  }

  CandDigitListHandle *cdlh = dynamic_cast<CandDigitListHandle *>
    (candrec->FindCandHandle("CandDigitListHandle",namelistin));

  if (cdlh == 0) {
    MSG("LIFilter", Msg::kWarning)
                << "No CandDigitListHandle named " << namelistin
                                           << " in CandRecord." << endl;
    result.SetWarning().SetFailed();
    return result;
  }

  //-- Getting an instance to a factory of LIPatternFinder algorithms

  MSG("LIFilter", Msg::kDebug)
                        << "Getting a LIFinderFactory instance" << endl;
                        
  LIFinderFactory * factory = LIFinderFactory::GetInstance();

  //-- Get the algorithm by name
    
  MSG("LIFilter", Msg::kDebug)
           << "Asking LIFinderFactory for ["
                                  << algname << "] algorithm " << endl;

  LIPatternFinder * li_finder = factory->LIFinder(algname);

  //-- Ask the algorithm object if we should pass this event
  
  if(!li_finder) {
    
      MSG("LIFilter", Msg::kWarning)
           << " ** Unspecified LIPatternFinder ** "
           << " Skipping Pattern Finding & passing event to next module"
           << endl;
            
      return JobCResult::kPassed;

  } else {

      li_finder->Configure(r);

      if ( li_finder->IsLightInjectionTrash(cdlh) ) {

         //-- print the rejected snarl number

         RawRecord * rawrec = dynamic_cast<RawRecord *>
                                        (mom->GetFragment("RawRecord"));
         if (rawrec != 0) {

           const RawDaqSnarlHeader * snrlHdr = 
                        dynamic_cast<const RawDaqSnarlHeader *> 
                                               (rawrec->GetRawHeader());
           if(snrlHdr != 0) {
              MAXMSG("LIFilter",Msg::kInfo, 200)
                 << "REJECTING SNARL = " << snrlHdr->GetSnarl() << endl;
           }
         }

         return JobCResult::kFailed;
      }
  }                  

  return JobCResult::kPassed;
}
//______________________________________________________________________

---