L1TCaloSummary Class Reference

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

Inheritance diagram for L1TCaloSummary:

Public Member Functions
	L1TCaloSummary (const edm::ParameterSet &)
	~L1TCaloSummary () override
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

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 72 of file L1TCaloSummary.cc.

Constructor & Destructor Documentation

L1TCaloSummary()

L1TCaloSummary::L1TCaloSummary ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 125 of file L1TCaloSummary.cc.

References caloScaleFactor, edm::ParameterSet::getParameter(), l1tTowerCalibrationProducer_cfi::iEta, loader, SiStripPI::max, model, nPumBins, and pumLUT.

    : nPumBins(iConfig.getParameter<unsigned int>("nPumBins")),
      pumLUT(nPumBins, std::vector<std::vector<uint32_t>>(2, std::vector<uint32_t>(13))),
      caloScaleFactor(iConfig.getParameter<double>("caloScaleFactor")),
      jetSeed(iConfig.getParameter<unsigned int>("jetSeed")),
      tauSeed(iConfig.getParameter<unsigned int>("tauSeed")),
      tauIsolationFactor(iConfig.getParameter<double>("tauIsolationFactor")),
      eGammaSeed(iConfig.getParameter<unsigned int>("eGammaSeed")),
      eGammaIsolationFactor(iConfig.getParameter<double>("eGammaIsolationFactor")),
      boostedJetPtFactor(iConfig.getParameter<double>("boostedJetPtFactor")),
      verbose(iConfig.getParameter<bool>("verbose")),
      fwVersion(iConfig.getParameter<int>("firmwareVersion")),
      regionToken(consumes<L1CaloRegionCollection>(edm::InputTag("simCaloStage2Layer1Digis"))),
      loader(hls4mlEmulator::ModelLoader(iConfig.getParameter<string>("CICADAModelVersion"))) {
  std::vector<double> pumLUTData;
  char pumLUTString[10];
  for (uint32_t pumBin = 0; pumBin < nPumBins; pumBin++) {
    for (uint32_t side = 0; side < 2; side++) {
      if (side == 0)
        sprintf(pumLUTString, "pumLUT%2.2dp", pumBin);
      else
        sprintf(pumLUTString, "pumLUT%2.2dn", pumBin);
      pumLUTData = iConfig.getParameter<std::vector<double>>(pumLUTString);
      for (uint32_t iEta = 0; iEta < std::max((uint32_t)pumLUTData.size(), MaxUCTRegionsEta); iEta++) {
        pumLUT[pumBin][side][iEta] = (uint32_t)round(pumLUTData[iEta] / caloScaleFactor);
      }
      if (pumLUTData.size() != (MaxUCTRegionsEta))
        edm::LogError("L1TCaloSummary") << "PUM LUT Data size integrity check failed; Expected size = "
                                        << MaxUCTRegionsEta << "; Provided size = " << pumLUTData.size()
                                        << "; Will use what is provided :(" << std::endl;
    }
  }
  produces<L1JetParticleCollection>("Boosted");

  //anomaly trigger loading
  model = loader.load_model();
  produces<float>("anomalyScore");
}

~L1TCaloSummary()

L1TCaloSummary::~L1TCaloSummary ( )

override

Definition at line 164 of file L1TCaloSummary.cc.

{}

Member Function Documentation

beginRun()

