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 noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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
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(), caHitNtupletGeneratorKernels::good, DetId::HGCalEE, DetId::HGCalHSi, l1ctLayer2EG_cff::id, info(), edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, HGCSiliconDetId::layer(), TrackingDataMCValidation_Standalone_cff::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(), false);
          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, ProducerED_cfi::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().