SiStripGaussianSmearingRecHitConverterAlgorithm.cc

Go to the documentation of this file.

// SiStrip Gaussian Smearing
#include "FastSimulation/TrackingRecHitProducer/interface/SiStripGaussianSmearingRecHitConverterAlgorithm.h"

// Famos
#include "FastSimulation/Utilities/interface/RandomEngineAndDistribution.h"

// STL

// #define FAMOS_DEBUG

SiStripGaussianSmearingRecHitConverterAlgorithm::SiStripGaussianSmearingRecHitConverterAlgorithm()
{ }

void 
SiStripGaussianSmearingRecHitConverterAlgorithm::smearHit(const PSimHit& simHit, 
                              double localPositionResolutionX, 
                              double localPositionResolutionY,  
                              double localPositionResolutionZ,
                                                          double boundX,
                                                          double boundY,
                                                          RandomEngineAndDistribution const* random)
{

  // Gaussian Smearing
  // x is smeared,
  // y fixed at the centre of the strip, 
  // z fixed at the centre of the active area
  // For the double-sided modules it will be a problem 
  // for the RecHit matcher (starting from these PSimHits)
  //
  do {
    thePosition = 
      Local3DPoint(random->gaussShoot((double)simHit.localPosition().x(), localPositionResolutionX), 
           0.,0.);
#ifdef FAMOS_DEBUG
    std::cout << " Detector bounds: "
              << "\t\tx = " << boundX
              << "\ty = " << boundY
              << std::endl;
    std::cout << " Generated local position "
              << "\tx = " << thePosition.x()
              << "\ty = " << thePosition.y()
              << std::endl;       
#endif  
  } while(fabs(thePosition.x()) > boundX);
  
  //
  thePositionX = thePosition.x();
  thePositionY = thePosition.y();
  thePositionZ = thePosition.z();
  //
  theErrorX = localPositionResolutionX;
  theErrorY = localPositionResolutionY;
  theErrorZ = localPositionResolutionZ;

  theError = LocalError( theErrorX * theErrorX,
             0.0,
             theErrorY * theErrorY ); 
  // Local Error is 2D: (xx,xy,yy), square of sigma 
  // in first an third position as for resolution matrix
  //
}