InDet::TruthPixelClusterSplitProbTool Class Reference

#include <TruthPixelClusterSplitProbTool.h>

Inheritance diagram for InDet::TruthPixelClusterSplitProbTool:

Collaboration diagram for InDet::TruthPixelClusterSplitProbTool:

Public Member Functions
	TruthPixelClusterSplitProbTool (const std::string &t, const std::string &n, const IInterface *p)
virtual	~TruthPixelClusterSplitProbTool ()=default
virtual StatusCode	initialize () override
virtual InDet::PixelClusterSplitProb	splitProbability (const InDet::PixelCluster &origCluster) const override
virtual InDet::PixelClusterSplitProb	splitProbability (const InDet::PixelCluster &origCluster, const Trk::TrackParameters &trackParameters) const override

Private Member Functions
InDet::PixelClusterSplitProb	compileSplitProbability (std::vector< double > &vectorOfProbs) const

Private Attributes
ToolHandle< TruthClusterizationFactory >	m_truthClusterizationFactory {this, "NnClusterizationFactory", "InDet::NnClusterizationFactory/TruthClusterizationFactory"}
DoubleArrayProperty	m_priorMultiplicityContent {this, "PriorMultiplicityContent", {2793337, 82056, 19944}}

Detailed Description

Definition at line 32 of file TruthPixelClusterSplitProbTool.h.

Constructor & Destructor Documentation

TruthPixelClusterSplitProbTool()

InDet::TruthPixelClusterSplitProbTool::TruthPixelClusterSplitProbTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 27 of file TruthPixelClusterSplitProbTool.cxx.

    :base_class(t,n,p)
  {
    declareInterface<IPixelClusterSplitProbTool>(this);
  }

~TruthPixelClusterSplitProbTool()

virtual InDet::TruthPixelClusterSplitProbTool::~TruthPixelClusterSplitProbTool ( )

virtualdefault

Member Function Documentation

compileSplitProbability()

InDet::PixelClusterSplitProb InDet::TruthPixelClusterSplitProbTool::compileSplitProbability ( std::vector< double > &  vectorOfProbs ) const

private

Definition at line 100 of file TruthPixelClusterSplitProbTool.cxx.

  {
 
 
    double sum=0;
 
    std::vector<double>::iterator begin=vectorOfProbs.begin();
    std::vector<double>::iterator end=vectorOfProbs.end();
 
    for (std::vector<double>::iterator iter=begin;iter!=end;++iter)
    {
      sum+=*iter;
    }
 
    
    ATH_MSG_VERBOSE(" Sum of cluster probabilities is: "<<sum);
 
    std::vector<double> vectorOfSplitProbs;
 
    for (std::vector<double>::iterator iter=begin;iter!=end;++iter)
    {
      (*iter)/=sum;
    }
 
    if (m_priorMultiplicityContent.size()<vectorOfProbs.size())
    {
      ATH_MSG_ERROR("Prior compatibilities count " <<  m_priorMultiplicityContent.size() << " is too small: please correct through job properties.");
      return InDet::PixelClusterSplitProb(std::vector<double>());
    }
 
    double psum=0;
    int count=0;
    for (std::vector<double>::iterator iter=begin;iter!=end;++iter,++count)
    {
      psum+=(*iter)/m_priorMultiplicityContent[count];
    }
 
    count=0;
    for (std::vector<double>::iterator iter=begin;iter!=end;++iter,++count)
    {
      (*iter)/=m_priorMultiplicityContent[count];
      (*iter)/=psum;
    }
 
    double sumTest=0;
 
    for (std::vector<double>::iterator iter=begin;iter!=end;++iter)
    {
      ATH_MSG_VERBOSE("After update prob is: " << *iter);
      sumTest+=*iter;
    }
 
    ATH_MSG_VERBOSE(" Sum of cluster probabilities is: "<<sumTest);
    for (std::vector<double>::iterator iter=begin;iter!=end;++iter)
    {
      if (iter!=begin)
      {
        vectorOfSplitProbs.push_back(*iter);
      }
    }
 
    ATH_MSG_VERBOSE(" normalized P(1->2): " << vectorOfSplitProbs[0] << " P(2->3): " << vectorOfSplitProbs[1] );
 
    PixelClusterSplitProb  clusterSplitProb(vectorOfSplitProbs);
 
    ATH_MSG_VERBOSE("SplitProb: " << clusterSplitProb.splitProbability(2) << " -->3 " << clusterSplitProb.splitProbability(3) );
    
    return clusterSplitProb;
  }

