VectorHitBuilderAlgorithmBase.cc

Go to the documentation of this file.

#include "RecoLocalTracker/SiPhase2VectorHitBuilder/interface/VectorHitBuilderAlgorithmBase.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "RecoLocalTracker/Records/interface/TkPhase2OTCPERecord.h"

#include "Geometry/CommonTopologies/interface/PixelTopology.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

VectorHitBuilderAlgorithmBase::VectorHitBuilderAlgorithmBase(
    const edm::ParameterSet& conf,
    const TrackerGeometry* tkGeomProd,
    const TrackerTopology* tkTopoProd,
    const ClusterParameterEstimator<Phase2TrackerCluster1D>* cpeProd)
    : tkGeom_(tkGeomProd),
      tkTopo_(tkTopoProd),
      cpe_(cpeProd),
      nMaxVHforeachStack_(conf.getParameter<int>("maxVectorHitsInAStack")),
      barrelCut_(conf.getParameter<std::vector<double> >("BarrelCut")),
      endcapCut_(conf.getParameter<std::vector<double> >("EndcapCut")),
      cpeTag_(conf.getParameter<edm::ESInputTag>("CPE")) {}

double VectorHitBuilderAlgorithmBase::computeParallaxCorrection(const PixelGeomDetUnit* geomDetUnit_low,
                                                                const Point3DBase<float, LocalTag>& lPosClu_low,
                                                                const PixelGeomDetUnit* geomDetUnit_upp,
                                                                const Point3DBase<float, LocalTag>& lPosClu_upp) const {
  double parallCorr = 0.0;
  Global3DPoint origin(0, 0, 0);
  Global3DPoint gPosClu_low = geomDetUnit_low->surface().toGlobal(lPosClu_low);
  GlobalVector gV = gPosClu_low - origin;
  LogTrace("VectorHitsBuilderValidation") << " global vector passing to the origin:" << gV;

  LocalVector lV = geomDetUnit_low->surface().toLocal(gV);
  LogTrace("VectorHitsBuilderValidation")
      << " local vector passing to the origin (in the lower detector system of reference):" << lV;
  LocalVector lV_norm = lV / lV.z();
  LogTrace("VectorHitsBuilderValidation")
      << " normalized local vector passing to the origin (in low the lower detector system of reference):" << lV_norm;

  Global3DPoint gPosClu_upp = geomDetUnit_upp->surface().toGlobal(lPosClu_upp);
  Local3DPoint lPosClu_uppInLow = geomDetUnit_low->surface().toLocal(gPosClu_upp);
  parallCorr = lV_norm.x() * lPosClu_uppInLow.z();

  return parallCorr;
}

void VectorHitBuilderAlgorithmBase::printClusters(const edmNew::DetSetVector<Phase2TrackerCluster1D>& clusters) const {
  int nCluster = 0;
  for (const auto& DSViter : clusters) {
    // Loop over the clusters in the detector unit
    for (const auto& clustIt : DSViter) {
      nCluster++;
      // get the detector unit's id
      const GeomDetUnit* geomDetUnit(tkGeom_->idToDetUnit(DSViter.detId()));
      if (!geomDetUnit)
        return;
      printCluster(geomDetUnit, &clustIt);
    }
  }
  LogDebug("VectorHitBuilder") << " Number of input clusters: " << nCluster << std::endl;
}

void VectorHitBuilderAlgorithmBase::printCluster(const GeomDet* geomDetUnit,
                                                 const Phase2TrackerCluster1D* clustIt) const {
  if (!geomDetUnit)
    return;
  const PixelGeomDetUnit* pixelGeomDetUnit = dynamic_cast<const PixelGeomDetUnit*>(geomDetUnit);
  const PixelTopology& topol = pixelGeomDetUnit->specificTopology();
  if (!pixelGeomDetUnit)
    return;

  unsigned int layer = tkTopo_->layer(geomDetUnit->geographicalId());
  unsigned int module = tkTopo_->module(geomDetUnit->geographicalId());
  LogTrace("VectorHitBuilder") << "Layer:" << layer << " and DetId: " << geomDetUnit->geographicalId().rawId()
                               << std::endl;
  TrackerGeometry::ModuleType mType = tkGeom_->getDetectorType(geomDetUnit->geographicalId());
  if (mType == TrackerGeometry::ModuleType::Ph2PSP)
    LogTrace("VectorHitBuilder") << "Pixel cluster (module:" << module << ") " << std::endl;
  else if (mType == TrackerGeometry::ModuleType::Ph2SS || mType == TrackerGeometry::ModuleType::Ph2PSS)
    LogTrace("VectorHitBuilder") << "Strip cluster (module:" << module << ") " << std::endl;
  else
    LogTrace("VectorHitBuilder") << "no module?!" << std::endl;
  LogTrace("VectorHitBuilder") << "with pitch:" << topol.pitch().first << " , " << topol.pitch().second << std::endl;
  LogTrace("VectorHitBuilder") << " and width:" << pixelGeomDetUnit->surface().bounds().width()
                               << " , lenght:" << pixelGeomDetUnit->surface().bounds().length() << std::endl;

  auto&& lparams = cpe_->localParameters(*clustIt, *pixelGeomDetUnit);
  Global3DPoint gparams = pixelGeomDetUnit->surface().toGlobal(lparams.first);

  LogTrace("VectorHitBuilder") << "\t global pos " << gparams << std::endl;
  LogTrace("VectorHitBuilder") << "\t local  pos " << lparams.first << "with err " << lparams.second << std::endl;
  LogTrace("VectorHitBuilder") << std::endl;

  return;
}