MuIsolatorResultProducer< BT > Class Template Reference

BT == base type. More...

#include <MuIsolatorResultProducer.h>

Inheritance diagram for MuIsolatorResultProducer< BT >:

Classes
struct	DepositConf
struct	VetoCuts

Public Member Functions
	MuIsolatorResultProducer (const edm::ParameterSet &)
	constructor with config More...
void	produce (edm::Event &, const edm::EventSetup &) override
	~MuIsolatorResultProducer () override
	destructor More...
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 muisorhelper::CandMap< BT >	CandMap
typedef muisorhelper::DepositContainer	DepositContainer
typedef muisorhelper::Isolator	Isolator
typedef muisorhelper::Result	Result
typedef muisorhelper::Results	Results
typedef muisorhelper::ResultType	ResultType

Private Member Functions
void	callWhatProduces ()
	declare what's going to be produced More...
unsigned int	initAssociation (edm::Event &event, CandMap &candMapT) const
void	initVetos (reco::TrackBase::Point const &theBeam, std::vector< reco::IsoDeposit::Vetos *> &vetos, CandMap &candMap) const
void	writeOut (edm::Event &event, const CandMap &candMap, const Results &results) const
	choose which result type to write here More...
template<typename RT >
void	writeOutImpl (edm::Event &event, const CandMap &candMapT, const Results &results) const
	actually do the writing here More...

Private Attributes
std::string	theBeamlineOption
	beam spot More...
edm::InputTag	theBeamSpotLabel
edm::ParameterSet	theConfig
std::vector< DepositConf >	theDepositConfs
Isolator *	theIsolator
	the isolator More...
bool	theRemoveOtherVetos
	choose which muon vetos should be removed from all deposits More...
ResultType	theResultType
VetoCuts	theVetoCuts

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<typename BT = reco::Candidate>
class MuIsolatorResultProducer< BT >

BT == base type.

Definition at line 53 of file MuIsolatorResultProducer.h.

Member Typedef Documentation

CandMap

template<typename BT = reco::Candidate>

typedef muisorhelper::CandMap<BT> MuIsolatorResultProducer< BT >::CandMap

private

Definition at line 68 of file MuIsolatorResultProducer.h.

DepositContainer

template<typename BT = reco::Candidate>

typedef muisorhelper::DepositContainer MuIsolatorResultProducer< BT >::DepositContainer

private

Definition at line 66 of file MuIsolatorResultProducer.h.

Isolator

template<typename BT = reco::Candidate>

typedef muisorhelper::Isolator MuIsolatorResultProducer< BT >::Isolator

private

Definition at line 62 of file MuIsolatorResultProducer.h.

Result

template<typename BT = reco::Candidate>

typedef muisorhelper::Result MuIsolatorResultProducer< BT >::Result

private

Definition at line 63 of file MuIsolatorResultProducer.h.

Results

template<typename BT = reco::Candidate>

typedef muisorhelper::Results MuIsolatorResultProducer< BT >::Results

private

Definition at line 65 of file MuIsolatorResultProducer.h.

ResultType

template<typename BT = reco::Candidate>

typedef muisorhelper::ResultType MuIsolatorResultProducer< BT >::ResultType

private

Definition at line 64 of file MuIsolatorResultProducer.h.

Constructor & Destructor Documentation

MuIsolatorResultProducer()

template<typename BT >

MuIsolatorResultProducer< BT >::MuIsolatorResultProducer

(

const edm::ParameterSet &

par

)

constructor with config

read input config for deposit types and weights and thresholds to apply to them

will switch to a factory at some point

FIXME

FIXME: need to get the file name here

"other vetoes" is limited to the same collection now for non-trivial choice an external map with pre-made selection flags can be a better choice

Definition at line 186 of file MuIsolatorResultProducer.h.

