#include <ImpactParameterAlgorithm.h>

Public Member Functions
	ImpactParameterAlgorithm ()

	ImpactParameterAlgorithm (const edm::ParameterSet &parameters)

void	setPrimaryVertex (reco::Vertex *pv)

void	setTransientTrackBuilder (const TransientTrackBuilder *)

std::pair< float, reco::TauImpactParameterInfo >	tag (const reco::IsolatedTauTagInfoRef &, const reco::Vertex &)

	~ImpactParameterAlgorithm ()

Private Attributes
double	ip_max

double	ip_min

reco::Vertex *	primaryVertex

double	sip_min

const TransientTrackBuilder *	transientTrackBuilder

bool	use3D

bool	use_sign

Detailed Description

Definition at line 19 of file ImpactParameterAlgorithm.h.

Constructor & Destructor Documentation

◆ ImpactParameterAlgorithm() [1/2]

ImpactParameterAlgorithm::ImpactParameterAlgorithm ( const edm::ParameterSet & parameters )

Definition at line 17 of file ImpactParameterAlgorithm.cc.

                                                                                     {
   ip_min = parameters.getParameter<double>("TauImpactParameterMin");
   ip_max = parameters.getParameter<double>("TauImpactParameterMax");
   sip_min = parameters.getParameter<double>("TauImpactParameterSignificanceMin");
   use_sign = parameters.getParameter<bool>("UseTauImpactParameterSign");
   use3D = parameters.getParameter<bool>("UseTau3DImpactParameter");
 }

References ip_max, ip_min, sip_min, use3D, and use_sign.

◆ ImpactParameterAlgorithm() [2/2]

ImpactParameterAlgorithm::ImpactParameterAlgorithm ( )

Definition at line 9 of file ImpactParameterAlgorithm.cc.

                                                    {
   ip_min = -9999;
   ip_max = 9999;
   sip_min = 0;
   use_sign = false;
   use3D = false;
 }

References ip_max, ip_min, sip_min, use3D, and use_sign.

◆ ~ImpactParameterAlgorithm()

ImpactParameterAlgorithm::~ImpactParameterAlgorithm ( )

inline

Definition at line 28 of file ImpactParameterAlgorithm.h.

28 {}

Member Function Documentation

◆ setPrimaryVertex()

void ImpactParameterAlgorithm::setPrimaryVertex ( reco::Vertex * pv )

inline

Definition at line 30 of file ImpactParameterAlgorithm.h.

30 { primaryVertex = pv; }

References primaryVertex, and MetAnalyzer::pv().

◆ setTransientTrackBuilder()

void ImpactParameterAlgorithm::setTransientTrackBuilder ( const TransientTrackBuilder * builder )

Definition at line 25 of file ImpactParameterAlgorithm.cc.

                                                                                             {
   transientTrackBuilder = builder;
 }

References transientTrackBuilder.

Referenced by ImpactParameter::produce().

◆ tag()

std::pair< float, reco::TauImpactParameterInfo > ImpactParameterAlgorithm::tag	(	const reco::IsolatedTauTagInfoRef &	tauRef,
		const reco::Vertex &	pv
	)

Definition at line 29 of file ImpactParameterAlgorithm.cc.

                                                                                                    {
   if (transientTrackBuilder == nullptr) {
     throw cms::Exception("NullTransientTrackBuilder") << "Transient track builder is 0. ";
   }
  
   reco::TauImpactParameterInfo resultExtended;
   resultExtended.setIsolatedTauTag(tauRef);
  
   const reco::Jet* jet = tauRef->jet().get();
   GlobalVector direction(jet->px(), jet->py(), jet->pz());
  
   const reco::TrackRefVector& tracks = tauRef->selectedTracks();
  
   edm::RefVector<reco::TrackCollection>::const_iterator iTrack;
   for (iTrack = tracks.begin(); iTrack != tracks.end(); iTrack++) {
     const reco::TransientTrack transientTrack = (transientTrackBuilder->build(&(**iTrack)));
  
     SignedTransverseImpactParameter stip;
     Measurement1D ip = stip.apply(transientTrack, direction, pv).second;
  
     SignedImpactParameter3D signed_ip3D;
     Measurement1D ip3D = signed_ip3D.apply(transientTrack, direction, pv).second;
     LogDebug("ImpactParameterAlgorithm::tag") << "check pv,ip3d " << pv.z() << " " << ip3D.value();
     if (!use_sign) {
       Measurement1D tmp2D(fabs(ip.value()), ip.error());
       ip = tmp2D;
  
       Measurement1D tmp3D(fabs(ip3D.value()), ip3D.error());
       ip3D = tmp3D;
     }
  
     reco::TauImpactParameterTrackData theData;
  
     theData.transverseIp = ip;
     theData.ip3D = ip3D;
     resultExtended.storeTrackData(*iTrack, theData);
   }
  
   float discriminator = resultExtended.discriminator(ip_min, ip_max, sip_min, use_sign, use3D);
  
   return std::make_pair(discriminator, resultExtended);
 }