ECAL SC parameters to EventSetup

CondCore/EcalPlugins/src/plugin.cc

@@ -161,6 +161,11 @@
 #include "CondFormats/EcalObjects/interface/EcalSimPulseShape.h"
 #include "CondFormats/DataRecord/interface/EcalSimPulseShapeRcd.h"
 
+#include "CondFormats/EcalObjects/interface/EcalMustacheSCParameters.h"
+#include "CondFormats/DataRecord/interface/EcalMustacheSCParametersRcd.h"
+#include "CondFormats/EcalObjects/interface/EcalSCDynamicDPhiParameters.h"
+#include "CondFormats/DataRecord/interface/EcalSCDynamicDPhiParametersRcd.h"
+
 REGISTER_PLUGIN(EcalPedestalsRcd, EcalCondObjectContainer<EcalPedestal>);
 REGISTER_PLUGIN(EcalWeightXtalGroupsRcd, EcalCondObjectContainer<EcalXtalGroupId>);
 REGISTER_PLUGIN(EcalTBWeightsRcd, EcalTBWeights);
@@ -222,3 +227,6 @@ REGISTER_PLUGIN(EcalSamplesCorrelationRcd, EcalSamplesCorrelation);
 REGISTER_PLUGIN(EcalPulseShapesRcd, EcalCondObjectContainer<EcalPulseShape>);
 REGISTER_PLUGIN(EcalPulseCovariancesRcd, EcalCondObjectContainer<EcalPulseCovariance>);
 REGISTER_PLUGIN(EcalPulseSymmCovariancesRcd, EcalCondObjectContainer<EcalPulseSymmCovariance>);
+
+REGISTER_PLUGIN(EcalMustacheSCParametersRcd, EcalMustacheSCParameters);
+REGISTER_PLUGIN(EcalSCDynamicDPhiParametersRcd, EcalSCDynamicDPhiParameters);

CondFormats/DataRecord/interface/EcalMustacheSCParametersRcd.h

@@ -0,0 +1,26 @@
+#ifndef CondFormats_DataRecord_EcalMustacheSCParametersRcd_h
+#define CondFormats_DataRecord_EcalMustacheSCParametersRcd_h
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     EcalMustacheSCParametersRcd
+//
+/**\class EcalMustacheSCParametersRcd EcalMustacheSCParametersRcd.h CondFormats/DataRecord/interface/EcalMustacheSCParametersRcd.h
+
+ Description: [one line class summary]
+
+ Usage:
+    <usage>
+
+*/
+//
+// Author:      Thomas Reis
+// Created:     Fri, 23 Oct 2020 09:52:36 GMT
+//
+
+#include "FWCore/Framework/interface/EventSetupRecordImplementation.h"
+
+class EcalMustacheSCParametersRcd
+    : public edm::eventsetup::EventSetupRecordImplementation<EcalMustacheSCParametersRcd> {};
+
+#endif

CondFormats/DataRecord/interface/EcalSCDynamicDPhiParametersRcd.h

@@ -0,0 +1,26 @@
+#ifndef CondFormats_DataRecord_EcalSCDynamicDPhiParametersRcd_h
+#define CondFormats_DataRecord_EcalSCDynamicDPhiParametersRcd_h
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     EcalSCDynamicDPhiParametersRcd
+//
+/**\class EcalSCDynamicDPhiParametersRcd EcalSCDynamicDPhiParametersRcd.h CondFormats/DataRecord/interface/EcalSCDynamicDPhiParametersRcd.h
+
+ Description: [one line class summary]
+
+ Usage:
+    <usage>
+
+*/
+//
+// Author:      Thomas Reis
+// Created:     Wed, 28 Oct 2020 12:56:00 GMT
+//
+
+#include "FWCore/Framework/interface/EventSetupRecordImplementation.h"
+
+class EcalSCDynamicDPhiParametersRcd
+    : public edm::eventsetup::EventSetupRecordImplementation<EcalSCDynamicDPhiParametersRcd> {};
+
+#endif

