SiPixelPlugins/test/testSiPixelPayloadInspector.cpp

0001 #include <iostream>
0002 #include <sstream>
0003 #include "CondCore/Utilities/interface/PayloadInspector.h"
0004 #include "CondCore/SiPixelPlugins/plugins/SiPixelLorentzAngle_PayloadInspector.cc"
0005 #include "CondCore/SiPixelPlugins/plugins/SiPixelQuality_PayloadInspector.cc"
0006 #include "CondCore/SiPixelPlugins/plugins/SiPixelGainCalibrationOffline_PayloadInspector.cc"
0007 #include "CondCore/SiPixelPlugins/plugins/SiPixelTemplateDBObject_PayloadInspector.cc"
0008 #include "CondCore/SiPixelPlugins/plugins/SiPixelGenErrorDBObject_PayloadInspector.cc"
0009 #include "CondCore/SiPixelPlugins/plugins/SiPixelVCal_PayloadInspector.cc"
0010 #include "CondCore/SiPixelPlugins/plugins/SiPixelQualityProbabilities_PayloadInspector.cc"
0011 #include "CondCore/SiPixelPlugins/plugins/SiPixelDynamicInefficiency_PayloadInspector.cc"
0012 #include "CondCore/SiPixelPlugins/plugins/SiPixelFEDChannelContainer_PayloadInspector.cc"
0013 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0014 #include "FWCore/PluginManager/interface/PluginManager.h"
0015 #include "FWCore/PluginManager/interface/standard.h"
0016 #include "FWCore/PluginManager/interface/SharedLibrary.h"
0017 #include "FWCore/ServiceRegistry/interface/ServiceRegistry.h"
0018
0019 int main(int argc, char** argv) {
0020   Py_Initialize();
0021
0022   edmplugin::PluginManager::Config config;
0023   edmplugin::PluginManager::configure(edmplugin::standard::config());
0024
0025   std::vector<edm::ParameterSet> psets;
0026   edm::ParameterSet pSet;
0027   pSet.addParameter("@service_type", std::string("SiteLocalConfigService"));
0028   psets.push_back(pSet);
0029   edm::ServiceToken servToken(edm::ServiceRegistry::createSet(psets));
0030   edm::ServiceRegistry::Operate operate(servToken);
0031
0032   std::string connectionString("frontier://FrontierProd/CMS_CONDITIONS");
0033
0034   // Lorentz Angle
0035
0036   std::string tag = "SiPixelLorentzAngle_v11_offline";
0037   cond::Time_t start = static_cast<unsigned long long>(303790);
0038   cond::Time_t end = static_cast<unsigned long long>(324245);
0039
0040   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising Lorentz Angle plots " << std::endl;
0041
0042   SiPixelLorentzAngleValues histo1;
0043   histo1.process(connectionString, PI::mk_input(tag, start, start));
0044   edm::LogPrint("testSiPixelPayloadInspector") << histo1.data() << std::endl;
0045
0046   SiPixelLorentzAngleValueComparisonSingleTag histo2;
0047   histo2.process(connectionString, PI::mk_input(tag, start, end));
0048   edm::LogPrint("testSiPixelPayloadInspector") << histo2.data() << std::endl;
0049
0050   SiPixelLorentzAngleByRegionComparisonSingleTag histo3;
0051   histo3.process(connectionString, PI::mk_input(tag, start, end));
0052   edm::LogPrint("testSiPixelPayloadInspector") << histo3.data() << std::endl;
0053
0054   SiPixelBPixLorentzAngleMap histo4;
0055   histo4.process(connectionString, PI::mk_input(tag, start, start));
0056   edm::LogPrint("testSiPixelPayloadInspector") << histo4.data() << std::endl;
0057
0058   SiPixelFPixLorentzAngleMap histo5;
0059   histo5.process(connectionString, PI::mk_input(tag, end, end));
0060   edm::LogPrint("testSiPixelPayloadInspector") << histo5.data() << std::endl;
0061
0062   std::string f_tagPhase2 = "SiPixelLorentzAngle_phase2_T15_v5_mc";
0063   std::string l_tagPhase2 = "SiPixelLorentzAngle_phase2_T19_v1_mc";
0064
0065   SiPixelLorentzAngleValuesBarrelCompareTwoTags histoPhase2;
0066   histoPhase2.process(connectionString, PI::mk_input(f_tagPhase2, 1, 1, l_tagPhase2, 1, 1));
0067   edm::LogPrint("testSiPixelPayloadInspector") << histoPhase2.data() << std::endl;
0068
0069   // 2 tags comparisons
0070
0071   std::string tag2 = "SiPixelLorentzAngle_2016_ultralegacymc_v2";
0072   cond::Time_t start2 = static_cast<unsigned long long>(1);
0073
0074   SiPixelLorentzAngleValueComparisonTwoTags histo6;
0075   histo6.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
0076   edm::LogPrint("testSiPixelPayloadInspector") << histo6.data() << std::endl;
0077
0078   SiPixelLorentzAngleByRegionComparisonTwoTags histo7;
0079   histo7.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
0080   edm::LogPrint("testSiPixelPayloadInspector") << histo7.data() << std::endl;
0081
0082   // SiPixelQuality
0083
0084   tag = "SiPixelQuality_forDigitizer_phase1_2018_permanentlyBad";
0085   start = static_cast<unsigned long long>(1);
0086   end = static_cast<unsigned long long>(1);
0087
0088   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelQuality plots " << std::endl;
0089
0090   SiPixelBPixQualityMap histo8;
0091   histo8.process(connectionString, PI::mk_input(tag, start, end));
0092   edm::LogPrint("testSiPixelPayloadInspector") << histo8.data() << std::endl;
0093
0094   SiPixelFPixQualityMap histo9;
0095   histo9.process(connectionString, PI::mk_input(tag, start, end));
0096   edm::LogPrint("testSiPixelPayloadInspector") << histo9.data() << std::endl;
0097
0098   SiPixelQualityBadFractionMap histo9bis;
0099   histo9bis.process(connectionString, PI::mk_input(tag, start, end));
0100   edm::LogPrint("testSiPixelPayloadInspector") << histo9bis.data() << std::endl;
0101
0102   // SiPixelGainCalibrationOffline
0103
0104   tag = "SiPixelGainCalibration_2009runs_express";
0105   start = static_cast<unsigned long long>(312203);
0106   end = static_cast<unsigned long long>(312203);
0107
0108   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelGainCalibrationOffline plots " << std::endl;
0109
0110   SiPixelGainCalibrationOfflineGainsValues histo10;
0111   histo10.process(connectionString, PI::mk_input(tag, start, end));
0112   edm::LogPrint("testSiPixelPayloadInspector") << histo10.data() << std::endl;
0113
0114   SiPixelGainCalibrationOfflinePedestalsValues histo11;
0115   histo11.process(connectionString, PI::mk_input(tag, start, end));
0116   edm::LogPrint("testSiPixelPayloadInspector") << histo11.data() << std::endl;
0117
0118   SiPixelGainCalibrationOfflineGainsByPart histo12;
0119   histo12.process(connectionString, PI::mk_input(tag, start, end));
0120   edm::LogPrint("testSiPixelPayloadInspector") << histo12.data() << std::endl;
0121
0122   SiPixelGainCalibrationOfflinePedestalsByPart histo13;
0123   histo13.process(connectionString, PI::mk_input(tag, start, end));
0124   edm::LogPrint("testSiPixelPayloadInspector") << histo13.data() << std::endl;
0125
0126   end = static_cast<unsigned long long>(326851);
0127
0128   SiPixelGainCalibOfflinePedestalComparisonSingleTag histo14;
0129   histo14.process(connectionString, PI::mk_input(tag, start, end));
0130   edm::LogPrint("testSiPixelPayloadInspector") << histo14.data() << std::endl;
0131
0132   SiPixelGainCalibOfflineGainByRegionComparisonSingleTag histo15;
0133   histo15.process(connectionString, PI::mk_input(tag, start, end));
0134   edm::LogPrint("testSiPixelPayloadInspector") << histo15.data() << std::endl;
0135
0136   SiPixelGainCalibrationOfflineCorrelations histo16;
0137   histo16.process(connectionString, PI::mk_input(tag, end, end));
0138   edm::LogPrint("testSiPixelPayloadInspector") << histo16.data() << std::endl;
0139
0140   py::dict inputs;
0141   inputs["SetLog"] = "True";  // sets to true, 1,True,Yes will work
0142
0143   SiPixelGainCalibrationOfflineGainsValuesBarrel histo17;
0144   histo17.setInputParamValues(inputs);
0145   histo17.process(connectionString, PI::mk_input(tag, end, end));
0146   edm::LogPrint("testSiPixelPayloadInspector") << histo17.data() << std::endl;
0147
0148   // SiPixelTemplates
0149
0150   tag = "SiPixelTemplateDBObject38Tv3_express";
0151   start = static_cast<unsigned long long>(326083);
0152   end = static_cast<unsigned long long>(326083);
0153
0154   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelTemplates plots " << std::endl;
0155
0156   SiPixelTemplateIDsBPixMap histo18;
0157   histo18.process(connectionString, PI::mk_input(tag, start, end));
0158   edm::LogPrint("testSiPixelPayloadInspector") << histo18.data() << std::endl;
0159
0160   SiPixelTemplateLAFPixMap histo19;
0161   histo19.process(connectionString, PI::mk_input(tag, start, end));
0162   edm::LogPrint("testSiPixelPayloadInspector") << histo19.data() << std::endl;
0163
0164   SiPixelTemplateQScaleMap histoQscale;
0165   histoQscale.process(connectionString, PI::mk_input(tag, start, end));
0166   edm::LogPrint("testSiPixelPayloadInspector") << histoQscale.data() << std::endl;
0167
0168   // SiPixelVCal
0169
0170   tag = "SiPixelVCal_v1";
0171
0172   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelVCal plots " << std::endl;
0173
0174   SiPixelVCalValues histo20;
0175   histo20.process("frontier://FrontierPrep/CMS_CONDITIONS", PI::mk_input(tag, 1, 1));
0176   edm::LogPrint("testSiPixelPayloadInspector") << histo20.data() << std::endl;
0177
0178   SiPixelVCalSlopeValuesBarrel histo21;
0179   histo21.process("frontier://FrontierPrep/CMS_CONDITIONS", PI::mk_input(tag, 1, 1));
0180   edm::LogPrint("testSiPixelPayloadInspector") << histo21.data() << std::endl;
0181
0182   SiPixelVCalOffsetValuesEndcap histo22;
0183   histo22.process("frontier://FrontierPrep/CMS_CONDITIONS", PI::mk_input(tag, 1, 1));
0184   edm::LogPrint("testSiPixelPayloadInspector") << histo22.data() << std::endl;
0185
0186   // SiPixelGenErrors
0187
0188   tag = "SiPixelGenErrorDBObject_phase1_BoR3_HV350_Tr2000";
0189   start = static_cast<unsigned long long>(1);
0190   end = static_cast<unsigned long long>(1);
0191
0192   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelGenErrors plots " << std::endl;
0193
0194   SiPixelGenErrorHeaderTable histo23;
0195   histo23.process(connectionString, PI::mk_input(tag, start, end));
0196   edm::LogPrint("testSiPixelPayloadInspector") << histo23.data() << std::endl;
0197
0198   SiPixelGenErrorIDsBPixMap histo24;
0199   histo24.process(connectionString, PI::mk_input(tag, start, end));
0200   edm::LogPrint("testSiPixelPayloadInspector") << histo24.data() << std::endl;
0201
0202   // SiPixelQualityProbabilities
0203
0204   tag = "SiPixelQualityProbabilities_UltraLegacy2018_v0_mc";
0205   start = static_cast<unsigned long long>(1);
0206   end = static_cast<unsigned long long>(1);
0207
0208   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelQualityProbabilities plots " << std::endl;
0209
0210   SiPixelQualityProbabilitiesScenariosCount histo25;
0211   histo25.process(connectionString, PI::mk_input(tag, start, end));
0212   edm::LogPrint("testSiPixelPayloadInspector") << histo25.data() << std::endl;
0213
0214   // SiPixelDynamicInefficiency
0215
0216   tag = "SiPixelDynamicInefficiency_PhaseI_v9";
0217   start = static_cast<unsigned long long>(1);
0218   end = static_cast<unsigned long long>(1);
0219
0220   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixeDynamicInefficiency plots " << std::endl;
0221
0222   SiPixelFullIneffROCfromDynIneffMap histo26;
0223   histo26.process(connectionString, PI::mk_input(tag, start, end));
0224   edm::LogPrint("testSiPixelPayloadInspector") << histo26.data() << std::endl;
0225
0226   tag2 = "SiPixelDynamicInefficiency_PhaseI_v6";
0227   start2 = static_cast<unsigned long long>(1);
0228
0229   SiPixelFullIneffROCsMapCompareTwoTags histo27;
0230   histo27.process(connectionString, PI::mk_input(tag, start, start, tag2, start2, start2));
0231   edm::LogPrint("testSiPixelPayloadInspector") << histo27.data() << std::endl;
0232
0233   SiPixelDynamicInefficiencyPUPixelMaps histo28;
0234   histo28.process(connectionString, PI::mk_input(tag2, start, end));
0235   edm::LogPrint("testSiPixelPayloadInspector") << histo28.data() << std::endl;
0236
0237   SiPixelDynamicInefficiencyPUParametrization histo29;
0238   histo29.process(connectionString, PI::mk_input(tag2, start, end));
0239   edm::LogPrint("testSiPixelPayloadInspector") << histo29.data() << std::endl;
0240
0241   SiPixelDynamicInefficiencyPUParamComparisonTwoTags histo30;
0242   histo30.process(connectionString, PI::mk_input(tag, start, end, tag2, start2, start2));
0243   edm::LogPrint("testSiPixelPayloadInspector") << histo30.data() << std::endl;
0244
0245   // SiPixelFEDChannelContainer
0246   tag = "SiPixelStatusScenarios_StuckTBMandOther_2023_v2_mc";
0247   tag2 = "SiPixelQualityProbabilities_2023_v2_mc";
0248   start = static_cast<unsigned long long>(1);
0249   end = static_cast<unsigned long long>(1);
0250
0251   edm::LogPrint("testSiPixelPayloadInspector") << "## Exercising SiPixelFEDChannelContainer plots " << std::endl;
0252
0253   inputs["SiPixelQualityProbabilitiesTag"] = tag2;  // Quality Probabilities tag to use
0254   SiPixelBPixFEDChannelContainerWeightedMap histo31;
0255   histo31.setInputParamValues(inputs);
0256   histo31.process(connectionString, PI::mk_input(tag, start, end));
0257   edm::LogPrint("testSiPixelPayloadInspector") << histo31.data() << std::endl;
0258
0259   inputs["Scenarios"] = "370097_302";
0260   SiPixelBPixFEDChannelContainerMap histo32;
0261   histo32.setInputParamValues(inputs);
0262   histo32.process(connectionString, PI::mk_input(tag, start, end));
0263   edm::LogPrint("testSiPixelPayloadInspector") << histo32.data() << std::endl;
0264
0265   inputs.clear();
0266 #if PY_MAJOR_VERSION >= 3
0267   // TODO I don't know why this Py_INCREF is necessary...
0268   Py_INCREF(inputs.ptr());
0269 #endif
0270
0271   Py_Finalize();
0272 }