References MuIsolatorResultProducer< BT >::callWhatProduces(), HLT_2023v12_cff::coneSize, edm::ParameterSet::getParameter(), LogDebug, MuIsolatorResultProducer< BT >::VetoCuts::muAbsEtaMax, MuIsolatorResultProducer< BT >::VetoCuts::muAbsZMax, MuIsolatorResultProducer< BT >::VetoCuts::muD0Max, MuIsolatorResultProducer< BT >::VetoCuts::muPtMin, muonisolation::MuIsoBaseIsolator::resultType(), MuIsolatorResultProducer< BT >::VetoCuts::selectAll, AlCaHLTBitMon_QueryRunRegistry::string, MuIsolatorResultProducer< BT >::DepositConf::tag, MuIsolatorResultProducer< BT >::theBeamlineOption, MuIsolatorResultProducer< BT >::theBeamSpotLabel, MuIsolatorResultProducer< BT >::theDepositConfs, MuIsolatorResultProducer< BT >::theIsolator, MuIsolatorResultProducer< BT >::theRemoveOtherVetos, MuIsolatorResultProducer< BT >::theResultType, MuIsolatorResultProducer< BT >::theVetoCuts, MuIsolatorResultProducer< BT >::DepositConf::threshold, and MuIsolatorResultProducer< BT >::DepositConf::weight.

    : theConfig(par), theRemoveOtherVetos(par.getParameter<bool>("RemoveOtherVetos")), theIsolator(nullptr) {
  LogDebug("RecoMuon|MuonIsolation") << " MuIsolatorResultProducer CTOR";

  std::vector<edm::ParameterSet> depositInputs = par.getParameter<std::vector<edm::ParameterSet>>("InputMuIsoDeposits");

  std::vector<double> dWeights(depositInputs.size());
  std::vector<double> dThresholds(depositInputs.size());

  for (unsigned int iDep = 0; iDep < depositInputs.size(); ++iDep) {
    DepositConf dConf;
    dConf.tag = depositInputs[iDep].getParameter<edm::InputTag>("DepositTag");
    dConf.weight = depositInputs[iDep].getParameter<double>("DepositWeight");
    dConf.threshold = depositInputs[iDep].getParameter<double>("DepositThreshold");

    dWeights[iDep] = dConf.weight;
    dThresholds[iDep] = dConf.threshold;

    theDepositConfs.push_back(dConf);
  }

  edm::ParameterSet isoPset = par.getParameter<edm::ParameterSet>("IsolatorPSet");
  std::string isolatorType = isoPset.getParameter<std::string>("ComponentName");
  if (isolatorType == "IsolatorByDeposit") {
    std::string coneSizeType = isoPset.getParameter<std::string>("ConeSizeType");
    if (coneSizeType == "FixedConeSize") {
      float coneSize(isoPset.getParameter<double>("coneSize"));

      theIsolator = new muonisolation::IsolatorByDeposit(coneSize, dWeights, dThresholds);

      //      theIsolator = new IsolatorByDeposit(isoPset);
    } else if (coneSizeType == "CutsConeSize") {
      //       Cuts cuts(isoPset.getParameter<edm::ParameterSet>("CutsPSet"));

      //       theIsolator = new IsolatorByDeposit(coneSize, dWeights, dThresholds);
    }
  } else if (isolatorType == "IsolatorByNominalEfficiency") {
    theIsolator =
        new muonisolation::IsolatorByNominalEfficiency("noname", std::vector<std::string>(1, "8:0.97"), dWeights);
  } else if (isolatorType == "IsolatorByDepositCount") {
    std::string coneSizeType = isoPset.getParameter<std::string>("ConeSizeType");
    if (coneSizeType == "FixedConeSize") {
      float coneSize(isoPset.getParameter<double>("coneSize"));

      theIsolator = new muonisolation::IsolatorByDepositCount(coneSize, dThresholds);

      //      theIsolator = new IsolatorByDeposit(isoPset);
    } else if (coneSizeType == "CutsConeSize") {
      //       Cuts cuts(isoPset.getParameter<edm::ParameterSet>("CutsPSet"));

      //       theIsolator = new IsolatorByDeposit(coneSize, dWeights, dThresholds);
    }
  }

  if (theIsolator == nullptr) {
    edm::LogError("MuonIsolationProducers") << "Failed to initialize an isolator";
  }
  theResultType = theIsolator->resultType();

  callWhatProduces();

  if (theRemoveOtherVetos) {
    edm::ParameterSet vetoPSet = par.getParameter<edm::ParameterSet>("VetoPSet");
    theVetoCuts.selectAll = vetoPSet.getParameter<bool>("SelectAll");

    if (!theVetoCuts.selectAll) {
      theVetoCuts.muAbsEtaMax = vetoPSet.getParameter<double>("MuAbsEtaMax");
      theVetoCuts.muPtMin = vetoPSet.getParameter<double>("MuPtMin");
      theVetoCuts.muAbsZMax = vetoPSet.getParameter<double>("MuAbsZMax");
      theVetoCuts.muD0Max = vetoPSet.getParameter<double>("MuD0Max");
      theBeamlineOption = par.getParameter<std::string>("BeamlineOption");
      theBeamSpotLabel = par.getParameter<edm::InputTag>("BeamSpotLabel");
    }
  }
}

