#include <Validation/HGCalValidation/plugins/HGCalWaferValidation.cc>

Inheritance diagram for HGCalWaferValidation:

Classes
struct	WaferInfo

Public Member Functions
	HGCalWaferValidation (const edm::ParameterSet &)

	~HGCalWaferValidation () override=default

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Types
using	WaferCoord = std::tuple< int, int, int >

using	WaferShapeMap = std::map< std::string, int >

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

bool	DDFindHGCal (DDCompactView::GraphWalker &walker, std::string targetName)

void	DDFindWafers (DDCompactView::GraphWalker &walker)

void	endJob () override

bool	isRotationMatched (const bool isNewFile, const int layer, const int fileShapeCode, const int geoRotCode, const int fileRotCode)

bool	isThicknessMatched (const int geoThickClass, const int fileThickness)

void	ProcessWaferLayer (DDCompactView::GraphWalker &walker)

std::string	strWaferCoord (const WaferCoord &coord)

Private Attributes
edm::FileInPath	geometryFileName_

unsigned int	layerCount_

std::vector< int >	layerTypes_

const std::string	logcat = "HGCalWaferValidation"

edm::ESGetToken< DDCompactView, IdealGeometryRecord >	viewToken_

std::map< WaferCoord, struct WaferInfo >	waferData_

const WaferShapeMap	waferShapeMapDD

const WaferShapeMap	waferShapeMapHD

const WaferShapeMap	waferShapeMapLD

std::map< WaferCoord, bool >	waferValidated_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Validates HGCal wafer data inside DD against specifications given in a flat text file.

