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"

namespace CLHEP {
  class HepRandomEngine;
}

class CSCDriftSim
{
public:
  CSCDriftSim();

  ~CSCDriftSim();

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

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

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(CLHEP::HepRandomEngine*);
  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;
};

#endif