#include <TwoTrackVrtBDTSelector.h>

Inheritance diagram for Rec::TwoTrackVrtBDTSelector:

Collaboration diagram for Rec::TwoTrackVrtBDTSelector:

Public Member Functions
	TwoTrackVrtBDTSelector (const std::string &type, const std::string &name, const IInterface *parent)

virtual	~TwoTrackVrtBDTSelector ()

StatusCode	initialize ()

StatusCode	finalize ()

bool	isgood (const std::pair< const xAOD::TrackParticle , const xAOD::TrackParticle > tracks, const xAOD::Vertex &candV, std::pair< ROOT::Math::XYZTVector, ROOT::Math::XYZTVector > moms, const xAOD::Vertex &tPV) const final

bool	isgood (const std::pair< const xAOD::TrackParticle , const xAOD::TrackParticle > tracks, const xAOD::Vertex &candV, std::pair< ROOT::Math::XYZTVector, ROOT::Math::XYZTVector > moms, const xAOD::Vertex &tPV, float &quality) const final

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Static Private Member Functions
static double	vrtRadiusError (const Amg::Vector3D &secVrt, const std::vector< float > &vrtErr)

static int	getIBLHit (const xAOD::TrackParticle *Part)

static int	getBLHit (const xAOD::TrackParticle *Part)

Private Attributes
Gaudi::Property< float >	m_vrt2TrMassLimit {this,"Vrt2TrMassLimit",4000., "Maximal allowed mass for 2-track vertices" }

Gaudi::Property< float >	m_vrt2TrPtMin {this, "Vrt2TrPtMin", 1000., "Minimal allowed Pt for 2-track vertices." }

Gaudi::Property< float >	m_vrt2TrPtMax {this, "Vrt2TrPtMax", 5.e5, "Maximal allowed Pt for 2-track vertices. Calibration limit" }

Gaudi::Property< float >	m_sel2VrtProbCut {this, "Sel2VrtProbCut", 0.02, "Cut on probability of 2-track vertex for initial selection" }

Gaudi::Property< float >	m_maxSVRadiusCut {this, "MaxSVRadiusCut", 140., "Cut on maximal radius of SV (def = Pixel detector size)" }

Gaudi::Property< float >	m_cosSVPVCut {this, "cosSVPVCut", 0., "Cut on cos of angle between SV-PV and full vertex momentum" }

Gaudi::Property< bool >	m_do2TrkIBLChecks {this, "do2TrkIBLChecks", true, "IBL and B-layer hit requrirements based on the position of 2-track DV." }

Gaudi::Property< bool >	m_useVertexCleaning {this, "useVertexCleaning", true, "Clean vertices by requiring pixel hit presence according to vertex position" }

Gaudi::Property< float >	m_firstPixelLayerR {this, "FirstPixelLayerR", 32.0,"Radius of the first Pixel layer" }

Gaudi::Property< float >	m_v2tBDTCut {this, "v2tBDTCut", -0.5, "BDT cut to select 2-track vertices" }

Gaudi::Property< std::string >	m_calibFileName {this, "CalibFileName", "Fake2TrVertexReject.MVA.v02.root", " MVA calibration file for 2-track fake vertices removal" }

std::unique_ptr< MVAUtils::BDT >	m_SV2T_BDT

ToolHandle< Trk::TrkVKalVrtFitter >	m_fitSvc {this, "VertexFitterTool", "Trk::TrkVKalVrtFitter/VertexFitterTool", "Vertex Fitter tool for 2-track selector"}

double	m_massPi {}

double	m_massP {}

double	m_massE {}

std::string	m_instanceName {}

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 36 of file TwoTrackVrtBDTSelector.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TwoTrackVrtBDTSelector()

