#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

uint32_t	nPumBins

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

edm::EDGetTokenT < L1CaloRegionCollection >	regionToken

UCTSummaryCard *	summaryCard

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

Constructor & Destructor Documentation

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

explicit

Definition at line 119 of file L1TCaloSummary.cc.

References caloScaleFactor, eGammaIsolationFactor, eGammaSeed, edm::ParameterSet::getParameter(), jetSeed, SiStripPI::max, nPumBins, pumLUT, summaryCard, tauIsolationFactor, and tauSeed.

     : 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"))) {
   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");
   summaryCard = new UCTSummaryCard(&pumLUT, jetSeed, tauSeed, tauIsolationFactor, eGammaSeed, eGammaIsolationFactor);
 }

L1TCaloSummary::~L1TCaloSummary ( )

override

Definition at line 154 of file L1TCaloSummary.cc.

References summaryCard.

                                 {
   if (summaryCard != nullptr)
     delete summaryCard;
 }

Member Function Documentation

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

overrideprivate

Definition at line 257 of file L1TCaloSummary.cc.

257 {}

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

static

Definition at line 284 of file L1TCaloSummary.cc.

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

                                                                                 {
   //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);
 }

void L1TCaloSummary::print ( void )

private

Definition at line 247 of file L1TCaloSummary.cc.

247 {}

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

overrideprivate

Definition at line 164 of file L1TCaloSummary.cc.

References funct::abs(), boostedJetPtFactor, caloScaleFactor, PVValHelper::eta, beamvalidation::exit(), fwVersion, g, edm::Event::getByToken(), mps_fire::i, ResonanceBuilder::mass, eostools::move(), phi, DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), alignCSCRings::r, HLT_FULL_cff::region, regionToken, summaryCard, submitPVValidationJobs::t, and ctpps_dqm_sourceclient-live_cfg::test.

                                                                           {
   using namespace edm;
 
   std::unique_ptr<L1JetParticleCollection> bJetCands(new L1JetParticleCollection);
 
   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.
   summaryCard->clearRegions();
   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);
   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);
   }
   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>();
     if (std::abs(eta) < 2.5 && (regionEta == "010" || regionPhi == "010" || regionEta == "110" || regionPhi == "110" ||
                                 regionEta == "011" || regionPhi == "011"))
       bJetCands->push_back(L1JetParticle(math::PtEtaPhiMLorentzVector(pt, eta, phi, mass), L1JetParticle::kCentral));
   }
 
   iEvent.put(std::move(bJetCands), "Boosted");
 }