GEMSignalModel Class Reference

#include <GEMSignalModel.h>

Inheritance diagram for GEMSignalModel:

Public Member Functions
	GEMSignalModel (const edm::ParameterSet &)
int	getSimHitBx (const PSimHit *, CLHEP::HepRandomEngine *)
void	simulate (const GEMEtaPartition *, const edm::PSimHitContainer &, CLHEP::HepRandomEngine *, Strips &, DetectorHitMap &) override
std::vector< std::pair< int, int > >	simulateClustering (const TrapezoidalStripTopology *, const PSimHit *, const int, CLHEP::HepRandomEngine *)
	~GEMSignalModel () override
Public Member Functions inherited from GEMDigiModel
void	setGeometry (const GEMGeometry *geom)
virtual	~GEMDigiModel ()

Private Attributes
double	averageEfficiency_
double	averageShapingTime_
const double	cspeed
bool	digitizeOnlyMuons_
const double	elecEffLowCoeff
const double	elecEffLowParam0
const double	elecEffMidCoeff
const double	elecEffMidParam0
const double	elecEffMidParam1
const double	elecMomCut1
const double	elecMomCut2
const double	momConvFact
const int	muonPdgId
double	resolutionX_
double	signalPropagationSpeed_
double	timeJitter_
double	timeResolution_

Additional Inherited Members
Protected Member Functions inherited from GEMDigiModel
	GEMDigiModel (const edm::ParameterSet &)
Protected Attributes inherited from GEMDigiModel
const GEMGeometry *	geometry_

Detailed Description

Class for the GEM strip intrinsic noise simulation based on a very simple model Originally comes from GEMSimpleModel

Author

Sven Dildick by Roumyana Hadjiiska by Yechan Kang

Class for the GEM strip response simulation based on a very simple model Originally comes from GEMSimpleModel

Author

Sven Dildick by Roumyana Hadjiiska by Yechan Kang

Definition at line 24 of file GEMSignalModel.h.

Constructor & Destructor Documentation

GEMSignalModel::GEMSignalModel

(

const edm::ParameterSet &

config

)

Definition at line 16 of file GEMSignalModel.cc.

    : GEMDigiModel(config),
      averageEfficiency_(config.getParameter<double>("averageEfficiency")),
      averageShapingTime_(config.getParameter<double>("averageShapingTime")),
      timeResolution_(config.getParameter<double>("timeResolution")),
      timeJitter_(config.getParameter<double>("timeJitter")),
      signalPropagationSpeed_(config.getParameter<double>("signalPropagationSpeed")),
      digitizeOnlyMuons_(config.getParameter<bool>("digitizeOnlyMuons")),
      resolutionX_(config.getParameter<double>("resolutionX")),
      muonPdgId(13),
      cspeed(geant_units::operators::convertMmToCm(CLHEP::c_light)),
      momConvFact(1000.),
      elecMomCut1(1.96e-03),
      elecMomCut2(10.e-03),
      elecEffLowCoeff(1.7e-05),
      elecEffLowParam0(2.1),
      elecEffMidCoeff(1.34),
      elecEffMidParam0(-5.75e-01),
      elecEffMidParam1(7.96e-01) {}

GEMSignalModel::~GEMSignalModel ( )

override

Definition at line 36 of file GEMSignalModel.cc.

{}

Member Function Documentation