~MuIsolatorResultProducer()

template<typename BT >

MuIsolatorResultProducer< BT >::~MuIsolatorResultProducer ( )

override

destructor

Definition at line 271 of file MuIsolatorResultProducer.h.

References LogDebug.

                                                        {
  if (theIsolator)
    delete theIsolator;
  LogDebug("RecoMuon|MuIsolatorResultProducer") << " MuIsolatorResultProducer DTOR";
}

Member Function Documentation

callWhatProduces()

template<typename BT >

void MuIsolatorResultProducer< BT >::callWhatProduces ( )

inlineprivate

declare what's going to be produced

Definition at line 153 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

                                                           {
  if (theResultType == Isolator::ISOL_FLOAT_TYPE)
    produces<edm::ValueMap<float>>();
  if (theResultType == Isolator::ISOL_INT_TYPE)
    produces<edm::ValueMap<int>>();
  if (theResultType == Isolator::ISOL_BOOL_TYPE)
    produces<edm::ValueMap<bool>>();
}

initAssociation()

template<typename BT >

unsigned int MuIsolatorResultProducer< BT >::initAssociation ( edm::Event &  event, CandMap &  candMapT  ) const

private

WARNING: the input ValueMaps are better be for a single key product ID no effort is done (FIXME) for more complex cases
init {muon, {[deposit,veto]_type}} container

Definition at line 345 of file MuIsolatorResultProducer.h.

References edm::ValueMap< T >::const_iterator::begin(), edm::ValueMap< T >::begin(), edm::ValueMap< T >::empty(), edm::ValueMap< T >::const_iterator::end(), muisorhelper::CandMap< BT >::get(), edm::ValueMap< T >::const_iterator::id(), LogDebug, metname, muisorhelper::CandMap< BT >::setHandle(), edm::ValueMap< T >::size(), AlCaHLTBitMon_QueryRunRegistry::string, and makeGlobalPositionRcd_cfg::tag.

                                                                                                   {
  std::string metname = "RecoMuon|MuonIsolationProducers";

  typedef reco::IsoDepositMap::container CT;

  for (unsigned int iMap = 0; iMap < theDepositConfs.size(); ++iMap) {
    edm::Handle<reco::IsoDepositMap> depH;
    event.getByLabel(theDepositConfs[iMap].tag, depH);
    LogDebug(metname) << "Got Deposits of size " << depH->size();
    if (depH->empty())
      continue;

    typename CandMap::handle_type keyH;
    event.get(depH->begin().id(), keyH);
    candMapT.setHandle(keyH);
    typename CT::const_iterator depHCI = depH->begin().begin();
    typename CT::const_iterator depEnd = depH->begin().end();
    unsigned int keyI = 0;
    for (; depHCI != depEnd; ++depHCI, ++keyI) {
      typename CandMap::key_type muPtr(keyH->refAt(keyI));
      if (iMap == 0)
        candMapT.get().push_back(typename CandMap::pair_type(muPtr, DepositContainer(theDepositConfs.size())));
      typename CandMap::iterator muI = candMapT.get().begin();
      for (; muI != candMapT.get().end(); ++muI) {
        if (muI->first == muPtr)
          break;
      }
      if (muI->first != muPtr) {
        edm::LogError("MuonIsolationProducers") << "Failed to align muon map";
      }
      muI->second[iMap].dep = &*depHCI;
    }
  }

  LogDebug(metname) << "Picked and aligned nDeps = " << candMapT.get().size();
  return candMapT.get().size();
}

initVetos()

template<typename BT >

void MuIsolatorResultProducer< BT >::initVetos ( reco::TrackBase::Point const &  theBeam, std::vector< reco::IsoDeposit::Vetos *> &  vetos, CandMap &  candMap  ) const

