Merge pull request #8890 from dildick/for-75X-GEM-CSC-Trigger

GEM-CSC integrated local trigger (part I)

DataFormats/CSCDigi/interface/GEMCSCLCTDigi.h

@@ -0,0 +1,32 @@
+#ifndef CSCDigi_GEMCSCLCTDigi_h
+#define CSCDigi_GEMCSCLCTDigi_h
+
+#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h"
+#include <boost/cstdint.hpp>
+#include <iosfwd>
+
+class GEMCSCLCTDigi 
+{
+ public:
+
+  /// Constructors
+  GEMCSCLCTDigi(const CSCCorrelatedLCTDigi, float);
+  GEMCSCLCTDigi();                               /// default
+
+  /// return track number
+  const CSCCorrelatedLCTDigi& getDigi() const { return digi_; }
+
+  /// return bend
+  int getBend()    const { return bend_; }
+
+  ///Comparison
+  bool operator == (const GEMCSCLCTDigi &) const;
+  bool operator != (const GEMCSCLCTDigi &rhs) const
+    { return !(this->operator==(rhs)); }
+
+ private:
+  CSCCorrelatedLCTDigi digi_;
+  float bend_;
+};
+
+#endif

DataFormats/CSCDigi/interface/GEMCSCLCTDigiCollection.h

@@ -0,0 +1,10 @@
+#ifndef GEMCSCLCTDigi_GEMCSCLCTDigiCollection_h
+#define GEMCSCLCTDigi_GEMCSCLCTDigiCollection_h
+
+#include <DataFormats/MuonDetId/interface/CSCDetId.h>
+#include <DataFormats/CSCDigi/interface/GEMCSCLCTDigi.h>
+#include <DataFormats/MuonData/interface/MuonDigiCollection.h>
+
+typedef MuonDigiCollection<CSCDetId,GEMCSCLCTDigi> GEMCSCLCTDigiCollection;
+
+#endif

DataFormats/CSCDigi/src/GEMCSCLCTDigi.cc

@@ -0,0 +1,18 @@
+#include "DataFormats/CSCDigi/interface/GEMCSCLCTDigi.h"
+#include "FWCore/MessageLogger/interface/MessageLogger.h"
+#include <iostream>
+
+/// Constructors
+GEMCSCLCTDigi::GEMCSCLCTDigi(const CSCCorrelatedLCTDigi digi, float bend) : 
+  digi_(digi),
+  bend_(bend)
+{}
+
+/// Default
+GEMCSCLCTDigi::GEMCSCLCTDigi() {
+}
+
+/// Comparison
+bool GEMCSCLCTDigi::operator==(const GEMCSCLCTDigi &rhs) const {
+  return ( digi_ == rhs.getDigi() && bend_ == rhs.getBend() );
+}

DataFormats/CSCDigi/src/classes.h

@@ -12,6 +12,8 @@
 #include <DataFormats/CSCDigi/interface/CSCALCTDigiCollection.h>
 #include <DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h>
 #include <DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h>
+#include <DataFormats/CSCDigi/interface/GEMCSCLCTDigi.h>
+#include <DataFormats/CSCDigi/interface/GEMCSCLCTDigiCollection.h>
 #include <DataFormats/CSCDigi/interface/CSCCFEBStatusDigi.h>
 #include <DataFormats/CSCDigi/interface/CSCCFEBStatusDigiCollection.h>
 #include <DataFormats/CSCDigi/interface/CSCDMBStatusDigi.h>
