RunInfoPayloadInspectoHelper.h

Go to the documentation of this file.

#ifndef CONDCORE_RUNINFOPLUGINS_RUNINFOPAYLOADINSPECTORHELPER_H
#define CONDCORE_RUNINFOPLUGINS_RUNINFOPAYLOADINSPECTORHELPER_H

#include <vector>
#include <string>
#include "TH1.h"
#include "TH2.h"
#include "TStyle.h"
#include "TPaveText.h"
#include "TStyle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

namespace RunInfoPI {

  // values are taken from https://github.com/cms-sw/cmssw/blob/master/MagneticField/GeomBuilder/plugins/VolumeBasedMagneticFieldESProducerFromDB.cc#L74-L75
  constexpr std::array<int, 7> nominalCurrents{{-1, 0, 9558, 14416, 16819, 18268, 19262}};
  constexpr std::array<float, 7> nominalFields{{3.8, 0., 2., 3., 3.5, 3.8, 4.}};

  // all parameter than can be displayed
  enum parameters {
    m_run,                   // int
    m_start_time_ll,         // long long;
    m_stop_time_ll,          // long long
    m_start_current,         // float
    m_stop_current,          // float
    m_avg_current,           // float
    m_max_current,           // float
    m_min_current,           // float
    m_run_intervall_micros,  // float
    m_fedIN,                 // unsigned int
    m_BField,                // float
    END_OF_TYPES
  };

  /************************************************/
  float theBField(const float current) {
    // logic is taken from https://github.com/cms-sw/cmssw/blob/master/MagneticField/GeomBuilder/plugins/VolumeBasedMagneticFieldESProducerFromDB.cc#L156

    int i = 0;
    for (; i < (int)nominalFields.size() - 1; i++) {
      if (2 * current < nominalCurrents[i] + nominalCurrents[i + 1]) {
        return nominalFields[i];
      }
    }
    return nominalFields[i];
  }

  /************************************************/
  std::string getStringFromTypeEnum(const parameters& parameter) {
    switch (parameter) {
      case m_run:
        return "run number";
      case m_start_time_ll:
        return "start time";
      case m_stop_time_ll:
        return "stop time";
      case m_start_current:
        return "start current [A]";
      case m_stop_current:
        return "stop current [A]";
      case m_avg_current:
        return "average current [A]";
      case m_max_current:
        return "max current [A]";
      case m_min_current:
        return "min current [A]";
      case m_run_intervall_micros:
        return "run duration [#mus]";
      case m_fedIN:
        return "n. FEDs";
      case m_BField:
        return "B-field intensity [T]";
      default:
        return "should never be here";
    }
  }

  /************************************************/
  void reportSummaryMapPalette(TH2* obj) {
    static int pcol[20];

    float rgb[20][3];

    for (int i = 0; i < 20; i++) {
      if (i < 17) {
        rgb[i][0] = 0.80 + 0.01 * i;
        rgb[i][1] = 0.00 + 0.03 * i;
        rgb[i][2] = 0.00;
      } else if (i < 19) {
        rgb[i][0] = 0.80 + 0.01 * i;
        rgb[i][1] = 0.00 + 0.03 * i + 0.15 + 0.10 * (i - 17);
        rgb[i][2] = 0.00;
      } else if (i == 19) {
        rgb[i][0] = 0.00;
        rgb[i][1] = 0.80;
        rgb[i][2] = 0.00;
      }
      pcol[i] = TColor::GetColor(rgb[i][0], rgb[i][1], rgb[i][2]);
    }

    gStyle->SetPalette(20, pcol);

    if (obj) {
      obj->SetMinimum(-1.e-15);
      obj->SetMaximum(+1.0);
      obj->SetOption("colz");
    }
  }

};  // namespace RunInfoPI
#endif