L1T ZDC and uGT Emulator to CMSSW_13_3_X

DataFormats/L1TGlobal/interface/GlobalObject.h

@@ -43,6 +43,8 @@ namespace l1t {
     gtCentrality6,
     gtCentrality7,
     gtExternal,
+    gtZDCP,
+    gtZDCM,
     ObjNull
   };
 

DataFormats/L1TGlobal/src/GlobalObject.cc

@@ -53,6 +53,8 @@ l1t::GlobalObject l1TGtObjectStringToEnum(const std::string& label) {
                                                                             {"CENT6", gtCentrality6},
                                                                             {"CENT7", gtCentrality7},
                                                                             {"External", gtExternal},
+                                                                            {"ZDCP", gtZDCP},
+                                                                            {"ZDCM", gtZDCM},
                                                                             {"ObjNull", ObjNull},
                                                                             {nullptr, (GlobalObject)-1}};
 
@@ -200,6 +202,14 @@ std::string l1t::l1TGtObjectEnumToString(const GlobalObject& gtObject) {
       gtObjectString = "External";
     } break;
 
+    case gtZDCP: {
+      gtObjectString = "ZDCP";
+    } break;
+
+    case gtZDCM: {
+      gtObjectString = "ZDCM";
+    } break;
+
     case ObjNull: {
       gtObjectString = "ObjNull";
       edm::LogInfo("L1TGlobal") << "\n  ObjNull means no valid GlobalObject defined!";

DataFormats/L1Trigger/interface/EtSum.h

@@ -52,6 +52,8 @@ namespace l1t {
       kAsymHt,
       kAsymEtHF,
       kAsymHtHF,
+      kZDCP,
+      kZDCM,
       kUninitialized
     };
 

L1Trigger/L1TCalorimeter/interface/CaloParamsHelper.h

@@ -75,7 +75,8 @@ namespace l1t {
       metHFPhiCalibration = 49,
       layer1HCalFBUpper = 50,
       layer1HCalFBLower = 51,
-      NUM_CALOPARAMNODES = 52
+      hiZDC = 52,
+      NUM_CALOPARAMNODES = 53
     };
 
     CaloParamsHelper() { pnode_.resize(NUM_CALOPARAMNODES); }
@@ -525,6 +526,11 @@ namespace l1t {
     l1t::LUT const* q2LUT() const { return &pnode_[hiQ2].LUT_; }
     void setQ2LUT(const l1t::LUT& lut) { pnode_[hiQ2].LUT_ = lut; }
 
+    // HI ZDC
+    l1t::LUT* zdcLUT() { return &pnode_[hiZDC].LUT_; }
+    l1t::LUT const* zdcLUT() const { return &pnode_[hiZDC].LUT_; }
+    void setZDCLUT(const l1t::LUT& lut) { pnode_[hiZDC].LUT_ = lut; }
+
     // HI parameters
     double etSumCentLower(unsigned centClass) const {
       if (pnode_[etSumCentralityLower].dparams_.size() > centClass)

L1Trigger/L1TCalorimeter/plugins/L1TCaloStage2ParamsESProducer.cc

@@ -344,6 +344,12 @@ L1TCaloStage2ParamsESProducer::L1TCaloStage2ParamsESProducer(const edm::Paramete
   std::shared_ptr<LUT> q2LUT(new LUT(q2LUTStream));
   m_params_helper.setQ2LUT(*q2LUT);
 
+  // HI ZDC calibration LUT for trigger
+  edm::FileInPath zdcLUTFile = conf.getParameter<edm::FileInPath>("zdcLUTFile");
+  std::ifstream zdcLUTStream(zdcLUTFile.fullPath());
+  std::shared_ptr<LUT> zdcLUT(new LUT(zdcLUTStream));
+  m_params_helper.setZDCLUT(*zdcLUT);
+
   // Layer 1 LUT specification
   m_params_helper.setLayer1ECalScaleFactors(conf.getParameter<std::vector<double>>("layer1ECalScaleFactors"));
   m_params_helper.setLayer1HCalScaleFactors(conf.getParameter<std::vector<double>>("layer1HCalScaleFactors"));

L1Trigger/L1TCalorimeter/python/caloParams_cfi.py

@@ -142,7 +142,8 @@
     minimumBiasThresholds = cms.vint32(0, 0, 0, 0),
     centralityLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/centralityLUT_stage1.txt"),
     q2LUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/q2LUT_stage1.txt"),
-
+    zdcLUTFile        = cms.FileInPath("L1Trigger/L1TZDC/data/zdcLUT_HI_v0_1.txt"),
+
     # HCal FB LUT
     layer1HCalFBLUTUpper = cms.vuint32([
     0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 0xBBBABBBA, 

L1Trigger/L1TGlobal/interface/EnergySumZdcCondition.h

@@ -0,0 +1,86 @@
+#ifndef L1Trigger_L1TGlobal_EnergySumZdcCondition_h
+#define L1Trigger_L1TGlobal_EnergySumZdcCondition_h
+
+/**
+ * \class EnergySumZdcCondition
+ * 
+ * 
+ * Description: evaluation of a CondEnergySumZdc condition.
+ * 
+ * Implementation:
+ *    <TODO: enter implementation details>
+ *   
+ * \author: Vasile Mihai Ghete   - HEPHY Vienna 
+ * 
+ *
+ */
+
+// system include files
+#include <iosfwd>
+#include <string>
+
+// user include files
+//   base classes
+#include "L1Trigger/L1TGlobal/interface/ConditionEvaluation.h"
+
+// forward declarations
+class GlobalCondition;
+class EnergySumZdcTemplate;
+
+namespace l1t {
+
+  class L1Candidate;
+
+  class GlobalBoard;
+
+  // class declaration
+  class EnergySumZdcCondition : public ConditionEvaluation {
+  public:
+    /// constructors
+    ///     default
+    EnergySumZdcCondition();
+
+    ///     from base template condition (from event setup usually)
+    EnergySumZdcCondition(const GlobalCondition*, const GlobalBoard*);
+
+    // copy constructor
+    EnergySumZdcCondition(const EnergySumZdcCondition&);
+
+    // destructor
+    ~EnergySumZdcCondition() override;
+
+    // assign operator
+    EnergySumZdcCondition& operator=(const EnergySumZdcCondition&);
+
+  public:
+    /// the core function to check if the condition matches
+    const bool evaluateCondition(const int bxEval) const override;
+
+    /// print condition
+    void print(std::ostream& myCout) const override;
+
+  public:
+    ///   get / set the pointer to a L1GtCondition
+    inline const EnergySumZdcTemplate* gtEnergySumZdcTemplate() const { return m_gtEnergySumZdcTemplate; }
+
+    void setGtEnergySumZdcTemplate(const EnergySumZdcTemplate*);
+
+    ///   get / set the pointer to uGt GlobalBoard
+    inline const GlobalBoard* getuGtB() const { return m_uGtB; }
+
+    void setuGtB(const GlobalBoard*);
+
+  private:
+    /// copy function for copy constructor and operator=
+    void copy(const EnergySumZdcCondition& cp);
+
+  private:
+    /// pointer to a EnergySumZdcTemplate
+    const EnergySumZdcTemplate* m_gtEnergySumZdcTemplate;
+
+    /// pointer to uGt GlobalBoard, to be able to get the trigger objects
+    const GlobalBoard* m_uGtB;
+  };
+
+}  // namespace l1t
+#endif

L1Trigger/L1TGlobal/interface/EnergySumZdcTemplate.h

@@ -0,0 +1,79 @@
+#ifndef L1Trigger_L1TGlobal_EnergySumZdcTemplate_h
+#define L1Trigger_L1TGlobal_EnergySumZdcTemplate_h
+
+/**
+ * \class EnergySumZdcTemplate
+ *
+ *
+ * Description: L1 Global Trigger energy-sum template.
+ *
+ * Implementation:
+ *    <TODO: enter implementation details>
+ *
+ * \author: Vasile Mihai Ghete - HEPHY Vienna
+ *
+ * $Date$
+ * $Revision$
+ *
+ */
+
+// system include files
+#include <string>
+#include <iosfwd>
+
+// user include files
+
+//   base class
+#include "L1Trigger/L1TGlobal/interface/GlobalCondition.h"
+
+// forward declarations
+
+// class declaration
+class EnergySumZdcTemplate : public GlobalCondition {
+public:
+  // constructor
+  EnergySumZdcTemplate();
+
+  // constructor
+  EnergySumZdcTemplate(const std::string&);
+
+  // constructor
+  EnergySumZdcTemplate(const std::string&, const l1t::GtConditionType&);
+
+  // copy constructor
+  EnergySumZdcTemplate(const EnergySumZdcTemplate&);
+
+  // destructor
+  ~EnergySumZdcTemplate() override;
+
+  // assign operator
+  EnergySumZdcTemplate& operator=(const EnergySumZdcTemplate&);
+
+public:
+  struct ObjectParameter {
+    unsigned int etLowThreshold;
+    unsigned int etHighThreshold;
+  };
+
+public:
+  inline const std::vector<ObjectParameter>* objectParameter() const { return &m_objectParameter; }
+
+  /// set functions
+  void setConditionParameter(const std::vector<ObjectParameter>&);
+
+  /// print the condition
+  void print(std::ostream& myCout) const override;
+
+  /// output stream operator
+  friend std::ostream& operator<<(std::ostream&, const EnergySumZdcTemplate&);
+
+private:
+  /// copy function for copy constructor and operator=
+  void copy(const EnergySumZdcTemplate& cp);
+
+private:
+  /// variables containing the parameters
+  std::vector<ObjectParameter> m_objectParameter;
+};
+
+#endif

L1Trigger/L1TGlobal/interface/GlobalBoard.h

@@ -68,13 +68,15 @@ namespace l1t {
                                const edm::EDGetTokenT<BXVector<l1t::Tau>>&,
                                const edm::EDGetTokenT<BXVector<l1t::Jet>>&,
                                const edm::EDGetTokenT<BXVector<l1t::EtSum>>&,
+                               const edm::EDGetTokenT<BXVector<l1t::EtSum>>&,
                                const bool receiveEG,
                                const int nrL1EG,
                                const bool receiveTau,
                                const int nrL1Tau,
                                const bool receiveJet,
                                const int nrL1Jet,
-                               const bool receiveEtSums);
+                               const bool receiveEtSums,
+                               const bool receiveEtSumsZdc);
 
     void receiveMuonObjectData(const edm::Event&,
                                const edm::EDGetTokenT<BXVector<l1t::Muon>>&,
@@ -161,10 +163,13 @@ namespace l1t {
     /// pointer to Tau data list
     inline const BXVector<const l1t::L1Candidate*>* getCandL1Tau() const { return m_candL1Tau; }
 
-    /// pointer to Tau data list
+    /// pointer to EtSum data list
     inline const BXVector<const l1t::EtSum*>* getCandL1EtSum() const { return m_candL1EtSum; }
 
-    /// pointer to Tau data list
+    /// pointer to ZDC EtSum data list
+    inline const BXVector<const l1t::EtSum*>* getCandL1EtSumZdc() const { return m_candL1EtSumZdc; }
+
+    /// pointer to External data list
     inline const BXVector<const GlobalExtBlk*>* getCandL1External() const { return m_candL1External; }
 
     /*  Drop individual EtSums for Now
@@ -226,6 +231,7 @@ namespace l1t {
     BXVector<const l1t::L1Candidate*>* m_candL1Tau;
     BXVector<const l1t::L1Candidate*>* m_candL1Jet;
     BXVector<const l1t::EtSum*>* m_candL1EtSum;
+    BXVector<const l1t::EtSum*>* m_candL1EtSumZdc;
     BXVector<const GlobalExtBlk*>* m_candL1External;
 
     //    BXVector<const l1t::EtSum*>* m_candETM;

L1Trigger/L1TGlobal/interface/GlobalDefinitions.h

@@ -79,6 +79,8 @@ namespace l1t {
     TypeAsymHt,
     TypeAsymEtHF,
     TypeAsymHtHF,
+    TypeZDCP,
+    TypeZDCM,
     GtConditionTypeInvalid = -1
   };
 
@@ -101,6 +103,7 @@ namespace l1t {
     CondCorrelationWithOverlapRemoval,
     CondCorrelationThreeBody,
     CondMuonShower,
+    CondEnergySumZdc,
     GtConditionCategoryInvalid = -1
   };
 

L1Trigger/L1TGlobal/interface/TriggerMenu.h

@@ -35,6 +35,7 @@
 #include "L1Trigger/L1TGlobal/interface/MuonShowerTemplate.h"
 #include "L1Trigger/L1TGlobal/interface/CaloTemplate.h"
 #include "L1Trigger/L1TGlobal/interface/EnergySumTemplate.h"
+#include "L1Trigger/L1TGlobal/interface/EnergySumZdcTemplate.h"
 #include "L1Trigger/L1TGlobal/interface/ExternalTemplate.h"
 #include "L1Trigger/L1TGlobal/interface/CorrelationTemplate.h"
 #include "L1Trigger/L1TGlobal/interface/CorrelationThreeBodyTemplate.h"
@@ -57,6 +58,7 @@ class TriggerMenu {
               const std::vector<std::vector<MuonShowerTemplate> >&,
               const std::vector<std::vector<CaloTemplate> >&,
               const std::vector<std::vector<EnergySumTemplate> >&,
+              const std::vector<std::vector<EnergySumZdcTemplate> >&,
               const std::vector<std::vector<ExternalTemplate> >&,
               const std::vector<std::vector<CorrelationTemplate> >&,
               const std::vector<std::vector<CorrelationThreeBodyTemplate> >&,
@@ -130,6 +132,14 @@ class TriggerMenu {
 
   void setVecEnergySumTemplate(const std::vector<std::vector<EnergySumTemplate> >&);
 
+  //
+  inline const std::vector<std::vector<EnergySumZdcTemplate> >& vecEnergySumZdcTemplate() const {
+    return m_vecEnergySumZdcTemplate;
+  }
+
+  void setVecEnergySumZdcTemplate(const std::vector<std::vector<EnergySumZdcTemplate> >&);
+
+  //
   inline const std::vector<std::vector<ExternalTemplate> >& vecExternalTemplate() const {
     return m_vecExternalTemplate;
   }
@@ -230,6 +240,7 @@ class TriggerMenu {
   std::vector<std::vector<MuonShowerTemplate> > m_vecMuonShowerTemplate;
   std::vector<std::vector<CaloTemplate> > m_vecCaloTemplate;
   std::vector<std::vector<EnergySumTemplate> > m_vecEnergySumTemplate;
+  std::vector<std::vector<EnergySumZdcTemplate> > m_vecEnergySumZdcTemplate;
 
   std::vector<std::vector<ExternalTemplate> > m_vecExternalTemplate;