|
| HGCalTimingAnalyzer (edm::ParameterSet const &) |
|
| ~HGCalTimingAnalyzer () override=default |
|
| 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 |
|
void | callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func) |
|
| EDAnalyzerBase () |
|
ModuleDescription const & | moduleDescription () const |
|
bool | wantsStreamLuminosityBlocks () const |
|
bool | wantsStreamRuns () const |
|
| ~EDAnalyzerBase () override |
|
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) |
|
|
typedef EDAnalyzerBase | ModuleType |
|
typedef ProductLabels | Labels |
|
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) |
|
Definition at line 44 of file HGCalTimingAnalyzer.cc.
Implements edm::one::EDAnalyzerBase.
Definition at line 161 of file HGCalTimingAnalyzer.cc.
References analyzeSimHits(), analyzeSimTracks(), edm::HandleBase::clear(), detectorBeam_, detectorEE_, doTree_, edm::HepMCProduct::GetEvent(), iEvent, edm::HandleBase::isValid(), dqmdumpme::k, AlCaHLTBitMon_ParallelJobs::p, simHitCellEnBeam_, simHitCellEnEE_, simHitCellIdBeam_, simHitCellIdEE_, simHitCellTmBeam_, simHitCellTmEE_, tok_hepMC_, tok_hitsBeam_, tok_hitsEE_, tok_simTk_, tok_simVtx_, and tree_.
170 for (HepMC::GenEvent::particle_const_iterator
p = myGenEvent->particles_begin();
p != myGenEvent->particles_end();
172 edm::LogVerbatim(
"HGCSim") <<
"Particle[" <<
k <<
"] with p " << (*p)->momentum().rho() <<
" theta " 173 << (*p)->momentum().theta() <<
" phi " << (*p)->momentum().phi();
190 std::vector<PCaloHit> caloHits;
192 if (theCaloHitContainers.
isValid()) {
198 caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
207 if (caloHitContainerBeam.
isValid()) {
213 caloHits.insert(caloHits.end(), caloHitContainerBeam->begin(), caloHitContainerBeam->end());
Log< level::Info, true > LogVerbatim
const edm::EDGetTokenT< edm::PCaloHitContainer > tok_hitsBeam_
std::vector< float > simHitCellEnEE_
const edm::EDGetTokenT< edm::HepMCProduct > tok_hepMC_
std::vector< float > simHitCellTmBeam_
const std::string detectorBeam_
std::vector< float > simHitCellEnBeam_
const edm::EDGetTokenT< edm::PCaloHitContainer > tok_hitsEE_
std::vector< uint32_t > simHitCellIdBeam_
std::vector< float > simHitCellTmEE_
const HepMC::GenEvent * GetEvent() const
std::vector< uint32_t > simHitCellIdEE_
const edm::EDGetTokenT< edm::SimVertexContainer > tok_simVtx_
const edm::EDGetTokenT< edm::SimTrackContainer > tok_simTk_
const std::string detectorEE_
Log< level::Warning, false > LogWarning
void analyzeSimHits(int type, std::vector< PCaloHit > const &hits)
void analyzeSimTracks(edm::Handle< edm::SimTrackContainer > const &SimTk, edm::Handle< edm::SimVertexContainer > const &SimVtx)
void HGCalTimingAnalyzer::analyzeSimHits |
( |
int |
type, |
|
|
std::vector< PCaloHit > const & |
hits |
|
) |
| |
|
private |
Definition at line 224 of file HGCalTimingAnalyzer.cc.
References TauDecayModes::dec, hcalRecHitTable_cff::energy, dqmdumpme::first, hgcons_, hfClusterShapes_cfi::hits, mps_fire::i, hit::id, submitPVResolutionJobs::key, nano_mu_digi_cff::layer, DetId::rawId(), edm::second(), nano_mu_digi_cff::sector, simHitCellEnBeam_, simHitCellEnEE_, simHitCellIdBeam_, simHitCellIdEE_, simHitCellTmBeam_, simHitCellTmEE_, HGCalTBDDDConstants::simToReco(), nano_mu_digi_cff::subsector, hcalRecHitTable_cff::time, timeUnit_, HGCalTestNumbering::unpackHexagonIndex(), HcalTestBeamNumbering::unpackIndex(), x, y, and ecaldqm::zside().
Referenced by analyze().
228 std::map<std::pair<uint32_t, uint64_t>, std::pair<double, double>> map_hits;
232 uint32_t
id =
hit.
id();
242 << recoLayerCell.second <<
" Energy " <<
energy <<
" Time " <<
time;
248 edm::LogVerbatim(
"HGCSim") <<
"SimHit:Hit[" <<
i <<
"] Beam Subdet " << subdet <<
" Layer " <<
layer <<
" x|y " 249 <<
x <<
":" <<
y <<
" Energy " <<
energy <<
" Time " <<
time;
253 std::pair<uint32_t, uint64_t>
key(
id, tid);
254 auto itr = map_hits.find(
key);
255 if (itr == map_hits.end()) {
256 map_hits[
key] = std::pair<double, double>(
time, 0.0);
257 itr = map_hits.find(
key);
260 map_hits[
key] = std::pair<double, double>((itr->second).
first,
energy);
267 edm::LogVerbatim(
"HGCSim") <<
"analyzeSimHits: Finds " << map_hits.size() <<
" hits " 268 <<
" from the Hit Vector of size " <<
hits.size() <<
" for type " <<
type;
270 for (
const auto& itr : map_hits) {
271 uint32_t
id = (itr.first).
first;
Log< level::Info, true > LogVerbatim
const HGCalTBDDDConstants * hgcons_
U second(std::pair< T, U > const &p)
std::vector< float > simHitCellEnEE_
key
prepare the HTCondor submission files and eventually submit them
std::vector< float > simHitCellTmBeam_
std::vector< float > simHitCellEnBeam_
std::vector< uint32_t > simHitCellIdBeam_
std::vector< float > simHitCellTmEE_
std::vector< uint32_t > simHitCellIdEE_
unsigned long long uint64_t
constexpr uint32_t rawId() const
get the raw id
static void unpackIndex(const uint32_t &idx, int &det, int &lay, int &x, int &y)
static void unpackHexagonIndex(const uint32_t &idx, int &subdet, int &z, int &lay, int &wafer, int &celltyp, int &cell)
std::pair< int, int > simToReco(int cell, int layer, int mod, bool half) const