FiltCosmicModule.cxx

Go to the documentation of this file.

#include "FiltCosmicModule.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro
#include "JobControl/JobCommand.h"  
#include "Conventions/Munits.h"

#include "CandData/CandHeader.h"
#include "RecoBase/CandEventListHandle.h"
#include "RecoBase/CandEventHandle.h"

// ROOT includes
#include "TSystem.h"
#include "TFile.h"
#include "TTree.h"
#include "TMath.h"

JOBMODULE(FiltCosmicModule, "FiltCosmicModule",
          "Filt");
CVSID("$Id: FiltCosmicModule.cxx,v 1.1 2008/09/17 02:56:19 tjyang Exp $");
//......................................................................

FiltCosmicModule::FiltCosmicModule() 
{

}
//......................................................................

FiltCosmicModule::~FiltCosmicModule()
{
}

//......................................................................

void FiltCosmicModule::BeginJob()
{
}

//......................................................................

void FiltCosmicModule::EndJob()
{
}

//......................................................................

void FiltCosmicModule::BeginFile()
{
}

//......................................................................

void FiltCosmicModule::EndFile()
{
}

//......................................................................

void FiltCosmicModule::BeginRun()
{
}

//......................................................................

void FiltCosmicModule::EndRun()
{
}

//......................................................................

JobCResult FiltCosmicModule::Ana(const MomNavigator* mom)
{
  JobCResult result = JobCResult::kFailed;

  CandRecord *record = dynamic_cast<CandRecord *>(mom->GetFragment("CandRecord", "PrimaryCandidateRecord"));
 
  if (record == 0) {
    MSG("Filt",Msg::kError) << "No PrimaryCandidateRecord CandRecord found." << endl;
    return result;
  }

  const CandHeader* header  = dynamic_cast<const CandHeader*>(record->GetHeader());
  if(header){
    MSG("Filt",Msg::kDebug) << " Run/Snarl: " << header->GetRun() << " / "<< header->GetSnarl() <<endl;
  }

  const CandEventListHandle * eventlist = dynamic_cast<CandEventListHandle*>(record->FindCandHandle("CandEventListHandle")); 

  if (eventlist==NULL) {
    MSG("Filt",Msg::kDebug) << " Bailing out of Event eventlist = " << eventlist << endl;
    return result;
  }

  TIter event_iter(eventlist->GetDaughterIterator());

  // for cosmic ray, we only test the first event of the snarl
  const CandEventHandle* event = dynamic_cast<const CandEventHandle*>(event_iter());  
  
  if (event == NULL) {
    MSG("Filt",Msg::kDebug) << " Bailing out of Event = " << event << endl;
    return result;
  } else {
    MSG("Filt",Msg::kDebug) << " New event " <<endl;
  
    const int nshowers = event->GetLastShower() + 1;
    const int ntracks = event->GetLastTrack() + 1;

    MSG("Filt",Msg::kDebug) << " ntrack : " << ntracks << " nshower : " << nshowers << endl;

    if (nshowers<1 || ntracks<1) {
      return result; 
    } else {
      MSG("Filt",Msg::kDebug) << "event passed" << endl;
      result = JobCResult::kPassed; 
    }
  }
  return result;
}

//......................................................................

JobCResult FiltCosmicModule::Reco(MomNavigator* /* mom */)
{
  return JobCResult::kPassed; // kNoDecision, kFailed, etc.
}

//......................................................................

void FiltCosmicModule::Report()
{
}

//......................................................................

void FiltCosmicModule::Reset()
{
}

//......................................................................

void FiltCosmicModule::HandleCommand(JobCommand* /* c */)
{
}

---