HGCalValidator Class Reference

#include <HGCalValidator.h>

Inheritance diagram for HGCalValidator:

Public Types
using	Histograms = HGCalValidatorHistograms
Public Types inherited from DQMGlobalEDAnalyzerBase< HGCalValidatorHistograms, Args...>
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &, Histograms &) const override
	Method called to book the DQM histograms. More...
void	cpParametersAndSelection (const Histograms &histograms, std::vector< CaloParticle > const &cPeff, std::vector< SimVertex > const &simVertices, std::vector< size_t > &selected_cPeff, unsigned int layers, std::unordered_map< DetId, const HGCRecHit * > const &) const
void	dqmAnalyze (const edm::Event &, const edm::EventSetup &, const Histograms &) const override
	Method called once per event. More...
	HGCalValidator (const edm::ParameterSet &pset)
	Constructor. More...
	~HGCalValidator () override
	Destructor. More...
Public Member Functions inherited from DQMGlobalEDAnalyzer< HGCalValidatorHistograms >
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &, HGCalValidatorHistogramsconst &) const
void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup) const final
Public Member Functions inherited from DQMGlobalEDAnalyzerBase< HGCalValidatorHistograms, Args...>
void	accumulate (edm::StreamID id, edm::Event const &event, edm::EventSetup const &setup) const final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &, HGCalValidatorHistograms &) const
	DQMGlobalEDAnalyzerBase ()
std::shared_ptr < HGCalValidatorHistograms >	globalBeginRun (edm::Run const &run, edm::EventSetup const &setup) const final
void	globalEndRun (edm::Run const &, edm::EventSetup const &) const final
Public Member Functions inherited from edm::global::EDProducer< T >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
EDProducer &	operator= (const EDProducer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
bool	wantsStreamLuminosityBlocks () const final
bool	wantsStreamRuns () const final
Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) 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)

Protected Attributes
edm::InputTag	associator_
edm::EDGetTokenT < hgcal::RecoToSimCollection >	associatorMapRtS
edm::EDGetTokenT < hgcal::RecoToSimCollectionWithSimClusters >	associatorMapRtSim
edm::EDGetTokenT < hgcal::SimToRecoCollectionWithSimClusters >	associatorMapSimtR
edm::EDGetTokenT < hgcal::SimToRecoCollection >	associatorMapStR
edm::InputTag	associatorSim_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_
std::vector< edm::EDGetTokenT < std::vector< float > > >	clustersMaskTokens_
const edm::FileInPath	cummatbudinxo_
edm::EDGetTokenT< Density >	density_
const bool	doCaloParticlePlots_
const bool	doCaloParticleSelection_
const bool	doLayerClustersPlots_
const bool	doSimClustersPlots_
const bool	doTrackstersPlots_
std::unique_ptr < HGVHistoProducerAlgo >	histoProducerAlgo_
edm::EDGetTokenT < std::unordered_map< DetId, const HGCRecHit * > >	hitMap_
std::vector< edm::InputTag >	label_clustersmask
edm::EDGetTokenT< std::vector < CaloParticle > >	label_cp_effic
edm::EDGetTokenT< std::vector < CaloParticle > >	label_cp_fake
edm::InputTag	label_layerClustersPlots_
edm::InputTag	label_lcl
edm::InputTag	label_LCToCPLinking_
edm::InputTag	label_SimClustersLevel_
edm::InputTag	label_SimClustersPlots_
edm::InputTag	label_simTSFromCP
std::vector< edm::InputTag >	label_tst
edm::InputTag	label_TSToCPLinking_
std::vector< edm::EDGetTokenT < ticl::TracksterCollection > >	label_tstTokens
std::vector< edm::EDGetTokenT < reco::CaloClusterCollection > >	labelToken
edm::EDGetTokenT < reco::CaloClusterCollection >	layerclusters_
const bool	SaveGeneralInfo_
edm::EDGetTokenT< std::vector < SimCluster > >	simClusters_
edm::EDGetTokenT < ticl::TracksterCollection >	simTrackstersFromCPs_
edm::EDGetTokenT< std::vector < SimVertex > >	simVertices_
Protected Attributes inherited from DQMGlobalEDAnalyzerBase< HGCalValidatorHistograms, Args...>
DQMStore *	dqmstore_
edm::EDPutTokenT< DQMToken >	runToken_

