CACell Class Reference

#include <CACell.h>

Public Types
using	CAColl = std::vector< CACell >
using	CAntuple = std::vector< unsigned int >
using	CAntuplet = std::vector< unsigned int >
using	CAStatusColl = std::vector< CACellStatus >
using	Hit = RecHitsSortedInPhi::Hit

Public Member Functions
int	areAlignedRZ (float r1, float z1, float ro, float zo, const float ptmin, const float thetaCut) const
	CACell (const HitDoublets *doublets, int doubletId, const int innerHitId, const int outerHitId)
void	checkAlignmentAndAct (CAColl &allCells, CAntuple &innerCells, const float ptmin, const float region_origin_x, const float region_origin_y, const float region_origin_radius, const CACut::CAValuesByInnerLayerIds &thetaCutByInnerLayer, const CACut::CAValuesByInnerLayerIds &phiCutByInnerLayer, const float hardPtCut, std::vector< CACell::CAntuplet > *foundTriplets)
void	checkAlignmentAndPushTriplet (CAColl &allCells, CAntuple &innerCells, std::vector< CACell::CAntuplet > &foundTriplets, const float ptmin, const float region_origin_x, const float region_origin_y, const float region_origin_radius, const CACut::CAValuesByInnerLayerIds &thetaCutByInnerLayer, const CACut::CAValuesByInnerLayerIds &phiCutByInnerLayer, const float hardPtCut)
void	checkAlignmentAndTag (CAColl &allCells, CAntuple &innerCells, const float ptmin, const float region_origin_x, const float region_origin_y, const float region_origin_radius, const CACut::CAValuesByInnerLayerIds &thetaCutByInnerLayer, const CACut::CAValuesByInnerLayerIds &phiCutByInnerLayer, const float hardPtCut)
void	evolve (unsigned int me, CAStatusColl &allStatus)
void	findNtuplets (CAColl &allCells, std::vector< CAntuplet > &foundNtuplets, CAntuplet &tmpNtuplet, const unsigned int minHitsPerNtuplet) const
Hit const &	getInnerHit () const
int	getInnerLayer () const
float	getInnerPhi () const
float	getInnerR () const
float	getInnerX () const
float	getInnerY () const
float	getInnerZ () const
Hit const &	getOuterHit () const
int	getOuterLayer () 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
void	tagAsOuterNeighbor (unsigned int otherCell)

Private Attributes
const int	theDoubletId
const HitDoublets *	theDoublets
const float	theInnerR
const float	theInnerZ
CAntuple	theOuterNeighbors

Detailed Description

Definition at line 27 of file CACell.h.

Member Typedef Documentation

CAColl

using CACell::CAColl = std::vector<CACell>

Definition at line 32 of file CACell.h.

CAntuple

using CACell::CAntuple = std::vector<unsigned int>

Definition at line 30 of file CACell.h.

CAntuplet

using CACell::CAntuplet = std::vector<unsigned int>

Definition at line 31 of file CACell.h.

CAStatusColl

using CACell::CAStatusColl = std::vector<CACellStatus>

Definition at line 33 of file CACell.h.

Hit

using CACell::Hit = RecHitsSortedInPhi::Hit

Definition at line 29 of file CACell.h.

Constructor & Destructor Documentation

CACell()

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

inline

Definition at line 35 of file CACell.h.

      : theDoublets(doublets),
        theDoubletId(doubletId),
        theInnerR(doublets->rv(doubletId, HitDoublets::inner)),
        theInnerZ(doublets->z(doubletId, HitDoublets::inner)) {}

Member Function Documentation

areAlignedRZ()

int CACell::areAlignedRZ ( float  r1, float  z1, float  ro, float  zo, const float  ptmin, const float  thetaCut  ) const

inline

Definition at line 175 of file CACell.h.

References funct::abs(), getInnerR(), getInnerZ(), ALCARECOTkAlMinBias_cff::pMin, ptmin, diffTwoXMLs::r1, and mathSSE::sqrt().

Referenced by checkAlignmentAndAct().

                                                                                                          {
    float radius_diff = std::abs(r1 - ro);
    float distance_13_squared = radius_diff * radius_diff + (z1 - zo) * (z1 - zo);

    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 * (getInnerR() - ro) + getInnerZ() * (ro - r1) + zo * (r1 - getInnerR()));
    return tan_12_13_half_mul_distance_13_squared * pMin <= thetaCut * distance_13_squared * radius_diff;
  }

checkAlignmentAndAct()

void CACell::checkAlignmentAndAct ( CAColl &  allCells, CAntuple &  innerCells, const float  ptmin, const float  region_origin_x, const float  region_origin_y, const float  region_origin_radius, const CACut::CAValuesByInnerLayerIds &  thetaCutByInnerLayer, const CACut::CAValuesByInnerLayerIds &  phiCutByInnerLayer, const float  hardPtCut, std::vector< CACell::CAntuplet > *  foundTriplets  )

