Inheritance diagram for MaterialBudgetVolumeAnalysis:

Public Member Functions
	MaterialBudgetVolumeAnalysis (edm::ParameterSet const &)

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	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)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	bookHisto ()

Private Attributes
const int	binEta_

const int	binPhi_

const double	etaHigh_

const double	etaLow_

std::vector< TProfile * >	meIntLEta_

std::vector< TProfile2D * >	meIntLEtaPhi_

std::vector< TProfile * >	meIntLPhi_

std::vector< TProfile * >	meRadLEta_

std::vector< TProfile2D * >	meRadLEtaPhi_

std::vector< TProfile * >	meRadLPhi_

std::vector< TProfile * >	meStepEta_

std::vector< TProfile2D * >	meStepEtaPhi_

std::vector< TProfile * >	meStepPhi_

const std::vector< std::string >	names_

const double	phiHigh_

const double	phiLow_

const edm::InputTag	tag_

edm::EDGetTokenT< edm::MaterialInformationContainer >	tok_info_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

Detailed Description

Definition at line 27 of file MaterialBudgetVolumeAnalysis.cc.

Constructor & Destructor Documentation

◆ MaterialBudgetVolumeAnalysis()

MaterialBudgetVolumeAnalysis::MaterialBudgetVolumeAnalysis ( edm::ParameterSet const & p )

explicit

Definition at line 48 of file MaterialBudgetVolumeAnalysis.cc.

References binEta_, binPhi_, bookHisto(), etaHigh_, etaLow_, dqmdumpme::k, TFileService::kSharedResource, names_, tag_, and tok_info_.

     : names_(p.getParameter<std::vector<std::string> >("names")),
       tag_(p.getParameter<edm::InputTag>("inputTag")),
       binEta_(p.getParameter<int>("nBinEta")),
       binPhi_(p.getParameter<int>("nBinPhi")),
       etaLow_(p.getParameter<double>("etaLow")),
       etaHigh_(p.getParameter<double>("etaHigh")),
       phiLow_(-1._pi),
       phiHigh_(1._pi) {
   usesResource(TFileService::kSharedResource);
   tok_info_ = consumes<edm::MaterialInformationContainer>(tag_);
 
   edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolumeAnalysis: Eta plot: NX " << binEta_ << " Range "
                                      << -etaLow_ << ":" << etaHigh_ << " Phi plot: NX " << binPhi_ << " Range "
                                      << -1._pi << ":" << 1._pi << " for " << names_.size() << " detectors from "
                                      << tag_;
   std::ostringstream st1;
   for (unsigned int k = 0; k < names_.size(); ++k)
     st1 << " [" << k << "] " << names_[k];
   edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume: " << st1.str();
   bookHisto();
 }

Member Function Documentation

◆ analyze()

void MaterialBudgetVolumeAnalysis::analyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 96 of file MaterialBudgetVolumeAnalysis.cc.