Implementation: Uses GraphWalker to follow DD hierarchy to find HGCal EE module and the HE modules. Search of wafer layers and iterates each wafer found. Extract x, y coordinate position from wafer positioning; thickness, u & v coords from copyNo. Wafer shape and rotation are extracted from given names of wafer logical volumes. All extracted wafer info saved into a map indexed by (layer#, u, v). Each line in flat text file are compared against wafer information in the map. Any errors are reported, counted and summarized at the end. Unaccounted wafers, which are in DD but not in the flat text file, are also reported and counted.

Definition at line 72 of file HGCalWaferValidation.cc.

Member Typedef Documentation

◆ WaferCoord

using HGCalWaferValidation::WaferCoord = std::tuple<int, int, int>

private

Definition at line 84 of file HGCalWaferValidation.cc.

◆ WaferShapeMap

using HGCalWaferValidation::WaferShapeMap = std::map<std::string, int>

private

Definition at line 89 of file HGCalWaferValidation.cc.

Constructor & Destructor Documentation

◆ HGCalWaferValidation()

HGCalWaferValidation::HGCalWaferValidation ( const edm::ParameterSet & iConfig )

explicit

Definition at line 167 of file HGCalWaferValidation.cc.

References viewToken_.

     : geometryFileName_(iConfig.getParameter<edm::FileInPath>("GeometryFileName")) {
   viewToken_ = esConsumes<DDCompactView, IdealGeometryRecord>();
   //now do what ever initialization is needed
 }

◆ ~HGCalWaferValidation()

HGCalWaferValidation::~HGCalWaferValidation ( )

overridedefault

Member Function Documentation

◆ analyze()

void HGCalWaferValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 297 of file HGCalWaferValidation.cc.

References visDQMUpload::buf, DDFindHGCal(), DDFindWafers(), MillePedeFileConverter_cfg::fileName, edm::FileInPath::fullPath(), geometryFileName_, edm::EventSetup::getHandle(), mps_fire::i, isRotationMatched(), isThicknessMatched(), layerCount_, layerTypes_, logcat, HGCalWaferValidation::WaferInfo::rotCode, HGCalWaferValidation::WaferInfo::shapeCode, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, strWaferCoord(), HGCalWaferValidation::WaferInfo::thickClass, viewToken_, waferData_, HGCalWaferIndex::waferLayer(), HGCalWaferValidation::WaferInfo::waferName, HGCalTypes::WaferOut, waferShapeMapDD, waferShapeMapHD, waferShapeMapLD, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), waferValidated_, HGCalWaferValidation::WaferInfo::x, and HGCalWaferValidation::WaferInfo::y.

                                                                                       {
   using namespace edm;
 
   // Get the CMS DD
   auto viewH = iSetup.getHandle(viewToken_);
 
   if (!viewH.isValid()) {
     edm::LogPrint(logcat) << "Error obtaining geometry handle!";
     return;
   }
 
   edm::LogVerbatim(logcat) << "Root is : " << viewH->root();
   edm::LogVerbatim(logcat) << std::endl;
 
   // find HGCalEE
   auto eeWalker = viewH->walker();
   const bool eeFound = DDFindHGCal(eeWalker, "HGCalEE");
   if (eeFound) {
     edm::LogVerbatim(logcat) << "HGCalEE found!";
     edm::LogVerbatim(logcat) << "name     = " << eeWalker.current().first.name().name();
     edm::LogVerbatim(logcat) << "fullname = " << eeWalker.current().first.name().fullname();
   } else {
     edm::LogPrint(logcat) << "HGCalEE not found!";
   }
   edm::LogVerbatim(logcat) << std::endl;
 
   // find HGCalHEsil
   auto hesilWalker = viewH->walker();
   const bool hesilFound = DDFindHGCal(hesilWalker, "HGCalHEsil");
   if (hesilFound) {
     edm::LogVerbatim(logcat) << "HGCalHEsil found!";
     edm::LogVerbatim(logcat) << "name     = " << hesilWalker.current().first.name().name();
     edm::LogVerbatim(logcat) << "fullname = " << hesilWalker.current().first.name().fullname();
   } else {
     edm::LogPrint(logcat) << "HGCalHEsil not found!";
   }
   edm::LogVerbatim(logcat) << std::endl;
 
   // find HGCalHEmix
   auto hemixWalker = viewH->walker();
   const bool hemixFound = DDFindHGCal(hemixWalker, "HGCalHEmix");
   if (hemixFound) {
     edm::LogVerbatim(logcat) << "HGCalHEmix found!";
     edm::LogVerbatim(logcat) << "name     = " << hemixWalker.current().first.name().name();
     edm::LogVerbatim(logcat) << "fullname = " << hemixWalker.current().first.name().fullname();
   } else {
     edm::LogPrint(logcat) << "HGCalHEmix not found!";
   }
   edm::LogVerbatim(logcat) << std::endl;
 
   // give up if no HGCal found at all
   if (!(eeFound || hesilFound || hemixFound)) {
     edm::LogPrint(logcat) << "Nothing found. Giving up.";
     return;
   }
 
   // Now walk the HGCalEE walker to find the first wafer on each layer and process them
   edm::LogVerbatim(logcat) << "Calling DDFindWafers(eeWalker);";
   DDFindWafers(eeWalker);
 
   // Walk the HGCalHEsilwalker to find the first wafer on each layer and process them
   edm::LogVerbatim(logcat) << "Calling DDFindWafers(hesilWalker);";
   DDFindWafers(hesilWalker);
 
   // Walk the HGCalHEmix walker to find the first wafer on each layer and process them
   edm::LogVerbatim(logcat) << "Calling DDFindWafers(hemixWalker);";
   DDFindWafers(hemixWalker);
 
   // Confirm all the DD wafers have been read
   edm::LogVerbatim(logcat) << "Number of wafers read from DD: " << waferData_.size();
 
   // Now open the geometry text file
   std::string fileName = geometryFileName_.fullPath();
   edm::LogVerbatim(logcat) << "Opening geometry text file: " << fileName;
   std::ifstream geoTxtFile(fileName);
 
   if (!geoTxtFile) {
     edm::LogPrint(logcat) << "Cannot open geometry text file.";
     return;
   }
 
   // total processed counter
   int nTotalProcessed = 0;
 
   // geometry error counters
   int nMissing = 0;
   int nThicknessError = 0;
   int nPosXError = 0;
   int nPosYError = 0;
   int nShapeError = 0;
   int nRotError = 0;
   int nUnaccounted = 0;
 
   std::string buf;
 
   // find out if this file is an old file or a new file
   std::getline(geoTxtFile, buf);
   std::stringstream ss(buf);
   std::vector<std::string> first_tokens;
   while (ss >> buf)
     if (!buf.empty())
       first_tokens.push_back(buf);
 
   const bool isNewFile(first_tokens.size() == 1);
 
   if (isNewFile) {
     edm::LogVerbatim(logcat) << "Text file is of newer version.";
     layerCount_ = std::stoi(buf);
     std::getline(geoTxtFile, buf);
     std::stringstream layerTypesSS(buf);
     while (layerTypesSS >> buf)
       if (!buf.empty())
         layerTypes_.push_back(std::stoi(buf));
     if (layerTypes_.size() != layerCount_)
       edm::LogWarning(logcat) << "Number of layer types does not tally with layer count.";
 
     // TEMP: make sure reading is correct
     edm::LogVerbatim(logcat) << "layerCount = " << layerCount_;
     for (unsigned i = 0; i < layerTypes_.size(); i++) {
       edm::LogVerbatim(logcat) << "  layerType " << i + 1 << " = " << layerTypes_[i];  // 1-based
     }
   } else {
     layerCount_ = 0;
     // rewind back
     geoTxtFile.clear();
     geoTxtFile.seekg(0);
   }
 
   // process each line on the text file
   while (std::getline(geoTxtFile, buf)) {
     std::stringstream ss(buf);
     std::vector<std::string> tokens;
     while (ss >> buf)
       if (!buf.empty())
         tokens.push_back(buf);
     if (tokens.size() != 8)
       continue;
 
     nTotalProcessed++;
 
     // extract wafer info from a textfile line
     const int waferLayer(std::stoi(tokens[0]));
     const std::string waferShapeStr(tokens[1]);
     const std::string waferDensityStr(isNewFile ? tokens[2].substr(0, 1) : "");
     const int waferThickness(isNewFile ? std::stoi(tokens[2].substr(1)) : std::stoi(tokens[2]));
     const double waferX(std::stod(tokens[3]));
     const double waferY(std::stod(tokens[4]));
     const int waferRotCode(std::stoi(tokens[5]));
     const int waferU(std::stoi(tokens[6]));
     const int waferV(std::stoi(tokens[7]));
     const int waferShapeCode(isNewFile ? (waferDensityStr == "l"   ? waferShapeMapLD.at(waferShapeStr)
                                           : waferDensityStr == "h" ? waferShapeMapHD.at(waferShapeStr)
                                                                    : HGCalTypes::WaferOut)
                                        : waferShapeMapDD.at(waferShapeStr));
 
     // map index for crosschecking with DD
     const WaferCoord waferCoord(waferLayer, waferU, waferV);
 
     // now check for (and report) wafer data disagreements
 
     if (waferData_.find(waferCoord) == waferData_.end()) {
       nMissing++;
       edm::LogVerbatim(logcat) << "MISSING: " << strWaferCoord(waferCoord);
       continue;
     }
 
     const struct WaferInfo waferInfo = waferData_[waferCoord];
     waferValidated_[waferCoord] = true;
 
     if (!isThicknessMatched(waferInfo.thickClass, waferThickness)) {
       nThicknessError++;
       edm::LogVerbatim(logcat) << "THICKNESS ERROR: " << strWaferCoord(waferCoord);
     }
 
     // it seems that wafer x-coords relative to their layer plane are mirrored...
     if (fabs(-waferInfo.x - waferX) > 0.015) {  // assuming this much tolerance
       nPosXError++;
       edm::LogVerbatim(logcat) << "POSITION x ERROR: " << strWaferCoord(waferCoord);
     }
 
     if (fabs(waferInfo.y - waferY) > 0.015) {  // assuming this much tolerance
       nPosYError++;
       edm::LogVerbatim(logcat) << "POSITION y ERROR: " << strWaferCoord(waferCoord);
     }
 
     if (waferInfo.shapeCode != waferShapeCode || waferShapeCode == HGCalTypes::WaferOut) {
       nShapeError++;
       edm::LogVerbatim(logcat) << "SHAPE ERROR: " << strWaferCoord(waferCoord) << "  ( " << waferInfo.shapeCode
                                << " != " << waferDensityStr << waferShapeCode << " )  name=" << waferInfo.waferName;
     }
 
     if (!isRotationMatched(isNewFile, std::get<0>(waferCoord), waferShapeCode, waferInfo.rotCode, waferRotCode)) {
       nRotError++;
       edm::LogVerbatim(logcat) << "ROTATION ERROR: " << strWaferCoord(waferCoord) << "  ( " << waferInfo.rotCode
                                << " != " << waferRotCode << " (" << waferShapeCode
                                << ") )  name=" << waferInfo.waferName;
     }
   }
 
   geoTxtFile.close();
 
   // Find unaccounted DD wafers
   for (auto const& accounted : waferValidated_) {
     if (!accounted.second) {
       nUnaccounted++;
       edm::LogVerbatim(logcat) << "UNACCOUNTED: " << strWaferCoord(accounted.first);
     }
   }
 
   // Print out error counts
   edm::LogVerbatim(logcat) << std::endl;
   edm::LogVerbatim(logcat) << "*** ERROR COUNTS ***";
   edm::LogVerbatim(logcat) << "Missing         :  " << nMissing;
   edm::LogVerbatim(logcat) << "Thickness error :  " << nThicknessError;
   edm::LogVerbatim(logcat) << "Pos-x error     :  " << nPosXError;
   edm::LogVerbatim(logcat) << "Pos-y error     :  " << nPosYError;
   edm::LogVerbatim(logcat) << "Shape error     :  " << nShapeError;
   edm::LogVerbatim(logcat) << "Rotation error  :  " << nRotError;
   edm::LogVerbatim(logcat) << "Unaccounted     :  " << nUnaccounted;
   edm::LogVerbatim(logcat) << std::endl;
   edm::LogVerbatim(logcat) << "Total wafers processed from geotxtfile = " << nTotalProcessed;
   edm::LogVerbatim(logcat) << std::endl;
 
   // Issue a LogPrint (warning) if there is at least one wafer errors
   if (nMissing > 0 || nThicknessError > 0 || nPosXError > 0 || nPosYError > 0 || nShapeError > 0 || nRotError > 0 ||
       nUnaccounted > 0) {
     edm::LogPrint(logcat) << "There are at least one wafer error.";
   }
 }

