d5/dd8/CSCDriftSim_8h_source.html

#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