BmnHistProf.cxx

Go to the documentation of this file.

#include "BmnHistProf.h"
 
#include "BmnSiProfDigit.h"
#include "BmnSiProfStationSet.h"
 
BmnHistProf::BmnHistProf(TString title, TString path, Int_t PeriodID, BmnSetup setup)
    : BmnHist(PeriodID, setup)
{
    sumMods = 0;
    maxLayers = 0;
    refPath = path;
    fTitle = title;
    fName = title + "_cl";
    TString name;
    unique_ptr<BmnSiProfStationSet> stationSet = BmnSiProfStationSet::Create(PeriodID, fSetup);
    for (Int_t iStation = 0; iStation < stationSet->GetNStations(); iStation++) {
        vector<vector<TH1F*>> rowGEM;
        auto* st = stationSet->GetStation(iStation);
        sumMods += st->GetNModules();
        for (Int_t iModule = 0; iModule < st->GetNModules(); iModule++) {
            vector<TH1F*> colGEM;
            auto* mod = st->GetModule(iModule);
            if (maxLayers < mod->GetNStripLayers())
                maxLayers = mod->GetNStripLayers();
            for (Int_t iLayer = 0; iLayer < mod->GetNStripLayers(); iLayer++) {
                auto lay = mod->GetStripLayer(iLayer);
                name = Form(fTitle + "_Station_%d_module_%d_layer_%d", iStation, iModule, iLayer);
                TH1F* h = new TH1F(name, name, lay.GetNStrips(), 0, lay.GetNStrips());
                h->SetTitleSize(0.06, "XY");
                h->SetLabelSize(0.08, "XY");
                h->GetXaxis()->SetTitle("Strip Number");
                h->GetXaxis()->SetTitleColor(kOrange + 10);
                h->GetYaxis()->SetTitle("Activation Count");
                h->GetYaxis()->SetTitleColor(kOrange + 10);
                colGEM.push_back(h);
            }
            rowGEM.push_back(colGEM);
        }
        histSiliconStrip.push_back(rowGEM);
    }
    maxLayers = 2;
    // Create canvas
    name = fTitle + "Canvas";
    canStrip = new TCanvas(name, name, PAD_WIDTH * maxLayers, PAD_HEIGHT * sumMods);
    canStrip->Divide(maxLayers, sumMods);
    Int_t modCtr = 0;   // filling GEM Canvas' pads
    canStripPads.resize(sumMods * maxLayers);
    for (auto& pad : canStripPads) {
        pad = nullptr;
    }
    Names.resize(sumMods * maxLayers);
    for (Int_t iStation = 0; iStation < stationSet->GetNStations(); iStation++) {
        auto* st = stationSet->GetStation(iStation);
        for (Int_t iModule = 0; iModule < st->GetNModules(); iModule++) {
            auto* mod = st->GetModule(iModule);
            for (Int_t iLayer = 0; iLayer < 2 /*mod->GetNStripLayers()*/; iLayer++) {
                Int_t iPad = modCtr * maxLayers + iLayer;
                PadInfo* p = new PadInfo();
                p->current = histSiliconStrip[iStation][iModule][iLayer];
                canStripPads[iPad] = p;
                canStrip->GetPad(iPad + 1)->SetGrid();
                Names[iPad] = canStripPads[iPad]->current->GetName();
            }
            modCtr++;
        }
    }
}
 
BmnHistProf::~BmnHistProf()
{
    //    delete gemStationSet;
    delete canStrip;
    if (fDir)
        return;
    for (auto pad : canStripPads)
        delete pad;
}
 
void BmnHistProf::Register(THttpServer* serv)
{
    fServer = serv;
    fServer->Register("/", this);
    TString path = "/" + fTitle + "/";
    fServer->Register(path, canStrip);
    fServer->SetItemField(path, "_monitoring", "2000");
    fServer->SetItemField(path, "_layout", "grid3x3");
    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 BmnHistProf::SetDir(TFile* outFile, TTree* recoTree)
{
    frecoTree = recoTree;
    fDir = NULL;
    if (outFile != NULL)
        fDir = outFile->mkdir(fTitle + "_hists");
    for (auto row : histSiliconStrip)
        for (auto col : row)
            for (auto el : col)
                el->SetDirectory(fDir);
}
 
void BmnHistProf::DrawBoth()
{
    BmnHist::DrawRef(canStrip, &canStripPads);
}
 
void BmnHistProf::FillFromDigi(DigiArrays* fDigiArrays)
{
    TClonesArray* gemDigits = fDigiArrays->siprof;
    if (!gemDigits)
        return;
    for (Int_t digIndex = 0; digIndex < gemDigits->GetEntriesFast(); digIndex++) {
        BmnSiProfDigit* gs = (BmnSiProfDigit*)gemDigits->At(digIndex);
        Int_t module = gs->GetModule();
        Int_t station = gs->GetStation();
        Int_t layer = gs->GetStripLayer();
        Int_t gemStrip = gs->GetStripNumber();
        histSiliconStrip[station][module][layer]->Fill(gemStrip);
    }
}
 
BmnStatus BmnHistProf::SetRefRun(Int_t id)
{
    if (refID != id) {
        TString FileName = Form("bmn_run%04d_hist.root", id);
        printf("SetRefRun: %s\n", FileName.Data());
        refRunName = FileName;
        refID = id;
        BmnHist::LoadRefRun(refID, refPath + FileName, fTitle, canStripPads, Names);
        DrawBoth();
    }
 
    return kBMNSUCCESS;
}
 
void BmnHistProf::ClearRefRun()
{
    for (auto& pad : canStripPads) {
        if (pad) {
            if (pad->ref)
                delete pad->ref;
            pad->ref = NULL;
        }
    }
    refID = 0;
}
 
void BmnHistProf::Reset()
{
    for (auto& row : histSiliconStrip)
        for (auto& col : row)
            for (auto& el : col)
                el->Reset();
}