inline

Definition at line 82 of file CACell.h.

References areAlignedRZ(), CACut::CAValuesByInnerLayerIds::at(), getOuterR(), getOuterZ(), haveSimilarCurvature(), mps_fire::i, dqmiolumiharvest::j, heppy_loop::loop, convertSQLiteXML::ok, ptmin, and diffTwoXMLs::r1.

Referenced by checkAlignmentAndPushTriplet(), and checkAlignmentAndTag().

                                                                       {
    int ncells = innerCells.size();
    int constexpr VSIZE = 16;
    int ok[VSIZE];
    float r1[VSIZE];
    float z1[VSIZE];
    float thetaCut[VSIZE];
    float phiCut[VSIZE];
    auto ro = getOuterR();
    auto zo = getOuterZ();
    unsigned int cellId = this - &allCells.front();
    auto loop = [&](int i, int vs) {
      for (int j = 0; j < vs; ++j) {
        auto koc = innerCells[i + j];
        auto& oc = allCells[koc];
        r1[j] = oc.getInnerR();
        z1[j] = oc.getInnerZ();
        thetaCut[j] = thetaCutByInnerLayer.at(oc.getInnerLayer());
        phiCut[j] = phiCutByInnerLayer.at(oc.getInnerLayer());
      }
      // this vectorize!
      for (int j = 0; j < vs; ++j)
        ok[j] = areAlignedRZ(r1[j], z1[j], ro, zo, ptmin, thetaCut[j]);
      for (int j = 0; j < vs; ++j) {
        auto koc = innerCells[i + j];
        auto& oc = allCells[koc];
        if (ok[j] && haveSimilarCurvature(
                         oc, ptmin, region_origin_x, region_origin_y, region_origin_radius, phiCut[j], hardPtCut)) {
          if (foundTriplets)
            foundTriplets->emplace_back(CACell::CAntuplet{koc, cellId});
          else {
            oc.tagAsOuterNeighbor(cellId);
          }
        }
      }
    };
    auto lim = VSIZE * (ncells / VSIZE);
    for (int i = 0; i < lim; i += VSIZE)
      loop(i, VSIZE);
    loop(lim, ncells - lim);
  }

checkAlignmentAndPushTriplet()

void CACell::checkAlignmentAndPushTriplet ( CAColl &  allCells, CAntuple &  innerCells, std::vector< CACell::CAntuplet > &  foundTriplets, const float  ptmin, const float  region_origin_x, const float  region_origin_y, const float  region_origin_radius, const CACut::CAValuesByInnerLayerIds &  thetaCutByInnerLayer, const CACut::CAValuesByInnerLayerIds &  phiCutByInnerLayer, const float  hardPtCut  )

inline

Definition at line 153 of file CACell.h.

References checkAlignmentAndAct(), and ptmin.

                                                           {
    checkAlignmentAndAct(allCells,
                         innerCells,
                         ptmin,
                         region_origin_x,
                         region_origin_y,
                         region_origin_radius,
                         thetaCutByInnerLayer,
                         phiCutByInnerLayer,
                         hardPtCut,
                         &foundTriplets);
  }

checkAlignmentAndTag()

void CACell::checkAlignmentAndTag ( CAColl &  allCells, CAntuple &  innerCells, const float  ptmin, const float  region_origin_x, const float  region_origin_y, const float  region_origin_radius, const CACut::CAValuesByInnerLayerIds &  thetaCutByInnerLayer, const CACut::CAValuesByInnerLayerIds &  phiCutByInnerLayer, const float  hardPtCut  )

inline

Definition at line 133 of file CACell.h.

References checkAlignmentAndAct(), and ptmin.

                                                   {
    checkAlignmentAndAct(allCells,
                         innerCells,
                         ptmin,
                         region_origin_x,
                         region_origin_y,
                         region_origin_radius,
                         thetaCutByInnerLayer,
                         phiCutByInnerLayer,
                         hardPtCut,
                         nullptr);
  }

evolve()

void CACell::evolve ( unsigned int  me, CAStatusColl &  allStatus  )

inline

Definition at line 69 of file CACell.h.

References hlt_dqm_clientPB-live_cfg::me, and theOuterNeighbors.

                                                        {
    allStatus[me].hasSameStateNeighbors = 0;
    auto mystate = allStatus[me].theCAState;

    for (auto oc : theOuterNeighbors) {
      if (allStatus[oc].getCAState() == mystate) {
        allStatus[me].hasSameStateNeighbors = 1;

        break;
      }
    }
  }

findNtuplets()

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

inline

Definition at line 260 of file CACell.h.

