#include <TemplatedSimpleSecondaryVertexComputer.h>

Inheritance diagram for TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >:

Public Types
typedef reco::TemplatedSecondaryVertexTagInfo< IPTI, VTX >	TagInfo

Public Member Functions
float	discriminator (const TagInfoHelper &tagInfos) const

	TemplatedSimpleSecondaryVertexComputer (const edm::ParameterSet &parameters)

Public Member Functions inherited from JetTagComputer
const std::vector< std::string > &	getInputLabels () const

virtual void	initialize (const JetTagComputerRecord &)

	JetTagComputer ()

	JetTagComputer (const edm::ParameterSet &configuration)

float	operator() (const reco::BaseTagInfo &info) const

float	operator() (const TagInfoHelper &helper) const

void	setupDone ()

virtual	~JetTagComputer ()

Private Attributes
unsigned int	minTracks

unsigned int	minVertices_

bool	unBoost

bool	use2d

bool	useSig

Additional Inherited Members
Protected Member Functions inherited from JetTagComputer
virtual float	discriminator (const reco::BaseTagInfo &) const

void	uses (unsigned int id, const std::string &label)

void	uses (const std::string &label)

Detailed Description

template<class IPTI, class VTX>
class TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >

Definition at line 17 of file TemplatedSimpleSecondaryVertexComputer.h.

Member Typedef Documentation

template<class IPTI , class VTX >

typedef reco::TemplatedSecondaryVertexTagInfo<IPTI,VTX> TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::TagInfo

Definition at line 19 of file TemplatedSimpleSecondaryVertexComputer.h.

Constructor & Destructor Documentation

template<class IPTI , class VTX >

TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::TemplatedSimpleSecondaryVertexComputer ( const edm::ParameterSet & parameters )

inline

Definition at line 21 of file TemplatedSimpleSecondaryVertexComputer.h.

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::minVertices_, and JetTagComputer::uses().

                                                                               :
             use2d(!parameters.getParameter<bool>("use3d")),
         useSig(parameters.getParameter<bool>("useSignificance")),
             unBoost(parameters.getParameter<bool>("unBoost")),
             minTracks(parameters.getParameter<unsigned int>("minTracks")),
             minVertices_(1)
           { 
             uses("svTagInfos"); 
             minVertices_    = parameters.existsAs<unsigned int>("minVertices") ?  parameters.getParameter<unsigned int>("minVertices") : 1 ;
           }

Member Function Documentation

template<class IPTI , class VTX >

float TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::discriminator ( const TagInfoHelper & tagInfos ) const

inlinevirtual

Reimplemented from JetTagComputer.

Definition at line 32 of file TemplatedSimpleSecondaryVertexComputer.h.

References reco::TemplatedSecondaryVertexTagInfo< IPTI, VTX >::flightDistance(), CustomPhysics_cfi::gamma, JetTagComputer::TagInfoHelper::get(), training_settings::idx, info(), cmsBatch::log, TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::minTracks, TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::minVertices_, reco::TemplatedSecondaryVertexTagInfo< IPTI, VTX >::nVertexTracks(), reco::TemplatedSecondaryVertexTagInfo< IPTI, VTX >::nVertices(), reco::TemplatedSecondaryVertexTagInfo< IPTI, VTX >::secondaryVertex(), Measurement1D::significance(), TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::unBoost, TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::use2d, TemplatedSimpleSecondaryVertexComputer< IPTI, VTX >::useSig, Measurement1D::value(), relativeConstraints::value, and reco::TrackKinematics::vectorSum().

     {
         const TagInfo &info =
                 tagInfos.get<TagInfo>();
         if(info.nVertices() < minVertices_) return -1;
                 unsigned int idx = 0;
         while(idx < info.nVertices()) {
             if (info.nVertexTracks(idx) >= minTracks)
                 break;
             idx++;
         }
         if (idx >= info.nVertices())
             return -1.0;
 
         double gamma;
         if (unBoost) {
             reco::TrackKinematics kinematics(
                         info.secondaryVertex(idx));
             gamma = kinematics.vectorSum().Gamma();
         } else
             gamma = 1.0;
 
         double value;
         if (useSig)
             value = info.flightDistance(idx, use2d).significance();
         else
             value = info.flightDistance(idx, use2d).value();
 
         value /= gamma;
 
         if (useSig)
             value = (value > 0) ? +std::log(1 + value)
                                 : -std::log(1 - value);
 
         return value;
     }