HLTDisplacedmumumuFilter Class Reference

#include <HLTDisplacedmumumuFilter.h>

Inheritance diagram for HLTDisplacedmumumuFilter:

Public Member Functions
void	beginJob () override
void	endJob () override
	HLTDisplacedmumumuFilter (const edm::ParameterSet &)
bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
	~HLTDisplacedmumumuFilter () override
Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)
int	module (edm::Event const &) const
const std::string *	moduleLabel () const
int	path (edm::Event const &) const
const std::string *	pathName (edm::Event const &) const
std::pair< int, int >	pmid (edm::Event const &) const
bool	saveTags () const
	~HLTFilter () override
Public Member Functions inherited from edm::global::EDFilter<>
	EDFilter ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsStreamLuminosityBlocks () const final
bool	wantsStreamRuns () const final
Public Member Functions inherited from edm::global::EDFilterBase
	EDFilterBase ()
ModuleDescription const &	moduleDescription () const
	~EDFilterBase () 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
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)
	~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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)
Static Public Member Functions inherited from edm::global::EDFilterBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::InputTag	beamSpotTag_
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
edm::InputTag	DisplacedVertexTag_
edm::EDGetTokenT< reco::VertexCollection >	DisplacedVertexToken_
bool	fastAccept_
double	maxLxySignificance_
double	maxNormalisedChi2_
double	minCosinePointingAngle_
double	minLxySignificance_
double	minVtxProbability_
edm::InputTag	MuonTag_
edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	MuonToken_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	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
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
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<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 12 of file HLTDisplacedmumumuFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file HLTDisplacedmumumuFilter.cc.

References ~HLTDisplacedmumumuFilter().

                                                                                 :
  HLTFilter(iConfig),
  fastAccept_ (iConfig.getParameter<bool>("FastAccept")),
  minLxySignificance_ (iConfig.getParameter<double>("MinLxySignificance")),
  maxLxySignificance_ (iConfig.getParameter<double>("MaxLxySignificance")),
  maxNormalisedChi2_ (iConfig.getParameter<double>("MaxNormalisedChi2")),
  minVtxProbability_ (iConfig.getParameter<double>("MinVtxProbability")),
  minCosinePointingAngle_ (iConfig.getParameter<double>("MinCosinePointingAngle")),
  DisplacedVertexTag_(iConfig.getParameter<edm::InputTag>("DisplacedVertexTag")),
  DisplacedVertexToken_(consumes<reco::VertexCollection>(DisplacedVertexTag_)),
  beamSpotTag_ (iConfig.getParameter<edm::InputTag> ("BeamSpotTag")),
  beamSpotToken_(consumes<reco::BeamSpot>(beamSpotTag_)),
  MuonTag_ (iConfig.getParameter<edm::InputTag>("MuonTag")),
  MuonToken_(consumes<reco::RecoChargedCandidateCollection>(MuonTag_))
{
}

HLTDisplacedmumumuFilter::~HLTDisplacedmumumuFilter ( )

overridedefault

Referenced by HLTDisplacedmumumuFilter().

Member Function Documentation

void HLTDisplacedmumumuFilter::beginJob ( void  )

overridevirtual

Reimplemented from edm::global::EDFilterBase.

Definition at line 72 of file HLTDisplacedmumumuFilter.cc.

{

}

void HLTDisplacedmumumuFilter::endJob ( void  )

overridevirtual

Reimplemented from edm::global::EDFilterBase.

Definition at line 78 of file HLTDisplacedmumumuFilter.cc.

{
    
}

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

static

Definition at line 56 of file HLTDisplacedmumumuFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLTFilter::makeHLTFilterDescription().

                                                                                          {
  edm::ParameterSetDescription desc;
  makeHLTFilterDescription(desc);
  desc.add<bool>("FastAccept",false);
  desc.add<double>("MinLxySignificance",0.0);
  desc.add<double>("MaxLxySignificance",0.0);
  desc.add<double>("MaxNormalisedChi2",10.0);
  desc.add<double>("MinVtxProbability",0.0);
  desc.add<double>("MinCosinePointingAngle",-2.0);
  desc.add<edm::InputTag>("DisplacedVertexTag",edm::InputTag("hltDisplacedmumumuVtx"));
  desc.add<edm::InputTag>("BeamSpotTag",edm::InputTag("hltOnlineBeamSpot"));
  desc.add<edm::InputTag>("MuonTag",edm::InputTag("hltL3MuonCandidates"));
  descriptions.add("hltDisplacedmumumuFilter", desc);
}

