#include <HLTDisplacedtktkVtxProducer.h>

Inheritance diagram for HLTDisplacedtktkVtxProducer:

Public Member Functions
virtual void	beginJob ()

virtual void	endJob ()

	HLTDisplacedtktkVtxProducer (const edm::ParameterSet &)

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

	~HLTDisplacedtktkVtxProducer ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Private Attributes
const int	chargeOpt_

const double	massParticle1_

const double	massParticle2_

const double	maxEta_

const double	maxInvMass_

const double	minInvMass_

const double	minPt_

const double	minPtPair_

const edm::InputTag	previousCandTag_

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

const edm::InputTag	srcTag_

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

const int	triggerTypeDaughters_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
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 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 35 of file HLTDisplacedtktkVtxProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 35 of file HLTDisplacedtktkVtxProducer.cc.

                                                                                       : 
     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>("MaxEta")),
     minPt_ (iConfig.getParameter<double>("MinPt")),
     minPtPair_ (iConfig.getParameter<double>("MinPtPair")),
     minInvMass_ (iConfig.getParameter<double>("MinInvMass")),
     maxInvMass_ (iConfig.getParameter<double>("MaxInvMass")),
     massParticle1_ (iConfig.getParameter<double>("massParticle1")),
     massParticle2_ (iConfig.getParameter<double>("massParticle2")),
     chargeOpt_ (iConfig.getParameter<int>("ChargeOpt")),
     triggerTypeDaughters_(iConfig.getParameter<int>("triggerTypeDaughters"))
 
 {
     produces<VertexCollection>();
 }

HLTDisplacedtktkVtxProducer::~HLTDisplacedtktkVtxProducer ( )

Definition at line 55 of file HLTDisplacedtktkVtxProducer.cc.

 {
 
 }

Member Function Documentation

void HLTDisplacedtktkVtxProducer::beginJob ( void )

virtual

Definition at line 79 of file HLTDisplacedtktkVtxProducer.cc.

 {
 
 }

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

private

Definition at line 200 of file HLTDisplacedtktkVtxProducer.cc.

References i, and convertSQLiteXML::ok.

Referenced by produce().

                                                                                                                       {
   bool ok=false;
   for (unsigned int i=0; i<refVect.size(); i++) {
     if ( refVect[i]->get<TrackRef>() == trackref ) {
       ok=true;
       break;
     }
   }
   return ok;
 }

void HLTDisplacedtktkVtxProducer::endJob ( void )

virtual

Definition at line 85 of file HLTDisplacedtktkVtxProducer.cc.

 {
     
 }

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

static

Definition at line 61 of file HLTDisplacedtktkVtxProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_25ns10e33_v2_cff::InputTag.

                                                                                         {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("Src",edm::InputTag("hltL3MuonCandidates"));
   desc.add<edm::InputTag>("PreviousCandTag",edm::InputTag(""));
   desc.add<double>("MaxEta",2.5);
   desc.add<double>("MinPt",0.0);
   desc.add<double>("MinPtPair",0.0);
   desc.add<double>("MinInvMass",1.0);
   desc.add<double>("MaxInvMass",20.0);
   desc.add<double>("massParticle1",0.1396);
   desc.add<double>("massParticle2",0.4937);
   desc.add<int>("ChargeOpt",-1);
   desc.add<int>("triggerTypeDaughters",0);
 
   descriptions.add("hltDisplacedtktkVtxProducer", desc);
 }

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

virtual

Implements edm::stream::EDProducerBase.

Definition at line 91 of file HLTDisplacedtktkVtxProducer.cc.

References funct::abs(), chargeOpt_, checkPreviousCand(), reco::e1, reco::e2, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), TransientVertex::isValid(), LogDebug, massParticle1_, massParticle2_, maxEta_, maxInvMass_, minInvMass_, minPt_, minPtPair_, AlCaHLTBitMon_ParallelJobs::p, p1, p2, previousCandToken_, edm::Event::put(), mathSSE::sqrt(), srcToken_, triggerTypeDaughters_, KalmanVertexFitter::vertex(), and GoodVertex_cfg::vertexCollection.

 {
     double const firstTrackMass = massParticle1_;
     double const firstTrackMass2 = firstTrackMass*firstTrackMass;
     double const secondTrackMass = massParticle2_;
     double const secondTrackMass2 = secondTrackMass*secondTrackMass;
 
     // get hold of track trks
     Handle<RecoChargedCandidateCollection> trackcands;
     iEvent.getByToken(srcToken_,trackcands);
     
     //get the transient track builder:
     edm::ESHandle<TransientTrackBuilder> theB;
     iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",theB);
 
     std::auto_ptr<VertexCollection> vertexCollection(new VertexCollection());
 
     // look at all trackcands,  check cuts and make vertices
     double e1,e2;
     Particle::LorentzVector p,p1,p2;
     
     RecoChargedCandidateCollection::const_iterator cand1;
     RecoChargedCandidateCollection::const_iterator cand2;
     
     // get the objects passing the previous filter
     Handle<TriggerFilterObjectWithRefs> previousCands;
     iEvent.getByToken(previousCandToken_,previousCands);
 
     vector<RecoChargedCandidateRef> vPrevCands;
     previousCands->getObjects(triggerTypeDaughters_,vPrevCands);
     
     for (cand1=trackcands->begin(); cand1!=trackcands->end(); cand1++) {
            TrackRef tk1 = cand1->get<TrackRef>();
            LogDebug("HLTDisplacedtktkVtxProducer") << " 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( ! checkPreviousCand( tk1, vPrevCands) ) continue;
     
            // cuts
            if (abs(cand1->eta())>maxEta_) continue;
            if (cand1->pt() < minPt_) continue;
         
         
         cand2=trackcands->begin();
         if(massParticle1_==massParticle2_){cand2 = cand1+1;}
 
         for (; cand2!=trackcands->end(); cand2++) {
              
              TrackRef tk2 = cand2->get<TrackRef>();
              if(tk1==tk2) continue;
          
              // eta cut
              LogDebug("HLTDisplacedtktkVtxProducer") << " 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( ! checkPreviousCand( tk2, vPrevCands) ) continue;
              
              // cuts
              if (abs(cand2->eta())>maxEta_) continue;
              if (cand2->pt() < minPt_) continue;
              
              // opposite sign or same sign
              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);
              p1 = Particle::LorentzVector(cand1->px(),cand1->py(),cand1->pz(),e1);
              p2 = Particle::LorentzVector(cand2->px(),cand2->py(),cand2->pz(),e2);
              p = p1+p2;
              
              
              if (p.pt()<minPtPair_) continue;
              
              double invmass = abs(p.mass());
              LogDebug("HLTDisplacedtktkVtxProducer") << " ... 1-2 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);
              t_tks.push_back(ttkp1);
              t_tks.push_back(ttkp2);
              
         
              if (t_tks.size()!=2) 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(vertexCollection);
 }