Ecal phase2 new digitization 11 2 x

CalibFormats/CaloObjects/src/CaloTSamples.cc

@@ -1,7 +1,9 @@
-#include "CalibFormats/CaloObjects/interface/CaloTSamples.icc"
 #include "CalibFormats/CaloObjects/interface/CaloTSamplesBase.icc"
+#include "CalibFormats/CaloObjects/interface/CaloTSamples.icc"
+#include "CondFormats/EcalObjects/interface/EcalConstants.h"
 
 template class CaloTSamplesBase<float>;
 
-template class CaloTSamples<float, 10>;
+template class CaloTSamples<float, ecalPh2::sampleSize>;
 template class CaloTSamples<float, 3>;
+template class CaloTSamples<float, 10>;

CondFormats/DataRecord/interface/EcalCATIAGainRatiosRcd.h

@@ -0,0 +1,6 @@
+#ifndef CondFormats_DataRecord_EcalCATIAGainRatiosRcd_H
+#define CondFormats_DataRecord_EcalCATIAGainRatiosRcd_H
+
+#include "FWCore/Framework/interface/EventSetupRecordImplementation.h"
+class EcalCATIAGainRatiosRcd : public edm::eventsetup::EventSetupRecordImplementation<EcalCATIAGainRatiosRcd> {};
+#endif

CondFormats/DataRecord/interface/EcalLiteDTUPedestalsRcd.h

@@ -0,0 +1,6 @@
+#ifndef CondFormats_DataRecord_EcalLiteDTUPedestalsRcd_h
+#define CondFormats_DataRecord_EcalLiteDTUPedestalsRcd_h
+
+#include "FWCore/Framework/interface/EventSetupRecordImplementation.h"
+class EcalLiteDTUPedestalsRcd : public edm::eventsetup::EventSetupRecordImplementation<EcalLiteDTUPedestalsRcd> {};
+#endif

CondFormats/DataRecord/src/EcalCATIAGainRatiosRcd.cc

@@ -0,0 +1,4 @@
+#include "CondFormats/DataRecord/interface/EcalCATIAGainRatiosRcd.h"
+#include "FWCore/Framework/interface/eventsetuprecord_registration_macro.h"
+
+EVENTSETUP_RECORD_REG(EcalCATIAGainRatiosRcd);

CondFormats/DataRecord/src/EcalLiteDTUPedestalsRcd.cc

@@ -0,0 +1,4 @@
+#include "CondFormats/DataRecord/interface/EcalLiteDTUPedestalsRcd.h"
+#include "FWCore/Framework/interface/eventsetuprecord_registration_macro.h"
+
+EVENTSETUP_RECORD_REG(EcalLiteDTUPedestalsRcd);

CondFormats/EcalObjects/interface/EcalCATIAGainRatios.h

@@ -0,0 +1,10 @@
+#ifndef CondFormats_EcalObjects_EcalCATIAGainRatios_h
+#define CondFormats_EcalObjects_EcalCATIAGainRatios_h
+
+#include "CondFormats/EcalObjects/interface/EcalCondObjectContainer.h"
+
+typedef float EcalCATIAGainRatio;
+typedef EcalFloatCondObjectContainer EcalCATIAGainRatioMap;
+typedef EcalCATIAGainRatioMap EcalCATIAGainRatios;
+
+#endif

CondFormats/EcalObjects/interface/EcalConstants.h

@@ -0,0 +1,31 @@
+//Namespaces for Phase1 and Phase2
+#ifndef CondFormats_EcalObject_EcalConstants_h
+#define CondFormats_EcalObject_EcalConstants_h
+
+class ecalPh2 {
+public:
+  static constexpr double Samp_Period = 6.25;
+  static constexpr unsigned int NGAINS = 2;
+  static constexpr float gains[NGAINS] = {10., 1.};
+  static constexpr unsigned int gainId1 = 1;
+  static constexpr unsigned int gainId10 = 0;
+  static constexpr unsigned int sampleSize = 16;
+  static constexpr unsigned int NBITS = 12;                 // number of available bits
+  static constexpr unsigned int MAXADC = (1 << NBITS) - 1;  // 2^12 -1,  adc max range
+  static constexpr unsigned int kEBChannels = 61200;
+  static constexpr double maxEneEB = 2000.;
+  static constexpr unsigned int kNOffsets = 2000;
+  static constexpr unsigned int kAdcMask = 0xFFF;
+  static constexpr unsigned int kGainIdMask = 0x3;
+
+};  // namespace ecalPh2
+
+class ecalPh1 {
+public:
+  static constexpr double Samp_Period = 25.;
+  static constexpr unsigned int NGAINS = 4;
+  static constexpr float gains[NGAINS] = {0., 12., 6., 1.};
+  static constexpr unsigned int sampleSize = 10;
+  static constexpr unsigned int kNOffsets = 2000;
+};  // namespace ecalPh1
+#endif