void L1TCaloSummary::beginRun ( edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 304 of file L1TCaloSummary.cc.

{}

fillDescriptions()

void L1TCaloSummary::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 331 of file L1TCaloSummary.cc.

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

                                                                                {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

print()

void L1TCaloSummary::print ( void  )

private

Definition at line 294 of file L1TCaloSummary.cc.

{}

produce()

void L1TCaloSummary::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 171 of file L1TCaloSummary.cc.

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

                                                                          {
  using namespace edm;

  std::unique_ptr<L1JetParticleCollection> bJetCands(new L1JetParticleCollection);

  std::unique_ptr<float> anomalyScore = std::make_unique<float>();

  UCTGeometry g;

  // Here we read region data from the region collection created by L1TCaloLayer1 instead of
  // independently creating regions from TPGs for processing by the summary card. This results
  // in a single region vector of size 252 whereas from independent creation we had 3*6 vectors
  // of size 7*2. Indices are mapped in UCTSummaryCard accordingly.
  UCTSummaryCard summaryCard =
      UCTSummaryCard(&pumLUT, jetSeed, tauSeed, tauIsolationFactor, eGammaSeed, eGammaIsolationFactor);
  std::vector<UCTRegion*> inputRegions;
  inputRegions.clear();
  edm::Handle<std::vector<L1CaloRegion>> regionCollection;
  if (!iEvent.getByToken(regionToken, regionCollection))
    edm::LogError("L1TCaloSummary") << "UCT: Failed to get regions from region collection!";
  iEvent.getByToken(regionToken, regionCollection);
  //Model input
  //This is done as a flat vector input, but future versions may involve 2D input
  //This will have to be handled later
  //Would also be good to be able to configure the precision of the ap_fixed type
  ap_ufixed<10, 10> modelInput[252];
  for (const L1CaloRegion& i : *regionCollection) {
    UCTRegionIndex r = g.getUCTRegionIndexFromL1CaloRegion(i.gctEta(), i.gctPhi());
    UCTTowerIndex t = g.getUCTTowerIndexFromL1CaloRegion(r, i.raw());
    uint32_t absCaloEta = std::abs(t.first);
    uint32_t absCaloPhi = std::abs(t.second);
    bool negativeEta = false;
    if (t.first < 0)
      negativeEta = true;
    uint32_t crate = g.getCrate(t.first, t.second);
    uint32_t card = g.getCard(t.first, t.second);
    uint32_t region = g.getRegion(absCaloEta, absCaloPhi);
    UCTRegion* test = new UCTRegion(crate, card, negativeEta, region, fwVersion);
    test->setRegionSummary(i.raw());
    inputRegions.push_back(test);
    //This *should* fill the tensor in the proper order to be fed to the anomaly model
    //We take 4 off of the GCT eta/iEta.
    //iEta taken from this ranges from 4-17, (I assume reserving lower and higher for forward regions)
    //So our first index, index 0, is technically iEta=4, and so-on.
    //CICADA v1 reads this as a flat vector
    modelInput[14 * i.gctPhi() + (i.gctEta() - 4)] = i.et();
  }
  //Extract model output
  //Would be good to be able to configure the precision of the result
  ap_fixed<11, 5> modelResult[1];
  model->prepare_input(modelInput);
  model->predict();
  model->read_result(modelResult);

  *anomalyScore = modelResult[0].to_float();

  summaryCard.setRegionData(inputRegions);

  if (!summaryCard.process()) {
    edm::LogError("L1TCaloSummary") << "UCT: Failed to process summary card" << std::endl;
    exit(1);
  }

  double pt = 0;
  double eta = -999.;
  double phi = -999.;
  double mass = 0;

  std::list<UCTObject*> boostedJetObjs = summaryCard.getBoostedJetObjs();
  for (std::list<UCTObject*>::const_iterator i = boostedJetObjs.begin(); i != boostedJetObjs.end(); i++) {
    const UCTObject* object = *i;
    pt = ((double)object->et()) * caloScaleFactor * boostedJetPtFactor;
    eta = g.getUCTTowerEta(object->iEta());
    phi = g.getUCTTowerPhi(object->iPhi());
    bitset<3> activeRegionEtaPattern = 0;
    for (uint32_t iEta = 0; iEta < 3; iEta++) {
      bool activeStrip = false;
      for (uint32_t iPhi = 0; iPhi < 3; iPhi++) {
        if (object->boostedJetRegionET()[3 * iEta + iPhi] > 30 &&
            object->boostedJetRegionET()[3 * iEta + iPhi] > object->et() * 0.0625)
          activeStrip = true;
      }
      if (activeStrip)
        activeRegionEtaPattern |= (0x1 << iEta);
    }
    bitset<3> activeRegionPhiPattern = 0;
    for (uint32_t iPhi = 0; iPhi < 3; iPhi++) {
      bool activeStrip = false;
      for (uint32_t iEta = 0; iEta < 3; iEta++) {
        if (object->boostedJetRegionET()[3 * iEta + iPhi] > 30 &&
            object->boostedJetRegionET()[3 * iEta + iPhi] > object->et() * 0.0625)
          activeStrip = true;
      }
      if (activeStrip)
        activeRegionPhiPattern |= (0x1 << iPhi);
    }
    string regionEta = activeRegionEtaPattern.to_string<char, std::string::traits_type, std::string::allocator_type>();
    string regionPhi = activeRegionPhiPattern.to_string<char, std::string::traits_type, std::string::allocator_type>();

    bool centralHighest = object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[0] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[1] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[2] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[3] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[5] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[6] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[7] &&
                          object->boostedJetRegionET()[4] >= object->boostedJetRegionET()[8];

    if (abs(eta) < 2.5 && ((regionEta == "101" && (regionPhi == "110" || regionPhi == "101" || regionPhi == "010")) ||
                           ((regionEta == "110" || regionEta == "101" || regionEta == "010") && regionPhi == "101") ||
                           (regionEta == "111" && (regionPhi == "110" || regionPhi == "010")) ||
                           ((regionEta == "110" || regionEta == "010") && regionPhi == "111") ||
                           ((regionEta == "010" || regionPhi == "010" || regionEta == "110" || regionPhi == "110" ||
                             regionEta == "011" || regionPhi == "011") &&
                            centralHighest)))
      bJetCands->push_back(L1JetParticle(math::PtEtaPhiMLorentzVector(pt, eta, phi, mass), L1JetParticle::kCentral));
  }

  iEvent.put(std::move(bJetCands), "Boosted");
  //Write out anomaly score
  iEvent.put(std::move(anomalyScore), "anomalyScore");
}

Member Data Documentation

boostedJetPtFactor

double L1TCaloSummary::boostedJetPtFactor

private

Definition at line 101 of file L1TCaloSummary.cc.

Referenced by produce().

caloScaleFactor

double L1TCaloSummary::caloScaleFactor

private

Definition at line 94 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary(), and produce().

eGammaIsolationFactor

double L1TCaloSummary::eGammaIsolationFactor

private

Definition at line 100 of file L1TCaloSummary.cc.

Referenced by produce().

eGammaSeed

uint32_t L1TCaloSummary::eGammaSeed

private

Definition at line 99 of file L1TCaloSummary.cc.

Referenced by produce().

fwVersion

int L1TCaloSummary::fwVersion

private

Definition at line 104 of file L1TCaloSummary.cc.

Referenced by produce().

jetSeed

uint32_t L1TCaloSummary::jetSeed

private

Definition at line 96 of file L1TCaloSummary.cc.

Referenced by produce().

layer1

UCTLayer1* L1TCaloSummary::layer1

private

Definition at line 108 of file L1TCaloSummary.cc.

loader

hls4mlEmulator::ModelLoader L1TCaloSummary::loader

private

Definition at line 110 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

model

std::shared_ptr<hls4mlEmulator::Model> L1TCaloSummary::model

private

Definition at line 111 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary(), produce(), and querying.connection::search_everything().

nPumBins

uint32_t L1TCaloSummary::nPumBins

private

Definition at line 90 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

pumLUT

std::vector<std::vector<std::vector<uint32_t> > > L1TCaloSummary::pumLUT

private

Definition at line 92 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary(), and produce().

regionToken

edm::EDGetTokenT<L1CaloRegionCollection> L1TCaloSummary::regionToken

private

Definition at line 106 of file L1TCaloSummary.cc.

Referenced by produce().

tauIsolationFactor

float L1TCaloSummary::tauIsolationFactor

private

Definition at line 98 of file L1TCaloSummary.cc.

Referenced by produce().

tauSeed

uint32_t L1TCaloSummary::tauSeed

private

Definition at line 97 of file L1TCaloSummary.cc.

Referenced by produce().

verbose

bool L1TCaloSummary::verbose

private

Definition at line 103 of file L1TCaloSummary.cc.

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