CondFormats/DataRecord/src/EcalMustacheSCParametersRcd.cc

@@ -0,0 +1,15 @@
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     EcalMustacheSCParametersRcd
+//
+// Implementation:
+//     [Notes on implementation]
+//
+// Author:      Thomas Reis
+// Created:     Fri, 23 Oct 2020 09:52:36 GMT
+
+#include "CondFormats/DataRecord/interface/EcalMustacheSCParametersRcd.h"
+#include "FWCore/Framework/interface/eventsetuprecord_registration_macro.h"
+
+EVENTSETUP_RECORD_REG(EcalMustacheSCParametersRcd);

CondFormats/DataRecord/src/EcalSCDynamicDPhiParametersRcd.cc

@@ -0,0 +1,15 @@
+// -*- C++ -*-
+//
+// Package:     CondFormats/DataRecord
+// Class  :     EcalSCDynamicDPhiParametersRcd
+//
+// Implementation:
+//     [Notes on implementation]
+//
+// Author:      Thomas Reis
+// Created:     Wed, 28 Oct 2020 12:56:00 GMT
+
+#include "CondFormats/DataRecord/interface/EcalSCDynamicDPhiParametersRcd.h"
+#include "FWCore/Framework/interface/eventsetuprecord_registration_macro.h"
+
+EVENTSETUP_RECORD_REG(EcalSCDynamicDPhiParametersRcd);

CondFormats/EcalObjects/interface/EcalMustacheSCParameters.h

@@ -0,0 +1,54 @@
+#ifndef CondFormats_EcalObjects_EcalMustacheSCParameters_h
+#define CondFormats_EcalObjects_EcalMustacheSCParameters_h
+
+#include <iostream>
+#include <vector>
+
+#include "CondFormats/Serialization/interface/Serializable.h"
+
+namespace reco {
+  class MustacheSCParametersHelper;
+}  // namespace reco
+
+class EcalMustacheSCParameters {
+public:
+  EcalMustacheSCParameters(){};
+  virtual ~EcalMustacheSCParameters() = default;
+
+  struct ParabolaParameters {
+    double log10EMin;
+    double etaMin;
+    std::vector<double> pUp;
+    std::vector<double> pLow;
+    std::vector<double> w0Up;
+    std::vector<double> w1Up;
+    std::vector<double> w0Low;
+    std::vector<double> w1Low;
+
+    COND_SERIALIZABLE;
+  };
+
+  float sqrtLogClustETuning() const;
+
+  const ParabolaParameters* parabolaParameters(float log10ClustE, float absSeedEta) const;
+
+  // helper class to set parameters
+  friend class reco::MustacheSCParametersHelper;
+
+  // print parameters to stream:
+  void print(std::ostream&) const;
+  friend std::ostream& operator<<(std::ostream& out, const EcalMustacheSCParameters& params) {
+    params.print(out);
+    return out;
+  }
+
+protected:
+  float sqrtLogClustETuning_;
+
+  // collection is expected to be sorted in ascending ParabolaParameters.log10EMin and ascending ParabolaParameters.etaMax
+  std::vector<ParabolaParameters> parabolaParametersCollection_;
+
+  COND_SERIALIZABLE;
+};
+
+#endif

CondFormats/EcalObjects/interface/EcalSCDynamicDPhiParameters.h

