pat::L1MuonMatcher Class Reference

Matcher of reconstructed objects to L1 Muons. More...

Inheritance diagram for pat::L1MuonMatcher:

Public Member Functions
	L1MuonMatcher (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~L1MuonMatcher () 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

Private Types
typedef pat::TriggerObjectStandAlone	PATPrimitive
typedef pat::TriggerObjectStandAloneCollection	PATPrimitiveCollection
typedef pat::TriggerObjectStandAloneMatch	PATTriggerAssociation

Private Member Functions
template<typename Hand , typename T >
void	storeExtraInfo (edm::Event &iEvent, const Hand &handle, const std::vector< T > &values, const std::string &label) const
	Store extra information in a ValueMap. More...

Private Attributes
int	firstBX_
	Skim stage2 BX vector. More...
edm::EDGetTokenT< std::vector< l1extra::L1MuonParticle > >	l1Token_
edm::EDGetTokenT< l1t::MuonBxCollection >	l1tToken_
std::string	labelL1_
	Labels to set as filter names in the output. More...
std::string	labelProp_
int	lastBX_
L1MuonMatcherAlgo	matcher_
edm::EDGetTokenT< edm::View< reco::Candidate > >	recoToken_
	Tokens for input collections. More...
bool	useStage2L1_
	Allow to run both on legacy or stage2 (2016) L1 Muon trigger output. More...
bool	writeExtraInfo_
	Write out additional info as ValueMaps. More...

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

Matcher of reconstructed objects to L1 Muons.

Author

Giovanni Petrucciani

Version

Id

HLTL1MuonMatcher.cc,v 1.3 2010/07/12 20:56:11 gpetrucc Exp

Author

Giovanni Petrucciani

Version

Id

L1MuonMatcher.cc,v 1.4 2011/03/31 09:59:33 gpetrucc Exp

Definition at line 32 of file L1MuonMatcher.cc.

Member Typedef Documentation

PATPrimitive

typedef pat::TriggerObjectStandAlone pat::L1MuonMatcher::PATPrimitive

private

Definition at line 40 of file L1MuonMatcher.cc.

PATPrimitiveCollection

typedef pat::TriggerObjectStandAloneCollection pat::L1MuonMatcher::PATPrimitiveCollection

private

Definition at line 41 of file L1MuonMatcher.cc.

PATTriggerAssociation

typedef pat::TriggerObjectStandAloneMatch pat::L1MuonMatcher::PATTriggerAssociation

private

Definition at line 42 of file L1MuonMatcher.cc.

Constructor & Destructor Documentation

L1MuonMatcher()

pat::L1MuonMatcher::L1MuonMatcher ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 74 of file L1MuonMatcher.cc.

References edm::ParameterSet::getParameter(), l1Token_, l1tToken_, useStage2L1_, and writeExtraInfo_.

    : matcher_(iConfig, consumesCollector()),
      recoToken_(consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("src"))),
      labelL1_(iConfig.getParameter<std::string>("setL1Label")),
      labelProp_(iConfig.getParameter<std::string>("setPropLabel")),
      writeExtraInfo_(iConfig.getParameter<bool>("writeExtraInfo")),
      useStage2L1_(iConfig.getParameter<bool>("useStage2L1")),
      firstBX_(iConfig.getParameter<int>("firstBX")),
      lastBX_(iConfig.getParameter<int>("lastBX")) {
  if (useStage2L1_) {
    l1tToken_ = consumes<l1t::MuonBxCollection>(iConfig.getParameter<edm::InputTag>("matched"));
  } else {
    l1Token_ = consumes<std::vector<l1extra::L1MuonParticle> >(iConfig.getParameter<edm::InputTag>("matched"));
  }
  produces<PATPrimitiveCollection>("l1muons");         // l1 in PAT format
  produces<PATPrimitiveCollection>("propagatedReco");  // reco to muon station 2
  produces<PATTriggerAssociation>("propagatedReco");   // asso reco to propagated reco
  produces<PATTriggerAssociation>();                   // asso reco to l1
  if (writeExtraInfo_) {
    produces<edm::ValueMap<float> >("deltaR");
    produces<edm::ValueMap<float> >("deltaPhi");
    produces<edm::ValueMap<int> >("quality");
    produces<edm::ValueMap<int> >("bx");
    if (useStage2L1_) {
      produces<edm::ValueMap<int> >("iPhi");
      produces<edm::ValueMap<int> >("tfIndex");
    }
    produces<edm::ValueMap<int> >("isolated");
    produces<edm::ValueMap<reco::CandidatePtr> >();
    produces<edm::ValueMap<reco::CandidatePtr> >("l1ToReco");
  }
}

~L1MuonMatcher()

pat::L1MuonMatcher::~L1MuonMatcher ( )

inlineoverride

Definition at line 35 of file L1MuonMatcher.cc.

{}

Member Function Documentation

produce()

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

override

Definition at line 107 of file L1MuonMatcher.cc.

