TrackingRecHitNoSmearingPlugin.cc

Go to the documentation of this file.

#include "FastSimulation/TrackingRecHitProducer/interface/TrackingRecHitAlgorithm.h"
#include "FastSimulation/TrackingRecHitProducer/interface/TrackingRecHitAlgorithmFactory.h"
#include "FastSimulation/TrackingRecHitProducer/interface/TrackingRecHitProduct.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"
#include "DataFormats/TrackerRecHit2D/interface/FastSingleTrackerRecHit.h"

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

#include <string>
#include <iostream>

class TrackingRecHitNoSmearingPlugin : public TrackingRecHitAlgorithm {
private:
  double _errorXX;
  double _errorXY;
  double _errorYY;

public:
  TrackingRecHitNoSmearingPlugin(const std::string& name,
                                 const edm::ParameterSet& config,
                                 edm::ConsumesCollector& consumesCollector)
      : TrackingRecHitAlgorithm(name, config, consumesCollector),
        _errorXX(0.001 * 0.001),
        _errorXY(0.0),
        _errorYY(0.001 * 0.001) {
    if (config.exists("errorXX")) {
      _errorXX = config.getParameter<double>("errorXX");
    }

    if (config.exists("errorXY")) {
      _errorXY = config.getParameter<double>("errorXY");
    }

    if (config.exists("errorYY")) {
      _errorYY = config.getParameter<double>("errorYY");
    }
  }

  TrackingRecHitProductPtr process(TrackingRecHitProductPtr product) const override {
    for (const std::pair<unsigned int, const PSimHit*>& simHitIdPair : product->getSimHitIdPairs()) {
      const PSimHit* simHit = simHitIdPair.second;
      const Local3DPoint& position = simHit->localPosition();
      LocalError error(_errorXX, _errorXY, _errorYY);
      const GeomDet* geomDet = getTrackerGeometry().idToDetUnit(product->getDetId());

      //TODO: this is only a minimal example
      FastSingleTrackerRecHit recHit(
          position,  //const LocalPoint &
          error,     //const LocalError &
          *geomDet,  //GeomDet const &idet
          fastTrackerRecHitType::
              siPixel  // since this is a dummy class anyway: pretend all hits are pixel hits (only effect: hits are defined in 2D (?))
      );
      product->addRecHit(recHit, {simHitIdPair});
    }

    return product;
  }
};

DEFINE_EDM_PLUGIN(TrackingRecHitAlgorithmFactory, TrackingRecHitNoSmearingPlugin, "TrackingRecHitNoSmearingPlugin");