int GEMSignalModel::getSimHitBx	(	const PSimHit *	simhit,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 78 of file GEMSignalModel.cc.

References averageShapingTime_, l1GtPatternGenerator_cfi::bx, cspeed, PSimHit::detUnitId(), GEMGeometry::etaPartition(), Exception, GEMDigiModel::geometry_, triggerObjects_cff::id, PSimHit::localPosition(), LogDebug, signalPropagationSpeed_, mathSSE::sqrt(), timeJitter_, PSimHit::timeOfFlight(), timeResolution_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by simulate().

                                                                                   {
  int bx = -999;
  const LocalPoint simHitPos(simhit->localPosition());
  const float tof(simhit->timeOfFlight());
  // random Gaussian time correction due to electronics jitter
  float randomJitterTime = CLHEP::RandGaussQ::shoot(engine, 0., timeJitter_);
  const GEMDetId id(simhit->detUnitId());
  const GEMEtaPartition* roll(geometry_->etaPartition(id));
  if (!roll) {
    throw cms::Exception("Geometry") << "GEMSignalModel::getSimHitBx() - GEM simhit id does not match any GEM roll id: "
                                     << id << "\n";
    return 999;
  }
  if (roll->id().region() == 0) {
    throw cms::Exception("Geometry")
        << "GEMSignalModel::getSimHitBx() - this GEM id is from barrel, which cannot happen: " << roll->id() << "\n";
  }
  const int nstrips = roll->nstrips();
  float middleStrip = nstrips / 2.;
  const LocalPoint& middleOfRoll = roll->centreOfStrip(middleStrip);
  const GlobalPoint& globMiddleRol = roll->toGlobal(middleOfRoll);
  double muRadius = sqrt(globMiddleRol.x() * globMiddleRol.x() + globMiddleRol.y() * globMiddleRol.y() +
                         globMiddleRol.z() * globMiddleRol.z());
  double timeCalibrationOffset_ = muRadius / cspeed;  //[ns]

  const TrapezoidalStripTopology* top(dynamic_cast<const TrapezoidalStripTopology*>(&(roll->topology())));
  const float halfStripLength(0.5 * top->stripLength());
  const float distanceFromEdge(halfStripLength - simHitPos.y());

  // signal propagation speed in material in [cm/ns]
  double signalPropagationSpeedTrue = signalPropagationSpeed_ * cspeed;

  // average time for the signal to propagate from the SimHit to the top of a strip
  const float averagePropagationTime(distanceFromEdge / signalPropagationSpeedTrue);
  // random Gaussian time correction due to the finite timing resolution of the detector
  float randomResolutionTime = CLHEP::RandGaussQ::shoot(engine, 0., timeResolution_);
  const float simhitTime(tof + averageShapingTime_ + randomResolutionTime + averagePropagationTime + randomJitterTime);
  float referenceTime = 0.;
  referenceTime = timeCalibrationOffset_ + halfStripLength / signalPropagationSpeedTrue + averageShapingTime_;
  const float timeDifference(simhitTime - referenceTime);
  // assign the bunch crossing
  bx = static_cast<int>(std::round((timeDifference) / 25.));

  // check time
  LogDebug("GEMDigiProducer") << "checktime "
                              << "bx = " << bx << "\tdeltaT = " << timeDifference << "\tsimT =  " << simhitTime
                              << "\trefT =  " << referenceTime << "\ttof = " << tof
                              << "\tavePropT =  " << averagePropagationTime
                              << "\taveRefPropT = " << halfStripLength / signalPropagationSpeedTrue << "\n";
  return bx;
}

void GEMSignalModel::simulate ( const GEMEtaPartition *  roll, const edm::PSimHitContainer &  simHits, CLHEP::HepRandomEngine *  engine, Strips &  strips_, DetectorHitMap &  detectorHitMap_  )

overridevirtual

Implements GEMDigiModel.

Definition at line 38 of file GEMSignalModel.cc.

References funct::abs(), averageEfficiency_, l1GtPatternGenerator_cfi::bx, digitizeOnlyMuons_, elecEffLowCoeff, elecEffLowParam0, elecEffMidCoeff, elecEffMidParam0, elecEffMidParam1, elecMomCut1, elecMomCut2, JetChargeProducer_cfi::exp, getSimHitBx(), dqm-mbProfile::log, momConvFact, muonPdgId, simulateClustering(), and GEMEtaPartition::topology().

                                                               {
  bool digiMuon = false;
  bool digiElec = false;
  const TrapezoidalStripTopology* top(dynamic_cast<const TrapezoidalStripTopology*>(&(roll->topology())));
  for (const auto& hit : simHits) {
    if (std::abs(hit.particleType()) != muonPdgId && digitizeOnlyMuons_)
      continue;
    double elecEff = 0.;
    double partMom = hit.pabs();
    double checkMuonEff = CLHEP::RandFlat::shoot(engine, 0., 1.);
    double checkElecEff = CLHEP::RandFlat::shoot(engine, 0., 1.);
    if (std::abs(hit.particleType()) == muonPdgId && checkMuonEff < averageEfficiency_)
      digiMuon = true;
    if (std::abs(hit.particleType()) != muonPdgId)  //consider all non muon particles with gem efficiency to electrons
    {
      if (partMom <= elecMomCut1)
        elecEff = elecEffLowCoeff * std::exp(elecEffLowParam0 * partMom * momConvFact);
      if (partMom > elecMomCut1 && partMom < elecMomCut2)
        elecEff = elecEffMidCoeff * log(elecEffMidParam1 * partMom * momConvFact - elecEffMidParam0) /
                  (elecEffMidCoeff + log(elecEffMidParam1 * partMom * momConvFact - elecEffMidParam0));
      if (partMom > elecMomCut2)
        elecEff = 1.;
      if (checkElecEff < elecEff)
        digiElec = true;
    }
    if (!(digiMuon || digiElec))
      continue;
    const int bx(getSimHitBx(&hit, engine));
    const std::vector<std::pair<int, int> >& cluster(simulateClustering(top, &hit, bx, engine));
    for (const auto& digi : cluster) {
      detectorHitMap_.emplace(digi, &hit);
      strips_.emplace(digi);
    }
  }
}

std::vector< std::pair< int, int > > GEMSignalModel::simulateClustering	(	const TrapezoidalStripTopology *	top,
		const PSimHit *	simHit,
		const int	bx,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 130 of file GEMSignalModel.cc.

References funct::abs(), TrapezoidalStripTopology::channel(), PSimHit::entryPoint(), PSimHit::exitPoint(), resolutionX_, digitizers_cfi::strip, and PV3DBase< T, PVType, FrameType >::x().

Referenced by simulate().

                                                                                                   {
  const LocalPoint& hit_entry(simHit->entryPoint());
  const LocalPoint& hit_exit(simHit->exitPoint());

  LocalPoint start_point, end_point;
  if (hit_entry.x() < hit_exit.x()) {
    start_point = hit_entry;
    end_point = hit_exit;
  } else {
    start_point = hit_exit;
    end_point = hit_entry;
  }

  // Add Gaussian noise to the points towards outside.
  float smeared_start_x = start_point.x() - std::abs(CLHEP::RandGaussQ::shoot(engine, 0, resolutionX_));
  float smeared_end_x = end_point.x() + std::abs(CLHEP::RandGaussQ::shoot(engine, 0, resolutionX_));

  LocalPoint smeared_start_point(smeared_start_x, start_point.y(), start_point.z());
  LocalPoint smeared_end_point(smeared_end_x, end_point.y(), end_point.z());

  int cluster_start = top->channel(smeared_start_point);
  int cluster_end = top->channel(smeared_end_point);

  std::vector<std::pair<int, int> > cluster;
  for (int strip = cluster_start; strip <= cluster_end; strip++) {
    cluster.emplace_back(strip, bx);
  }

  return cluster;
}

Member Data Documentation

double GEMSignalModel::averageEfficiency_

private

Definition at line 41 of file GEMSignalModel.h.

Referenced by simulate().

double GEMSignalModel::averageShapingTime_

private

Definition at line 42 of file GEMSignalModel.h.

Referenced by getSimHitBx().

const double GEMSignalModel::cspeed

private

Definition at line 50 of file GEMSignalModel.h.

Referenced by getSimHitBx().

bool GEMSignalModel::digitizeOnlyMuons_

private

Definition at line 46 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecEffLowCoeff

private

Definition at line 54 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecEffLowParam0

private

Definition at line 55 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecEffMidCoeff

private

Definition at line 56 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecEffMidParam0

private

Definition at line 57 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecEffMidParam1

private

Definition at line 58 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecMomCut1

private

Definition at line 52 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::elecMomCut2

private

Definition at line 53 of file GEMSignalModel.h.

Referenced by simulate().

const double GEMSignalModel::momConvFact

private

Definition at line 51 of file GEMSignalModel.h.

Referenced by simulate().

const int GEMSignalModel::muonPdgId

private

Definition at line 49 of file GEMSignalModel.h.

Referenced by simulate().

double GEMSignalModel::resolutionX_

private

Definition at line 47 of file GEMSignalModel.h.

Referenced by simulateClustering().

double GEMSignalModel::signalPropagationSpeed_

private

Definition at line 45 of file GEMSignalModel.h.

Referenced by getSimHitBx().

double GEMSignalModel::timeJitter_

private

Definition at line 44 of file GEMSignalModel.h.

Referenced by getSimHitBx().

double GEMSignalModel::timeResolution_

private

Definition at line 43 of file GEMSignalModel.h.

Referenced by getSimHitBx().