#include <SimpleSecondaryVertexComputer.h>

Inheritance diagram for SimpleSecondaryVertexComputer:

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

	SimpleSecondaryVertexComputer (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

Definition at line 16 of file SimpleSecondaryVertexComputer.h.

Constructor & Destructor Documentation

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

inline

Definition at line 18 of file SimpleSecondaryVertexComputer.h.

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), 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

float SimpleSecondaryVertexComputer::discriminator ( const TagInfoHelper & tagInfos ) const

inlinevirtual

Reimplemented from JetTagComputer.

Definition at line 29 of file SimpleSecondaryVertexComputer.h.

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

     {
         const reco::SecondaryVertexTagInfo &info =
                 tagInfos.get<reco::SecondaryVertexTagInfo>();
         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;
     }