Private Attributes
CaloParticleSelector	cpSelector
std::map< double, double >	cummatbudg
std::string	dirName_
std::vector< int >	particles_to_monitor_
std::vector< int >	thicknesses_to_monitor_
std::shared_ptr < hgcal::RecHitTools >	tools_
unsigned	totallayers_to_monitor_

Additional Inherited Members
Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from DQMGlobalEDAnalyzerBase< HGCalValidatorHistograms, Args...>
uint64_t	meId (edm::Run const &run) const
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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 () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
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

Class that produces histograms to validate HGCal Reconstruction performances

Author

HGCal

Definition at line 39 of file HGCalValidator.h.

Member Typedef Documentation

using HGCalValidator::Histograms = HGCalValidatorHistograms

Definition at line 41 of file HGCalValidator.h.

Constructor & Destructor Documentation

HGCalValidator::HGCalValidator

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 11 of file HGCalValidator.cc.

References associator_, associatorMapRtS, associatorMapRtSim, associatorMapSimtR, associatorMapStR, associatorSim_, clustersMaskTokens_, edm::EDConsumerBase::consumes(), cpSelector, cummatbudg, cummatbudinxo_, density_, dirName_, fmb, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), histoProducerAlgo_, hitMap_, HLT_FULL_cff::InputTag, label_clustersmask, label_cp_effic, label_cp_fake, label_lcl, label_simTSFromCP, label_tst, label_tstTokens, layerclusters_, particles_to_monitor_, simClusters_, simTrackstersFromCPs_, simVertices_, AlCaHLTBitMon_QueryRunRegistry::string, thicknesses_to_monitor_, tools_, and totallayers_to_monitor_.

    : caloGeomToken_(esConsumes<CaloGeometry, CaloGeometryRecord>()),
      label_lcl(pset.getParameter<edm::InputTag>("label_lcl")),
      label_tst(pset.getParameter<std::vector<edm::InputTag>>("label_tst")),
      label_simTSFromCP(pset.getParameter<edm::InputTag>("label_simTSFromCP")),
      associator_(pset.getUntrackedParameter<edm::InputTag>("associator")),
      associatorSim_(pset.getUntrackedParameter<edm::InputTag>("associatorSim")),
      SaveGeneralInfo_(pset.getUntrackedParameter<bool>("SaveGeneralInfo")),
      doCaloParticlePlots_(pset.getUntrackedParameter<bool>("doCaloParticlePlots")),
      doCaloParticleSelection_(pset.getUntrackedParameter<bool>("doCaloParticleSelection")),
      doSimClustersPlots_(pset.getUntrackedParameter<bool>("doSimClustersPlots")),
      label_SimClustersPlots_(pset.getParameter<edm::InputTag>("label_SimClusters")),
      label_SimClustersLevel_(pset.getParameter<edm::InputTag>("label_SimClustersLevel")),
      doLayerClustersPlots_(pset.getUntrackedParameter<bool>("doLayerClustersPlots")),
      label_layerClustersPlots_(pset.getParameter<edm::InputTag>("label_layerClusterPlots")),
      label_LCToCPLinking_(pset.getParameter<edm::InputTag>("label_LCToCPLinking")),
      doTrackstersPlots_(pset.getUntrackedParameter<bool>("doTrackstersPlots")),
      label_TSToCPLinking_(pset.getParameter<edm::InputTag>("label_TSToCPLinking")),
      label_clustersmask(pset.getParameter<std::vector<edm::InputTag>>("LayerClustersInputMask")),
      cummatbudinxo_(pset.getParameter<edm::FileInPath>("cummatbudinxo")) {
  //In this way we can easily generalize to associations between other objects also.
  const edm::InputTag& label_cp_effic_tag = pset.getParameter<edm::InputTag>("label_cp_effic");
  const edm::InputTag& label_cp_fake_tag = pset.getParameter<edm::InputTag>("label_cp_fake");

  label_cp_effic = consumes<std::vector<CaloParticle>>(label_cp_effic_tag);
  label_cp_fake = consumes<std::vector<CaloParticle>>(label_cp_fake_tag);

  simVertices_ = consumes<std::vector<SimVertex>>(pset.getParameter<edm::InputTag>("simVertices"));

  for (auto& itag : label_clustersmask) {
    clustersMaskTokens_.push_back(consumes<std::vector<float>>(itag));
  }

  associatorMapSimtR = consumes<hgcal::SimToRecoCollectionWithSimClusters>(associatorSim_);
  associatorMapRtSim = consumes<hgcal::RecoToSimCollectionWithSimClusters>(associatorSim_);

  hitMap_ = consumes<std::unordered_map<DetId, const HGCRecHit*>>(edm::InputTag("hgcalRecHitMapProducer"));

  density_ = consumes<Density>(edm::InputTag("hgcalLayerClusters"));

  simClusters_ = consumes<std::vector<SimCluster>>(pset.getParameter<edm::InputTag>("label_scl"));

  layerclusters_ = consumes<reco::CaloClusterCollection>(label_lcl);

  for (auto& itag : label_tst) {
    label_tstTokens.push_back(consumes<ticl::TracksterCollection>(itag));
  }

  simTrackstersFromCPs_ = consumes<ticl::TracksterCollection>(label_simTSFromCP);

  associatorMapRtS = consumes<hgcal::RecoToSimCollection>(associator_);
  associatorMapStR = consumes<hgcal::SimToRecoCollection>(associator_);

  cpSelector = CaloParticleSelector(pset.getParameter<double>("ptMinCP"),
                                    pset.getParameter<double>("ptMaxCP"),
                                    pset.getParameter<double>("minRapidityCP"),
                                    pset.getParameter<double>("maxRapidityCP"),
                                    pset.getParameter<double>("lipCP"),
                                    pset.getParameter<double>("tipCP"),
                                    pset.getParameter<int>("minHitCP"),
                                    pset.getParameter<int>("maxSimClustersCP"),
                                    pset.getParameter<bool>("signalOnlyCP"),
                                    pset.getParameter<bool>("intimeOnlyCP"),
                                    pset.getParameter<bool>("chargedOnlyCP"),
                                    pset.getParameter<bool>("stableOnlyCP"),
                                    pset.getParameter<bool>("notConvertedOnlyCP"),
                                    pset.getParameter<std::vector<int>>("pdgIdCP"));

  tools_.reset(new hgcal::RecHitTools());

  particles_to_monitor_ = pset.getParameter<std::vector<int>>("pdgIdCP");
  totallayers_to_monitor_ = pset.getParameter<int>("totallayers_to_monitor");
  thicknesses_to_monitor_ = pset.getParameter<std::vector<int>>("thicknesses_to_monitor");

  //For the material budget file here
  std::ifstream fmb(cummatbudinxo_.fullPath().c_str());
  double thelay = 0.;
  double mbg = 0.;
  for (unsigned ilayer = 1; ilayer <= totallayers_to_monitor_; ++ilayer) {
    fmb >> thelay >> mbg;
    cummatbudg.insert(std::pair<double, double>(thelay, mbg));
  }

  fmb.close();

  ParameterSet psetForHistoProducerAlgo = pset.getParameter<ParameterSet>("histoProducerAlgoBlock");
  histoProducerAlgo_ = std::make_unique<HGVHistoProducerAlgo>(psetForHistoProducerAlgo);

  dirName_ = pset.getParameter<std::string>("dirName");
}

