#include <HLTDisplacedtktktkVtxProducer.h>

Inheritance diagram for HLTDisplacedtktktkVtxProducer:

Public Member Functions
	HLTDisplacedtktktkVtxProducer (const edm::ParameterSet &)

void	produce (edm::Event &, const edm::EventSetup &) override

	~HLTDisplacedtktktkVtxProducer () override

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

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

Private Attributes
const int	chargeOpt_

double	firstTrackMass

double	firstTrackMass2

double	firstTrackPt

const double	massParticle1_

const double	massParticle2_

const double	massParticle3_

const double	maxEta_

const double	maxInvMass_

const double	maxInvMassRes_

const double	minInvMass_

const double	minInvMassRes_

const double	minPtRes_

const double	minPtTk1_

const double	minPtTk2_

const double	minPtTk3_

const double	minPtTri_

const edm::InputTag	previousCandTag_

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

const int	resOpt_

double	secondTrackMass

double	secondTrackMass2

double	secondTrackPt

const edm::InputTag	srcTag_

const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	srcToken_

double	thirdTrackMass

double	thirdTrackMass2

double	thirdTrackPt

const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	transientTrackRecordToken_

const int	triggerTypeDaughters_

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 35 of file HLTDisplacedtktktkVtxProducer.h.

Constructor & Destructor Documentation

◆ HLTDisplacedtktktkVtxProducer()

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

explicit

Definition at line 29 of file HLTDisplacedtktktkVtxProducer.cc.

References firstTrackMass, firstTrackMass2, firstTrackPt, massParticle1_, massParticle2_, massParticle3_, minPtTk1_, minPtTk2_, minPtTk3_, resOpt_, secondTrackMass, secondTrackMass2, secondTrackPt, edm::swap(), thirdTrackMass, thirdTrackMass2, and thirdTrackPt.

     : transientTrackRecordToken_(esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))),
       srcTag_(iConfig.getParameter<edm::InputTag>("Src")),
       srcToken_(consumes<reco::RecoChargedCandidateCollection>(srcTag_)),
       previousCandTag_(iConfig.getParameter<edm::InputTag>("PreviousCandTag")),
       previousCandToken_(consumes<trigger::TriggerFilterObjectWithRefs>(previousCandTag_)),
       maxEta_(iConfig.getParameter<double>("MaxEtaTk")),
       minPtTk1_(iConfig.getParameter<double>("MinPtResTk1")),
       minPtTk2_(iConfig.getParameter<double>("MinPtResTk2")),
       minPtTk3_(iConfig.getParameter<double>("MinPtThirdTk")),
       minPtRes_(iConfig.getParameter<double>("MinPtRes")),
       minPtTri_(iConfig.getParameter<double>("MinPtTri")),
       minInvMassRes_(iConfig.getParameter<double>("MinInvMassRes")),
       maxInvMassRes_(iConfig.getParameter<double>("MaxInvMassRes")),
       minInvMass_(iConfig.getParameter<double>("MinInvMass")),
       maxInvMass_(iConfig.getParameter<double>("MaxInvMass")),
       massParticle1_(iConfig.getParameter<double>("massParticleRes1")),
       massParticle2_(iConfig.getParameter<double>("massParticleRes2")),
       massParticle3_(iConfig.getParameter<double>("massParticle3")),
       chargeOpt_(iConfig.getParameter<int>("ChargeOpt")),
       resOpt_(iConfig.getParameter<int>("ResOpt")),
       triggerTypeDaughters_(iConfig.getParameter<int>("triggerTypeDaughters")) {
   produces<VertexCollection>();
 
   firstTrackMass = massParticle1_;
   secondTrackMass = massParticle2_;
   thirdTrackMass = massParticle3_;
   firstTrackPt = minPtTk1_;
   secondTrackPt = minPtTk2_;
   thirdTrackPt = minPtTk3_;
   if (resOpt_ <= 0 && massParticle1_ != massParticle2_) {
     if (massParticle1_ == massParticle3_) {
       std::swap(secondTrackMass, thirdTrackMass);
       std::swap(secondTrackPt, thirdTrackPt);
     }
     if (massParticle2_ == massParticle3_) {
       std::swap(firstTrackMass, thirdTrackMass);
       std::swap(firstTrackPt, thirdTrackPt);
     }
   }
   firstTrackMass2 = firstTrackMass * firstTrackMass;
   secondTrackMass2 = secondTrackMass * secondTrackMass;
   thirdTrackMass2 = thirdTrackMass * thirdTrackMass;
 }

