Phase2EndcapLayer Class Reference

#include <Phase2EndcapLayer.h>

Inheritance diagram for Phase2EndcapLayer:

Classes
struct	RingPar

Public Member Functions
const std::vector< const GeomDet * > &	basicComponents () const override
const std::vector< const GeometricSearchDet * > &	components () const override __attribute__((cold))
void	groupedCompatibleDetsV (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetGroup > &result) const override __attribute__((hot))
Phase2EndcapLayer &	operator= (const Phase2EndcapLayer &)=delete
	Phase2EndcapLayer (std::vector< const Phase2EndcapRing * > &rings, const bool isOT) __attribute__((cold))
	Phase2EndcapLayer (const Phase2EndcapLayer &)=delete
SubDetector	subDetector () const override
	~Phase2EndcapLayer () override __attribute__((cold))

Private Member Functions
BoundDisk *	computeDisk (const std::vector< const Phase2EndcapRing * > &rings) const __attribute__((cold))
float	computeWindowSize (const GeomDet *det, const TrajectoryStateOnSurface &tsos, const MeasurementEstimator &est) const __attribute__((hot))
void	fillRingPars (int i) __attribute__((cold))
std::array< int, 3 >	findThreeClosest (std::vector< GlobalPoint >) const __attribute__((hot))
bool	overlapInR (const TrajectoryStateOnSurface &tsos, int i, double ymax) const __attribute__((hot))
std::array< int, 3 >	ringIndicesByCrossingProximity (const TrajectoryStateOnSurface &startingState, const Propagator &prop) const

Private Attributes
const bool	isOuterTracker
std::vector< RingPar >	ringPars
std::vector< GeomDet const * >	theBasicComps
std::atomic< std::vector< const GeometricSearchDet * > * >	theComponents
std::vector< const Phase2EndcapRing * >	theComps
int	theRingSize

Detailed Description

A concrete implementation for Phase 2 Endcap/Forward layer built out of Phase2EndcapRings this classs is used for both OT and Pixel detector

Definition at line 15 of file Phase2EndcapLayer.h.

Constructor & Destructor Documentation

Phase2EndcapLayer::Phase2EndcapLayer	(	std::vector< const Phase2EndcapRing * > &	rings,
		const bool	isOT
	)

Definition at line 51 of file Phase2EndcapLayer.cc.

References basicComponents(), computeDisk(), fillRingPars(), mps_fire::i, LogDebug, position, theBasicComps, theComps, and theRingSize.

                                                                                           :
  RingedForwardLayer(true),
  isOuterTracker(isOT),
  theComponents{nullptr}
{
  //They should be already R-ordered. TO BE CHECKED!!
  //sort( theRings.begin(), theRings.end(), DetLessR());

  theRingSize = rings.size();
  LogDebug("TkDetLayers") << "Number of rings in Phase2 OT EC layer is " << theRingSize << std::endl;
  setSurface( computeDisk( rings ) );

  for(unsigned int i=0; i!=rings.size(); ++i) {
    theComps.push_back(rings[i]);
    fillRingPars(i);
    theBasicComps.insert(theBasicComps.end(),   
             (*rings[i]).basicComponents().begin(),
             (*rings[i]).basicComponents().end());
  }

 
  LogDebug("TkDetLayers") << "==== DEBUG Phase2EndcapLayer =====" ; 
  LogDebug("TkDetLayers") << "r,zed pos  , thickness, innerR, outerR: " 
              << this->position().perp() << " , "
              << this->position().z() << " , "
              << this->specificSurface().bounds().thickness() << " , "
              << this->specificSurface().innerRadius() << " , "
              << this->specificSurface().outerRadius() ;
}

Phase2EndcapLayer::~Phase2EndcapLayer ( )

override

Definition at line 113 of file Phase2EndcapLayer.cc.

References EnergyCorrector::c, theComponents, and theComps.

                                     {
  for (auto c : theComps) delete c;

  delete theComponents.load();
} 

Phase2EndcapLayer::Phase2EndcapLayer ( const Phase2EndcapLayer &  )

delete

Member Function Documentation

const std::vector<const GeomDet*>& Phase2EndcapLayer::basicComponents ( ) const

inlineoverride

Definition at line 26 of file Phase2EndcapLayer.h.

