MuIsoDepositProducer Class Reference

#include <MuIsoDepositProducer.h>

Inheritance diagram for MuIsoDepositProducer:

Public Member Functions
	MuIsoDepositProducer (const edm::ParameterSet &)
	constructor More...
void	produce (edm::Event &, const edm::EventSetup &) override
	data making method More...
	~MuIsoDepositProducer () 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 Attributes
std::vector< std::string >	theDepositNames
bool	theExtractForCandidate
std::unique_ptr< reco::isodeposit::IsoDepositExtractor >	theExtractor
std::string	theInputType
bool	theMultipleDepositsFlag
edm::EDGetToken	theMuonCollectionTag
std::string	theMuonTrackRefType

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

Definition at line 17 of file MuIsoDepositProducer.h.

Constructor & Destructor Documentation

MuIsoDepositProducer()

MuIsoDepositProducer::MuIsoDepositProducer

(

const edm::ParameterSet &

par

)

constructor

constructor with config

Definition at line 31 of file MuIsoDepositProducer.cc.

References SiStripOfflineCRack_cfg::alias, relativeConstraints::empty, Exception, get, edm::ParameterSet::getParameter(), mps_fire::i, LogDebug, metname, AlCaHLTBitMon_QueryRunRegistry::string, theDepositNames, theExtractForCandidate, theExtractor, theInputType, theMultipleDepositsFlag, theMuonCollectionTag, and theMuonTrackRefType.

    : theDepositNames(std::vector<std::string>(1, std::string())) {
  static const std::string metname = "RecoMuon|MuonIsolationProducers|MuIsoDepositProducer";
  LogDebug(metname) << " MuIsoDepositProducer CTOR";

  edm::ParameterSet ioPSet = par.getParameter<edm::ParameterSet>("IOPSet");

  theInputType = ioPSet.getParameter<std::string>("InputType");
  theExtractForCandidate = ioPSet.getParameter<bool>("ExtractForCandidate");
  theMuonTrackRefType = ioPSet.getParameter<std::string>("MuonTrackRefType");

  bool readFromRecoTrack = theInputType == "TrackCollection";
  bool readFromRecoMuon = theInputType == "MuonCollection";
  bool readFromCandidateView = theInputType == "CandidateView";
  if (readFromRecoTrack) {
    theMuonCollectionTag = consumes<View<Track>>(ioPSet.getParameter<edm::InputTag>("inputMuonCollection"));
  } else if (readFromRecoMuon) {
    theMuonCollectionTag = consumes<View<RecoCandidate>>(ioPSet.getParameter<edm::InputTag>("inputMuonCollection"));
  } else if (readFromCandidateView) {
    theMuonCollectionTag = consumes<View<Candidate>>(ioPSet.getParameter<edm::InputTag>("inputMuonCollection"));
  } else {
    throw cms::Exception("Configuration") << "Inconsistent configuration or failure to read Candidate-muon view";
  }

  theMultipleDepositsFlag = ioPSet.getParameter<bool>("MultipleDepositsFlag");

  if (theMultipleDepositsFlag) {
    theDepositNames = par.getParameter<edm::ParameterSet>("ExtractorPSet")
                          .getParameter<std::vector<std::string>>("DepositInstanceLabels");
  }

  for (unsigned int i = 0; i < theDepositNames.size(); ++i) {
    std::string alias = par.getParameter<std::string>("@module_label");
    if (!theDepositNames[i].empty())
      alias += "_" + theDepositNames[i];
    produces<reco::IsoDepositMap>(theDepositNames[i]).setBranchAlias(alias);
  }

  edm::ParameterSet extractorPSet = par.getParameter<edm::ParameterSet>("ExtractorPSet");
  std::string extractorName = extractorPSet.getParameter<std::string>("ComponentName");
  theExtractor = IsoDepositExtractorFactory::get()->create(extractorName, extractorPSet, consumesCollector());
  LogDebug(metname) << " Load extractor..." << extractorName;
}

~MuIsoDepositProducer()

MuIsoDepositProducer::~MuIsoDepositProducer ( )

override

destructor

Definition at line 76 of file MuIsoDepositProducer.cc.

