GEMBkgModel Class Reference

#include <GEMBkgModel.h>

Inheritance diagram for GEMBkgModel:

Public Member Functions
	GEMBkgModel (const edm::ParameterSet &)
void	simulate (const GEMEtaPartition *, const edm::PSimHitContainer &, CLHEP::HepRandomEngine *, Strips &, DetectorHitMap &) override
	~GEMBkgModel () override
Public Member Functions inherited from GEMDigiModel
void	setGeometry (const GEMGeometry *geom)
virtual	~GEMDigiModel ()

Private Attributes
double	averageEfficiency_
double	bxWidth_
const double	clusterSizeCut
bool	digitizeOnlyMuons_
bool	fixedRollRadius_
double	GE11ElecBkgParam0_
double	GE11ElecBkgParam1_
double	GE11ElecBkgParam2_
double	GE11ModNeuBkgParam0_
double	GE11ModNeuBkgParam1_
double	GE11ModNeuBkgParam2_
double	GE21ElecBkgParam0_
double	GE21ElecBkgParam1_
double	GE21ElecBkgParam2_
double	GE21ModNeuBkgParam0_
double	GE21ModNeuBkgParam1_
double	GE21ModNeuBkgParam2_
double	instLumi_
int	maxBunch_
int	minBunch_
double	rateFact_
const double	referenceInstLumi_
double	resolutionX_
bool	simulateElectronBkg_
bool	simulateNoiseCLS_

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 response to background simulation based on a very simple model Originally comes from GEMSimpleModel

Author

Sven Dildick by Roumyana Hadjiiska by Yechan Kang

Definition at line 23 of file GEMBkgModel.h.

Constructor & Destructor Documentation

GEMBkgModel::GEMBkgModel

(

const edm::ParameterSet &

config

)

Definition at line 13 of file GEMBkgModel.cc.

    : GEMDigiModel(config),
      clusterSizeCut(0.53),
      averageEfficiency_(config.getParameter<double>("averageEfficiency")),
      minBunch_(config.getParameter<int>("minBunch")),
      maxBunch_(config.getParameter<int>("maxBunch")),
      simulateNoiseCLS_(config.getParameter<bool>("simulateNoiseCLS")),
      fixedRollRadius_(config.getParameter<bool>("fixedRollRadius")),
      simulateElectronBkg_(config.getParameter<bool>("simulateElectronBkg")),
      instLumi_(config.getParameter<double>("instLumi")),
      rateFact_(config.getParameter<double>("rateFact")),
      bxWidth_(config.getParameter<double>("bxWidth")),
      referenceInstLumi_(config.getParameter<double>("referenceInstLumi")),
      GE11ElecBkgParam0_(config.getParameter<double>("GE11ElecBkgParam0")),
      GE11ElecBkgParam1_(config.getParameter<double>("GE11ElecBkgParam1")),
      GE11ElecBkgParam2_(config.getParameter<double>("GE11ElecBkgParam2")),
      GE21ElecBkgParam0_(config.getParameter<double>("GE21ElecBkgParam0")),
      GE21ElecBkgParam1_(config.getParameter<double>("GE21ElecBkgParam1")),
      GE21ElecBkgParam2_(config.getParameter<double>("GE21ElecBkgParam2")),
      GE11ModNeuBkgParam0_(config.getParameter<double>("GE11ModNeuBkgParam0")),
      GE11ModNeuBkgParam1_(config.getParameter<double>("GE11ModNeuBkgParam1")),
      GE11ModNeuBkgParam2_(config.getParameter<double>("GE11ModNeuBkgParam2")),
      GE21ModNeuBkgParam0_(config.getParameter<double>("GE11ModNeuBkgParam0")),
      GE21ModNeuBkgParam1_(config.getParameter<double>("GE11ModNeuBkgParam1")),
      GE21ModNeuBkgParam2_(config.getParameter<double>("GE11ModNeuBkgParam2")) {}

GEMBkgModel::~GEMBkgModel ( )

