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 37 of file ME0ReDigiProducer.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file ME0ReDigiProducer.cc.

References buildTopo(), relativeConstraints::chamber, ME0Geometry::chambers(), chambers, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), Exception, getPartCenter(), createfilelist::int, LogDebug, middleDistanceFromBeam, MuonTCMETValueMapProducer_cff::nLayers, ME0ReDigiProducer::numberOfPartitions, ME0ReDigiProducer::numberOfStrips, 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 116 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 120 of file ME0ReDigiProducer.cc.

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

Referenced by 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 106 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;
 }