References components(), groupedCompatibleDetsV(), mps_fire::result, and theBasicComps.

Referenced by Phase2EndcapLayer().

{return theBasicComps;}

const std::vector< const GeometricSearchDet * > & Phase2EndcapLayer::components ( ) const

override

Definition at line 22 of file Phase2EndcapLayer.cc.

References EnergyCorrector::c, and groupFilesInBlocks::temp.

Referenced by basicComponents().

                                                                               {
  if (not theComponents) {
    auto temp = std::make_unique<std::vector<const GeometricSearchDet*>>();
    temp->reserve(15);   // This number is just an upper bound
    for ( auto c: theComps) temp->push_back(c);
    std::vector<const GeometricSearchDet*>* expected = nullptr;
    if(theComponents.compare_exchange_strong(expected,temp.get())) {
      //this thread set the value
      temp.release();
    }
  }
  return *theComponents;
 }

BoundDisk * Phase2EndcapLayer::computeDisk ( const std::vector< const Phase2EndcapRing * > &  rings ) const

private

Definition at line 83 of file Phase2EndcapLayer.cc.

References mps_fire::i, SiStripPI::max, min(), and makeMuonMisalignmentScenario::rot.

Referenced by Phase2EndcapLayer(), and subDetector().

{
  float theRmin = rings.front()->specificSurface().innerRadius();
  float theRmax = rings.front()->specificSurface().outerRadius();
  float theZmin = rings.front()->position().z() -
    rings.front()->surface().bounds().thickness()/2;
  float theZmax = rings.front()->position().z() +
    rings.front()->surface().bounds().thickness()/2;
  
  for (vector<const Phase2EndcapRing*>::const_iterator i = rings.begin(); i != rings.end(); i++) {
    float rmin = (**i).specificSurface().innerRadius();
    float rmax = (**i).specificSurface().outerRadius();
    float zmin = (**i).position().z() - (**i).surface().bounds().thickness()/2.;
    float zmax = (**i).position().z() + (**i).surface().bounds().thickness()/2.;
    theRmin = min( theRmin, rmin);
    theRmax = max( theRmax, rmax);
    theZmin = min( theZmin, zmin);
    theZmax = max( theZmax, zmax);
  }
  
  float zPos = (theZmax+theZmin)/2.;
  PositionType pos(0.,0.,zPos);
  RotationType rot;

  return new BoundDisk( pos, rot, new SimpleDiskBounds(theRmin, theRmax,    
                      theZmin-zPos, theZmax-zPos));

}

float Phase2EndcapLayer::computeWindowSize ( const GeomDet *  det, const TrajectoryStateOnSurface &  tsos, const MeasurementEstimator &  est  ) const

private

Definition at line 298 of file Phase2EndcapLayer.cc.

References particleFlowClusterHGC_cfi::maxDistance, MeasurementEstimator::maximalLocalDisplacement(), GeomDet::surface(), and PV2DBase< T, PVType, FrameType >::y().

Referenced by groupedCompatibleDetsV(), and subDetector().

{
  const Plane& startPlane = det->surface();  
  MeasurementEstimator::Local2DVector maxDistance = 
    est.maximalLocalDisplacement( tsos, startPlane);
  return maxDistance.y();
}

void Phase2EndcapLayer::fillRingPars ( int  i )

private

Definition at line 38 of file Phase2EndcapLayer.cc.

References funct::abs(), mps_fire::i, Phase2EndcapLayer::RingPar::theRingR, Phase2EndcapLayer::RingPar::thetaRingMax, and Phase2EndcapLayer::RingPar::thetaRingMin.

Referenced by Phase2EndcapLayer(), and subDetector().

                                     {
  const BoundDisk& ringDisk = static_cast<const BoundDisk&>(theComps[i]->surface());
  float ringMinZ = std::abs( ringDisk.position().z()) - ringDisk.bounds().thickness()/2.;
  float ringMaxZ = std::abs( ringDisk.position().z()) + ringDisk.bounds().thickness()/2.; 
  RingPar tempPar; 
  tempPar.thetaRingMin =  ringDisk.innerRadius()/ ringMaxZ;
  tempPar.thetaRingMax =  ringDisk.outerRadius()/ ringMinZ;
  tempPar.theRingR=( ringDisk.innerRadius() +
             ringDisk.outerRadius())/2.;
  ringPars.push_back(tempPar); 
}

