CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
L1RCTInputProducer Class Reference

#include <L1RCTInputProducer.h>

Inheritance diagram for L1RCTInputProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 L1RCTInputProducer (const edm::ParameterSet &ps)
 
void produce (edm::Event &e, const edm::EventSetup &c) override
 
 ~L1RCTInputProducer () override
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducer () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

edm::InputTag ecalDigisLabel
 
edm::InputTag hcalDigisLabel
 
L1RCT * rct
 
L1RCTLookupTables * rctLookupTables
 
bool useEcal
 
bool useHcal
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 19 of file L1RCTInputProducer.h.

Constructor & Destructor Documentation

L1RCTInputProducer::L1RCTInputProducer ( const edm::ParameterSet &  ps)
explicit

Definition at line 30 of file L1RCTInputProducer.cc.

31  : rctLookupTables(new L1RCTLookupTables), rct(new L1RCT(rctLookupTables)),
32  useEcal(conf.getParameter<bool>("useEcal")),
33  useHcal(conf.getParameter<bool>("useHcal")),
34  ecalDigisLabel(conf.getParameter<edm::InputTag>("ecalDigisLabel")),
35  hcalDigisLabel(conf.getParameter<edm::InputTag>("hcalDigisLabel")) {
36  produces<std::vector<unsigned short>>("rctCrate");
37  produces<std::vector<unsigned short>>("rctCard");
38  produces<std::vector<unsigned short>>("rctTower");
39  produces<std::vector<unsigned int>>("rctEGammaET");
40  produces<std::vector<bool>>("rctHoEFGVetoBit");
41  produces<std::vector<unsigned int>>("rctJetMETET");
42  produces<std::vector<bool>>("rctTowerActivityBit");
43  produces<std::vector<bool>>("rctTowerMIPBit");
44  produces<std::vector<unsigned short>>("rctHFCrate");
45  produces<std::vector<unsigned short>>("rctHFRegion");
46  produces<std::vector<unsigned int>>("rctHFET");
47  produces<std::vector<bool>>("rctHFFG");
48 }
L1RCTInputProducer::ecalDigisLabel
edm::InputTag ecalDigisLabel
Definition: L1RCTInputProducer.h:30
L1RCTInputProducer::hcalDigisLabel
edm::InputTag hcalDigisLabel
Definition: L1RCTInputProducer.h:31
L1RCTInputProducer::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTInputProducer.h:26
L1RCT
Definition: L1RCT.h:20
L1RCTInputProducer::~L1RCTInputProducer ( )
override

Definition at line 50 of file L1RCTInputProducer.cc.

References rct, and rctLookupTables.

50  {
51  if (rct != nullptr)
52  delete rct;
53  if (rctLookupTables != nullptr)
54  delete rctLookupTables;
55 }
L1RCTInputProducer::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTInputProducer.h:26

Member Function Documentation

void L1RCTInputProducer::produce ( edm::Event &  e,
const edm::EventSetup &  c 
)
override

Definition at line 57 of file L1RCTInputProducer.cc.

References EnergyCorrector::c, L1RCT::digiInput(), MillePedeFileConverter_cfg::e, digitizers_cfi::ecal, L1RCT::ecalCompressedET(), ecalDigisLabel, L1RCT::ecalFineGrainBit(), edm::EventSetup::get(), h, digitizers_cfi::hcal, L1RCT::hcalCompressedET(), hcalDigisLabel, L1RCT::hfCompressedET(), edm::HandleBase::isValid(), L1RCTLookupTables::lookup(), eostools::move(), edm::ESHandle< T >::product(), alignCSCRings::r, rct, rctLookupTables, alignCSCRings::s, L1RCTLookupTables::setChannelMask(), L1RCTLookupTables::setEcalScale(), L1RCTLookupTables::setHcalScale(), L1RCTLookupTables::setL1CaloEtScale(), L1RCTLookupTables::setRCTParameters(), useEcal, and useHcal.