◆ beginJob()

void HGCalWaferValidation::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 528 of file HGCalWaferValidation.cc.

                                     {
   // please remove this method if not needed
 }

◆ DDFindHGCal()

bool HGCalWaferValidation::DDFindHGCal	(	DDCompactView::GraphWalker &	walker,
		std::string	targetName
	)

private

Definition at line 185 of file HGCalWaferValidation.cc.

References math::GraphWalker< N, E >::current(), math::GraphWalker< N, E >::firstChild(), math::GraphWalker< N, E >::nextSibling(), and math::GraphWalker< N, E >::parent().

Referenced by analyze().

                                                                                              {
   if (walker.current().first.name().name() == targetName) {
     // target found
     return true;
   }
   if (walker.firstChild()) {
     do {
       if (DDFindHGCal(walker, targetName))
         // target inside child
         return true;
     } while (walker.nextSibling());
     walker.parent();
   }
   return false;
 }

◆ DDFindWafers()

void HGCalWaferValidation::DDFindWafers ( DDCompactView::GraphWalker & walker )

private

Definition at line 202 of file HGCalWaferValidation.cc.

References math::GraphWalker< N, E >::current(), math::GraphWalker< N, E >::firstChild(), math::GraphWalker< N, E >::nextSibling(), math::GraphWalker< N, E >::parent(), and ProcessWaferLayer().

