PulserSummaryModule.cxx

Go to the documentation of this file.

// $Id: PulserSummaryModule.cxx,v 1.13 2005/09/19 18:28:12 murgia Exp $
//
// PulserSummaryModule.cxx
//
// Makes List of averages of pulser summaries
//
// Author: Phil Adamson 4/2001
//
// $Log: PulserSummaryModule.cxx,v $
// Revision 1.13  2005/09/19 18:28:12  murgia
// Remove fSummaryList if no triggers are found, before call to Finish.
//
// Revision 1.12  2005/08/16 15:51:50  tjyang
// Save the last point in the last subrun if the run stopped between two LED pulses. This bug was pointed out by Pat Ward. Thanks Robert H. for showing how to pass an object between two modules.
//
//
//
// Revision 1.11  2004/11/08 22:02:42  murgia
// If empty summaries are encountered, do not calculate trigfrac and print a warning message.
//
// Revision 1.10  2004/10/19 20:33:50  murgia
// Fix bug in Reco: when matching rawLI to fSummaryList, check if it is the same pulser box and whether it is drift point or gain curve to avoid adding pieces of gain curve to drift point.
//
// Revision 1.9  2004/10/05 21:39:48  murgia
// Fixed code to write points with low light levels. Added call to EndFile.
//
// Revision 1.8  2004/09/08 21:50:51  pa
// Checkin of new PulserCalibration version. MIght even work
//
// Revision 1.7  2003/09/03 09:31:16  rjn
// Fixed a couple of little bugs in the code that stopped things working.
//
// Revision 1.6  2003/02/05 22:56:24  rjn
// Changed Pulser Summary List To Output Box As Well As LEd (cosmetic change
// -- not sure if I have CVS write access)
//
// Revision 1.5  2002/10/17 05:55:53  rhatcher
// various mods to reduce the number of complaints about comparing
// signed and unsigned values
//
// Revision 1.4  2002/10/01 16:58:32  phil.adamson
// PulserTimePlotter isn't workimng, but it compiles, so it won't break
// anything else. Cosmetic changes.
//
// Revision 1.3  2002/08/07 16:16:52  phil.adamson
// Minor cosmetic details, plus provide ZeroCorr() methods in PulserSummary.
// Phil.
//
//

#include "JobControl/JobCModuleRegistry.h"
#include "MinosObjectMap/MomNavigator.h"
#include "PulserCalibration/PulserSummaryModule.h"
#include "MessageService/MsgService.h"
#include "RawData/RawLIAdcSummaryBlock.h"
#include "RawData/RawRecord.h"
#include "RawData/RawHeader.h"
#include "RawData/RawDaqHeaderBlock.h"

#include "TROOT.h"
#include "TFolder.h"

#include <cstdio>
#include <iostream>
#include <fstream>
ClassImp(PulserSummaryModule)

  // ......

CVSID("$Id: PulserSummaryModule.cxx,v 1.13 2005/09/19 18:28:12 murgia Exp $");
JOBMODULE(PulserSummaryModule,"PulserSummaryModule","Average over subaverages of pulser data from online");

PulserSummaryModule::PulserSummaryModule()
{
  MSG("Pulser",Msg::kVerbose) << "PulserSummaryModule ctor\n";
  fSummaryList=0;
  fFirst = true;
}

PulserSummaryModule::~PulserSummaryModule()
{
  MSG("Pulser",Msg::kVerbose) << "PulserSummaryModule dtor\n";
}

JobCResult PulserSummaryModule::Reco(MomNavigator *mom)
{
  MSG("Pulser",Msg::kVerbose) << "PulserSummaryModule::Reco()\n";
  
  // Find RawRecord fragment in MOM.
  TIter iter = mom->FragmentIter();
  while (TObject *obj = iter.Next()) {
    RawRecord *rawrec = dynamic_cast<RawRecord *>(obj);
    if (rawrec) {
      MSG("Pulser",Msg::kVerbose) << "PSM::Reco() got vld: " <<(rawrec->GetRawHeader())->GetVldContext() <<endl;
      int sec = rawrec->GetHeader()->GetVldContext().GetTimeStamp().GetSec();
      // look for a RawLIAdcSummaryBlock in the RawRecord
      VldTimeStamp ts = rawrec->GetHeader()->GetVldContext().GetTimeStamp();
      const RawLIAdcSummaryBlock *rawsum = 
        dynamic_cast<const RawLIAdcSummaryBlock *>
        (rawrec->FindRawBlock("RawLIAdcSummaryBlock"));

      if (rawsum) {
        // Got a LI summary
        MSG("Pulser",Msg::kVerbose) << 
          "PSM::Reco() : Got a LI summary" << sec<<endl;
        fLastSec = sec;


        // Code to use liMaster.log obsolete so now deleted

        //OK - now rawsum definitely contains its header info.
        // if (this rawLI doesn't match fSUmmaryList, dispatch list
        // elsewhere - if time -fLastSec >5, dispatch...
        if (fSummaryList) {
          //      if (rawsum->GetLed()==fSummaryList->GetLed() && 
          //          rawsum->GetCalibPoint()==fSummaryList->GetPointIndex()) {
          // Match pulser box and whether it is a drift point or gain curve
          if (rawsum->GetLed()==fSummaryList->GetLed() && 
              rawsum->GetCalibPoint()==fSummaryList->GetPointIndex() &&
              rawsum->GetCalibType()==fSummaryList->GetExpectedPoints() &&
              rawsum->GetPulserBox()==fSummaryList->GetPulserBox()) {
            fSummaryList->Add(*rawsum);
            
            MSG("Pulser",Msg::kDebug) << 
              "PSM::Reco() : Existing Point. Pulser Box, Led" << rawsum->GetPulserBox() << ", " << rawsum->GetLed()<<endl;

          }
          else {
            // We're on a new point now!
            MSG("Pulser",Msg::kDebug)<<"New point at "<<sec<<endl;
            MSG("Pulser",Msg::kDebug) << 
              "PSM::Reco() : New Point. Pulser Box, Led" << rawsum->GetPulserBox() << ", " << rawsum->GetLed()<<endl;
            FinishList(mom,ts);
            fSummaryList = new PulserSummaryList(*rawsum);
          }
        } 
        else {
          fSummaryList = new PulserSummaryList(*rawsum);
        }
      }
      
      if (fSummaryList&&(sec-fLastSec>10)) {
         // 4 secs without an LI summary? Must be end of point...
         MSG("Pulser",Msg::kInfo)<<"Point timed out "<<sec<<endl;
         FinishList(mom,ts);
      }
      
      
    }
  }
  return JobCResult::kAOK;
}


