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

Constructor & Destructor Documentation

L1TCaloSummary()

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

explicit

Definition at line 121 of file L1TCaloSummary.cc.

    : 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);
}

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

~L1TCaloSummary()

L1TCaloSummary::~L1TCaloSummary ( )

override

Definition at line 156 of file L1TCaloSummary.cc.

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

References summaryCard.

Member Function Documentation

beginRun()

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

overrideprivate

Definition at line 259 of file L1TCaloSummary.cc.

{}

fillDescriptions()

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

static

Definition at line 286 of file L1TCaloSummary.cc.

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

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

print()

void L1TCaloSummary::print ( void  )

private

Definition at line 249 of file L1TCaloSummary.cc.

{}

produce()

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

overrideprivate

Definition at line 166 of file L1TCaloSummary.cc.

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

References funct::abs(), boostedJetPtFactor, caloScaleFactor, PVValHelper::eta, beamvalidation::exit(), fwVersion, g, mps_fire::i, L1TowerCalibrationProducer_cfi::iEta, iEvent, EgHLTOffHistBins_cfi::mass, eostools::move(), resolutioncreator_cfi::object, phi, DiDispStaMuonMonitor_cfi::pt, alignCSCRings::r, HLT_FULL_cff::region, regionToken, summaryCard, and submitPVValidationJobs::t.

Member Data Documentation

boostedJetPtFactor

double L1TCaloSummary::boostedJetPtFactor

private

Definition at line 99 of file L1TCaloSummary.cc.

Referenced by produce().

caloScaleFactor

double L1TCaloSummary::caloScaleFactor

private

Definition at line 92 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary(), and produce().

eGammaIsolationFactor

double L1TCaloSummary::eGammaIsolationFactor

private

Definition at line 98 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

eGammaSeed

uint32_t L1TCaloSummary::eGammaSeed

private

Definition at line 97 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

fwVersion

int L1TCaloSummary::fwVersion

private

Definition at line 102 of file L1TCaloSummary.cc.

Referenced by produce().

jetSeed

uint32_t L1TCaloSummary::jetSeed

private

Definition at line 94 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

layer1

UCTLayer1* L1TCaloSummary::layer1

private

Definition at line 106 of file L1TCaloSummary.cc.

nPumBins

uint32_t L1TCaloSummary::nPumBins

private

Definition at line 88 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

pumLUT

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

private

Definition at line 90 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

regionToken

edm::EDGetTokenT<L1CaloRegionCollection> L1TCaloSummary::regionToken

private

Definition at line 104 of file L1TCaloSummary.cc.

Referenced by produce().

summaryCard

UCTSummaryCard* L1TCaloSummary::summaryCard

private

Definition at line 107 of file L1TCaloSummary.cc.

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

tauIsolationFactor

float L1TCaloSummary::tauIsolationFactor

private

Definition at line 96 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

tauSeed

uint32_t L1TCaloSummary::tauSeed

private

Definition at line 95 of file L1TCaloSummary.cc.

Referenced by L1TCaloSummary().

verbose

bool L1TCaloSummary::verbose

private

Definition at line 101 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().