HGCalValidator::~HGCalValidator ( )

override

Destructor.

Definition at line 102 of file HGCalValidator.cc.

{}

Member Function Documentation

void HGCalValidator::bookHistograms ( DQMStore::IBooker &  ibook, edm::Run const &  , edm::EventSetup const &  setup, Histograms &  histograms  ) const

overridevirtual

Method called to book the DQM histograms.

Implements DQMGlobalEDAnalyzerBase< HGCalValidatorHistograms, Args...>.

Definition at line 104 of file HGCalValidator.cc.

References dqm::implementation::NavigatorBase::cd(), cummatbudinxo_, TrackerOfflineValidation_Dqm_cff::dirName, dirName_, doCaloParticlePlots_, doLayerClustersPlots_, doSimClustersPlots_, doTrackstersPlots_, edm::FileInPath::fullPath(), HGCalValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, edm::InputTag::instance(), edm::InputTag::label(), label_clustersmask, label_layerClustersPlots_, label_LCToCPLinking_, label_SimClustersLevel_, label_SimClustersPlots_, label_tst, label_TSToCPLinking_, LogDebug, particles_to_monitor_, edm::InputTag::process(), SaveGeneralInfo_, dqm::implementation::NavigatorBase::setCurrentFolder(), thicknesses_to_monitor_, totallayers_to_monitor_, and gpuVertexFinder::ws.

                                                                  {
  if (SaveGeneralInfo_) {
    ibook.cd();
    ibook.setCurrentFolder(dirName_ + "GeneralInfo");
    histoProducerAlgo_->bookInfo(ibook, histograms.histoProducerAlgo);
  }

  if (doCaloParticlePlots_) {
    ibook.cd();

    for (auto const particle : particles_to_monitor_) {
      ibook.setCurrentFolder(dirName_ + "SelectedCaloParticles/" + std::to_string(particle));
      histoProducerAlgo_->bookCaloParticleHistos(
          ibook, histograms.histoProducerAlgo, particle, totallayers_to_monitor_);
    }
    ibook.cd();
    ibook.setCurrentFolder(dirName_);
  }

  //Booking histograms concerning with simClusters
  if (doSimClustersPlots_) {
    ibook.cd();
    ibook.setCurrentFolder(dirName_ + label_SimClustersPlots_.label() + "/" + label_SimClustersLevel_.label());
    histoProducerAlgo_->bookSimClusterHistos(
        ibook, histograms.histoProducerAlgo, totallayers_to_monitor_, thicknesses_to_monitor_);

    for (unsigned int ws = 0; ws < label_clustersmask.size(); ws++) {
      ibook.cd();
      InputTag algo = label_clustersmask[ws];
      string dirName = dirName_ + label_SimClustersPlots_.label() + "/";
      if (!algo.process().empty())
        dirName += algo.process() + "_";
      LogDebug("HGCalValidator") << dirName << "\n";
      if (!algo.label().empty())
        dirName += algo.label() + "_";
      LogDebug("HGCalValidator") << dirName << "\n";
      if (!algo.instance().empty())
        dirName += algo.instance() + "_";
      LogDebug("HGCalValidator") << dirName << "\n";

      if (!dirName.empty()) {
        dirName.resize(dirName.size() - 1);
      }

      LogDebug("HGCalValidator") << dirName << "\n";

      ibook.setCurrentFolder(dirName);

      histoProducerAlgo_->bookSimClusterAssociationHistos(
          ibook, histograms.histoProducerAlgo, totallayers_to_monitor_, thicknesses_to_monitor_);
    }  //end of loop over masks
  }    //if for simCluster plots

  //Booking histograms concerning with hgcal layer clusters
  if (doLayerClustersPlots_) {
    ibook.cd();
    ibook.setCurrentFolder(dirName_ + label_layerClustersPlots_.label() + "/ClusterLevel");
    histoProducerAlgo_->bookClusterHistos_ClusterLevel(ibook,
                                                       histograms.histoProducerAlgo,
                                                       totallayers_to_monitor_,
                                                       thicknesses_to_monitor_,
                                                       cummatbudinxo_.fullPath());
    ibook.cd();
    ibook.setCurrentFolder(dirName_ + label_layerClustersPlots_.label() + "/" + label_LCToCPLinking_.label());
    histoProducerAlgo_->bookClusterHistos_LCtoCP_association(
        ibook, histograms.histoProducerAlgo, totallayers_to_monitor_);

    ibook.cd();
    ibook.setCurrentFolder(dirName_ + label_layerClustersPlots_.label() + "/CellLevel");
    histoProducerAlgo_->bookClusterHistos_CellLevel(
        ibook, histograms.histoProducerAlgo, totallayers_to_monitor_, thicknesses_to_monitor_);
  }

  //Booking histograms for Tracksters
  for (unsigned int www = 0; www < label_tst.size(); www++) {
    ibook.cd();
    InputTag algo = label_tst[www];
    string dirName = dirName_;
    if (!algo.process().empty())
      dirName += algo.process() + "_";
    LogDebug("HGCalValidator") << dirName << "\n";
    if (!algo.label().empty())
      dirName += algo.label() + "_";
    LogDebug("HGCalValidator") << dirName << "\n";
    if (!algo.instance().empty())
      dirName += algo.instance() + "_";
    LogDebug("HGCalValidator") << dirName << "\n";

    if (!dirName.empty()) {
      dirName.resize(dirName.size() - 1);
    }

    LogDebug("HGCalValidator") << dirName << "\n";

    ibook.setCurrentFolder(dirName);

    //Booking histograms concerning HGCal tracksters
    if (doTrackstersPlots_) {
      histoProducerAlgo_->bookTracksterHistos(ibook, histograms.histoProducerAlgo, totallayers_to_monitor_);
      ibook.setCurrentFolder(dirName + "/" + label_TSToCPLinking_.label());
      histoProducerAlgo_->bookTracksterCPLinkingHistos(ibook, histograms.histoProducerAlgo);
    }
  }  //end of booking Tracksters loop
}