@@ -0,0 +1,50 @@
+#ifndef CondFormats_EcalObjects_EcalSCDynamicDPhiParameters_h
+#define CondFormats_EcalObjects_EcalSCDynamicDPhiParameters_h
+
+#include <iostream>
+#include <vector>
+
+#include "CondFormats/Serialization/interface/Serializable.h"
+
+namespace reco {
+  class SCDynamicDPhiParametersHelper;
+}  // namespace reco
+
+class EcalSCDynamicDPhiParameters {
+public:
+  EcalSCDynamicDPhiParameters(){};
+  virtual ~EcalSCDynamicDPhiParameters() = default;
+
+  struct DynamicDPhiParameters {
+    double eMin;
+    double etaMin;
+    double yoffset;
+    double scale;
+    double xoffset;
+    double width;
+    double saturation;
+    double cutoff;
+
+    COND_SERIALIZABLE;
+  };
+
+  const DynamicDPhiParameters* dynamicDPhiParameters(double clustE, double absSeedEta) const;
+
+  // helper class to set parameters
+  friend class reco::SCDynamicDPhiParametersHelper;
+
+  // print parameters to stream:
+  void print(std::ostream&) const;
+  friend std::ostream& operator<<(std::ostream& out, const EcalSCDynamicDPhiParameters& params) {
+    params.print(out);
+    return out;
+  }
+
+protected:
+  // collection is expected to be sorted in ascending DynamicDPhiParameters.eMin and ascending DynamicDPhiParameters.etaMax
+  std::vector<DynamicDPhiParameters> dynamicDPhiParametersCollection_;
+
+  COND_SERIALIZABLE;
+};
+
+#endif

CondFormats/EcalObjects/src/EcalMustacheSCParameters.cc

@@ -0,0 +1,71 @@
+#include "CondFormats/EcalObjects/interface/EcalMustacheSCParameters.h"
+
+float EcalMustacheSCParameters::sqrtLogClustETuning() const { return sqrtLogClustETuning_; }
+
+const EcalMustacheSCParameters::ParabolaParameters* EcalMustacheSCParameters::parabolaParameters(
+    float log10ClustE, float absSeedEta) const {
+  // assume the collection is lexicographically sorted in ascending ParabolaParameters.log10EMin and ascending ParabolaParameters.etaMin
+  // find the matching log10EMin value
+  auto it1 = std::lower_bound(parabolaParametersCollection_.begin(),
+                              parabolaParametersCollection_.end(),
+                              log10ClustE,
+                              [](const EcalMustacheSCParameters::ParabolaParameters& params, const double var) {
+                                return params.log10EMin < var;
+                              });
+  if (it1 != parabolaParametersCollection_.begin()) {
+    --it1;
+  }
+
+  // find the matching log10EMin and etaMin entry going only up to the sets matching for log10ClustE
+  const auto vars = std::make_pair(it1->log10EMin, absSeedEta);
+  auto it2 = std::lower_bound(
+      parabolaParametersCollection_.begin(),
+      it1 + 1,
+      vars,
+      [](const EcalMustacheSCParameters::ParabolaParameters& params, const std::pair<double, double> vars) {
+        return params.log10EMin < vars.first || params.etaMin < vars.second;
+      });
+
+  return (it2 != parabolaParametersCollection_.begin()) ? &*(it2 - 1) : nullptr;
+}
+
+void EcalMustacheSCParameters::print(std::ostream& out) const {
+  out << "Mustache SC parameters:" << std::endl;
+  out << " sqrtLogClustETuning: " << sqrtLogClustETuning_ << std::endl;
+  out << " Parabola parameters are binned in " << parabolaParametersCollection_.size() << " (log10(E), |eta|) regions."
+      << std::endl;
+  for (const auto& params : parabolaParametersCollection_) {
+    out << " Parameters for log10(E_min)=" << params.log10EMin << " and |eta_min|=" << params.etaMin << ":"
+        << std::endl;
+
+    out << "  pUp:" << std::endl;
+    for (size_t i = 0; i < params.pUp.size(); ++i) {
+      out << "   [" << i << "]: " << params.pUp[i] << std::endl;
+    }
+
+    out << "  pLow:" << std::endl;
+    for (size_t i = 0; i < params.pLow.size(); ++i) {
+      out << "   [" << i << "]: " << params.pLow[i] << std::endl;
+    }
+
+    out << "  w0Up:" << std::endl;
+    for (size_t i = 0; i < params.w0Up.size(); ++i) {
+      out << "   [" << i << "]: " << params.w0Up[i] << std::endl;
+    }
+
+    out << "  w1Up:" << std::endl;
+    for (size_t i = 0; i < params.w1Up.size(); ++i) {
+      out << "   [" << i << "]: " << params.w1Up[i] << std::endl;
+    }
+
+    out << "  w0Low:" << std::endl;
+    for (size_t i = 0; i < params.w0Low.size(); ++i) {
+      out << "   [" << i << "]: " << params.w0Low[i] << std::endl;
+    }
+
+    out << "  w1Low:" << std::endl;
+    for (size_t i = 0; i < params.w1Low.size(); ++i) {
+      out << "   [" << i << "]: " << params.w1Low[i] << std::endl;
+    }
+  }
+}

