ZToMuMuSelector Class Reference

Inheritance diagram for ZToMuMuSelector:

Public Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override
	ZToMuMuSelector (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDFilter () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
double	etaCut_
edm::EDGetTokenT< edm::ValueMap< bool > >	isoToken_
edm::EDGetTokenT< edm::ValueMap< bool > >	isoTrackerToken_
double	massZMax_
double	massZMin_
int	minTrackerHits_
edm::EDGetTokenT< reco::TrackCollection >	muonToken_
bool	onlyGlobalMuons_
double	ptCut_
edm::EDGetTokenT< reco::TrackCollection >	trackerToken_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 12 of file ZToMuMuSelector.cc.

Constructor & Destructor Documentation

ZToMuMuSelector()

ZToMuMuSelector::ZToMuMuSelector

(

const edm::ParameterSet &

cfg

)

Definition at line 40 of file ZToMuMuSelector.cc.

    : muonToken_(consumes<TrackCollection>(cfg.getParameter<edm::InputTag>("MuonTag"))),
      isoToken_(consumes<edm::ValueMap<bool> >(cfg.getParameter<edm::InputTag>("IsolationTag"))),
      ptCut_(cfg.getParameter<double>("PtCut")),
      etaCut_(cfg.getParameter<double>("EtaCut")),
      massZMin_(cfg.getParameter<double>("MassZMin")),
      massZMax_(cfg.getParameter<double>("MassZMax")),
 
      onlyGlobalMuons_(cfg.getParameter<bool>("OnlyGlobalMuons")),
      trackerToken_(mayConsume<TrackCollection>(
          cfg.getUntrackedParameter<edm::InputTag>("TrackerTag", edm::InputTag("ctfWithMaterialTracks")))),
      isoTrackerToken_(mayConsume<edm::ValueMap<bool> >(
          cfg.getUntrackedParameter<edm::InputTag>("TrackerIsolationTag", edm::InputTag("zMuMuTrackerIsolations")))),
      minTrackerHits_(cfg.getUntrackedParameter<int>("MinTrackerHits", 7)) {}

Member Function Documentation

filter()

bool ZToMuMuSelector::filter ( edm::Event &  ev, const edm::EventSetup &    )

overridevirtual

Implements edm::EDFilter.

Definition at line 55 of file ZToMuMuSelector.cc.

                                                          {
  Handle<TrackCollection> muonCollection;
  ev.getByToken(muonToken_, muonCollection);
  if (!muonCollection.isValid()) {
    LogTrace("") << ">>> Muon collection does not exist !!!";
    return false;
  }
 
  Handle<edm::ValueMap<bool> > isoMap;
  ev.getByToken(isoToken_, isoMap);
  if (!isoMap.isValid()) {
    LogTrace("") << ">>> ISO Muon collection does not exist !!!";
    return false;
  }
 
  Handle<TrackCollection> trackerCollection;
  Handle<edm::ValueMap<bool> > isoTrackerMap;
  if (!onlyGlobalMuons_) {
    ev.getByToken(trackerToken_, trackerCollection);
    if (!trackerCollection.isValid()) {
      LogTrace("") << ">>> Tracker collection does not exist !!!";
      return false;
    }
 
    ev.getByToken(isoTrackerToken_, isoTrackerMap);
    if (!isoTrackerMap.isValid()) {
      LogTrace("") << ">>> ISO Tracker collection does not exist !!!";
      return false;
    }
  }
 
  unsigned int npairs = 0;
  bool globalCombinationFound = false;
  for (unsigned int i = 0; i < muonCollection->size(); i++) {
    TrackRef mu(muonCollection, i);
    LogTrace("") << "> Processing muon number " << i << "...";
    double pt = mu->pt();
    LogTrace("") << "\t... pt= " << pt << " GeV";
    if (pt < ptCut_)
      continue;
    double eta = mu->eta();
    LogTrace("") << "\t... eta= " << eta;
    if (fabs(eta) > etaCut_)
      continue;
    bool iso = (*isoMap)[mu];
    LogTrace("") << "\t... isolated? " << iso;
    if (!iso)
      continue;
 
    for (unsigned int j = i + 1; j < muonCollection->size(); j++) {
      TrackRef mu2(muonCollection, j);
      LogTrace("") << "> Processing second muon number " << j << "...";
      double pt2 = mu2->pt();
      LogTrace("") << "\t... pt2= " << pt2 << " GeV";
      if (pt2 < ptCut_)
        continue;
      double eta2 = mu2->eta();
      LogTrace("") << "\t... eta2= " << eta2;
      if (fabs(eta2) > etaCut_)
        continue;
      bool iso2 = (*isoMap)[mu2];
      LogTrace("") << "\t... isolated2? " << iso2;
      if (!iso2)
        continue;
 
      double z_en = mu->p() + mu2->p();
      double z_px = mu->px() + mu2->px();
      double z_py = mu->py() + mu2->py();
      double z_pz = mu->pz() + mu2->pz();
      double massZ = z_en * z_en - z_px * z_px - z_py * z_py - z_pz * z_pz;
      massZ = (massZ > 0) ? sqrt(massZ) : 0;
      LogTrace("") << "\t... Z_en, Z_px, Z_py, Z_pz= " << z_en << ", " << z_px << ", " << z_py << ", " << z_pz
                   << " GeV";
      LogTrace("") << "\t... (GM-GM) Invariant reconstructed mass= " << massZ << " GeV";
      if (massZ < massZMin_)
        continue;
      if (massZ > massZMax_)
        continue;
      globalCombinationFound = true;
      npairs++;
    }
 
    if (onlyGlobalMuons_ || globalCombinationFound)
      continue;
 
    for (unsigned int j = 0; j < trackerCollection->size(); j++) {
      TrackRef mu2(trackerCollection, j);
      LogTrace("") << "> Processing track number " << j << "...";
      double pt2 = mu2->pt();
      LogTrace("") << "\t... pt3= " << pt2 << " GeV";
      if (pt2 < ptCut_)
        continue;
      double eta2 = mu2->eta();
      LogTrace("") << "\t... eta3= " << eta2;
      if (fabs(eta2) > etaCut_)
        continue;
      int nhits2 = mu2->numberOfValidHits();
      LogTrace("") << "\t... nhits3= " << nhits2;
      if (nhits2 < minTrackerHits_)
        continue;
      bool iso2 = (*isoTrackerMap)[mu2];
      LogTrace("") << "\t... isolated3? " << iso2;
      if (!iso2)
        continue;
 
      double z_en = mu->p() + mu2->p();
      double z_px = mu->px() + mu2->px();
      double z_py = mu->py() + mu2->py();
      double z_pz = mu->pz() + mu2->pz();
      double massZ = z_en * z_en - z_px * z_px - z_py * z_py - z_pz * z_pz;
      massZ = (massZ > 0) ? sqrt(massZ) : 0;
      LogTrace("") << "\t... Z_en, Z_px, Z_py, Z_pz= " << z_en << ", " << z_px << ", " << z_py << ", " << z_pz
                   << " GeV";
      LogTrace("") << "\t... (GM-TK) Invariant reconstructed mass= " << massZ << " GeV";
      if (massZ < massZMin_)
        continue;
      if (massZ > massZMax_)
        continue;
      npairs++;
    }
  }
 
  LogTrace("") << "> Number of Z pairs found= " << npairs;
  if (npairs < 1) {
    LogTrace("") << ">>>> Event REJECTED";
    return false;
  }
  LogTrace("") << ">>>> Event SELECTED!!!";
 
  return true;
}

References PVValHelper::eta, HLT_FULL_cff::eta2, etaCut_, ev, mps_fire::i, isoToken_, isoTrackerToken_, edm::HandleBase::isValid(), dqmiolumiharvest::j, LogTrace, massZMax_, massZMin_, minTrackerHits_, amptDefaultParameters_cff::mu, pdwgLeptonRecoSkim_cfi::muonCollection, muonToken_, onlyGlobalMuons_, DiDispStaMuonMonitor_cfi::pt, HLT_FULL_cff::pt2, ptCut_, mathSSE::sqrt(), and trackerToken_.

Member Data Documentation

etaCut_

double ZToMuMuSelector::etaCut_

private

Definition at line 21 of file ZToMuMuSelector.cc.

Referenced by filter().

isoToken_

edm::EDGetTokenT<edm::ValueMap<bool> > ZToMuMuSelector::isoToken_

private

Definition at line 19 of file ZToMuMuSelector.cc.

Referenced by filter().

isoTrackerToken_

edm::EDGetTokenT<edm::ValueMap<bool> > ZToMuMuSelector::isoTrackerToken_

private

Definition at line 27 of file ZToMuMuSelector.cc.

Referenced by filter().

massZMax_

double ZToMuMuSelector::massZMax_

private

Definition at line 23 of file ZToMuMuSelector.cc.

Referenced by filter().

massZMin_

double ZToMuMuSelector::massZMin_

private

Definition at line 22 of file ZToMuMuSelector.cc.

Referenced by filter().

minTrackerHits_

int ZToMuMuSelector::minTrackerHits_

private

Definition at line 28 of file ZToMuMuSelector.cc.

Referenced by filter().

muonToken_

edm::EDGetTokenT<reco::TrackCollection> ZToMuMuSelector::muonToken_

private

Definition at line 18 of file ZToMuMuSelector.cc.

Referenced by filter().

onlyGlobalMuons_

bool ZToMuMuSelector::onlyGlobalMuons_

private

Definition at line 25 of file ZToMuMuSelector.cc.

Referenced by filter().

ptCut_

double ZToMuMuSelector::ptCut_

private

Definition at line 20 of file ZToMuMuSelector.cc.

Referenced by badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons(), and filter().

trackerToken_

edm::EDGetTokenT<reco::TrackCollection> ZToMuMuSelector::trackerToken_

private

Definition at line 26 of file ZToMuMuSelector.cc.

Referenced by filter().