void HGCalValidator::cpParametersAndSelection	(	const Histograms &	histograms,
		std::vector< CaloParticle > const &	cPeff,
		std::vector< SimVertex > const &	simVertices,
		std::vector< size_t > &	selected_cPeff,
		unsigned int	layers,
		std::unordered_map< DetId, const HGCRecHit * > const &	hitMap
	)		const

Definition at line 212 of file HGCalValidator.cc.

References cpSelector, doCaloParticlePlots_, doCaloParticleSelection_, HGCalValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, and dqmiolumiharvest::j.

Referenced by dqmAnalyze().

                                                                                                             {
  selected_cPeff.reserve(cPeff.size());

  size_t j = 0;
  for (auto const& caloParticle : cPeff) {
    int id = caloParticle.pdgId();

    if (!doCaloParticleSelection_ || (doCaloParticleSelection_ && cpSelector(caloParticle, simVertices))) {
      selected_cPeff.push_back(j);
      if (doCaloParticlePlots_) {
        histoProducerAlgo_->fill_caloparticle_histos(
            histograms.histoProducerAlgo, id, caloParticle, simVertices, layers, hitMap);
      }
    }
    ++j;
  }  //end of loop over caloparticles
}

void HGCalValidator::dqmAnalyze ( const edm::Event &  event, const edm::EventSetup &  setup, const Histograms &  histograms  ) const

