HGCalTestPartialWaferRecHits Class Reference

Inheritance diagram for HGCalTestPartialWaferRecHits:

Public Member Functions
	HGCalTestPartialWaferRecHits (const edm::ParameterSet &ps)
	~HGCalTestPartialWaferRecHits () 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)

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

Private Attributes
std::vector< int >	dumpDets_
const std::string	missingFile_
const std::string	nameDetector_
const edm::EDGetTokenT< HGCRecHitCollection >	recHitSource_
const edm::InputTag	source_
const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	tok_hgcGeom_
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 30 of file HGCalTestPartialWaferRecHits.cc.

Constructor & Destructor Documentation

HGCalTestPartialWaferRecHits()

HGCalTestPartialWaferRecHits::HGCalTestPartialWaferRecHits

(

const edm::ParameterSet &

ps

)

Definition at line 50 of file HGCalTestPartialWaferRecHits.cc.

References nameDetector_.

    : nameDetector_(ps.getParameter<std::string>("detectorName")),
      missingFile_(ps.getParameter<std::string>("missingFile")),
      source_(ps.getParameter<edm::InputTag>("source")),
      recHitSource_(consumes<HGCRecHitCollection>(source_)),
      tok_hgcGeom_(esConsumes<HGCalGeometry, IdealGeometryRecord>(edm::ESInputTag{"", nameDetector_})) {
  edm::LogVerbatim("HGCalSim") << "Test Hit ID using RecHits for " << nameDetector_ << " with module Label: " << source_
                               << " 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") << nameDetector_ << " Dump detector " << dumpdet;
        }
      }
      edm::LogVerbatim("HGCalSim") << "HGCalTestPartialWaferRecHits::Reads in " << wafers_.size() << ":"
                                   << dumpDets_.size() << " wafer|detector information from " << fileName;
      fInput.close();
    }
  }
}

~HGCalTestPartialWaferRecHits()

HGCalTestPartialWaferRecHits::~HGCalTestPartialWaferRecHits ( )

overridedefault

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 94 of file HGCalTestPartialWaferRecHits.cc.

References python.cmstools::all(), HGCSiliconDetId::cellU(), HGCSiliconDetId::cellV(), dumpDets_, MillePedeFileConverter_cfg::e, spr::find(), relativeConstraints::geom, edm::EventSetup::getData(), good, DetId::HGCalEE, DetId::HGCalHSi, l1ctLayer2EG_cff::id, info(), edm::HandleBase::isValid(), HGCSiliconDetId::layer(), nhits, edm::Handle< T >::product(), recHitSource_, edm::SortedCollection< T, SORT >::size(), tok_hgcGeom_, HGCalTypes::WaferFull, HGCalWaferIndex::waferIndex(), wafers_, HGCSiliconDetId::waferU(), and HGCSiliconDetId::waferV().

                                                                                     {
  // get hgcalGeometry
  const HGCalGeometry* geom = &iS.getData(tok_hgcGeom_);
  const HGCalDDDConstants& hgc = geom->topology().dddConstants();
  int firstLayer = hgc.getLayerOffset();
  // get the hit collection
  const edm::Handle<HGCRecHitCollection>& theRecHitContainers = e.getHandle(recHitSource_);
  bool getHits = (theRecHitContainers.isValid());
  uint32_t nhits = (getHits) ? theRecHitContainers->size() : 0;
  uint32_t good(0), allSi(0), all(0);
  edm::LogVerbatim("HGCalSim") << "HGCalTestPartialWaferRecHits: Input flags Hits " << getHits << " with " << nhits
                               << " hits: Layer Offset " << firstLayer;

  if (getHits) {
    // Loop over all hits
    for (const auto& it : *(theRecHitContainers.product())) {
      ++all;
      DetId id(it.id());
      if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
        ++allSi;
        HGCSiliconDetId hid(id);
        const auto& info = hgc.waferInfo(hid.layer(), hid.waferU(), hid.waferV());
        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());
          edm::LogVerbatim("HGCalSim") << "Hit[" << all << ":" << allSi << ":" << good << "]" << hid
                                       << " Wafer Type:Part:Orient:Cassette " << info.type << ":" << info.part << ":"
                                       << info.orient << ":" << info.cassette << " at (" << pos.x() << ", " << pos.y()
                                       << ", " << pos.z() << ") Validity " << valid1 << ":" << valid2;
        }
      }
    }
  }
  edm::LogVerbatim("HGCalSim") << "Total hits = " << all << ":" << nhits << " Good DetIds = " << allSi << ":" << good;
}

beginJob()

void HGCalTestPartialWaferRecHits::beginJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 38 of file HGCalTestPartialWaferRecHits.cc.

{}

endJob()

void HGCalTestPartialWaferRecHits::endJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 39 of file HGCalTestPartialWaferRecHits.cc.

{}

fillDescriptions()

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

static

Definition at line 86 of file HGCalTestPartialWaferRecHits.cc.

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

                                                                                              {
  edm::ParameterSetDescription desc;
  desc.add<std::string>("detectorName", "HGCalEESensitive");
  desc.add<std::string>("missingFile", "");
  desc.add<edm::InputTag>("source", edm::InputTag("HGCalRecHit", "HGCEERecHits"));
  descriptions.add("hgcalRecHitPartialEE", desc);
}

Member Data Documentation

dumpDets_

std::vector<int> HGCalTestPartialWaferRecHits::dumpDets_

private

Definition at line 47 of file HGCalTestPartialWaferRecHits.cc.

Referenced by analyze().

missingFile_

const std::string HGCalTestPartialWaferRecHits::missingFile_

private

Definition at line 42 of file HGCalTestPartialWaferRecHits.cc.

nameDetector_

const std::string HGCalTestPartialWaferRecHits::nameDetector_

private

Definition at line 42 of file HGCalTestPartialWaferRecHits.cc.

Referenced by HGCalTestPartialWaferRecHits().

recHitSource_

const edm::EDGetTokenT<HGCRecHitCollection> HGCalTestPartialWaferRecHits::recHitSource_

private

Definition at line 44 of file HGCalTestPartialWaferRecHits.cc.

Referenced by analyze().

source_

const edm::InputTag HGCalTestPartialWaferRecHits::source_

private

Definition at line 43 of file HGCalTestPartialWaferRecHits.cc.

Referenced by Config.Process::dumpConfig(), Config.Process::dumpPython(), and Config.Process::splitPython().

tok_hgcGeom_

const edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalTestPartialWaferRecHits::tok_hgcGeom_

private

Definition at line 45 of file HGCalTestPartialWaferRecHits.cc.

Referenced by analyze().

wafers_

std::vector<int> HGCalTestPartialWaferRecHits::wafers_

private

Definition at line 46 of file HGCalTestPartialWaferRecHits.cc.

Referenced by analyze().