Simulation of the pixel bad components on the FED channel basis

CondCore/SiPixelPlugins/src/plugin.cc

@@ -19,6 +19,10 @@
 #include "CondFormats/SiPixelObjects/interface/SiPixelQuality.h"
 #include "CondFormats/DataRecord/interface/SiPixelQualityFromDbRcd.h"
 #include "CondFormats/DataRecord/interface/SiPixelQualityRcd.h"
+#include "CondFormats/SiPixelObjects/interface/SiPixelFEDChannelContainer.h"
+#include "CondFormats/DataRecord/interface/SiPixelStatusScenariosRcd.h"
+#include "CondFormats/SiPixelObjects/interface/SiPixelQualityProbabilities.h"
+#include "CondFormats/DataRecord/interface/SiPixelStatusScenarioProbabilityRcd.h"
 #include "CondFormats/SiPixelObjects/interface/SiPixelCPEGenericErrorParm.h"
 #include "CondFormats/DataRecord/interface/SiPixelCPEGenericErrorParmRcd.h"
 
@@ -61,6 +65,8 @@ REGISTER_PLUGIN(SiPixelDynamicInefficiencyRcd,SiPixelDynamicInefficiency);
 REGISTER_PLUGIN(SiPixelCalibConfigurationRcd,SiPixelCalibConfiguration);
 REGISTER_PLUGIN(SiPixelPerformanceSummaryRcd,SiPixelPerformanceSummary);
 REGISTER_PLUGIN(SiPixelQualityFromDbRcd,SiPixelQuality);
+REGISTER_PLUGIN(SiPixelStatusScenariosRcd,SiPixelFEDChannelContainer);
+REGISTER_PLUGIN(SiPixelStatusScenarioProbabilityRcd,SiPixelQualityProbabilities);
 REGISTER_PLUGIN(SiPixelCPEGenericErrorParmRcd,SiPixelCPEGenericErrorParm);
 REGISTER_PLUGIN(SiPixelTemplateDBObjectRcd,SiPixelTemplateDBObject);
 REGISTER_PLUGIN(SiPixelTemplateDBObject38TRcd,SiPixelTemplateDBObject);

CondCore/Utilities/plugins/Module_2XML.cc

