BDUniquifyTestModule.cxx

Go to the documentation of this file.

#include "BDUniquifyTestModule.h"
#include <BeamDataUtil/BDUniquify.h>
#include <BeamDataUtil/BDEarliest.h>

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

using namespace std;

#include <MinosObjectMap/MomNavigator.h>
#include <JobControl/JobCModuleRegistry.h>
#include <JobControl/JobCResult.h>
#include <MessageService/MsgService.h>
CVSID("$Id: BDUniquifyTestModule.cxx,v 1.2 2005/11/07 23:59:31 bv Exp $");
JOBMODULE(BDUniquifyTestModule,
          "BDUniquifyTest", "Test BDUniquify class.");


BDUniquifyTestModule::BDUniquifyTestModule()
{
}
BDUniquifyTestModule::~BDUniquifyTestModule()
{
}
static double get_scalar(const RawBeamMonBlock& rbmb, const char* name)
{
    const RawBeamData* data = rbmb[name];
    if (!data) return 0.0;
    const double* ddata = data->GetData();
    if (!ddata) return 0.0;
    return ddata[0];
}


static void dump(const char* prefix,
                 const RawBeamMonHeaderBlock& rbmhb,
                 const RawBeamMonBlock& rbmb)
{
    static BDEarliest earliest;
    earliest.SetSpill(rbmhb,rbmb);
    double dae=0,vme=0;
    earliest.GetTimestamps(dae, vme);
    VldTimeStamp vts_dae(dae), vts_vme(vme);

    double tortgt = get_scalar(rbmb,"E:TORTGT");
    double tor101 = get_scalar(rbmb,"E:TOR101");

    int bad = 0;
    vector<string> names = rbmb.GetNames();
    size_t siz = names.size();
    VldTimeStamp vts = rbmhb.GetTimeStamp();
    for (size_t ind=0; ind<siz; ++ind) {
        const RawBeamData* rbd = rbmb[names[ind]];
        if (!rbd) continue;

        VldTimeStamp tmp(rbd->GetSeconds(),rbd->GetMsecs()*1000000);
        if (fabs(vts - tmp) > 0.5) ++bad;
    }

    cerr << prefix << rbmhb.GetTimeStamp() << " "
         << "dae: " << vts_dae << " "
         << "vme: " << vts_vme << " "
         << bad << "/" << siz
         << " "
         << rbmhb.GetSpillCountNum()
         << " "
         << tor101 << " " << tortgt
         << endl;
}

JobCResult BDUniquifyTestModule::Ana(const MomNavigator *mom)
{
    static BDUniquify uniqifier;
    vector<BDUniquify::BlockPair> blocks = uniqifier.GetUniqueBlocks(*mom);

    cerr << "Got " << blocks.size() << " unique blocks\n";
    for (size_t ind=0; ind<blocks.size(); ++ind) {
        dump("u:\t",*blocks[ind].first,*blocks[ind].second);
    }


    // Loop through all RawRecords in mom:
    TObject* tobj=0;
    TIter fragiter = mom->FragmentIter();
    while( ( tobj = fragiter.Next() ) ) {
        const RawRecord* rawrec = dynamic_cast<const RawRecord*>(tobj);
        if (!rawrec) continue;

        const RawBeamMonHeaderBlock* rbmhb=0;
        const RawBeamMonBlock* rbmb=0;

        TIter blockiter = rawrec->GetRawBlockIter();
        while( ( tobj = blockiter.Next() ) ) {
            if (!rbmhb)
                rbmhb = dynamic_cast<const RawBeamMonHeaderBlock*>(tobj);
            if (!rbmb)
                rbmb = dynamic_cast<const RawBeamMonBlock*>(tobj);
        }

        if (!rbmhb || !rbmb) continue;
        dump("a:\t",*rbmhb,*rbmb);
    }


    return JobCResult::kAOK;
}
void BDUniquifyTestModule::BeginJob()
{
}
void BDUniquifyTestModule::EndJob()
{
}

---