bool HLTDisplacedmumumuFilter::hltFilter ( edm::Event &  iEvent, const edm::EventSetup &  iSetup, trigger::TriggerFilterObjectWithRefs &  filterproduct  ) const

overridevirtual

Implements HLTFilter.

Definition at line 84 of file HLTDisplacedmumumuFilter.cc.

References trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), beamSpotToken_, DisplacedVertexToken_, SoftLeptonByDistance_cfi::distance, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), Exception, edm::Event::getByToken(), LogDebug, maxLxySignificance_, maxNormalisedChi2_, minCosinePointingAngle_, minLxySignificance_, minVtxProbability_, MuonTag_, MuonToken_, HLTFilter::saveTags(), mathSSE::sqrt(), trigger::TriggerMuon, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

{


  // get beam spot
  reco::BeamSpot vertexBeamSpot;
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  iEvent.getByToken(beamSpotToken_,recoBeamSpotHandle);
  vertexBeamSpot = *recoBeamSpotHandle;


  // get displaced vertices
  reco::VertexCollection displacedVertexColl;
  edm::Handle<reco::VertexCollection> displacedVertexCollHandle;
  bool foundVertexColl = iEvent.getByToken(DisplacedVertexToken_, displacedVertexCollHandle);
  if(foundVertexColl) displacedVertexColl = *displacedVertexCollHandle;


  // get muon collection
  edm::Handle<reco::RecoChargedCandidateCollection> mucands;
  iEvent.getByToken(MuonToken_,mucands);

  // All HLT filters must create and fill an HLT filter object,
  // recording any reconstructed physics objects satisfying (or not)
  // this HLT filter, and place it in the Event.


  // Ref to Candidate object to be recorded in filter object
  reco::RecoChargedCandidateRef ref1;
  reco::RecoChargedCandidateRef ref2;
  reco::RecoChargedCandidateRef ref3;

  if (saveTags()) filterproduct.addCollectionTag(MuonTag_);


  bool triggered = false;

  // loop over vertex collection
  for(auto displacedVertex : displacedVertexColl){
          // check if the vertex actually consists of exactly two muon tracks, throw exception if not
          if(displacedVertex.tracksSize() != 3)  throw cms::Exception("BadLogic") << "HLTDisplacedmumumuFilter: ERROR: the Jpsi vertex must have exactly three muons by definition. It now has n muons = "
                                                << displacedVertex.tracksSize() << std::endl;

          float normChi2 = displacedVertex.normalizedChi2();
      if (normChi2 > maxNormalisedChi2_) continue;

      double vtxProb = 0.0;
      if( (displacedVertex.chi2()>=0.0) && (displacedVertex.ndof()>0) ) vtxProb = TMath::Prob(displacedVertex.chi2(), displacedVertex.ndof() );
      if (vtxProb < minVtxProbability_) continue;

          // get the two muons from the vertex
          auto trackIt =  displacedVertex.tracks_begin();
          reco::TrackRef vertextkRef1 =  (*trackIt).castTo<reco::TrackRef>() ;
          // the second one
          trackIt++;
          reco::TrackRef vertextkRef2 =  (*trackIt).castTo<reco::TrackRef>();
          // the second one
          trackIt++;
          reco::TrackRef vertextkRef3 =  (*trackIt).castTo<reco::TrackRef>();

      reco::RecoChargedCandidateCollection::const_iterator cand1;
      reco::RecoChargedCandidateCollection::const_iterator cand2;   
      reco::RecoChargedCandidateCollection::const_iterator cand3;   

      // first find these two tracks in the muon collection
      int iFoundRefs = 0;
      for (auto cand=mucands->begin(); cand!=mucands->end(); cand++) {
        reco::TrackRef tkRef = cand->get<reco::TrackRef>();
        if(tkRef == vertextkRef1) {cand1 = cand; iFoundRefs++;}
        if(tkRef == vertextkRef2) {cand2 = cand; iFoundRefs++;}
        if(tkRef == vertextkRef3) {cand3 = cand; iFoundRefs++;}
      }
      if(iFoundRefs < 3) throw cms::Exception("BadLogic") << "HLTDisplacedmumumuFilter: ERROR: the muons matched with the Jpsi vertex tracks should be at least three by definition."  << std::endl;

          // calculate two-track transverse momentum
          math::XYZVector pperp(cand1->px() + cand2->px() + cand3->px(),
                cand1->py() + cand2->py() + cand3->py(),
                      0.);


      const reco::Vertex::Point& vpoint=displacedVertex.position();
      //translate to global point, should be improved
      GlobalPoint secondaryVertex (vpoint.x(), vpoint.y(), vpoint.z());

      reco::Vertex::Error verr = displacedVertex.error();
      // translate to global error, should be improved
      GlobalError err(verr.At(0,0), verr.At(1,0), verr.At(1,1), verr.At(2,0), verr.At(2,1), verr.At(2,2) );

      GlobalPoint displacementFromBeamspot( -1*((vertexBeamSpot.x0() -  secondaryVertex.x()) +  (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dxdz()),
                              -1*((vertexBeamSpot.y0() - secondaryVertex.y())+  (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dydz()), 0);

          float lxy = displacementFromBeamspot.perp();
          float lxyerr = sqrt(err.rerr(displacementFromBeamspot));


          //calculate the angle between the decay length and the mumumu momentum
      reco::Vertex::Point vperp(displacementFromBeamspot.x(),displacementFromBeamspot.y(),0.);

          float cosAlpha = vperp.Dot(pperp)/(vperp.R()*pperp.R());

          // check thresholds
          if (cosAlpha < minCosinePointingAngle_) continue;
          if (minLxySignificance_ > 0. && lxy/lxyerr < minLxySignificance_) continue;
      if (maxLxySignificance_ > 0. && lxy/lxyerr > maxLxySignificance_) continue;
      triggered = true;

      // now add the muons that passed to the filter object
    
      ref1=reco::RecoChargedCandidateRef( edm::Ref<reco::RecoChargedCandidateCollection>    (mucands,distance(mucands->begin(), cand1)));
      filterproduct.addObject(trigger::TriggerMuon,ref1);
    
      ref2=reco::RecoChargedCandidateRef( edm::Ref<reco::RecoChargedCandidateCollection> (mucands,distance(mucands->begin(),cand2 )));
      filterproduct.addObject(trigger::TriggerMuon,ref2);

      ref3=reco::RecoChargedCandidateRef( edm::Ref<reco::RecoChargedCandidateCollection> (mucands,distance(mucands->begin(),cand3 )));
      filterproduct.addObject(trigger::TriggerMuon,ref3);
  }

  LogDebug("HLTDisplacedMumumuFilter") << " >>>>> Result of HLTDisplacedMumumuFilter is "<< triggered <<std::endl;

  return triggered;
}

Member Data Documentation

edm::InputTag HLTDisplacedmumumuFilter::beamSpotTag_

private

Definition at line 31 of file HLTDisplacedmumumuFilter.h.

edm::EDGetTokenT<reco::BeamSpot> HLTDisplacedmumumuFilter::beamSpotToken_

private

Definition at line 32 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

edm::InputTag HLTDisplacedmumumuFilter::DisplacedVertexTag_

private

Definition at line 29 of file HLTDisplacedmumumuFilter.h.

edm::EDGetTokenT<reco::VertexCollection> HLTDisplacedmumumuFilter::DisplacedVertexToken_

private

Definition at line 30 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

bool HLTDisplacedmumumuFilter::fastAccept_

private

Definition at line 23 of file HLTDisplacedmumumuFilter.h.

double HLTDisplacedmumumuFilter::maxLxySignificance_

private

Definition at line 25 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

double HLTDisplacedmumumuFilter::maxNormalisedChi2_

private

Definition at line 26 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

double HLTDisplacedmumumuFilter::minCosinePointingAngle_

private

Definition at line 28 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

double HLTDisplacedmumumuFilter::minLxySignificance_

private

Definition at line 24 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

double HLTDisplacedmumumuFilter::minVtxProbability_

private

Definition at line 27 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

edm::InputTag HLTDisplacedmumumuFilter::MuonTag_

private

Definition at line 33 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().

edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTDisplacedmumumuFilter::MuonToken_

private

Definition at line 34 of file HLTDisplacedmumumuFilter.h.

Referenced by hltFilter().