CondFormats/EcalObjects/interface/EcalLiteDTUPedestals.h

@@ -0,0 +1,39 @@
+#ifndef CondFormats_EcalObject_EcalLiteDTUPedestals_h
+#define CondFormats_EcalObject_EcalLiteDTUPedestals_h
+
+#include "CondFormats/Serialization/interface/Serializable.h"
+#include "CondFormats/EcalObjects/interface/EcalCondObjectContainer.h"
+#include "CondFormats/EcalObjects/interface/EcalConstants.h"
+
+class EcalLiteDTUPedestals {
+public:
+  int setMean(unsigned int i, float value) {
+    if (i >= ecalPh2::NGAINS)
+      return -1;
+    else
+      meanarray[i] = value;
+    return 1;
+  }
+
+  int setRMS(unsigned int i, float value) {
+    if (i >= ecalPh2::NGAINS)
+      return -1;
+    else
+      rmsarray[i] = value;
+    return 1;
+  }
+
+  float mean(unsigned int i) const { return meanarray[i]; }
+
+  float rms(unsigned int i) const { return rmsarray[i]; }
+
+private:
+  float meanarray[ecalPh2::NGAINS] = {13., 8.};
+  float rmsarray[ecalPh2::NGAINS] = {2.8, 1.2};
+  COND_SERIALIZABLE;
+};
+
+typedef EcalCondObjectContainer<EcalLiteDTUPedestals> EcalLiteDTUPedestalsMap;
+typedef EcalLiteDTUPedestalsMap::const_iterator EcalLiteDTUPedestalsMapIterator;
+
+#endif

CondFormats/EcalObjects/src/SerializationManual.h

@@ -4,6 +4,7 @@ COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalLaserAPDPNRatios::Eca
 COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalMappingElement>);
 COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalMGPAGainRatio>);
 COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalPedestal>);
+COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalLiteDTUPedestals>);
 COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalPulseShape>);
 COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalPulseCovariance>);
 COND_SERIALIZATION_INSTANTIATE(EcalCondObjectContainer<EcalPulseSymmCovariance>);

CondFormats/EcalObjects/src/T_EventSetup_EcalLiteDTUPedestals.cc

@@ -0,0 +1,4 @@
+#include "CondFormats/EcalObjects/interface/EcalLiteDTUPedestals.h"
+#include "FWCore/Utilities/interface/typelookup.h"
+
+TYPELOOKUP_DATA_REG(EcalLiteDTUPedestalsMap);

CondFormats/EcalObjects/src/classes.h

@@ -1,4 +1,4 @@
-
+#include <cstdint>
 #include "CondFormats/EcalObjects/interface/EcalCondObjectContainer.h"
 #include "CondFormats/EcalObjects/interface/EcalCondTowerObjectContainer.h"
 #include "CondFormats/EcalObjects/interface/EcalPedestals.h"
@@ -64,7 +64,9 @@
 #include "CondFormats/EcalObjects/interface/EcalTPGSpike.h"
 #include "CondFormats/EcalObjects/interface/EcalSRSettings.h"
 #include "CondFormats/EcalObjects/interface/EcalSimPulseShape.h"