overridevirtual

Method called once per event.

Implements DQMGlobalEDAnalyzerBase< HGCalValidatorHistograms, Args...>.

Definition at line 235 of file HGCalValidator.cc.

References associatorMapRtS, associatorMapRtSim, associatorMapSimtR, associatorMapStR, caloGeomToken_, caloTruthProducer_cfi::caloParticles, HLT_FULL_cff::clusters, clustersMaskTokens_, cpParametersAndSelection(), cummatbudg, density_, doLayerClustersPlots_, doSimClustersPlots_, doTrackstersPlots_, edmPickEvents::event, relativeConstraints::geom, edm::EventSetup::getHandle(), HGCalValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, hitMap_, edm::InputTag::instance(), edm::InputTag::label(), label_clustersmask, label_cp_effic, label_lcl, label_tst, label_tstTokens, layerclusters_, LogDebug, LogTrace, or, edm::InputTag::process(), dt_dqm_sourceclient_common_cff::reco, removeCPFromPU(), SaveGeneralInfo_, simClusters_, simTrackstersFromCPs_, simVertices_, thicknesses_to_monitor_, tools_, totallayers_to_monitor_, w, and gpuVertexFinder::ws.

                                                                    {
  using namespace reco;

  LogDebug("HGCalValidator") << "\n===================================================="
                             << "\n"
                             << "Analyzing new event"
                             << "\n"
                             << "====================================================\n"
                             << "\n";

  edm::Handle<std::vector<SimVertex>> simVerticesHandle;
  event.getByToken(simVertices_, simVerticesHandle);
  std::vector<SimVertex> const& simVertices = *simVerticesHandle;

  edm::Handle<std::vector<CaloParticle>> caloParticleHandle;
  event.getByToken(label_cp_effic, caloParticleHandle);
  std::vector<CaloParticle> const& caloParticles = *caloParticleHandle;

  edm::Handle<ticl::TracksterCollection> simTracksterFromCPHandle;
  event.getByToken(simTrackstersFromCPs_, simTracksterFromCPHandle);
  ticl::TracksterCollection const& simTrackstersFromCPs = *simTracksterFromCPHandle;

  edm::ESHandle<CaloGeometry> geom = setup.getHandle(caloGeomToken_);
  tools_->setGeometry(*geom);
  histoProducerAlgo_->setRecHitTools(tools_);

  edm::Handle<hgcal::SimToRecoCollection> simtorecoCollectionH;
  event.getByToken(associatorMapStR, simtorecoCollectionH);
  auto simRecColl = *simtorecoCollectionH;
  edm::Handle<hgcal::RecoToSimCollection> recotosimCollectionH;
  event.getByToken(associatorMapRtS, recotosimCollectionH);
  auto recSimColl = *recotosimCollectionH;

  edm::Handle<std::unordered_map<DetId, const HGCRecHit*>> hitMapHandle;
  event.getByToken(hitMap_, hitMapHandle);
  const std::unordered_map<DetId, const HGCRecHit*>* hitMap = &*hitMapHandle;

  //Some general info on layers etc.
  if (SaveGeneralInfo_) {
    histoProducerAlgo_->fill_info_histos(histograms.histoProducerAlgo, totallayers_to_monitor_);
  }

  std::vector<size_t> cPIndices;
  //Consider CaloParticles coming from the hard scatterer
  //excluding the PU contribution and save the indices.
  removeCPFromPU(caloParticles, cPIndices);

  // ##############################################
  // fill caloparticles histograms
  // ##############################################
  // HGCRecHit are given to select the SimHits which are also reconstructed
  LogTrace("HGCalValidator") << "\n# of CaloParticles: " << caloParticles.size() << "\n" << std::endl;
  std::vector<size_t> selected_cPeff;
  cpParametersAndSelection(histograms, caloParticles, simVertices, selected_cPeff, totallayers_to_monitor_, *hitMap);

  //get collections from the event
  //simClusters
  edm::Handle<std::vector<SimCluster>> simClustersHandle;
  event.getByToken(simClusters_, simClustersHandle);
  std::vector<SimCluster> const& simClusters = *simClustersHandle;

  //Layer clusters
  edm::Handle<reco::CaloClusterCollection> clusterHandle;
  event.getByToken(layerclusters_, clusterHandle);
  const reco::CaloClusterCollection& clusters = *clusterHandle;

  //Density
  edm::Handle<Density> densityHandle;
  event.getByToken(density_, densityHandle);
  const Density& densities = *densityHandle;

  auto nSimClusters = simClusters.size();
  std::vector<size_t> sCIndices;
  //There shouldn't be any SimTracks from different crossings, but maybe they will be added later.
  //At the moment there should be one SimTrack in each SimCluster.
  for (unsigned int scId = 0; scId < nSimClusters; ++scId) {
    if (simClusters[scId].g4Tracks()[0].eventId().event() != 0 or
        simClusters[scId].g4Tracks()[0].eventId().bunchCrossing() != 0) {
      LogDebug("HGCalValidator") << "Excluding SimClusters from event: "
                                 << simClusters[scId].g4Tracks()[0].eventId().event()
                                 << " with BX: " << simClusters[scId].g4Tracks()[0].eventId().bunchCrossing()
                                 << std::endl;
      continue;
    }
    sCIndices.emplace_back(scId);
  }

  // ##############################################
  // fill simCluster histograms
  // ##############################################
  if (doSimClustersPlots_) {
    histoProducerAlgo_->fill_simCluster_histos(
        histograms.histoProducerAlgo, simClusters, totallayers_to_monitor_, thicknesses_to_monitor_);

    for (unsigned int ws = 0; ws < label_clustersmask.size(); ws++) {
      const auto& inputClusterMask = event.get(clustersMaskTokens_[ws]);

      edm::Handle<hgcal::SimToRecoCollectionWithSimClusters> simtorecoCollectionH;
      event.getByToken(associatorMapSimtR, simtorecoCollectionH);
      auto simRecColl = *simtorecoCollectionH;
      edm::Handle<hgcal::RecoToSimCollectionWithSimClusters> recotosimCollectionH;
      event.getByToken(associatorMapRtSim, recotosimCollectionH);
      auto recSimColl = *recotosimCollectionH;

      histoProducerAlgo_->fill_simClusterAssociation_histos(histograms.histoProducerAlgo,
                                                            ws,
                                                            clusterHandle,
                                                            clusters,
                                                            simClustersHandle,
                                                            simClusters,
                                                            sCIndices,
                                                            inputClusterMask,
                                                            *hitMap,
                                                            totallayers_to_monitor_,
                                                            recSimColl,
                                                            simRecColl);

      //General Info on simClusters
      LogTrace("HGCalValidator") << "\n# of SimClusters: " << nSimClusters
                                 << ", layerClusters mask label: " << label_clustersmask[ws].label() << "\n";
    }  //end of loop overs masks
  }

  // ##############################################
  // fill layercluster histograms
  // ##############################################
  int w = 0;  //counter counting the number of sets of histograms
  if (doLayerClustersPlots_) {
    histoProducerAlgo_->fill_generic_cluster_histos(histograms.histoProducerAlgo,
                                                    w,
                                                    clusterHandle,
                                                    clusters,
                                                    densities,
                                                    caloParticleHandle,
                                                    caloParticles,
                                                    cPIndices,
                                                    selected_cPeff,
                                                    *hitMap,
                                                    cummatbudg,
                                                    totallayers_to_monitor_,
                                                    thicknesses_to_monitor_,
                                                    recSimColl,
                                                    simRecColl);

    for (unsigned int layerclusterIndex = 0; layerclusterIndex < clusters.size(); layerclusterIndex++) {
      histoProducerAlgo_->fill_cluster_histos(histograms.histoProducerAlgo, w, clusters[layerclusterIndex]);
    }

    //General Info on hgcalLayerClusters
    LogTrace("HGCalValidator") << "\n# of layer clusters with " << label_lcl.process() << ":" << label_lcl.label()
                               << ":" << label_lcl.instance() << ": " << clusters.size() << "\n";
  }

  // ##############################################
  // fill Trackster histograms
  // ##############################################
  for (unsigned int wml = 0; wml < label_tstTokens.size(); wml++) {
    if (doTrackstersPlots_) {
      edm::Handle<ticl::TracksterCollection> tracksterHandle;
      event.getByToken(label_tstTokens[wml], tracksterHandle);
      const ticl::TracksterCollection& tracksters = *tracksterHandle;

      histoProducerAlgo_->fill_trackster_histos(histograms.histoProducerAlgo,
                                                wml,
                                                tracksters,
                                                clusters,
                                                simTrackstersFromCPs,
                                                caloParticles,
                                                cPIndices,
                                                selected_cPeff,
                                                *hitMap,
                                                totallayers_to_monitor_);

      //General Info on Tracksters
      LogTrace("HGCalValidator") << "\n# of Tracksters with " << label_tst[wml].process() << ":"
                                 << label_tst[wml].label() << ":" << label_tst[wml].instance() << ": "
                                 << tracksters.size() << "\n"
                                 << std::endl;
    }
  }  //end of loop over Trackster input labels
}

