Inheritance diagram for PFRecoTauDiscriminationByFlightPathSignificance:

Public Member Functions
void	beginEvent (const edm::Event &, const edm::EventSetup &) override

double	discriminate (const reco::PFTauRef &) const override

	PFRecoTauDiscriminationByFlightPathSignificance (const ParameterSet &iConfig)

	~PFRecoTauDiscriminationByFlightPathSignificance () override

Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
virtual double	discriminate (const TauRef &tau) const =0

virtual void	endEvent (edm::Event &)

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

	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)

	TauDiscriminationProducerBase ()

	~TauDiscriminationProducerBase () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Member Functions
double	threeProngFlightPathSig (const PFTauRef &) const

double	vertexSignificance (reco::Vertex const &, reco::Vertex const &, GlobalVector const &) const

Private Attributes
bool	booleanOutput

double	flightPathSig

const TransientTrackBuilder *	transientTrackBuilder

reco::tau::RecoTauVertexAssociator *	vertexAssociator_

bool	withPVError

Additional Inherited Members
Public Types inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
typedef std::vector< TauType >	TauCollection

typedef edm::Ref< TauCollection >	TauRef

typedef edm::RefProd< TauCollection >	TauRefProd

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

Protected Attributes inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
std::string	moduleLabel_

double	prediscriminantFailValue_

edm::EDGetTokenT< TauCollection >	Tau_token

size_t	tauIndex_

edm::InputTag	TauProducer_

Detailed Description

Definition at line 24 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

Constructor & Destructor Documentation

PFRecoTauDiscriminationByFlightPathSignificance::PFRecoTauDiscriminationByFlightPathSignificance ( const ParameterSet & iConfig )

inlineexplicit

Definition at line 27 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References CaloRecoTauDiscriminationByFlightPathSignificance_cfi::flightPathSig, and edm::ParameterSet::getParameter().

       :PFTauDiscriminationProducerBase(iConfig){
       flightPathSig     = iConfig.getParameter<double>("flightPathSig");
       withPVError       = iConfig.getParameter<bool>("UsePVerror");
       booleanOutput         = iConfig.getParameter<bool>("BooleanOutput");
       //      edm::ConsumesCollector iC(consumesCollector());
       vertexAssociator_ = new reco::tau::RecoTauVertexAssociator(iConfig.getParameter<ParameterSet>("qualityCuts"),consumesCollector());
     }

PFRecoTauDiscriminationByFlightPathSignificance::~PFRecoTauDiscriminationByFlightPathSignificance ( )

inlineoverride

Definition at line 36 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

36 {}

Member Function Documentation

void PFRecoTauDiscriminationByFlightPathSignificance::beginEvent	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 54 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References edm::EventSetup::get(), and edm::ESHandle< T >::product().

                                                   {
 
    vertexAssociator_->setEvent(iEvent);
 
   // Transient Tracks
   edm::ESHandle<TransientTrackBuilder> builder;
   iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",builder);
   transientTrackBuilder = builder.product();
 
 }

double PFRecoTauDiscriminationByFlightPathSignificance::discriminate ( const reco::PFTauRef & tau ) const

override

Definition at line 66 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References CaloRecoTauDiscriminationByFlightPathSignificance_cfi::flightPathSig.

                                                                                              {
 
   if(booleanOutput) return ( threeProngFlightPathSig(tau) > flightPathSig ? 1. : 0. );
   return threeProngFlightPathSig(tau);
 }

double PFRecoTauDiscriminationByFlightPathSignificance::threeProngFlightPathSig ( const PFTauRef & tau ) const

private

Definition at line 72 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References reco::PFCandidate::gsfTrackRef(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::isNull(), TransientVertex::isValid(), impactParameterTagInfos_cfi::primaryVertex, reco::PFCandidate::trackRef(), and KalmanVertexFitter::vertex().

                                {
   double flightPathSignificance = 0;
 
   reco::VertexRef primaryVertex = vertexAssociator_->associatedVertex(*tau);
 
   if (primaryVertex.isNull()) {
     edm::LogError("FlightPathSignficance") << "Could not get vertex associated"
       << " to tau, returning -999!" << std::endl;
     return -999;
   }
 
   //Secondary vertex
   const vector<PFCandidatePtr>& pfSignalCandidates = tau->signalPFChargedHadrCands();
   vector<TransientTrack> transientTracks;
   vector<PFCandidatePtr>::const_iterator iTrack;
   for(iTrack = pfSignalCandidates.begin(); iTrack!= pfSignalCandidates.end(); iTrack++){
     const PFCandidate& pfCand = *(iTrack->get());
     if(pfCand.trackRef().isNonnull()){
       const TransientTrack transientTrack = transientTrackBuilder->build(pfCand.trackRef());
       transientTracks.push_back(transientTrack);
     }
     else if(pfCand.gsfTrackRef().isNonnull()){
       const TransientTrack transientTrack = transientTrackBuilder->build(pfCand.gsfTrackRef());
       transientTracks.push_back(transientTrack);
     }
   }
   if(transientTracks.size() > 1){
     KalmanVertexFitter kvf(true);
     TransientVertex tv = kvf.vertex(transientTracks);
 
     if(tv.isValid()){
       GlobalVector tauDir(tau->px(),
           tau->py(),
           tau->pz());
       Vertex secVer = tv;
       // We have to un-const the PV for some reason
       reco::Vertex primaryVertexNonConst = *primaryVertex;
       flightPathSignificance = vertexSignificance(primaryVertexNonConst,secVer,tauDir);
     }
   }
   return flightPathSignificance;
 }

double PFRecoTauDiscriminationByFlightPathSignificance::vertexSignificance	(	reco::Vertex const &	pv,
		reco::Vertex const &	sv,
		GlobalVector const &	direction
	)		const

private

Definition at line 116 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References reco::TemplatedSecondaryVertex< SV >::computeDist3d(), DEFINE_FWK_MODULE, and Measurement1D::significance().

                                                                                    {
   return SecondaryVertex::computeDist3d(pv,sv,direction,withPVError).significance();
 }