Skip to content

Commit

Permalink
Merge pull request #22720 from mmusich/runInfoPI
Browse files Browse the repository at this point in the history 
introduce RunInfo Payload Inspector
  • Loading branch information
@cmsbuild
cmsbuild committed Mar 30, 2018
2 parents f26ea7c + 6612356 commit a2ac776
Show file tree
Hide file tree
Showing 6 changed files with 421 additions and 0 deletions.
99 changes: 99 additions & 0 deletions CondCore/RunInfoPlugins/interface/RunInfoPayloadInspectoHelper.h
View file
@@ -0,0 +1,99 @@
#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
8 changes: 8 additions & 0 deletions CondCore/RunInfoPlugins/plugins/BuildFile.xml
View file
@@ -0,0 +1,8 @@
<flags CXX_FLAGS="-Wno-unused-variable"/>

<library file="RunInfo_PayloadInspector.cc" name="RunInfo_PayloadInspector">
<use name="CondCore/Utilities"/>
<use name="CondCore/CondDB"/>
<use name="CondFormats/Common"/>
<use name="boost_python"/>
</library>
233 changes: 233 additions & 0 deletions CondCore/RunInfoPlugins/plugins/RunInfo_PayloadInspector.cc
View file
@@ -0,0 +1,233 @@
/*!
\file RunInfo_PayloadInspector
\Payload Inspector Plugin for RunInfo
\author M. Musich
\version $Revision: 1.0 $
\date $Date: 2018/03/18 10:01:00 $
*/

#include "CondCore/Utilities/interface/PayloadInspectorModule.h"
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/CondDB/interface/Time.h"

// the data format of the condition to be inspected
#include "CondFormats/RunInfo/interface/RunInfo.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

// helper
#include "CondCore/RunInfoPlugins/interface/RunInfoPayloadInspectoHelper.h"

// system includes
#include <memory>
#include <sstream>
#include <iostream>

// include ROOT
#include "TProfile.h"
#include "TH2F.h"
#include "TLegend.h"
#include "TCanvas.h"
#include "TLine.h"
#include "TStyle.h"
#include "TLatex.h"
#include "TPave.h"
#include "TPaveStats.h"
#include "TPaletteAxis.h"

