NtpBDLiteModule.cxx

Go to the documentation of this file.

#include "BeamDataNtuple/Module/NtpBDLiteModule.h"
#include "BeamDataNtuple/NtpBDLiteRecord.h"
#include "BeamDataNtuple/BeamDataLiteHeader.h"

#include <BeamDataUtil/BDSpillAccessor.h>

#include <Record/RecMinos.h>
#include <Record/RecMinosHdr.h>
#include <Record/RecRecord.h>
#include <RawData/RawRecord.h>

#include <MessageService/MsgService.h>
#include <MinosObjectMap/MomNavigator.h>

#include <Validity/VldContext.h>
#include <DataUtil/GetVldContext.h>

#include <SpillTiming/SpillTimeFinder.h>

#include <JobControl/JobCModuleRegistry.h> // For JOBMODULE macro

#include <vector>

ClassImp(NtpBDLiteModule)

JOBMODULE(NtpBDLiteModule, "NtpBDLite",
          "Fill the beam monitoring data ntuple");
CVSID("$Id: NtpBDLiteModule.cxx,v 1.2 2010/01/22 20:38:50 rhatcher Exp $");
//......................................................................


NtpBDLiteModule::NtpBDLiteModule()
    :fInputRecordsFilter("!ConfigRecord")
{}

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

NtpBDLiteModule::~NtpBDLiteModule()
{}

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

const Registry& NtpBDLiteModule::DefaultConfig() const
{
    MSG("NtpBDLite",Msg::kDebug) <<
        "NtpBDLiteModule::DefaultConfig" << endl;

    static Registry r;
    std::string name = this->JobCModule::GetName();
    name += ".config.default";
    r.SetName(name.c_str());

    r.UnLockValues();
    r.Set("InputRecordsFilter","!ConfigRecord");
    r.LockValues();

    return r;
}

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

void NtpBDLiteModule::Config(const Registry& r)
{
    MSG("NtpBDLite",Msg::kDebug) << "NtpBDLiteModule::Config" << endl;

    const char* tmps;
    if ( r.Get("InputRecordsFilter",tmps) ) fInputRecordsFilter = tmps;   
}

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

JobCResult NtpBDLiteModule::Reco(MomNavigator* mom)
{
    MSG("NtpBDLite",Msg::kDebug)<<"In NtpBDLite::Reco"<<endl;
    
    // try to get the vld context from somewhere

    VldContext vldc;
    
    RawRecord* daqsnarl =
      dynamic_cast<RawRecord *>(mom->GetFragment("RawRecord",0,"DaqSnarl"));
    if ( daqsnarl ) {
      // prefer a DaqSnarl record if one is available
      vldc = daqsnarl->GetHeader()->GetVldContext();
    } else {
      // no DaqSnarl
      if (fInputRecordsFilter.empty()) {
        MAXMSG("NtpBDLite",Msg::kWarning,5)
          << "No input record name set. Might pick up wrong VldContext." 
          << endl;
      }
    
      MSG("NtpBDLite",Msg::kDebug) << "Record filter = " 
                                   << fInputRecordsFilter << endl;

      vector<VldContext> vldcv = 
        DataUtil::GetVldContext(mom,fInputRecordsFilter.c_str());
      if ( vldcv.empty() ){
        MSG("NtpBDLite",Msg::kDebug)
          << "In NtpBDLite::Reco: "
          << "Nothing interesting found in mom." <<  endl;
        return JobCResult::kFailed;
      }

      vldc = vldcv[0];
      // Look for earliest time in the streams
      for (UInt_t i=1; i<vldcv.size(); ++i) {
        if (vldcv[i].GetTimeStamp()<vldc.GetTimeStamp()) vldc = vldcv[i] ;
      }
    }

    NtpBDLiteRecord* ntpbdr = FillNtpBDLite(vldc);

    //give the ntpbdr object to mom, she'll write it to the file
    mom->AdoptFragment(ntpbdr);
        
    return JobCResult::kPassed;

}

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

NtpBDLiteRecord* NtpBDLiteModule::FillNtpBDLite(const VldContext& vldc)
{

    MSG("NtpBDLite",Msg::kDebug)<< "vldc = "  << vldc <<endl;
    BeamDataLiteHeader bdh(vldc);

    const BeamMonSpill* spill = BDSpillAccessor::Get().LoadSpill(vldc.GetTimeStamp());
    //BDSpillAccessor B;
    //const BeamMonSpill* spill = B.LoadSpill(vldc.GetTimeStamp());
    if (!spill){ 
        MSG("NtpBDLite",Msg::kDebug)<< "No spill found!" <<endl;
        NtpBDLiteRecord* ntpbdr = new NtpBDLiteRecord(bdh);
        return ntpbdr;
    }
    
    // fill members of the BeamDataLiteHeader
    
    VldTimeStamp spill_time = spill->SpillTime();
    bdh.SetEarliestTimeStamp(spill_time);
    VldTimeStamp vldd = spill_time-vldc.GetTimeStamp();
    bdh.SetTimeDiffStreamSpill(vldd.GetSeconds()); 
    bdh.SetFoundBD(1);
    bdh.SetStatus(spill->GetStatusInt());
    
    NtpBDLiteRecord* ntpbdr = new NtpBDLiteRecord(bdh);
    MSG("NtpBDLite",Msg::kDebug)<< spill->fTor101  <<endl;
    ntpbdr->tor101 = spill->fTor101;           
    ntpbdr->tr101d = spill->fTr101d;           
    ntpbdr->tortgt = spill->fTortgt;           
    ntpbdr->trtgtd = spill->fTrtgtd;
    ntpbdr->horncur = spill->fHornCur;     
    for (Int_t i=0;i<6;++i){
        ntpbdr->bposx[i]=spill->fTargBpmX[i];
        ntpbdr->bposy[i]=spill->fTargBpmY[i];
        ntpbdr->bpmint[i]=spill->fBpmInt[i];
    }
    ntpbdr->bwidx = spill->fProfWidX;  
    ntpbdr->bwidy = spill->fProfWidY;  
    ntpbdr->hadint = spill->fHadInt; 
    ntpbdr->muint1 = spill->fMuInt1;  
    ntpbdr->muint2 = spill->fMuInt2;  
    ntpbdr->muint3 = spill->fMuInt3;  

    SpillTimeFinder &stf = SpillTimeFinder::Instance();
    ntpbdr->dt_nearest = stf.GetTimeToNearestSpill(vldc);
    VldTimeStamp tons = stf.GetTimeOfNearestSpill(vldc);
    ntpbdr->nearest_sec = tons.GetSec();
    ntpbdr->nearest_nsec = tons.GetNanoSec();

    return ntpbdr;
}

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

---