CondFormats/EcalObjects/src/EcalSCDynamicDPhiParameters.cc

@@ -0,0 +1,42 @@
+#include "CondFormats/EcalObjects/interface/EcalSCDynamicDPhiParameters.h"
+
+const EcalSCDynamicDPhiParameters::DynamicDPhiParameters* EcalSCDynamicDPhiParameters::dynamicDPhiParameters(
+    double clustE, double absSeedEta) const {
+  // assume the collection is lexicographically sorted in ascending DynamicDPhiParams.eMin and ascending DynamicDPhiParams.etaMin
+  // find the matching eMin value
+  auto it1 = std::lower_bound(dynamicDPhiParametersCollection_.begin(),
+                              dynamicDPhiParametersCollection_.end(),
+                              clustE,
+                              [](const EcalSCDynamicDPhiParameters::DynamicDPhiParameters& params, const double var) {
+                                return params.eMin < var;
+                              });
+  if (it1 != dynamicDPhiParametersCollection_.begin()) {
+    --it1;
+  }
+
+  // find the matching eMin and etaMin entry going only up to the sets matching for clustE
+  const auto vars = std::make_pair(it1->eMin, absSeedEta);
+  auto it2 = std::lower_bound(
+      dynamicDPhiParametersCollection_.begin(),
+      it1 + 1,
+      vars,
+      [](const EcalSCDynamicDPhiParameters::DynamicDPhiParameters& params, const std::pair<double, double> vars) {
+        return params.eMin < vars.first || params.etaMin < vars.second;
+      });
+
+  return (it2 != dynamicDPhiParametersCollection_.begin()) ? &*(it2 - 1) : nullptr;
+}
+
+void EcalSCDynamicDPhiParameters::print(std::ostream& out) const {
+  out << "SC dynamic dPhi parameters:" << std::endl;
+  out << " Parameters are binned in " << dynamicDPhiParametersCollection_.size() << " (E, |eta|) regions." << std::endl;
+  for (const auto& params : dynamicDPhiParametersCollection_) {
+    out << " Parameters for E_min=" << params.eMin << " and |eta_min|=" << params.etaMin << ":" << std::endl;
+    out << "  yoffset:    " << params.yoffset << std::endl;
+    out << "  scale:      " << params.scale << std::endl;
+    out << "  xoffset:    " << params.xoffset << std::endl;
+    out << "  width:      " << params.width << std::endl;
+    out << "  saturation: " << params.saturation << std::endl;
+    out << "  cutoff:     " << params.cutoff << std::endl;
+  }
+}

CondFormats/EcalObjects/src/T_EventSetup_EcalMustacheSCParameters.cc

@@ -0,0 +1,5 @@
+// user include files
+#include "CondFormats/EcalObjects/interface/EcalMustacheSCParameters.h"
+#include "FWCore/Utilities/interface/typelookup.h"
+
+TYPELOOKUP_DATA_REG(EcalMustacheSCParameters);

CondFormats/EcalObjects/src/T_EventSetup_EcalSCDynamicDPhiParameters.cc

@@ -0,0 +1,5 @@
+// user include files
+#include "CondFormats/EcalObjects/interface/EcalSCDynamicDPhiParameters.h"
+#include "FWCore/Utilities/interface/typelookup.h"
+
+TYPELOOKUP_DATA_REG(EcalSCDynamicDPhiParameters);