AtNuFindModule.cxx

Go to the documentation of this file.

#include  "AtNuFindModule.h"

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"

#include "CandData/CandRecord.h"
#include "Candidate/CandContext.h"
#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "Validity/VldContext.h"
#include "Plex/PlexPlaneId.h"

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

#include "CandAtNuRecoHandle.h"
#include "CandShowerAtNuHandle.h"
#include "CandShowerAtNuListHandle.h"
#include "CandTrackAtNuHandle.h"
#include "CandTrackAtNuListHandle.h"

//
// $Log $
//

ClassImp(AtNuFindModule)

CVSID("$Id: AtNuFindModule.cxx,v 1.12 2006/03/28 13:54:11 blake Exp $");
JOBMODULE(AtNuFindModule, "AtNuFindModule","AtNu Track/Shower Finder Module");

AtNuFindModule::AtNuFindModule() :
  fListIn("CandSliceListHandle"),
  fMCTruthFlag(0)
{
  
}

AtNuFindModule::~AtNuFindModule()
{
 
}

void AtNuFindModule::BeginJob()
{

  MSG("AtNuFindModule", Msg::kInfo) << " *** AtNuFindModule::BeginJob() *** " << endl;

  /*
  AlgFactory &af = AlgFactory::GetInstance();
  af.Register("AlgAtNuReco","default","libAtNuReco.so","AlgConfig");
  af.Register("AlgAtNuRecoMCTruth","default","libAtNuReco.so","AlgConfig");
  af.Register("AlgTrackAtNuList", "default", "libAtNuReco.so", "AlgConfig");
  af.Register("AlgTrackAtNu", "default", "libAtNuReco.so", "AlgConfig");
  af.Register("AlgShowerAtNuList", "default", "libAtNuReco.so", "AlgConfig");
  af.Register("AlgShowerAtNu", "default", "libAtNuReco.so", "AlgConfig");
  */
}

JobCResult AtNuFindModule::Reco(MomNavigator *mom)
{

  // AtNuFindModule::Reco()
  MSG("AtNuFindModule", Msg::kInfo) << " *** AtNuFindModule::Reco() *** " << endl;
  JobCResult result(JobCResult::kPassed);

  // Find PrimaryCandidateRecord fragment in mom.
  CandRecord *candrec = dynamic_cast<CandRecord *> (mom->GetFragment("CandRecord", "PrimaryCandidateRecord"));
  if(!candrec){ 
    MSG("AtNuFindModule", Msg::kWarning) << " Failed to Find CandRecord " << endl;
    return result.SetFailed();
  }

  // Find CandSliceList in PrimaryCandidateRecord.   
  CandSliceListHandle *cslh = dynamic_cast<CandSliceListHandle*>(candrec->FindCandHandle("CandSliceListHandle",fListIn.Data()));
  if(!cslh){
    MSG("AtNuFindModule", Msg::kWarning) << " Failed to Find " << fListIn.Data() << endl;
    return result.SetFailed();
  }

  // Make CandAtNuRecoHandle
  MSG("AtNuFindModule", Msg::kDebug) << " Make CandAtNuRecoHandle " << endl;
  TString alg_atnu("AlgAtNuReco"); if(fMCTruthFlag) alg_atnu.Append("MCTruth");
  MSG("AtNuFindModule", Msg::kDebug) << "     ... using algorithm " << alg_atnu.Data() << endl;

  AlgFactory &af = AlgFactory::GetInstance();
  AlgHandle ah_atnu = af.GetAlgHandle(alg_atnu.Data(), "default");
  
  CandContext cx_atnu(this, mom);
  cx_atnu.SetCandRecord(candrec);
  cx_atnu.SetDataIn(cslh);
  CandAtNuRecoHandle candatnu = CandAtNuReco::MakeCandidate(ah_atnu, cx_atnu);
  candatnu.SetName("CandAtNuRecoHandle");
  candatnu.SetTitle(TString("Created by AtNuFindModule from ").Append(cslh->GetName()));
  
  candrec->SecureCandHandle(candatnu);

  MSG("AtNuFindModule", Msg::kInfo) << " *** AtNuFindModule::Reco() FINISHED *** " << endl;

  return result;

}

JobCResult AtNuFindModule::Ana(const MomNavigator * /*mom*/)
{

  JobCResult result(JobCResult::kPassed);


  return result;

}

