NtpBDFullModule.cxx

Go to the documentation of this file.

#include "BeamDataNtuple/Module/NtpBDFullModule.h"
#include "BeamDataNtuple/NtpBDFullRecord.h"
#include "BeamDataNtuple/BeamDataLiteHeader.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro

#include <RawData/RawRecord.h>
#include <RawData/RawBeamMonHeaderBlock.h>
#include <RawData/RawBeamMonBlock.h>

#include <Validity/VldContext.h>

#include <DataUtil/GetRecords.h>
#include <DataUtil/GetRawBlock.h>

#include <BeamDataUtil/BDScalar.h>

#include <vector>

ClassImp(NtpBDFullModule)

JOBMODULE(NtpBDFullModule, "NtpBDFull",
          "Fill the beam monitoring data ntuple");
CVSID("$Id: NtpBDFullModule.cxx,v 1.2 2007/03/01 17:38:35 rhatcher Exp $");
//......................................................................


NtpBDFullModule::NtpBDFullModule()
{}

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

NtpBDFullModule::~NtpBDFullModule()
{}

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

void NtpBDFullModule::BeginJob()
{
//======================================================================
// 
//======================================================================
    MSG("NtpBDFull",Msg::kDebug)<<"In NtpBDFull::BeginJob"<<endl;
}

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

void NtpBDFullModule::BeginFile()
{
//======================================================================
// Get the "run" and "subrun" number from the filename
//======================================================================
    MSG("NtpBDFull",Msg::kDebug)<<"In NtpBDFull::BeginFile"<<endl;
}

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

void NtpBDFullModule::EndJob()
{
//======================================================================
// Print out some summary at the end of the job
//======================================================================
    MSG("NtpBDFull",Msg::kDebug)<<"In NtpBDFull::EndJob"<<endl;
}

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


JobCResult NtpBDFullModule::Reco(MomNavigator* mom)
{
//======================================================================
// Read in the raw beam monitoring data files and fill the ntuples
//======================================================================
    MSG("NtpBDFull",Msg::kDebug)<<"In NtpBDFull::Reco"<<endl;

    JobCResult result(JobCResult::kPassed); // kNoDecision, kFailed, etc.

    vector<const RawBeamMonBlock*> rbmbs = DataUtil::GetRawBlocks<RawBeamMonBlock>(mom);
    vector<const RawBeamMonHeaderBlock*> rbmhbs = DataUtil::GetRawBlocks<RawBeamMonHeaderBlock>(mom);

    if (rbmhbs.size() != 1 && rbmbs.size() != 1) {
        MSG("BDFull",Msg::kWarning)
            << "No beam monitoring data in mom\n";
        return JobCResult::kFailed;
    }
    
    const RawBeamMonHeaderBlock& rbmhb = *rbmhbs[0];
    const RawBeamMonBlock& rbmb = *rbmbs[0];

    VldContext vc = rbmhb.GetVldContext();
    
    // make a header
    BeamDataLiteHeader bdh(vc);

    // TODO: fill the other members of the BeamDataLiteHeader
    
    NtpBDFullRecord* ntpbdr = new NtpBDFullRecord(bdh);

    NtpBDFullRecord& ntpbdr_r = *(ntpbdr);
    
    this->FillToroids(ntpbdr_r,rbmhb, rbmb);
        
    //give the ntpbdr object to mom, she'll write it to the file
    mom->AdoptFragment(ntpbdr);

    return result;

}

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


void NtpBDFullModule::FillToroids(NtpBDFullRecord& ntpbdr_r,
                                  const RawBeamMonHeaderBlock& rbmhb,
                                  const RawBeamMonBlock& rbmb)
{
    BDScalar bdpi1(5,10,"E:TRTGTD");
    bdpi1.SetSpill(rbmhb, rbmb);
    ntpbdr_r.protonIntensity = bdpi1.GetValue();
    MSG("NtpBDFull",Msg::kDebug) << "In NtpBDFullModule::FillToroids: PoT = "
                                 << ntpbdr_r.protonIntensity << endl;
}

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

---