HGCalTestScintHits Class Reference

Inheritance diagram for HGCalTestScintHits:

Public Member Functions
	HGCalTestScintHits (const edm::ParameterSet &ps)
	~HGCalTestScintHits () 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_
const std::string	g4Label_
const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geomToken_
const std::string	nameSense_
const std::string	tileFileName_
std::vector< int >	tiles_
const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_calo_

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 HGCalTestScintHits.cc.

Constructor & Destructor Documentation

HGCalTestScintHits()

HGCalTestScintHits::HGCalTestScintHits

(

const edm::ParameterSet &

ps

)

Definition at line 48 of file HGCalTestScintHits.cc.

References nameSense_.

    : g4Label_(ps.getParameter<std::string>("moduleLabel")),
      caloHitSource_(ps.getParameter<std::string>("caloHitSource")),
      nameSense_(ps.getParameter<std::string>("nameSense")),
      tileFileName_(ps.getParameter<std::string>("tileFileName")),
      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_ << " Tile file " << tileFileName_;
  if (!tileFileName_.empty()) {
    edm::FileInPath filetmp("SimG4CMS/Calo/data/" + tileFileName_);
    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 ring = std::atoi(items[1].c_str());
          int phi = std::atoi(items[2].c_str());
          tiles_.emplace_back(HGCalTileIndex::tileIndex(layer, ring, phi));
        }
      }
      edm::LogVerbatim("HGCalSim") << "Reads in " << tiles_.size() << " tile information from " << fileName;
      fInput.close();
    }
  }
}

~HGCalTestScintHits()

HGCalTestScintHits::~HGCalTestScintHits ( )

overridedefault

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 89 of file HGCalTestScintHits.cc.

References python.cmstools::all(), MillePedeFileConverter_cfg::e, spr::find(), relativeConstraints::geom, geomToken_, edm::EventSetup::getData(), caHitNtupletGeneratorKernels::good, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, hit::id, edm::HandleBase::isValid(), nhits, HGCalTileIndex::tileIndex(), tiles_, and tok_calo_.

                                                                           {
  // get hcalGeometry
  const HGCalGeometry* geom = &iS.getData(geomToken_);
  const HGCalDDDConstants& hgc = geom->topology().dddConstants();
  int firstLayer = hgc.firstLayer() - 1;
  // 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), all(0);
  edm::LogVerbatim("HGCalSim") << "HGCalTestScintHits: Input flags Hits " << getHits << " with " << nhits
                               << " hits first Layer " << 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());
        HGCScintillatorDetId hid(id);
        std::pair<int, int> typm = hgc.tileType(hid.layer(), hid.ring(), 0);
        if (typm.first >= 0) {
          hid.setType(typm.first);
          hid.setSiPM(typm.second);
          id = static_cast<DetId>(hid);
        }
        bool toCheck(true);
        if (!tiles_.empty()) {
          int indx = HGCalTileIndex::tileIndex(firstLayer + hid.layer(), hid.ring(), hid.iphi());
          if (std::find(tiles_.begin(), tiles_.end(), indx) != tiles_.end())
            toCheck = true;
        }
        if (toCheck) {
          ++good;
          GlobalPoint pos = geom->getPosition(id);
          bool valid1 = geom->topology().valid(id);
          bool valid2 = hgc.isValidTrap(hid.zside(), hid.layer(), hid.ring(), hid.iphi());
          edm::LogVerbatim("HGCalSim") << "Hit[" << all << ":" << good << "]" << hid << " at (" << pos.x() << ", "
                                       << pos.y() << ", " << pos.z() << ") Validity " << valid1 << ":" << valid2
                                       << " Energy " << hit.energy() << " Time " << hit.time();
        }
      }
    }
  }
  edm::LogVerbatim("HGCalSim") << "Total hits = " << all << ":" << nhits << " Good DetIds = " << good;
}

beginJob()

void HGCalTestScintHits::beginJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 38 of file HGCalTestScintHits.cc.

{}

endJob()

void HGCalTestScintHits::endJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 39 of file HGCalTestScintHits.cc.

{}

fillDescriptions()

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

static

Definition at line 80 of file HGCalTestScintHits.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", "HGCHitsHEback");
  desc.add<std::string>("nameSense", "HGCalHEScintillatorSensitive");
  desc.add<std::string>("tileFileName", "");
  descriptions.add("hgcalHitScintillator", desc);
}

Member Data Documentation

caloHitSource_

const std::string HGCalTestScintHits::caloHitSource_

private

Definition at line 42 of file HGCalTestScintHits.cc.

g4Label_

const std::string HGCalTestScintHits::g4Label_

private

Definition at line 42 of file HGCalTestScintHits.cc.

geomToken_

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

private

Definition at line 44 of file HGCalTestScintHits.cc.

Referenced by analyze().

nameSense_

const std::string HGCalTestScintHits::nameSense_

private

Definition at line 42 of file HGCalTestScintHits.cc.

Referenced by HGCalTestScintHits().

tileFileName_

const std::string HGCalTestScintHits::tileFileName_

private

Definition at line 42 of file HGCalTestScintHits.cc.

tiles_

std::vector<int> HGCalTestScintHits::tiles_

private

Definition at line 45 of file HGCalTestScintHits.cc.

Referenced by analyze().

tok_calo_

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

private

Definition at line 43 of file HGCalTestScintHits.cc.

Referenced by analyze().