-#include <cstdint>
+//ECAL PH2:
+#include "CondFormats/EcalObjects/interface/EcalLiteDTUPedestals.h"
+#include "CondFormats/EcalObjects/interface/EcalCATIAGainRatios.h"
 
 namespace CondFormats_EcalObjects {
   struct dictionary {
@@ -93,6 +95,11 @@ namespace CondFormats_EcalObjects {
     EcalContainer<EBDetId, EcalPedestal> ec_ebDetId_ecalPedestal;
     EcalCondObjectContainer<EcalPedestal> pedmap;  //typedef EcalPedestals
 
+    //ECAL PH2:
+    std::vector<EcalLiteDTUPedestals> v_ecalDTUPedestals;
+    EcalContainer<EBDetId, EcalLiteDTUPedestals> ec_ebDetId_ecalLiteDTUPedestals;
+    EcalCondObjectContainer<EcalLiteDTUPedestals> pedmapDTU;  //typedef EcalPedestals
+
     std::vector<EcalTPGCrystalStatusCode> v_ecalTPGCrystalStatusCode;
     EcalContainer<EEDetId, EcalTPGCrystalStatusCode> ec_eeDetId_ecalTPGCrystalStatusCode;
     EcalContainer<EBDetId, EcalTPGCrystalStatusCode> ec_ebDetId_ecalTPGCrystalStatusCode;

CondFormats/EcalObjects/src/classes_def.xml

@@ -82,6 +82,14 @@
    <field name = "ee_" mapping = "blob"/>
 </class>
 
+<class name="EcalLiteDTUPedestals"/>
+<!-- templates for EcalCondObjectContainer having the class as template parameter -->
+<class name="std::vector<EcalLiteDTUPedestals>"/>
+<class name="EcalContainer<EBDetId,EcalLiteDTUPedestals>"/>
+<class name="EcalCondObjectContainer<EcalLiteDTUPedestals>">
+   <field name = "eb_" mapping = "blob"/>
+</class>
+
 <!-- base class -->
 <class name="EcalTPGCrystalStatusCode"/>
 <!-- templates for EcalCondObjectContainer having the class as template parameter -->

Configuration/Eras/python/Era_Phase2C11_Ecal_Devel_cff.py

@@ -0,0 +1,7 @@
+import FWCore.ParameterSet.Config as cms
+
+from Configuration.Eras.Era_Phase2C11_cff import Phase2C11
+from Configuration.Eras.Modifier_phase2_ecal_devel_cff import phase2_ecal_devel
+
+Phase2C11_Ecal_Devel = cms.ModifierChain(Phase2C11,phase2_ecal_devel)
+

Configuration/Eras/python/Modifier_phase2_ecal_devel_cff.py

@@ -0,0 +1,4 @@
+import FWCore.ParameterSet.Config as cms
+
+phase2_ecal_devel =  cms.Modifier()
+

Configuration/StandardSequences/python/DigiToRaw_cff.py

@@ -55,3 +55,6 @@
 
 from Configuration.Eras.Modifier_fastSim_cff import fastSim
 fastSim.toReplaceWith(DigiToRawTask, DigiToRawTask.copyAndExclude([siPixelRawData,SiStripDigiToRaw,castorRawData,ctppsRawData]))
+
+from Configuration.Eras.Modifier_phase2_ecal_devel_cff import phase2_ecal_devel
+phase2_ecal_devel.toReplaceWith(DigiToRawTask, DigiToRawTask.copyAndExclude([esDigiToRaw]))

Configuration/StandardSequences/python/Digi_cff.py

@@ -83,3 +83,24 @@ def _modifyEnableHcalHardcode( theProcess ):
 
 from Configuration.Eras.Modifier_hcalHardcodeConditions_cff import hcalHardcodeConditions
 modifyEnableHcalHardcode_ = hcalHardcodeConditions.makeProcessModifier( _modifyEnableHcalHardcode )
+
+
+#phase 2 ecal pedestals
+def _modifyEcalPedestals( process ):
+    process.load("SimCalorimetry.EcalSimProducers.esEcalLiteDTUPedestalsProducer_cfi")
+
+    from CondCore.CondDB.CondDB_cfi import CondDB
+    CondDB.connect = cms.string('sqlite_file:SimCalorimetry/EcalSimProducers/data/simPulseShapePhaseII.db')
+    process.ecalConditions = cms.ESSource("PoolDBESSource", CondDB,
+      toGet = cms.VPSet(         
+                  cms.PSet(
+                      record = cms.string('EcalSimPulseShapeRcd') ,
+                      tag = cms.string('EcalSimPulseShape_default_mc')
+                  )
+              )
+        )
+    process.es_prefer_ecalPulseShape = cms.ESPrefer("PoolDBESSource","ecalConditions")
+
+from Configuration.Eras.Modifier_phase2_ecal_devel_cff import phase2_ecal_devel
+modifyDigi_Phase2EcalPed = phase2_ecal_devel.makeProcessModifier(_modifyEcalPedestals)
+

Configuration/StandardSequences/python/Eras.py

@@ -44,6 +44,7 @@ def __init__(self):
                  'Phase2C11_dd4hep',
                  'Phase2C12_dd4hep',
                  'Phase2C11M9',
+                 'Phase2C11_Ecal_Devel'
         ]
 
         internalUseMods = ['run2_common', 'run2_25ns_specific',
@@ -58,7 +59,7 @@ def __init__(self):
                            'phase2_common', 'phase2_tracker',
                            'phase2_muon', 'phase2_GEM', 'phase2_GE0',
                            'phase2_hgcal', 'phase2_timing', 'phase2_hfnose', 'phase2_hgcalV10', 'phase2_hgcalV11', 'phase2_hgcalV12',
-                           'phase2_timing_layer', 'phase2_hcal', 'phase2_ecal',
+                           'phase2_timing_layer', 'phase2_hcal', 'phase2_ecal','phase2_ecal_devel',
                            'phase2_trigger',
                            'trackingLowPU', 'trackingPhase1', 'ctpps_2016', 'ctpps_2017', 'ctpps_2018', 'ctpps_2021', 'trackingPhase2PU140','highBetaStar_2018',
                            'tracker_apv_vfp30_2016', 'pf_badHcalMitigation', 'run2_miniAOD_80XLegacy','run2_miniAOD_94XFall17', 'run2_nanoAOD_92X',

DataFormats/EcalDigi/interface/EBDataFrame.h

@@ -3,6 +3,7 @@
 
 #include "DataFormats/EcalDetId/interface/EBDetId.h"
 #include "DataFormats/EcalDigi/interface/EcalDataFrame.h"
+#include "DataFormats/EcalDigi/interface/EcalDataFrame_Ph2.h"
 #include <iosfwd>
 
 /** \class EBDataFrame

DataFormats/EcalDigi/interface/EcalDataFrame_Ph2.h

@@ -0,0 +1,38 @@
+#ifndef DataFormats_EcalDigi_EcalDataFrame_Ph2_h
+#define DataFormats_EcalDigi_EcalDataFrame_Ph2_h
+#include "DataFormats/EcalDigi/interface/EcalLiteDTUSample.h"
+#include "DataFormats/DetId/interface/DetId.h"
+#include "DataFormats/Common/interface/DataFrame.h"
+#include "CondFormats/EcalObjects/interface/EcalConstants.h"
+
+/** \class EcalDataFrame_Ph2
+
+*/
+class EcalDataFrame_Ph2 {
+public:
+  EcalDataFrame_Ph2() {}
+  EcalDataFrame_Ph2(edm::DataFrame const& iframe) : m_data(iframe) {}
+
+  virtual ~EcalDataFrame_Ph2() {}
+
+  DetId id() const { return m_data.id(); }
+
+  int size() const { return m_data.size(); }
+
+  EcalLiteDTUSample operator[](int i) const { return m_data[i]; }
+  EcalLiteDTUSample sample(int i) const { return m_data[i]; }
+
+  // FIXME (shall we throw??)
+  void setSize(int) {}
+  void setSample(int i, EcalLiteDTUSample sam) { m_data[i] = sam; }
+
+  static constexpr int MAXSAMPLES = ecalPh2::sampleSize;
+
+  edm::DataFrame const& frame() const { return m_data; }
+  edm::DataFrame& frame() { return m_data; }
+
+private:
+  edm::DataFrame m_data;
+};
+
+#endif

DataFormats/EcalDigi/interface/EcalDigiCollections.h

@@ -18,6 +18,29 @@
 #include "DataFormats/DetId/interface/DetId.h"
 #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
 #include "DataFormats/Common/interface/DataFrameContainer.h"
+#include "CondFormats/EcalObjects/interface/EcalConstants.h"
+
+class EcalDigiCollectionPh2 : public edm::DataFrameContainer {
+public:
+  typedef edm::DataFrameContainer::size_type size_type;
+  static const size_type MAXSAMPLES = ecalPh2::sampleSize;
+  explicit EcalDigiCollectionPh2(size_type istride = MAXSAMPLES, int isubdet = 0)
+      : edm::DataFrameContainer(istride, isubdet) {}
+  void swap(DataFrameContainer& other) { this->DataFrameContainer::swap(other); }
+};
+
+class EBDigiCollectionPh2 : public EcalDigiCollectionPh2 {
+public:
+  typedef edm::DataFrameContainer::size_type size_type;
+  typedef EBDataFrame Digi;
+  typedef Digi::key_type DetId;
+
+  EBDigiCollectionPh2(size_type istride = MAXSAMPLES) : EcalDigiCollectionPh2(istride, EcalBarrel) {}
+  void swap(EBDigiCollectionPh2& other) { this->EcalDigiCollectionPh2::swap(other); }
+  void push_back(const Digi& digi) { DataFrameContainer::push_back(digi.id(), digi.frame().begin()); }
+  void push_back(id_type iid) { DataFrameContainer::push_back(iid); }
+  void push_back(id_type iid, data_type const* idata) { DataFrameContainer::push_back(iid, idata); }
+};
 
 class EcalDigiCollection : public edm::DataFrameContainer {
 public:
@@ -89,6 +112,10 @@ inline void swap(EEDigiCollection& lhs, EEDigiCollection& rhs) { lhs.swap(rhs);
 
 inline void swap(ESDigiCollection& lhs, ESDigiCollection& rhs) { lhs.swap(rhs); }
 
+inline void swap(EcalDigiCollectionPh2& lhs, EcalDigiCollectionPh2& rhs) { lhs.swap(rhs); }
+
+inline void swap(EBDigiCollectionPh2& lhs, EBDigiCollectionPh2& rhs) { lhs.swap(rhs); }
+
 typedef edm::SortedCollection<EcalTimeDigi> EcalTimeDigiCollection;
 typedef edm::SortedCollection<EcalTriggerPrimitiveDigi> EcalTrigPrimDigiCollection;
 typedef edm::SortedCollection<EcalEBTriggerPrimitiveDigi> EcalEBTrigPrimDigiCollection;

DataFormats/EcalDigi/interface/EcalLiteDTUSample.h

@@ -0,0 +1,47 @@
+#ifndef DataFormats_EcalDigi_EcalLiteDTUSample_h
+#define DataFormats_EcalDigi_EcalLiteDTUSample_h
+
+#include <iosfwd>
+#include <cstdint>
+#include "CondFormats/EcalObjects/interface/EcalConstants.h"
+
+namespace ecalLiteDTU {
+  typedef uint16_t sample_type;
+
+  /// get the ADC sample (12 bits)
+  constexpr int adc(sample_type sample) { return sample & ecalPh2::kAdcMask; }
+  /// get the gainId (2 bits)
+  constexpr int gainId(sample_type sample) { return (sample >> ecalPh2::NBITS) & ecalPh2::kGainIdMask; }
+  constexpr sample_type pack(int adc, int gainId) {
+    return (adc & ecalPh2::kAdcMask) | ((gainId & ecalPh2::kGainIdMask) << ecalPh2::NBITS);
+  }
+}  // namespace ecalLiteDTU
+
+/** \class EcalLiteDTUSample
+ *  Simple container packer/unpacker for a single sample from the Lite_CATIA electronics
+ *
+ *
+ */
+class EcalLiteDTUSample {
+public:
+  EcalLiteDTUSample() { theSample = 0; }
+  EcalLiteDTUSample(uint16_t data) { theSample = data; }
+  EcalLiteDTUSample(int adc, int gainId);
+
+  /// get the raw word
+  uint16_t raw() const { return theSample; }
+  /// get the ADC sample (12 bits)
+  int adc() const { return theSample & ecalPh2::kAdcMask; }
+  /// get the gainId (2 bits)
+  int gainId() const { return (theSample >> ecalPh2::NBITS) & ecalPh2::kGainIdMask; }
+  /// for streaming
+  uint16_t operator()() const { return theSample; }
+  operator uint16_t() const { return theSample; }
+
+private:
+  uint16_t theSample;
+};
+
+std::ostream& operator<<(std::ostream&, const EcalLiteDTUSample&);
+
+#endif

DataFormats/EcalDigi/src/EcalLiteDTUSample.cc

@@ -0,0 +1,11 @@
+#include "DataFormats/EcalDigi/interface/EcalLiteDTUSample.h"
+#include <iostream>
+
+EcalLiteDTUSample::EcalLiteDTUSample(int adc, int gainId) {
+  theSample = (adc & ecalPh2::kAdcMask) | ((gainId & ecalPh2::kGainIdMask) << ecalPh2::NBITS);
+}
+
+std::ostream& operator<<(std::ostream& s, const EcalLiteDTUSample& samp) {
+  s << "ADC=" << samp.adc() << ", gainId=" << samp.gainId();
+  return s;
+}