CACell Class Reference

#include <CACell.h>

Public Types
using	CAntuplet = std::vector< CACell * >
using	Hit = RecHitsSortedInPhi::Hit

Public Member Functions
bool	areAlignedRZ (const CACell *otherCell, const float ptmin, const float thetaCut) const
	CACell (const HitDoublets *doublets, int doubletId, const unsigned int cellId, const int innerHitId, const int outerHitId)
void	checkAlignmentAndPushTriplet (CACell *innerCell, std::vector< CACell::CAntuplet > &foundTriplets, const float ptmin, const float region_origin_x, const float region_origin_y, const float region_origin_radius, const float thetaCut, const float phiCut, const float hardPtCut)
void	checkAlignmentAndTag (CACell *innerCell, const float ptmin, const float region_origin_x, const float region_origin_y, const float region_origin_radius, const float thetaCut, const float phiCut, const float hardPtCut)
void	evolve ()
void	findNtuplets (std::vector< CAntuplet > &foundNtuplets, CAntuplet &tmpNtuplet, const unsigned int minHitsPerNtuplet) const
unsigned int	getCAState () const
unsigned int	getCellId () const
Hit const &	getInnerHit () const
unsigned int	getInnerHitId () const
float	getInnerPhi () const
float	getInnerR () const
float	getInnerX () const
float	getInnerY () const
float	getInnerZ () const
Hit const &	getOuterHit () const
unsigned int	getOuterHitId () const
float	getOuterPhi () const
float	getOuterR () const
float	getOuterX () const
float	getOuterY () const
float	getOuterZ () const
bool	haveSimilarCurvature (const CACell *otherCell, const float ptmin, const float region_origin_x, const float region_origin_y, const float region_origin_radius, const float phiCut, const float hardPtCut) const
bool	isRootCell (const unsigned int minimumCAState) const
void	tagAsInnerNeighbor (CACell *otherCell)
void	tagAsOuterNeighbor (CACell *otherCell)
void	updateState ()

Private Attributes
unsigned int	hasSameStateNeighbors
unsigned int	theCAState
const unsigned int	theCellId
const int	theDoubletId
const HitDoublets *	theDoublets
const unsigned int	theInnerHitId
std::vector< CACell * >	theInnerNeighbors
const float	theInnerR
const float	theInnerZ
const unsigned int	theOuterHitId
std::vector< CACell * >	theOuterNeighbors
const float	theOuterR
const float	theOuterZ

Detailed Description

Definition at line 13 of file CACell.h.

Member Typedef Documentation

using CACell::CAntuplet = std::vector<CACell*>

Definition at line 16 of file CACell.h.

using CACell::Hit = RecHitsSortedInPhi::Hit

Definition at line 15 of file CACell.h.

Constructor & Destructor Documentation

CACell::CACell ( const HitDoublets *  doublets, int  doubletId, const unsigned int  cellId, const int  innerHitId, const int  outerHitId  )

inline

Definition at line 18 of file CACell.h.

                                                                                                                              :
    theCAState(0), theInnerHitId(innerHitId), theOuterHitId(outerHitId), theCellId(cellId), hasSameStateNeighbors(0), theDoublets(doublets), theDoubletId(doubletId),
    theInnerR(doublets->r(doubletId, HitDoublets::inner)), theOuterR(doublets->r(doubletId, HitDoublets::outer)),
    theInnerZ(doublets->z(doubletId, HitDoublets::inner)), theOuterZ(doublets->z(doubletId, HitDoublets::outer)) {
    }

Member Function Documentation

bool CACell::areAlignedRZ ( const CACell *  otherCell, const float  ptmin, const float  thetaCut  ) const

inline

Definition at line 129 of file CACell.h.

References funct::abs(), getInnerR(), getInnerZ(), HLT_FULL_cff::pMin, diffTwoXMLs::r1, mathSSE::sqrt(), theInnerR, theInnerZ, theOuterR, and theOuterZ.

Referenced by checkAlignmentAndPushTriplet(), and checkAlignmentAndTag().

    {

        float r1 = otherCell->getInnerR();
        float z1 = otherCell->getInnerZ();
        float radius_diff = std::abs(r1 - theOuterR);
        float distance_13_squared = radius_diff*radius_diff + (z1 - theOuterZ)*(z1 - theOuterZ);

        float pMin = ptmin*std::sqrt(distance_13_squared); //this needs to be divided by radius_diff later

        float tan_12_13_half_mul_distance_13_squared = fabs(z1 * (theInnerR - theOuterR) + theInnerZ * (theOuterR - r1) + theOuterZ * (r1 - theInnerR)) ;
        return tan_12_13_half_mul_distance_13_squared * pMin <= thetaCut * distance_13_squared * radius_diff;
    }

