HLTmumutkVtxProducer Class Reference

#include <HLTmumutkVtxProducer.h>

Inheritance diagram for HLTmumutkVtxProducer:

Public Member Functions
	HLTmumutkVtxProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~HLTmumutkVtxProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
bool	checkPreviousCand (const reco::TrackRef &trackref, const std::vector< reco::RecoChargedCandidateRef > &ref2) const
bool	overlap (const reco::TrackRef &trackref1, const reco::TrackRef &trackref2)

Static Private Member Functions
static FreeTrajectoryState	initialFreeState (const reco::Track &, const MagneticField *)

Private Attributes
const edm::InputTag	beamSpotTag_
const edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
const double	maxEta_
const double	maxInvMass_
const std::string	mfName_
const double	minD0Significance_
const double	minInvMass_
const double	minPt_
const edm::InputTag	muCandTag_
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	muCandToken_
const double	overlapDR_
const edm::InputTag	previousCandTag_
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	previousCandToken_
const double	thirdTrackMass_
const edm::InputTag	trkCandTag_
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	trkCandToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

HLT Filter for b to (mumu) + X

Implementation: <Notes on="" implementation>="">

Definition at line 43 of file HLTmumutkVtxProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 31 of file HLTmumutkVtxProducer.cc.

References ~HLTmumutkVtxProducer().

    : muCandTag_(iConfig.getParameter<edm::InputTag>("MuCand")),
      muCandToken_(consumes<reco::RecoChargedCandidateCollection>(muCandTag_)),
      trkCandTag_(iConfig.getParameter<edm::InputTag>("TrackCand")),
      trkCandToken_(consumes<reco::RecoChargedCandidateCollection>(trkCandTag_)),
      previousCandTag_(iConfig.getParameter<edm::InputTag>("PreviousCandTag")),
      previousCandToken_(consumes<trigger::TriggerFilterObjectWithRefs>(previousCandTag_)),
      mfName_(iConfig.getParameter<std::string>("SimpleMagneticField")),
      thirdTrackMass_(iConfig.getParameter<double>("ThirdTrackMass")),
      maxEta_(iConfig.getParameter<double>("MaxEta")),
      minPt_(iConfig.getParameter<double>("MinPt")),
      minInvMass_(iConfig.getParameter<double>("MinInvMass")),
      maxInvMass_(iConfig.getParameter<double>("MaxInvMass")),
      minD0Significance_(iConfig.getParameter<double>("MinD0Significance")),
      overlapDR_(iConfig.getParameter<double>("OverlapDR")),
      beamSpotTag_(iConfig.getParameter<edm::InputTag>("BeamSpotTag")),
      beamSpotToken_(consumes<reco::BeamSpot>(beamSpotTag_)) {
  produces<VertexCollection>();
}

HLTmumutkVtxProducer::~HLTmumutkVtxProducer ( )

overridedefault

Referenced by HLTmumutkVtxProducer().

Member Function Documentation

bool HLTmumutkVtxProducer::checkPreviousCand ( const reco::TrackRef &  trackref, const std::vector< reco::RecoChargedCandidateRef > &  ref2  ) const

private

Definition at line 237 of file HLTmumutkVtxProducer.cc.

References mps_fire::i, and convertSQLiteXML::ok.

Referenced by produce().

                                                                                                   {
  bool ok = false;
  for (auto& i : refVect) {
    if (i->get<TrackRef>() == trackref) {
      ok = true;
      break;
    }
  }
  return ok;
}

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

static

Definition at line 54 of file HLTmumutkVtxProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                      {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("MuCand", edm::InputTag("hltMuTracks"));
  desc.add<edm::InputTag>("TrackCand", edm::InputTag("hltMumukAllConeTracks"));
  desc.add<edm::InputTag>("PreviousCandTag", edm::InputTag("hltDisplacedmumuFilterDoubleMu4Jpsi"));
  desc.add<std::string>("SimpleMagneticField", "");
  desc.add<double>("ThirdTrackMass", 0.493677);
  desc.add<double>("MaxEta", 2.5);
  desc.add<double>("MinPt", 3.0);
  desc.add<double>("MinInvMass", 0.0);
  desc.add<double>("MaxInvMass", 99999.);
  desc.add<double>("MinD0Significance", 0.0);
  desc.add<double>("OverlapDR", 1.44e-4);
  desc.add<edm::InputTag>("BeamSpotTag", edm::InputTag("hltOfflineBeamSpot"));
  descriptions.add("HLTmumutkVtxProducer", desc);
}

FreeTrajectoryState HLTmumutkVtxProducer::initialFreeState ( const reco::Track &  tk, const MagneticField *  field  )

staticprivate

