---

SiPixelGaussianSmearingRecHitConverterAlgorithm.h

Go to the documentation of this file.

#ifndef FastSimulation_TrackingRecHitProducer_SiPixelGaussianSmearingRecHitConverterAlgorithm_h
#define FastSimulation_TrackingRecHitProducer_SiPixelGaussianSmearingRecHitConverterAlgorithm_h

//---------------------------------------------------------------------------
//---------------------------------------------------------------------------


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

// PSimHit
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"

// Geometry
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"

// Vectors
#include "DataFormats/GeometryVector/interface/Point3DBase.h"
#include "DataFormats/GeometrySurface/interface/LocalError.h"

// STL
#include <vector>
#include <string>

class TFile;
class TH1F;
class PixelErrorParametrization;
class RandomEngine;
class HistogramGenerator;

class SiPixelGaussianSmearingRecHitConverterAlgorithm {
public:
  //--- Constructor, virtual destructor (just in case)
  explicit SiPixelGaussianSmearingRecHitConverterAlgorithm(                
   const edm::ParameterSet& pset,
   GeomDetType::SubDetector pixelPart,
   std::vector<TH1F*>& alphaMultiplicityCumulativeProbabilities,
   std::vector<TH1F*>& betaMultiplicityCumulativeProbabilities,
   TFile* pixelResolutionFile,
   const RandomEngine* engine);

  // destructor
  virtual ~SiPixelGaussianSmearingRecHitConverterAlgorithm();
  
  // return results
  Local3DPoint getPosition()  {return thePosition;}
  double       getPositionX() {return thePositionX;}
  double       getPositionY() {return thePositionY;}
  double       getPositionZ() {return thePositionZ;}
  LocalError   getError()     {return theError;}
  double       getErrorX()    {return theErrorX;}
  double       getErrorY()    {return theErrorY;}
  double       getErrorZ()    {return theErrorZ;}
  unsigned int getPixelMultiplicityAlpha() {return thePixelMultiplicityAlpha;}
  unsigned int getPixelMultiplicityBeta()  {return thePixelMultiplicityBeta;}
  //
  
  //
  void smearHit( const PSimHit& simHit, const PixelGeomDetUnit* detUnit, const double boundX, const double boundY);

private:
  // Switch between old (ORCA) and new (CMSSW) pixel parameterization
  bool useCMSSWPixelParameterization;
  //
  bool isFlipped(const PixelGeomDetUnit* theDet) const;
  //
  bool negativeErrorProtection; // in case it is true protect against PixelErrorParametrization negative variance
  //
  // resolution bins
  double resAlpha_binMin , resAlpha_binWidth;
  unsigned int resAlpha_binN;
  double resBeta_binMin  , resBeta_binWidth;
  unsigned int resBeta_binN;
  //
  edm::ParameterSet pset_;
  // Useful private members
  GeomDetType::SubDetector thePixelPart;

  std::vector<TH1F*> theAlphaMultiplicityCumulativeProbabilities;
  std::vector<TH1F*> theBetaMultiplicityCumulativeProbabilities;

  std::map<unsigned,const HistogramGenerator*> theAlphaHistos;
  std::map<unsigned,const HistogramGenerator*> theBetaHistos;

  TFile* thePixelResolutionFile;

  unsigned int theLayer;
  // output
  Local3DPoint thePosition;
  double       thePositionX;
  double       thePositionY;
  double       thePositionZ;
  LocalError   theError;
  double       theErrorX;
  double       theErrorY;
  double       theErrorZ;
  unsigned int thePixelMultiplicityAlpha;
  unsigned int thePixelMultiplicityBeta;
  //
  PixelErrorParametrization* pixelError;

  // The random engine
  const RandomEngine* random;
  
};

#endif

---