std::array< int, 3 > Phase2EndcapLayer::findThreeClosest ( std::vector< GlobalPoint >  ringCrossing ) const

private

Definition at line 309 of file Phase2EndcapLayer.cc.

References funct::abs(), mps_fire::i, perp(), ringPars, and theRingSize.

Referenced by ringIndicesByCrossingProximity(), and subDetector().

{
  std::array<int,3> theBins={{-1,-1,-1}};
  theBins[0] = 0;
  float initialR =  ringPars[0].theRingR;
  float rDiff0 = std::abs( ringCrossing[0].perp() - initialR);
  float rDiff1 = -1.;
  float rDiff2 = -1.;
  for (int i = 1; i < theRingSize ; i++){
    float ringR =  ringPars[i].theRingR;
    float testDiff = std::abs( ringCrossing[i].perp() - ringR);
    if ( testDiff<rDiff0 ) {
      rDiff2 = rDiff1;
      rDiff1 = rDiff0;
      rDiff0 = testDiff;
      theBins[2] = theBins[1];
      theBins[1] = theBins[0];
      theBins[0] = i;
    } 
    else 
    if ( rDiff1 < 0 || testDiff<rDiff1 ) {
      rDiff2 = rDiff1;
      rDiff1 = testDiff;
      theBins[2] = theBins[1];
      theBins[1] = i;
    } 
    else 
    if ( rDiff2 < 0 || testDiff<rDiff2 ) {
      rDiff2 = testDiff;
      theBins[2] = i;
    } 
  }

  return theBins;
}

void Phase2EndcapLayer::groupedCompatibleDetsV ( const TrajectoryStateOnSurface &  tsos, const Propagator &  prop, const MeasurementEstimator &  est, std::vector< DetGroup > &  result  ) const

override

Definition at line 122 of file Phase2EndcapLayer.cc.

References DetGroupMerger::addSameLevel(), alongMomentum, computeWindowSize(), lumiContext::fill, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, diffTreeTool::index, eostools::move(), DetGroupMerger::orderAndMergeTwoLevels(), overlapInR(), position, Propagator::propagationDirection(), ringIndicesByCrossingProximity(), theComps, theRingSize, z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by basicComponents().