Member Data Documentation

edm::InputTag HGCalValidator::associator_

protected

Definition at line 66 of file HGCalValidator.h.

Referenced by HGCalValidator().

edm::EDGetTokenT<hgcal::RecoToSimCollection> HGCalValidator::associatorMapRtS

protected

Definition at line 91 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<hgcal::RecoToSimCollectionWithSimClusters> HGCalValidator::associatorMapRtSim

protected

Definition at line 94 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<hgcal::SimToRecoCollectionWithSimClusters> HGCalValidator::associatorMapSimtR

protected

Definition at line 93 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<hgcal::SimToRecoCollection> HGCalValidator::associatorMapStR

protected

Definition at line 92 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::InputTag HGCalValidator::associatorSim_

protected

Definition at line 67 of file HGCalValidator.h.

Referenced by HGCalValidator().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HGCalValidator::caloGeomToken_

protected

Definition at line 62 of file HGCalValidator.h.

Referenced by dqmAnalyze().

std::vector<edm::EDGetTokenT<std::vector<float> > > HGCalValidator::clustersMaskTokens_

protected

Definition at line 88 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

CaloParticleSelector HGCalValidator::cpSelector

private

Definition at line 98 of file HGCalValidator.h.

Referenced by cpParametersAndSelection(), and HGCalValidator().

