tango_fhcal_hv.C

Go to the documentation of this file.

// Macro for drawing the high voltage of the FHCal from the Tango database as line graphs or as a surface
//
// The FHCal high voltage is stored in Tango separately for every section of every module:
// fhcal/hv/hv_1/vch_<module>_<section>, module = 2, 4, ..., 116 (except 64, 78, 92, 114), section = 1, ..., 10
//
// Parameters:
//  period, start_run, end_run - period and run range (or start_date, end_date - time interval)
//  module - FHCal module number to show all its sections, 0 - show the given section of all the modules
//  section - section number to show for all the modules (used only if module = 0)
//  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_fhcal_hv.C(9, 9696, 9730)'                     - all sections of module 2 as line graphs
//  root -l 'tango_fhcal_hv.C(9, 9696, 9730, 0, 1, "SURF1")'      - section 1 of all the modules as a surface
//  root -l 'tango_fhcal_hv.C("2026-03-01 10:00:00", "2026-03-01 16:00:00", 2, 1, "LEGO2")'
 
// FHCal modules whose high voltage is stored in Tango
vector<int> GetFHCalModules()
{
    vector<int> modules;
    for (int module = 2; module <= 116; module += 2)
        if ((module != 64) && (module != 78) && (module != 92) && (module != 114))
            modules.push_back(module);
    return modules;
}
 
// read the high voltage of the selected FHCal channels and draw them,
// create_tango - function creating TangoData object for the given Tango attribute name
void DrawFHCalHV(function<TangoData*(const char*)> create_tango, int module, int section, TString draw_option)
{
    vector<int> modules, sections;
    if (module > 0) {
        modules.push_back(module);
        for (int s = 1; s <= 10; s++)
            sections.push_back(s);
    } else {
        if ((section < 1) || (section > 10)) {
            cout << "ERROR: section number should be from 1 to 10" << endl;
            return;
        }
        modules = GetFHCalModules();
        sections.push_back(section);
    }
 
    // TangoData object for every channel, the channels without data are skipped
    vector<TangoData*> vecTangoData;
    cout << "Channel indices on the graph:" << endl;
    for (int m : modules) {
        for (int s : sections) {
            TString attribute_name = TString::Format("vch_%d_%d", m, s);
            TangoData* tango_data = create_tango(attribute_name.Data());
            if ((tango_data->GetTangoData() == nullptr) || (tango_data->GetTangoData()->IsEmpty())) {
                cout << "WARNING: there are no Tango data for " << attribute_name << ", the channel is skipped" << endl;
                delete tango_data;
                continue;
            }
            cout << "   " << vecTangoData.size() << " - module " << m << ", section " << s << endl;
            vecTangoData.push_back(tango_data);
        }
    }
    if (vecTangoData.empty()) {
        cout << "Macro finished with errors: there are no Tango data for the selected FHCal 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",
                                                   module > 0 ? "section index" : "module 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_fhcal_hv(int period = 9,
                    int start_run = 9696,
                    int end_run = 9730,
                    int module = 2,
                    int section = 1,
                    TString draw_option = "lines")
{
    DrawFHCalHV(
        [&](const char* attribute_name) {
            return new TangoData(period, start_run, end_run, "fhcal", attribute_name, "hv_1");
        },
        module, section, draw_option);
}
 
// for the time interval
void tango_fhcal_hv(const char* start_date,
                    const char* end_date,
                    int module = 2,
                    int section = 1,
                    TString draw_option = "lines")
{
    DrawFHCalHV(
        [&](const char* attribute_name) {
            return new TangoData(start_date, end_date, "fhcal", attribute_name, "hv_1");
        },
        module, section, draw_option);
}