@@ -44,6 +46,7 @@ namespace DataFormats_CSCDigi {
   CSCCLCTDigi cCLCTD_;
   CSCALCTDigi cALCTD_;
   CSCCorrelatedLCTDigi cCorLCTD_;
+  GEMCSCLCTDigi gcLCTD_;
   CSCCFEBStatusDigi cCSD_;
   CSCTMBStatusDigi cTMBSD_;
   CSCDCCFormatStatusDigi cDFSD_;
@@ -61,6 +64,7 @@ namespace DataFormats_CSCDigi {
   std::vector<CSCCLCTDigi> vCLCTD_;
   std::vector<CSCALCTDigi> vALCTD_;
   std::vector<CSCCorrelatedLCTDigi> vCorLCTD_;
+  std::vector<GEMCSCLCTDigi> vgcLCTD_;
   std::vector<CSCCFEBStatusDigi>  vCSD_;
   std::vector<CSCTMBStatusDigi>  vTMBSD_;
   std::vector<CSCDCCFormatStatusDigi>  vDFSD_;
@@ -77,6 +81,7 @@ namespace DataFormats_CSCDigi {
   std::vector<std::vector<CSCCLCTDigi> > vvCLCTD_;
   std::vector<std::vector<CSCALCTDigi> > vvALCTD_;
   std::vector<std::vector<CSCCorrelatedLCTDigi> > vvCorLCTD_;
+  std::vector<std::vector<GEMCSCLCTDigi> > vvgcLCTD_;
   std::vector<std::vector<CSCCFEBStatusDigi> >  vvCSD_;
   std::vector<std::vector<CSCTMBStatusDigi> >  vvTMBSD_;
   std::vector<std::vector<CSCDMBStatusDigi> >  vvDMBSD_;
@@ -93,6 +98,7 @@ namespace DataFormats_CSCDigi {
   CSCCLCTDigiCollection clCLCTD_;
   CSCALCTDigiCollection clALCTD_;
   CSCCorrelatedLCTDigiCollection clCorLCTD_;
+  GEMCSCLCTDigiCollection clgcLCTD_;  
   CSCCFEBStatusDigiCollection clCSD_;
   CSCTMBStatusDigiCollection clTMBSD_;
   CSCDCCFormatStatusDigiCollection clDFSD_;
@@ -109,6 +115,7 @@ namespace DataFormats_CSCDigi {
   edm::Wrapper<CSCCLCTDigiCollection> wCLCTD_;
   edm::Wrapper<CSCALCTDigiCollection> wALCTD_;
   edm::Wrapper<CSCCorrelatedLCTDigiCollection> wCorLCTD_;
+  edm::Wrapper<GEMCSCLCTDigiCollection> wgcLCTD_;
   edm::Wrapper<CSCCFEBStatusDigiCollection> wCSD_;
   edm::Wrapper<CSCTMBStatusDigiCollection> wTMBSD_;
   edm::Wrapper<CSCDCCFormatStatusDigiCollection> wDFSD_;

DataFormats/CSCDigi/src/classes_def.xml

@@ -14,6 +14,9 @@
    <class name="CSCCorrelatedLCTDigi" ClassVersion="10">
     <version ClassVersion="10" checksum="1301481852"/>
    </class>
+   <class name="GEMCSCLCTDigi" ClassVersion="10">
+    <version ClassVersion="10" checksum="967689346"/>
+   </class>
    <class name="CSCCLCTDigi" ClassVersion="10">
     <version ClassVersion="10" checksum="3910057547"/>
    </class>
@@ -47,6 +50,7 @@
    <class name="std::vector<CSCStripDigi>"/>
    <class name="std::vector<CSCComparatorDigi>"/>
    <class name="std::vector<CSCCorrelatedLCTDigi>"/>
+   <class name="std::vector<GEMCSCLCTDigi>"/>
    <class name="std::vector<CSCCLCTDigi>"/>
    <class name="std::vector<CSCALCTDigi>"/>
    <class name="std::vector<CSCCFEBStatusDigi>"/>
@@ -62,6 +66,7 @@
    <class name="std::map<CSCDetId,std::vector<CSCStripDigi> >"/>
    <class name="std::map<CSCDetId,std::vector<CSCComparatorDigi> >"/>
    <class name="std::map<CSCDetId,std::vector<CSCCorrelatedLCTDigi> >"/>
+   <class name="std::map<CSCDetId,std::vector<GEMCSCLCTDigi> >"/>
    <class name="std::map<CSCDetId,std::vector<CSCCLCTDigi> >"/>
    <class name="std::map<CSCDetId,std::vector<CSCALCTDigi> >"/>
    <class name="std::map<CSCDetId,std::vector<CSCCFEBStatusDigi> >"/>
@@ -78,6 +83,7 @@
    <class name="std::pair<CSCDetId,std::vector<CSCStripDigi> >"/>
    <class name="std::pair<CSCDetId,std::vector<CSCComparatorDigi> >"/>
    <class name="std::pair<CSCDetId,std::vector<CSCCorrelatedLCTDigi> >"/>
+   <class name="std::pair<CSCDetId,std::vector<GEMCSCLCTDigi> >"/>
    <class name="std::pair<CSCDetId,std::vector<CSCCLCTDigi> >"/>
    <class name="std::pair<CSCDetId,std::vector<CSCALCTDigi> >"/>
    <class name="std::pair<CSCDetId,std::vector<CSCCFEBStatusDigi> >"/>
@@ -95,6 +101,7 @@
    <class name="MuonDigiCollection<CSCDetId,CSCStripDigi>"/>
    <class name="MuonDigiCollection<CSCDetId,CSCComparatorDigi>"/>
    <class name="MuonDigiCollection<CSCDetId,CSCCorrelatedLCTDigi>"/>
+   <class name="MuonDigiCollection<CSCDetId,GEMCSCLCTDigi>"/>
    <class name="MuonDigiCollection<CSCDetId,CSCCLCTDigi>"/>
    <class name="MuonDigiCollection<CSCDetId,CSCALCTDigi>"/>
    <class name="MuonDigiCollection<CSCDetId,CSCCFEBStatusDigi>"/>
@@ -111,6 +118,7 @@
    <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,CSCStripDigi>>"  splitLevel="0"/>
    <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,CSCComparatorDigi>>"  splitLevel="0"/>
    <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,CSCCorrelatedLCTDigi> >" splitLevel="0"/>
+   <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,GEMCSCLCTDigi> >" splitLevel="0"/>
    <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,CSCCLCTDigi> >" splitLevel="0"/>
    <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,CSCALCTDigi> >" splitLevel="0"/>
    <class name="edm::Wrapper<MuonDigiCollection<CSCDetId,CSCCFEBStatusDigi> >" splitLevel="0"/>

DataFormats/GEMDigi/interface/GEMCoPadDigi.h

@@ -16,12 +16,13 @@
 class GEMCoPadDigi{
 
 public:
-  explicit GEMCoPadDigi(GEMPadDigi pad1, GEMPadDigi pad2);
+  explicit GEMCoPadDigi(uint8_t roll, GEMPadDigi pad1, GEMPadDigi pad2);
   GEMCoPadDigi();
 
   bool operator==(const GEMCoPadDigi& digi) const;
   bool operator!=(const GEMCoPadDigi& digi) const;
 
+  int roll() const {return roll_;}
   int pad(int l) const;
   int bx(int l) const;
 
@@ -31,6 +32,7 @@ class GEMCoPadDigi{
   void print() const;
 
 private:
+  uint8_t roll_;
   GEMPadDigi first_;
   GEMPadDigi second_;
 };

DataFormats/GEMDigi/src/GEMCoPadDigi.cc

@@ -1,13 +1,15 @@
 #include "DataFormats/GEMDigi/interface/GEMCoPadDigi.h"
 #include <iostream>
 
-GEMCoPadDigi::GEMCoPadDigi(GEMPadDigi f, GEMPadDigi s):
+GEMCoPadDigi::GEMCoPadDigi(uint8_t roll, GEMPadDigi f, GEMPadDigi s):
+  roll_(roll),
   first_(f),
   second_(s)
 {}
 
 
 GEMCoPadDigi::GEMCoPadDigi():
+  roll_(0),
   first_(GEMPadDigi()),
   second_(GEMPadDigi())
 {}
@@ -16,14 +18,14 @@ GEMCoPadDigi::GEMCoPadDigi():
 // Comparison
 bool GEMCoPadDigi::operator == (const GEMCoPadDigi& digi) const
 {
-  return digi.first() == first_ and digi.second() == second_;
+  return digi.first() == first_ and digi.second() == second_ and digi.roll() == roll_;
 }
 
 
 // Comparison
 bool GEMCoPadDigi::operator != (const GEMCoPadDigi& digi) const
 {
-  return digi.first() != first_ or digi.second() != second_;
+  return digi.first() != first_ or digi.second() != second_ or digi.roll() != roll_;
 }
 
 
@@ -45,12 +47,12 @@ int GEMCoPadDigi::bx(int l) const
 
 void GEMCoPadDigi::print() const
 {
-  std::cout << "Pad1 " << first_.pad() << " bx1 " << first_.bx() 
+  std::cout << "Roll " << roll_ << ", pad1 " << first_.pad() << " bx1 " << first_.bx() 
             << ", Pad2 " << second_.pad() << " bx2 " << second_.bx() << std::endl;
 }
 
 
 std::ostream & operator<<(std::ostream & o, const GEMCoPadDigi& digi)
 {
-  return o << " 1:" << digi.first() << ", 2:" << digi.second();
+  return o << "Roll: " << digi.roll() << " 1:" << digi.first() << ", 2:" << digi.second();
 }

DataFormats/GEMDigi/src/classes_def.xml

@@ -16,8 +16,8 @@
 <class name="MuonDigiCollection<GEMDetId,GEMPadDigi>"/> 
 <class name="edm::Wrapper<MuonDigiCollection<GEMDetId,GEMPadDigi> >" splitLevel="0"/> 
 
-<class name="GEMCoPadDigi" ClassVersion="11">
- <version ClassVersion="11" checksum="1150395405"/>
+<class name="GEMCoPadDigi" ClassVersion="12">
+ <version ClassVersion="12" checksum="2678638706"/>
 </class>
 <class name="std::vector<GEMCoPadDigi>"/>
 <class name="std::map<GEMDetId,std::vector<GEMCoPadDigi> >"/>

L1Trigger/CSCTriggerPrimitives/BuildFile.xml

@@ -1,10 +1,14 @@
 <use   name="CondFormats/CSCObjects"/>
 <use   name="DataFormats/CSCDigi"/>
+<use   name="DataFormats/GEMDigi"/>
+<use   name="DataFormats/RPCDigi"/>
 <use   name="DataFormats/L1CSCTrackFinder"/>
 <use   name="DataFormats/MuonDetId"/>
 <use   name="FWCore/MessageLogger"/>
 <use   name="FWCore/ParameterSet"/>
 <use   name="FWCore/Utilities"/>
+<use   name="Geometry/GEMGeometry"/>
+<use   name="Geometry/RPCGeometry"/>
 <use   name="L1Trigger/CSCCommonTrigger"/>
 <export>
   <lib   name="1"/>

L1Trigger/CSCTriggerPrimitives/plugins/CSCTriggerPrimitivesProducer.cc

@@ -22,13 +22,18 @@
 #include "Geometry/Records/interface/MuonGeometryRecord.h"
 #include "L1Trigger/CSCCommonTrigger/interface/CSCTriggerGeometry.h"
 #include "CondFormats/DataRecord/interface/CSCBadChambersRcd.h"
+#include "Geometry/GEMGeometry/interface/GEMGeometry.h"
+#include "Geometry/RPCGeometry/interface/RPCGeometry.h"
 
 //#include "DataFormats/CSCDigi/interface/CSCComparatorDigiCollection.h"
 //#include "DataFormats/CSCDigi/interface/CSCWireDigiCollection.h"
 #include "DataFormats/CSCDigi/interface/CSCALCTDigiCollection.h"
 #include "DataFormats/CSCDigi/interface/CSCCLCTDigiCollection.h"
 #include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h"
 
+#include "DataFormats/CSCDigi/interface/GEMCSCLCTDigiCollection.h"
+#include "DataFormats/GEMDigi/interface/GEMCoPadDigiCollection.h"
+
 // Configuration via EventSetup
 #include "CondFormats/CSCObjects/interface/CSCDBL1TPParameters.h"
 #include "CondFormats/DataRecord/interface/CSCDBL1TPParametersRcd.h"
@@ -38,32 +43,38 @@ CSCTriggerPrimitivesProducer::CSCTriggerPrimitivesProducer(const edm::ParameterS
 
   // if false, parameters will be read in from DB using EventSetup mechanism
   // else will use all parameters from the config file
-  debugParameters_ = conf.getUntrackedParameter<bool>("debugParameters",false);
+  debugParameters_ = conf.getParameter<bool>("debugParameters");
 
   wireDigiProducer_ = conf.getParameter<edm::InputTag>("CSCWireDigiProducer");
   compDigiProducer_ = conf.getParameter<edm::InputTag>("CSCComparatorDigiProducer");
-  checkBadChambers_ = conf.getUntrackedParameter<bool>("checkBadChambers", true);
-
-  lctBuilder_ = new CSCTriggerPrimitivesBuilder(conf); // pass on the conf
-
+  gemPadDigiProducer_ = conf.existsAs<edm::InputTag>("GEMPadDigiProducer")?conf.getParameter<edm::InputTag>("GEMPadDigiProducer"):edm::InputTag("");
+  rpcDigiProducer_ = conf.existsAs<edm::InputTag>("RPCDigiProducer")?conf.getParameter<edm::InputTag>("RPCDigiProducer"):edm::InputTag("");
+  checkBadChambers_ = conf.getParameter<bool>("checkBadChambers");
+  lctBuilder_.reset( new CSCTriggerPrimitivesBuilder(conf) ); // pass on the conf
+
   wire_token_ = consumes<CSCWireDigiCollection>(wireDigiProducer_);
   comp_token_ = consumes<CSCComparatorDigiCollection>(compDigiProducer_);
+  gem_pad_token_ = consumes<GEMPadDigiCollection>(gemPadDigiProducer_);
+  rpc_digi_token_ = consumes<RPCDigiCollection>(rpcDigiProducer_);
 
   // register what this produces
   produces<CSCALCTDigiCollection>();
   produces<CSCCLCTDigiCollection>();
   produces<CSCCLCTPreTriggerCollection>();
   produces<CSCCorrelatedLCTDigiCollection>();
   produces<CSCCorrelatedLCTDigiCollection>("MPCSORTED");
+  produces<GEMCoPadDigiCollection>();
+  produces<GEMCSCLCTDigiCollection>();
   usesResource("CSCTriggerGeometry");
   consumes<CSCComparatorDigiCollection>(compDigiProducer_);
   consumes<CSCWireDigiCollection>(wireDigiProducer_);
+  consumes<GEMPadDigiCollection>(gemPadDigiProducer_);
+  consumes<RPCDigiCollection>(rpcDigiProducer_);
 }
 
 CSCTriggerPrimitivesProducer::~CSCTriggerPrimitivesProducer() {
   LogDebug("L1CSCTrigger")
     << "deleting trigger primitives after " << iev << " events.";
-  delete lctBuilder_;
 }
 
 //void CSCTriggerPrimitivesProducer::beginRun(const edm::EventSetup& setup) {
@@ -80,6 +91,25 @@ void CSCTriggerPrimitivesProducer::produce(edm::Event& ev,
     edm::ESHandle<CSCGeometry> h;
     setup.get<MuonGeometryRecord>().get(h);
     CSCTriggerGeometry::setGeometry(h);
+    lctBuilder_->setCSCGeometry(&*h);
+
+    edm::ESHandle<GEMGeometry> h_gem;
+    try {
+      setup.get<MuonGeometryRecord>().get(h_gem);
+      lctBuilder_->setGEMGeometry(&*h_gem);
+    } catch (edm::eventsetup::NoProxyException<GEMGeometry>& e) {
+      edm::LogInfo("L1CSCTPEmulatorNoGEMGeometry") 
+	<< "+++ Info: GEM geometry is unavailable. Running CSC-only trigger algorithm. +++\n";
+    }
+
+    edm::ESHandle<RPCGeometry> h_rpc;
+    try {
+      setup.get<MuonGeometryRecord>().get(h_rpc);
+      lctBuilder_->setRPCGeometry(&*h_rpc);
+    } catch (edm::eventsetup::NoProxyException<RPCGeometry>& e) {
+      edm::LogInfo("L1CSCTPEmulatorNoRPCGeometry") 
+	<< "+++ Info: RPC geometry is unavailable. Running CSC-only trigger algorithm. +++\n";
+    }
   }
 
   // Find conditions data for bad chambers.
@@ -109,13 +139,30 @@ void CSCTriggerPrimitivesProducer::produce(edm::Event& ev,
   ev.getByToken(comp_token_, compDigis);
   ev.getByToken(wire_token_, wireDigis);
 
-  // Create empty collections of ALCTs, CLCTs, and correlated LCTs upstream
+
+  const GEMPadDigiCollection *gemPads = nullptr;
+  if (!gemPadDigiProducer_.label().empty()) {
+    edm::Handle<GEMPadDigiCollection> gemPadDigis; 
+    ev.getByToken(gem_pad_token_, gemPadDigis);
+    gemPads = gemPadDigis.product();
+  }
+
+  const RPCDigiCollection *rpcDigis = nullptr;
+  if (!rpcDigiProducer_.label().empty()) {
+    edm::Handle<RPCDigiCollection> rpcs; 
+    ev.getByToken(rpc_digi_token_, rpcs);
+    rpcDigis = rpcs.product();
+  }
+
+ // Create empty collections of ALCTs, CLCTs, and correlated LCTs upstream
   // and downstream of MPC.
   std::auto_ptr<CSCALCTDigiCollection> oc_alct(new CSCALCTDigiCollection);
   std::auto_ptr<CSCCLCTDigiCollection> oc_clct(new CSCCLCTDigiCollection);
   std::auto_ptr<CSCCLCTPreTriggerCollection> oc_pretrig(new CSCCLCTPreTriggerCollection);
   std::auto_ptr<CSCCorrelatedLCTDigiCollection> oc_lct(new CSCCorrelatedLCTDigiCollection);
   std::auto_ptr<CSCCorrelatedLCTDigiCollection> oc_sorted_lct(new CSCCorrelatedLCTDigiCollection);
+  std::auto_ptr<GEMCoPadDigiCollection> oc_gemcopad(new GEMCoPadDigiCollection);
+  std::auto_ptr<GEMCSCLCTDigiCollection> oc_gemcsclct(new GEMCSCLCTDigiCollection);
 
   if (!wireDigis.isValid()) {
     edm::LogWarning("L1CSCTPEmulatorNoInputCollection")
@@ -135,8 +182,8 @@ void CSCTriggerPrimitivesProducer::produce(edm::Event& ev,
   if (wireDigis.isValid() && compDigis.isValid()) {   
     const CSCBadChambers* temp = checkBadChambers_ ? pBadChambers.product() : new CSCBadChambers;
     lctBuilder_->build(temp,
-		       wireDigis.product(), compDigis.product(),
-		       *oc_alct, *oc_clct, *oc_pretrig, *oc_lct, *oc_sorted_lct);
+		       wireDigis.product(), compDigis.product(), gemPads, rpcDigis,
+		       *oc_alct, *oc_clct, *oc_pretrig, *oc_lct, *oc_sorted_lct, *oc_gemcopad, *oc_gemcsclct);
     if (!checkBadChambers_)
       delete temp;
   }
@@ -147,4 +194,6 @@ void CSCTriggerPrimitivesProducer::produce(edm::Event& ev,
   ev.put(oc_pretrig);
   ev.put(oc_lct);
   ev.put(oc_sorted_lct,"MPCSORTED");
+  ev.put(oc_gemcopad);
+  ev.put(oc_gemcsclct);
 }