TestCalibModule.cxx

Go to the documentation of this file.

#include <BeamDataUtil/test/TestCalibModule.h>
#include <BeamDataUtil/BDSwicCalibrator.h>
#include <BeamDataUtil/BDHadMuMon.h>
#include <BeamDataUtil/BDDevices.h>

#include <Conventions/Munits.h>
#include <DataUtil/GetRawBlock.h>

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


#include <JobControl/JobCModuleRegistry.h>
#include <JobControl/JobCResult.h>
#include <MessageService/MsgService.h>
CVSID("$Id: TestCalibModule.cxx,v 1.1 2005/05/09 22:12:26 minoscvs Exp $");
JOBMODULE(TestCalibModule,
          "BeamMonTestCalib", "Test calibrating the SWIC data.");

#include <vector>
#include <string>

using namespace std;
using namespace DataUtil;

TestCalibModule::TestCalibModule()
{
    for (int ind=0; ind<4; ++ind) {
        fHadMu[ind] = new BDHadMuMon;
    }
}
TestCalibModule::~TestCalibModule()
{
}

static void dump_hadmus(BDHadMuMon* hm[])
{    
    for (int ind=0; ind<4; ++ind) {
        double x,y,w,h;
        BDHadMuMon* hmm = hm[ind];

        hmm->GetStats(x,y,w,h);

        cerr << "\t" << ind << ": "
             << hmm->GetTotalCharge()
             << " (" << x<<", " << y << "),["<< w<<", "<<h<<"]\n";
        continue;
        int nrowcol = hmm->GetNrowcol();
        cerr << "pC:";
        for (int xind=0; xind<nrowcol; ++xind) {
            cerr << "\t\t";
            for (int yind=0; yind<nrowcol; ++yind) {
                double q = hmm->GetCharge(hmm->Index(hmm->Channel(xind+1,yind+1)));
                cerr << Form("%7.4e ",q/Munits::picocoloumb);
            }
            cerr << endl;
        }

    }
}

JobCResult TestCalibModule::Ana(const MomNavigator *mom)
{
    vector<const RawBeamMonBlock*> rbmbs = GetRawBlocks<RawBeamMonBlock>(mom);
    vector<const RawBeamMonHeaderBlock*> rbmhbs = GetRawBlocks<RawBeamMonHeaderBlock>(mom);

    if (rbmhbs.size() != 1 && rbmbs.size() != 1) {
        MSG("BD",Msg::kWarning)
            << "No beam monitoring data in mom\n";
        return JobCResult::kFailed;
    }

    const RawBeamMonHeaderBlock& rbmhb = *rbmhbs[0];
    const RawBeamMonBlock& rbmb = *rbmbs[0];

    vector<string> hadmu_name = BDDevices::HadMuMonitors();
    for (int ind=0; ind<4; ++ind) {
        const RawBeamData* d = rbmb[hadmu_name[ind].c_str()];
        if (!d) {
            MSG("BD",Msg::kDebug)
                << "No device data for " << hadmu_name[ind] << ", skipping\n";
            continue;
        }
        fHadMu[ind]->SetData(*d);
    }
    cerr << "Before calibration\n";
    dump_hadmus(fHadMu);

    BDSwicCalibrator::Get().Calibrate(rbmhb,rbmb);

    cerr << "After calibration\n";
    dump_hadmus(fHadMu);

    cerr << endl;

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

---