BmnHistMSC.cxx

Go to the documentation of this file.

#include "BmnHistMSC.h"
 
#include <algorithm>
#include <numeric>
// BmnRoot
#include <BmnMSCZSDigit.h>
 
BmnHistMSC::BmnHistMSC(TString title, TString path, Int_t PeriodID, BmnSetup stp, BmnTrigRaw2Digit* trigMapper)
    : BmnHist(PeriodID, stp)
    , fTrigMapper(trigMapper)
{
    sumMods = 0;
    maxLayers = 0;
    refPath = path;
    fTitle = title;
    fName = title + "_cl";
 
    scalers2index_map = fTrigMapper->GetScaler2IdxMap();
    vector<string>& scaler_idx2name = fTrigMapper->GetScalerIdx2NameMap();
    uint16_t scalers_ar_size = scaler_idx2name.size();
 
    TString name;
    name = fTitle + "_Scalers_Total";
    TH1I* h = new TH1I(name, name, scalers_ar_size, 0, scalers_ar_size);
    h->GetXaxis()->SetTitle("Scaler");
    h->GetYaxis()->SetTitle("Activation Count");
    TAxis* xaTotal = h->GetXaxis();
    for (uint16_t i = 0; i < scalers_ar_size; i++) {
        string& sname = scaler_idx2name[i];
        xaTotal->SetBinLabel(i + 1, sname.data());
    }
    BmnHist::SetHistStyleTH1(h);
    h->SetLineWidth(2);
    histMscAll.push_back(h);
    maxLayers = 1;
    sumMods = 1;
    // Create canvas
    name = fTitle + "Canvas";
    canGemStrip = new TCanvas(name, name, PAD_WIDTH * maxLayers, PAD_HEIGHT * sumMods);
    canGemStrip->Divide(maxLayers, sumMods, 0.01 / maxLayers, 0.01 / sumMods);
    //    Int_t modCtr = 0; // filling GEM Canvas' pads
    canStripPads.resize(sumMods * maxLayers);
    Names.resize(sumMods * maxLayers);
    PadInfo* pad = new PadInfo();
    pad->opt = "hist";
    pad->current = h;
    pad->logy = true;
    Int_t iPad = 0;   // modCtr * maxLayers + iLayer;
    canStripPads[iPad] = pad;
    canGemStrip->GetPad(iPad + 1)->SetGrid();
    Names[iPad] = canStripPads[iPad]->current->GetName();
}
 
BmnHistMSC::~BmnHistMSC()
{
    delete canGemStrip;
    if (fDir)
        return;
    for (auto pad : canStripPads)
        delete pad;
}
 
void BmnHistMSC::Register(THttpServer* serv)
{
    fServer = serv;
    fServer->Register("/", this);
    TString path = "/" + fTitle + "/";
    fServer->Register(path, canGemStrip);
    TString cmd = "/" + fName + "/->Reset()";
    TString cmdTitle = path + "Reset";
    fServer->RegisterCommand(cmdTitle.Data(), cmd.Data(), "button;");
    //    fServer->Restrict(cmdTitle.Data(), "visible=shift");
    //    fServer->Restrict(cmdTitle.Data(), "allow=shift");
    //    fServer->Restrict(cmdTitle.Data(), "deny=guest");
    cmd = "/" + fName + "/->SetRefRun(%arg1%)";
    cmdTitle = path + "SetRefRun";
    fServer->RegisterCommand(cmdTitle.Data(), cmd.Data(), "button;");
}
 
void BmnHistMSC::SetDir(TFile* outFile, TTree* recoTree)
{
    frecoTree = recoTree;
    fDir = NULL;
    if (outFile != NULL)
        fDir = outFile->mkdir(fTitle + "_hists");
    for (auto el : histMscAll)
        el->SetDirectory(fDir);
}
 
void BmnHistMSC::DrawBoth()
{
    BmnHist::DrawRef(canGemStrip, &canStripPads);
}
 
void BmnHistMSC::FillFromDigi(DigiArrays* fDigiArrays)
{
    TClonesArray* digs = fDigiArrays->msc;
    //            printf("evtype %u\n", fDigiArrays->header->GetEventType());
    //            printf("entries %d\n", digs->GetEntriesFast());
    if (!digs)
        return;
    if (digs->GetClass() == BmnMSCZSDigit<uint8_t>::Class()) {
        for (Int_t digIndex = 0; digIndex < digs->GetEntriesFast(); digIndex++) {
            BmnMSCZSDigit<uint8_t>* gs = static_cast<BmnMSCZSDigit<uint8_t>*>(digs->At(digIndex));
            auto it = scalers2index_map.find(make_pair(gs->GetSerial(), gs->GetChannel()));
            if (it == scalers2index_map.end())
                continue;
            vector<uint8_t>& vec = gs->GetValues();
            uint64_t sum = std::accumulate(vec.begin(), vec.end(), 0);
            //            printf("MSC[%2u] = %lu\n", it->second, sum);
            histMscAll[0]->Fill(it->second, sum);
        }
    }
}
 
BmnStatus BmnHistMSC::SetRefRun(Int_t id)
{
    if (refID != id) {
        refID = id;
        refRunName = Form("bmn_run%04u%s_hist.root", refID, fHistFileNameAdd.Data());
        BmnHist::LoadRefRun(refID, refPath + refRunName, fTitle, canStripPads, Names);
        DrawBoth();
    }
 
    return kBMNSUCCESS;
}
 
void BmnHistMSC::ClearRefRun()
{
    for (auto pad : canStripPads) {
        if (pad)
            if (pad->ref) {
                delete pad->ref;
                pad->ref = nullptr;
            }
    }
    refID = 0;
}
 
void BmnHistMSC::Reset()
{
    for (auto el : histMscAll)
        el->Reset();
}