References caHitNtupletGeneratorKernels::foundNtuplets, mps_fire::i, HLT_2023v12_cff::minHitsPerNtuplet, 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 {
      unsigned int numberOfOuterNeighbors = theOuterNeighbors.size();
      for (unsigned int i = 0; i < numberOfOuterNeighbors; ++i) {
        tmpNtuplet.push_back((theOuterNeighbors[i]));
        allCells[theOuterNeighbors[i]].findNtuplets(allCells, foundNtuplets, tmpNtuplet, minHitsPerNtuplet);
        tmpNtuplet.pop_back();
      }
    }
  }

getInnerHit()

Hit const& CACell::getInnerHit ( ) const

inline

Definition at line 41 of file CACell.h.

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

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

getInnerLayer()

int CACell::getInnerLayer ( ) const

inline

Definition at line 45 of file CACell.h.

References HitDoublets::detLayer(), HitDoublets::inner, DetLayer::seqNum(), and theDoublets.

{ return theDoublets->detLayer(HitDoublets::inner)->seqNum(); }

getInnerPhi()

float CACell::getInnerPhi ( ) const

inline

Definition at line 65 of file CACell.h.

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

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

getInnerR()

float CACell::getInnerR ( ) const

inline

Definition at line 61 of file CACell.h.

References theInnerR.

Referenced by areAlignedRZ().

{ return theInnerR; }

getInnerX()

float CACell::getInnerX ( ) const

inline

Definition at line 49 of file CACell.h.

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

Referenced by haveSimilarCurvature().

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

getInnerY()

float CACell::getInnerY ( ) const

inline

Definition at line 53 of file CACell.h.

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

Referenced by haveSimilarCurvature().

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

getInnerZ()

float CACell::getInnerZ ( ) const

inline

Definition at line 57 of file CACell.h.

References theInnerZ.

Referenced by areAlignedRZ().

{ return theInnerZ; }

getOuterHit()

Hit const& CACell::getOuterHit ( ) const

inline

Definition at line 43 of file CACell.h.

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

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

getOuterLayer()

int CACell::getOuterLayer ( ) const

inline

Definition at line 47 of file CACell.h.

References HitDoublets::detLayer(), HitDoublets::outer, DetLayer::seqNum(), and theDoublets.

{ return theDoublets->detLayer(HitDoublets::outer)->seqNum(); }

getOuterPhi()

float CACell::getOuterPhi ( ) const

inline

Definition at line 67 of file CACell.h.

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

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

getOuterR()

float CACell::getOuterR ( ) const

inline

Definition at line 63 of file CACell.h.

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

Referenced by checkAlignmentAndAct().

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

getOuterX()

float CACell::getOuterX ( ) const

inline

Definition at line 51 of file CACell.h.

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

Referenced by haveSimilarCurvature().

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

getOuterY()

float CACell::getOuterY ( ) const

inline

Definition at line 55 of file CACell.h.

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

Referenced by haveSimilarCurvature().

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

getOuterZ()

float CACell::getOuterZ ( ) const

inline

Definition at line 59 of file CACell.h.

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

Referenced by checkAlignmentAndAct().

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

haveSimilarCurvature()

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 188 of file CACell.h.

References funct::abs(), hippyaddtobaddatafiles::cd(), MillePedeFileConverter_cfg::e, f, getInnerX(), getInnerY(), getOuterX(), getOuterY(), hltrates_dqm_sourceclient-live_cfg::offset, ptmin, CosmicsPD_Skims::radius, mathSSE::sqrt(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by checkAlignmentAndAct().

                                                         {
    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));
    // high pt : just straight
    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;

      return distance_13_beamspot_squared < (region_origin_radius + phiCut) * (region_origin_radius + phiCut);
    }

    //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;  // FIXME move out and use real MagField

    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;  // hard cut on pt

    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;
  }

tagAsOuterNeighbor()

void CACell::tagAsOuterNeighbor ( unsigned int  otherCell )

inline

Definition at line 186 of file CACell.h.

References theOuterNeighbors.

{ theOuterNeighbors.push_back(otherCell); }

Member Data Documentation

theDoubletId

const int CACell::theDoubletId

private

Definition at line 285 of file CACell.h.

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

theDoublets

const HitDoublets* CACell::theDoublets

private

Definition at line 284 of file CACell.h.

Referenced by getInnerHit(), getInnerLayer(), getInnerPhi(), getInnerX(), getInnerY(), getOuterHit(), getOuterLayer(), getOuterPhi(), getOuterR(), getOuterX(), getOuterY(), and getOuterZ().

theInnerR

const float CACell::theInnerR

private

Definition at line 287 of file CACell.h.

Referenced by getInnerR().

theInnerZ

const float CACell::theInnerZ

private

Definition at line 288 of file CACell.h.

Referenced by getInnerZ().

theOuterNeighbors

CAntuple CACell::theOuterNeighbors

private

Definition at line 282 of file CACell.h.

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