SiLinearChargeDivider.h

Go to the documentation of this file.

#ifndef Tracker_SiLinearChargeDivider_H
#define Tracker_SiLinearChargeDivider_H

#include <memory>

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

#include "SiChargeDivider.h"
#include "SimTracker/Common/interface/SiG4UniversalFluctuation.h"
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"

namespace CLHEP{
  class HepRandomEngine;
}

class SiLinearChargeDivider : public SiChargeDivider{
 public:

  // constructor
  SiLinearChargeDivider(const edm::ParameterSet& conf);

  // destructor
  virtual ~SiLinearChargeDivider();

  // main method: divide the charge (from the PSimHit) into several energy deposits in the bulk
  SiChargeDivider::ionization_type divide(const PSimHit*, const LocalVector&, double, const StripGeomDetUnit& det, CLHEP::HepRandomEngine*);

  // set the ParticleDataTable (used to fluctuate the charge properly)
  void setParticleDataTable(const ParticleDataTable * pdt) { theParticleDataTable = pdt; }
  
 private:
  // configuration data
  const bool   peakMode;
  const bool   fluctuateCharge;
  const int    chargedivisionsPerStrip;
  const double deltaCut ;
  const double cosmicShift;

  const ParticleDataTable * theParticleDataTable;

  // Geant4 engine used by fluctuateEloss()
  std::unique_ptr<SiG4UniversalFluctuation> fluctuate; 
  // utility: drifts the charge to the surface to estimate the number of relevant strips
  inline float driftXPos(const Local3DPoint& pos, const LocalVector& drift, double thickness) { 
    return pos.x()+(thickness/2.-pos.z())*drift.x()/drift.z();
  }
  // fluctuate the Eloss
  void fluctuateEloss(double const particleMass, float momentum, float eloss, float length, int NumberOfSegmentation, float elossVector[], CLHEP::HepRandomEngine*);
  // time response (from the pulse shape)
  inline float TimeResponse( const PSimHit* hit, const StripGeomDetUnit& det) {
    return (peakMode ? PeakShape(hit,det) : DeconvolutionShape(hit,det));
  } 
  // pulse shape in peak mode
  float PeakShape(const PSimHit*, const StripGeomDetUnit& det);
  // pulse shape in deconvolution mode
  float DeconvolutionShape( const PSimHit*, const StripGeomDetUnit& det);
  // data table for pulse shape in peak mode
  static float const peakValues[921];
  // data table for pulse shape in deconvolution mode
  static float const decoValues[651];

};

#endif