Definition at line 221 of file HLTmumutkVtxProducer.cc.

References reco::TrackBase::charge(), reco::TrackBase::covariance(), runTheMatrix::err, reco::TrackBase::momentum(), and reco::TrackBase::vertex().

Referenced by produce().

                                                                                                          {
  Basic3DVector<float> pos(tk.vertex());
  GlobalPoint gpos(pos);
  Basic3DVector<float> mom(tk.momentum());
  GlobalVector gmom(mom);
  GlobalTrajectoryParameters par(gpos, gmom, tk.charge(), field);
  CurvilinearTrajectoryError err(tk.covariance());
  return FreeTrajectoryState(par, err);
}

bool HLTmumutkVtxProducer::overlap ( const reco::TrackRef &  trackref1, const reco::TrackRef &  trackref2  )

private

Definition at line 231 of file HLTmumutkVtxProducer.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, and overlapDR_.

Referenced by produce().

                                                                                       {
  if (deltaR(trackref1->eta(), trackref1->phi(), trackref2->eta(), trackref2->phi()) < overlapDR_)
    return true;
  return false;
}

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

override

Definition at line 72 of file HLTmumutkVtxProducer.cc.

References funct::abs(), beamSpotToken_, checkPreviousCand(), StorageManager_cfg::e1, SiPixelPhase1Clusters_cfi::e3, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), initialFreeState(), TransientVertex::isValid(), LogDebug, maxEta_, maxInvMass_, mfName_, minD0Significance_, minInvMass_, minPt_, eostools::move(), muCandToken_, overlap(), AlCaHLTBitMon_ParallelJobs::p, p1, p2, p3, previousCandToken_, edm::ESHandle< T >::product(), edm::Event::put(), Measurement1D::significance(), mathSSE::sqrt(), thirdTrackMass_, TrajectoryStateClosestToBeamLine::transverseImpactParameter(), trigger::TriggerMuon, trkCandToken_, bphysicsOniaDQM_cfi::vertex, KalmanVertexFitter::vertex(), and spclusmultinvestigator_cfi::vertexCollection.

                                                                                {
  const double MuMass(0.106);
  const double MuMass2(MuMass * MuMass);
  const double thirdTrackMass2(thirdTrackMass_ * thirdTrackMass_);

  // get hold of muon trks
  Handle<RecoChargedCandidateCollection> mucands;
  iEvent.getByToken(muCandToken_, mucands);

  //get the transient track builder:
  edm::ESHandle<TransientTrackBuilder> theB;
  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", theB);

  //get the beamspot position
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  iEvent.getByToken(beamSpotToken_, recoBeamSpotHandle);

  //get the b field
  ESHandle<MagneticField> bFieldHandle;
  iSetup.get<IdealMagneticFieldRecord>().get(mfName_, bFieldHandle);
  const MagneticField* magField = bFieldHandle.product();
  TSCBLBuilderNoMaterial blsBuilder;

  // get track candidates around displaced muons
  Handle<RecoChargedCandidateCollection> trkcands;
  iEvent.getByToken(trkCandToken_, trkcands);

  unique_ptr<VertexCollection> vertexCollection(new VertexCollection());

  // Ref to Candidate object to be recorded in filter object
  RecoChargedCandidateRef refMu1;
  RecoChargedCandidateRef refMu2;
  RecoChargedCandidateRef refTrk;

  double e1, e2, e3;
  Particle::LorentzVector p, p1, p2, p3;

  if (mucands->size() < 2)
    return;
  if (trkcands->empty())
    return;

  RecoChargedCandidateCollection::const_iterator mucand1;
  RecoChargedCandidateCollection::const_iterator mucand2;
  RecoChargedCandidateCollection::const_iterator trkcand;

  // get the objects passing the previous filter
  Handle<TriggerFilterObjectWithRefs> previousCands;
  iEvent.getByToken(previousCandToken_, previousCands);

  vector<RecoChargedCandidateRef> vPrevCands;
  previousCands->getObjects(TriggerMuon, vPrevCands);

  for (mucand1 = mucands->begin(); mucand1 != mucands->end(); ++mucand1) {
    TrackRef trk1 = mucand1->get<TrackRef>();
    LogDebug("HLTmumutkVtxProducer") << " 1st muon: q*pt= " << trk1->charge() * trk1->pt() << ", eta= " << trk1->eta()
                                     << ", hits= " << trk1->numberOfValidHits();

    //first check if this muon passed the previous filter
    if (!checkPreviousCand(trk1, vPrevCands))
      continue;

    // eta and pt cut
    if (fabs(trk1->eta()) > maxEta_)
      continue;
    if (trk1->pt() < minPt_)
      continue;

    mucand2 = mucand1;
    ++mucand2;
    for (; mucand2 != mucands->end(); mucand2++) {
      TrackRef trk2 = mucand2->get<TrackRef>();

      LogDebug("HLTDisplacedMumukFilter") << " 2nd muon: q*pt= " << trk2->charge() * trk2->pt()
                                          << ", eta= " << trk2->eta() << ", hits= " << trk2->numberOfValidHits();

      //first check if this muon passed the previous filter
      if (!checkPreviousCand(trk2, vPrevCands))
        continue;
      // eta and pt cut
      if (fabs(trk2->eta()) > maxEta_)
        continue;
      if (trk2->pt() < minPt_)
        continue;

      //loop on track collection
      for (trkcand = trkcands->begin(); trkcand != trkcands->end(); ++trkcand) {
        TrackRef trk3 = trkcand->get<TrackRef>();
        if (overlap(trk1, trk3))
          continue;
        if (overlap(trk2, trk3))
          continue;

        LogDebug("HLTDisplacedMumukFilter") << " 3rd track: q*pt= " << trk3->charge() * trk3->pt()
                                            << ", eta= " << trk3->eta() << ", hits= " << trk3->numberOfValidHits();

        // eta and pt cut
        if (fabs(trk3->eta()) > maxEta_)
          continue;
        if (trk3->pt() < minPt_)
          continue;

        // Combined system
        e1 = sqrt(trk1->momentum().Mag2() + MuMass2);
        e2 = sqrt(trk2->momentum().Mag2() + MuMass2);
        e3 = sqrt(trk3->momentum().Mag2() + thirdTrackMass2);

        p1 = Particle::LorentzVector(trk1->px(), trk1->py(), trk1->pz(), e1);
        p2 = Particle::LorentzVector(trk2->px(), trk2->py(), trk2->pz(), e2);
        p3 = Particle::LorentzVector(trk3->px(), trk3->py(), trk3->pz(), e3);

        p = p1 + p2 + p3;

        //invariant mass cut
        double invmass = abs(p.mass());
        LogDebug("HLTDisplacedMumukFilter") << " Invmass= " << invmass;
        if (invmass < minInvMass_)
          continue;
        if (invmass > maxInvMass_)
          continue;

        // do the vertex fit
        vector<TransientTrack> t_tks;
        t_tks.push_back((*theB).build(&trk1));
        t_tks.push_back((*theB).build(&trk2));
        t_tks.push_back((*theB).build(&trk3));
        if (t_tks.size() != 3)
          continue;

        FreeTrajectoryState InitialFTS = initialFreeState(*trk3, magField);
        TrajectoryStateClosestToBeamLine tscb(blsBuilder(InitialFTS, *recoBeamSpotHandle));
        double d0sig = tscb.transverseImpactParameter().significance();
        if (d0sig < minD0Significance_)
          continue;

        KalmanVertexFitter kvf;
        TransientVertex tv = kvf.vertex(t_tks);
        if (!tv.isValid())
          continue;
        Vertex vertex = tv;

        // put vertex in the event
        vertexCollection->push_back(vertex);
      }
    }
  }
  iEvent.put(std::move(vertexCollection));
}

