CSCDriftSim.h

Go to the documentation of this file.

#ifndef MU_END_DRIFT_SIM_H
#define MU_END_DRIFT_SIM_H

#include <vector>
class CSCLayer;
class CSCDetectorHit;
class PSimHit;
class MagneticField;
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "CLHEP/Random/RandFlat.h"
#include "CLHEP/Random/RandGaussQ.h"

class CSCDriftSim
{
public:
  CSCDriftSim();

  ~CSCDriftSim();

  CSCDetectorHit getWireHit(const Local3DPoint & ionClusterPosition,
                  const CSCLayer *, int wire, const PSimHit & simHit);

  void setMagneticField(const MagneticField * field) {theMagneticField = field;}

  void setRandomEngine(CLHEP::HepRandomEngine& engine);

private:
  // helper functions
  double avgPathLengthLowB();
  double pathSigmaLowB();
  double avgDriftTimeLowB();
  double driftTimeSigmaLowB();
  double avgPathLengthHighB();
  double pathSigmaHighB();
  double avgDriftTimeHighB();
  double driftTimeSigmaHighB();
  double avgDrift() const;
  double driftSigma() const;
  double avalancheCharge();
  double gasGain(const CSCDetId & id) const;

  // local magnetic field
  float bz;
  // distances from wire
  double ycell, zcell;
  // dN/dE for the avalanche sim.
  std::vector<double> dNdEIntegral;
  const double STEP_SIZE;

  const double ELECTRON_DIFFUSION_COEFF;

  const MagneticField * theMagneticField;

  CLHEP::RandGaussQ * theRandGaussQ;
  CLHEP::RandFlat *   theRandFlat;
};

#endif