Rec::TwoTrackVrtBDTSelector::TwoTrackVrtBDTSelector	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 26 of file TwoTrackVrtBDTSelector.cxx.

                                                                     :
     AthAlgTool(type,name,parent),
     m_SV2T_BDT(nullptr),
     m_instanceName(name)
    {
 //
 // Declare additional interface
 //
     declareInterface< ITwoTrackVertexSelector >(this);
 //
     m_massPi  =  Trk::ParticleMasses::mass[Trk::pion];
     m_massP   =  Trk::ParticleMasses::mass[Trk::proton];
     m_massE   =  Trk::ParticleMasses::mass[Trk::electron];
   }

◆ ~TwoTrackVrtBDTSelector()

Rec::TwoTrackVrtBDTSelector::~TwoTrackVrtBDTSelector ( )

virtual

Definition at line 44 of file TwoTrackVrtBDTSelector.cxx.

                                                  {
      ATH_MSG_DEBUG("TwoTrackVertBDTSelector destructor called");
   }

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                      {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ finalize()

StatusCode Rec::TwoTrackVrtBDTSelector::finalize ( )

Definition at line 70 of file TwoTrackVrtBDTSelector.cxx.

   {
     ATH_MSG_DEBUG("TwoTrackVrtBDTSelector finalize()");
     return StatusCode::SUCCESS; 
   }

◆ getBLHit()

int Rec::TwoTrackVrtBDTSelector::getBLHit ( const xAOD::TrackParticle * Part )

staticprivate

Definition at line 183 of file TwoTrackVrtBDTSelector.cxx.

   {
         uint8_t BLhit,BLexp;
         if(!Part->summaryValue( BLexp,  xAOD::expectNextToInnermostPixelLayerHit) )     BLexp = 0;
         if( BLexp==0 ) return -1;
         if(!Part->summaryValue( BLhit,  xAOD::numberOfNextToInnermostPixelLayerHits) )  BLhit = 0;
         if(BLhit) return 1;
         else      return 0;
   }

◆ getIBLHit()

int Rec::TwoTrackVrtBDTSelector::getIBLHit ( const xAOD::TrackParticle * Part )

staticprivate

Definition at line 174 of file TwoTrackVrtBDTSelector.cxx.

   {
         uint8_t IBLhit,IBLexp;
         if(!Part->summaryValue( IBLexp,  xAOD::expectInnermostPixelLayerHit) )    IBLexp = 0;
         if( IBLexp==0 ) return -1;
         if(!Part->summaryValue( IBLhit,  xAOD::numberOfInnermostPixelLayerHits) ) IBLhit = 0;
         if(IBLhit) return 1;
         else       return 0;
   }

◆ initialize()

StatusCode Rec::TwoTrackVrtBDTSelector::initialize ( )

Definition at line 49 of file TwoTrackVrtBDTSelector.cxx.

                                                {
     ATH_MSG_DEBUG( "Initialising TwoTrackVrtBDTSelector" );
  
     ATH_CHECK( m_fitSvc.retrieve() );
 //--------------------------------------------------------
     //std::string fileName="NewVrtSecInclusiveTool/Fake2TrVertexReject.MVA.v01.root";   ///For local calibration file
     //std::string rootFilePath = PathResolver::find_file(fileName, "DATAPATH");         ///
     std::string rootFilePath = PathResolver::find_calib_file("NewVrtSecInclusiveTool/"+m_calibFileName);
     std::unique_ptr<TFile> rootFile(TFile::Open(rootFilePath.c_str(), "READ")); 
     if (!rootFile) {
         ATH_MSG_FATAL("Could not retrieve root file: " << m_calibFileName);
         return StatusCode::FAILURE;
     }
     std::unique_ptr<TTree> training((TTree*)rootFile->Get("BDT"));
     m_SV2T_BDT = std::make_unique<MVAUtils::BDT>(training.get());
 //--------------------------------------------------------
     return StatusCode::SUCCESS;
   }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& Rec::ITwoTrackVertexSelector::interfaceID ( )

inlinestaticinherited

Definition at line 35 of file ITwoTrackVertexSelector.h.

35 { return IID_ITwoTrackVertexSelector;}

◆ isgood() [1/2]

bool Rec::TwoTrackVrtBDTSelector::isgood	(	const std::pair< const xAOD::TrackParticle , const xAOD::TrackParticle >	tracks,
		const xAOD::Vertex &	candV,
		std::pair< ROOT::Math::XYZTVector, ROOT::Math::XYZTVector >	moms,
		const xAOD::Vertex &	tPV
	)		const

finalvirtual

Implements Rec::ITwoTrackVertexSelector.

Definition at line 77 of file TwoTrackVrtBDTSelector.cxx.

   {
       float quality;
       return isgood(iTrks, candV, moms,tPV,quality);
   }

◆ isgood() [2/2]

bool Rec::TwoTrackVrtBDTSelector::isgood	(	const std::pair< const xAOD::TrackParticle , const xAOD::TrackParticle >	tracks,
		const xAOD::Vertex &	candV,
		std::pair< ROOT::Math::XYZTVector, ROOT::Math::XYZTVector >	moms,
		const xAOD::Vertex &	tPV,
		float &	quality
	)		const

finalvirtual

Implements Rec::ITwoTrackVertexSelector.

Definition at line 86 of file TwoTrackVrtBDTSelector.cxx.

   {
       quality=-99.;  
       double Prob2v=TMath::Prob(candV.chiSquared(),1);
       if( Prob2v       < m_sel2VrtProbCut )   return false;    // Vertex probability check
       double vrtR=candV.position().perp();
       auto sumMom=moms.first+moms.second;
       if( sumMom.M()   > m_vrt2TrMassLimit ) return false;    // Invariant mass check 
       if( vrtR         > m_maxSVRadiusCut)   return false;    // Too far from interaction point
       if( sumMom.Pt()  < m_vrt2TrPtMin)      return false;    // Summary momentum > minimal allowed 2-track vertex pt
       if( sumMom.Pt()  > m_vrt2TrPtMax)      return false;    // Maximal allowed 2-track vertex Pt
  
       ROOT::Math::XYZVector vSVPV(candV.x()-tPV.x(),candV.y()-tPV.y(),candV.z()-tPV.z());
       double cosSVPV=vSVPV.Unit().Dot(sumMom.Vect().Unit());
       if(cosSVPV  < m_cosSVPVCut)             return false;    // Angle between tV-PV direction and summary momentum
          
  
       double vrtRErr=vrtRadiusError(candV.position(),candV.covariance() );
 //
 // Check pixel hits vs vertex positions.
       int ihitIBL  = getIBLHit(iTrks.first);
       int jhitIBL  = getIBLHit(iTrks.second);
       if( m_do2TrkIBLChecks && ( (ihitIBL==0&&jhitIBL>0) || (ihitIBL>0&&jhitIBL==0) ) ) return false;
       int ihitBL   = getBLHit (iTrks.first);
       int jhitBL   = getBLHit (iTrks.second);
 //--Very general cleaning cuts based on ID geometry and applicable to all processes
       if( m_do2TrkIBLChecks && vrtR<m_firstPixelLayerR-2.*vrtRErr ){
         if( ihitIBL<1 && ihitBL<1) return false;
         if( jhitIBL<1 && jhitBL<1) return false;
       }
       float ihitR  = iTrks.first->radiusOfFirstHit();
       float jhitR  = iTrks.second->radiusOfFirstHit();
       if(std::abs(ihitR-jhitR)>50.)return false;             //- FMPs are in very different layers
       if( vrtR-std::min(ihitR,jhitR) > 50.) return false;    //- FMP is closer to (0,0) than SV itself
       if(ihitR-vrtR > 180.+2.*vrtRErr)return false;  //- Distance FMP-vertex should be less then SCT-Pixel gap
       if(jhitR-vrtR > 180.+2.*vrtRErr)return false;  //- Distance FMP-vertex should be less then SCT-Pixel gap
 //-------------------------------------------------------
       if(m_useVertexCleaning){ //More agressive cleaning 
         if(std::abs(ihitR-jhitR)>12.) return false;
         if( ihitR-vrtR > 36.) return false; // Too big dR between vertex and hit in pixel
         if( jhitR-vrtR > 36.) return false; // Should be another layer in between 
         if( ihitR-vrtR <-2.*vrtRErr) return false; // Vertex is behind hit in pixel 
         if( jhitR-vrtR <-2.*vrtRErr) return false; // Vertex is behind hit in pixel 
       }
 //-------------------BDT based rejection
       std::vector<double> impact,impactError;
       m_fitSvc->VKalGetImpact( iTrks.first, tPV.position(), 1, impact, impactError);
       float trk1Signif = sqrt( impact[0]*impact[0]/impactError[0] + impact[1]*impact[1]/impactError[2]);
       m_fitSvc->VKalGetImpact( iTrks.second, tPV.position(), 1, impact, impactError);
       float trk2Signif = sqrt( impact[0]*impact[0]/impactError[0] + impact[1]*impact[1]/impactError[2]);
       float minPtT = std::min(iTrks.first->pt(),iTrks.second->pt());
       std::vector<float> VARS(10);
       VARS[0]=Prob2v;
       VARS[1]=log(sumMom.Pt());
       VARS[2]=log(minPtT);
       VARS[3]=log(vrtR<20. ? vSVPV.Rho() : vrtR);
       VARS[4]=log(std::min(trk1Signif,trk2Signif));
       VARS[5]=log(std::max(trk1Signif,trk2Signif));
       VARS[6]=sumMom.M();
       VARS[7]=sqrt(std::abs(1.-cosSVPV*cosSVPV));
       VARS[8]=vSVPV.Eta();
       VARS[9]=std::max(ihitR,jhitR);
       quality=m_SV2T_BDT->GetGradBoostMVA(VARS);
       if(quality<m_v2tBDTCut) return false;         // BDT rejection
  
       return true;   
   }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

◆ vrtRadiusError()

double Rec::TwoTrackVrtBDTSelector::vrtRadiusError	(	const Amg::Vector3D &	secVrt,
		const std::vector< float > &	vrtErr
	)

staticprivate

Definition at line 161 of file TwoTrackVrtBDTSelector.cxx.

   {
     double DirX=SecVrt.x(), DirY=SecVrt.y(); 
     double Covar =    DirX*VrtErr[0]*DirX
                   +2.*DirX*VrtErr[1]*DirY
                      +DirY*VrtErr[2]*DirY;
     Covar /= DirX*DirX + DirY*DirY;
     Covar=std::sqrt(std::abs(Covar));
     if(Covar != Covar)  Covar = 0.;
     return Covar;
   }

Member Data Documentation

◆ m_calibFileName

Gaudi::Property<std::string> Rec::TwoTrackVrtBDTSelector::m_calibFileName {this, "CalibFileName", "Fake2TrVertexReject.MVA.v02.root", " MVA calibration file for 2-track fake vertices removal" }

private

Definition at line 78 of file TwoTrackVrtBDTSelector.h.

◆ m_cosSVPVCut

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_cosSVPVCut {this, "cosSVPVCut", 0., "Cut on cos of angle between SV-PV and full vertex momentum" }

private

Definition at line 69 of file TwoTrackVrtBDTSelector.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_do2TrkIBLChecks

Gaudi::Property<bool> Rec::TwoTrackVrtBDTSelector::m_do2TrkIBLChecks {this, "do2TrkIBLChecks", true, "IBL and B-layer hit requrirements based on the position of 2-track DV." }

private

Definition at line 72 of file TwoTrackVrtBDTSelector.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_firstPixelLayerR

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_firstPixelLayerR {this, "FirstPixelLayerR", 32.0,"Radius of the first Pixel layer" }

private

Definition at line 74 of file TwoTrackVrtBDTSelector.h.

◆ m_fitSvc

ToolHandle<Trk::TrkVKalVrtFitter> Rec::TwoTrackVrtBDTSelector::m_fitSvc {this, "VertexFitterTool", "Trk::TrkVKalVrtFitter/VertexFitterTool", "Vertex Fitter tool for 2-track selector"}

private

Definition at line 82 of file TwoTrackVrtBDTSelector.h.

◆ m_instanceName

std::string Rec::TwoTrackVrtBDTSelector::m_instanceName {}

private

Definition at line 91 of file TwoTrackVrtBDTSelector.h.

◆ m_massE

double Rec::TwoTrackVrtBDTSelector::m_massE {}

private

Definition at line 90 of file TwoTrackVrtBDTSelector.h.

◆ m_massP

double Rec::TwoTrackVrtBDTSelector::m_massP {}

private

Definition at line 89 of file TwoTrackVrtBDTSelector.h.

◆ m_massPi

double Rec::TwoTrackVrtBDTSelector::m_massPi {}

private

Definition at line 88 of file TwoTrackVrtBDTSelector.h.

◆ m_maxSVRadiusCut

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_maxSVRadiusCut {this, "MaxSVRadiusCut", 140., "Cut on maximal radius of SV (def = Pixel detector size)" }

private

Definition at line 68 of file TwoTrackVrtBDTSelector.h.

◆ m_sel2VrtProbCut

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_sel2VrtProbCut {this, "Sel2VrtProbCut", 0.02, "Cut on probability of 2-track vertex for initial selection" }

private

Definition at line 67 of file TwoTrackVrtBDTSelector.h.

◆ m_SV2T_BDT

std::unique_ptr<MVAUtils::BDT> Rec::TwoTrackVrtBDTSelector::m_SV2T_BDT

private

Definition at line 80 of file TwoTrackVrtBDTSelector.h.

◆ m_useVertexCleaning

Gaudi::Property<bool> Rec::TwoTrackVrtBDTSelector::m_useVertexCleaning {this, "useVertexCleaning", true, "Clean vertices by requiring pixel hit presence according to vertex position" }

private

Definition at line 73 of file TwoTrackVrtBDTSelector.h.

◆ m_v2tBDTCut

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_v2tBDTCut {this, "v2tBDTCut", -0.5, "BDT cut to select 2-track vertices" }

private

Definition at line 77 of file TwoTrackVrtBDTSelector.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_vrt2TrMassLimit

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_vrt2TrMassLimit {this,"Vrt2TrMassLimit",4000., "Maximal allowed mass for 2-track vertices" }

private

Definition at line 64 of file TwoTrackVrtBDTSelector.h.

◆ m_vrt2TrPtMax

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_vrt2TrPtMax {this, "Vrt2TrPtMax", 5.e5, "Maximal allowed Pt for 2-track vertices. Calibration limit" }

private

Definition at line 66 of file TwoTrackVrtBDTSelector.h.

◆ m_vrt2TrPtMin

Gaudi::Property<float> Rec::TwoTrackVrtBDTSelector::m_vrt2TrPtMin {this, "Vrt2TrPtMin", 1000., "Minimal allowed Pt for 2-track vertices." }

private

Definition at line 65 of file TwoTrackVrtBDTSelector.h.

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ TwoTrackVrtBDTSelector()

◆ ~TwoTrackVrtBDTSelector()

Member Function Documentation

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ finalize()

◆ getBLHit()

◆ getIBLHit()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ isgood() [1/2]

◆ isgood() [2/2]

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

◆ vrtRadiusError()

Member Data Documentation

◆ m_calibFileName

◆ m_cosSVPVCut

◆ m_detStore

◆ m_do2TrkIBLChecks

◆ m_evtStore

◆ m_firstPixelLayerR

◆ m_fitSvc

◆ m_instanceName

◆ m_massE

◆ m_massP

◆ m_massPi

◆ m_maxSVRadiusCut

◆ m_sel2VrtProbCut

◆ m_SV2T_BDT

◆ m_useVertexCleaning

◆ m_v2tBDTCut

◆ m_varHandleArraysDeclared

◆ m_vhka

◆ m_vrt2TrMassLimit

◆ m_vrt2TrPtMax

◆ m_vrt2TrPtMin