private

Definition at line 387 of file MuIsolatorResultProducer.h.

References d0, muisorhelper::CandMap< BT >::get(), LogDebug, metname, amptDefaultParameters_cff::mu, AlCaHLTBitMon_QueryRunRegistry::string, and boostedElectronIsolation_cff::vetos.

                                                                      {
  std::string metname = "RecoMuon|MuonIsolationProducers";

  if (theRemoveOtherVetos) {
    LogDebug(metname) << "Start checking for vetos based on input/expected vetos.size of " << vetos.size()
                      << " passed at " << &vetos << " and an input map.size of " << candMapT.get().size();

    unsigned int muI = 0;
    for (; muI < candMapT.get().size(); ++muI) {
      typename CandMap::key_type mu = candMapT.get()[muI].first;
      double d0 = ((mu->vx() - theBeam.x()) * mu->py() - (mu->vy() - theBeam.y()) * mu->px()) / mu->pt();
      LogDebug(metname) << "Muon at index " << muI;
      if (theVetoCuts.selectAll || (fabs(mu->eta()) < theVetoCuts.muAbsEtaMax && mu->pt() > theVetoCuts.muPtMin &&
                                    fabs(mu->vz()) < theVetoCuts.muAbsZMax && fabs(d0) < theVetoCuts.muD0Max)) {
        LogDebug(metname) << "muon passes the cuts";
        for (unsigned int iDep = 0; iDep < candMapT.get()[muI].second.size(); ++iDep) {
          if (vetos[iDep] == nullptr)
            vetos[iDep] = new reco::IsoDeposit::Vetos();

          vetos[iDep]->push_back(candMapT.get()[muI].second[iDep].dep->veto());
        }
      }
    }

    LogDebug(metname) << "Assigning vetos";
    muI = 0;
    for (; muI < candMapT.get().size(); ++muI) {
      for (unsigned int iDep = 0; iDep < candMapT.get()[muI].second.size(); ++iDep) {
        candMapT.get()[muI].second[iDep].vetos = vetos[iDep];
      }
    }
    LogDebug(metname) << "Done with vetos";
  }
}

produce()

template<typename BT >

void MuIsolatorResultProducer< BT >::produce ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

override

do it only if needed

"smart" container used to repackage deposits_type_candIndex into deposits_candIndex_type have to have it for veto removal (could do away without it otherwise) IMPORTANT: ALL THE REFERENCING BUSINESS IS DONE THROUGH POINTERS Access to the mapped values as reference type HAS TO BE AVAILABLE



isolator results will be here

extra vetos will be filled here

call the isolator result, passing {[deposit,vetos]_type} set and the candidate

do cleanup

Definition at line 278 of file MuIsolatorResultProducer.h.

References pwdgSkimBPark_cfi::beamSpot, muisorhelper::CandMap< BT >::get(), edm::HandleBase::isValid(), LogDebug, LogTrace, metname, reco::BeamSpot::position(), mysort::results, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                           {
  std::string metname = "RecoMuon|MuonIsolationProducers";
  LogDebug(metname) << " Muon Deposit producing..."
                    << " BEGINING OF EVENT "
                    << "================================";

  reco::TrackBase::Point theBeam = reco::TrackBase::Point(0, 0, 0);

  if (theRemoveOtherVetos && !theVetoCuts.selectAll) {
    if (theBeamlineOption == "BeamSpotFromEvent") {
      //pick beamSpot
      reco::BeamSpot beamSpot;
      edm::Handle<reco::BeamSpot> beamSpotH;

      event.getByLabel(theBeamSpotLabel, beamSpotH);

      if (beamSpotH.isValid()) {
        theBeam = beamSpotH->position();
        LogTrace(metname) << "Extracted beam point at " << theBeam << std::endl;
      }
    }
  }

  CandMap candMapT;

  unsigned int colSize = initAssociation(event, candMapT);

  Results results(colSize);

  std::vector<reco::IsoDeposit::Vetos*> vetoDeps(theDepositConfs.size(), nullptr);

  if (colSize != 0) {
    if (theRemoveOtherVetos) {
      initVetos(theBeam, vetoDeps, candMapT);
    }

    for (unsigned int muI = 0; muI < colSize; ++muI) {
      results[muI] = theIsolator->result(candMapT.get()[muI].second, *(candMapT.get()[muI].first));

      if (results[muI].typeF() != theIsolator->resultType()) {
        edm::LogError("MuonIsolationProducers") << "Failed to get result from the isolator";
      }
    }
  }

  LogDebug(metname) << "Ready to write out results of size " << results.size();
  writeOut(event, candMapT, results);

  for (unsigned int iDep = 0; iDep < vetoDeps.size(); ++iDep) {
    if (vetoDeps[iDep]) {
      delete vetoDeps[iDep];
      vetoDeps[iDep] = nullptr;
    }
  }
}

