MeuCalModule.cxx

Go to the documentation of this file.

// Coded by Jeff Hartnell Apr/2007 onwards
//
// Contact: j.j.hartnell@rl.ac.uk

#include <cassert>

#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "BeamDataNtuple/NtpBDLiteRecord.h"
#include "StandardNtuple/NtpStRecord.h"

#include "MeuCal/MeuAnalysis.h"
#include "MeuCal/MeuCalModule.h"

JOBMODULE(MeuCalModule,
          "MeuCalModule",
          "A module to produce MeuSummary ntuples");

CVSID("$Id: MeuCalModule.cxx,v 1.2 2007/04/04 09:08:57 hartnell Exp $");

ClassImp(MeuCalModule)

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

MeuCalModule::MeuCalModule() 
{
  MSG("MeuCalModule", Msg::kDebug) 
   <<"Running MeuCalModule constructor..."<<endl;

  //initialise data members
  fMeuAnalysis=0;

  MSG("MeuCalModule", Msg::kDebug) 
   <<"Finished constructor"<<endl;
}

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

MeuCalModule::~MeuCalModule()
{
  MSG("MeuCalModule", Msg::kDebug) 
   <<"Running destructor..."<<endl;

  MSG("MeuCalModule", Msg::kDebug) 
   <<"Finished destructor"<<endl;
}

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

JobCResult MeuCalModule::Ana(const MomNavigator *mom)
{
  //to test number of entries
  //static Int_t counter=0;
  //counter++;
  //cout<<"counter="<<counter<<endl;

  //set default result to kPassed
  JobCResult result(JobCResult::kPassed);
  
  //check that mom exists.
  assert(mom);
  
  const NtpStRecord* pntp=dynamic_cast<NtpStRecord*>
    (mom->GetFragment("NtpStRecord"));
  if (!pntp){
    MSG("MeuCalModule",Msg::kWarning) 
      << "No NtpStRecord in Mom" << endl;
  }
  
  const NtpBDLiteRecord* pntpBD=dynamic_cast<NtpBDLiteRecord*>
    (mom->GetFragment("NtpBDLiteRecord"));
  if (!pntpBD){
    MAXMSG("MeuCalModule",Msg::kInfo,1)
      <<"No NtpBDLiteRecord in Mom (note: not present for MC)"<<endl;
  }
  
  //run the code to calc the meu values for events in this snarl
  fMeuAnalysis->MakeSummaryTreeWithNtpStOneSnarl(pntp,pntpBD);
  
  return result;
}

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

void MeuCalModule::BeginJob() 
{
  MSG("MeuCalModule",Msg::kInfo)
    <<"Running MeuCalModule::BeginJob()..."<<endl;

  //construct the MeuAnalysis object
  fMeuAnalysis=new MeuAnalysis();

  MSG("MeuCalModule",Msg::kInfo)
    <<"Finished running MeuCalModule::BeginJob()"<<endl;
}

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

void MeuCalModule::EndJob() 
{  
  MSG("MeuCalModule",Msg::kInfo)
    <<"Running MeuCalModule::EndJob()..."<<endl;

  //finish the job
  fMeuAnalysis->StoreOrFinishSummaryTree(NULL,NULL,1);

  //write out any histos if a file was opened
  fMeuAnalysis->WriteOutHistos();
  
  MSG("MeuCalModule",Msg::kInfo)
    <<"Finished running MeuCalModule::EndJob()"<<endl;
}

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

const Registry& MeuCalModule::DefaultConfig() const
{

  MSG("MeuCalModule", Msg::kDebug)
    << "Running MeuCalModule::DefaultConfig..." << endl;
  
  static Registry r;
  
  string name=this->JobCModule::GetName();
  name+=".config.default";
  r.SetName(name.c_str());
  r.UnLockValues();
  
  //r.Set("LIRunSeqRequired",-1);
  //r.Set("WriteOnlyGainCurves",0);
  
  r.LockValues();
  
  return r;
}

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

void MeuCalModule::Config(const Registry& /*r*/)
{

  MSG("MeuCalModule",Msg::kDebug)<<"MeuCalModule::Config"<<endl;

  /*
    Int_t runSeqReq=-1;
    if (r.Get("LIRunSeqRequired",runSeqReq)) fLIRunSeqReq=runSeqReq;
  
  Int_t writeOnlyGainCurves=0;
  if (r.Get("WriteOnlyGainCurves",
  writeOnlyGainCurves)) fWriteOnlyGainCurves=writeOnlyGainCurves;
  */
}

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

---