58  {
59 
60  // Refresh configuration information every event
61  // Hopefully, this does not take too much time
62  // There should be a call back function in future to
63  // handle changes in configuration
64 
65  edm::ESHandle<L1RCTParameters> rctParameters;
66  eventSetup.get<L1RCTParametersRcd>().get(rctParameters);
67  const L1RCTParameters *r = rctParameters.product();
68  edm::ESHandle<L1RCTChannelMask> channelMask;
69  eventSetup.get<L1RCTChannelMaskRcd>().get(channelMask);
70  const L1RCTChannelMask *c = channelMask.product();
71  edm::ESHandle<L1CaloEcalScale> ecalScale;
72  eventSetup.get<L1CaloEcalScaleRcd>().get(ecalScale);
73  const L1CaloEcalScale *e = ecalScale.product();
74  edm::ESHandle<L1CaloHcalScale> hcalScale;
75  eventSetup.get<L1CaloHcalScaleRcd>().get(hcalScale);
76  const L1CaloHcalScale *h = hcalScale.product();
77  edm::ESHandle<L1CaloEtScale> emScale;
78  eventSetup.get<L1EmEtScaleRcd>().get(emScale);
79  const L1CaloEtScale *s = emScale.product();
80 
81  rctLookupTables->setRCTParameters(r);
82  rctLookupTables->setChannelMask(c);
83  rctLookupTables->setHcalScale(h);
84  rctLookupTables->setEcalScale(e);
85  rctLookupTables->setL1CaloEtScale(s);
86 
87  edm::Handle<EcalTrigPrimDigiCollection> ecal;
88  edm::Handle<HcalTrigPrimDigiCollection> hcal;
89 
90  if (useEcal) {
91  event.getByLabel(ecalDigisLabel, ecal);
92  }
93  if (useHcal) {
94  event.getByLabel(hcalDigisLabel, hcal);
95  }
96 
97  EcalTrigPrimDigiCollection ecalColl;
98  HcalTrigPrimDigiCollection hcalColl;
99  if (ecal.isValid()) {
100  ecalColl = *ecal;
101  }
102  if (hcal.isValid()) {
103  hcalColl = *hcal;
104  }
105 
106  rct->digiInput(ecalColl, hcalColl);
107 
108  // Stuff to create
109 
110  std::unique_ptr<std::vector<unsigned short>> rctCrate(
111  new std::vector<unsigned short>);
112  std::unique_ptr<std::vector<unsigned short>> rctCard(
113  new std::vector<unsigned short>);
114  std::unique_ptr<std::vector<unsigned short>> rctTower(
115  new std::vector<unsigned short>);
116  std::unique_ptr<std::vector<unsigned int>> rctEGammaET(
117  new std::vector<unsigned int>);
118  std::unique_ptr<std::vector<bool>> rctHoEFGVetoBit(new std::vector<bool>);
119  std::unique_ptr<std::vector<unsigned int>> rctJetMETET(
120  new std::vector<unsigned int>);
121  std::unique_ptr<std::vector<bool>> rctTowerActivityBit(new std::vector<bool>);
122  std::unique_ptr<std::vector<bool>> rctTowerMIPBit(new std::vector<bool>);
123 
124  for (int crate = 0; crate < 18; crate++) {
125  for (int card = 0; card < 7; card++) {
126  for (int tower = 0; tower < 32; tower++) {
127  unsigned short ecalCompressedET =
128  rct->ecalCompressedET(crate, card, tower);
129  unsigned short ecalFineGrainBit =
130  rct->ecalFineGrainBit(crate, card, tower);
131  unsigned short hcalCompressedET =
132  rct->hcalCompressedET(crate, card, tower);
133  unsigned int lutBits =
134  rctLookupTables->lookup(ecalCompressedET, hcalCompressedET,
135  ecalFineGrainBit, crate, card, tower);
136  unsigned int eGammaETCode = lutBits & 0x0000007F;
137  bool hOeFGVetoBit = (lutBits >> 7) & 0x00000001;
138  unsigned int jetMETETCode = (lutBits >> 8) & 0x000001FF;
139  bool activityBit = (lutBits >> 17) & 0x00000001;
140  if (eGammaETCode > 0 || jetMETETCode > 0 || hOeFGVetoBit ||
141  activityBit) {
142  rctCrate->push_back(crate);
143  rctCard->push_back(card);
144  rctTower->push_back(tower);
145  rctEGammaET->push_back(eGammaETCode);
146  rctHoEFGVetoBit->push_back(hOeFGVetoBit);
147  rctJetMETET->push_back(jetMETETCode);
148  rctTowerActivityBit->push_back(activityBit);
149  rctTowerMIPBit->push_back(false); // FIXME: MIP bit is not yet defined
150  }
151  }
152  }
153  }
154 
155  std::unique_ptr<std::vector<unsigned short>> rctHFCrate(
156  new std::vector<unsigned short>);
157  std::unique_ptr<std::vector<unsigned short>> rctHFRegion(
158  new std::vector<unsigned short>);
159  std::unique_ptr<std::vector<unsigned int>> rctHFET(
160  new std::vector<unsigned int>);
161  std::unique_ptr<std::vector<bool>> rctHFFG(new std::vector<bool>);
162  for (int crate = 0; crate < 18; crate++) {
163  for (int hfRegion = 0; hfRegion < 8; hfRegion++) {
164  unsigned short hfCompressedET = rct->hfCompressedET(crate, hfRegion);
165  unsigned int hfETCode =
166  rctLookupTables->lookup(hfCompressedET, crate, 999, hfRegion);
167  if (hfETCode > 0) {
168  rctHFCrate->push_back(crate);
169  rctHFRegion->push_back(hfRegion);
170  rctHFET->push_back(hfETCode);
171  rctHFFG->push_back(false); // FIXME: HF FG is not yet defined
172  }
173  }
174  }
175 
176  // putting stuff back into event
177  event.put(std::move(rctCrate), "rctCrate");
178  event.put(std::move(rctCard), "rctCard");
179  event.put(std::move(rctTower), "rctTower");
180  event.put(std::move(rctEGammaET), "rctEGammaET");
181  event.put(std::move(rctHoEFGVetoBit), "rctHoEFGVetoBit");
182  event.put(std::move(rctJetMETET), "rctJetMETET");
183  event.put(std::move(rctTowerActivityBit), "rctTowerActivityBit");
184  event.put(std::move(rctTowerMIPBit), "rctTowerMIPBit");
185  event.put(std::move(rctHFCrate), "rctHFCrate");
186  event.put(std::move(rctHFRegion), "rctHFRegion");
187  event.put(std::move(rctHFET), "rctHFET");
188  event.put(std::move(rctHFFG), "rctHFFG");
189 }
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
L1RCT::hcalCompressedET
unsigned short hcalCompressedET(int crate, int card, int tower)
Definition: L1RCT.h:74
L1RCT::hfCompressedET
unsigned short hfCompressedET(int crate, int tower)
Definition: L1RCT.h:80
L1RCTInputProducer::ecalDigisLabel
edm::InputTag ecalDigisLabel
Definition: L1RCTInputProducer.h:30
L1RCTLookupTables::setHcalScale
void setHcalScale(const L1CaloHcalScale *hcalScale)
Definition: L1RCTLookupTables.h:35
L1RCT::ecalFineGrainBit
unsigned short ecalFineGrainBit(int crate, int card, int tower)
Definition: L1RCT.h:71
L1RCT::ecalCompressedET
unsigned short ecalCompressedET(int crate, int card, int tower)
Definition: L1RCT.h:68
digitizers_cfi.hcal
hcal
Definition: digitizers_cfi.py:25
digitizers_cfi.ecal
ecal
Definition: digitizers_cfi.py:22
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
L1RCTLookupTables::setL1CaloEtScale
void setL1CaloEtScale(const L1CaloEtScale *etScale)
Definition: L1RCTLookupTables.h:39
L1RCTInputProducer::hcalDigisLabel
edm::InputTag hcalDigisLabel
Definition: L1RCTInputProducer.h:31
L1RCTLookupTables::setEcalScale
void setEcalScale(const L1CaloEcalScale *ecalScale)
Definition: L1RCTLookupTables.h:41
L1RCTLookupTables::setRCTParameters
void setRCTParameters(const L1RCTParameters *rctParameters)
Definition: L1RCTLookupTables.h:23
L1RCT::digiInput
void digiInput(const EcalTrigPrimDigiCollection &ecalCollection, const HcalTrigPrimDigiCollection &hcalCollection)
Definition: L1RCT.cc:109
L1RCTInputProducer::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTInputProducer.h:26
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
eostools.move
def move(src, dest)
Definition: eostools.py:511
L1RCTLookupTables::setChannelMask
void setChannelMask(const L1RCTChannelMask *channelMask)
Definition: L1RCTLookupTables.h:27
L1RCTLookupTables::lookup
unsigned int lookup(unsigned short ecalInput, unsigned short hcalInput, unsigned short fgbit, unsigned short crtNo, unsigned short crdNo, unsigned short twrNo) const
Definition: L1RCTLookupTables.cc:23

Member Data Documentation

edm::InputTag L1RCTInputProducer::ecalDigisLabel
private

Definition at line 30 of file L1RCTInputProducer.h.

Referenced by produce().

edm::InputTag L1RCTInputProducer::hcalDigisLabel
private

Definition at line 31 of file L1RCTInputProducer.h.

Referenced by produce().

L1RCT* L1RCTInputProducer::rct
private

Definition at line 27 of file L1RCTInputProducer.h.

Referenced by produce(), and ~L1RCTInputProducer().

L1RCTLookupTables* L1RCTInputProducer::rctLookupTables
private

Definition at line 26 of file L1RCTInputProducer.h.

Referenced by produce(), and ~L1RCTInputProducer().

bool L1RCTInputProducer::useEcal
private

Definition at line 28 of file L1RCTInputProducer.h.

Referenced by produce().

bool L1RCTInputProducer::useHcal
private

Definition at line 29 of file L1RCTInputProducer.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11