Inheritance diagram for HGCalTestGuardRing:

Public Member Functions
	HGCalTestGuardRing (const edm::ParameterSet &ps)

	~HGCalTestGuardRing () 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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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)

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

void	beginJob () override

void	endJob () override

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 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)

Private Attributes
const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geomToken_

const double	guardRingOffset_

const std::string	nameSense_

const std::string	waferFile_

std::map< HGCSiliconDetId, int >	waferID_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 37 of file HGCalTestGuardRing.cc.

Constructor & Destructor Documentation

◆ HGCalTestGuardRing()

HGCalTestGuardRing::HGCalTestGuardRing ( const edm::ParameterSet & ps )

Definition at line 55 of file HGCalTestGuardRing.cc.

References nameSense_.

     : nameSense_(ps.getParameter<std::string>("nameSense")),
       waferFile_(ps.getParameter<std::string>("waferFile")),
       guardRingOffset_(ps.getParameter<double>("guardRingOffset")),
       geomToken_(esConsumes<HGCalGeometry, IdealGeometryRecord>(edm::ESInputTag{"", nameSense_})) {
   DetId::Detector det = (nameSense_ != "HGCalHESiliconSensitive") ? DetId::HGCalEE : DetId::HGCalHSi;
   edm::LogVerbatim("HGCalSim") << "Test Guard Ring Offset " << guardRingOffset_ << " for " << nameSense_ << ":" << det
                                << " for wafers read from file " << waferFile_;
   if (!waferFile_.empty()) {
     std::string thick[4] = {"h120", "l200", "l300", "h200"};
     int addType[4] = {HGCalTypes::WaferFineThin,
                       HGCalTypes::WaferCoarseThin,
                       HGCalTypes::WaferCoarseThick,
                       HGCalTypes::WaferFineThick};
     const int partTypeH[6] = {HGCalTypes::WaferFull,
                               HGCalTypes::WaferHalf2,
                               HGCalTypes::WaferChopTwoM,
                               HGCalTypes::WaferSemi2,
                               HGCalTypes::WaferSemi2,
                               HGCalTypes::WaferFive2};
     const int partTypeL[7] = {HGCalTypes::WaferFull,
                               HGCalTypes::WaferHalf,
                               HGCalTypes::WaferHalf,
                               HGCalTypes::WaferSemi,
                               HGCalTypes::WaferSemi,
                               HGCalTypes::WaferFive,
                               HGCalTypes::WaferThree};
     edm::FileInPath filetmp("SimG4CMS/Calo/data/" + waferFile_);
     std::string fileName = filetmp.fullPath();
     std::ifstream fInput(fileName.c_str());
     if (!fInput.good()) {
       edm::LogVerbatim("HGCalSim") << "Cannot open file " << fileName;
     } else {
       char buffer[80];
       while (fInput.getline(buffer, 80)) {
         std::vector<std::string> items = CaloSimUtils::splitString(std::string(buffer));
         if (items.size() > 6) {
           int layer = std::atoi(items[0].c_str());
           int waferU = std::atoi(items[4].c_str());
           int waferV = std::atoi(items[5].c_str());
           int thck = static_cast<int>(std::find(thick, thick + 4, items[2]) - thick);
           int type = (thck < 4) ? addType[thck] : 0;
           HGCSiliconDetId id(det, -1, type, layer, waferU, waferV, 0, 0);
           int orient = std::atoi(items[5].c_str());
           int part = std::atoi(items[1].c_str());
           if (part >= 0) {
             if (type == HGCalTypes::WaferFineThin)
               part = partTypeH[part];
             else
               part = partTypeL[part];
           }
           waferID_[id] = orient * 100 + part;
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("HGCalSim") << "HGCalTestGuardRing::Reads " << id << " Orientation:Partial " << orient << ":"
                                        << part;
 #endif
         }
       }
       edm::LogVerbatim("HGCalSim") << "HGCalTestGuardRing::Reads in " << waferID_.size() << " wafer information from "
                                    << fileName;
       fInput.close();
     }
   }
 }

◆ ~HGCalTestGuardRing()

HGCalTestGuardRing::~HGCalTestGuardRing ( )

overridedefault

Member Function Documentation

◆ analyze()

void HGCalTestGuardRing::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	iS
	)

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 128 of file HGCalTestGuardRing.cc.