@@ -245,6 +245,8 @@ PAYLOAD_2XML_MODULE( pluginUtilities_payload2xml ){
   PAYLOAD_2XML_CLASS( SiPixelGenErrorDBObject );
   PAYLOAD_2XML_CLASS( SiPixelLorentzAngle );
   PAYLOAD_2XML_CLASS( SiPixelQuality );
+  PAYLOAD_2XML_CLASS( SiPixelFEDChannelContainer );
+  PAYLOAD_2XML_CLASS( SiPixelQualityProbabilities );
   PAYLOAD_2XML_CLASS( SiPixelTemplateDBObject );
   PAYLOAD_2XML_CLASS( SiStripApvGain );
   PAYLOAD_2XML_CLASS( SiStripBackPlaneCorrection );

CondCore/Utilities/src/CondDBFetch.cc

@@ -276,6 +276,8 @@ namespace cond {
       FETCH_PAYLOAD_CASE( SiPixelLorentzAngle )
       FETCH_PAYLOAD_CASE( SiPixelDynamicInefficiency )
       FETCH_PAYLOAD_CASE( SiPixelQuality )
+      FETCH_PAYLOAD_CASE( SiPixelFEDChannelContainer )
+      FETCH_PAYLOAD_CASE( SiPixelQualityProbabilities )
       FETCH_PAYLOAD_CASE( SiPixelTemplateDBObject )
       FETCH_PAYLOAD_CASE( SiPixel2DTemplateDBObject )
       FETCH_PAYLOAD_CASE( SiStripApvGain )

CondCore/Utilities/src/CondDBImport.cc

@@ -296,6 +296,8 @@ namespace cond {
       IMPORT_PAYLOAD_CASE( SiPixelLorentzAngle )
       IMPORT_PAYLOAD_CASE( SiPixelDynamicInefficiency )
       IMPORT_PAYLOAD_CASE( SiPixelQuality )
+      IMPORT_PAYLOAD_CASE( SiPixelFEDChannelContainer )
+      IMPORT_PAYLOAD_CASE( SiPixelQualityProbabilities )
       IMPORT_PAYLOAD_CASE( SiPixelTemplateDBObject )
       IMPORT_PAYLOAD_CASE( SiPixel2DTemplateDBObject )
       IMPORT_PAYLOAD_CASE( SiStripApvGain )

CondCore/Utilities/src/CondFormats.h

@@ -138,6 +138,8 @@
 #include "CondFormats/SiPixelObjects/interface/SiPixelLorentzAngle.h"
 #include "CondFormats/SiPixelObjects/interface/SiPixelDynamicInefficiency.h"
 #include "CondFormats/SiPixelObjects/interface/SiPixelQuality.h"
+#include "CondFormats/SiPixelObjects/interface/SiPixelFEDChannelContainer.h"
+#include "CondFormats/SiPixelObjects/interface/SiPixelQualityProbabilities.h"
 #include "CondFormats/SiPixelObjects/interface/SiPixelTemplateDBObject.h"
 #include "CondFormats/SiPixelObjects/interface/SiPixel2DTemplateDBObject.h"
 #include "CondFormats/SiPixelObjects/interface/SiPixelGenErrorDBObject.h"

CondFormats/DataRecord/interface/SiPixelStatusScenarioProbabilityRcd.h

@@ -0,0 +1,19 @@
+#ifndef SiPixelStatusScenarioProbabilityRcd_SiPixelStatusScenarioProbabilityRcd_h
+#define SiPixelStatusScenarioProbabilityRcd_SiPixelStatusScenarioProbabilityRcd_h
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     SiPixelStatusScenarioProbabilityRcd
+// 
+/**\class SiPixelStatusScenarioProbabilityRcd SiPixelStatusScenarioProbabilityRcd.h CondFormats/DataRecord/interface/SiPixelStatusScenarioProbabilityRcd.h
+ Description: [one line class summary]
+ Usage:
+    <usage>
+*/
+//
+// Author:      Marco Musich
+// Created:     Wed, 30 Nov 2018 11:44:30 GMT
+//
+#include "FWCore/Framework/interface/EventSetupRecordImplementation.h"
+class SiPixelStatusScenarioProbabilityRcd : public edm::eventsetup::EventSetupRecordImplementation<SiPixelStatusScenarioProbabilityRcd> {};
+#endif

CondFormats/DataRecord/interface/SiPixelStatusScenariosRcd.h

@@ -0,0 +1,19 @@
+#ifndef SiPixelStatusScenariosRcd_SiPixelStatusScenariosRcd_h
+#define SiPixelStatusScenariosRcd_SiPixelStatusScenariosRcd_h
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     SiPixelStatusScenariosRcd
+// 
+/**\class SiPixelStatusScenariosRcd SiPixelStatusScenariosRcd.h CondFormats/DataRecord/interface/SiPixelStatusScenariosRcd.h
+ Description: [one line class summary]
+ Usage:
+    <usage>
+*/
+//
+// Author:      Marco Musich
+// Created:     Wed, 28 Nov 2018 10:31:30 GMT
+//
+#include "FWCore/Framework/interface/EventSetupRecordImplementation.h"
+class SiPixelStatusScenariosRcd : public edm::eventsetup::EventSetupRecordImplementation<SiPixelStatusScenariosRcd> {};
+#endif

CondFormats/DataRecord/src/SiPixelStatusScenarioProbabilityRcd.cc

@@ -0,0 +1,13 @@
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     SiPixelStatusScenarioProbabilityRcd
+// 
+// Implementation:
+//     [Notes on implementation]
+//
+// Author:      Marco Musich
+// Created:     Wed, 22 Feb 2017 10:31:30 GMT
+#include "CondFormats/DataRecord/interface/SiPixelStatusScenarioProbabilityRcd.h"
+#include "FWCore/Framework/interface/eventsetuprecord_registration_macro.h"
+EVENTSETUP_RECORD_REG(SiPixelStatusScenarioProbabilityRcd);

CondFormats/DataRecord/src/SiPixelStatusScenariosRcd.cc

@@ -0,0 +1,13 @@
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     SiPixelStatusScenariosRcd
+// 
+// Implementation:
+//     [Notes on implementation]
+//
+// Author:      Marco Musich
+// Created:     Wed, 22 Feb 2017 10:31:30 GMT
+#include "CondFormats/DataRecord/interface/SiPixelStatusScenariosRcd.h"
+#include "FWCore/Framework/interface/eventsetuprecord_registration_macro.h"
+EVENTSETUP_RECORD_REG(SiPixelStatusScenariosRcd);

CondFormats/External/interface/PixelFEDChannel.h

@@ -0,0 +1,69 @@
+#ifndef CondFormats_External_PIXELFEDCHANNEL_H
+#define CondFormats_External_PIXELFEDCHANNEL_H
+
+#include <boost/serialization/base_object.hpp>
+#include <boost/serialization/nvp.hpp>
+#include <boost/serialization/split_free.hpp>
+
+#include <boost/serialization/vector.hpp>
+#include <boost/serialization/string.hpp>
+#include <boost/serialization/map.hpp>
+
+#include "DataFormats/SiPixelDetId/interface/PixelFEDChannel.h"
+
+// struct PixelFEDChannel {
+//   unsigned int fed, link, roc_first, roc_last;
+// };
+
+namespace boost {
+  namespace serialization {
+
+    /*
+     * Note regarding object tracking: all autos used here
+     * must resolve to untracked types, since we use local
+     * variables in the stack which could end up with the same
+     * address. For the moment, all types resolved by auto here
+     * are primitive types, which are untracked by default
+     * by Boost Serialization.
+     */
+
+    template<class Archive>
+    void save(Archive & ar, const PixelFEDChannel & obj, const unsigned int)
+    {
+      auto fed       = obj.fed;
+      auto link      = obj.link;
+      auto roc_first = obj.roc_first;
+      auto roc_last  = obj.roc_last;
+      ar & boost::serialization::make_nvp("fed_"       , fed );
+      ar & boost::serialization::make_nvp("link_"      , link );
+      ar & boost::serialization::make_nvp("roc_first_" , roc_first );  
+      ar & boost::serialization::make_nvp("roc_last_"  , roc_last );  
+    }
+
+    template<class Archive>
+    void load(Archive & ar, PixelFEDChannel & obj, const unsigned int)
+    {
+      unsigned int fed_;
+      unsigned int link_;
+      unsigned int roc_first_;
+      unsigned int roc_last_;
+
+      ar & boost::serialization::make_nvp("fed_"       , fed_ );
+      ar & boost::serialization::make_nvp("link_"      , link_ );
+      ar & boost::serialization::make_nvp("roc_first_" , roc_first_ );  
+      ar & boost::serialization::make_nvp("roc_last_"  , roc_last_ );  
+      PixelFEDChannel tmp{fed_, link_, roc_first_, roc_last_};
+      obj = tmp;
+
+    }
+
+    template<class Archive>
+    void serialize(Archive & ar, PixelFEDChannel & obj, const unsigned int v)
+    {
+      split_free(ar, obj, v);
+    }
+
+  } // namespace serialization
+} // namespace boost
+
+#endif

CondFormats/External/test/testSerializationExternal.cpp

@@ -7,6 +7,7 @@
 #include "CondFormats/External/interface/L1GtLogicParser.h"
 #include "CondFormats/External/interface/SMatrix.h"
 #include "CondFormats/External/interface/Timestamp.h"
+#include "CondFormats/External/interface/PixelFEDChannel.h"
 
 int main()
 {
@@ -19,6 +20,7 @@ int main()
     testSerialization<CLHEP::Hep3Vector>();
     testSerialization<CLHEP::HepEulerAngles>();
     testSerialization<ROOT::Math::SMatrix<double, 2, 3>>();
+    testSerialization<PixelFEDChannel>();
 
     return 0;
 }

CondFormats/SiPixelObjects/interface/SiPixelFEDChannelContainer.h

@@ -0,0 +1,50 @@
+#ifndef CondFormats_SiPixelObjects_SiPixelFEDChannelContainer_h 
+#define CondFormats_SiPixelObjects_SiPixelFEDChannelContainer_h
+
+#include "CondFormats/Serialization/interface/Serializable.h"
+#include "DataFormats/DetId/interface/DetId.h"
+#include "DataFormats/SiPixelDetId/interface/PixelFEDChannel.h"   // N.B. a DataFormat is serialized here (need for dedicated serialization rules, see CondFormats/External/interface/PixelFEDChannel.h)
+
+#include <map>
+#include <string>
+#include <vector>
+
+class SiPixelFEDChannelContainer{
+public:
+  typedef std::map<DetId,std::vector<PixelFEDChannel> > SiPixelFEDChannelCollection;
+  typedef std::unordered_map<std::string, SiPixelFEDChannelCollection> SiPixelBadFEDChannelsScenarioMap;
+
+  SiPixelFEDChannelContainer(){}
+  SiPixelFEDChannelContainer( const SiPixelFEDChannelContainer& rhs ){ m_scenarioMap = rhs.getScenarioMap(); };
+  virtual ~SiPixelFEDChannelContainer(){}
+
+  void setScenario(const std::string &theScenarioId, const SiPixelFEDChannelCollection &theBadFEDChannels);
+
+  const SiPixelBadFEDChannelsScenarioMap& getScenarioMap () const  {return m_scenarioMap;}
+
+  SiPixelFEDChannelCollection getSiPixelBadFedChannels(const std::string &ScenarioId) const;
+  const SiPixelFEDChannelCollection & getSiPixelBadFedChannels(const std::string &ScenarioId);
+
+  const std::vector<PixelFEDChannel>& getSiPixelBadFedChannelsInDetId(const std::string &theScenarioId,DetId theDetId);
+
+  std::unique_ptr<PixelFEDChannelCollection> getDetSetBadPixelFedChannels(const std::string &ScenarioId) const;
+
+  double size()const {return m_scenarioMap.size();}
+  std::vector<std::string> getScenarioList() const;
+
+  void printAll() const;
+
+  //dumping values on output stream
+  void print(std::ostream & os) const;
+
+private:
+
+  SiPixelBadFEDChannelsScenarioMap m_scenarioMap;
+
+  COND_SERIALIZABLE;
+
+};
+
+std::ostream & operator<<( std::ostream &, SiPixelFEDChannelContainer FEDChannels);
+
+#endif //CondFormats_SiPixelObjects_SiPixelFEDChannelContainer_h

CondFormats/SiPixelObjects/interface/SiPixelQualityProbabilities.h

@@ -0,0 +1,45 @@
+#ifndef CondFormats_SiPixelObjects_SiPixelQualityProbabilities_h 
+#define CondFormats_SiPixelObjects_SiPixelQualityProbabilities_h
+
+#include "CondFormats/Serialization/interface/Serializable.h"
+
+#include <map>
+#include <string>
+#include <vector>
+
+class SiPixelQualityProbabilities{
+public:
+  typedef std::vector<std::pair<std::string,float> > probabilityVec; 
+  typedef std::map<unsigned int,probabilityVec> probabilityMap; 
+
+  SiPixelQualityProbabilities(){}
+  SiPixelQualityProbabilities( const SiPixelQualityProbabilities& rhs ){ m_probabilities = rhs.getProbability_Map(); };
+  virtual ~SiPixelQualityProbabilities(){}
+
+  void setProbabilities(const unsigned int puBin, const probabilityVec &theProbabilities);
+
+  const probabilityMap& getProbability_Map () const  {return m_probabilities;}
+
+  probabilityVec   getProbabilities(const unsigned int puBin) const;
+  const probabilityVec & getProbabilities(const unsigned int puBin);
+
+  double size()const {return m_probabilities.size();}
+  double nelements(const int puBin)const {return m_probabilities.at(puBin).size();}
+  std::vector<unsigned int> getPileUpBins() const;
+
+  void printAll() const;
+
+  //dumping values on output stream
+  void print(std::ostream & os) const;
+
+private:
+
+  probabilityMap m_probabilities;
+
+  COND_SERIALIZABLE;
+
+};
+
+std::ostream & operator<<( std::ostream &, SiPixelQualityProbabilities theProbabilities);
+
+#endif //CondFormats_SiPixelObjects_SiPixelQualityProbabilities_h