SiStripPlugins/test/testSiStripPayloadInspector.cpp

0001 #include <iostream>
0002 #include <sstream>
0003 #include "CondCore/Utilities/interface/PayloadInspector.h"
0004 #include "CondCore/SiStripPlugins/plugins/SiStripApvGain_PayloadInspector.cc"
0005 #include "CondCore/SiStripPlugins/plugins/SiStripNoises_PayloadInspector.cc"
0006 #include "CondCore/SiStripPlugins/plugins/SiStripPedestals_PayloadInspector.cc"
0007 #include "CondCore/SiStripPlugins/plugins/SiStripThreshold_PayloadInspector.cc"
0008 #include "CondCore/SiStripPlugins/plugins/SiStripLatency_PayloadInspector.cc"
0009 #include "CondCore/SiStripPlugins/plugins/SiStripFedCabling_PayloadInspector.cc"
0010 #include "CondCore/SiStripPlugins/plugins/SiStripBadStrip_PayloadInspector.cc"
0011
0012 #include "FWCore/PluginManager/interface/PluginManager.h"
0013 #include "FWCore/PluginManager/interface/standard.h"
0014 #include "FWCore/ServiceRegistry/interface/ServiceRegistry.h"
0015
0016 int main(int argc, char** argv) {
0017   Py_Initialize();
0018
0019   edmplugin::PluginManager::Config config;
0020   edmplugin::PluginManager::configure(edmplugin::standard::config());
0021
0022   std::vector<edm::ParameterSet> psets;
0023   edm::ParameterSet pSet;
0024   pSet.addParameter("@service_type", std::string("SiteLocalConfigService"));
0025   psets.push_back(pSet);
0026   edm::ServiceToken servToken(edm::ServiceRegistry::createSet(psets));
0027   edm::ServiceRegistry::Operate operate(servToken);
0028
0029   std::string connectionString("frontier://FrontierProd/CMS_CONDITIONS");
0030
0031   // Gains
0032
0033   std::string tag = "SiStripApvGain_FromParticles_GR10_v11_offline";
0034
0035   cond::Time_t start = static_cast<unsigned long long>(132440);
0036   cond::Time_t end = static_cast<unsigned long long>(285368);
0037   py::dict inputs;
0038
0039   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising Gains plots " << std::endl;
0040
0041   SiStripApvGainsAverageTrackerMap histo1;
0042   histo1.process(connectionString, PI::mk_input(tag, start, start));
0043   edm::LogPrint("testSiStripPayloadInspector") << histo1.data() << std::endl;
0044
0045   SiStripApvGainsAvgDeviationRatioWithPreviousIOVTrackerMap histo2;
0046   inputs["nsigma"] = "1";
0047   histo2.setInputParamValues(inputs);
0048   histo2.process(connectionString, PI::mk_input(tag, start, end));
0049   edm::LogPrint("testSiStripPayloadInspector") << histo2.data() << std::endl;
0050
0051   SiStripApvGainsMaxDeviationRatioWithPreviousIOVTrackerMap histo3;
0052   inputs["nsigma"] = "1";
0053   histo3.setInputParamValues(inputs);
0054   histo3.process(connectionString, PI::mk_input(tag, start, end));
0055   edm::LogPrint("testSiStripPayloadInspector") << histo3.data() << std::endl;
0056
0057   SiStripApvGainsValuesComparatorSingleTag histo4;
0058   histo4.process(connectionString, PI::mk_input(tag, start, end));
0059   edm::LogPrint("testSiStripPayloadInspector") << histo4.data() << std::endl;
0060
0061   SiStripApvGainsComparatorSingleTag histo5;
0062   histo5.process(connectionString, PI::mk_input(tag, start, end));
0063   edm::LogPrint("testSiStripPayloadInspector") << histo5.data() << std::endl;
0064
0065   SiStripApvGainsComparatorByRegionSingleTag histo6;
0066   histo6.process(connectionString, PI::mk_input(tag, start, end));
0067   edm::LogPrint("testSiStripPayloadInspector") << histo6.data() << std::endl;
0068
0069   SiStripApvGainsRatioComparatorByRegionSingleTag histo7;
0070   histo7.process(connectionString, PI::mk_input(tag, start, end));
0071   edm::LogPrint("testSiStripPayloadInspector") << histo7.data() << std::endl;
0072
0073   SiStripApvGainCompareByPartitionSingleTag histoGainByPart;
0074   histoGainByPart.process(connectionString, PI::mk_input(tag, start, end));
0075   edm::LogPrint("testSiStripPayloadInspector") << histoGainByPart.data() << std::endl;
0076
0077   SiStripApvGainRatioByPartitionSingleTag histoGainRatioByPart;
0078   histoGainRatioByPart.process(connectionString, PI::mk_input(tag, start, end));
0079   edm::LogPrint("testSiStripPayloadInspector") << histoGainRatioByPart.data() << std::endl;
0080
0081   // Noise
0082
0083   tag = "SiStripNoise_GR10_v1_hlt";
0084   start = static_cast<unsigned long long>(312968);
0085   end = static_cast<unsigned long long>(313120);
0086
0087   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising Noise plots " << std::endl;
0088
0089   SiStripNoiseValuePerAPV histo8;
0090   histo8.process(connectionString, PI::mk_input(tag, start, start));
0091   edm::LogPrint("testSiStripPayloadInspector") << histo8.data() << std::endl;
0092
0093   SiStripNoiseValueComparisonPerAPVSingleTag histo9;
0094   histo9.process(connectionString, PI::mk_input(tag, start, end));
0095   edm::LogPrint("testSiStripPayloadInspector") << histo9.data() << std::endl;
0096
0097   SiStripNoiseComparatorMeanByRegionSingleTag histoCompareMeanByRegion;
0098   histoCompareMeanByRegion.process(connectionString, PI::mk_input(tag, start, start));
0099   edm::LogPrint("testSiStripPayloadInspector") << histoCompareMeanByRegion.data() << std::endl;
0100
0101   SiStripNoisePerDetId histoNoiseForDetId;
0102   inputs["DetIds"] = "470148232";
0103   histoNoiseForDetId.setInputParamValues(inputs);
0104   histoNoiseForDetId.process(connectionString, PI::mk_input(tag, start, start));
0105   edm::LogPrint("testSiStripPayloadInspector") << histoNoiseForDetId.data() << std::endl;
0106
0107   SiStripNoiseCorrelationByPartition histoNoiseCorrelationByPartition;
0108   histoNoiseCorrelationByPartition.process(connectionString, PI::mk_input(tag, start, start));
0109   edm::LogPrint("testSiStripPayloadInspector") << histoNoiseCorrelationByPartition.data() << std::endl;
0110
0111   SiStripNoiseDiffByPartitionSingleTag histoNoiseDiffByPart;
0112   histoNoiseDiffByPart.process(connectionString, PI::mk_input(tag, start, end));
0113   edm::LogPrint("testSiStripPayloadInspector") << histoNoiseDiffByPart.data() << std::endl;
0114
0115   // Pedestals
0116
0117   tag = "SiStripPedestals_v2_prompt";
0118   start = static_cast<unsigned long long>(303420);
0119   end = static_cast<unsigned long long>(313120);
0120
0121   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising Pedestal plots " << std::endl;
0122
0123   SiStripPedestalValuePerStrip histo10;
0124   histo10.process(connectionString, PI::mk_input(tag, start, start));
0125   edm::LogPrint("testSiStripPayloadInspector") << histo10.data() << std::endl;
0126
0127   SiStripPedestalValueComparisonPerModuleSingleTag histo11;
0128   histo11.process(connectionString, PI::mk_input(tag, start, end));
0129   edm::LogPrint("testSiStripPayloadInspector") << histo11.data() << std::endl;
0130
0131   SiStripPedestalPerDetId histoPedestalForDetId;
0132   inputs["DetIds"] = py::cast<std::string>(inputs["DetIds"]) +
0133                      ",470065830,369121594,369124670,470177668";  // add a bunch of other DetIds
0134   histoPedestalForDetId.setInputParamValues(inputs);
0135   histoPedestalForDetId.process(connectionString, PI::mk_input(tag, start, start));
0136   edm::LogPrint("testSiStripPayloadInspector") << histoPedestalForDetId.data() << std::endl;
0137
0138   SiStripPedestalCorrelationByPartition histoPedestalCorrelationByPartition;
0139   histoPedestalCorrelationByPartition.process(connectionString, PI::mk_input(tag, start, start));
0140   edm::LogPrint("testSiStripPayloadInspector") << histoPedestalCorrelationByPartition.data() << std::endl;
0141
0142   SiStripPedestalDiffByPartitionSingleTag histoPedDiffByPart;
0143   histoPedDiffByPart.process(connectionString, PI::mk_input(tag, start, end));
0144   edm::LogPrint("testSiStripPayloadInspector") << histoPedDiffByPart.data() << std::endl;
0145
0146   // Latency
0147
0148   tag = "SiStripLatency_v2_prompt";
0149   start = static_cast<unsigned long long>(315347);
0150   end = static_cast<unsigned long long>(316675);
0151
0152   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising Latency plots " << std::endl;
0153
0154   SiStripLatencyMode histo12;
0155   histo12.process(connectionString, PI::mk_input(tag, start, start));
0156   edm::LogPrint("testSiStripPayloadInspector") << histo12.data() << std::endl;
0157
0158   SiStripLatencyModeHistory histo13;
0159   histo13.process(connectionString, PI::mk_input(tag, start, end));
0160   edm::LogPrint("testSiStripPayloadInspector") << histo13.data() << std::endl;
0161
0162   // Threshold
0163   tag = "SiStripThreshold_v1_prompt";
0164   start = static_cast<unsigned long long>(315352);
0165   end = static_cast<unsigned long long>(315460);
0166
0167   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising Threshold plots " << std::endl;
0168
0169   SiStripThresholdValueHigh histo14;
0170   histo14.process(connectionString, PI::mk_input(tag, start, end));
0171   edm::LogPrint("testSiStripPayloadInspector") << histo14.data() << std::endl;
0172
0173   // test SiStripTkMaps
0174   tag = "SiStripBadComponents_startupMC_for2017_v1_mc";
0175   start = static_cast<unsigned long long>(1);
0176   end = static_cast<unsigned long long>(1);
0177
0178   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising SiStripTkMaps plots " << std::endl;
0179
0180   SiStripBadStripFractionTH2PolyTkMap histoTkMap;
0181   histoTkMap.process(connectionString, PI::mk_input(tag, start, end));
0182   edm::LogPrint("testSiStripPayloadInspector") << histoTkMap.data() << std::endl;
0183
0184   // test FED cabling
0185   tag = "SiStripFedCabling_GR10_v1_hlt";
0186   start = static_cast<unsigned long long>(1);
0187   end = static_cast<unsigned long long>(400000);
0188
0189   edm::LogPrint("testSiStripPayloadInspector") << "## Exercising FEDCabling plots " << std::endl;
0190
0191   SiStripFedCablingComparisonTrackerMapSingleTag histoFEDCablingCompare;
0192   histoFEDCablingCompare.process(connectionString, PI::mk_input(tag, start, end));
0193   edm::LogPrint("testSiStripPayloadInspector") << histoFEDCablingCompare.data() << std::endl;
0194
0195   inputs.clear();
0196 #if PY_MAJOR_VERSION >= 3
0197   // TODO I don't know why this Py_INCREF is necessary...
0198   Py_INCREF(inputs.ptr());
0199 #endif
0200
0201   Py_Finalize();
0202 }