void CACell::checkAlignmentAndPushTriplet ( CACell *  innerCell, std::vector< CACell::CAntuplet > &  foundTriplets, const float  ptmin, const float  region_origin_x, const float  region_origin_y, const float  region_origin_radius, const float  thetaCut, const float  phiCut, const float  hardPtCut  )

inline

Definition at line 114 of file CACell.h.

References areAlignedRZ(), and haveSimilarCurvature().

                                   {

        if (areAlignedRZ(innerCell, ptmin, thetaCut) &&
                haveSimilarCurvature(innerCell,ptmin, region_origin_x, region_origin_y,
                        region_origin_radius, phiCut, hardPtCut)) {
            foundTriplets.emplace_back(CACell::CAntuplet{innerCell,this});

        }
    }

void CACell::checkAlignmentAndTag ( CACell *  innerCell, const float  ptmin, const float  region_origin_x, const float  region_origin_y, const float  region_origin_radius, const float  thetaCut, const float  phiCut, const float  hardPtCut  )

inline

Definition at line 102 of file CACell.h.

References areAlignedRZ(), haveSimilarCurvature(), tagAsInnerNeighbor(), and tagAsOuterNeighbor().

                                                       {

        if (areAlignedRZ(innerCell, ptmin, thetaCut) &&
                haveSimilarCurvature(innerCell,ptmin, region_origin_x, region_origin_y,
                        region_origin_radius, phiCut, hardPtCut)) {
            tagAsInnerNeighbor(innerCell);
            innerCell->tagAsOuterNeighbor(this);
        }
    }

void CACell::evolve ( )

inline

Definition at line 84 of file CACell.h.

References getCAState(), hasSameStateNeighbors, i, theCAState, and theOuterNeighbors.

                  {

        hasSameStateNeighbors = 0;
        unsigned int numberOfNeighbors = theOuterNeighbors.size();

        for (unsigned int i = 0; i < numberOfNeighbors; ++i) {

            if (theOuterNeighbors[i]->getCAState() == theCAState) {

                hasSameStateNeighbors = 1;

                break;
            }
        }

    }

void CACell::findNtuplets ( std::vector< CAntuplet > &  foundNtuplets, CAntuplet &  tmpNtuplet, const unsigned int  minHitsPerNtuplet  ) const

inline

Definition at line 242 of file CACell.h.

References i, and theOuterNeighbors.

                                                                                                                              {

        // the building process for a track ends if:
        // it has no outer neighbor
        // it has no compatible neighbor
        // the ntuplets is then saved if the number of hits it contains is greater than a threshold

        if (tmpNtuplet.size() == minHitsPerNtuplet - 1)
        {
            foundNtuplets.push_back(tmpNtuplet);
        }
        else
        {
            if(theOuterNeighbors.size() == 0)
            {
                return;
            }
            else
            {
                unsigned int numberOfOuterNeighbors = theOuterNeighbors.size();
                for (unsigned int i = 0; i < numberOfOuterNeighbors; ++i) {
                    tmpNtuplet.push_back((theOuterNeighbors[i]));
                    theOuterNeighbors[i]->findNtuplets(foundNtuplets, tmpNtuplet, minHitsPerNtuplet);
                    tmpNtuplet.pop_back();
                }
            }

        }

    }

unsigned int CACell::getCAState ( ) const

inline

Definition at line 226 of file CACell.h.

References theCAState.

Referenced by evolve().

                                    {
        return theCAState;
    }

unsigned int CACell::getCellId ( ) const

inline

Definition at line 24 of file CACell.h.

References theCellId.

                                   {
        return theCellId;
    }

Hit const& CACell::getInnerHit ( ) const

inline

Definition at line 28 of file CACell.h.

References HitDoublets::hit(), HitDoublets::inner, theDoubletId, and theDoublets.

                                    {
        return theDoublets->hit(theDoubletId, HitDoublets::inner);
    }

unsigned int CACell::getInnerHitId ( ) const

inline

Definition at line 76 of file CACell.h.

References theInnerHitId.

                                       {
        return theInnerHitId;
    }

float CACell::getInnerPhi ( ) const

inline

Definition at line 68 of file CACell.h.

References HitDoublets::inner, HitDoublets::phi(), theDoubletId, and theDoublets.

                              {
        return theDoublets->phi(theDoubletId, HitDoublets::inner);
    }

float CACell::getInnerR ( ) const

inline

Definition at line 60 of file CACell.h.

