PixelCPEGenericBase.cc

Go to the documentation of this file.

#include "RecoLocalTracker/SiPixelRecHits/interface/PixelCPEGenericBase.h"
#include "Geometry/TrackerGeometryBuilder/interface/RectangularPixelTopology.h"

namespace {
  constexpr float micronsToCm = 1.0e-4;
  const auto convertDoubleVecToFloatVec = [](std::vector<double> const& iIn) {
    return std::vector<float>(iIn.begin(), iIn.end());
  };
}  // namespace

PixelCPEGenericBase::PixelCPEGenericBase(edm::ParameterSet const& conf,
                                         const MagneticField* mag,
                                         const TrackerGeometry& geom,
                                         const TrackerTopology& ttopo,
                                         const SiPixelLorentzAngle* lorentzAngle,
                                         const SiPixelGenErrorDBObject* genErrorDBObject,
                                         const SiPixelLorentzAngle* lorentzAngleWidth = nullptr)
    : PixelCPEBase(conf, mag, geom, ttopo, lorentzAngle, genErrorDBObject, nullptr, lorentzAngleWidth, 0),
      edgeClusterErrorX_{static_cast<float>(conf.getParameter<double>("EdgeClusterErrorX"))},
      edgeClusterErrorY_{static_cast<float>(conf.getParameter<double>("EdgeClusterErrorY"))},
      useErrorsFromTemplates_{conf.getParameter<bool>("UseErrorsFromTemplates")},
      truncatePixelCharge_{conf.getParameter<bool>("TruncatePixelCharge")},
      xerr_barrel_l1_{convertDoubleVecToFloatVec(conf.getParameter<std::vector<double>>("xerr_barrel_l1"))},
      yerr_barrel_l1_{convertDoubleVecToFloatVec(conf.getParameter<std::vector<double>>("yerr_barrel_l1"))},
      xerr_barrel_ln_{convertDoubleVecToFloatVec(conf.getParameter<std::vector<double>>("xerr_barrel_ln"))},
      yerr_barrel_ln_{convertDoubleVecToFloatVec(conf.getParameter<std::vector<double>>("yerr_barrel_ln"))},
      xerr_endcap_{convertDoubleVecToFloatVec(conf.getParameter<std::vector<double>>("xerr_endcap"))},
      yerr_endcap_{convertDoubleVecToFloatVec(conf.getParameter<std::vector<double>>("yerr_endcap"))},
      xerr_barrel_l1_def_{static_cast<float>(conf.getParameter<double>("xerr_barrel_l1_def"))},
      yerr_barrel_l1_def_{static_cast<float>(conf.getParameter<double>("yerr_barrel_l1_def"))},
      xerr_barrel_ln_def_{static_cast<float>(conf.getParameter<double>("xerr_barrel_ln_def"))},
      yerr_barrel_ln_def_{static_cast<float>(conf.getParameter<double>("yerr_barrel_ln_def"))},
      xerr_endcap_def_{static_cast<float>(conf.getParameter<double>("xerr_endcap_def"))},
      yerr_endcap_def_{static_cast<float>(conf.getParameter<double>("yerr_endcap_def"))} {};

std::unique_ptr<PixelCPEBase::ClusterParam> PixelCPEGenericBase::createClusterParam(const SiPixelCluster& cl) const {
  return std::make_unique<ClusterParamGeneric>(cl);
}

//-----------------------------------------------------------------------------
//-----------------------------------------------------------------------------
void PixelCPEGenericBase::collect_edge_charges(ClusterParam& theClusterParamBase,  
                                               int& q_f_X,                         
                                               int& q_l_X,                         
                                               int& q_f_Y,                         
                                               int& q_l_Y,                         
                                               bool truncate) {
  ClusterParamGeneric& theClusterParam = static_cast<ClusterParamGeneric&>(theClusterParamBase);

  // Initialize return variables.
  q_f_X = q_l_X = 0;
  q_f_Y = q_l_Y = 0;

  // Obtain boundaries in index units
  int xmin = theClusterParam.theCluster->minPixelRow();
  int xmax = theClusterParam.theCluster->maxPixelRow();
  int ymin = theClusterParam.theCluster->minPixelCol();
  int ymax = theClusterParam.theCluster->maxPixelCol();

  // Iterate over the pixels.
  int isize = theClusterParam.theCluster->size();
  for (int i = 0; i != isize; ++i) {
    auto const& pixel = theClusterParam.theCluster->pixel(i);
    // ggiurgiu@fnal.gov: add pixel charge truncation
    int pix_adc = pixel.adc;
    if (truncate)
      pix_adc = std::min(pix_adc, theClusterParam.pixmx);

    //
    // X projection
    if (pixel.x == xmin)
      q_f_X += pix_adc;
    if (pixel.x == xmax)
      q_l_X += pix_adc;
    //
    // Y projection
    if (pixel.y == ymin)
      q_f_Y += pix_adc;
    if (pixel.y == ymax)
      q_l_Y += pix_adc;
  }
}

