Inheritance diagram for HGCalTestPartialWaferHits:

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

	~HGCalTestPartialWaferHits () 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 std::string	caloHitSource_

std::vector< int >	dumpDets_

const std::string	g4Label_

const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geomToken_

const std::string	missingFile_

const std::string	nameSense_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_calo_

std::vector< int >	wafers_

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 33 of file HGCalTestPartialWaferHits.cc.

Constructor & Destructor Documentation

◆ HGCalTestPartialWaferHits()

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

Definition at line 52 of file HGCalTestPartialWaferHits.cc.

References nameSense_.

     : g4Label_(ps.getParameter<std::string>("moduleLabel")),
       caloHitSource_(ps.getParameter<std::string>("caloHitSource")),
       nameSense_(ps.getParameter<std::string>("nameSense")),
       missingFile_(ps.getParameter<std::string>("missingFile")),
       tok_calo_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, caloHitSource_))),
       geomToken_(esConsumes<HGCalGeometry, IdealGeometryRecord>(edm::ESInputTag{"", nameSense_})) {
   edm::LogVerbatim("HGCalSim") << "Test Hit ID using SimHits for " << nameSense_ << " with module Label: " << g4Label_
                                << "   Hits: " << caloHitSource_ << " Missing Wafer file " << missingFile_;
   if (!missingFile_.empty()) {
     edm::FileInPath filetmp("SimG4CMS/Calo/data/" + missingFile_);
     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() > 2) {
           int layer = std::atoi(items[0].c_str());
           int waferU = std::atoi(items[1].c_str());
           int waferV = std::atoi(items[2].c_str());
           wafers_.emplace_back(HGCalWaferIndex::waferIndex(layer, waferU, waferV, false));
         } else if (items.size() == 1) {
           int dumpdet = std::atoi(items[0].c_str());
           dumpDets_.emplace_back(dumpdet);
           edm::LogVerbatim("HGCalSim") << nameSense_ << " Dump detector " << dumpdet;
         }
       }
       edm::LogVerbatim("HGCalSim") << "HGCalTestPartialWaferHits::Reads in " << wafers_.size() << ":"
                                    << dumpDets_.size() << " wafer|detector information from " << fileName;
       fInput.close();
     }
   }
 }

◆ ~HGCalTestPartialWaferHits()

HGCalTestPartialWaferHits::~HGCalTestPartialWaferHits ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 98 of file HGCalTestPartialWaferHits.cc.

