MultiValuedTransverseTrackHitsTool.cc

Go to the documentation of this file.

 
#include "Pandora/AlgorithmHeaders.h"
 
#include "larpandoracontent/LArThreeDReco/LArHitCreation/MultiValuedTransverseTrackHitsTool.h"
 
using namespace pandora;
 
namespace lar_content
{
 
void MultiValuedTransverseTrackHitsTool::GetTransverseTrackHit3D(const MatchedSlidingFitMap &matchedSlidingFitMap, ProtoHit &protoHit) const
{
    const CaloHit *const pCaloHit2D(protoHit.GetParentCaloHit2D());
    const HitType hitType(pCaloHit2D->GetHitType());
    const HitType hitType1((TPC_VIEW_U == hitType) ? TPC_VIEW_V : (TPC_VIEW_V == hitType) ? TPC_VIEW_W : TPC_VIEW_U);
    const HitType hitType2((TPC_VIEW_U == hitType) ? TPC_VIEW_W : (TPC_VIEW_V == hitType) ? TPC_VIEW_U : TPC_VIEW_V);
 
    CartesianPointVector fitPositionList1, fitPositionList2;
 
    MatchedSlidingFitMap::const_iterator iter1 = matchedSlidingFitMap.find(hitType1);
 
    if (matchedSlidingFitMap.end() != iter1)
    {
        const TwoDSlidingFitResult &fitResult1 = iter1->second;
        PANDORA_THROW_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
            fitResult1.GetGlobalFitPositionListAtX(pCaloHit2D->GetPositionVector().GetX(), fitPositionList1));
    }
 
    MatchedSlidingFitMap::const_iterator iter2 = matchedSlidingFitMap.find(hitType2);
 
    if (matchedSlidingFitMap.end() != iter2)
    {
        const TwoDSlidingFitResult &fitResult2 = iter2->second;
        PANDORA_THROW_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
            fitResult2.GetGlobalFitPositionListAtX(pCaloHit2D->GetPositionVector().GetX(), fitPositionList2));
    }
 
    unsigned int nViews(1);
    if (fitPositionList1.size() > 0)
        ++nViews;
    if (fitPositionList2.size() > 0)
        ++nViews;
 
    if (nViews < m_minViews)
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
    if (nViews <= 2)
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
    this->GetBestPosition3D(hitType1, hitType2, fitPositionList1, fitPositionList2, protoHit);
}
 
} // namespace lar_content