Referenced by analyze().

                                                                         {
   if (walker.current().first.name().fullname().rfind("hgcalwafer:", 0) == 0) {
     // first wafer found. Now process the entire layer of wafers.
     ProcessWaferLayer(walker);
     return;
   }
   if (walker.firstChild()) {
     do {
       DDFindWafers(walker);
     } while (walker.nextSibling());
     walker.parent();
   }
 }

◆ endJob()

void HGCalWaferValidation::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 533 of file HGCalWaferValidation.cc.

                                   {
   // please remove this method if not needed
 }

◆ fillDescriptions()

void HGCalWaferValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 538 of file HGCalWaferValidation.cc.

References edm::ConfigurationDescriptions::add(), and submitPVResolutionJobs::desc.

                                                                                       {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.add<edm::FileInPath>("GeometryFileName",
                             edm::FileInPath("Validation/HGCalValidation/data/geomnew_corrected_360.txt"));
   descriptions.add("hgcalWaferValidation", desc);
 }

◆ isRotationMatched()

bool HGCalWaferValidation::isRotationMatched	(	const bool	isNewFile,
		const int	layer,
		const int	fileShapeCode,
		const int	geoRotCode,
		const int	fileRotCode
	)

private

Definition at line 280 of file HGCalWaferValidation.cc.

References pixelTopology::layer, layerTypes_, and HGCalTypes::WaferFull.

Referenced by analyze().

                                                                                                                  {
   if (fileShapeCode != HGCalTypes::WaferFull && geoRotCode == fileRotCode)
     return true;
   if (fileShapeCode == HGCalTypes::WaferFull) {
     if (isNewFile && layerTypes_[layer - 1] == 3) {  // this array is index-0 based
       if ((geoRotCode + 1) % 2 == fileRotCode % 2)
         return true;
     } else {
       if (geoRotCode % 2 == fileRotCode % 2)
         return true;
     }
   }
   return false;
 }

