TTStubAlgorithm_official.cc

Go to the documentation of this file.

#include "L1Trigger/TrackTrigger/interface/TTStubAlgorithm_official.h"

template< >
void TTStubAlgorithm_official< Ref_Phase2TrackerDigi_ >::PatternHitCorrelation( bool &aConfirmation,
                                                                       int &aDisplacement, 
                                                                       int &anOffset, 
                                                                       const TTStub< Ref_Phase2TrackerDigi_ > &aTTStub ) const
{ 
  // These are already corrected for being at the center of each pixel
  MeasurementPoint mp0 = aTTStub.getClusterRef(0)->findAverageLocalCoordinates();
  MeasurementPoint mp1 = aTTStub.getClusterRef(1)->findAverageLocalCoordinates();

  bool isPS = (theTrackerGeom_->getDetectorType(aTTStub.getDetId())==TrackerGeometry::ModuleType::Ph2PSP);
  // TODO temporary: should use a method from the topology
  DetId stDetId( aTTStub.getDetId() );
  const GeomDetUnit* det0 = theTrackerGeom_->idToDetUnit( stDetId+1 );
  const GeomDetUnit* det1 = theTrackerGeom_->idToDetUnit( stDetId+2 );

  const PixelGeomDetUnit* pix0 = dynamic_cast< const PixelGeomDetUnit* >( det0 );
  const PixelGeomDetUnit* pix1 = dynamic_cast< const PixelGeomDetUnit* >( det1 );
  const PixelTopology* top0 = dynamic_cast< const PixelTopology* >( &(pix0->specificTopology()) );
  const PixelTopology* top1 = dynamic_cast< const PixelTopology* >( &(pix1->specificTopology()) );
  std::pair< float, float > pitch0 = top0->pitch();
  std::pair< float, float > pitch1 = top1->pitch();

  int cols0 = top0->ncolumns();
  int cols1 = top1->ncolumns();
  int ratio = cols0/cols1; 
  int segment0 = floor( mp0.y() / ratio );

//  if ( ratio == 1 ) /// 2S Modules
  if (!isPS)
  {
    if ( mPerformZMatching2S && ( segment0 != floor( mp1.y() ) ) )
      return;
  }
  else 
  {
    if ( mPerformZMatchingPS && ( segment0 != floor( mp1.y() ) ) )
      return;
  }

  double R0 = det0->position().perp();
  double R1 = det1->position().perp();
  double Z0 = det0->position().z();
  double Z1 = det1->position().z();

  double DR = R1-R0;
  double DZ = Z1-Z0;

  double alpha = atan2(DR,DZ);
  double delta = sqrt(DR*DR+DZ*DZ)/(R0*sin(alpha)+Z0*cos(alpha));

  int window=0;



  double dispD = 2 * (mp1.x() - mp0.x()) * (pitch0.first / pitch1.first); 
  int dispI = ((dispD>0)-(dispD<0))*floor(std::abs(dispD)); 
  
  double offsetD = 2 * delta * ( mp0.x() - (top0->nrows()/2 - 0.5) ) * (pitch0.first / pitch1.first); 
  int offsetI = ((offsetD>0)-(offsetD<0))*floor(std::abs(offsetD)); 

  if (stDetId.subdetId()==StripSubdetector::TOB)
  {
    int layer  = theTrackerTopo_->layer(stDetId);
    int ladder = theTrackerTopo_->tobRod(stDetId);
    int type   = 2*theTrackerTopo_->tobSide(stDetId)-3; // -1 for tilted-, 1 for tilted+, 3 for flat
    double corr=0;

    if (type<3) // Only for tilted modules
    {
      corr   = (barrelNTilt.at(layer)+1)/2.;
      ladder = corr-(corr-ladder)*type; // Corrected ring number, bet 0 and barrelNTilt.at(layer), in ascending |z|
      window = 2*(tiltedCut.at(layer)).at(ladder);
    }
    else // Classis barrel window otherwise
    {
      window = 2*barrelCut.at( layer );
    }
 
  }
  else if (stDetId.subdetId()==StripSubdetector::TID)
  {
    window = 2*(ringCut.at( theTrackerTopo_->tidWheel(stDetId))).at(theTrackerTopo_->tidRing(stDetId));
  }

  if ( std::abs(dispI - offsetI) <= window ) 
  {
    aConfirmation = true;
    aDisplacement = dispI; 
    anOffset = offsetI; 
  } 

}