PulserFromRawModule.cxx

Go to the documentation of this file.

#include <cstdio>
#include <fstream>
// ROOT includes
#include "TFile.h"
// MINOS includes
#include "MessageService/MsgService.h"     // MSG text output
#include "MinosObjectMap/MomNavigator.h"   // Data access
#include "CandData/CandRecord.h"           // DigitList handling
#include "JobControl/JobCommand.h"         // JobCommand handling
#include "JobControl/JobCModuleRegistry.h" // JOBMODULE registration macro
#include "RawData/RawDaqSnarlHeader.h"
#include "RawData/RawSnarlHeaderBlock.h"
#include "RawData/RawRecord.h"
#include "RawData/RawDigitDataBlock.h"
#include "RawData/RawLIAdcSummaryBlock.h"
#include "PulserCalibration/PulserSummaryList.h"
#include "PulserCalibration/PulserSummary.h"
#include "PulserCalibration/PulserFromRawModule.h"
#include "Plex/PlexStripEndId.h"
#include "Plex/PlexHandle.h"
#include "RawData/RawChannelId.h"
#include "RawData/RawDigit.h"
#include "Validity/VldContext.h"
#include "Validity/VldRange.h"
#include "Conventions/StripEnd.h"
#include "DatabaseInterface/DbiWriter.h"
#include "Plex/PlexStripEndToLed.h"
#include "Plex/PlexPinDiodeToLed.h"
//include "PEGain/PEGainFitter.h"
ClassImp(PulserFromRawModule)

CVSID("$Id: PulserFromRawModule.cxx,v 1.5 2006/04/19 20:00:55 rhatcher Exp $");
JOBMODULE(PulserFromRawModule,"PulserFromRawModule",
          "Make Pulser Summaries From Raw Data\n");

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

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

PulserFromRawModule::PulserFromRawModule()
   //:fDoingThisSum(0)
{
//======================================================================
// Set the initial state of your module. The code above initializes the
// variables listed to the values in parentheses. In this case that's
// all I need to do so the method itself is empty...
//======================================================================
    fDoingThisSum=0;
    fSummaryList=0;

}

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

JobCResult PulserFromRawModule::Reco(MomNavigator *mom) 
{
   JobCResult result(JobCResult::kPassed);
   
   RawRecord *rawrec = dynamic_cast<RawRecord *>(mom->GetFragment("RawRecord"));
   if (rawrec == 0) {
      MSG("Pulser", Msg::kWarning) << "No RawRecord in MOM." << endl;
      return result;
   }

   MSG("Pulser",Msg::kVerbose) << "Reco() got vld: " 
                             << (rawrec->GetRawHeader())->GetVldContext() 
                             << endl;
   int sec = rawrec->GetHeader()->GetVldContext().GetTimeStamp().GetSec();
   UInt_t trigword = 0;
   UInt_t frun = 0;
   UInt_t fsnarl = 0;
   

   const RawDaqSnarlHeader* snarlHdr =
      dynamic_cast<const RawDaqSnarlHeader*>(rawrec->GetRawHeader());
   if (snarlHdr) {
      frun = snarlHdr->GetRun();
      fsnarl = snarlHdr->GetSnarl();
      trigword = snarlHdr->GetTrigSrc();
   }
   const RawLIAdcSummaryBlock *rawsum = 
      dynamic_cast<const RawLIAdcSummaryBlock *>
      (rawrec->FindRawBlock("RawLIAdcSummaryBlock"));   
   if (rawsum) {
      MSG("Pulser",Msg::kDebug) 
         << "Got LI Summary for Led :" 
         << rawsum->GetLed()
         << " H " << rawsum->GetPulseHeight() << endl;
      if(!fDoingThisSum) {
         fDoingThisSum=new PulserFromRawLiBlock(*rawsum);
      }
      else fDoingThisSum->New(*rawsum);
//      cerr << "Hello: " << fDoingThisSum << " "
//                 << fDoingThisSum->GetPulses() << endl;
   }
  //  cerr << "Here " << fDoingThisSum << " "
//      << fDoingThisSum->GetPulses() << endl;
   const RawDigitDataBlock *rawddb =
      dynamic_cast<const RawDigitDataBlock *>
      (rawrec->FindRawBlock("RawDigitDataBlock"));
   
   if (rawddb && fDoingThisSum) {
//       cerr << "Here also " << fDoingThisSum << endl;
//       cerr <<  fDoingThisSum->GetPulses() << endl;
      if (fSummaryList) {
         if (fDoingThisSum->GetLed()==fSummaryList->GetLed() && 
             fDoingThisSum->GetCalibPoint()==fSummaryList->GetPointIndex()) {
//           cerr << "Here !!!" << endl;
            fSummaryList->Add(*rawddb,*fDoingThisSum);
         }
         else {
            // We're on a new point now!
//           cerr << "Here !!!!" << endl;
            MSG("Pulser",Msg::kDebug)<< "New point at " << sec <<endl;
            FinishList(mom);
            fSummaryList = new PulserFromRawSummaryList(*rawddb,*fDoingThisSum);
         }
      } 
      else { 
//        cerr << "Here !!!!! " << fDoingThisSum->GetPulses() << endl;
         
         fSummaryList = new PulserFromRawSummaryList(*rawddb,*fDoingThisSum);
      }
   }
   
   return result;
   
}

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

void PulserFromRawModule::EndJob() 
{


}

void PulserFromRawModule::FinishList(MomNavigator* mom)
{
    MSG("Pulser",Msg::kInfo) <<"Summary Complete! Box " << fSummaryList->GetPulserBox() << " Led "<< fSummaryList->GetLed()<<", point "<<fSummaryList->GetPointIndex()<<" of "<<fSummaryList->GetExpectedPoints()<<endl;
   fSummaryList->Finish();
   mom->AdoptFragment(fSummaryList);
   fSummaryList=0;
}

---