#include <BTLDeviceSim.h>

Public Member Functions
	BTLDeviceSim (const edm::ParameterSet &pset, edm::ConsumesCollector iC)

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 edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord >	geomToken_

const float	LCEpositionSlope_

const float	LightCollEff_

const float	LightCollSlope_

const float	LightYield_

const float	PDE_

const MTDTopology *	topo_

const edm::ESGetToken< MTDTopology, MTDTopologyRcd >	topoToken_

Detailed Description

Definition at line 23 of file BTLDeviceSim.h.

Constructor & Destructor Documentation

◆ BTLDeviceSim()

BTLDeviceSim::BTLDeviceSim	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector	iC
	)

Definition at line 15 of file BTLDeviceSim.cc.

     : geomToken_(iC.esConsumes()),
       topoToken_(iC.esConsumes()),
       geom_(nullptr),
       topo_(nullptr),
       bxTime_(pset.getParameter<double>("bxTime")),
       LightYield_(pset.getParameter<double>("LightYield")),
       LightCollEff_(pset.getParameter<double>("LightCollectionEff")),
       LightCollSlope_(pset.getParameter<double>("LightCollectionSlope")),
       PDE_(pset.getParameter<double>("PhotonDetectionEff")),
       LCEpositionSlope_(pset.getParameter<double>("LCEpositionSlope")) {}

Member Function Documentation

◆ getEvent()

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

inline

Definition at line 27 of file BTLDeviceSim.h.

27 {}

◆ getEventSetup()

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

Definition at line 27 of file BTLDeviceSim.cc.

References geom_, geomToken_, edm::EventSetup::getData(), topo_, and topoToken_.

                                                          {
   geom_ = &evs.getData(geomToken_);
   topo_ = &evs.getData(topoToken_);
 }

◆ getHitsResponse()

void BTLDeviceSim::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 32 of file BTLDeviceSim.cc.

References bxTime_, cuy::col, BTLDetId::column(), angle_units::operators::convertGeVToMeV(), angle_units::operators::convertMmToCm(), MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, hcalRecHitTable_cff::detId, Exception, DetId::Forward, BTLDetId::geographicalId(), geom_, MTDTopology::getMTDTopologyMode(), hfClusterShapes_cfi::hits, MTDGeometry::idToDet(), mtd_digitizer::kInTimeBX, mtd_digitizer::kNumberOfBX, LCEpositionSlope_, LightCollEff_, LightCollSlope_, LightYield_, RectangularMTDTopology::nrows(), PDE_, RectangularMTDTopology::pitch(), RectangularMTDTopology::pixel(), RectangularMTDTopology::pixelToModuleLocalPoint(), position, DetId::rawId(), BTLDetId::row(), ProxyMTDTopology::specificTopology(), topo_, GeomDet::topology(), and hit::x.

                                                               {
   using namespace geant_units::operators;
 
   //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);
     DetId geoId = btlid.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topo_->getMTDTopologyMode()));
     const MTDGeomDet* thedet = geom_->idToDet(geoId);
 
     if (thedet == nullptr) {
       throw cms::Exception("BTLDeviceSim") << "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& position = hit.localPosition();
     Local3DPoint simscaled(convertMmToCm(position.x()), convertMmToCm(position.y()), convertMmToCm(position.z()));
     // 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("BTLDeviceSim") << "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 = convertGeVToMeV(hit.energyLoss()) * LightYield_ * LightCollEff_ * PDE_;
 
     // --- Calculate the light propagation time to the crystal bases (labeled L and R)
     double distR = 0.5 * topo.pitch().first - convertMmToCm(hit.localPosition().x());
     double distL = 0.5 * topo.pitch().first + convertMmToCm(hit.localPosition().x());
 
     double tR = std::get<2>(hitRef) + LightCollSlope_ * distR;
     double tL = std::get<2>(hitRef) + LightCollSlope_ * distL;
 
     // --- Accumulate in 15 buckets of 25ns (9 pre-samples, 1 in-time, 5 post-samples)
     const int iBXR = std::floor(tR / bxTime_) + mtd_digitizer::kInTimeBX;
     const int iBXL = std::floor(tL / bxTime_) + mtd_digitizer::kInTimeBX;
 
     // --- Right side
     if (iBXR > 0 && iBXR < mtd_digitizer::kNumberOfBX) {
       // Accumulate the energy of simHits in the same crystal
       (simHitIt->second).hit_info[0][iBXR] += Npe * (1. + LCEpositionSlope_ * convertMmToCm(hit.localPosition().x()));
 
       // Store the time of the first SimHit in the i-th BX
       if ((simHitIt->second).hit_info[1][iBXR] == 0 || tR < (simHitIt->second).hit_info[1][iBXR])
         (simHitIt->second).hit_info[1][iBXR] = tR - (iBXR - mtd_digitizer::kInTimeBX) * bxTime_;
     }
 
     // --- Left side
     if (iBXL > 0 && iBXL < mtd_digitizer::kNumberOfBX) {
       // Accumulate the energy of simHits in the same crystal
       (simHitIt->second).hit_info[2][iBXL] += Npe * (1. - LCEpositionSlope_ * convertMmToCm(hit.localPosition().x()));
 
       // Store the time of the first SimHit in the i-th BX
       if ((simHitIt->second).hit_info[3][iBXL] == 0 || tL < (simHitIt->second).hit_info[3][iBXL])
         (simHitIt->second).hit_info[3][iBXL] = tL - (iBXL - mtd_digitizer::kInTimeBX) * bxTime_;
     }
 
   }  // hitRef loop
 }