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

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

Static Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
static void	fillProducerDescriptions (edm::ParameterSetDescription &desc)

static std::string	getTauTypeString ()
	helper method to retrieve tau type name, e.g. to build correct cfi getter More...

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 27 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

Constructor & Destructor Documentation

PFRecoTauDiscriminationByFlightPathSignificance::PFRecoTauDiscriminationByFlightPathSignificance ( const ParameterSet & iConfig )

inlineexplicit

Definition at line 30 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 39 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References fillDescriptions().

39 {}

Member Function Documentation

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

overridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 58 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 70 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References CaloRecoTauDiscriminationByFlightPathSignificance_cfi::flightPathSig.

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

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

static

Definition at line 123 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addOptional(), DEFINE_FWK_MODULE, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                             {
   // pfRecoTauDiscriminationByFlightPathSignificance
   edm::ParameterSetDescription desc;
   desc.add<double>("flightPathSig", 1.5);
   desc.add<edm::InputTag>("PVProducer", edm::InputTag("offlinePrimaryVertices"));
   desc.add<bool>("BooleanOutput", true);
   desc.add<edm::InputTag>("PFTauProducer", edm::InputTag("pfRecoTauProducer"));
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("BooleanOperator", "and");
     {
       edm::ParameterSetDescription psd1;
       psd1.add<double>("cut");
       psd1.add<edm::InputTag>("Producer");
       psd0.addOptional<edm::ParameterSetDescription>("leadTrack", psd1);
     }
     desc.add<edm::ParameterSetDescription>("Prediscriminants", psd0);
   }
   desc.add<bool>("UsePVerror", true);
   descriptions.add("pfRecoTauDiscriminationByFlightPathSignificance", desc);
 }

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

private

Definition at line 76 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

References edm::Ref< C, T, F >::isNull(), TransientVertex::isValid(), impactParameterTagInfos_cfi::primaryVertex, 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
   vector<TransientTrack> transientTracks;
   for(const auto& pfSignalCand : tau->signalPFChargedHadrCands()){
     if(pfSignalCand->trackRef().isNonnull()){
       const TransientTrack transientTrack = transientTrackBuilder->build(pfSignalCand->trackRef());
       transientTracks.push_back(transientTrack);
     }
     else if(pfSignalCand->gsfTrackRef().isNonnull()){
       const TransientTrack transientTrack = transientTrackBuilder->build(pfSignalCand->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 117 of file PFRecoTauDiscriminationByFlightPathSignificance.cc.

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

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