Public Member Functions
unsigned int	findEtaPartition (float locY) const

float	getCentralTOF (const ME0DetId &me0Id, unsigned int partIdx) const

float	getPartCenter (const unsigned int partIdx) const

const TrapezoidalStripTopology *	getTopo (const unsigned int partIdx) const

unsigned int	numLayers () const

	TemporaryGeometry (const ME0Geometry *geometry, const unsigned int numberOfStrips, const unsigned int numberOfPartitions)

	~TemporaryGeometry ()

Private Member Functions
TrapezoidalStripTopology *	buildTopo (const std::vector< float > &_p) const

Private Attributes
float	middleDistanceFromBeam

std::vector< float >	partitionTops

std::vector< TrapezoidalStripTopology * >	stripTopos

std::vector< std::vector< double > >	tofs

Detailed Description

Definition at line 47 of file ME0ReDigiProducer.cc.

Constructor & Destructor Documentation

ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry	(	const ME0Geometry *	geometry,
		const unsigned int	numberOfStrips,
		const unsigned int	numberOfPartitions
	)

Definition at line 126 of file ME0ReDigiProducer.cc.

References buildTopo(), relativeConstraints::chamber, ME0Geometry::chambers(), chambers, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), Exception, getPartCenter(), createfilelist::int, LogDebug, middleDistanceFromBeam, ticl::constants::nLayers, ME0ReDigiProducer::numberOfPartitions, ME0ReDigiProducer::numberOfStrips, CalibrationSummaryClient_cfi::params, partitionTops, TrapezoidalStripTopology::radius(), TrapezoidalStripTopology::stripLength(), stripTopos, tofs, and PV3DBase< T, PVType, FrameType >::y().

                                                                                                {
   //First test geometry to make sure that it is compatible with our assumptions
   const auto& chambers = geometry->chambers();
   if (chambers.empty())
     throw cms::Exception("Setup")
         << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - No ME0Chambers in geometry.";
   const auto* mainChamber = chambers.front();
   const unsigned int nLayers = chambers.front()->nLayers();
   if (!nLayers)
     throw cms::Exception("Setup")
         << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - ME0Chamber has no layers.";
   const auto* mainLayer = mainChamber->layers()[0];
   if (!mainLayer->nEtaPartitions())
     throw cms::Exception("Setup")
         << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - ME0Layer has no partitions.";
   if (mainLayer->nEtaPartitions() != 1)
     throw cms::Exception("Setup") << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - This module is only "
                                      "compatitble with geometries that contain only one partition per ME0Layer.";
 
   const auto* mainPartition = mainLayer->etaPartitions()[0];
   const TrapezoidalStripTopology* mainTopo = dynamic_cast<const TrapezoidalStripTopology*>(&mainPartition->topology());
   if (!mainTopo)
     throw cms::Exception("Setup") << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - ME0 strip topology "
                                      "must be of type TrapezoidalStripTopology. This module cannot be used";
 
   for (const auto& chamber : geometry->chambers()) {
     if (chamber->nLayers() != int(nLayers))
       throw cms::Exception("Setup") << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - Not all "
                                        "ME0Chambers have the same number of layers. This module cannot be used.";
     for (unsigned int iL = 0; iL < nLayers; ++iL) {
       if (chamber->layers()[iL]->nEtaPartitions() != mainLayer->nEtaPartitions())
         throw cms::Exception("Setup") << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - Not all "
                                          "ME0Layers have the same number of partitions. This module cannot be used.";
       if (chamber->layers()[iL]->etaPartitions()[0]->specs()->parameters() != mainPartition->specs()->parameters())
         throw cms::Exception("Setup") << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - Not all ME0 ETA "
                                          "partitions have the same properties. This module cannot be used.";
       if (std::fabs(chamber->layers()[iL]->etaPartitions()[0]->position().z()) !=
           std::fabs(mainChamber->layers()[iL]->etaPartitions()[0]->position().z()))
         throw cms::Exception("Setup")
             << "ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry() - Not all ME0 ETA partitions in a single "
                "layer have the same Z position. This module cannot be used.";
     }
   }
 
   //Calculate radius to center of partition
   middleDistanceFromBeam = mainTopo->radius();
 
   //calculate the top of each eta partition, assuming equal distance in eta between partitions
   const auto localTop = LocalPoint(0, mainTopo->stripLength() / 2);
   const auto localBottom = LocalPoint(0, -1 * mainTopo->stripLength() / 2);
   const auto globalTop = mainPartition->toGlobal(localTop);
   const auto globalBottom = mainPartition->toGlobal(localBottom);
   const double etaTop = globalTop.eta();
   const double etaBottom = globalBottom.eta();
   const double zBottom = globalBottom.z();
 
   //Build topologies
   partitionTops.reserve(numberOfPartitions);
   stripTopos.reserve(numberOfPartitions);
   const auto& mainPars = mainPartition->specs()->parameters();
   for (unsigned int iP = 0; iP < numberOfPartitions; ++iP) {
     const double eta = (etaTop - etaBottom) * double(iP + 1) / double(numberOfPartitions) + etaBottom;
     const double distFromBeam = std::fabs(zBottom / std::sinh(eta));
     partitionTops.push_back(distFromBeam - middleDistanceFromBeam);
     LogDebug("ME0ReDigiProducer::TemporaryGeometry") << "Top of new partition: " << partitionTops.back() << std::endl;
 
     std::vector<float> params(4, 0);
 
     //half width of trapezoid at local coordinate Y
     auto getWidth = [&](float locY) -> float {
       return (mainPars[2] * (mainPars[1] + mainPars[0]) + locY * (mainPars[1] - mainPars[0])) / (2 * mainPars[2]);
     };
 
     params[0] = iP == 0 ? mainPars[0] : getWidth(partitionTops[iP - 1]);  // Half width of bottom of chamber
     params[1] =
         iP + 1 == numberOfPartitions ? mainPars[1] : getWidth(partitionTops[iP]);  // Half width of top of chamber
     params[2] = ((iP + 1 == numberOfPartitions ? localTop.y() : partitionTops[iP]) -
                  (iP == 0 ? localBottom.y() : partitionTops[iP - 1])) /
                 2;  // Half width of height of chamber
     params[3] = numberOfStrips;
 
     stripTopos.push_back(buildTopo(params));
   }
 
   //Get TOF at center of each partition
   tofs.resize(nLayers);
   LogDebug("ME0ReDigiProducer::TemporaryGeometry") << "TOF numbers [layer][partition]: ";
   for (unsigned int iL = 0; iL < nLayers; ++iL) {
     tofs[iL].resize(numberOfPartitions);
     for (unsigned int iP = 0; iP < numberOfPartitions; ++iP) {
       const LocalPoint partCenter(0., getPartCenter(iP), 0.);
       const GlobalPoint centralGP(mainChamber->layers()[iL]->etaPartitions()[0]->toGlobal(partCenter));
       tofs[iL][iP] = (centralGP.mag() / (CLHEP::c_light / CLHEP::cm));  //speed of light [cm/ns]
       LogDebug("ME0ReDigiProducer::TemporaryGeometry")
           << "[" << iL << "][" << iP << "]=" << tofs[iL][iP] << " " << std::endl;
     }
   }
 }

