#include <ME0PreRecoGaussianModel.h>

Inheritance diagram for ME0PreRecoGaussianModel:

Public Member Functions
	ME0PreRecoGaussianModel (const edm::ParameterSet &)

void	setup ()

void	simulateNoise (const ME0EtaPartition , CLHEP::HepRandomEngine ) override

void	simulateSignal (const ME0EtaPartition , const edm::PSimHitContainer &, CLHEP::HepRandomEngine ) override

	~ME0PreRecoGaussianModel ()

Public Member Functions inherited from ME0DigiPreRecoModel
void	fillDigis (int rollDetId, ME0DigiPreRecoCollection &)

const ME0Geometry *	getGeometry ()

void	setGeometry (const ME0Geometry *geom)

virtual	~ME0DigiPreRecoModel ()

Private Attributes
double	averageEfficiency_

double	averageNoiseRate_

int	bxwidth_

bool	corr

bool	digitizeOnlyMuons_

bool	doBkgNoise_

bool	etaproj

int	maxBunch_

double	ME0ModElecBkgParam0

double	ME0ModElecBkgParam1

double	ME0ModElecBkgParam2

double	ME0ModElecBkgParam3

double	ME0ModElecBkgParam4

double	ME0ModElecBkgParam5

double	ME0ModElecBkgParam6

double	ME0ModElecBkgParam7

double	ME0ModNeuBkgParam0

double	ME0ModNeuBkgParam1

double	ME0ModNeuBkgParam2

double	ME0ModNeuBkgParam3

double	ME0ModNeuBkgParam4

double	ME0ModNeuBkgParam5

double	ME0ModNeuBkgParam6

int	minBunch_

double	sigma_t

double	sigma_u

double	sigma_v

bool	simulateElectronBkg_

bool	simulateIntrinsicNoise_

Additional Inherited Members
Protected Member Functions inherited from ME0DigiPreRecoModel
	ME0DigiPreRecoModel (const edm::ParameterSet &)

Protected Attributes inherited from ME0DigiPreRecoModel
std::set< ME0DigiPreReco >	digi_

const ME0Geometry *	geometry_

Detailed Description

Class for the ME0 Gaussian response simulation as pre-reco step

Definition at line 18 of file ME0PreRecoGaussianModel.h.

Constructor & Destructor Documentation

ME0PreRecoGaussianModel::ME0PreRecoGaussianModel ( const edm::ParameterSet & config )

Definition at line 14 of file ME0PreRecoGaussianModel.cc.

References ME0ModElecBkgParam0, ME0ModElecBkgParam1, ME0ModElecBkgParam2, ME0ModElecBkgParam3, ME0ModElecBkgParam4, ME0ModElecBkgParam5, ME0ModElecBkgParam6, ME0ModNeuBkgParam0, ME0ModNeuBkgParam1, ME0ModNeuBkgParam2, ME0ModNeuBkgParam3, ME0ModNeuBkgParam4, ME0ModNeuBkgParam5, and ME0ModNeuBkgParam6.

                                                                               :
   ME0DigiPreRecoModel(config)
 , sigma_t(config.getParameter<double> ("timeResolution"))
 , sigma_u(config.getParameter<double> ("phiResolution"))
 , sigma_v(config.getParameter<double> ("etaResolution"))
 , corr(config.getParameter<bool> ("useCorrelation"))
 , etaproj(config.getParameter<bool> ("useEtaProjectiveGEO"))
 , digitizeOnlyMuons_(config.getParameter<bool> ("digitizeOnlyMuons"))
 , averageEfficiency_(config.getParameter<double> ("averageEfficiency"))
 , doBkgNoise_(config.getParameter<bool> ("doBkgNoise"))
 , simulateIntrinsicNoise_(config.getParameter<bool> ("simulateIntrinsicNoise"))
 , simulateElectronBkg_(config.getParameter<bool>("simulateElectronBkg"))
 , averageNoiseRate_(config.getParameter<double> ("averageNoiseRate"))
 , bxwidth_(config.getParameter<int> ("bxwidth"))
 , minBunch_(config.getParameter<int> ("minBunch"))
 , maxBunch_(config.getParameter<int> ("maxBunch"))
 
 {
   //params for the simple pol6 model of neutral and electron bkg for ME0:
   ME0ModNeuBkgParam0 = 5.69e+06;
   ME0ModNeuBkgParam1 = -293334;
   ME0ModNeuBkgParam2 = 6279.6;
   ME0ModNeuBkgParam3 = -71.2928;
   ME0ModNeuBkgParam4 = 0.452244;
   ME0ModNeuBkgParam5 = -0.0015191;
   ME0ModNeuBkgParam6 = 2.1106e-06;
 
   ME0ModElecBkgParam0 = 3.77712e+06;
   ME0ModElecBkgParam1 = -199280;
   ME0ModElecBkgParam2 = 4340.69;
   ME0ModElecBkgParam3 = -49.922;
   ME0ModElecBkgParam4 = 0.319699;
   ME0ModElecBkgParam5 = -0.00108113;
   ME0ModElecBkgParam6 = 1.50889e-06;
 
 }