◆ isThicknessMatched()

bool HGCalWaferValidation::isThicknessMatched	(	const int	geoThickClass,
		const int	fileThickness
	)

private

Definition at line 269 of file HGCalWaferValidation.cc.

Referenced by analyze().

                                                                                               {
   if (geoThickClass == 0 && fileThickness == 120)
     return true;
   if (geoThickClass == 1 && fileThickness == 200)
     return true;
   if (geoThickClass == 2 && fileThickness == 300)
     return true;
   return false;
 }

◆ ProcessWaferLayer()

void HGCalWaferValidation::ProcessWaferLayer ( DDCompactView::GraphWalker & walker )

private

Definition at line 217 of file HGCalWaferValidation.cc.

References math::GraphWalker< N, E >::current(), HGCalTypes::getUnpackedType(), HGCalTypes::getUnpackedU(), HGCalTypes::getUnpackedV(), logcat, math::GraphWalker< N, E >::nextSibling(), HGCalWaferValidation::WaferInfo::rotCode, HGCalWaferValidation::WaferInfo::shapeCode, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, HGCalWaferValidation::WaferInfo::thickClass, waferData_, HGCalWaferIndex::waferLayer(), HGCalWaferValidation::WaferInfo::waferName, waferShapeMapDD, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), waferValidated_, HGCalWaferValidation::WaferInfo::x, and HGCalWaferValidation::WaferInfo::y.

Referenced by DDFindWafers().

                                                                              {
   int waferLayer = 0;  // layer numbers in DD are assumed to be sequential from 1
   waferLayer++;
   edm::LogVerbatim(logcat) << "ProcessWaferLayer: Processing layer " << waferLayer;
   do {
     if (walker.current().first.name().fullname().rfind("hgcalwafer:", 0) == 0) {
       auto wafer = walker.current();
       const std::string waferName(walker.current().first.name().fullname());
       //edm::LogVerbatim(logcat) << "  " << waferName; // DEBUG: in case wafer name info is needed
       const int copyNo = wafer.second->copyno();
       // extract DD layer properties
       const int waferType = HGCalTypes::getUnpackedType(copyNo);
       const int waferU = HGCalTypes::getUnpackedU(copyNo);
       const int waferV = HGCalTypes::getUnpackedV(copyNo);
       const WaferCoord waferCoord(waferLayer, waferU, waferV);  // map index
       // build struct of DD wafer properties
       struct WaferInfo waferInfo;
       waferInfo.waferName = waferName;  //TEMPORARY
       waferInfo.thickClass = waferType;
       waferInfo.x = wafer.second->translation().x();
       waferInfo.y = wafer.second->translation().y();
       const std::string waferNameData =
           std::regex_replace(waferName,
                              std::regex("(HGCal[EH]E)(Wafer[01])(Fine|Coarse[12])([a-z]*)([0-9]*)"),
                              "$1 $2-$3 $4 $5",
                              std::regex_constants::format_no_copy);
       std::stringstream ss(waferNameData);
       std::string EEorHE;
       std::string typeStr;
       std::string shapeStr;
       std::string rotStr;
       ss >> EEorHE >> typeStr >> shapeStr >> rotStr;
       // assume rotational symmetry of full-sized wafers
       if (shapeStr.empty())
         shapeStr = "F";
       if (rotStr.empty())
         rotStr = "0";
       const int rotCode(std::stoi(rotStr));
       //edm::LogVerbatim(logcat) << "rotStr " << rotStr << " rotCode " << rotCode;
 
       // convert shape code to wafer types defined in HGCalTypes.h
       waferInfo.shapeCode = waferShapeMapDD.at(shapeStr);
 
       waferInfo.rotCode = rotCode;
       // populate the map
       waferData_[waferCoord] = waferInfo;
       waferValidated_[waferCoord] = false;
     }
   } while (walker.nextSibling());
 }

◆ strWaferCoord()

std::string HGCalWaferValidation::strWaferCoord ( const WaferCoord & coord )

private

Definition at line 178 of file HGCalWaferValidation.cc.

References contentValuesCheck::ss.