References theInnerR.

Referenced by areAlignedRZ().

                            {
        return theInnerR;
    }

float CACell::getInnerX ( ) const

inline

Definition at line 36 of file CACell.h.

References HitDoublets::inner, theDoubletId, theDoublets, and HitDoublets::x().

Referenced by haveSimilarCurvature().

                            {
        return theDoublets->x(theDoubletId, HitDoublets::inner);
    }

float CACell::getInnerY ( ) const

inline

Definition at line 44 of file CACell.h.

References HitDoublets::inner, theDoubletId, theDoublets, and HitDoublets::y().

Referenced by haveSimilarCurvature().

                            {
        return theDoublets->y(theDoubletId, HitDoublets::inner);
    }

float CACell::getInnerZ ( ) const

inline

Definition at line 52 of file CACell.h.

References theInnerZ.

Referenced by areAlignedRZ().

                            {
        return theInnerZ;
    }

Hit const& CACell::getOuterHit ( ) const

inline

Definition at line 32 of file CACell.h.

References HitDoublets::hit(), HitDoublets::outer, theDoubletId, and theDoublets.

                                    {
        return theDoublets->hit(theDoubletId, HitDoublets::outer);
    }

unsigned int CACell::getOuterHitId ( ) const

inline

Definition at line 80 of file CACell.h.

References theOuterHitId.

                                       {
        return theOuterHitId;
    }

float CACell::getOuterPhi ( ) const

inline

Definition at line 72 of file CACell.h.

References HitDoublets::outer, HitDoublets::phi(), theDoubletId, and theDoublets.

                              {
        return theDoublets->phi(theDoubletId, HitDoublets::outer);
    }

float CACell::getOuterR ( ) const

inline

Definition at line 64 of file CACell.h.

References theOuterR.

                            {
        return theOuterR;
    }

float CACell::getOuterX ( ) const

inline

Definition at line 40 of file CACell.h.

References HitDoublets::outer, theDoubletId, theDoublets, and HitDoublets::x().

Referenced by haveSimilarCurvature().

                            {
        return theDoublets->x(theDoubletId, HitDoublets::outer);
    }

float CACell::getOuterY ( ) const

inline

Definition at line 48 of file CACell.h.

References HitDoublets::outer, theDoubletId, theDoublets, and HitDoublets::y().

Referenced by haveSimilarCurvature().

                            {
        return theDoublets->y(theDoubletId, HitDoublets::outer);
    }

float CACell::getOuterZ ( ) const

inline

Definition at line 56 of file CACell.h.

References theOuterZ.

                            {
        return theOuterZ;
    }

bool CACell::haveSimilarCurvature ( const CACell *  otherCell, const float  ptmin, const float  region_origin_x, const float  region_origin_y, const float  region_origin_radius, const float  phiCut, const float  hardPtCut  ) const

inline

Definition at line 153 of file CACell.h.

References funct::abs(), alignCSCRings::e, f, getInnerX(), getInnerY(), getOuterX(), getOuterY(), hltrates_dqm_sourceclient-live_cfg::offset, CosmicsPD_Skims::radius, and mathSSE::sqrt().

