DeltaT.cxx

Go to the documentation of this file.

#include "DeltaT.h"

#include <HistMan/HistMan.h>

#include <BeamData/ana/bmnt/Knot.h>
#include <BeamData/ana/bmnt/Swic.h>
#include <BeamData/ana/bmnt/AcnetDevice.h>

#include <TTree.h>
#include <TFile.h>
#include <TH2D.h>

#include <cmath>
#include <cassert>
#include <vector>
#include <iostream>
using namespace std;

DeltaT::DeltaT()
    : fFile(0)
{
}

DeltaT::~DeltaT()
{
    if (fFile) delete fFile; fFile = 0;
}

void DeltaT::NewFile(const char* filename)
{
    if (fFile) delete fFile; fFile = 0;
    fFile = new TFile(filename);

    TObject* o = fFile->FindObjectAny("bd");
    TTree* bd = dynamic_cast<TTree*>(o);
    assert(bd);
    fKnot.SetTree(*bd);
}

void DeltaT::Book(PlotterManager& pm)
{

    vector<const char*> devs = pm.GetSwicDevNames();
    size_t siz = devs.size();

    HistMan& hm = pm.GetHistMan();

    TH2D* devhist[3];

    int i=0;
    devhist[i++] = hm.Book<TH2D>("dae_vme_dt","DAE-VME time difference",
                                 siz,0,siz, 400,0,4);
    devhist[i++] = hm.Book<TH2D>("dae_dt_prev","DAE time to previous spill",
                                 siz,0,siz, 240,0,240);
    devhist[i++] = hm.Book<TH2D>("vme_dt_prev","VME time to previous spill",
                                 siz,0,siz, 240,0,240);
    devhist[i++] = hm.Book<TH2D>("dae_vme_dt_prev","DAE-VME-(DAE time to prev spill)",
                                 siz,0,siz, 400,0,4);

    TH1D* missing =
        hm.Book<TH1D>("missing","Number of times the device is missing",
                      siz,0,siz);

    hm.Book<TH1D>("earliest_dae_vme_dt","Delta between earliest DAE and VME",
                  1000,-5,5);

    for (size_t bin=0; bin<siz; ++bin) {
        for (int ind=0; ind<i; ++ind) {
            devhist[ind]->GetXaxis()->SetBinLabel(bin+1,devs[bin]);
        }
        missing->GetXaxis()->SetBinLabel(bin+1,devs[bin]);
    }


}
bool DeltaT::Fill(PlotterManager& pm)
{
    Knot& k = pm.GetKnot();
    HistMan& hm = pm.GetHistMan();

    vector<const char*> devs = pm.GetSwicDevNames();
    size_t siz = devs.size();
    
    int entry = k.GetEntryNumber();
    if (entry > 0)  fKnot.GetEntry(entry-1);

    double early_dae = -1.0, early_vme = -1.0;

    Swic swic1, swic2;
    for (size_t dev=0; dev<siz; ++dev) {
        const AcnetDevice* ad1 = k.GetDevice(devs[dev]);
        if (!ad1->ndata) {
            hm.Fill1d("missing",dev,1);
            continue;
        }

        swic1.SetAcnetDevice(*ad1);

        double diff0 = ad1->timestamp-swic1.GetVmeTimestamp();
        if (fabs(diff0)<1e-6) diff0=0.0;
        hm.Fill2d("dae_vme_dt",dev,diff0);

        // Track earliest DAE & VME timestamp.
        if (ad1->timestamp > 1.0)
            if (early_dae < 0.0 || ad1->timestamp < early_dae)
                early_dae = ad1->timestamp;
        double vme_timestamp = swic1.GetVmeTimestamp();
        if (vme_timestamp > 1.0)
            if (early_vme < 0.0 || vme_timestamp < early_vme)
                early_vme = vme_timestamp;
        
        // Need to pass the first spill before comparing to previous....
        if (entry < 1) continue;

        const AcnetDevice* ad2 = fKnot.GetDevice(devs[dev]);
        if (!ad2->ndata) continue;

        swic2.SetAcnetDevice(*ad2);

        double diff1 = ad1->timestamp-ad2->timestamp;
        if (fabs(diff1)<1e-6) diff1=0.0;
        hm.Fill2d("dae_dt_prev",dev,diff1);

        double diff2 = swic1.GetVmeTimestamp()-swic2.GetVmeTimestamp();
        if (fabs(diff2)<1e-6) diff2=0.0;
        hm.Fill2d("vme_dt_prev",dev,diff2);

        double diff3 = diff0-diff1;
        if (diff0 < 0 || diff0 > 4) {
            hm.Fill2d("dae_vme_dt_prev",dev,diff3);
        }
    }
    hm.Fill1d("earliest_dae_vme_dt",early_dae-early_vme);
    return true;
}

---