References spr::find(), iEvent, edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmdumpme::k, meIntLEta_, meIntLEtaPhi_, meIntLPhi_, meRadLEta_, meRadLEtaPhi_, meRadLPhi_, meStepEta_, meStepEtaPhi_, meStepPhi_, names_, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::ntot, edm::Handle< T >::product(), and tok_info_.

                                                                                        {
   edm::Handle<edm::MaterialInformationContainer> materialInformationContainer;
   iEvent.getByToken(tok_info_, materialInformationContainer);
 #ifdef EDM_ML_DEBUG
   unsigned int nsize(0), ntot(0), nused(0);
 #endif
   if (materialInformationContainer.isValid()) {
 #ifdef EDM_ML_DEBUG
     nsize = materialInformationContainer->size();
 #endif
     for (const auto& it : *(materialInformationContainer.product())) {
 #ifdef EDM_ML_DEBUG
       ntot++;
 #endif
       if (std::find(names_.begin(), names_.end(), it.vname()) != names_.end()) {
 #ifdef EDM_ML_DEBUG
         nused++;
 #endif
         unsigned int k =
             static_cast<unsigned int>(std::find(names_.begin(), names_.end(), it.vname()) - names_.begin());
         meStepEta_[k]->Fill(it.trackEta(), it.stepLength());
         meRadLEta_[k]->Fill(it.trackEta(), it.radiationLength());
         meIntLEta_[k]->Fill(it.trackEta(), it.interactionLength());
         meStepPhi_[k]->Fill(it.trackPhi(), it.stepLength());
         meRadLPhi_[k]->Fill(it.trackPhi(), it.radiationLength());
         meIntLPhi_[k]->Fill(it.trackPhi(), it.interactionLength());
         meStepEtaPhi_[k]->Fill(it.trackEta(), it.trackPhi(), it.stepLength());
         meRadLEtaPhi_[k]->Fill(it.trackEta(), it.trackPhi(), it.radiationLength());
         meIntLEtaPhi_[k]->Fill(it.trackEta(), it.trackPhi(), it.interactionLength());
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MaterialBudget") << "MaterialInformation with " << nsize << ":" << ntot << " elements of which "
                                      << nused << " are used";
 #endif
 }

◆ bookHisto()

void MaterialBudgetVolumeAnalysis::bookHisto ( )

private

Definition at line 134 of file MaterialBudgetVolumeAnalysis.cc.

References binEta_, binPhi_, etaHigh_, etaLow_, compareTotals::fs, dqmdumpme::k, meIntLEta_, meIntLEtaPhi_, meIntLPhi_, meRadLEta_, meRadLEtaPhi_, meRadLPhi_, meStepEta_, meStepEtaPhi_, meStepPhi_, Skims_PA_cff::name, names_, phiHigh_, phiLow_, and runGCPTkAlMap::title.

Referenced by MaterialBudgetVolumeAnalysis().

                                              {
   edm::Service<TFileService> fs;
   char name[40], title[100];
   for (unsigned int k = 0; k < names_.size(); ++k) {
     sprintf(name, "stepEta%s", names_[k].c_str());
     sprintf(title, "MB(Step) vs #eta for %s", names_[k].c_str());
     meStepEta_.emplace_back(fs->make<TProfile>(name, title, binEta_, etaLow_, etaHigh_));
     sprintf(name, "radlEta%s", names_[k].c_str());
     sprintf(title, "MB(X0) vs #eta for %s", names_[k].c_str());
     meRadLEta_.emplace_back(fs->make<TProfile>(name, title, binEta_, etaLow_, etaHigh_));
     sprintf(name, "intlEta%s", names_[k].c_str());
     sprintf(title, "MB(L0) vs #eta for %s", names_[k].c_str());
     meIntLEta_.emplace_back(fs->make<TProfile>(name, title, binEta_, etaLow_, etaHigh_));
     sprintf(name, "stepPhi%s", names_[k].c_str());
     sprintf(title, "MB(Step) vs #phi for %s", names_[k].c_str());
     meStepPhi_.emplace_back(fs->make<TProfile>(name, title, binPhi_, phiLow_, phiHigh_));
     sprintf(name, "radlPhi%s", names_[k].c_str());
     sprintf(title, "MB(X0) vs #phi for %s", names_[k].c_str());
     meRadLPhi_.emplace_back(fs->make<TProfile>(name, title, binPhi_, phiLow_, phiHigh_));
     sprintf(name, "intlPhi%s", names_[k].c_str());
     sprintf(title, "MB(L0) vs #phi for %s", names_[k].c_str());
     meIntLPhi_.emplace_back(fs->make<TProfile>(name, title, binPhi_, phiLow_, phiHigh_));
     sprintf(name, "stepEtaPhi%s", names_[k].c_str());
     sprintf(title, "MB(Step) vs #eta and #phi for %s", names_[k].c_str());
     meStepEtaPhi_.emplace_back(
         fs->make<TProfile2D>(name, title, binEta_ / 2, etaLow_, etaHigh_, binPhi_ / 2, phiLow_, phiHigh_));
     sprintf(name, "radlEtaPhi%s", names_[k].c_str());
     sprintf(title, "MB(X0) vs #eta and #phi for %s", names_[k].c_str());
     meRadLEtaPhi_.emplace_back(
         fs->make<TProfile2D>(name, title, binEta_ / 2, etaLow_, etaHigh_, binPhi_ / 2, phiLow_, phiHigh_));
     sprintf(name, "intlEtaPhi%s", names_[k].c_str());
     sprintf(title, "MB(L0) vs #eta and #phi for %s", names_[k].c_str());
     meIntLEtaPhi_.emplace_back(
         fs->make<TProfile2D>(name, title, binEta_ / 2, etaLow_, etaHigh_, binPhi_ / 2, phiLow_, phiHigh_));
   }
 }

◆ fillDescriptions()

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

static

Definition at line 71 of file MaterialBudgetVolumeAnalysis.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and names.

                                                                                               {
   edm::ParameterSetDescription desc;
   std::vector<std::string> names = {"BEAM",
                                     "BEAM1",
                                     "BEAM2",
                                     "BEAM3",
                                     "BEAM4",
                                     "Tracker",
                                     "ECAL",
                                     "HCal",
                                     "MUON",
                                     "VCAL",
                                     "MGNT",
                                     "OQUA",
                                     "CALOEC",
                                     "HFNoseVol"};
   desc.add<std::vector<std::string> >("names", names);
   desc.add<edm::InputTag>("inputTag", edm::InputTag("g4SimHits", "MaterialInformation"));
   desc.add<int>("nBinEta", 300);
   desc.add<int>("nBinPhi", 180);
   desc.add<double>("etaLow", -6.0);
   desc.add<double>("etaHigh", 6.0);
   descriptions.add("materialBudgetVolumeAnalysis", desc);
 }