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