override

Definition at line 39 of file GEMBkgModel.cc.

{}

Member Function Documentation

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

overridevirtual

Implements GEMDigiModel.

Definition at line 41 of file GEMBkgModel.cc.

References averageEfficiency_, bxWidth_, clusterSizeCut, Exception, fixedRollRadius_, GE11ElecBkgParam0_, GE11ElecBkgParam1_, GE11ElecBkgParam2_, GE11ModNeuBkgParam0_, GE11ModNeuBkgParam1_, GE11ModNeuBkgParam2_, GE21ElecBkgParam0_, GE21ElecBkgParam1_, GE21ElecBkgParam2_, GE21ModNeuBkgParam0_, GE21ModNeuBkgParam1_, GE21ModNeuBkgParam2_, mps_fire::i, GEMEtaPartition::id(), instLumi_, maxBunch_, minBunch_, GEMEtaPartition::nstrips(), GEMEtaPartition::pitch(), rateFact_, referenceInstLumi_, simulateElectronBkg_, simulateNoiseCLS_, and GEMEtaPartition::topology().

                                                            {
  const GEMDetId& gemId(roll->id());
  const int nstrips(roll->nstrips());
  double trArea(0.0);
  double trStripArea(0.0);
  if (gemId.region() == 0) {
    throw cms::Exception("Geometry") << "GEMBkgModel::simulate() - this GEM id is from barrel, which cannot happen.";
  }
  const GEMStripTopology* top_(dynamic_cast<const GEMStripTopology*>(&(roll->topology())));
  const float striplength(top_->stripLength());
  trStripArea = (roll->pitch()) * striplength;
  trArea = trStripArea * nstrips;
  const int nBxing(maxBunch_ - minBunch_ + 1);
  const float rollRadius(
      fixedRollRadius_
          ? top_->radius()
          : top_->radius() + CLHEP::RandFlat::shoot(engine, -1. * top_->stripLength() / 2., top_->stripLength() / 2.));

  //calculate noise from model
  double averageNeutralNoiseRatePerRoll = 0.;
  double averageNoiseElectronRatePerRoll = 0.;
  double averageNoiseRatePerRoll = 0.;
  if (gemId.station() == 1) {
    //simulate neutral background for GE1/1
    averageNeutralNoiseRatePerRoll =
        (GE11ModNeuBkgParam0_ + GE11ModNeuBkgParam1_ * rollRadius +
         GE11ModNeuBkgParam2_ * rollRadius * rollRadius);  //simulate electron background for GE1/1
    if (simulateElectronBkg_)
      averageNoiseElectronRatePerRoll =
          (GE11ElecBkgParam0_ + GE11ElecBkgParam1_ * rollRadius + GE11ElecBkgParam2_ * rollRadius * rollRadius);

    // Scale up/down for desired instantaneous lumi (reference is 5E34, double from config is in units of 1E34)
    averageNoiseRatePerRoll = averageNeutralNoiseRatePerRoll + averageNoiseElectronRatePerRoll;
    averageNoiseRatePerRoll *= instLumi_ * rateFact_ / referenceInstLumi_;
  } else if (gemId.station() == 2) {
    //simulate neutral background for GE2/1
    averageNeutralNoiseRatePerRoll =
        (GE21ModNeuBkgParam0_ + GE21ModNeuBkgParam1_ * rollRadius + GE21ModNeuBkgParam2_ * rollRadius * rollRadius);
    //simulate electron background for GE2/1
    if (simulateElectronBkg_)
      averageNoiseElectronRatePerRoll =
          (GE21ElecBkgParam0_ + GE21ElecBkgParam1_ * rollRadius + GE21ElecBkgParam2_ * rollRadius * rollRadius);

    // Scale up/down for desired instantaneous lumi (reference is 5E34, double from config is in units of 1E34)
    averageNoiseRatePerRoll = averageNeutralNoiseRatePerRoll + averageNoiseElectronRatePerRoll;
    averageNoiseRatePerRoll *= instLumi_ * rateFact_ / referenceInstLumi_;
  }

  //simulate bkg contribution
  const double averageNoise(averageNoiseRatePerRoll * nBxing * trArea * bxWidth_);
  CLHEP::RandPoissonQ randPoissonQ(*engine, averageNoise);
  const int n_hits(randPoissonQ.fire());
  for (int i = 0; i < n_hits; ++i) {
    const int centralStrip(static_cast<int>(CLHEP::RandFlat::shoot(engine, 1, nstrips)));
    const int time_hit(static_cast<int>(CLHEP::RandFlat::shoot(engine, nBxing)) + minBunch_);
    if (simulateNoiseCLS_) {
      std::vector<std::pair<int, int> > cluster_;
      cluster_.clear();
      cluster_.emplace_back(centralStrip, time_hit);
      int clusterSize((CLHEP::RandFlat::shoot(engine)) <= clusterSizeCut ? 1 : 2);
      if (clusterSize == 2) {
        if (CLHEP::RandFlat::shoot(engine) < 0.5) {
          if (CLHEP::RandFlat::shoot(engine) < averageEfficiency_ && (centralStrip - 1 > 0))
            cluster_.emplace_back(centralStrip - 1, time_hit);
        } else {
          if (CLHEP::RandFlat::shoot(engine) < averageEfficiency_ && (centralStrip + 1 <= nstrips))
            cluster_.emplace_back(centralStrip + 1, time_hit);
        }
      }
      for (const auto& digi : cluster_) {
        strips_.emplace(digi);
      }
    }  //end simulateNoiseCLS_
    else {
      strips_.emplace(centralStrip, time_hit);
    }
  }
  return;
}