References python.cmstools::all(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), TauDecayModes::dec, change_name::diff, dumpDets_, PVValHelper::dx, PVValHelper::dy, MillePedeFileConverter_cfg::e, spr::find(), relativeConstraints::geom, geomToken_, edm::EventSetup::getData(), caHitNtupletGeneratorKernels::good, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, hit::id, info(), edm::HandleBase::isValid(), nameSense_, TrackingDataMCValidation_Standalone_cff::nhits, mathSSE::sqrt(), tok_calo_, validateGeometry_cfg::valid, HGCalTypes::WaferFull, HGCalWaferIndex::waferIndex(), wafers_, geometryCSVtoXML::xx, and geometryCSVtoXML::xy.

                                                                                   {
   // get HGCalGeometry
   const HGCalGeometry* geom = &iS.getData(geomToken_);
   const HGCalDDDConstants& hgc = geom->topology().dddConstants();
   int firstLayer = hgc.getLayerOffset();
   // get the hit collection
   const edm::Handle<edm::PCaloHitContainer>& hitsCalo = e.getHandle(tok_calo_);
   bool getHits = (hitsCalo.isValid());
   uint32_t nhits = (getHits) ? hitsCalo->size() : 0;
   uint32_t good(0), allSi(0), all(0), allSc(0), bad(0);
   constexpr double tol = 2.0;
   edm::LogVerbatim("HGCalSim") << "HGCalTestPartialWaferHits: Input flags Hits " << getHits << " with " << nhits
                                << " hits: Layer Offset " << firstLayer;
 
   if (getHits) {
     std::vector<PCaloHit> hits;
     hits.insert(hits.end(), hitsCalo->begin(), hitsCalo->end());
     if (!hits.empty()) {
       // Loop over all hits
       for (auto hit : hits) {
         ++all;
         DetId id(hit.id());
         bool valid = (geom->topology()).valid(id);
         std::ostringstream st1;
         if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
           ++allSi;
           HGCSiliconDetId hid(id);
           st1 << hid;
           if (((id.det() == DetId::HGCalEE) && (nameSense_ == "HGCalEESensitive")) ||
               ((id.det() == DetId::HGCalHSi) && (nameSense_ == "HGCalHESiliconSensitive"))) {
             std::string_view pid =
                 ((hgc.cassetteShiftSilicon(hid.zside(), hid.layer(), hid.waferU(), hid.waferV())) ? "HGCSim"
                                                                                                   : "HGCalSim");
             const auto& info = hgc.waferInfo(hid.layer(), hid.waferU(), hid.waferV());
             if (!valid)
               st1 << " Wafer Type:Part:Orient:Cassette " << info.type << ":" << info.part << ":" << info.orient << ":"
                   << info.cassette;
             bool toCheck(false);
             if (!wafers_.empty()) {
               int indx = HGCalWaferIndex::waferIndex(firstLayer + hid.layer(), hid.waferU(), hid.waferV(), false);
               if (std::find(wafers_.begin(), wafers_.end(), indx) != wafers_.end())
                 toCheck = true;
             } else if (!dumpDets_.empty()) {
               if ((std::find(dumpDets_.begin(), dumpDets_.end(), static_cast<int>(id.det())) != dumpDets_.end()) &&
                   (info.part != HGCalTypes::WaferFull))
                 toCheck = true;
             } else {
               // Only partial wafers
               toCheck = (info.part != HGCalTypes::WaferFull);
             }
             if (toCheck) {
               ++good;
               GlobalPoint pos = geom->getPosition(id);
               bool valid1 = geom->topology().valid(id);
               bool valid2 = hgc.isValidHex8(hid.layer(), hid.waferU(), hid.waferV(), hid.cellU(), hid.cellV(), false);
               auto xy = hgc.locateCell(hid, false);
               double xx = (hid.zside() > 0) ? xy.first : -xy.first;
               double dx = xx - pos.x();
               double dy = xy.second - pos.y();
               double diff = std::sqrt(dx * dx + dy * dy);
               if ((diff > tol) || (!valid1) || (!valid2))
                 pid = "HGCalError";
               edm::LogVerbatim(pid) << "Hit[" << all << ":" << allSi << ":" << good << "]" << hid
                                     << " Wafer Type:Part:Orient:Cassette " << info.type << ":" << info.part << ":"
                                     << info.orient << ":" << info.cassette << " at (" << pos.x() << ":" << xx << ":"
                                     << dx << ", " << pos.y() << ":" << xy.second << ":" << dy << ", " << pos.z()
                                     << ") Valid " << valid1 << ":" << valid2 << " Distance " << diff;
             }
           }
         } else if (id.det() == DetId::HGCalHSc) {
           ++allSc;
           HGCScintillatorDetId hid(id);
           st1 << hid;
           if ((id.det() == DetId::HGCalHSc) && (nameSense_ == "HGCalHEScintillatorSensitive")) {
             std::string_view pid =
                 ((hgc.cassetteShiftScintillator(hid.zside(), hid.layer(), hid.iphi())) ? "HGCSim" : "HGCalSim");
             GlobalPoint pos = geom->getPosition(id);
             bool valid1 = geom->topology().valid(id);
             bool valid2 = hgc.isValidTrap(hid.zside(), hid.layer(), hid.ring(), hid.iphi());
             if ((!valid1) || (!valid2))
               pid = "HGCalError";
             int cassette = hgc.cassetteTile(hid.iphi());
             edm::LogVerbatim(pid) << "Hit[" << all << ":" << allSc << "] " << hid << " Cassette " << cassette << " at ("
                                   << pos.x() << ", " << pos.y() << ", " << pos.z() << ") Valid " << valid1 << ":"
                                   << valid2;
           }
         } else {
           st1 << std::hex << id.rawId() << std::dec;
         }
         if (!valid) {
           edm::LogVerbatim("HGCalError") << "Invalid ID " << st1.str();
           ++bad;
         }
       }
     }
   }
   edm::LogVerbatim("HGCalSim") << "Total hits = " << all << ":" << nhits << " Good Silicon DetIds = " << allSi << ":"
                                << good << " Scintitllator = " << allSc << " Invalid = " << bad;
 }

◆ beginJob()

void HGCalTestPartialWaferHits::beginJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 41 of file HGCalTestPartialWaferHits.cc.

41 {}

◆ endJob()

void HGCalTestPartialWaferHits::endJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 42 of file HGCalTestPartialWaferHits.cc.

42 {}

◆ fillDescriptions()

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

static

Definition at line 89 of file HGCalTestPartialWaferHits.cc.

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

                                                                                            {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("moduleLabel", "g4SimHits");
   desc.add<std::string>("caloHitSource", "HGCHitsEE");
   desc.add<std::string>("nameSense", "HGCalEESensitive");
   desc.add<std::string>("missingFile", "");
   descriptions.add("hgcalHitPartialEE", desc);
 }