PrintStat.cxx

Go to the documentation of this file.

// $Id: PrintStat.cxx,v 1.3 2004/09/11 17:35:43 weber Exp $
//
// A small bit of program to print out some status line for each event
//
// A.Weber@rl.ac.uk
#include "PrintStat.h"
#include <cstdio>
// MINOS includes
#include "MessageService/MsgService.h"     // MSG text output
#include "MinosObjectMap/MomNavigator.h"   // Data access
#include "RawData/RawRecord.h"
#include "RawData/RawDaqSnarlHeader.h"
//#include "JobControl/JobCommand.h"         // JobCommand handling
#include "JobControl/JobCModuleRegistry.h" // JOBMODULE registration macro
#include "Conventions/Munits.h"
#include "Validity/VldContext.h"
#include "TString.h"

CVSID("$Id: PrintStat.cxx,v 1.3 2004/09/11 17:35:43 weber Exp $");
JOBMODULE(PrintStat,"PrintStat","Analyse Raw Hits\n");

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

PrintStat::PrintStat() :
  fActOutLine(-1)      // event counter
{
//======================================================================
// Set the initial state of your module. The code above initializes the
// variables listed to the values in parentheses. 
//======================================================================

  MSG("Stat",Msg::kDebug) << "Constructing PrintStat.\n";

}

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


PrintStat::~PrintStat()
{
//======================================================================
// This is my destructor, which currently does nothing.
//======================================================================

  MSG("Stat",Msg::kDebug) << "Destructing PrintStat.\n";

}

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

void PrintStat::BeginJob()
{
//======================================================================
// This is my destructor, which currently does nothing.
//======================================================================

  MSG("Stat",Msg::kDebug) << "PrintStat::BeginJob starts here.\n";

}

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

const Registry& PrintStat::DefaultConfig() const
{
//======================================================================
// Create a registry which holds the default configuration and return it
//======================================================================
  static Registry r;

  MSG("Stat",Msg::kDebug) << "PrintStat::DefaultConfig starts here.\n";

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

  // Set values of config
  r.UnLockValues();
  r.Set("OutputRate",  10);
  r.Set("MaxOutLine", 100000);
  r.LockValues();

  return r;
}

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

void PrintStat::Config(const Registry& r)
{
//======================================================================
// Configure the module given the registry r
// Do all configuration here.
// Beware: this methode may be called several times!
//======================================================================

  int    tmpi;

  if (r.Get("OutputRate", tmpi)) { fOutputRate  = tmpi; }
  if (fOutputRate<1) fOutputRate = 1;

  if (r.Get("MaxOutLine", tmpi)) { fMaxOutLine  = tmpi; }

}

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

JobCResult PrintStat::Ana(const MomNavigator *mom) 
{
//======================================================================
// Given the object to which the data is attached perform some
// analysis of the data
//======================================================================

  MSG("Stat", Msg::kDebug) << "Starting PrintStat::Ana.\n";

  // Accumulate statistics

  fActOutLine++;
  if (fActOutLine%fOutputRate!=0)           return JobCResult::kPassed;
  if (fActOutLine/fOutputRate>fMaxOutLine)  return JobCResult::kPassed;

  // find some pointers

  TObject* tobj;
  TIter    fragiter = mom->FragmentIter();

  RawRecord*                 rr   = 0;
  const RawDaqSnarlHeader*   rdsh = 0;

  // loop over all fragments of mom and get pointers

  while( ( tobj = fragiter.Next() ) ) { 
    if( ( rr = dynamic_cast<RawRecord*>(tobj) ) ) {

      MSG("Stat", Msg::kDebug) << "Found Next RawRecord. rr= "
                                 << rr << endl;

      if ( !rdsh )  rdsh  = dynamic_cast<const RawDaqSnarlHeader*>
                      (rr->GetRawHeader());
    } else {
      MSG("Stat", Msg::kVerbose) << "Not a RawRecord.\n";

    }  // if rawrec
  } // while


  if ( rdsh ) {

    Int_t        Run   = rdsh->GetRun();
    Int_t        Snarl = rdsh->GetSnarl();
    VldContext   vldc  = rdsh->GetVldContext();
    VldTimeStamp vldts = vldc.GetTimeStamp();

    MSG("Stat",Msg::kInfo) << "===>> " 
                           << SimFlag::AsString(vldc.GetSimFlag())
                           << " Snarl(" << setw(5) << fActOutLine << "):"
                           << " -- run = "    << setw(6) << Run
                           << " -- snarl = "  << setw(6) << Snarl
                           << " -- time = "   << vldts.AsString("c")
                           << " <<===" << endl;
  } else {

    MSG("Stat",Msg::kDebug) << "No RawDaqSnarlHeader.";

  }
  
  return JobCResult::kPassed;

}

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

void PrintStat::EndJob() 
{
//======================================================================
// At the end of the job print some stuff out and save the histogram
// to a file
//======================================================================

  MSG("Stat",Msg::kInfo) << "PrintStat: Total Events processed = " 
                         << fActOutLine << "\n";

}

---