FiltTriggerPrescale.cxx

Go to the documentation of this file.

// $Id: FiltTriggerPrescale.cxx,v 1.2 2005/06/02 06:50:29 rhatcher Exp $
//
// A filter module that accepts/rejects DaqSnarl records based on
// prescaling of thier trigger source.
//
// rhatcher@fnal.gov
#include "Filtration/FiltTriggerPrescale.h"
#include <cstdio>
using namespace std;

// ROOT includes
#include "TFile.h"
#include "TStyle.h"
#include "TMath.h"
#include "TString.h"  // for Form()

// MINOS includes
#include "MessageService/MsgService.h"     // MSG text output
#include "MinosObjectMap/MomNavigator.h"   // Data access
#include "Util/UtilString.h"

#include "RawData/RawRecord.h"
#include "RawData/RawDaqSnarlHeader.h"
#include "RawData/RawTriggerCodes.h"

//#include "Conventions/Munits.h"            // Unit conversions
#include "JobControl/JobCommand.h"         // JobCommand handling
#include "JobControl/JobCModuleRegistry.h" // JOBMODULE registration macro

CVSID("$Id: FiltTriggerPrescale.cxx,v 1.2 2005/06/02 06:50:29 rhatcher Exp $");
JOBMODULE(FiltTriggerPrescale,"FiltTriggerPrescale","Select events on prescaled trigger counts.\n");

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

FiltTriggerPrescale::FiltTriggerPrescale()
{
    fRecordSets = fSnarls = fSnarlsPassed = fNoSnarlRecSets = 0;
    for (int i=0; i<32; ++i) {
        fN[i] = fM[i] = fTrigSeen[i] = fTrigPass[i] = 0;
    }
}

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

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

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

  // Set values of config
  r.UnLockValues();

  for (int ibit=0; ibit<32; ++ibit) {
      UInt_t oneBit = 1 << ibit;
      RawTriggerCodes::TrigBits_t tbit = 
          static_cast<RawTriggerCodes::TrigBits_t>(oneBit);
      std::string trigName = RawTriggerCodes::AsString(tbit);
      int defaultN = 1;
      int defaultM = 1;

      // handle unnamed trigger bits (default for this is "ignore")
      // but allow users to configure with "TrigBit<n>" naming
      if ( trigName.find("TrigBit") == 0 || trigName.find("0x") == 0 ) {
          trigName = "TrigBit" + UtilString::ToString(ibit);
          defaultN = 0;
          defaultM = 0;
      }

      std::string nName = trigName + "_N";
      std::string mName = trigName + "_M";
      r.Set(nName.c_str(),defaultN);
      r.Set(mName.c_str(),defaultM);
  }

  r.LockValues();

  return r;
}

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

void FiltTriggerPrescale::Config(const Registry& r)
{
//======================================================================
// Configure the module given the registry r
//======================================================================
//  char   tmpb;
//  double tmpd;
    int    tmpn, tmpm;

  for (int ibit=0; ibit<32; ++ibit) {
      UInt_t oneBit = 1 << ibit;
      RawTriggerCodes::TrigBits_t tbit = 
          static_cast<RawTriggerCodes::TrigBits_t>(oneBit);
      std::string trigName = RawTriggerCodes::AsString(tbit);

      // handle unnamed trigger bits
      if ( trigName.find("TrigBit") == 0 || trigName.find("0x") == 0 ) {
          trigName = "TrigBit" + UtilString::ToString(ibit);
      }
      
      std::string nName = trigName + "_N";
      std::string mName = trigName + "_M";
      if (r.Get(nName.c_str(),tmpn)) { fN[ibit] = tmpn; }
      if (r.Get(mName.c_str(),tmpm)) { fM[ibit] = tmpm; }
  }

}

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