ME0PreRecoGaussianModel::~ME0PreRecoGaussianModel ( )

Definition at line 51 of file ME0PreRecoGaussianModel.cc.

52 {

53 }

Member Function Documentation

void ME0PreRecoGaussianModel::setup ( void )

inlinevirtual

Implements ME0DigiPreRecoModel.

Definition at line 30 of file ME0PreRecoGaussianModel.h.

30 {}

void ME0PreRecoGaussianModel::simulateNoise	(	const ME0EtaPartition *	roll,
		CLHEP::HepRandomEngine *	engine
	)

overridevirtual

Implements ME0DigiPreRecoModel.

Definition at line 77 of file ME0PreRecoGaussianModel.cc.

References bxwidth_, corr, ME0DigiPreRecoModel::digi_, doBkgNoise_, alignCSCRings::e, Exception, i, ME0EtaPartition::id(), maxBunch_, ME0ModElecBkgParam0, ME0ModElecBkgParam1, ME0ModElecBkgParam2, ME0ModElecBkgParam3, ME0ModElecBkgParam4, ME0ModElecBkgParam5, ME0ModElecBkgParam6, ME0ModNeuBkgParam0, ME0ModNeuBkgParam1, ME0ModNeuBkgParam2, ME0ModNeuBkgParam3, ME0ModNeuBkgParam4, ME0ModNeuBkgParam5, ME0ModNeuBkgParam6, minBunch_, Parameters::parameters, ME0EtaPartitionSpecs::parameters(), sigma_t, sigma_u, sigma_v, simulateElectronBkg_, ME0EtaPartition::specs(), mathSSE::sqrt(), GeomDet::toGlobal(), ME0EtaPartition::topology(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

 {
   const double cspeed = 299792458;
   double trArea(0.0);
   const ME0DetId me0Id(roll->id());
 
   if (me0Id.region() == 0)
   {
     throw cms::Exception("Geometry")
         << "GEMSynchronizer::simulateNoise() - this GEM id is from barrel, which cannot happen.";
   }
   const TrapezoidalStripTopology* top_(dynamic_cast<const TrapezoidalStripTopology*> (&(roll->topology())));
 
   // base_bottom, base_top, height, strips, pads (all half length)
   auto& parameters(roll->specs()->parameters());
   float semiBottomEdge(parameters[0]);
   float semiTopEdge(parameters[1]);
   float semiHeight(parameters[2]);
   float myTanPhi = (semiTopEdge - semiBottomEdge) / (semiHeight * 2);
   double rollRadius = top_->radius();
   const int nBxing(maxBunch_ - minBunch_ + 1);
   trArea = 2 * semiHeight * (semiTopEdge + semiBottomEdge);
 
   //simulate bkg contribution
   if (!doBkgNoise_)
     return;
 
   double aveNeutrRateBotRoll = 0.;
   double averageNoiseElectronRatePerRoll = 0.;
   int pdgid = 0;
 
   float myRand = CLHEP::RandFlat::shoot(engine);
   float yy_rand = 2 * semiHeight * (myRand - 0.5);
 
   double radius_rand = rollRadius + yy_rand;
 
     if(simulateElectronBkg_)
     averageNoiseElectronRatePerRoll = ME0ModElecBkgParam0
                       + ME0ModElecBkgParam1 * radius_rand
                       + ME0ModElecBkgParam2 * radius_rand * radius_rand
                       + ME0ModElecBkgParam3 * radius_rand * radius_rand * radius_rand
                       + ME0ModElecBkgParam4 * radius_rand * radius_rand * radius_rand * radius_rand
                       + ME0ModElecBkgParam5 * radius_rand * radius_rand * radius_rand * radius_rand * radius_rand
                       + ME0ModElecBkgParam6 * radius_rand * radius_rand * radius_rand * radius_rand * radius_rand * radius_rand;
 
   const double averageNoiseElec(averageNoiseElectronRatePerRoll * nBxing * bxwidth_ * trArea * 1.0e-9);
   const int n_elechits(CLHEP::RandPoissonQ::shoot(engine, averageNoiseElec));
 
   double xMax = semiTopEdge - (semiHeight - yy_rand) * myTanPhi;
   for (int i = 0; i < n_elechits; ++i)
   {
     //calculate xx_rand at a given yy_rand
     float myRandX = CLHEP::RandFlat::shoot(engine);
     float xx_rand = 2 * xMax * (myRandX - 0.5);
 
     float ex = sigma_u;
     float ey = sigma_v;
     float corr = 0.;
 
     GlobalPoint pointDigiHit = roll->toGlobal(LocalPoint(xx_rand, yy_rand));
     //calc tof to the random estimated point
     double stripRadius = sqrt(pointDigiHit.x() * pointDigiHit.x() + pointDigiHit.y() * pointDigiHit.y()
         + pointDigiHit.z() * pointDigiHit.z());
     double timeCalibrationOffset_ = (stripRadius * 1e+9) / (cspeed * 1e+2); //[ns]
     float tof = CLHEP::RandGaussQ::shoot(engine, timeCalibrationOffset_, sigma_t);
 
     float myrand = CLHEP::RandFlat::shoot(engine);
     if (myrand <= 0.5)
       pdgid = -11;
     else
       pdgid = 11;
 
     ME0DigiPreReco digi(xx_rand, yy_rand, ex, ey, corr, tof, pdgid);
     digi_.insert(digi);
   }
 
   aveNeutrRateBotRoll = ME0ModNeuBkgParam0
                       + ME0ModNeuBkgParam1 * radius_rand
                       + ME0ModNeuBkgParam2 * radius_rand * radius_rand
                       + ME0ModNeuBkgParam3 * radius_rand * radius_rand * radius_rand
                       + ME0ModNeuBkgParam4 * radius_rand * radius_rand * radius_rand * radius_rand
                       + ME0ModNeuBkgParam5 * radius_rand * radius_rand * radius_rand * radius_rand * radius_rand
                       + ME0ModNeuBkgParam6 * radius_rand * radius_rand * radius_rand * radius_rand * radius_rand * radius_rand;
 
   const double averageNoiseNeutral(aveNeutrRateBotRoll * nBxing * bxwidth_ * trArea * 1.0e-9);
   const int n_hits(CLHEP::RandPoissonQ::shoot(engine, averageNoiseNeutral));
 
   for (int i = 0; i < n_hits; ++i)
   {
     //calculate xx_rand at a given yy_rand
     float myRandX = CLHEP::RandFlat::shoot(engine);
     float xx_rand = 2 * xMax * (myRandX - 0.5);
 
     float ex = sigma_u;
     float ey = sigma_v;
     float corr = 0.;
 
     GlobalPoint pointDigiHit = roll->toGlobal(LocalPoint(xx_rand, yy_rand));
     //calc tof to the random estimated point
     double stripRadius = sqrt(pointDigiHit.x() * pointDigiHit.x() + pointDigiHit.y() * pointDigiHit.y()
         + pointDigiHit.z() * pointDigiHit.z());
     double timeCalibrationOffset_ = (stripRadius * 1e+9) / (cspeed * 1e+2); //[ns]
     float tof = CLHEP::RandGaussQ::shoot(engine, timeCalibrationOffset_, sigma_t);
 
     //distribute bkg between neutrons and gammas
     float myrand = CLHEP::RandFlat::shoot(engine);
     if (myrand <= 0.1)
       pdgid = 2112; // neutrons
     else
       pdgid = 22;
 
     ME0DigiPreReco digi(xx_rand, yy_rand, ex, ey, corr, tof, pdgid);
     digi_.insert(digi);
   }
 }

void ME0PreRecoGaussianModel::simulateSignal	(	const ME0EtaPartition *	roll,
		const edm::PSimHitContainer &	simHits,
		CLHEP::HepRandomEngine *	engine
	)

overridevirtual

Implements ME0DigiPreRecoModel.

Definition at line 55 of file ME0PreRecoGaussianModel.cc.

References funct::abs(), averageEfficiency_, corr, ME0DigiPreRecoModel::digi_, digitizeOnlyMuons_, sigma_t, sigma_u, sigma_v, x, and y.

 {
 for (const auto & hit: simHits)
 {
   if (std::abs(hit.particleType()) != 13 && digitizeOnlyMuons_) continue;
   // GEM efficiency
   if (CLHEP::RandFlat::shoot(engine) > averageEfficiency_) continue;
 
   auto entry = hit.entryPoint();
   float x=CLHEP::RandGaussQ::shoot(engine, entry.x(), sigma_u);
   float y=CLHEP::RandGaussQ::shoot(engine, entry.y(), sigma_v);
   float ex=sigma_u;
   float ey=sigma_v;
   float corr=0.;
   float tof=CLHEP::RandGaussQ::shoot(engine, hit.timeOfFlight(), sigma_t);
   int pdgid = hit.particleType();
   // please keep hit time always 0 for this model
   ME0DigiPreReco digi(x,y,ex,ey,corr,tof,pdgid);
   digi_.insert(digi);
 }
 }