Referenced by analyze().

                                                                      {
   std::stringstream ss;
   ss << "(" << std::get<0>(coord) << "," << std::get<1>(coord) << "," << std::get<2>(coord) << ")";
   return ss.str();
 }

Member Data Documentation

◆ geometryFileName_

edm::FileInPath HGCalWaferValidation::geometryFileName_

private

Definition at line 130 of file HGCalWaferValidation.cc.

Referenced by analyze().

◆ layerCount_

unsigned int HGCalWaferValidation::layerCount_

private

Definition at line 133 of file HGCalWaferValidation.cc.

Referenced by analyze().

◆ layerTypes_

std::vector<int> HGCalWaferValidation::layerTypes_

private

Definition at line 134 of file HGCalWaferValidation.cc.

Referenced by analyze(), and isRotationMatched().

◆ logcat

const std::string HGCalWaferValidation::logcat = "HGCalWaferValidation"

private

Definition at line 81 of file HGCalWaferValidation.cc.

Referenced by analyze(), and ProcessWaferLayer().

◆ viewToken_

edm::ESGetToken<DDCompactView, IdealGeometryRecord> HGCalWaferValidation::viewToken_

private

Definition at line 147 of file HGCalWaferValidation.cc.

Referenced by analyze(), and HGCalWaferValidation().

◆ waferData_

std::map<WaferCoord, struct WaferInfo> HGCalWaferValidation::waferData_

private

Definition at line 150 of file HGCalWaferValidation.cc.

Referenced by analyze(), and ProcessWaferLayer().

◆ waferShapeMapDD

const WaferShapeMap HGCalWaferValidation::waferShapeMapDD

private

Initial value:

= {{"F", HGCalTypes::WaferFull},
                                         {"a", HGCalTypes::WaferHalf},
                                         {"am", HGCalTypes::WaferHalf2},
                                         {"b", HGCalTypes::WaferFive},
                                         {"bm", HGCalTypes::WaferFive2},
                                         {"c", HGCalTypes::WaferThree},
                                         {"d", HGCalTypes::WaferSemi},
                                         {"dm", HGCalTypes::WaferSemi2},
                                         {"g", HGCalTypes::WaferChopTwo},
                                         {"gm", HGCalTypes::WaferChopTwoM}}

Definition at line 91 of file HGCalWaferValidation.cc.

Referenced by analyze(), and ProcessWaferLayer().

◆ waferShapeMapHD

const WaferShapeMap HGCalWaferValidation::waferShapeMapHD

private

Initial value:

= {{"0", HGCalTypes::WaferFull},
                                         {"1", HGCalTypes::WaferHalf2},
                                         {"2", HGCalTypes::WaferChopTwoM},
                                         {"3", HGCalTypes::WaferSemi2},
                                         {"4", HGCalTypes::WaferSemi2},
                                         {"5", HGCalTypes::WaferFive2}}

Definition at line 110 of file HGCalWaferValidation.cc.

Referenced by analyze().

◆ waferShapeMapLD

const WaferShapeMap HGCalWaferValidation::waferShapeMapLD

private

Initial value:

= {{"0", HGCalTypes::WaferFull},
                                         {"1", HGCalTypes::WaferHalf},
                                         {"2", HGCalTypes::WaferHalf},
                                         {"3", HGCalTypes::WaferSemi},
                                         {"4", HGCalTypes::WaferSemi},
                                         {"5", HGCalTypes::WaferFive},
                                         {"6", HGCalTypes::WaferThree}}

Definition at line 102 of file HGCalWaferValidation.cc.

Referenced by analyze().

◆ waferValidated_

std::map<WaferCoord, bool> HGCalWaferValidation::waferValidated_

private

Definition at line 153 of file HGCalWaferValidation.cc.

Referenced by analyze(), and ProcessWaferLayer().

Classes

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ WaferCoord

◆ WaferShapeMap

Constructor & Destructor Documentation

◆ HGCalWaferValidation()

◆ ~HGCalWaferValidation()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ DDFindHGCal()

◆ DDFindWafers()

◆ endJob()

◆ fillDescriptions()

◆ isRotationMatched()

◆ isThicknessMatched()

◆ ProcessWaferLayer()

◆ strWaferCoord()

Member Data Documentation

◆ geometryFileName_

◆ layerCount_

◆ layerTypes_

◆ logcat

◆ viewToken_

◆ waferData_

◆ waferShapeMapDD

◆ waferShapeMapHD

◆ waferShapeMapLD

◆ waferValidated_