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");
    }
  }

};
#endif