CondCore/SiPixelPlugins/test/testSiPixelPayloadInspector.cpp - master - githubmirror/cmssw

CondCore/SiPixelPlugins/test/testSiPixelPayloadInspector.cpp
277 строк12 KB
mmusich
add class to display a pixel tracker map of the difference of two SiPixelQuality payloads
19 сен 2024, 15:01
19 сен 2024, 15:01afb887a
#include <iostream>
#include <sstream>
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/SiPixelPlugins/plugins/SiPixelLorentzAngle_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelQuality_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelGainCalibrationOffline_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelTemplateDBObject_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelGenErrorDBObject_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelVCal_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelQualityProbabilities_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelDynamicInefficiency_PayloadInspector.cc"
#include "CondCore/SiPixelPlugins/plugins/SiPixelFEDChannelContainer_PayloadInspector.cc"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#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) {
  Py_Initialize();

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

  // Lorentz Angle

  std::string tag = "SiPixelLorentzAngle_v11_offline";
  cond::Time_t start = static_cast<unsigned long long>(303790);
  cond::Time_t end = static_cast<unsigned long long>(324245);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising Lorentz Angle plots " << std::endl;

  SiPixelLorentzAngleValues histo1;
  histo1.process(connectionString, PI::mk_input(tag, start, start));
  edm::LogPrint("testSiPixelPayloadInspector") << histo1.data() << std::endl;

  SiPixelLorentzAngleValueComparisonSingleTag histo2;
  histo2.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo2.data() << std::endl;

  SiPixelLorentzAngleByRegionComparisonSingleTag histo3;
  histo3.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo3.data() << std::endl;

  SiPixelBPixLorentzAngleMap histo4;
  histo4.process(connectionString, PI::mk_input(tag, start, start));
  edm::LogPrint("testSiPixelPayloadInspector") << histo4.data() << std::endl;

  SiPixelFPixLorentzAngleMap histo5;
  histo5.process(connectionString, PI::mk_input(tag, end, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo5.data() << std::endl;

  std::string f_tagPhase2 = "SiPixelLorentzAngle_phase2_T15_v5_mc";
  std::string l_tagPhase2 = "SiPixelLorentzAngle_phase2_T19_v1_mc";

  SiPixelLorentzAngleValuesBarrelCompareTwoTags histoPhase2;
  histoPhase2.process(connectionString, PI::mk_input(f_tagPhase2, 1, 1, l_tagPhase2, 1, 1));
  edm::LogPrint("testSiPixelPayloadInspector") << histoPhase2.data() << std::endl;

  // 2 tags comparisons

  std::string tag2 = "SiPixelLorentzAngle_2016_ultralegacymc_v2";
  cond::Time_t start2 = static_cast<unsigned long long>(1);

  SiPixelLorentzAngleValueComparisonTwoTags histo6;
  histo6.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
  edm::LogPrint("testSiPixelPayloadInspector") << histo6.data() << std::endl;

  SiPixelLorentzAngleByRegionComparisonTwoTags histo7;
  histo7.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
  edm::LogPrint("testSiPixelPayloadInspector") << histo7.data() << std::endl;

  // SiPixelQuality

  tag = "SiPixelQuality_forDigitizer_phase1_2018_permanentlyBad";
  tag2 = "SiPixelQuality_phase1_2023_v6_mc";
  start = static_cast<unsigned long long>(1);
  end = static_cast<unsigned long long>(1);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelQuality plots " << std::endl;

  SiPixelBPixQualityMap histo8;
  histo8.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo8.data() << std::endl;

  SiPixelFPixQualityMap histo9;
  histo9.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo9.data() << std::endl;

  SiPixelQualityBadFractionMap histo9bis;
  histo9bis.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo9bis.data() << std::endl;

  SiPixelQualityBadFracCompareTwoTags histo9ter;
  histo9ter.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
  edm::LogPrint("testSiPixelPayloadInspector") << histo9ter.data() << std::endl;

  // SiPixelGainCalibrationOffline

  tag = "SiPixelGainCalibration_2009runs_express";
  start = static_cast<unsigned long long>(312203);
  end = static_cast<unsigned long long>(312203);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelGainCalibrationOffline plots " << std::endl;

  SiPixelGainCalibrationOfflineGainsValues histo10;
  histo10.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo10.data() << std::endl;

  SiPixelGainCalibrationOfflinePedestalsValues histo11;
  histo11.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo11.data() << std::endl;

  SiPixelGainCalibrationOfflineGainsByPart histo12;
  histo12.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo12.data() << std::endl;

  SiPixelGainCalibrationOfflinePedestalsByPart histo13;
  histo13.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo13.data() << std::endl;

  end = static_cast<unsigned long long>(326851);

  SiPixelGainCalibOfflinePedestalComparisonSingleTag histo14;
  histo14.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo14.data() << std::endl;

  SiPixelGainCalibOfflineGainByRegionComparisonSingleTag histo15;
  histo15.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo15.data() << std::endl;

  SiPixelGainCalibrationOfflineCorrelations histo16;
  histo16.process(connectionString, PI::mk_input(tag, end, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo16.data() << std::endl;

  py::dict inputs;
  inputs["SetLog"] = "True";  // sets to true, 1,True,Yes will work

  SiPixelGainCalibrationOfflineGainsValuesBarrel histo17;
  histo17.setInputParamValues(inputs);
  histo17.process(connectionString, PI::mk_input(tag, end, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo17.data() << std::endl;

  // SiPixelTemplates

  tag = "SiPixelTemplateDBObject38Tv3_express";
  start = static_cast<unsigned long long>(326083);
  end = static_cast<unsigned long long>(326083);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelTemplates plots " << std::endl;

  SiPixelTemplateIDsBPixMap histo18;
  histo18.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo18.data() << std::endl;

  SiPixelTemplateLAFPixMap histo19;
  histo19.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo19.data() << std::endl;

  SiPixelTemplateQScaleMap histoQscale;
  histoQscale.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histoQscale.data() << std::endl;

  // SiPixelVCal

  tag = "SiPixelVCal_v1";

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelVCal plots " << std::endl;

  SiPixelVCalValues histo20;
  histo20.process("frontier://FrontierPrep/CMS_CONDITIONS", PI::mk_input(tag, 1, 1));
  edm::LogPrint("testSiPixelPayloadInspector") << histo20.data() << std::endl;

  SiPixelVCalSlopeValuesBarrel histo21;
  histo21.process("frontier://FrontierPrep/CMS_CONDITIONS", PI::mk_input(tag, 1, 1));
  edm::LogPrint("testSiPixelPayloadInspector") << histo21.data() << std::endl;

  SiPixelVCalOffsetValuesEndcap histo22;
  histo22.process("frontier://FrontierPrep/CMS_CONDITIONS", PI::mk_input(tag, 1, 1));
  edm::LogPrint("testSiPixelPayloadInspector") << histo22.data() << std::endl;

  // SiPixelGenErrors

  tag = "SiPixelGenErrorDBObject_phase1_BoR3_HV350_Tr2000";
  start = static_cast<unsigned long long>(1);
  end = static_cast<unsigned long long>(1);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelGenErrors plots " << std::endl;

  SiPixelGenErrorHeaderTable histo23;
  histo23.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo23.data() << std::endl;

  SiPixelGenErrorIDsBPixMap histo24;
  histo24.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo24.data() << std::endl;

  // SiPixelQualityProbabilities

  tag = "SiPixelQualityProbabilities_UltraLegacy2018_v0_mc";
  start = static_cast<unsigned long long>(1);
  end = static_cast<unsigned long long>(1);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelQualityProbabilities plots " << std::endl;

  SiPixelQualityProbabilitiesScenariosCount histo25;
  histo25.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo25.data() << std::endl;

  // SiPixelDynamicInefficiency

  tag = "SiPixelDynamicInefficiency_PhaseI_v9";
  start = static_cast<unsigned long long>(1);
  end = static_cast<unsigned long long>(1);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixeDynamicInefficiency plots " << std::endl;

  SiPixelFullIneffROCfromDynIneffMap histo26;
  histo26.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo26.data() << std::endl;

  tag2 = "SiPixelDynamicInefficiency_PhaseI_v6";
  start2 = static_cast<unsigned long long>(1);

  SiPixelFullIneffROCsMapCompareTwoTags histo27;
  histo27.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
  edm::LogPrint("testSiPixelPayloadInspector") << histo27.data() << std::endl;

  SiPixelDynamicInefficiencyPUPixelMaps histo28;
  histo28.process(connectionString, PI::mk_input(tag2, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo28.data() << std::endl;

  SiPixelDynamicInefficiencyPUParametrization histo29;
  histo29.process(connectionString, PI::mk_input(tag2, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo29.data() << std::endl;

  SiPixelDynamicInefficiencyPUParamComparisonTwoTags histo30;
  histo30.process(connectionString, PI::mk_input(tag, start, end, tag2, start2, start2));
  edm::LogPrint("testSiPixelPayloadInspector") << histo30.data() << std::endl;

  // SiPixelFEDChannelContainer
  tag = "SiPixelStatusScenarios_StuckTBMandOther_2023_v2_mc";
  tag2 = "SiPixelQualityProbabilities_2023_v2_mc";
  start = static_cast<unsigned long long>(1);
  end = static_cast<unsigned long long>(1);

  edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelFEDChannelContainer plots " << std::endl;

  inputs["SiPixelQualityProbabilitiesTag"] = tag2;  // Quality Probabilities tag to use
  SiPixelBPixFEDChannelContainerWeightedMap histo31;
  histo31.setInputParamValues(inputs);
  histo31.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo31.data() << std::endl;

  inputs["Scenarios"] = "370097_302";
  SiPixelBPixFEDChannelContainerMap histo32;
  histo32.setInputParamValues(inputs);
  histo32.process(connectionString, PI::mk_input(tag, start, end));
  edm::LogPrint("testSiPixelPayloadInspector") << histo32.data() << std::endl;

  inputs.clear();
#if PY_MAJOR_VERSION >= 3
  // TODO I don't know why this Py_INCREF is necessary...
  Py_INCREF(inputs.ptr());
#endif

  Py_Finalize();
}