Member Data Documentation

double GEMBkgModel::averageEfficiency_

private

Definition at line 34 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::bxWidth_

private

Definition at line 43 of file GEMBkgModel.h.

Referenced by simulate().

const double GEMBkgModel::clusterSizeCut

private

Definition at line 33 of file GEMBkgModel.h.

Referenced by simulate().

bool GEMBkgModel::digitizeOnlyMuons_

private

Definition at line 37 of file GEMBkgModel.h.

bool GEMBkgModel::fixedRollRadius_

private

Definition at line 39 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE11ElecBkgParam0_

private

Definition at line 48 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE11ElecBkgParam1_

private

Definition at line 49 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE11ElecBkgParam2_

private

Definition at line 50 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE11ModNeuBkgParam0_

private

Definition at line 55 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE11ModNeuBkgParam1_

private

Definition at line 56 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE11ModNeuBkgParam2_

private

Definition at line 57 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE21ElecBkgParam0_

private

Definition at line 51 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE21ElecBkgParam1_

private

Definition at line 52 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE21ElecBkgParam2_

private

Definition at line 53 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE21ModNeuBkgParam0_

private

Definition at line 58 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE21ModNeuBkgParam1_

private

Definition at line 59 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::GE21ModNeuBkgParam2_

private

Definition at line 60 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::instLumi_

private

Definition at line 41 of file GEMBkgModel.h.

Referenced by simulate().

int GEMBkgModel::maxBunch_

private

Definition at line 36 of file GEMBkgModel.h.

Referenced by simulate().

int GEMBkgModel::minBunch_

private

Definition at line 35 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::rateFact_

private

Definition at line 42 of file GEMBkgModel.h.

Referenced by simulate().

const double GEMBkgModel::referenceInstLumi_

private

Definition at line 44 of file GEMBkgModel.h.

Referenced by simulate().

double GEMBkgModel::resolutionX_

private

Definition at line 45 of file GEMBkgModel.h.

bool GEMBkgModel::simulateElectronBkg_

private

Definition at line 40 of file GEMBkgModel.h.

Referenced by simulate().

bool GEMBkgModel::simulateNoiseCLS_

private

Definition at line 38 of file GEMBkgModel.h.

Referenced by simulate().