◆ ~HLTDisplacedtktktkVtxProducer()

HLTDisplacedtktktkVtxProducer::~HLTDisplacedtktktkVtxProducer ( )

overridedefault

Member Function Documentation

◆ checkPreviousCand()

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

private

Definition at line 268 of file HLTDisplacedtktktkVtxProducer.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;
 }

◆ fillDescriptions()

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

static

Definition at line 76 of file HLTDisplacedtktktkVtxProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("Src", edm::InputTag("hltL3MuonCandidates"));
   desc.add<edm::InputTag>("PreviousCandTag", edm::InputTag(""));
   desc.add<double>("MaxEtaTk", 2.5);
   desc.add<double>("MinPtResTk1", 0.0);
   desc.add<double>("MinPtResTk2", 0.0);
   desc.add<double>("MinPtThirdTk", 0.0);
   desc.add<double>("MinPtRes", 0.0);
   desc.add<double>("MinPtTri", 0.0);
   desc.add<double>("MinInvMassRes", 1.0);
   desc.add<double>("MaxInvMassRes", 20.0);
   desc.add<double>("MinInvMass", 1.0);
   desc.add<double>("MaxInvMass", 20.0);
   desc.add<double>("massParticleRes1", 0.4937);
   desc.add<double>("massParticleRes2", 0.4937);
   desc.add<double>("massParticle3", 0.1396);
   desc.add<int>("ChargeOpt", -1);
   desc.add<int>("ResOpt", 1);
   desc.add<int>("triggerTypeDaughters", 0);
 
   descriptions.add("hltDisplacedtktktkVtxProducer", desc);
 }

◆ produce()

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

override

Definition at line 101 of file HLTDisplacedtktktkVtxProducer.cc.

