Trk::CaloCellSelectorMinPerp Class Reference

#include <CaloCellSelectorMinPerp.h>

Inheritance diagram for Trk::CaloCellSelectorMinPerp:

Collaboration diagram for Trk::CaloCellSelectorMinPerp:

Public Member Functions
	CaloCellSelectorMinPerp (double coneSize)
	~CaloCellSelectorMinPerp ()=default
bool	preSelectAction (const Trk::CaloExtension &caloExtension)
	Initialize selector wtih CaloExtension, returns false if the selection cannot be performed for the given input. More...
bool	select (const CaloCell &cell) const
	Method to determine if a cell is to be selected or not. More...
virtual bool	preSelectActionCells (const Rec::ParticleCellAssociation &cellAssociation)
	Initialize selector wtih ParticleCellAssociation, returns false if the selection cannot be performed for the given input More...

Private Attributes
const Trk::CaloExtension *	m_caloExtension
double	m_coneSize
Amg::Vector3D	m_meanPos
double	m_perp2cut

Detailed Description

Definition at line 22 of file CaloCellSelectorMinPerp.h.

Constructor & Destructor Documentation

CaloCellSelectorMinPerp()

Trk::CaloCellSelectorMinPerp::CaloCellSelectorMinPerp

(

double

coneSize

)

Definition at line 17 of file CaloCellSelectorMinPerp.cxx.

: m_caloExtension(nullptr), m_coneSize(coneSize), m_perp2cut(0) {}

~CaloCellSelectorMinPerp()

Trk::CaloCellSelectorMinPerp::~CaloCellSelectorMinPerp ( )

default

Member Function Documentation

preSelectAction()

bool Trk::CaloCellSelectorMinPerp::preSelectAction ( const Trk::CaloExtension &  caloExtension )

virtual

Initialize selector wtih CaloExtension, returns false if the selection cannot be performed for the given input.

Reimplemented from Trk::ICaloCellSelector.

Definition at line 19 of file CaloCellSelectorMinPerp.cxx.

                                                                                       {
        m_caloExtension = &caloExtension;
 
        // find a circular cone (defined by cut in perpendicular distance form track i.e. perp2)
        // with approx area as the dR cone
        // the dR cone is roughly thought as an elipse
        if (m_caloExtension->caloLayerIntersections().empty()) return false;
 
        Amg::Vector3D pos = m_caloExtension->caloLayerIntersections()[0].position();
 
        double a = fabs(pos.eta());
        double b = exp(-a);                // = tan(theta/2)
        double c = 2. * b / (1. - b * b);  // = tan(theta)
        double r = c / sqrt(1. + c * c);   // distance from beam axis
        double t1 = 2. * atan(exp(-(a - m_coneSize)));
        double t2 = 2. * atan(exp(-(a + m_coneSize)));
        double dRad1 = fabs(t1 - t2) / 2.;
        // double dRad1 = 2.*b/(1.+b*b) * m_coneSize; //semi-major axis of the dR cone elipse
        double dRad2 = r * m_coneSize;  // semi minor axis of the dR cone elipse
        m_perp2cut = dRad1 * dRad2;     // force area of circle = area of elipse, i.e. pi*r*r = pi*dRad1*dRad2
        return true;
    }

preSelectActionCells()

bool ICaloCellSelector::preSelectActionCells ( const Rec::ParticleCellAssociation &  cellAssociation )

inlinevirtualinherited

Initialize selector wtih ParticleCellAssociation, returns false if the selection cannot be performed for the given input

Definition at line 44 of file Reconstruction/RecoTools/RecoToolInterfaces/RecoToolInterfaces/ICaloCellSelector.h.

                                                                                                         {
    return preSelectAction(cellAssociation.caloExtension());
  }

select()

bool Trk::CaloCellSelectorMinPerp::select ( const CaloCell &  cell ) const

virtual

Method to determine if a cell is to be selected or not.

Parameters

CaloCell

Returns

bool

Implements Trk::ICaloCellSelector.

Definition at line 42 of file CaloCellSelectorMinPerp.cxx.

                                                                   {
        if (!m_caloExtension || m_caloExtension->caloLayerIntersections().empty()) return false;
 
        const CaloDetDescrElement* dde = cell.caloDDE();
        if (!dde) return false;
 
        Amg::Vector3D cellPos(dde->x(), dde->y(), dde->z());
 
        int nearestIdx;
        Amg::Vector3D nearestPos, nearestMom;
        Utils::findNearestPoint(cellPos, m_caloExtension, nearestIdx, nearestPos, nearestMom);
 
        // get the perp2 from track
        Amg::Vector3D dPos = cellPos - nearestPos;
        float perp2 = dPos.perp2(nearestMom);
 
        // get the total track length from IP to point of cloest approach to the cell
        // scale the perp2 cut with this length
        float totTrkLen = sqrt(dPos.mag2() - perp2);
        if (dPos.dot(nearestMom) < 0) { totTrkLen = -totTrkLen; }
 
        Amg::Vector3D pos;
        Amg::Vector3D oldPos(0., 0., 0.);
        const std::vector<Trk::CurvilinearParameters>& intersections = m_caloExtension->caloLayerIntersections();
        for (int i = 0; i <= nearestIdx; ++i) {
            pos = intersections[i].position();
            totTrkLen += (pos - oldPos).mag();
            std::swap(oldPos, pos);
        }
 
        // prevent cell at exact opposite of the track from being selected..
        if (totTrkLen < 0) return false;
 
        if (perp2 < (m_perp2cut * totTrkLen * totTrkLen)) { return true; }  // IF
 
        return false;
    }

Member Data Documentation

m_caloExtension

const Trk::CaloExtension* Trk::CaloCellSelectorMinPerp::m_caloExtension

private

Definition at line 32 of file CaloCellSelectorMinPerp.h.

m_coneSize

double Trk::CaloCellSelectorMinPerp::m_coneSize

private

Definition at line 33 of file CaloCellSelectorMinPerp.h.

m_meanPos

Amg::Vector3D Trk::CaloCellSelectorMinPerp::m_meanPos

private

Definition at line 34 of file CaloCellSelectorMinPerp.h.

m_perp2cut

double Trk::CaloCellSelectorMinPerp::m_perp2cut

private

Definition at line 35 of file CaloCellSelectorMinPerp.h.

---

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

• CaloCellSelectorMinPerp.h
• CaloCellSelectorMinPerp.cxx