{
  std::array<int,3> const & ringIndices = ringIndicesByCrossingProximity(startingState,prop);
  if ( ringIndices[0]==-1 || ringIndices[1] ==-1 || ringIndices[2] == -1 ) {
    edm::LogError("TkDetLayers") << "TkRingedForwardLayer::groupedCompatibleDets : error in CrossingProximity";
    return;
  }

  //order is: rings in front = 0; rings in back = 1
  //rings should be already ordered in r
  //if the layer has 1 ring, it does not matter
  //FIXME: to be optimized once the geometry is stable
  std::vector<int> ringOrder(theRingSize);
  std::fill(ringOrder.begin(), ringOrder.end(), 1);
  if (theRingSize > 1){
    if(fabs(theComps[0]->position().z()) < fabs(theComps[1]->position().z())){
      for (int i=0; i<theRingSize; i++) {
        if(i % 2 == 0) ringOrder[i] = 0;
      }
    } else if(fabs(theComps[0]->position().z()) > fabs(theComps[1]->position().z())){
      std::fill(ringOrder.begin(), ringOrder.end(), 0);
      for (int i=0; i<theRingSize; i++) {
        if(i % 2 == 0) ringOrder[i] = 1;
      }
    } else {
      throw DetLayerException("Rings in Endcap Layer have same z position, no idea how to order them!");
    }
  }

  auto index = [&ringIndices,& ringOrder](int i) { return ringOrder[ringIndices[i]];};

  std::vector<DetGroup> closestResult;
  theComps[ringIndices[0]]->groupedCompatibleDetsV( startingState, prop, est, closestResult);       
  // if the closest is empty, use the next one and exit: inherited from TID !
  if ( closestResult.empty() ){
    theComps[ringIndices[1]]->groupedCompatibleDetsV( startingState, prop, est, result); 
    return;
  }

  DetGroupElement closestGel( closestResult.front().front());  
  float rWindow = computeWindowSize( closestGel.det(), closestGel.trajectoryState(), est); 

  // check if next ring and next next ring are found and if there is overlap

  bool ring1ok = ringIndices[1] != -1 && overlapInR(closestGel.trajectoryState(),ringIndices[1],rWindow);
  bool ring2ok = ringIndices[2] != -1 && overlapInR(closestGel.trajectoryState(),ringIndices[2],rWindow);

  // look for the two rings in the same plane (are they only two?)

  // determine if we are propagating from in to out (0) or from out to in (1)

  int direction = 0;
  if(startingState.globalPosition().z()*startingState.globalMomentum().z()>0) {
    if(prop.propagationDirection() == alongMomentum) direction=0;
    else direction=1;
  }
  else{
    if(prop.propagationDirection() == alongMomentum) direction=1;
    else direction=0;
  }

  if((index(0) == index(1)) && (index(0) == index(2))) {
    edm::LogInfo("AllRingsInOnePlane") << " All rings: " 
                       << ringIndices[0] << " " 
                       << ringIndices[1] << " " 
                       << ringIndices[2] << " in one plane. Only the first two will be considered";
    ring2ok=false;
  }

  if(index(0) == index(1)) {
    if(ring1ok) {
      std::vector<DetGroup> ring1res;
      theComps[ringIndices[1]]->groupedCompatibleDetsV( startingState, prop, est, ring1res);
      DetGroupMerger::addSameLevel(std::move(ring1res),closestResult);
    }
    if(ring2ok) {
      std::vector<DetGroup> ring2res;
      theComps[ringIndices[2]]->groupedCompatibleDetsV( startingState, prop, est, ring2res);
      DetGroupMerger::orderAndMergeTwoLevels(std::move(closestResult),std::move(ring2res),result,index(0),direction);
      return;
    }
    else {
      result.swap(closestResult);
      return;
    }
  }
  else if(index(0) == index(2)) {
    if(ring2ok) {
      std::vector<DetGroup> ring2res;
      theComps[ringIndices[2]]->groupedCompatibleDetsV( startingState, prop, est, ring2res);
      DetGroupMerger::addSameLevel(std::move(ring2res),closestResult);
    }
    if(ring1ok) {
      std::vector<DetGroup> ring1res;
      theComps[ringIndices[1]]->groupedCompatibleDetsV( startingState, prop, est, ring1res);
      DetGroupMerger::orderAndMergeTwoLevels(std::move(closestResult),std::move(ring1res),result,index(0),direction);
      return;
    }
    else {
      result.swap(closestResult);
      return;
    }
  }
  else {
    std::vector<DetGroup> ring12res;
    if(ring1ok) {
      std::vector<DetGroup> ring1res;
      theComps[ringIndices[1]]->groupedCompatibleDetsV( startingState, prop, est, ring1res);
      ring12res.swap(ring1res);
    } 
    if(ring2ok) {
      std::vector<DetGroup> ring2res;
      theComps[ringIndices[2]]->groupedCompatibleDetsV( startingState, prop, est, ring2res);
      DetGroupMerger::addSameLevel(std::move(ring2res),ring12res);
    }
    if(!ring12res.empty()) {
      DetGroupMerger::orderAndMergeTwoLevels(std::move(closestResult),std::move(ring12res),result,index(0),direction);
      return;
    }
    else {
      result.swap(closestResult);
      return;
    }
  }
}

Phase2EndcapLayer& Phase2EndcapLayer::operator= ( const Phase2EndcapLayer &  )

delete

bool Phase2EndcapLayer::overlapInR ( const TrajectoryStateOnSurface &  tsos, int  i, double  ymax  ) const

private

Definition at line 346 of file Phase2EndcapLayer.cc.

References funct::abs(), f, TrajectoryStateOnSurface::globalPosition(), diffTreeTool::index, SiStripPI::max, PV3DBase< T, PVType, FrameType >::perp(), ringPars, Phase2TrackerMonitorDigi_cff::ymax, and PV3DBase< T, PVType, FrameType >::z().

Referenced by groupedCompatibleDetsV(), and subDetector().

{
  // assume "fixed theta window", i.e. margin in local y = r is changing linearly with z
  float tsRadius = tsos.globalPosition().perp();
  float thetamin = ( max(0.,tsRadius-ymax))/(std::abs(tsos.globalPosition().z())+10.f); // add 10 cm contingency 
  float thetamax = ( tsRadius + ymax)/(std::abs(tsos.globalPosition().z())-10.f);
  
  // do the theta regions overlap ?

  return !( thetamin > ringPars[index].thetaRingMax || ringPars[index].thetaRingMin > thetamax);
}