References funct::abs(), chargeOpt_, checkPreviousCand(), StorageManager_cfg::e1, SiPixelPhase1Clusters_cfi::e3, firstTrackMass, firstTrackMass2, firstTrackPt, edm::Ref< C, T, F >::get(), edm::EventSetup::getHandle(), trigger::TriggerRefsCollections::getObjects(), iEvent, TransientVertex::isValid(), LogDebug, maxEta_, maxInvMass_, maxInvMassRes_, minInvMass_, minInvMassRes_, minPtRes_, minPtTri_, eostools::move(), AlCaHLTBitMon_ParallelJobs::p, previousCandToken_, resOpt_, secondTrackMass, secondTrackMass2, secondTrackPt, mathSSE::sqrt(), srcToken_, thirdTrackMass, thirdTrackMass2, thirdTrackPt, transientTrackRecordToken_, triggerTypeDaughters_, bphysicsOniaDQM_cfi::vertex, KalmanVertexFitter::vertex(), and spclusmultinvestigator_cfi::vertexCollection.

                                                                                          {
   // get hold of track trks
   Handle<RecoChargedCandidateCollection> trackcands;
   iEvent.getByToken(srcToken_, trackcands);
   if (trackcands->size() < 3)
     return;
 
   // get the transient track builder
   auto const& theB = iSetup.getHandle(transientTrackRecordToken_);
 
   std::unique_ptr<VertexCollection> vertexCollection(new VertexCollection());
 
   // look at all trackcands,  check cuts and make vertices
   double e1, e2, e3;
   Particle::LorentzVector p, pres;
 
   RecoChargedCandidateCollection::const_iterator cand1;
   RecoChargedCandidateCollection::const_iterator cand2;
   RecoChargedCandidateCollection::const_iterator cand3;
 
   // get the objects passing the previous filter
   Handle<TriggerFilterObjectWithRefs> previousCands;
   iEvent.getByToken(previousCandToken_, previousCands);
 
   vector<RecoChargedCandidateRef> vPrevCands;
   previousCands->getObjects(triggerTypeDaughters_, vPrevCands);
 
   std::vector<bool> candComp;
   for (cand1 = trackcands->begin(); cand1 != trackcands->end(); cand1++)
     candComp.push_back(checkPreviousCand(cand1->get<TrackRef>(), vPrevCands));
 
   for (cand1 = trackcands->begin(); cand1 != trackcands->end(); cand1++) {
     TrackRef tk1 = cand1->get<TrackRef>();
     LogDebug("HLTDisplacedtktktkVtxProducer") << " 1st track in loop: q*pt= " << cand1->charge() * cand1->pt()
                                               << ", eta= " << cand1->eta() << ", hits= " << tk1->numberOfValidHits();
 
     //first check if this track passed the previous filter
     if (!candComp[cand1 - trackcands->begin()])
       continue;
     // if( ! checkPreviousCand( tk1, vPrevCands) ) continue;
 
     // cuts
     if (std::abs(cand1->eta()) > maxEta_)
       continue;
     if (cand1->pt() < firstTrackPt)
       continue;
 
     cand2 = trackcands->begin();
     if (firstTrackMass == secondTrackMass) {
       cand2 = cand1 + 1;
     }
 
     for (; cand2 != trackcands->end(); cand2++) {
       TrackRef tk2 = cand2->get<TrackRef>();
       if (tk1 == tk2)
         continue;
       LogDebug("HLTDisplacedtktktkVtxProducer")
           << " 2nd track in loop: q*pt= " << cand2->charge() * cand2->pt() << ", eta= " << cand2->eta()
           << ", hits= " << tk2->numberOfValidHits() << ", d0= " << tk2->d0();
 
       //first check if this track passed the previous filter
       if (!candComp[cand2 - trackcands->begin()])
         continue;
       // if( ! checkPreviousCand( tk2, vPrevCands) ) continue;
 
       // cuts
       if (std::abs(cand2->eta()) > maxEta_)
         continue;
       if (cand2->pt() < secondTrackPt)
         continue;
 
       // opposite sign or same sign for resonance
       if (resOpt_ > 0) {
         if (chargeOpt_ < 0) {
           if (cand1->charge() * cand2->charge() > 0)
             continue;
         } else if (chargeOpt_ > 0) {
           if (cand1->charge() * cand2->charge() < 0)
             continue;
         }
       }
 
       //
       // Combined ditrack system
       e1 = sqrt(cand1->momentum().Mag2() + firstTrackMass2);
       e2 = sqrt(cand2->momentum().Mag2() + secondTrackMass2);
       pres = Particle::LorentzVector(cand1->px(), cand1->py(), cand1->pz(), e1) +
              Particle::LorentzVector(cand2->px(), cand2->py(), cand2->pz(), e2);
 
       if (resOpt_ > 0) {
         if (pres.pt() < minPtRes_)
           continue;
         double invmassRes = std::abs(pres.mass());
         LogDebug("HLTDisplacedtktktkVtxProducer") << " ... 1-2 invmass= " << invmassRes;
         if (invmassRes < minInvMassRes_)
           continue;
         if (invmassRes > maxInvMassRes_)
           continue;
       }
 
       cand3 = trackcands->begin();
       if (firstTrackMass == secondTrackMass && firstTrackMass == thirdTrackMass && resOpt_ <= 0) {
         cand3 = cand2 + 1;
       }
 
       for (; cand3 != trackcands->end(); cand3++) {
         TrackRef tk3 = cand3->get<TrackRef>();
         if (tk1 == tk3 || tk2 == tk3)
           continue;
         LogDebug("HLTDisplacedtktktkVtxProducer")
             << " 3rd track in loop: q*pt= " << cand3->charge() * cand3->pt() << ", eta= " << cand3->eta()
             << ", hits= " << tk3->numberOfValidHits();
 
         //first check if this track passed the previous filter
         if (!candComp[cand3 - trackcands->begin()])
           continue;
         // if( ! checkPreviousCand( tk3, vPrevCands) ) continue;
 
         // cuts
         if (std::abs(cand3->eta()) > maxEta_)
           continue;
         if (cand3->pt() < thirdTrackPt)
           continue;
 
         e3 = sqrt(cand3->momentum().Mag2() + thirdTrackMass2);
         p = Particle::LorentzVector(cand1->px(), cand1->py(), cand1->pz(), e1) +
             Particle::LorentzVector(cand2->px(), cand2->py(), cand2->pz(), e2) +
             Particle::LorentzVector(cand3->px(), cand3->py(), cand3->pz(), e3);
 
         if (p.pt() < minPtTri_)
           continue;
         double invmass = std::abs(p.mass());
         LogDebug("HLTDisplacedtktktkVtxProducer") << " ... 1-2-3 invmass= " << invmass;
         if (invmass < minInvMass_)
           continue;
         if (invmass > maxInvMass_)
           continue;
 
         // do the vertex fit
         vector<TransientTrack> t_tks;
         TransientTrack ttkp1 = (*theB).build(&tk1);
         TransientTrack ttkp2 = (*theB).build(&tk2);
         TransientTrack ttkp3 = (*theB).build(&tk3);
 
         t_tks.push_back(ttkp1);
         t_tks.push_back(ttkp2);
         t_tks.push_back(ttkp3);
 
         if (t_tks.size() != 3)
           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));
 }