Member Data Documentation

const edm::InputTag HLTmumutkVtxProducer::beamSpotTag_

private

Definition at line 72 of file HLTmumutkVtxProducer.h.

const edm::EDGetTokenT<reco::BeamSpot> HLTmumutkVtxProducer::beamSpotToken_

private

Definition at line 73 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::maxEta_

private

Definition at line 65 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::maxInvMass_

private

Definition at line 68 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const std::string HLTmumutkVtxProducer::mfName_

private

Definition at line 62 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::minD0Significance_

private

Definition at line 69 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::minInvMass_

private

Definition at line 67 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::minPt_

private

Definition at line 66 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const edm::InputTag HLTmumutkVtxProducer::muCandTag_

private

Definition at line 55 of file HLTmumutkVtxProducer.h.

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTmumutkVtxProducer::muCandToken_

private

Definition at line 56 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::overlapDR_

private

Definition at line 70 of file HLTmumutkVtxProducer.h.

Referenced by overlap().

const edm::InputTag HLTmumutkVtxProducer::previousCandTag_

private

Definition at line 59 of file HLTmumutkVtxProducer.h.

const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTmumutkVtxProducer::previousCandToken_

private

Definition at line 60 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::thirdTrackMass_

private

Definition at line 64 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const edm::InputTag HLTmumutkVtxProducer::trkCandTag_

private

Definition at line 57 of file HLTmumutkVtxProducer.h.

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTmumutkVtxProducer::trkCandToken_

private

Definition at line 58 of file HLTmumutkVtxProducer.h.

Referenced by produce().