namespace {

/************************************************
RunInfo Payload Inspector of 1 IOV
*************************************************/
class RunInfoTest : public cond::payloadInspector::Histogram1D<RunInfo> {

public:
RunInfoTest() : cond::payloadInspector::Histogram1D<RunInfo>( "Test RunInfo", "Test RunInfo",10,0.0,10.0)
{
Base::setSingleIov( true );
}

bool fill( const std::vector<std::tuple<cond::Time_t,cond::Hash> >& iovs ) override{
auto iov = iovs.front();
std::shared_ptr<RunInfo> payload = Base::fetchPayload( std::get<1>(iov) );

if(payload.get() ) {
payload->printAllValues();
}
return true;
}
};

/************************************************
Summary of RunInfo of 1 IOV
*************************************************/
class RunInfoParameters : public cond::payloadInspector::PlotImage<RunInfo> {
public:
RunInfoParameters() : cond::payloadInspector::PlotImage<RunInfo>( "Display of RunInfo parameters" ){
setSingleIov( true );
}

bool fill( const std::vector<std::tuple<cond::Time_t,cond::Hash> >& iovs ) override{
auto iov = iovs.front();
std::shared_ptr<RunInfo> payload = fetchPayload( std::get<1>(iov) );

TCanvas canvas("Beam Spot Parameters Summary","RunInfo Parameters summary",1000,1000);
canvas.cd();

gStyle->SetHistMinimumZero();

canvas.SetTopMargin(0.08);
canvas.SetBottomMargin(0.06);
canvas.SetLeftMargin(0.3);
canvas.SetRightMargin(0.02);
canvas.Modified();
canvas.SetGrid();

auto h2_RunInfoParameters = std::unique_ptr<TH2F>(new TH2F("Parameters","",1,0.0,1.0,11,0,11.));
auto h2_RunInfoState = std::unique_ptr<TH2F>(new TH2F("State","",1,0.0,1.0,11,0,11.));
h2_RunInfoParameters->SetStats(false);
h2_RunInfoState->SetStats(false);

float fieldIntensity = RunInfoPI::theBField(payload->m_avg_current);

std::function<float(RunInfoPI::parameters)> cutFunctor = [&payload,fieldIntensity](RunInfoPI::parameters my_param) {
float ret(-999.);
switch(my_param){
case RunInfoPI::m_run : return float(payload->m_run);
case RunInfoPI::m_start_time_ll : return float(payload->m_start_time_ll);
case RunInfoPI::m_stop_time_ll : return float(payload->m_stop_time_ll);
case RunInfoPI::m_start_current : return payload->m_start_current;
case RunInfoPI::m_stop_current : return payload->m_stop_current;
case RunInfoPI::m_avg_current : return payload->m_avg_current;
case RunInfoPI::m_max_current : return payload->m_max_current;
case RunInfoPI::m_min_current : return payload->m_min_current;
case RunInfoPI::m_run_intervall_micros : return payload->m_run_intervall_micros;
case RunInfoPI::m_BField : return fieldIntensity;
case RunInfoPI::m_fedIN : return float((payload->m_fed_in).size());
case RunInfoPI::END_OF_TYPES : return ret;
default : return ret;
}
};

h2_RunInfoParameters->GetXaxis()->SetBinLabel(1,"Value");
h2_RunInfoState->GetXaxis()->SetBinLabel(1,"Value");

unsigned int yBin=11;
for(int foo = RunInfoPI::m_run; foo != RunInfoPI::END_OF_TYPES; foo++ ){
RunInfoPI::parameters param = static_cast<RunInfoPI::parameters>(foo);
std::string theLabel = RunInfoPI::getStringFromTypeEnum(param);
h2_RunInfoState->GetYaxis()->SetBinLabel(yBin,theLabel.c_str());
h2_RunInfoParameters->GetYaxis()->SetBinLabel(yBin,theLabel.c_str());
h2_RunInfoParameters->SetBinContent(1,yBin,cutFunctor(param));
// non-fake payload
if((payload->m_run)!=-1){
if ((payload->m_avg_current)<=-1){
// go in error state
h2_RunInfoState->SetBinContent(1,yBin,0.);
} else {
// all is OK
h2_RunInfoState->SetBinContent(1,yBin,1.);
}
} else {
// this is a fake payload
h2_RunInfoState->SetBinContent(1,yBin,0.9);
}
yBin--;

}

h2_RunInfoParameters->GetXaxis()->LabelsOption("h");
h2_RunInfoParameters->GetYaxis()->SetLabelSize(0.05);
h2_RunInfoParameters->GetXaxis()->SetLabelSize(0.05);
h2_RunInfoParameters->SetMarkerSize(1.5);

h2_RunInfoState->GetXaxis()->LabelsOption("h");
h2_RunInfoState->GetYaxis()->SetLabelSize(0.05);
h2_RunInfoState->GetXaxis()->SetLabelSize(0.05);
h2_RunInfoState->SetMarkerSize(1.5);

RunInfoPI::reportSummaryMapPalette(h2_RunInfoState.get());
h2_RunInfoState->Draw("col");

h2_RunInfoParameters->Draw("TEXTsame");

TLatex t1;
t1.SetNDC();
t1.SetTextAlign(12);
t1.SetTextSize(0.03);
t1.DrawLatex(0.1, 0.98,"RunInfo parameters:");
t1.DrawLatex(0.1, 0.95,"payload:");

t1.SetTextFont(42);
t1.SetTextColor(4);
t1.DrawLatex(0.37, 0.982,Form("IOV %s",std::to_string(+std::get<0>(iov)).c_str()));
t1.DrawLatex(0.21, 0.952,Form(" %s",(std::get<1>(iov)).c_str()));

std::string fileName(m_imageFileName);
canvas.SaveAs(fileName.c_str());

return true;
}


};

/************************************************
time history of Magnet currents from RunInfo
*************************************************/

template<RunInfoPI::parameters param> class RunInfoCurrentHistory : public cond::payloadInspector::HistoryPlot<RunInfo,float> {
public:
RunInfoCurrentHistory() : cond::payloadInspector::HistoryPlot<RunInfo,float>(getStringFromTypeEnum(param),getStringFromTypeEnum(param)+" value"){}
~RunInfoCurrentHistory() override = default;

float getFromPayload( RunInfo& payload ) override{

float fieldIntensity = RunInfoPI::theBField(payload.m_avg_current);

switch(param){
case RunInfoPI::m_start_current : return payload.m_start_current;
case RunInfoPI::m_stop_current : return payload.m_stop_current;
case RunInfoPI::m_avg_current : return payload.m_avg_current;
case RunInfoPI::m_max_current : return payload.m_max_current;
case RunInfoPI::m_min_current : return payload.m_min_current;
case RunInfoPI::m_BField : return fieldIntensity;
default:
edm::LogWarning("LogicError") << "Unknown parameter: " << param;
break;
}

} // payload

/************************************************/
std::string getStringFromTypeEnum (const RunInfoPI::parameters &parameter){
switch(parameter){
case RunInfoPI::m_start_current : return "Magent start current [A]";
case RunInfoPI::m_stop_current : return "Magnet stop current [A]";
case RunInfoPI::m_avg_current : return "Magnet average current [A]";
case RunInfoPI::m_max_current : return "Magnet max current [A]";
case RunInfoPI::m_min_current : return "Magnet min current [A]";
case RunInfoPI::m_BField : return "B-field intensity [T]";
default: return "should never be here";
}
}
};

typedef RunInfoCurrentHistory<RunInfoPI::m_start_current> RunInfoStartCurrentHistory;
typedef RunInfoCurrentHistory<RunInfoPI::m_stop_current> RunInfoStopCurrentHistory;
typedef RunInfoCurrentHistory<RunInfoPI::m_avg_current> RunInfoAverageCurrentHistory;
typedef RunInfoCurrentHistory<RunInfoPI::m_max_current> RunInfoMaxCurrentHistory;
typedef RunInfoCurrentHistory<RunInfoPI::m_min_current> RunInfoMinCurrentHistory;
typedef RunInfoCurrentHistory<RunInfoPI::m_BField> RunInfoBFieldHistory;

} // close namespace