writeOut()

template<typename BT >

void MuIsolatorResultProducer< BT >::writeOut ( edm::Event &  event, const CandMap &  candMap, const Results &  results  ) const

inlineprivate

choose which result type to write here

Definition at line 137 of file MuIsolatorResultProducer.h.

References edmPickEvents::event, LogDebug, metname, mysort::results, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
  std::string metname = "RecoMuon|MuonIsolationProducers";
  LogDebug(metname) << "Before calling writeOutMap  with result type " << theIsolator->resultType();

  if (theResultType == Isolator::ISOL_INT_TYPE)
    writeOutImpl<int>(event, candMapT, results);
  if (theResultType == Isolator::ISOL_FLOAT_TYPE)
    writeOutImpl<float>(event, candMapT, results);
  if (theResultType == Isolator::ISOL_BOOL_TYPE)
    writeOutImpl<bool>(event, candMapT, results);
}

writeOutImpl()

template<typename BT >

template<typename RT >

void MuIsolatorResultProducer< BT >::writeOutImpl ( edm::Event &  event, const CandMap &  candMapT, const Results &  results  ) const

inlineprivate

actually do the writing here

make an output vec of what's to be written with a concrete type

fill/insert of non-empty values only

Definition at line 116 of file MuIsolatorResultProducer.h.

References trigObjTnPSource_cfi::filler, muisorhelper::CandMap< BT >::get(), muisorhelper::CandMap< BT >::handle(), mps_fire::i, eostools::move(), and mysort::results.

                                                                                     {
  std::vector<RT> resV(results.size());
  for (unsigned int i = 0; i < resV.size(); ++i)
    resV[i] = results[i].val<RT>();
  auto outMap = std::make_unique<edm::ValueMap<RT>>();
  typename edm::ValueMap<RT>::Filler filler(*outMap);

  if (!candMapT.get().empty()) {
    filler.insert(candMapT.handle(), resV.begin(), resV.end());
    filler.fill();
  }

  event.put(std::move(outMap));
}

Member Data Documentation

theBeamlineOption

template<typename BT = reco::Candidate>

std::string MuIsolatorResultProducer< BT >::theBeamlineOption

private

beam spot

Definition at line 109 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

theBeamSpotLabel

template<typename BT = reco::Candidate>

edm::InputTag MuIsolatorResultProducer< BT >::theBeamSpotLabel

private

Definition at line 110 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

theConfig

template<typename BT = reco::Candidate>

edm::ParameterSet MuIsolatorResultProducer< BT >::theConfig

private

Definition at line 97 of file MuIsolatorResultProducer.h.

theDepositConfs

template<typename BT = reco::Candidate>

std::vector<DepositConf> MuIsolatorResultProducer< BT >::theDepositConfs

private

Definition at line 98 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

theIsolator

template<typename BT = reco::Candidate>

Isolator* MuIsolatorResultProducer< BT >::theIsolator

private

the isolator

Definition at line 105 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

theRemoveOtherVetos

template<typename BT = reco::Candidate>

bool MuIsolatorResultProducer< BT >::theRemoveOtherVetos

private

choose which muon vetos should be removed from all deposits

Definition at line 101 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

theResultType

template<typename BT = reco::Candidate>

ResultType MuIsolatorResultProducer< BT >::theResultType

private

Definition at line 106 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().

theVetoCuts

template<typename BT = reco::Candidate>

VetoCuts MuIsolatorResultProducer< BT >::theVetoCuts

private

Definition at line 102 of file MuIsolatorResultProducer.h.

Referenced by MuIsolatorResultProducer< BT >::MuIsolatorResultProducer().