Cuts.cxx

Go to the documentation of this file.

#include <cmath>

#include "Cuts.h"

#include <HistMan/HistMan.h>

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

#include <TH1D.h>
#include <TH2D.h>

void Cuts::Book(PlotterManager& )
{
    HistMan hm("cuts");
    hm.Book<TH1D>("tortgt","Proton intensity at E:TORTGT (x10^12 protons)",
                  100,0,10);
    hm.Book<TH1D>("tortgt_time","Time of E:TORTGT (from start of run 6067) (sec)",
                  108,-21600,21600);
    hm.Book<TH1D>("horncur","Horn Current, all quadrants (KAmp)",
                  40,-200,200);
}

// Check if timestamp is in the two runs we care about
static bool tortgt_ok(PlotterManager& pm)
{
    const int ti_6067 = 1106355930-1; // 2005-01-22 01:05:30 UCT
    const int tf_6067 = 1106362789+1; // 2005-01-22 02:59:49 UTC
    const int ti_6068 = 1106409473-1; // 2005-01-22 15:57:53 UTC
    const int tf_6068 = 1106420416+1; // 2005-01-22 19:00:16 UTC

    const AcnetDevice* tortgt = pm.GetKnot().GetDevice("E:TORTGT");
    if (!tortgt->ndata) return false;

    HistMan hm("cuts");

    hm.Fill1d("tortgt",tortgt->data[0]);
    if (tortgt->data[0] < 2.0) return false;

    int sec = tortgt->sec;
    hm.Fill1d("tortgt_time",sec-ti_6067+1);

    bool in_runs = ((sec>=ti_6067 && sec<=tf_6067) || (sec>=ti_6068 && sec<=tf_6068));
    return in_runs;
}

// Check if the horn current is decent
static bool horn_off(PlotterManager& pm)
{
    const char* hc[] = {
        "E:NSLINA",
        "E:NSLINB",
        "E:NSLINC",
        "E:NSLIND",
        0
    };
    double horn_current = 0;
    for (int ind = 0; hc[ind]; ++ind) {
        const AcnetDevice* ad = pm.GetKnot().GetDevice(hc[ind]);
        if (!ad->ndata) continue;
        horn_current += ad->data[0];
    }
    HistMan hm("cuts");
    hm.Fill1d("horncur",horn_current);
    return fabs(horn_current) < 175;
}

bool Cuts::Fill(PlotterManager& pm)
{
    if (!tortgt_ok(pm)) return false;

    if (horn_off(pm)) return false;
     
    return true;
}

---