References BXVector< T >::begin(), simKBmtfDigis_cfi::bx, L1MuGMTCand::bx(), ALCARECOTkAlJpsiMuMu_cff::charge, reco::LeafCandidate::charge(), filterCSVwithJSON::copy, HLT_2022v12_cff::deltaPhis, BXVector< T >::end(), edm::helper::Filler< Map >::fill(), BXVector< T >::getFirstBX(), BXVector< T >::getLastBX(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, iEvent, edm::helper::Filler< Map >::insert(), L1MuGMTCand::isol(), isSelected(), TrajectoryStateOnSurface::isValid(), L1Scalers_cfi::l1s, match(), SiStripPI::max, hltL1TEGammaFilteredCollectionProducer_cfi::maxBX, SiStripPI::min, hltL1TEGammaFilteredCollectionProducer_cfi::minBX, eostools::move(), amptDefaultParameters_cff::mu, dqmiodumpmetadata::n, reco::LeafCandidate::polarP4(), quality, L1MuGMTCand::quality(), dt_dqm_sourceclient_common_cff::reco, BXVector< T >::size(), and pfDeepBoostedJetPreprocessParams_cfi::upper_bound.

                                                                            {
  using namespace edm;
  using namespace std;

  matcher_.init(iSetup);

  Handle<View<reco::Candidate> > reco;
  Handle<vector<l1extra::L1MuonParticle> > l1s;
  Handle<l1t::MuonBxCollection> l1tBX;

  std::vector<l1t::Muon> l1ts;
  std::vector<size_t> bxIdxs;

  int minBxIdx = 0;
  size_t l1size = 0;

  iEvent.getByToken(recoToken_, reco);

  if (useStage2L1_) {
    iEvent.getByToken(l1tToken_, l1tBX);
    l1size = l1tBX->size();

    int minBX = max(firstBX_, l1tBX->getFirstBX());
    int maxBX = min(lastBX_, l1tBX->getLastBX());

    minBxIdx = l1tBX->begin(minBX) - l1tBX->begin();
    std::copy(l1tBX->begin(minBX), l1tBX->end(maxBX), std::back_inserter(l1ts));

    for (int ibx = l1tBX->getFirstBX(); ibx <= l1tBX->getLastBX(); ++ibx) {
      bxIdxs.push_back(l1tBX->end(ibx) - l1tBX->begin());
    }
  } else {
    iEvent.getByToken(l1Token_, l1s);
    l1size = l1s->size();
  }

  unique_ptr<PATPrimitiveCollection> propOut(new PATPrimitiveCollection());
  unique_ptr<PATPrimitiveCollection> l1Out(new PATPrimitiveCollection());
  std::vector<edm::Ptr<reco::Candidate> > l1rawMatches(reco->size());
  vector<int> isSelected(l1size, -1);
  std::vector<edm::Ptr<reco::Candidate> > whichRecoMatch(l1size);
  vector<int> propMatches(reco->size(), -1);
  vector<int> fullMatches(reco->size(), -1);
  vector<float> deltaRs(reco->size(), 999), deltaPhis(reco->size(), 999);
  vector<int> quality(reco->size(), 0), bx(reco->size(), -999), isolated(reco->size(), -999);
  vector<int> iPhi(reco->size(), 0), tfIndex(reco->size(), -999);
  for (int i = 0, n = reco->size(); i < n; ++i) {
    TrajectoryStateOnSurface propagated;
    const reco::Candidate &mu = (*reco)[i];
    int match = useStage2L1_ ? matcher_.match(mu, l1ts, deltaRs[i], deltaPhis[i], propagated)
                             : matcher_.match(mu, *l1s, deltaRs[i], deltaPhis[i], propagated);
    if (propagated.isValid()) {
      GlobalPoint pos = propagated.globalPosition();
      propMatches[i] = propOut->size();
      propOut->push_back(PATPrimitive(math::PtEtaPhiMLorentzVector(mu.pt(), pos.eta(), pos.phi(), mu.mass())));
      propOut->back().addFilterLabel(labelProp_);
      propOut->back().setCharge(mu.charge());
    }
    if (match != -1) {
      if (useStage2L1_) {
        match += minBxIdx;
      }

      whichRecoMatch[match] = reco->ptrAt(i);

      int charge = 0;
      math::PtEtaPhiMLorentzVector p4;

      if (useStage2L1_) {
        const l1t::Muon &l1t = (*l1tBX)[match];
        charge = l1t.charge();
        p4 = l1t.polarP4();
      } else {
        const l1extra::L1MuonParticle &l1 = (*l1s)[match];
        charge = l1.charge();
        p4 = l1.polarP4();
      }

      if (isSelected[match] == -1) {  // copy to output if needed
        isSelected[match] = l1Out->size();
        l1Out->push_back(PATPrimitive(p4));
        l1Out->back().addFilterLabel(labelL1_);
        l1Out->back().setCharge(charge);
      }

      fullMatches[i] = isSelected[match];  // index in the output collection

      if (useStage2L1_) {
        const l1t::Muon &l1t = (*l1tBX)[match];
        quality[i] = l1t.hwQual();
        bx[i] = l1tBX->getFirstBX() + (std::upper_bound(bxIdxs.begin(), bxIdxs.end(), match) - bxIdxs.begin());
        isolated[i] = l1t.hwIso();
        l1rawMatches[i] = edm::Ptr<reco::Candidate>(l1tBX, size_t(match));
        iPhi[i] = l1t.hwPhi();
        tfIndex[i] = l1t.tfMuonIndex();
      } else {
        const L1MuGMTCand &gmt = (*l1s)[match].gmtMuonCand();
        quality[i] = gmt.quality();
        bx[i] = gmt.bx();
        isolated[i] = gmt.isol();
        l1rawMatches[i] = edm::Ptr<reco::Candidate>(l1s, size_t(match));
      }
    }
  }

  OrphanHandle<PATPrimitiveCollection> l1Done = iEvent.put(std::move(l1Out), "l1muons");
  OrphanHandle<PATPrimitiveCollection> propDone = iEvent.put(std::move(propOut), "propagatedReco");

  unique_ptr<PATTriggerAssociation> propAss(new PATTriggerAssociation(propDone));
  PATTriggerAssociation::Filler propFiller(*propAss);
  propFiller.insert(reco, propMatches.begin(), propMatches.end());
  propFiller.fill();
  iEvent.put(std::move(propAss), "propagatedReco");

  unique_ptr<PATTriggerAssociation> fullAss(new PATTriggerAssociation(l1Done));
  PATTriggerAssociation::Filler fullFiller(*fullAss);
  fullFiller.insert(reco, fullMatches.begin(), fullMatches.end());
  fullFiller.fill();
  iEvent.put(std::move(fullAss));

  if (writeExtraInfo_) {
    storeExtraInfo(iEvent, reco, deltaRs, "deltaR");
    storeExtraInfo(iEvent, reco, deltaPhis, "deltaPhi");
    storeExtraInfo(iEvent, reco, bx, "bx");
    storeExtraInfo(iEvent, reco, isolated, "isolated");
    storeExtraInfo(iEvent, reco, quality, "quality");
    storeExtraInfo(iEvent, reco, l1rawMatches, "");
    if (useStage2L1_) {
      storeExtraInfo(iEvent, l1tBX, whichRecoMatch, "l1ToReco");
      storeExtraInfo(iEvent, reco, tfIndex, "tfIndex");
      storeExtraInfo(iEvent, reco, iPhi, "iPhi");
    } else {
      storeExtraInfo(iEvent, l1s, whichRecoMatch, "l1ToReco");
    }
  }
}

storeExtraInfo()

template<typename Hand , typename T >

void pat::L1MuonMatcher::storeExtraInfo ( edm::Event &  iEvent, const Hand &  handle, const std::vector< T > &  values, const std::string &  label  ) const

private

Store extra information in a ValueMap.

Definition at line 245 of file L1MuonMatcher.cc.

References trigObjTnPSource_cfi::filler, patZpeak::handle, iEvent, label, eostools::move(), and contentValuesCheck::values.

                                                                      {
  using namespace edm;
  using namespace std;
  unique_ptr<ValueMap<T> > valMap(new ValueMap<T>());
  typename edm::ValueMap<T>::Filler filler(*valMap);
  filler.insert(handle, values.begin(), values.end());
  filler.fill();
  iEvent.put(std::move(valMap), label);
}

Member Data Documentation

firstBX_

int pat::L1MuonMatcher::firstBX_

private

Skim stage2 BX vector.

Definition at line 61 of file L1MuonMatcher.cc.

l1Token_

edm::EDGetTokenT<std::vector<l1extra::L1MuonParticle> > pat::L1MuonMatcher::l1Token_

private

Definition at line 48 of file L1MuonMatcher.cc.

Referenced by L1MuonMatcher().

l1tToken_

edm::EDGetTokenT<l1t::MuonBxCollection> pat::L1MuonMatcher::l1tToken_

private

Definition at line 49 of file L1MuonMatcher.cc.

Referenced by L1MuonMatcher().

labelL1_

std::string pat::L1MuonMatcher::labelL1_

private

Labels to set as filter names in the output.

Definition at line 52 of file L1MuonMatcher.cc.

labelProp_

std::string pat::L1MuonMatcher::labelProp_

private

Definition at line 52 of file L1MuonMatcher.cc.

lastBX_

int pat::L1MuonMatcher::lastBX_

private

Definition at line 62 of file L1MuonMatcher.cc.

matcher_

L1MuonMatcherAlgo pat::L1MuonMatcher::matcher_

private

Definition at line 44 of file L1MuonMatcher.cc.

recoToken_

edm::EDGetTokenT<edm::View<reco::Candidate> > pat::L1MuonMatcher::recoToken_

private

Tokens for input collections.

Definition at line 47 of file L1MuonMatcher.cc.

useStage2L1_

bool pat::L1MuonMatcher::useStage2L1_

private

Allow to run both on legacy or stage2 (2016) L1 Muon trigger output.

Definition at line 58 of file L1MuonMatcher.cc.

Referenced by L1MuonMatcher().

writeExtraInfo_

bool pat::L1MuonMatcher::writeExtraInfo_

private

Write out additional info as ValueMaps.

Definition at line 55 of file L1MuonMatcher.cc.

Referenced by L1MuonMatcher().