#include <GEMSimAverage.h>

Inheritance diagram for GEMSimAverage:

Public Member Functions
	GEMSimAverage (const edm::ParameterSet &config)

void	setRandomEngine (CLHEP::HepRandomEngine &eng)

void	simulate (const GEMEtaPartition *roll, const edm::PSimHitContainer &rpcHits)

void	simulateNoise (const GEMEtaPartition *roll)

	~GEMSimAverage ()

Public Member Functions inherited from GEMSim
virtual void	fillDigis (int rollDetId, GEMDigiCollection &digis)

GEMSimSetUp *	getGEMSimSetUp ()

void	setGEMSimSetUp (GEMSimSetUp *setup)

const StripDigiSimLinks &	stripDigiSimLinks () const

virtual	~GEMSim ()

Private Member Functions
void	init ()

Private Attributes
double	averageEfficiency_

double	averageNoiseRate_

double	averageShapingTime_

double	bxwidth_

CLHEP::RandFlat *	flatDistr1_

CLHEP::RandFlat *	flatDistr2_

int	maxBunch_

int	minBunch_

CLHEP::RandPoissonQ *	poissonDistr_

GEMSynchronizer *	sync_

Additional Inherited Members
Public Types inherited from GEMSim
typedef edm::DetSet < StripDigiSimLink >	StripDigiSimLinks

Protected Types inherited from GEMSim
typedef std::multimap < std::pair< unsigned int, int > , const PSimHit *, std::less < std::pair< unsigned int, int > > >	DetectorHitMap

Protected Member Functions inherited from GEMSim
virtual void	addLinks (unsigned int strip, int bx)

	GEMSim (const edm::ParameterSet &config)

Protected Attributes inherited from GEMSim
DetectorHitMap	detectorHitMap_

GEMSimSetUp *	simSetUp_

StripDigiSimLinks	stripDigiSimLinks_

std::set< std::pair< int, int > >	strips_

Detailed Description

Class for the GEM strip response simulation that extends the Trivial model with average parameters for the GEM efficiency, timing and noise

Author: Sven Dildick

Definition at line 24 of file GEMSimAverage.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file GEMSimAverage.cc.

References averageEfficiency_, averageNoiseRate_, averageShapingTime_, bxwidth_, edm::ParameterSet::getParameter(), maxBunch_, minBunch_, and sync_.

                                                           :
   GEMSim(config)
 {
   averageEfficiency_ = config.getParameter<double>("averageEfficiency");
   averageShapingTime_ = config.getParameter<double>("averageShapingTime");
   averageNoiseRate_ = config.getParameter<double>("averageNoiseRate");
   bxwidth_ = config.getParameter<double>("bxwidth");
   minBunch_ = config.getParameter<int>("minBunch");
   maxBunch_ = config.getParameter<int>("maxBunch");
 
   sync_ = new GEMSynchronizer(config);
 }

GEMSimAverage::~GEMSimAverage ( )

Definition at line 44 of file GEMSimAverage.cc.

References flatDistr1_, flatDistr2_, poissonDistr_, and sync_.

 {
   if (flatDistr1_) delete flatDistr1_;
   if (flatDistr2_) delete flatDistr2_;
   if (poissonDistr_) delete poissonDistr_;
   delete sync_;
 }

Member Function Documentation

void GEMSimAverage::init ( void )

inlineprivatevirtual

Implements GEMSim.

Definition at line 40 of file GEMSimAverage.h.

40 {}

void GEMSimAverage::setRandomEngine ( CLHEP::HepRandomEngine & eng )

virtual

Implements GEMSim.

Definition at line 35 of file GEMSimAverage.cc.

References flatDistr1_, flatDistr2_, poissonDistr_, GEMSynchronizer::setRandomEngine(), and sync_.

 {
   flatDistr1_ = new CLHEP::RandFlat(eng);
   flatDistr2_ = new CLHEP::RandFlat(eng);
   poissonDistr_ = new CLHEP::RandPoissonQ(eng);
   sync_->setRandomEngine(eng);
 }

void GEMSimAverage::simulate	(	const GEMEtaPartition *	roll,
		const edm::PSimHitContainer &	rpcHits
	)

virtual

Implements GEMSim.

Definition at line 52 of file GEMSimAverage.cc.

References funct::abs(), averageEfficiency_, edm::DetSet< T >::clear(), GEMSim::detectorHitMap_, flatDistr1_, GEMSim::getGEMSimSetUp(), GEMSynchronizer::getSimHitBx(), GEMEtaPartition::id(), DetId::rawId(), GEMSynchronizer::setGEMSimSetUp(), GEMEtaPartition::specs(), GEMSim::stripDigiSimLinks_, GEMSim::strips_, sync_, and GEMEtaPartitionSpecs::topology().

 {
   sync_->setGEMSimSetUp(getGEMSimSetUp());
   stripDigiSimLinks_.clear();
   detectorHitMap_.clear();
   stripDigiSimLinks_ = StripDigiSimLinks(roll->id().rawId());
 
   const Topology& topology = roll->specs()->topology();
 
   for (const auto & hit: simHits)
   {
     if (std::abs(hit.particleType()) != 13) continue;
     // Check GEM efficiency
     if (flatDistr1_->fire(1) > averageEfficiency_) continue;
     auto entry = hit.entryPoint();
     
     int time_hit = sync_->getSimHitBx(&hit);
     std::pair<int, int> digi(topology.channel(entry) + 1, time_hit);
       
     detectorHitMap_.insert(DetectorHitMap::value_type(digi, &hit));
     strips_.insert(digi);
   }
 }

void GEMSimAverage::simulateNoise ( const GEMEtaPartition * roll )

virtual

Implements GEMSim.

Definition at line 78 of file GEMSimAverage.cc.

References averageNoiseRate_, bxwidth_, edm::hlt::Exception, flatDistr1_, flatDistr2_, i, GEMEtaPartition::id(), maxBunch_, minBunch_, GEMEtaPartition::nstrips(), poissonDistr_, GEMDetId::region(), GEMSim::strips_, GEMEtaPartition::topology(), x, SiStripMonitorClusterAlca_cfi::xmax, and SiStripMonitorClusterAlca_cfi::xmin.

 {
   GEMDetId gemId = roll->id();
   int nstrips = roll->nstrips();
   double area = 0.0;
   
   if ( gemId.region() == 0 )
     {
       throw cms::Exception("Geometry")
         << "GEMSynchronizer::simulateNoise() - this GEM id is from barrel, which cannot happen.";
     }
   else
     {
       const TrapezoidalStripTopology* top_=dynamic_cast<const TrapezoidalStripTopology*>(&(roll->topology()));
       float xmin = (top_->localPosition(0.)).x();
       float xmax = (top_->localPosition((float)roll->nstrips())).x();
       float striplength = (top_->stripLength());
       area = striplength*(xmax-xmin);
     }
   
   const int nBxing = maxBunch_ - minBunch_ + 1;
   double averageNoise = averageNoiseRate_ * nBxing * bxwidth_ * area * 1.0e-9;
 
   int n_hits = poissonDistr_->fire(averageNoise);
 
   for (int i = 0; i < n_hits; i++ ){
     int strip  = static_cast<int>(flatDistr1_->fire(1,nstrips));
     int time_hit = static_cast<int>(flatDistr2_->fire(nBxing)) + minBunch_;
     std::pair<int, int> digi(strip,time_hit);
     strips_.insert(digi);
   }
 }