BDEarliest.cxx

Go to the documentation of this file.

#include "BDEarliest.h"
#include "BDDevices.h"

#include <TString.h>            // for Form

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

#include <cmath>
using namespace std;

#include "MessageService/MsgService.h"
CVSID("$Id: BDEarliest.cxx,v 1.6 2005/11/07 23:59:31 bv Exp $");


BDEarliest::BDEarliest()
    : fDae(0), fVme(0), fLastDae(0), fLastVme(0), fBefore(0), fAfter(0)
{
    this->SetAllowedVmeDeviation();
}

static double find_earliest(const RawBeamMonBlock& rbmb, const vector<string>& dev)
{
    double earliest=-1;
    for (size_t ind=0; ind < dev.size(); ++ind) {
        const RawBeamData* rbd = rbmb[dev[ind]];
        if (!rbd) continue;

        double dae = rbd->GetSeconds() + rbd->GetMsecs()/1.0e3;
        if (earliest < 0 || dae < earliest)
            earliest = dae;
    }
    return earliest;
}
void BDEarliest::CalculateEarliest(const RawBeamMonBlock& rbmb,
                                   double& dae_out, double& vme_out,
                                   float before, float after)
{
    
    double earliest_vme=-1, earliest_dae=-1;

    RawBeamSwicData swic;

    vector<string> swic_names = BDDevices::SwicDevices();

    const float delay = rbmb.TclkTriggerDelay()/1000.0;
    MSG("BD",Msg::kVerbose)
        << "TCLK trigger delay is " << delay << endl;
    

    for (size_t ind=0; ind<swic_names.size(); ++ind) {
        const RawBeamData* rbd = rbmb[swic_names[ind]];
        if (!rbd) continue;

        double dae = rbd->GetSeconds() + rbd->GetMsecs()/1.0e3 - delay;

        if (earliest_dae < 0 || dae < earliest_dae)
            earliest_dae = dae;

        if (!swic.SetData(*rbd)) continue; // skip bogus data

        double vme = swic.VmeSeconds() + swic.VmeNanoseconds()/1.0e9;

        if (vme < 1.0) continue; // skip bogus timestamps

        double diff = dae-vme;

        if (diff < before || diff > after) {
            MSG("BD",Msg::kDebug)
                << swic_names[ind] << " has too large dae-vme=" << diff << endl;
            continue; // skip bogus timestamps
        }

        // Get here, then vme timestamp is non-bogus 

        if (earliest_vme < 0 || vme < earliest_vme) earliest_vme = vme;
    }

    if (earliest_vme < 0) earliest_vme = 0.0;
    if (earliest_dae < 0) {
        double tmp = find_earliest(rbmb,BDDevices::Toroids());
        if (tmp>1 && tmp<earliest_dae) earliest_dae = tmp;
    }
    if (earliest_dae < 0) {
        double tmp = find_earliest(rbmb,BDDevices::BpmPositions());
        if (tmp>1 && tmp<earliest_dae) earliest_dae = tmp;
    }
    if (earliest_dae < 0)
        earliest_dae = 0;
    
    dae_out = earliest_dae;
    vme_out = earliest_vme;

    MSG("BD",Msg::kVerbose)
        << Form("dae = %16.3f, vme = %16.3f, diff = %.3f",
                earliest_dae, earliest_vme, earliest_dae - earliest_vme)
        << endl;
}
    
void BDEarliest::SetSpill(const RawBeamMonHeaderBlock& /*rbmhb*/,
                          const RawBeamMonBlock& rbmb)
{
    fLastDae = fDae;
    fLastVme = fVme;
    
    fDae = fVme = 0.0;

    BDEarliest::CalculateEarliest(rbmb,fDae,fVme,fBefore,fAfter);
}

---