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 float thetaCut, const float phiCut, 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 float thetaCut, const float phiCut, 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 float thetaCut, const float phiCut, 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
float	getInnerPhi () const
float	getInnerR () const
float	getInnerX () const
float	getInnerY () const
float	getInnerZ () const
Hit const &	getOuterHit () 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 36 of file CACell.h.

Member Typedef Documentation

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

Definition at line 41 of file CACell.h.

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

Definition at line 39 of file CACell.h.

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

Definition at line 40 of file CACell.h.

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

Definition at line 42 of file CACell.h.

using CACell::Hit = RecHitsSortedInPhi::Hit

Definition at line 38 of file CACell.h.

Constructor & Destructor Documentation

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

inline

Definition at line 45 of file CACell.h.

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

Member Function Documentation

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

inline

Definition at line 174 of file CACell.h.

References funct::abs(), ALCARECOTkAlMinBias_cff::pMin, and mathSSE::sqrt().

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

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 float  thetaCut, const float  phiCut, const float  hardPtCut, std::vector< CACell::CAntuplet > *  foundTriplets  )

inline

Definition at line 120 of file CACell.h.

References constexpr, mps_fire::i, heppy_loop::loop, convertSQLiteXML::ok, and diffTwoXMLs::r1.

                                                                                                       {
    int ncells = innerCells.size();
    int constexpr VSIZE = 16;
    int ok[VSIZE];
    float r1[VSIZE];
    float z1[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();
      }
      // this vectorize!
      for (int j=0;j<vs; ++j) ok[j] = areAlignedRZ(r1[j], z1[j], ro, zo, ptmin, thetaCut);
      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, 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);
    
  }

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 float  thetaCut, const float  phiCut, const float  hardPtCut  )

inline

Definition at line 165 of file CACell.h.

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

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 float  thetaCut, const float  phiCut, const float  hardPtCut  )

inline

Definition at line 158 of file CACell.h.

                                                           {
    checkAlignmentAndAct(allCells, innerCells, ptmin, region_origin_x, region_origin_y, region_origin_radius, thetaCut,
             phiCut, hardPtCut, nullptr);
    
  }

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

inline

Definition at line 102 of file CACell.h.

References CACellStatus::getCAState().

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

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

inline

Definition at line 264 of file CACell.h.

References mps_fire::i.

                                                                                                                                              {
    
    // 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();
    }
      }
    
  }

Hit const& CACell::getInnerHit ( ) const

inline

Definition at line 53 of file CACell.h.

References HitDoublets::inner.

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

float CACell::getInnerPhi ( ) const

inline

Definition at line 94 of file CACell.h.

References HitDoublets::inner.

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

float CACell::getInnerR ( ) const

inline

Definition at line 86 of file CACell.h.

                          {
    return theInnerR;
  }

float CACell::getInnerX ( ) const

inline

Definition at line 62 of file CACell.h.

References HitDoublets::inner.

Referenced by haveSimilarCurvature().

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

float CACell::getInnerY ( ) const

inline

Definition at line 70 of file CACell.h.

References HitDoublets::inner.

Referenced by haveSimilarCurvature().

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

float CACell::getInnerZ ( ) const

inline

Definition at line 78 of file CACell.h.

                          {
    return theInnerZ;
  }

Hit const& CACell::getOuterHit ( ) const

inline

Definition at line 57 of file CACell.h.

References HitDoublets::outer.

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

float CACell::getOuterPhi ( ) const

inline

Definition at line 98 of file CACell.h.

References HitDoublets::outer.

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

float CACell::getOuterR ( ) const

inline

Definition at line 90 of file CACell.h.

References HitDoublets::outer.

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

float CACell::getOuterX ( ) const

inline

Definition at line 66 of file CACell.h.

References HitDoublets::outer.

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

float CACell::getOuterY ( ) const

inline

Definition at line 74 of file CACell.h.

References HitDoublets::outer.

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

float CACell::getOuterZ ( ) const

inline

Definition at line 82 of file CACell.h.

References HitDoublets::outer.

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

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

References funct::abs(), MillePedeFileConverter_cfg::e, f, getInnerX(), getInnerY(), PFRecoTauDiscriminationByIsolation_cfi::offset, TCMET_cfi::radius, mathSSE::sqrt(), globals_cff::x1, and globals_cff::x2.

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

void CACell::tagAsOuterNeighbor ( unsigned int  otherCell )

inline

Definition at line 186 of file CACell.h.

  {
    theOuterNeighbors.push_back(otherCell);
  }

Member Data Documentation

const int CACell::theDoubletId

private

Definition at line 293 of file CACell.h.

const HitDoublets* CACell::theDoublets

private

Definition at line 292 of file CACell.h.

const float CACell::theInnerR

private

Definition at line 295 of file CACell.h.

const float CACell::theInnerZ

private

Definition at line 296 of file CACell.h.

CAntuple CACell::theOuterNeighbors

private

Definition at line 290 of file CACell.h.