ME0ReDigiProducer::TemporaryGeometry::~TemporaryGeometry ( )

Definition at line 242 of file ME0ReDigiProducer.cc.

References AlCaHLTBitMon_ParallelJobs::p, and stripTopos.

                                                        {
   for (auto* p : stripTopos) {
     delete p;
   }
 }

Member Function Documentation

TrapezoidalStripTopology * ME0ReDigiProducer::TemporaryGeometry::buildTopo ( const std::vector< float > & _p ) const

private

Definition at line 248 of file ME0ReDigiProducer.cc.

References TtFullHadDaughter::B, b, h, LogDebug, and DigiDM_cff::strips.

Referenced by numLayers(), and TemporaryGeometry().

                                                                                                         {
   float b = _p[0];
   float B = _p[1];
   float h = _p[2];
   float r0 = h * (B + b) / (B - b);
   float striplength = h * 2;
   float strips = _p[3];
   float pitch = (b + B) / strips;
   int nstrip = static_cast<int>(strips);
 
   LogDebug("ME0ReDigiProducer::TemporaryGeometry")
       << "New partition parameters: "
       << "bottom width(" << 2 * b << ") top width(" << 2 * B << ") height(" << 2 * h << ") radius to center(" << r0
       << ") nStrips(" << strips << ") pitch(" << pitch << ")" << std::endl;
 
   return new TrapezoidalStripTopology(nstrip, pitch, striplength, r0);
 }

unsigned int ME0ReDigiProducer::TemporaryGeometry::findEtaPartition ( float locY ) const

Definition at line 227 of file ME0ReDigiProducer.cc.

References partitionTops, and stripTopos.

Referenced by ME0ReDigiProducer::getCustomStripProperties().

                                                                                   {
   unsigned int etaPart = stripTopos.size() - 1;
   for (unsigned int iP = 0; iP < stripTopos.size(); ++iP) {
     if (locY < partitionTops[iP]) {
       etaPart = iP;
       break;
     }
   }
   return etaPart;
 }

float ME0ReDigiProducer::TemporaryGeometry::getCentralTOF	(	const ME0DetId &	me0Id,
		unsigned int	partIdx
	)		const

inline

Definition at line 56 of file ME0ReDigiProducer.cc.

References ME0DetId::layer(), partIdx(), and tofs.

Referenced by ME0ReDigiProducer::getCustomStripProperties().

                                                                            {
       return tofs[me0Id.layer() - 1][partIdx];
     }  //in detId layer numbers stat at 1

float ME0ReDigiProducer::TemporaryGeometry::getPartCenter ( const unsigned int partIdx ) const

Definition at line 238 of file ME0ReDigiProducer.cc.

References middleDistanceFromBeam, partIdx(), and stripTopos.

Referenced by ME0ReDigiProducer::getCustomStripProperties(), getTopo(), and TemporaryGeometry().

                                                                                         {
   return stripTopos[partIdx]->radius() - middleDistanceFromBeam;
 }