BeamMonSwicPedsDbuModule.cxx

Go to the documentation of this file.

#include "BeamMonSwicPedsDbuModule.h"

#include <BeamDataUtil/BeamMonSwicPeds.h>
#include <BeamDataUtil/BDUniquify.h>

#include <BeamDataUtil/BDSwicPeds.h>
#include <BeamDataUtil/BDPipeline.h>

#include <MessageService/MsgService.h>
#include <JobControl/JobCModuleRegistry.h>

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

CVSID("$Id: BeamMonSwicPedsDbuModule.cxx,v 1.11 2006/04/19 20:58:05 rhatcher Exp $");

#include <DatabaseInterface/DbiWriter.h>
//  Instantiate associated Result Pointer class.
#include <DatabaseInterface/DbiWriter.tpl>
template class  DbiWriter<BeamMonSwicPeds>;

#include <Registry/Registry.h>

#include <DataUtil/GetRawBlock.h>
using namespace DataUtil;


JOBMODULE(BeamMonSwicPedsDbuModule,
          "BeamMonSwicPedsDbu", 
          "Update the database with Beam Monitoring device pedestals.");

using namespace std;

BeamMonSwicPedsDbuModule::BeamMonSwicPedsDbuModule()
    : fPipe()
    , fPedsProc(new BDSwicPeds)
    , fCount(0)
{
    fPipe.GetProcessPipeline().push_back(fPedsProc);
}

BeamMonSwicPedsDbuModule::~BeamMonSwicPedsDbuModule()
{
}

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

    size_t nblocks = blocks.size();

    if (nblocks == 0) {
        MSG("BD",Msg::kWarning)
            << "No good beam monitoring data in MOM\n";
        return JobCResult::kFailed;
    }

    const int trigger_event = 0xa9; // FIXME: needs to be a configurable param.
    for (size_t iblock=0; iblock<nblocks; ++iblock) {
        const RawBeamMonHeaderBlock& rbmhb = *blocks[iblock].first;
        const RawBeamMonBlock& rbmb = *blocks[iblock].second;

        if (rbmb.TclkTriggerEvent() != trigger_event) {
            MSG("BD",Msg::kWarning)
                << "Wrong trigger event: " << rbmb.TclkTriggerEvent()
                << " != " << trigger_event << endl;
            continue;
        }


        fPipe.SetSpill(rbmhb,rbmb);
    }

    return JobCResult::kAOK;
}
void BeamMonSwicPedsDbuModule::BeginFile()
{
    MSG("BD",Msg::kDebug) << "BeginFile: " << this->GetCurrentFile() << endl;
}


void BeamMonSwicPedsDbuModule::EndFile()
{
    MSG("BD",Msg::kDebug) << "EOF\n";
    this->Dump();
}

void BeamMonSwicPedsDbuModule::EndJob()
{
    MSG("BD",Msg::kDebug) << "EOJ\n";
}

void BeamMonSwicPedsDbuModule::Dump()
{
    VldTimeStamp vts_start = fPedsProc->GetStart();
    VldTimeStamp vts_stop = fPedsProc->GetStop();

    if (vts_stop < vts_start) {
        MSG("BD",Msg::kWarning)
            << "No Peds for any devices in this database dump\n";
        return;
    }

    const BDSwicPeds::SwicMap& chmap = fPedsProc->GetAllChannels();
    BDSwicPeds::SwicMap::const_iterator it, done = chmap.end();
    
    // Do an initial run through to see if there is at least one
    // pedestal set ready.  This is done to avoid creating VLD entries
    // with no associated main table rows.
    bool got_one = false;
    for (it=chmap.begin(); it != done; ++it) {
        assert(it->second);
        if (it->second->GetPopulationSize() > 0) {
            got_one = true;
            break;
        }
    }
    if (!got_one) {
        MSG("BD",Msg::kWarning) << "Aborting DB dump, no data\n";
        return;
    }

    VldRange range(Detector::kNear|Detector::kFar,SimFlag::kData,
                   vts_start,vts_stop,"Beam");
    
    VldTimeStamp now;
    DbiWriter<BeamMonSwicPeds> writer(range,-1,0,now);


    for (it=chmap.begin(); it != done; ++it) {
        const UtilRunningAverageVector& urav = *(it->second);

        int nsamples = urav.GetPopulationSize();
        if (!nsamples) continue;

        BeamMonSwicPeds ped(it->first.c_str());
        ped.SetNsamples(nsamples);

        vector<double> means = urav.PopulationMean();

        if (means.size() != 96) {
            MSG("BD",Msg::kFatal)
                << nsamples << " samples but size (" << means.size()
                << ") doesn't match for " << it->first << endl;
            abort();
        }

        double tot = 0;
        for (int ind=0; ind<96; ++ind) tot += means[ind];

        ped.SetChannelMeans(means);
        ped.SetChannelSigmas(urav.PopulationSigma());

        writer << ped;
        MSG("BD",Msg::kDebug)
            << it->first  << ", " << nsamples << " samples, tot peds = " << tot << endl;
        if (tot == 0.0) {
            MSG("BD",Msg::kWarning)
                << it->first  << ", " << nsamples << " samples, tot peds = " << tot << endl;
        }
            
    }

    writer.Close();
    fPedsProc->Clear();
}

---