void PixelCPEGenericBase::initializeLocalErrorVariables(
    float& xerr,
    float& yerr,
    bool& edgex,
    bool& edgey,
    bool& bigInX,
    bool& bigInY,
    int& maxPixelCol,
    int& maxPixelRow,
    int& minPixelCol,
    int& minPixelRow,
    uint& sizex,
    uint& sizey,
    DetParam const& theDetParam,
    ClusterParamGeneric const& theClusterParam) const {  // Default errors are the maximum error used for edge clusters.
  // These are determined by looking at residuals for edge clusters
  xerr = edgeClusterErrorX_ * micronsToCm;
  yerr = edgeClusterErrorY_ * micronsToCm;

  // Find if cluster is at the module edge.
  maxPixelCol = theClusterParam.theCluster->maxPixelCol();
  maxPixelRow = theClusterParam.theCluster->maxPixelRow();
  minPixelCol = theClusterParam.theCluster->minPixelCol();
  minPixelRow = theClusterParam.theCluster->minPixelRow();

  edgex = (theDetParam.theRecTopol->isItEdgePixelInX(minPixelRow)) ||
          (theDetParam.theRecTopol->isItEdgePixelInX(maxPixelRow));
  edgey = (theDetParam.theRecTopol->isItEdgePixelInY(minPixelCol)) ||
          (theDetParam.theRecTopol->isItEdgePixelInY(maxPixelCol));

  sizex = theClusterParam.theCluster->sizeX();
  sizey = theClusterParam.theCluster->sizeY();

  // Find if cluster contains double (big) pixels.
  bigInX = theDetParam.theRecTopol->containsBigPixelInX(minPixelRow, maxPixelRow);
  bigInY = theDetParam.theRecTopol->containsBigPixelInY(minPixelCol, maxPixelCol);
};

void PixelCPEGenericBase::setXYErrors(float& xerr,
                                      float& yerr,
                                      const bool edgex,
                                      const bool edgey,
                                      const unsigned int sizex,
                                      const unsigned int sizey,
                                      const bool bigInX,
                                      const bool bigInY,
                                      const bool useTemplateErrors,
                                      DetParam const& theDetParam,
                                      ClusterParamGeneric const& theClusterParam) const {
  if (useTemplateErrors) {
    if (!edgex) {  // Only use this for non-edge clusters
      if (sizex == 1) {
        if (!bigInX) {
          xerr = theClusterParam.sx1;
        } else {
          xerr = theClusterParam.sx2;
        }
      } else {
        xerr = theClusterParam.sigmax;
      }
    }

    if (!edgey) {  // Only use for non-edge clusters
      if (sizey == 1) {
        if (!bigInY) {
          yerr = theClusterParam.sy1;
        } else {
          yerr = theClusterParam.sy2;
        }
      } else {
        yerr = theClusterParam.sigmay;
      }
    }

  } else {  // simple errors

    if (GeomDetEnumerators::isTrackerPixel(theDetParam.thePart)) {
      if (GeomDetEnumerators::isBarrel(theDetParam.thePart)) {
        DetId id = (theDetParam.theDet->geographicalId());
        int layer = ttopo_.layer(id);
        if (layer == 1) {
          if (!edgex) {
            if (sizex <= xerr_barrel_l1_.size())
              xerr = xerr_barrel_l1_[sizex - 1];
            else
              xerr = xerr_barrel_l1_def_;
          }

          if (!edgey) {
            if (sizey <= yerr_barrel_l1_.size())
              yerr = yerr_barrel_l1_[sizey - 1];
            else
              yerr = yerr_barrel_l1_def_;
          }
        } else {  // layer 2,3
          if (!edgex) {
            if (sizex <= xerr_barrel_ln_.size())
              xerr = xerr_barrel_ln_[sizex - 1];
            else
              xerr = xerr_barrel_ln_def_;
          }

          if (!edgey) {
            if (sizey <= yerr_barrel_ln_.size())
              yerr = yerr_barrel_ln_[sizey - 1];
            else
              yerr = yerr_barrel_ln_def_;
          }
        }

      } else {  // EndCap

        if (!edgex) {
          if (sizex <= xerr_endcap_.size())
            xerr = xerr_endcap_[sizex - 1];
          else
            xerr = xerr_endcap_def_;
        }

        if (!edgey) {
          if (sizey <= yerr_endcap_.size())
            yerr = yerr_endcap_[sizey - 1];
          else
            yerr = yerr_endcap_def_;
        }
      }  // end endcap
    }
  }
}

void PixelCPEGenericBase::fillPSetDescription(edm::ParameterSetDescription& desc) {
  desc.add<std::vector<double>>("xerr_barrel_l1", {0.00115, 0.00120, 0.00088});
  desc.add<std::vector<double>>("yerr_barrel_l1",
                                {0.00375, 0.00230, 0.00250, 0.00250, 0.00230, 0.00230, 0.00210, 0.00210, 0.00240});
  desc.add<std::vector<double>>("xerr_barrel_ln", {0.00115, 0.00120, 0.00088});
  desc.add<std::vector<double>>("yerr_barrel_ln",
                                {0.00375, 0.00230, 0.00250, 0.00250, 0.00230, 0.00230, 0.00210, 0.00210, 0.00240});
  desc.add<std::vector<double>>("xerr_endcap", {0.0020, 0.0020});
  desc.add<std::vector<double>>("yerr_endcap", {0.00210});
  desc.add<double>("xerr_barrel_l1_def", 0.01030);
  desc.add<double>("yerr_barrel_l1_def", 0.00210);
  desc.add<double>("xerr_barrel_ln_def", 0.01030);
  desc.add<double>("yerr_barrel_ln_def", 0.00210);
  desc.add<double>("xerr_endcap_def", 0.0020);
  desc.add<double>("yerr_endcap_def", 0.00075);
}