TTStubAlgorithm_cbc3.cc

Go to the documentation of this file.

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

template <>
void TTStubAlgorithm_cbc3<Ref_Phase2TrackerDigi_>::PatternHitCorrelation(
    bool &aConfirmation,
    int &aDisplacement,
    int &anOffset,
    float &anHardBend,
    const TTStub<Ref_Phase2TrackerDigi_> &aTTStub) const

{
  /*
  MeasurementPoint mp0 = aTTStub.getClusterRef(0)->findAverageLocalCoordinates();
  MeasurementPoint mp1 = aTTStub.getClusterRef(1)->findAverageLocalCoordinates();

  StackedTrackerDetId stDetId( aTTStub.getDetId() );

  //bool isPS = TTStubAlgorithm< Ref_Phase2TrackerDigi_ >::theStackedTracker->isPSModule( stDetId );

  const GeomDetUnit* det0 = TTStubAlgorithm< Ref_Phase2TrackerDigi_ >::theStackedTracker->idToDetUnit( stDetId, 0 );

  const PixelGeomDetUnit* pix0 = dynamic_cast< const PixelGeomDetUnit* >( det0 );
  const PixelTopology* top0 = dynamic_cast< const PixelTopology* >( &(pix0->specificTopology()) );

  int myPosition = mp0.x();

  int chipSize = top0->rowsperroc();
  int asicNumber = myPosition / chipSize; /// ASIC in module
  int partitionSize = ceil( float(chipSize) / float( TTStubAlgorithm< Ref_Phase2TrackerDigi_ >::theStackedTracker->getPartitionsPerRoc() ) );
  int partitionNumber = (myPosition % chipSize) / partitionSize; /// Partition in ASIC

  anOffset = TTStubAlgorithm< Ref_Phase2TrackerDigi_ >::theStackedTracker->getASICOffset( stDetId, asicNumber, partitionNumber );

  int aPosition = 2 * mp0.x();
  int aBend = (2 * mp1.x()) - aPosition - anOffset;

  aDisplacement = (2 * mp1.x()) - aPosition;

  aConfirmation =  ( (abs(4 * aBend - 1)) <= (2 * TTStubAlgorithm< Ref_Phase2TrackerDigi_ >::theStackedTracker->getDetUnitWindow(stDetId) ) );

  if ( !mPerformZMatching2S ) // && !isPS
    return;

  //if ( !mPerformZMatchingPS && isPS )
  //  return;

  const GeomDetUnit* det1 = TTStubAlgorithm< Ref_Phase2TrackerDigi_ >::theStackedTracker->idToDetUnit( stDetId, 1 );
  const PixelGeomDetUnit* pix1 = dynamic_cast< const PixelGeomDetUnit* >( det1 );
  const PixelTopology* top1 = dynamic_cast< const PixelTopology* >( &(pix1->specificTopology()) );
  int cols0 = top0->ncolumns();
  int cols1 = top1->ncolumns();
  int ratio = cols0/cols1; /// This assumes the ratio is integer!
  int segment0 = floor( mp0.y() / ratio );
  if ( segment0 != floor( mp1.y() ) )
    aConfirmation = false;
  */
}