Referenced by checkAlignmentAndPushTriplet(), and checkAlignmentAndTag().

    {


        auto x1 = otherCell->getInnerX();
        auto y1 = otherCell->getInnerY();

        auto x2 = getInnerX();
        auto y2 = getInnerY();

        auto x3 = getOuterX();
        auto y3 = getOuterY();

        float distance_13_squared = (x1 - x3)*(x1 - x3) + (y1 - y3)*(y1 - y3);
        float tan_12_13_half_mul_distance_13_squared = std::abs(y1 * (x2 - x3) + y2 * (x3 - x1) + y3 * (x1 - x2)) ;
        if(tan_12_13_half_mul_distance_13_squared * ptmin <= 1.0e-4f*distance_13_squared)
        {

            float distance_3_beamspot_squared = (x3-region_origin_x) * (x3-region_origin_x) + (y3-region_origin_y) * (y3-region_origin_y);

            float dot_bs3_13 = ((x1 - x3)*( region_origin_x - x3) + (y1 - y3) * (region_origin_y-y3));
            float proj_bs3_on_13_squared = dot_bs3_13*dot_bs3_13/distance_13_squared;

            float distance_13_beamspot_squared  = distance_3_beamspot_squared -  proj_bs3_on_13_squared;

            if(distance_13_beamspot_squared > (region_origin_radius+phiCut)*(region_origin_radius+phiCut) )
            { 
                return false;
            }
            return true;
        }

        //87 cm/GeV = 1/(3.8T * 0.3)

        //take less than radius given by the hardPtCut and reject everything below
        float minRadius = hardPtCut*87.f;

        auto det = (x1 - x2) * (y2 - y3) - (x2 - x3) * (y1 - y2);


        auto offset = x2 * x2 + y2*y2;

        auto bc = (x1 * x1 + y1 * y1 - offset)*0.5f;

        auto cd = (offset - x3 * x3 - y3 * y3)*0.5f;



        auto idet = 1.f / det;

        auto x_center = (bc * (y2 - y3) - cd * (y1 - y2)) * idet;
        auto y_center = (cd * (x1 - x2) - bc * (x2 - x3)) * idet;

        auto radius = std::sqrt((x2 - x_center)*(x2 - x_center) + (y2 - y_center)*(y2 - y_center));

        if(radius < minRadius)
            return false;
        auto centers_distance_squared = (x_center - region_origin_x)*(x_center - region_origin_x) + (y_center - region_origin_y)*(y_center - region_origin_y);
        auto region_origin_radius_plus_tolerance = region_origin_radius + phiCut;
        auto minimumOfIntersectionRange = (radius - region_origin_radius_plus_tolerance)*(radius - region_origin_radius_plus_tolerance);

        if (centers_distance_squared >= minimumOfIntersectionRange) {
            auto maximumOfIntersectionRange = (radius + region_origin_radius_plus_tolerance)*(radius + region_origin_radius_plus_tolerance);
            return centers_distance_squared <= maximumOfIntersectionRange;
        } 

        return false;
        


    }

bool CACell::isRootCell ( const unsigned int  minimumCAState ) const

inline

Definition at line 235 of file CACell.h.

References theCAState.

                                                             {
        return (theCAState >= minimumCAState);
    }

void CACell::tagAsInnerNeighbor ( CACell *  otherCell )

inline

Definition at line 148 of file CACell.h.

References theInnerNeighbors.

Referenced by checkAlignmentAndTag().

    {
        theInnerNeighbors.push_back(otherCell);
    }

void CACell::tagAsOuterNeighbor ( CACell *  otherCell )

inline

Definition at line 143 of file CACell.h.

References theOuterNeighbors.

Referenced by checkAlignmentAndTag().

    {
        theOuterNeighbors.push_back(otherCell);
    }

void CACell::updateState ( )

inline

Definition at line 231 of file CACell.h.

References hasSameStateNeighbors, and theCAState.

                       {
        theCAState += hasSameStateNeighbors;
    }

Member Data Documentation

unsigned int CACell::hasSameStateNeighbors

private

Definition at line 283 of file CACell.h.

Referenced by evolve(), and updateState().

unsigned int CACell::theCAState

private

Definition at line 278 of file CACell.h.

Referenced by evolve(), getCAState(), isRootCell(), and updateState().

const unsigned int CACell::theCellId

private

Definition at line 282 of file CACell.h.

Referenced by getCellId().

const int CACell::theDoubletId

private

Definition at line 286 of file CACell.h.

Referenced by getInnerHit(), getInnerPhi(), getInnerX(), getInnerY(), getOuterHit(), getOuterPhi(), getOuterX(), and getOuterY().

const HitDoublets* CACell::theDoublets

private

Definition at line 285 of file CACell.h.

Referenced by getInnerHit(), getInnerPhi(), getInnerX(), getInnerY(), getOuterHit(), getOuterPhi(), getOuterX(), and getOuterY().

const unsigned int CACell::theInnerHitId

private

Definition at line 280 of file CACell.h.

Referenced by getInnerHitId().

std::vector<CACell*> CACell::theInnerNeighbors

private

Definition at line 275 of file CACell.h.

Referenced by tagAsInnerNeighbor().

const float CACell::theInnerR

private

Definition at line 288 of file CACell.h.

Referenced by areAlignedRZ(), and getInnerR().

const float CACell::theInnerZ

private

Definition at line 290 of file CACell.h.

Referenced by areAlignedRZ(), and getInnerZ().

const unsigned int CACell::theOuterHitId

private

Definition at line 281 of file CACell.h.

Referenced by getOuterHitId().

std::vector<CACell*> CACell::theOuterNeighbors

private

Definition at line 276 of file CACell.h.

Referenced by evolve(), findNtuplets(), and tagAsOuterNeighbor().

const float CACell::theOuterR

private

Definition at line 289 of file CACell.h.

Referenced by areAlignedRZ(), and getOuterR().

const float CACell::theOuterZ

private

Definition at line 291 of file CACell.h.

Referenced by areAlignedRZ(), and getOuterZ().