void PulserSummaryModule::FinishList(MomNavigator* mom, VldTimeStamp ts)
{
  MSG("Pulser",Msg::kInfo) <<ts<<" Summary Complete! Box " << fSummaryList->GetPulserBox() << " Led "<< fSummaryList->GetLed()<<", point "<<fSummaryList->GetPointIndex()<<" of "<<fSummaryList->GetExpectedPoints()<<endl;

  // Check if there are triggers:
  if (fSummaryList->GetTriggers()==0)  {
    MSG("Pulser",Msg::kDebug)<<"Found summary with no triggers. Triggers, Pulses: "
                             <<fSummaryList->GetTriggers() << ", "
                             <<fSummaryList->GetNumPulses() << endl;
    MSG("Pulser",Msg::kDebug)<<"LED, PB, PH, PW: "
                             <<fSummaryList->GetLed() << ", "
                             <<fSummaryList->GetPulserBox() << ", "
                             <<fSummaryList->GetHeight() << ", "
                             <<fSummaryList->GetWidth() << endl;

    fSummaryList=0;
  }

  if (fSummaryList) {
    fSummaryList->Finish();
    
    // Test to see if SummaryList is OK here - squeal if it isn't.
    
    float trigfrac; 
    if (fSummaryList->GetExpectedPoints()!=0) {
      // Do we have enough triggers?
      trigfrac = float(fSummaryList->GetTriggers())/float(fSummaryList->GetNumPulses());
    } else {
      MSG("Pulser",Msg::kWarning)<<"Found empty summary. Triggers, ExpectedPoints, Pulses: "
                                 <<fSummaryList->GetTriggers() << ", "<< fSummaryList->GetExpectedPoints()<<", "
                                 <<fSummaryList->GetNumPulses() << endl;
      trigfrac = 99999.;
    }
    
    if (trigfrac < 0.95) {
      MSG("Pulser",Msg::kError)<<"Only saw "<<trigfrac<<" of expected "<<fSummaryList->GetNumPulses()<<" triggers. Dumping point.\n";
      delete fSummaryList;
      fSummaryList=0;
      return;
    }
  }

  mom->AdoptFragment(fSummaryList);
  fSummaryList=0;
}

void PulserSummaryModule::BeginRun()
{
  MSG("Pulser",Msg::kDebug) << "PulserSummaryModule::BeginRun()"<<endl;
  fSummaryList=0;
}

void PulserSummaryModule::EndRun()
{
//  // Drop part-completed point at end of run, as liMaster will repeat it.
//  delete fSummaryList;
//  fSummaryList = 0;
  if (!fSummaryList) return;

  fSummaryList->Finish();

  // Test to see if SummaryList is OK here - squeal if it isn't.
  
  float trigfrac; 
  
  if (fSummaryList->GetExpectedPoints()!=0) {
    // Do we have enough triggers?
    trigfrac = float(fSummaryList->GetTriggers())/float(fSummaryList->GetNumPulses());
  } else {
    delete fSummaryList;
    fSummaryList=0;
    return;
  }
  
  if (trigfrac < 0.95) {
    MSG("Pulser",Msg::kInfo)<<"Only saw "<<trigfrac<<" of expected "<<fSummaryList->GetNumPulses()<<" triggers. Dumping point.\n";
    delete fSummaryList;
    fSummaryList=0;
    return;
  }
  
  MSG("Pulser",Msg::kDebug) <<"PulserSummaryModule::EndRun(): trying to save the last point"<<endl;
  
  TFolder *lf = dynamic_cast<TFolder *>
    (gROOT->GetRootFolder()->FindObject("Loon"));
  if (lf==0) {
    lf = gROOT->GetRootFolder()->AddFolder("Loon", "Loon analysis");
    gROOT->GetListOfBrowsables()->Add(lf, "Loon");
  }

  TFolder *fTFolder = lf->AddFolder("Pulser","Pulser analysis");
  fTFolder->Add(fSummaryList);
  MSG("Pulser",Msg::kDebug) <<"PulserSummaryModule::EndRun(): last point added to a temporary folder: Pulser"<<endl;


}

---