std::map<double, double> HGCalValidator::cummatbudg

private

Definition at line 100 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

const edm::FileInPath HGCalValidator::cummatbudinxo_

protected

Definition at line 78 of file HGCalValidator.h.

Referenced by bookHistograms(), and HGCalValidator().

edm::EDGetTokenT<Density> HGCalValidator::density_

protected

Definition at line 90 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

std::string HGCalValidator::dirName_

private

Definition at line 104 of file HGCalValidator.h.

Referenced by bookHistograms(), and HGCalValidator().

const bool HGCalValidator::doCaloParticlePlots_

protected

Definition at line 69 of file HGCalValidator.h.

Referenced by bookHistograms(), and cpParametersAndSelection().

const bool HGCalValidator::doCaloParticleSelection_

protected

Definition at line 70 of file HGCalValidator.h.

Referenced by cpParametersAndSelection().

const bool HGCalValidator::doLayerClustersPlots_

protected

Definition at line 73 of file HGCalValidator.h.

Referenced by bookHistograms(), and dqmAnalyze().

const bool HGCalValidator::doSimClustersPlots_

protected

Definition at line 71 of file HGCalValidator.h.

Referenced by bookHistograms(), and dqmAnalyze().

const bool HGCalValidator::doTrackstersPlots_

protected