void FiltTriggerPrescale::Report() 
{
  MSG("Filt",Msg::kInfo) 
    << "Report For FiltTriggerPrescale Module:" << endl;

  //Registry &r = GetConfig();
  //r.Print();

  MSG("Filt",Msg::kInfo) 
      << " Processed " << fRecordSets << " record-sets "
      << "(" << fNoSnarlRecSets << " without snarls)." << endl;

  MSG("Filt",Msg::kInfo) 
      << "  bit#  triggerName      passed/    seen                 N/   M" 
      << endl;
  for ( int ibit=0; ibit<32; ++ibit) {
      UInt_t nPass = fTrigPass[ibit];
      UInt_t mSeen = fTrigSeen[ibit];
      UInt_t n     = fN[ibit];
      UInt_t m     = fM[ibit];
      if ( nPass == 0 && mSeen == 0 && n == 0 && m == 0 ) continue;

      UInt_t oneBit = 1 << ibit;
      RawTriggerCodes::TrigBits_t tbit = 
          static_cast<RawTriggerCodes::TrigBits_t>(oneBit);
      std::string trigName = RawTriggerCodes::AsString(tbit);

      double passFrac = (mSeen>0) ? (double)nPass/(double)mSeen : -1;
      double targFrac = (m    >0) ? (double)n    /(double)m     : -1;

      MSG("Filt",Msg::kInfo)
          << "  [" << setw(2) << right << ibit << left << "]  "
          << setw(12) << trigName << "  "
          << right
          << " " << setw(8) << nPass << "/" << setw(8) << mSeen
          << " (" << setw(10) << passFrac << ")"
          << " " << setw(4) << n     << "/" << setw(4) << m
          << " (" << setw(10) << targFrac << ")"
          << left
          << endl;
  }
  double totPassFrac = 
      (fSnarls>0) ? (double)fSnarlsPassed/(double)fSnarls : -1;
  MSG("Filt",Msg::kInfo)
      << "        "
      << setw(12) << "-all triggers-" 
      << right
      << " " << setw(8) << fSnarlsPassed << "/" << setw(8) << fSnarls
      << " (" << setw(10) << totPassFrac << ")"
      << endl;

  MSG("Filt",Msg::kInfo) << endl;
}

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

void FiltTriggerPrescale::EndJob() 
{
    Report();
}

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

JobCResult FiltTriggerPrescale::Ana(const MomNavigator *mom) 
{

  JobCResult result = JobCResult::kFailed;
  int snarlsInRecSet = 0;

  // count record sets that have been processed
  fRecordSets++;
  
  TObject* recObj = 0;
  TIter recItr = mom->FragmentIter();
  while ( ( recObj = recItr.Next() ) ) {
    RawRecord* rawrec = dynamic_cast<RawRecord*>(recObj);
    if (!rawrec) continue;  // skip things in mom that aren't RawRecords

    // look for the DaqSnarl records
    
    const RawDaqSnarlHeader* snarlHdr = 
      dynamic_cast<const RawDaqSnarlHeader*>(rawrec->GetHeader());

    if ( snarlHdr ) {
      fSnarls++;
      snarlsInRecSet++;  // count # of snarls in recset
      bool passSnarl = false;

      int snarlNum  = snarlHdr->GetSnarl();
      int trigSrc   = snarlHdr->GetTrigSrc();
      const VldContext& vldc = snarlHdr->GetVldContext();
      std::string detName =  Detector::AsString(vldc.GetDetector());

      for (int ibit=0; ibit<32; ++ibit) {
          UInt_t oneBit = 1 << ibit;
          if ( oneBit & trigSrc ) {
              fTrigSeen[ibit]++;
              bool passBit = false;
              if ( fM[ibit] > 0 ) {
                  UInt_t modSnarlNum = snarlNum%fM[ibit];
                  if ( modSnarlNum < fN[ibit] ) passBit = true;
              }

              RawTriggerCodes::TrigBits_t tbit = 
                  static_cast<RawTriggerCodes::TrigBits_t>(oneBit);
              std::string trigName = RawTriggerCodes::AsString(tbit);
              MSG("Filt",Msg::kDebug )
                  << " " << detName
                  << " Run " << snarlHdr->GetRun() 
                  << " Snarl " << snarlNum
                  << " had TrigBit \"" << trigName << "\""
                  << " : "
                  << (passBit?"Passed":"Failed")
                  << endl;

              if ( passBit ) {
                  passSnarl = true;
                  fTrigPass[ibit]++;
              }
          }
      }
      if ( passSnarl ) {
          fSnarlsPassed++;
          result = JobCResult::kPassed;
      }

      MSG("Filt",Msg::kSynopsis )
          << " " << detName
          << " Run " << snarlHdr->GetRun() 
          << " SubRun " << snarlHdr->GetSubRun()
          << " Snarl " << snarlNum
          << " TrigSrc " << Form("0x%08x",trigSrc)
          << " : "
          << (passSnarl?"Passed":"Failed")
          << endl;

    } // snarl header
  }

// pass record sets that contain no DaqSnarl
  if ( snarlsInRecSet == 0 ) {
      fNoSnarlRecSets++;
      MSG("Filt",Msg::kSynopsis )
          << GetName() 
          << " RecordSet contained no DaqSnarl "
          << endl;
      return JobCResult::kPassed;
  }
  if ( snarlsInRecSet > 1 ) {
      MSG("Filt",Msg::kWarning )
          << GetName() 
          << " RecordSet contained more than one DaqSnarl "
          << endl;
  }
  return result;
}

---