References LogDebug.

                                            {
  LogDebug("RecoMuon/MuIsoDepositProducer") << " MuIsoDepositProducer DTOR";
}

Member Function Documentation

produce()

void MuIsoDepositProducer::produce ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

override

data making method

build deposits

read them as RecoCandidates: need to have track() standAloneMuon() etc in the interface

get rid of this at some point and use the cands

expect nMuons set already

OK, now we know how many deps for how many muons each we will create might linearize this at some point (lazy) do it in case some muons are there only






end for (nMuons)

now fill in selectively

some debugging stuff

fill the maps here

now figure out the source handle (see getByToken above)

now actually fill

end if (nMuons>0)

Definition at line 81 of file MuIsoDepositProducer.cc.

References HLT_2022v15_cff::cands, symbols::deps, mps_fire::end, options_cfi::eventSetup, trigObjTnPSource_cfi::filler, mps_fire::i, edm::EDGetToken::index(), LogDebug, LogTrace, metname, eostools::move(), PDWG_BPHSkim_cff::muons, custom_jme_cff::nMuons, susybsm::HSCParticleType::standAloneMuon, AlCaHLTBitMon_QueryRunRegistry::string, theDepositNames, theExtractForCandidate, theExtractor, theInputType, theMultipleDepositsFlag, theMuonCollectionTag, theMuonTrackRefType, HLT_2022v15_cff::track, and pwdgSkimBPark_cfi::tracks.

                                                                             {
  static const std::string metname = "RecoMuon|MuonIsolationProducers|MuIsoDepositProducer";

  LogDebug(metname) << " Muon Deposit producing..."
                    << " BEGINING OF EVENT "
                    << "================================";

  unsigned int nDeps = theMultipleDepositsFlag ? theDepositNames.size() : 1;

  // Take the muon container
  LogTrace(metname) << " Taking the muons: "
                    << theMuonCollectionTag.index();  //a more friendly print would use "inputMuonCollection"
  Handle<View<Track>> tracks;
  Handle<View<RecoCandidate>> muons;  
  Handle<View<Candidate>> cands;

  unsigned int nMuons = 0;

  bool readFromRecoTrack = theInputType == "TrackCollection";
  bool readFromRecoMuon = theInputType == "MuonCollection";
  bool readFromCandidateView = theInputType == "CandidateView";

  if (readFromRecoMuon) {
    event.getByToken(theMuonCollectionTag, muons);
    nMuons = muons->size();
    LogDebug(metname) << "Got Muons of size " << nMuons;
  }
  if (readFromRecoTrack) {
    event.getByToken(theMuonCollectionTag, tracks);
    nMuons = tracks->size();
    LogDebug(metname) << "Got MuonTracks of size " << nMuons;
  }
  if (readFromCandidateView || theExtractForCandidate) {
    event.getByToken(theMuonCollectionTag, cands);
    unsigned int nCands = cands->size();
    if (readFromRecoMuon && theExtractForCandidate) {
      if (nMuons != nCands)
        edm::LogError(metname) << "Inconsistent configuration or failure to read Candidate-muon view";
    }
    nMuons = nCands;
    LogDebug(metname) << "Got candidate view with size " << nMuons;
  }

  static const unsigned int MAX_DEPS = 10;
  std::unique_ptr<reco::IsoDepositMap> depMaps[MAX_DEPS];

  if (nDeps > 10)
    LogError(metname) << "Unable to handle more than 10 input deposits";
  for (unsigned int i = 0; i < nDeps; ++i) {
    depMaps[i] = std::make_unique<reco::IsoDepositMap>();
  }

  if (nMuons > 0) {
    std::vector<std::vector<IsoDeposit>> deps2D(nDeps, std::vector<IsoDeposit>(nMuons));

    for (unsigned int i = 0; i < nMuons; ++i) {
      TrackBaseRef muRef;
      if (readFromRecoMuon) {
        if (theMuonTrackRefType == "track") {
          muRef = TrackBaseRef((*muons)[i].track());
        } else if (theMuonTrackRefType == "standAloneMuon") {
          muRef = TrackBaseRef((*muons)[i].standAloneMuon());
        } else if (theMuonTrackRefType == "combinedMuon") {
          muRef = TrackBaseRef((*muons)[i].combinedMuon());
        } else if (theMuonTrackRefType == "bestGlbTrkSta") {
          if (!(*muons)[i].combinedMuon().isNull()) {
            muRef = TrackBaseRef((*muons)[i].combinedMuon());
          } else if (!(*muons)[i].track().isNull()) {
            muRef = TrackBaseRef((*muons)[i].track());
          } else {
            muRef = TrackBaseRef((*muons)[i].standAloneMuon());
          }
        } else if (theMuonTrackRefType == "bestTrkSta") {
          if (!(*muons)[i].track().isNull()) {
            muRef = TrackBaseRef((*muons)[i].track());
          } else {
            muRef = TrackBaseRef((*muons)[i].standAloneMuon());
          }
        } else {
          edm::LogWarning(metname) << "Wrong track type is supplied: breaking";
          break;
        }
      } else if (readFromRecoTrack) {
        muRef = TrackBaseRef(tracks, i);
      }

      if (!theMultipleDepositsFlag) {
        if (readFromCandidateView || theExtractForCandidate)
          deps2D[0][i] = theExtractor->deposit(event, eventSetup, (*cands)[i]);
        else
          deps2D[0][i] = theExtractor->deposit(event, eventSetup, muRef);

      } else {
        std::vector<IsoDeposit> deps(nDeps);
        if (readFromCandidateView || theExtractForCandidate)
          deps = theExtractor->deposits(event, eventSetup, (*cands)[i]);
        else
          deps = theExtractor->deposits(event, eventSetup, muRef);
        for (unsigned int iDep = 0; iDep < nDeps; ++iDep) {
          deps2D[iDep][i] = deps[iDep];
        }
      }
    }  

    for (unsigned int iDep = 0; iDep < nDeps; ++iDep) {
      for (unsigned int iMu = 0; iMu < nMuons; ++iMu) {
        LogTrace(metname) << "Contents of " << theDepositNames[iDep] << " for a muon at index " << iMu;
        LogTrace(metname) << deps2D[iDep][iMu].print();
      }

      reco::IsoDepositMap::Filler filler(*depMaps[iDep]);

      if (readFromRecoMuon) {
        filler.insert(muons, deps2D[iDep].begin(), deps2D[iDep].end());
      } else if (readFromRecoTrack) {
        filler.insert(tracks, deps2D[iDep].begin(), deps2D[iDep].end());
      } else if (readFromCandidateView) {
        filler.insert(cands, deps2D[iDep].begin(), deps2D[iDep].end());
      } else {
        edm::LogError(metname) << "Inconsistent configuration: unknown type requested";
      }

      filler.fill();
    }
  }  

  for (unsigned int iMap = 0; iMap < nDeps; ++iMap) {
    LogTrace(metname) << "About to put a deposit named " << theDepositNames[iMap] << " of size "
                      << depMaps[iMap]->size() << " into edm::Event";
    event.put(std::move(depMaps[iMap]), theDepositNames[iMap]);
  }

  LogTrace(metname) << " END OF EVENT "
                    << "================================";
}

Member Data Documentation

theDepositNames

std::vector<std::string> MuIsoDepositProducer::theDepositNames

private

Definition at line 37 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().

theExtractForCandidate

bool MuIsoDepositProducer::theExtractForCandidate

private

Definition at line 33 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().

theExtractor

std::unique_ptr<reco::isodeposit::IsoDepositExtractor> MuIsoDepositProducer::theExtractor

private

Definition at line 39 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().

theInputType

std::string MuIsoDepositProducer::theInputType

private

input type. Choose from:

Definition at line 31 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().

theMultipleDepositsFlag

bool MuIsoDepositProducer::theMultipleDepositsFlag

private

Definition at line 38 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().

theMuonCollectionTag

edm::EDGetToken MuIsoDepositProducer::theMuonCollectionTag

private

Definition at line 36 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().

theMuonTrackRefType

std::string MuIsoDepositProducer::theMuonTrackRefType

private

Definition at line 35 of file MuIsoDepositProducer.h.

Referenced by MuIsoDepositProducer(), and produce().