Definition at line 75 of file HGCalValidator.h.

Referenced by bookHistograms(), and dqmAnalyze().

std::unique_ptr<HGVHistoProducerAlgo> HGCalValidator::histoProducerAlgo_

protected

Definition at line 95 of file HGCalValidator.h.

Referenced by bookHistograms(), cpParametersAndSelection(), dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<std::unordered_map<DetId, const HGCRecHit*> > HGCalValidator::hitMap_

protected

Definition at line 89 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

std::vector<edm::InputTag> HGCalValidator::label_clustersmask

protected

Definition at line 77 of file HGCalValidator.h.

Referenced by bookHistograms(), dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<std::vector<CaloParticle> > HGCalValidator::label_cp_effic

protected

Definition at line 85 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<std::vector<CaloParticle> > HGCalValidator::label_cp_fake

protected

Definition at line 86 of file HGCalValidator.h.

Referenced by HGCalValidator().

edm::InputTag HGCalValidator::label_layerClustersPlots_

protected

Definition at line 74 of file HGCalValidator.h.

Referenced by bookHistograms().

edm::InputTag HGCalValidator::label_lcl

protected

Definition at line 63 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::InputTag HGCalValidator::label_LCToCPLinking_

protected

Definition at line 74 of file HGCalValidator.h.

Referenced by bookHistograms().

edm::InputTag HGCalValidator::label_SimClustersLevel_

protected

Definition at line 72 of file HGCalValidator.h.

Referenced by bookHistograms().

edm::InputTag HGCalValidator::label_SimClustersPlots_

protected

Definition at line 72 of file HGCalValidator.h.

Referenced by bookHistograms().

edm::InputTag HGCalValidator::label_simTSFromCP

protected

Definition at line 65 of file HGCalValidator.h.

Referenced by HGCalValidator().

std::vector<edm::InputTag> HGCalValidator::label_tst

protected

Definition at line 64 of file HGCalValidator.h.

Referenced by bookHistograms(), dqmAnalyze(), and HGCalValidator().

edm::InputTag HGCalValidator::label_TSToCPLinking_

protected

Definition at line 76 of file HGCalValidator.h.

Referenced by bookHistograms().

std::vector<edm::EDGetTokenT<ticl::TracksterCollection> > HGCalValidator::label_tstTokens

protected

Definition at line 83 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

std::vector<edm::EDGetTokenT<reco::CaloClusterCollection> > HGCalValidator::labelToken

protected

Definition at line 80 of file HGCalValidator.h.

edm::EDGetTokenT<reco::CaloClusterCollection> HGCalValidator::layerclusters_

protected

Definition at line 82 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

std::vector<int> HGCalValidator::particles_to_monitor_

private

Definition at line 101 of file HGCalValidator.h.

Referenced by bookHistograms(), and HGCalValidator().

const bool HGCalValidator::SaveGeneralInfo_

protected

Definition at line 68 of file HGCalValidator.h.

Referenced by bookHistograms(), and dqmAnalyze().

edm::EDGetTokenT<std::vector<SimCluster> > HGCalValidator::simClusters_

protected

Definition at line 81 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<ticl::TracksterCollection> HGCalValidator::simTrackstersFromCPs_

protected

Definition at line 84 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

edm::EDGetTokenT<std::vector<SimVertex> > HGCalValidator::simVertices_

protected

Definition at line 87 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

std::vector<int> HGCalValidator::thicknesses_to_monitor_

private

Definition at line 103 of file HGCalValidator.h.

Referenced by bookHistograms(), dqmAnalyze(), and HGCalValidator().

std::shared_ptr<hgcal::RecHitTools> HGCalValidator::tools_

private

Definition at line 99 of file HGCalValidator.h.

Referenced by dqmAnalyze(), and HGCalValidator().

unsigned HGCalValidator::totallayers_to_monitor_

private

Definition at line 102 of file HGCalValidator.h.

Referenced by bookHistograms(), dqmAnalyze(), and HGCalValidator().