CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
L1TCaloSummary< INPUT, OUTPUT > Class Template Reference

#include <L1Trigger/L1TCaloSummary/plugins/L1TCaloSummary.cc>

Inheritance diagram for L1TCaloSummary< INPUT, OUTPUT >:
edm::stream::EDProducer<>

Public Member Functions

 L1TCaloSummary (const edm::ParameterSet &)
 
 ~L1TCaloSummary () override=default
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Private Member Functions

void beginRun (edm::Run const &, edm::EventSetup const &) override
 
void print ()
 
void produce (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

double boostedJetPtFactor
 
double caloScaleFactor
 
double eGammaIsolationFactor
 
uint32_t eGammaSeed
 
int fwVersion
 
uint32_t jetSeed
 
UCTLayer1 * layer1
 
hls4mlEmulator::ModelLoader loader
 
std::shared_ptr< hls4mlEmulator::Model > model
 
uint32_t nPumBins
 
std::vector< std::vector< std::vector< uint32_t > > > pumLUT
 
edm::EDGetTokenT< L1CaloRegionCollection > regionToken
 
float tauIsolationFactor
 
uint32_t tauSeed
 
bool verbose
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
using CacheTypes = CacheContexts< T... >
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T... >
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 

Detailed Description

template<class INPUT, class OUTPUT>
class L1TCaloSummary< INPUT, OUTPUT >

Description: The package L1Trigger/L1TCaloSummary is prepared for monitoring the CMS Layer-1 Calorimeter Trigger.

Implementation: It prepares region objects and puts them in the event

Definition at line 73 of file L1TCaloSummary.cc.

Constructor & Destructor Documentation

◆ L1TCaloSummary()

template<class INPUT , class OUTPUT >
L1TCaloSummary< INPUT, OUTPUT >::L1TCaloSummary ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 127 of file L1TCaloSummary.cc.

References L1TCaloSummary< INPUT, OUTPUT >::caloScaleFactor, edm::ParameterSet::getParameter(), l1tTowerCalibrationProducer_cfi::iEta, L1TCaloSummary< INPUT, OUTPUT >::loader, SiStripPI::max, L1TCaloSummary< INPUT, OUTPUT >::model, L1TCaloSummary< INPUT, OUTPUT >::nPumBins, and L1TCaloSummary< INPUT, OUTPUT >::pumLUT.

128  : nPumBins(iConfig.getParameter<unsigned int>("nPumBins")),
129  pumLUT(nPumBins, std::vector<std::vector<uint32_t>>(2, std::vector<uint32_t>(13))),
130  caloScaleFactor(iConfig.getParameter<double>("caloScaleFactor")),
131  jetSeed(iConfig.getParameter<unsigned int>("jetSeed")),
132  tauSeed(iConfig.getParameter<unsigned int>("tauSeed")),
133  tauIsolationFactor(iConfig.getParameter<double>("tauIsolationFactor")),
134  eGammaSeed(iConfig.getParameter<unsigned int>("eGammaSeed")),
135  eGammaIsolationFactor(iConfig.getParameter<double>("eGammaIsolationFactor")),
136  boostedJetPtFactor(iConfig.getParameter<double>("boostedJetPtFactor")),
137  verbose(iConfig.getParameter<bool>("verbose")),
138  fwVersion(iConfig.getParameter<int>("firmwareVersion")),
139  regionToken(consumes<L1CaloRegionCollection>(edm::InputTag("simCaloStage2Layer1Digis"))),
140  loader(hls4mlEmulator::ModelLoader(iConfig.getParameter<string>("CICADAModelVersion"))) {
141  std::vector<double> pumLUTData;
142  char pumLUTString[10];
143  for (uint32_t pumBin = 0; pumBin < nPumBins; pumBin++) {
144  for (uint32_t side = 0; side < 2; side++) {
145  if (side == 0)
146  sprintf(pumLUTString, "pumLUT%2.2dp", pumBin);
147  else
148  sprintf(pumLUTString, "pumLUT%2.2dn", pumBin);
149  pumLUTData = iConfig.getParameter<std::vector<double>>(pumLUTString);
150  for (uint32_t iEta = 0; iEta < std::max((uint32_t)pumLUTData.size(), MaxUCTRegionsEta); iEta++) {
151  pumLUT[pumBin][side][iEta] = (uint32_t)round(pumLUTData[iEta] / caloScaleFactor);
152  }
153  if (pumLUTData.size() != (MaxUCTRegionsEta))
154  edm::LogError("L1TCaloSummary") << "PUM LUT Data size integrity check failed; Expected size = "
155  << MaxUCTRegionsEta << "; Provided size = " << pumLUTData.size()
156  << "; Will use what is provided :(" << std::endl;
157  }
158  }
159  produces<L1JetParticleCollection>("Boosted");
160 
161  //anomaly trigger loading
162  model = loader.load_model();
163  produces<float>("CICADAScore");
164 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
L1TCaloSummary::tauIsolationFactor
float tauIsolationFactor
Definition: L1TCaloSummary.cc:99
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
l1tTowerCalibrationProducer_cfi.iEta
iEta
Definition: l1tTowerCalibrationProducer_cfi.py:60
L1TCaloSummary::regionToken
edm::EDGetTokenT< L1CaloRegionCollection > regionToken
Definition: L1TCaloSummary.cc:107
L1TCaloSummary::model
std::shared_ptr< hls4mlEmulator::Model > model
Definition: L1TCaloSummary.cc:112
L1TCaloSummary::boostedJetPtFactor
double boostedJetPtFactor
Definition: L1TCaloSummary.cc:102
L1TCaloSummary::eGammaIsolationFactor
double eGammaIsolationFactor
Definition: L1TCaloSummary.cc:101
L1TCaloSummary::loader
hls4mlEmulator::ModelLoader loader
Definition: L1TCaloSummary.cc:111
L1TCaloSummary::pumLUT
std::vector< std::vector< std::vector< uint32_t > > > pumLUT
Definition: L1TCaloSummary.cc:93
L1TCaloSummary::caloScaleFactor
double caloScaleFactor
Definition: L1TCaloSummary.cc:95

◆ ~L1TCaloSummary()

template<class INPUT , class OUTPUT >
L1TCaloSummary< INPUT, OUTPUT >::~L1TCaloSummary ( )
overridedefault

Member Function Documentation

◆ beginRun()

template<class INPUT , class OUTPUT >
void L1TCaloSummary< INPUT, OUTPUT >::beginRun ( edm::Run const &  ,
edm::EventSetup const &   
)
inlineoverrideprivate

Definition at line 85 of file L1TCaloSummary.cc.

85 {};

◆ fillDescriptions()

template<class INPUT , class OUTPUT >
void L1TCaloSummary< INPUT, OUTPUT >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 296 of file L1TCaloSummary.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

296  {
297  //The following says we do not know what parameters are allowed so do no validation
298  // Please change this to state exactly what you do use, even if it is no parameters
299  edm::ParameterSetDescription desc;
300  desc.setUnknown();
301  descriptions.addDefault(desc);
302 }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99

◆ print()

template<class INPUT , class OUTPUT >
void L1TCaloSummary< INPUT, OUTPUT >::print ( )
private

◆ produce()

template<class INPUT , class OUTPUT >
void L1TCaloSummary< INPUT, OUTPUT >::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivate

Definition at line 172 of file L1TCaloSummary.cc.

References funct::abs(), simCaloStage2Layer1Summary_cfi::boostedJetPtFactor, simCaloStage2Layer1Summary_cfi::caloScaleFactor, simCaloStage2Layer1Summary_cfi::eGammaIsolationFactor, simCaloStage2Layer1Summary_cfi::eGammaSeed, PVValHelper::eta, beamvalidation::exit(), amcDumpToRaw_cfi::fwVersion, g, mps_fire::i, l1tTowerCalibrationProducer_cfi::iEta, iEvent, simCaloStage2Layer1Summary_cfi::jetSeed, EgHLTOffHistBins_cfi::mass, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, eostools::move(), resolutioncreator_cfi::object, DiDispStaMuonMonitor_cfi::pt, nano_mu_digi_cff::region, l1CaloSummaryTree_cfi::regionToken, submitPVValidationJobs::t, simCaloStage2Layer1Summary_cfi::tauIsolationFactor, and simCaloStage2Layer1Summary_cfi::tauSeed.

172  {
173  using namespace edm;
174 
175  std::unique_ptr<L1JetParticleCollection> bJetCands(new L1JetParticleCollection);
176 
177  std::unique_ptr<float> CICADAScore = std::make_unique<float>();
178 
179  UCTGeometry g;
180 
181  // Here we read region data from the region collection created by L1TCaloLayer1 instead of
182  // independently creating regions from TPGs for processing by the summary card. This results
183  // in a single region vector of size 252 whereas from independent creation we had 3*6 vectors
184  // of size 7*2. Indices are mapped in UCTSummaryCard accordingly.
185  UCTSummaryCard summaryCard =
186  UCTSummaryCard(&pumLUT, jetSeed, tauSeed, tauIsolationFactor, eGammaSeed, eGammaIsolationFactor);
187  std::vector<UCTRegion*> inputRegions;
188  inputRegions.clear();
189  edm::Handle<std::vector<L1CaloRegion>> regionCollection;
190  if (!iEvent.getByToken(regionToken, regionCollection))
191  edm::LogError("L1TCaloSummary") << "UCT: Failed to get regions from region collection!";
192  iEvent.getByToken(regionToken, regionCollection);
193  //Model input
194  //This is done as a flat vector input, but future versions may involve 2D input
195  //This will have to be handled later
196  INPUT modelInput[252];
197  for (const L1CaloRegion& i : *regionCollection) {
198  UCTRegionIndex r = g.getUCTRegionIndexFromL1CaloRegion(i.gctEta(), i.gctPhi());
199  UCTTowerIndex t = g.getUCTTowerIndexFromL1CaloRegion(r, i.raw());
200  uint32_t absCaloEta = std::abs(t.first);
201  uint32_t absCaloPhi = std::abs(t.second);
202  bool negativeEta = false;
203  if (t.first < 0)
204  negativeEta = true;
205  uint32_t crate = g.getCrate(t.first, t.second);
206  uint32_t card = g.getCard(t.first, t.second);
207  uint32_t region = g.getRegion(absCaloEta, absCaloPhi);
208  UCTRegion* test = new UCTRegion(crate, card, negativeEta, region, fwVersion);
209  test->setRegionSummary(i.raw());
210  inputRegions.push_back(test);
211  //This *should* fill the tensor in the proper order to be fed to the anomaly model
212  //We take 4 off of the GCT eta/iEta.
213  //iEta taken from this ranges from 4-17, (I assume reserving lower and higher for forward regions)
214  //So our first index, index 0, is technically iEta=4, and so-on.
215  //CICADA reads this as a flat vector
216  modelInput[14 * i.gctPhi() + (i.gctEta() - 4)] = i.et();
217  }
218  //Extract model output
219  OUTPUT modelResult[1] = {
220  OUTPUT("0.0", 10)}; //the 10 here refers to the fact that we read in "0.0" as a decimal number
221  model->prepare_input(modelInput);
222  model->predict();
223  model->read_result(modelResult);
224 
225  *CICADAScore = modelResult[0].to_float();
226 
227  summaryCard.setRegionData(inputRegions);
228 
229  if (!summaryCard.process()) {
230  edm::LogError("L1TCaloSummary") << "UCT: Failed to process summary card" << std::endl;
231  exit(1);
232  }
233 
234  double pt = 0;
235  double eta = -999.;
236  double phi = -999.;
237  double mass = 0;
238 
239  std::list<UCTObject*> boostedJetObjs = summaryCard.getBoostedJetObjs();
240  for (std::list<UCTObject*>::const_iterator i = boostedJetObjs.begin(); i != boostedJetObjs.end(); i++) {
241  const UCTObject* object = *i;
242  pt = ((double)object->et()) * caloScaleFactor * boostedJetPtFactor;
243  eta = g.getUCTTowerEta(object->iEta());
244  phi = g.getUCTTowerPhi(object->iPhi());
245  bitset<3> activeRegionEtaPattern = 0;
246  for (uint32_t iEta = 0; iEta < 3; iEta++) {
247  bool activeStrip = false;
248  for (uint32_t iPhi = 0; iPhi < 3; iPhi++) {
249  if (object->boostedJetRegionET()[3 * iEta + iPhi] > 30 &&
250  object->boostedJetRegionET()[3 * iEta + iPhi] > object->et() * 0.0625)
251  activeStrip = true;
252  }
253  if (activeStrip)
254  activeRegionEtaPattern |= (0x1 << iEta);
255  }
256  bitset<3> activeRegionPhiPattern = 0;
257  for (uint32_t iPhi = 0; iPhi < 3; iPhi++) {
258  bool activeStrip = false;
259  for (uint32_t iEta = 0; iEta < 3; iEta++) {
260  if (object->boostedJetRegionET()[3 * iEta + iPhi] > 30 &&
261  object->boostedJetRegionET()[3 * iEta + iPhi] > object->et() * 0.0625)
262  activeStrip = true;
263  }
264  if (activeStrip)
265  activeRegionPhiPattern |= (0x1 << iPhi);
266  }
267  string regionEta = activeRegionEtaPattern.to_string<char, std::string::traits_type, std::string::allocator_type>();
268  string regionPhi = activeRegionPhiPattern.to_string<char, std::string::traits_type, std::string::allocator_type>();
269 
270  bool centralHighest = object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[0] &&
271  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[1] &&
272  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[2] &&
273  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[3] &&
274  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[5] &&
275  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[6] &&
276  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[7] &&
277  object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[8];
278 
279  if (abs(eta) < 2.5 && ((regionEta == "101" && (regionPhi == "110" || regionPhi == "101" || regionPhi == "010")) ||
280  ((regionEta == "110" || regionEta == "101" || regionEta == "010") && regionPhi == "101") ||
281  (regionEta == "111" && (regionPhi == "110" || regionPhi == "010")) ||
282  ((regionEta == "110" || regionEta == "010") && regionPhi == "111") ||
283  ((regionEta == "010" || regionPhi == "010" || regionEta == "110" || regionPhi == "110" ||
284  regionEta == "011" || regionPhi == "011") &&
285  centralHighest)))
286  bJetCands->push_back(L1JetParticle(math::PtEtaPhiMLorentzVector(pt, eta, phi, mass), L1JetParticle::kCentral));
287  }
288 
289  iEvent.put(std::move(bJetCands), "Boosted");
290  //Write out anomaly score
291  iEvent.put(std::move(CICADAScore), "CICADAScore");
292 }
mps_fire.i
i
Definition: mps_fire.py:429
L1TCaloSummary::tauIsolationFactor
float tauIsolationFactor
Definition: L1TCaloSummary.cc:99
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
g
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
Definition: Activities.doc:4
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
l1tTowerCalibrationProducer_cfi.iEta
iEta
Definition: l1tTowerCalibrationProducer_cfi.py:60
math::PtEtaPhiMLorentzVector
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
iEvent
int iEvent
Definition: GenABIO.cc:224
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
L1TCaloSummary::regionToken
edm::EDGetTokenT< L1CaloRegionCollection > regionToken
Definition: L1TCaloSummary.cc:107
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
L1TCaloSummary::model
std::shared_ptr< hls4mlEmulator::Model > model
Definition: L1TCaloSummary.cc:112
L1TCaloSummary::boostedJetPtFactor
double boostedJetPtFactor
Definition: L1TCaloSummary.cc:102
L1TCaloSummary::eGammaIsolationFactor
double eGammaIsolationFactor
Definition: L1TCaloSummary.cc:101
edm
HLT enums.
Definition: AlignableModifier.h:19
L1TCaloSummary::pumLUT
std::vector< std::vector< std::vector< uint32_t > > > pumLUT
Definition: L1TCaloSummary.cc:93
L1CaloRegion
A calorimeter trigger region (sum of 4x4 trigger towers)
Definition: L1CaloRegion.h:21
nano_mu_digi_cff.region
region
Definition: nano_mu_digi_cff.py:43
L1TCaloSummary::caloScaleFactor
double caloScaleFactor
Definition: L1TCaloSummary.cc:95
resolutioncreator_cfi.object
object
Definition: resolutioncreator_cfi.py:4
eostools.move
def move(src, dest)
Definition: eostools.py:511
beamvalidation.exit
def exit(msg="")
Definition: beamvalidation.py:52

Member Data Documentation

◆ boostedJetPtFactor

template<class INPUT , class OUTPUT >
double L1TCaloSummary< INPUT, OUTPUT >::boostedJetPtFactor
private

Definition at line 102 of file L1TCaloSummary.cc.

◆ caloScaleFactor

template<class INPUT , class OUTPUT >
double L1TCaloSummary< INPUT, OUTPUT >::caloScaleFactor
private

Definition at line 95 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary< INPUT, OUTPUT >::L1TCaloSummary().

◆ eGammaIsolationFactor

template<class INPUT , class OUTPUT >
double L1TCaloSummary< INPUT, OUTPUT >::eGammaIsolationFactor
private

Definition at line 101 of file L1TCaloSummary.cc.

◆ eGammaSeed

template<class INPUT , class OUTPUT >
uint32_t L1TCaloSummary< INPUT, OUTPUT >::eGammaSeed
private

Definition at line 100 of file L1TCaloSummary.cc.

◆ fwVersion

template<class INPUT , class OUTPUT >
int L1TCaloSummary< INPUT, OUTPUT >::fwVersion
private

Definition at line 105 of file L1TCaloSummary.cc.

◆ jetSeed

template<class INPUT , class OUTPUT >
uint32_t L1TCaloSummary< INPUT, OUTPUT >::jetSeed
private

Definition at line 97 of file L1TCaloSummary.cc.

◆ layer1

template<class INPUT , class OUTPUT >
UCTLayer1* L1TCaloSummary< INPUT, OUTPUT >::layer1
private

Definition at line 109 of file L1TCaloSummary.cc.

◆ loader

template<class INPUT , class OUTPUT >
hls4mlEmulator::ModelLoader L1TCaloSummary< INPUT, OUTPUT >::loader
private

Definition at line 111 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary< INPUT, OUTPUT >::L1TCaloSummary().

◆ model

template<class INPUT , class OUTPUT >
std::shared_ptr<hls4mlEmulator::Model> L1TCaloSummary< INPUT, OUTPUT >::model
private

Definition at line 112 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary< INPUT, OUTPUT >::L1TCaloSummary(), and querying.connection::search_everything().

◆ nPumBins

template<class INPUT , class OUTPUT >
uint32_t L1TCaloSummary< INPUT, OUTPUT >::nPumBins
private

Definition at line 91 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary< INPUT, OUTPUT >::L1TCaloSummary().

◆ pumLUT

template<class INPUT , class OUTPUT >
std::vector<std::vector<std::vector<uint32_t> > > L1TCaloSummary< INPUT, OUTPUT >::pumLUT
private

Definition at line 93 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary< INPUT, OUTPUT >::L1TCaloSummary().

◆ regionToken

template<class INPUT , class OUTPUT >
edm::EDGetTokenT<L1CaloRegionCollection> L1TCaloSummary< INPUT, OUTPUT >::regionToken
private

Definition at line 107 of file L1TCaloSummary.cc.

◆ tauIsolationFactor

template<class INPUT , class OUTPUT >
float L1TCaloSummary< INPUT, OUTPUT >::tauIsolationFactor
private

Definition at line 99 of file L1TCaloSummary.cc.

◆ tauSeed

template<class INPUT , class OUTPUT >
uint32_t L1TCaloSummary< INPUT, OUTPUT >::tauSeed
private

Definition at line 98 of file L1TCaloSummary.cc.

◆ verbose

template<class INPUT , class OUTPUT >
bool L1TCaloSummary< INPUT, OUTPUT >::verbose
private

Definition at line 104 of file L1TCaloSummary.cc.

Referenced by python.diffProv.difference::list_diff(), python.diffProv.difference::module_diff(), python.diffProv.difference::onefilemodules(), core.TriggerMatchAnalyzer.TriggerMatchAnalyzer::process(), core.SkimAnalyzerCount.SkimAnalyzerCount::process(), and objects.VertexAnalyzer.VertexAnalyzer::process().

Generated for CMSSW Reference Manual by doxygen 1.8.14