tango_gem_hv.C

Go to the documentation of this file.

// Macro for drawing the high voltage of the GEM detectors from the Tango database as line graphs or as a surface
//
// The GEM high voltage (monitored values) is stored in Tango separately for every channel of the CAEN modules:
// gem/caen_240/mod_0/vmon_<N> and gem/caen_240/mod_2/vmon_<N>, N = 0, ..., 11
//
// Parameters:
//  period, start_run, end_run - period and run range (or start_date, end_date - time interval)
//  module - CAEN module ("mod_0" or "mod_2"), empty string - all the modules
//  draw_option - "lines" - line graphs, "lines3d" - line graphs in 3D view,
//                "SURF1" - surface, "LEGO2" - columns, "COLZ" - color map or another option of the 2D histogram
//
// Examples:
//  root -l 'tango_gem_hv.C(9, 9696, 9730)'                           - channels of mod_0 as line graphs
//  root -l 'tango_gem_hv.C(9, 9696, 9730, "", "SURF1")'              - channels of all the modules as a surface
//  root -l 'tango_gem_hv.C("2026-03-01 10:00:00", "2026-03-01 16:00:00", "mod_2", "LEGO2")'
 
// CAEN modules of the GEM high voltage stored in Tango
vector<TString> gem_hv_modules = {"mod_0", "mod_2"};
const int gem_hv_channel_count = 12;
 
// read the high voltage of the selected GEM channels and draw them,
// create_tango - function creating TangoData object for the given Tango attribute and member names
void DrawGemHV(function<TangoData*(const char*, const char*)> create_tango, TString module, TString draw_option)
{
    vector<TString> modules;
    if (module == "")
        modules = gem_hv_modules;
    else
        modules.push_back(module);
 
    // TangoData object for every channel, the channels without data are skipped
    vector<TangoData*> vecTangoData;
    cout << "Channel indices on the graph:" << endl;
    for (const auto& m : modules) {
        for (int ch = 0; ch < gem_hv_channel_count; ch++) {
            TString attribute_name = TString::Format("vmon_%d", ch);
            TangoData* tango_data = create_tango(attribute_name.Data(), m.Data());
            if ((tango_data->GetTangoData() == nullptr) || (tango_data->GetTangoData()->IsEmpty())) {
                cout << "WARNING: there are no Tango data for " << m << "/" << attribute_name
                     << ", the channel is skipped" << endl;
                delete tango_data;
                continue;
            }
            cout << "   " << vecTangoData.size() << " - " << m << ", " << attribute_name << endl;
            vecTangoData.push_back(tango_data);
        }
    }
    if (vecTangoData.empty()) {
        cout << "Macro finished with errors: there are no Tango data for the selected GEM channels" << endl;
        return;
    }
 
    TString option(draw_option);
    option.ToLower();
    int result;
    if ((option == "lines") || (option == "lines3d"))
        result = TangoData::PrintTangoArrayMultiGraph(vecTangoData, "voltage", option == "lines3d");
    else
        result = TangoData::PrintTangoArraySurface(vecTangoData, "voltage", "channel index", draw_option);
 
    for (auto tango_data : vecTangoData)
        delete tango_data;
 
    if (result < 0)
        cout << "Macro finished with errors" << endl;
    else
        cout << "Macro finished successfully" << endl;
}
 
// for the run range of the period
void tango_gem_hv(int period = 9,
                  int start_run = 9696,
                  int end_run = 9730,
                  TString module = "mod_0",
                  TString draw_option = "lines")
{
    DrawGemHV(
        [&](const char* attribute_name, const char* member_name) {
            return new TangoData(period, start_run, end_run, "gem", attribute_name, member_name);
        },
        module, draw_option);
}
 
// for the time interval
void tango_gem_hv(const char* start_date, const char* end_date, TString module = "mod_0", TString draw_option = "lines")
{
    DrawGemHV(
        [&](const char* attribute_name, const char* member_name) {
            return new TangoData(start_date, end_date, "gem", attribute_name, member_name);
        },
        module, draw_option);
}