#include <BTLTileDeviceSim.h>

Public Member Functions
	BTLTileDeviceSim (const edm::ParameterSet &pset)

void	getEvent (const edm::Event &evt)

void	getEventSetup (const edm::EventSetup &evt)

void	getHitsResponse (const std::vector< std::tuple< int, uint32_t, float > > &hitRefs, const edm::Handle< edm::PSimHitContainer > &hits, mtd_digitizer::MTDSimHitDataAccumulator simHitAccumulator, CLHEP::HepRandomEngine hre)

Private Attributes
const float	bxTime_

const MTDGeometry *	geom_

const float	LightCollEff_

const float	LightCollTime_

const float	LightYield_

const float	PDE_

const float	smearLightCollTime_

Detailed Description

Definition at line 21 of file BTLTileDeviceSim.h.

Constructor & Destructor Documentation

BTLTileDeviceSim::BTLTileDeviceSim ( const edm::ParameterSet & pset )

Definition at line 13 of file BTLTileDeviceSim.cc.

     : geom_(nullptr),
       bxTime_(pset.getParameter<double>("bxTime")),
       LightYield_(pset.getParameter<double>("LightYield")),
       LightCollEff_(pset.getParameter<double>("LightCollectionEff")),
       LightCollTime_(pset.getParameter<double>("LightCollectionTime")),
       smearLightCollTime_(pset.getParameter<double>("smearLightCollectionTime")),
       PDE_(pset.getParameter<double>("PhotonDetectionEff")) {}

Member Function Documentation

void BTLTileDeviceSim::getEvent ( const edm::Event & evt )

inline

Definition at line 25 of file BTLTileDeviceSim.h.

References hfClusterShapes_cfi::hits.

25 {}

void BTLTileDeviceSim::getEventSetup ( const edm::EventSetup & evt )

Definition at line 22 of file BTLTileDeviceSim.cc.

References relativeConstraints::geom, geom_, edm::EventSetup::get(), and edm::ESHandle< T >::product().

                                                              {
   edm::ESHandle<MTDGeometry> geom;
   evs.get<MTDDigiGeometryRecord>().get(geom);
   geom_ = geom.product();
 }

void BTLTileDeviceSim::getHitsResponse	(	const std::vector< std::tuple< int, uint32_t, float > > &	hitRefs,
		const edm::Handle< edm::PSimHitContainer > &	hits,
		mtd_digitizer::MTDSimHitDataAccumulator *	simHitAccumulator,
		CLHEP::HepRandomEngine *	hre
	)

Definition at line 28 of file BTLTileDeviceSim.cc.

References bxTime_, cuy::col, BTLDetId::column(), TauDecayModes::dec, PSimHit::energyLoss(), PSimHit::entryPoint(), Exception, DetId::Forward, geom_, MTDGeometry::idToDet(), BTLDetId::kTypeBoundariesReference, LightCollEff_, LightCollTime_, LightYield_, BTLDetId::modType(), BTLDetId::module(), MTDDetId::mtdRR(), MTDDetId::mtdSide(), RectangularMTDTopology::nrows(), PDE_, RectangularMTDTopology::pixel(), RectangularMTDTopology::pixelToModuleLocalPoint(), DetId::rawId(), BTLDetId::row(), smearLightCollTime_, ProxyMTDTopology::specificTopology(), and GeomDet::topology().

                                                                   {
   //loop over sorted simHits
   for (auto const& hitRef : hitRefs) {
     const int hitidx = std::get<0>(hitRef);
     const uint32_t id = std::get<1>(hitRef);
     const MTDDetId detId(id);
     const PSimHit& hit = hits->at(hitidx);
 
     // --- Safety check on the detector ID
     if (detId.det() != DetId::Forward || detId.mtdSubDetector() != 1)
       continue;
 
     if (id == 0)
       continue;  // to be ignored at RECO level
 
     BTLDetId btlid(detId);
     const int boundRef = BTLDetId::kTypeBoundariesReference[1];
     DetId geoId = BTLDetId(btlid.mtdSide(), btlid.mtdRR(), btlid.module() + boundRef * (btlid.modType() - 1), 0, 1);
     const MTDGeomDet* thedet = geom_->idToDet(geoId);
 
     if (thedet == nullptr) {
       throw cms::Exception("BTLTileDeviceSim") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                                << detId.rawId() << ") is invalid!" << std::dec << std::endl;
     }
     const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
     const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
     // calculate the simhit row and column
     const auto& pentry = hit.entryPoint();
     Local3DPoint simscaled(0.1 * pentry.x(), 0.1 * pentry.y(), 0.1 * pentry.z());  // mm -> cm here is the switch
     // translate from crystal-local coordinates to module-local coordinates to get the row and column
     simscaled = topo.pixelToModuleLocalPoint(simscaled, btlid.row(topo.nrows()), btlid.column(topo.nrows()));
     const auto& thepixel = topo.pixel(simscaled);
     uint8_t row(thepixel.first), col(thepixel.second);
 
     if (btlid.row(topo.nrows()) != row || btlid.column(topo.nrows()) != col) {
       edm::LogWarning("BTLTileDeviceSim")
           << "BTLDetId (row,column): (" << btlid.row(topo.nrows()) << ',' << btlid.column(topo.nrows())
           << ") is not equal to "
           << "topology (row,column): (" << uint32_t(row) << ',' << uint32_t(col) << "), overriding to detid";
       row = btlid.row(topo.nrows());
       col = btlid.column(topo.nrows());
     }
 
     // --- Store the detector element ID as a key of the MTDSimHitDataAccumulator map
     auto simHitIt =
         simHitAccumulator->emplace(mtd_digitizer::MTDCellId(id, row, col), mtd_digitizer::MTDCellInfo()).first;
 
     // --- Get the simHit energy and convert it from MeV to photo-electrons
     float Npe = 1000. * hit.energyLoss() * LightYield_ * LightCollEff_ * PDE_;
 
     // --- Get the simHit time of arrival and add the light collection time
     float toa = std::get<2>(hitRef) + LightCollTime_;
 
     if (smearLightCollTime_ > 0.)
       toa += CLHEP::RandGaussQ::shoot(hre, 0., smearLightCollTime_);
 
     if (toa > bxTime_ || toa < 0)  //just consider BX==0
       continue;
 
     (simHitIt->second).hit_info[0][0] += Npe;
 
     // --- Store the time of the first SimHit
     if ((simHitIt->second).hit_info[1][0] == 0 || toa < (simHitIt->second).hit_info[1][0])
       (simHitIt->second).hit_info[1][0] = toa;
 
   }  // hitRef loop
 }