initialize()

StatusCode InDet::TruthPixelClusterSplitProbTool::initialize ( )

overridevirtual

Definition at line 34 of file TruthPixelClusterSplitProbTool.cxx.

  {
 
    if (m_truthClusterizationFactory.retrieve().isFailure())
    {
      ATH_MSG_ERROR(" Unable to retrieve "<< m_truthClusterizationFactory );
      return StatusCode::FAILURE;
    }
 
    ATH_MSG_DEBUG(" Cluster split prob tool initialized successfully "<< m_truthClusterizationFactory );
    return StatusCode::SUCCESS;
  }

splitProbability() [1/2]

InDet::PixelClusterSplitProb InDet::TruthPixelClusterSplitProbTool::splitProbability ( const InDet::PixelCluster &  origCluster ) const

overridevirtual

Definition at line 48 of file TruthPixelClusterSplitProbTool.cxx.

  {
    std::vector<double> vectorOfProbs=m_truthClusterizationFactory->estimateNumberOfParticles(origCluster);
 
    ATH_MSG_VERBOSE(" Got splitProbability, size of vector: " << vectorOfProbs.size() );
 
    if (vectorOfProbs.empty())
    {
      std::vector<double> vectorOfSplitProbs;
      vectorOfSplitProbs.push_back(-100);
      PixelClusterSplitProb  clusterSplitProb(vectorOfSplitProbs);
      ATH_MSG_VERBOSE(" Returning single split prob equal to -100 " );
      return clusterSplitProb;
    }
     
 
    ATH_MSG_VERBOSE( 
        " P(1): " << vectorOfProbs[0] << 
        " P(2): " << vectorOfProbs[1] << 
        " P(>=3): " << vectorOfProbs[2] );
 
 
    return compileSplitProbability(vectorOfProbs); 
  }

splitProbability() [2/2]

InDet::PixelClusterSplitProb InDet::TruthPixelClusterSplitProbTool::splitProbability ( const InDet::PixelCluster &  origCluster, const Trk::TrackParameters &  trackParameters  ) const

overridevirtual

Definition at line 73 of file TruthPixelClusterSplitProbTool.cxx.

  {
    std::vector<double> vectorOfProbs=m_truthClusterizationFactory->estimateNumberOfParticles(origCluster);
 
    ATH_MSG_VERBOSE(" Got splitProbability, size of vector: " << vectorOfProbs.size() );
 
    if (vectorOfProbs.empty())
    {
      std::vector<double> vectorOfSplitProbs;
      vectorOfSplitProbs.push_back(-100);
      PixelClusterSplitProb  clusterSplitProb(vectorOfSplitProbs);
      ATH_MSG_VERBOSE(" Returning single split prob equal to -100 " );
      return clusterSplitProb;
    }
     
 
    ATH_MSG_VERBOSE( 
        " P(1): " << vectorOfProbs[0] << 
        " P(2): " << vectorOfProbs[1] << 
        " P(>=3): " << vectorOfProbs[2] );
 
 
    return compileSplitProbability(vectorOfProbs); 
  }

Member Data Documentation

m_priorMultiplicityContent

DoubleArrayProperty InDet::TruthPixelClusterSplitProbTool::m_priorMultiplicityContent {this, "PriorMultiplicityContent", {2793337, 82056, 19944}}

private

Definition at line 55 of file TruthPixelClusterSplitProbTool.h.

m_truthClusterizationFactory

ToolHandle<TruthClusterizationFactory> InDet::TruthPixelClusterSplitProbTool::m_truthClusterizationFactory {this, "NnClusterizationFactory", "InDet::NnClusterizationFactory/TruthClusterizationFactory"}

private

Definition at line 53 of file TruthPixelClusterSplitProbTool.h.

---

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

• TruthPixelClusterSplitProbTool.h
• TruthPixelClusterSplitProbTool.cxx