std::array< int, 3 > Phase2EndcapLayer::ringIndicesByCrossingProximity ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop  ) const

private

Definition at line 253 of file Phase2EndcapLayer.cc.

References findThreeClosest(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, Propagator::propagationDirection(), rho, theComps, theRingSize, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by groupedCompatibleDetsV(), and subDetector().

{
  typedef HelixForwardPlaneCrossing Crossing; 
  typedef MeasurementEstimator::Local2DVector Local2DVector;

  HelixPlaneCrossing::PositionType startPos( startingState.globalPosition());
  HelixPlaneCrossing::DirectionType startDir( startingState.globalMomentum());
  PropagationDirection propDir( prop.propagationDirection());
  float rho( startingState.transverseCurvature());

  // calculate the crossings with the ring surfaces
  // rings are assumed to be sorted in R !
  
  Crossing myXing(  startPos, startDir, rho, propDir );

  std::vector<GlobalPoint> ringCrossings;
  ringCrossings.reserve(theRingSize);
  // vector<GlobalVector>  ringXDirections;

  for (int i = 0; i < theRingSize ; i++ ) {
    const BoundDisk & theRing  = static_cast<const BoundDisk &>(theComps[i]->surface());
    pair<bool,double> pathlen = myXing.pathLength( theRing);
    if ( pathlen.first ) { 
      ringCrossings.push_back(GlobalPoint( myXing.position(pathlen.second )));
      // ringXDirections.push_back( GlobalVector( myXing.direction(pathlen.second )));
    } else {
      // TO FIX.... perhaps there is something smarter to do
      //throw DetLayerException("trajectory doesn't cross TID rings");
      ringCrossings.push_back(GlobalPoint( 0.,0.,0.));
      //  ringXDirections.push_back( GlobalVector( 0.,0.,0.));
    }
  }

  //find three closest rings to the crossing

  std::array<int,3> closests = findThreeClosest(ringCrossings);

  return closests;
}

SubDetector Phase2EndcapLayer::subDetector ( ) const

inlineoverride

Definition at line 36 of file Phase2EndcapLayer.h.

References computeDisk(), computeWindowSize(), fillRingPars(), findThreeClosest(), mps_fire::i, isOuterTracker, overlapInR(), GeomDetEnumerators::P2OTEC, GeomDetEnumerators::P2PXEC, ringIndicesByCrossingProximity(), GeomDetEnumerators::subDetGeom, and Phase2TrackerMonitorDigi_cff::ymax.

                                           { if(isOuterTracker) return GeomDetEnumerators::subDetGeom[GeomDetEnumerators::P2OTEC];
                                            else return GeomDetEnumerators::subDetGeom[GeomDetEnumerators::P2PXEC];}

Member Data Documentation

const bool Phase2EndcapLayer::isOuterTracker

private

Definition at line 63 of file Phase2EndcapLayer.h.

Referenced by subDetector().

std::vector<RingPar> Phase2EndcapLayer::ringPars

private

Definition at line 67 of file Phase2EndcapLayer.h.

Referenced by findThreeClosest(), and overlapInR().

std::vector<GeomDet const*> Phase2EndcapLayer::theBasicComps

private

Definition at line 62 of file Phase2EndcapLayer.h.

Referenced by basicComponents(), and Phase2EndcapLayer().

std::atomic<std::vector<const GeometricSearchDet*>*> Phase2EndcapLayer::theComponents

mutableprivate

Definition at line 64 of file Phase2EndcapLayer.h.

Referenced by ~Phase2EndcapLayer().

std::vector<const Phase2EndcapRing*> Phase2EndcapLayer::theComps

private

Definition at line 65 of file Phase2EndcapLayer.h.

Referenced by groupedCompatibleDetsV(), Phase2EndcapLayer(), ringIndicesByCrossingProximity(), and ~Phase2EndcapLayer().

int Phase2EndcapLayer::theRingSize

private

Definition at line 68 of file Phase2EndcapLayer.h.

Referenced by findThreeClosest(), groupedCompatibleDetsV(), Phase2EndcapLayer(), and ringIndicesByCrossingProximity().