const Registry& AtNuFindModule::DefaultConfig() const
{

  MSG("AtNuFindModule", Msg::kInfo) << " *** AtNuFindModule::DefaultConfig() *** " << endl;

  static Registry r;
  r.SetName("AtNuFindModule.config.default");
  r.UnLockValues();
  r.Set("ListIn",fListIn.Data());
  r.Set("MCTruthFlag",fMCTruthFlag);
  r.Set("ListOutTrk","CandTrackAtNuListHandle");
  r.Set("ListOutShw","CandShowerAtNuListHandle");
  r.Set("VtxFlag",1);
  r.Set("PeCut",1.0);
  r.Set("FibreIndex",1.77);
  r.Set("AtNuAnaOnOff",1);
  r.LockValues();

  AlgFactory &af = AlgFactory::GetInstance();

  af.Register("AlgAtNuReco","default");
  AlgHandle ah_atnu = af.GetAlgHandle("AlgAtNuReco","default");
  AlgConfig &ac_atnu = ah_atnu.GetAlgConfig();
  ac_atnu.UnLockValues();
  ac_atnu.Set("ListOutTrk","CandTrackAtNuListHandle");
  ac_atnu.Set("ListOutShw","CandShowerAtNuListHandle");
  ac_atnu.Set("VtxFlag",1);
  ac_atnu.Set("PeCut",1.0);
  ac_atnu.LockValues();

  af.Register("AlgAtNuRecoMCTruth","default");
  AlgHandle ah_atnumc = af.GetAlgHandle("AlgAtNuRecoMCTruth","default");
  AlgConfig &ac_atnumc = ah_atnumc.GetAlgConfig();
  ac_atnumc.UnLockValues();
  ac_atnumc.Set("ListOutTrk","CandMCTrackAtNuListHandle");
  ac_atnumc.Set("ListOutShw","CandMCShowerAtNuListHandle");
  ac_atnumc.LockValues();

  af.Register("AlgTrackAtNu", "default");
  AlgHandle ah_trk = af.GetAlgHandle("AlgTrackAtNu","default");
  AlgConfig &ac_trk = ah_trk.GetAlgConfig();
  ac_trk.UnLockValues();
  ac_trk.Set("FibreIndex",1.77);
  ac_trk.Set("AtNuAnaOnOff",1);
  ac_trk.LockValues();

  af.Register("AlgShowerAtNu", "default");
  AlgHandle ah_shw = af.GetAlgHandle("AlgShowerAtNu","default");
  AlgConfig &ac_shw = ah_shw.GetAlgConfig();
  ac_shw.UnLockValues();
  ac_shw.Set("FibreIndex",1.77);
  ac_shw.Set("AtNuAnaOnOff",1);
  ac_shw.LockValues();

  return r;
}

void AtNuFindModule::Config(const Registry &r)
{

  MSG("AtNuFindModule", Msg::kInfo) << " *** AtNuFindModule::Config() *** " << endl;

  Int_t tmpint; Double_t tmpdub; const char* tmpchar = 0;

  if(r.Get("ListIn",tmpchar)) fListIn = tmpchar;
  if(r.Get("MCTruthFlag",tmpint)) fMCTruthFlag = tmpint;

  AlgFactory &af = AlgFactory::GetInstance();

  AlgHandle ah_atnu = af.GetAlgHandle("AlgAtNuReco","default");
  AlgConfig &ac_atnu = ah_atnu.GetAlgConfig();
  ac_atnu.UnLockValues();
  if(r.Get("ListOutTrk",tmpchar)) ac_atnu.Set("ListOutTrk",tmpchar);
  if(r.Get("ListOutShw",tmpchar)) ac_atnu.Set("ListOutShw",tmpchar);
  if(r.Get("VtxFlag",tmpint)) ac_atnu.Set("VtxFlag",tmpint);
  if(r.Get("PeCut",tmpdub)) ac_atnu.Set("PeCut",tmpdub);
  ac_atnu.LockValues();

  AlgHandle ah_atnumc = af.GetAlgHandle("AlgAtNuRecoMCTruth","default");
  AlgConfig &ac_atnumc = ah_atnumc.GetAlgConfig();
  ac_atnumc.UnLockValues();
  if(r.Get("ListOutTrk",tmpchar)) ac_atnumc.Set("ListOutTrk",tmpchar);
  if(r.Get("ListOutShw",tmpchar)) ac_atnumc.Set("ListOutShw",tmpchar);
  ac_atnumc.LockValues();

  AlgHandle ah_trk = af.GetAlgHandle("AlgTrackAtNu","default");
  AlgConfig &ac_trk = ah_trk.GetAlgConfig();
  ac_trk.UnLockValues();
  if(r.Get("FibreIndex",tmpdub)) ac_trk.Set("FibreIndex",tmpdub);
  if(r.Get("AtNuAnaOnOff",tmpint)) ac_trk.Set("AtNuAnaOnOff",tmpint);
  ac_trk.LockValues();

  AlgHandle ah_shw = af.GetAlgHandle("AlgShowerAtNu","default");
  AlgConfig &ac_shw = ah_shw.GetAlgConfig();
  ac_shw.UnLockValues();
  if(r.Get("FibreIndex",tmpdub)) ac_shw.Set("FibreIndex",tmpdub);
  if(r.Get("AtNuAnaOnOff",tmpint)) ac_shw.Set("AtNuAnaOnOff",tmpint);
  ac_shw.LockValues();

  return;
}

void AtNuFindModule::HandleCommand(JobCommand *command)
{
  TString cmd = command->PopCmd();
  if(cmd=="Set"){
    TString opt = command->PopOpt();

  }
}

void AtNuFindModule::EndJob()
{   

}     

---