References python.cmstools::all(), dqmdumpme::first, relativeConstraints::geom, geomToken_, edm::EventSetup::getData(), caHitNtupletGeneratorKernels::good, guardRingOffset_, DetId::HGCalEE, DetId::HGCalHSi, HGCGuardRing::insidePolygon(), nano_mu_digi_cff::layer, nameSense_, caHitNtupletGeneratorKernels::nCells, hgcalPerformanceValidation::orient, edm::second(), findQualityFiles::v, validateGeometry_cfg::valid, HGCalTypes::waferFrontBack(), waferID_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), HGCalWaferMask::waferXY(), geometryCSVtoXML::xy, and ecaldqm::zside().

                                                                            {
   // get HGCalGeometry
   const HGCalGeometry* geom = &iS.getData(geomToken_);
   const HGCalDDDConstants& hgc = geom->topology().dddConstants();
   double waferSize = hgc.waferSize(false);
   HGCalCell wafer(waferSize, hgc.getUVMax(0), hgc.getUVMax(1));
   // get the hit collection
   edm::LogVerbatim("HGCalSim") << "HGCalTestGuardRing: Wafer Szie " << waferSize;
 
   // Loop over all IDs
   int all(0), allSi(0), good(0);
   for (std::map<HGCSiliconDetId, int>::const_iterator itr = waferID_.begin(); itr != waferID_.end(); ++itr) {
     HGCSiliconDetId id = itr->first;
     ++all;
     if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
       ++allSi;
       if (((id.det() == DetId::HGCalEE) && (nameSense_ == "HGCalEESensitive")) ||
           ((id.det() == DetId::HGCalHSi) && (nameSense_ == "HGCalHESiliconSensitive"))) {
         int partial = ((itr->second) % 100);
         int orient = (((itr->second) / 100) % 100);
         int type = id.type();
         int nCells = hgc.getUVMax(type);
         for (int u = 0; u < 2 * nCells; ++u) {
           for (int v = 0; v < 2 * nCells; ++v) {
             if (((v - u) < nCells) && ((u - v) <= nCells)) {
               HGCSiliconDetId hid(id.det(), id.zside(), id.type(), id.layer(), id.waferU(), id.waferV(), u, v);
               bool valid = (geom->topology()).valid(static_cast<DetId>(hid));
               if (valid) {
                 ++good;
                 int placeIndex = wafer.cellPlacementIndex(1, HGCalTypes::waferFrontBack(0), orient);
                 std::pair<double, double> xy = wafer.cellUV2XY1(u, v, placeIndex, type);
                 std::vector<std::pair<double, double> > wxy1 =
                     HGCalWaferMask::waferXY(partial, orient, -1, waferSize, 0.0, 0.0, 0.0);
                 bool check1 = HGCGuardRing::insidePolygon(xy.first, xy.second, wxy1);
                 std::ostringstream st1;
                 for (unsigned int k1 = 0; k1 < wxy1.size(); ++k1)
                   st1 << " (" << wxy1[k1].first << ", " << wxy1[k1].second << ")";
                 edm::LogVerbatim("HGCSim")
                     << "First " << hid << " Type:Partial:Orient:Place " << type << ":" << partial << ":" << orient
                     << ":" << placeIndex << " Boundary with " << wxy1.size() << " points: " << st1.str() << " check "
                     << check1 << " for (" << xy.first << ", " << xy.second << ")";
 
                 std::vector<std::pair<double, double> > wxy2 =
                     HGCalWaferMask::waferXY(partial, orient, -1, waferSize, guardRingOffset_, 0.0, 0.0);
                 bool check2 = HGCGuardRing::insidePolygon(xy.first, xy.second, wxy2);
                 std::ostringstream st2;
                 for (unsigned int k1 = 0; k1 < wxy2.size(); ++k1)
                   st2 << " (" << wxy2[k1].first << ", " << wxy2[k1].second << ")";
                 edm::LogVerbatim("HGCSim") << "Second Offset " << guardRingOffset_ << " Boundary with " << wxy2.size()
                                            << " points: " << st2.str() << " check " << check2 << " for (" << xy.first
                                            << ", " << xy.second << ")";
               }
             }
           }
         }
       }
     }
   }
   edm::LogVerbatim("HGCalSim") << "Total hits = " << all << " Good Silicon DetIds = " << allSi << ":" << good;
 }

◆ beginJob()

void HGCalTestGuardRing::beginJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 45 of file HGCalTestGuardRing.cc.

45 {}

◆ endJob()

void HGCalTestGuardRing::endJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 46 of file HGCalTestGuardRing.cc.

46 {}

◆ fillDescriptions()

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

static

Definition at line 120 of file HGCalTestGuardRing.cc.

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

                                                                                     {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("nameSense", "HGCalEESensitive");
   desc.add<std::string>("waferFile", "testWafersEE.txt");
   desc.add<double>("guardRingOffset", 1.0);
   descriptions.add("hgcalTestGuardRingEE", desc);
 }

Member Data Documentation

◆ geomToken_

const edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalTestGuardRing::geomToken_

private

Definition at line 51 of file HGCalTestGuardRing.cc.

Referenced by analyze().

◆ guardRingOffset_

const double HGCalTestGuardRing::guardRingOffset_

private

Definition at line 50 of file HGCalTestGuardRing.cc.

Referenced by analyze().

◆ nameSense_

const std::string HGCalTestGuardRing::nameSense_

private

Definition at line 49 of file HGCalTestGuardRing.cc.

Referenced by analyze(), and HGCalTestGuardRing().

◆ waferFile_

const std::string HGCalTestGuardRing::waferFile_

private

Definition at line 49 of file HGCalTestGuardRing.cc.

◆ waferID_

std::map<HGCSiliconDetId, int> HGCalTestGuardRing::waferID_

private

Definition at line 52 of file HGCalTestGuardRing.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalTestGuardRing()

◆ ~HGCalTestGuardRing()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

Member Data Documentation

◆ geomToken_

◆ guardRingOffset_

◆ nameSense_

◆ waferFile_

◆ waferID_