PAYLOAD_INSPECTOR_MODULE( RunInfo ){
PAYLOAD_INSPECTOR_CLASS( RunInfoTest );
PAYLOAD_INSPECTOR_CLASS( RunInfoParameters ) ;
PAYLOAD_INSPECTOR_CLASS( RunInfoStopCurrentHistory ) ;
PAYLOAD_INSPECTOR_CLASS( RunInfoAverageCurrentHistory) ;
PAYLOAD_INSPECTOR_CLASS( RunInfoMaxCurrentHistory ) ;
PAYLOAD_INSPECTOR_CLASS( RunInfoMinCurrentHistory ) ;
PAYLOAD_INSPECTOR_CLASS( RunInfoBFieldHistory ) ;
}
8 changes: 8 additions & 0 deletions CondCore/RunInfoPlugins/test/BuildFile.xml
View file
@@ -0,0 +1,8 @@
<use name="CondCore/CondDB"/>
<use name="CondCore/Utilities"/>
<use name="CondFormats/Common"/>
<use name="FWCore/PluginManager"/>
<use name="FWCore/Utilities"/>
<bin
file="testRunInfoPayloadInspector.cpp" name="testRunInfoPayloadInspector">
</bin>
44 changes: 44 additions & 0 deletions CondCore/RunInfoPlugins/test/testRunInfoPayloadInspector.cpp
View file
@@ -0,0 +1,44 @@
#include<iostream>
#include<sstream>
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/RunInfoPlugins/plugins/RunInfo_PayloadInspector.cc"

#include "FWCore/PluginManager/interface/PluginManager.h"
#include "FWCore/PluginManager/interface/standard.h"
#include "FWCore/PluginManager/interface/SharedLibrary.h"
#include "FWCore/ServiceRegistry/interface/ServiceRegistry.h"

int main(int argc, char** argv) {

edmplugin::PluginManager::Config config;
edmplugin::PluginManager::configure(edmplugin::standard::config());

std::vector<edm::ParameterSet> psets;
edm::ParameterSet pSet;
pSet.addParameter("@service_type",std::string("SiteLocalConfigService"));
psets.push_back(pSet);
edm::ServiceToken servToken(edm::ServiceRegistry::createSet(psets));
edm::ServiceRegistry::Operate operate(servToken);

std::string connectionString("frontier://FrontierProd/CMS_CONDITIONS");

std::string tag = "runinfo_31X_hlt";
std::string runTimeType = cond::time::timeTypeName(cond::runnumber);
cond::Time_t start = boost::lexical_cast<unsigned long long>(311950);
cond::Time_t end = boost::lexical_cast<unsigned long long>(312237);

std::cout <<"## RunInfo testing"<<std::endl;

RunInfoTest histo0;
histo0.process( connectionString, tag, runTimeType,end,end);
std::cout <<histo0.data()<<std::endl;

RunInfoParameters histo1;
histo1.process( connectionString, tag, runTimeType,end,end);
std::cout <<histo1.data()<<std::endl;

RunInfoBFieldHistory histo2;
histo2.process( connectionString, tag, runTimeType, start, end );
std::cout <<histo2.data()<<std::endl;

}

0 comments on commit a2ac776

Please sign in to comment.