HGCalShowerSeparation Class Reference

Inheritance diagram for HGCalShowerSeparation:

Public Member Functions
	HGCalShowerSeparation (const edm::ParameterSet &)
	~HGCalShowerSeparation () override=default
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	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
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	fillWithRecHits (std::unordered_map< DetId, const HGCRecHit *> &, DetId, unsigned int, float, int &, float &)

Private Attributes
edm::EDGetTokenT< std::vector< CaloParticle > >	caloParticles_
std::vector< MonitorElement * >	centers_
int	debug_
MonitorElement *	deltaEtaPhi_
std::vector< MonitorElement * >	distanceOnLayer_
MonitorElement *	energy1_
MonitorElement *	energy2_
MonitorElement *	energytot_
MonitorElement *	eta1_
MonitorElement *	eta2_
MonitorElement *	etaPhi_
bool	filterOnEnergyAndCaloP_
std::vector< MonitorElement * >	globalProfileOnLayer_
edm::EDGetTokenT< std::unordered_map< DetId, const HGCRecHit * > >	hitMap_
std::vector< MonitorElement * >	idealDeltaXY_
std::vector< MonitorElement * >	idealDistanceOnLayer_
MonitorElement *	layerDistance_
MonitorElement *	layerEnergy_
std::vector< MonitorElement * >	profileOnLayer_
hgcal::RecHitTools	recHitTools_
MonitorElement *	scEnergy_
MonitorElement *	showerProfile_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

Static Private Attributes
static constexpr int	layers_ = 52

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 36 of file HGCalShowerSeparation.cc.

Constructor & Destructor Documentation

HGCalShowerSeparation()

HGCalShowerSeparation::HGCalShowerSeparation ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 78 of file HGCalShowerSeparation.cc.

References caloTruthCellsProducer_cfi::caloParticles, caloParticles_, edm::ParameterSet::getParameter(), and hitMap_.

    : tok_geom_(esConsumes<CaloGeometry, CaloGeometryRecord>()),
      debug_(iConfig.getParameter<int>("debug")),
      filterOnEnergyAndCaloP_(iConfig.getParameter<bool>("filterOnEnergyAndCaloP")) {
  auto hitMapInputTag = iConfig.getParameter<edm::InputTag>("hitMapTag");
  auto caloParticles = iConfig.getParameter<edm::InputTag>("caloParticles");
  hitMap_ = consumes<std::unordered_map<DetId, const HGCRecHit*>>(hitMapInputTag);
  caloParticles_ = consumes<std::vector<CaloParticle>>(caloParticles);
}

~HGCalShowerSeparation()

HGCalShowerSeparation::~HGCalShowerSeparation ( )

overridedefault

Member Function Documentation

analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 148 of file HGCalShowerSeparation.cc.

References funct::abs(), caloTruthCellsProducer_cfi::caloParticles, caloParticles_, centers_, funct::cos(), submitPVResolutionJobs::count, debug_, deltaEtaPhi_, HLT_2023v12_cff::distance, distanceOnLayer_, HCALHighEnergyHPDFilter_cfi::energy, SimCluster::energy(), energy1_, energy2_, energytot_, HLT_2023v12_cff::eta1, eta1_, HLT_2023v12_cff::eta2, eta2_, etaPhi_, JetChargeProducer_cfi::exp, dqm::impl::MonitorElement::Fill(), filterOnEnergyAndCaloP_, edm::EventSetup::getData(), hgcal::RecHitTools::getLayerWithOffset(), hgcal::RecHitTools::getPosition(), globalProfileOnLayer_, hitMap_, SimCluster::hits_and_fractions(), idealDeltaXY_, idealDistanceOnLayer_, iEvent, IfLogTrace, layerDistance_, layerEnergy_, or, edm::Handle< T >::product(), profileOnLayer_, recHitTools_, scEnergy_, hgcal::RecHitTools::setGeometry(), showerProfile_, SimCluster::simEnergy(), funct::sin(), edm::RefVector< C, T, F >::size(), findQualityFiles::size, mathSSE::sqrt(), funct::tan(), tok_geom_, PV3DBase< T, PVType, FrameType >::x(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

                                                                                       {
  recHitTools_.setGeometry(iSetup.getData(tok_geom_));

  const edm::Handle<std::vector<CaloParticle>>& caloParticleHandle = iEvent.getHandle(caloParticles_);
  const std::vector<CaloParticle>& caloParticles = *(caloParticleHandle.product());

  edm::Handle<std::unordered_map<DetId, const HGCRecHit*>> hitMapHandle;
  iEvent.getByToken(hitMap_, hitMapHandle);
  const auto hitmap = *hitMapHandle;

  // loop over caloParticles
  IfLogTrace(debug_ > 0, "HGCalShowerSeparation") << "Number of caloParticles: " << caloParticles.size() << std::endl;
  if (caloParticles.size() == 2) {
    auto eta1 = caloParticles[0].eta();
    auto phi1 = caloParticles[0].phi();
    auto theta1 = 2. * std::atan(exp(-eta1));
    auto eta2 = caloParticles[1].eta();
    auto phi2 = caloParticles[1].phi();
    auto theta2 = 2. * std::atan(exp(-eta2));
    eta1_->Fill(eta1);
    eta2_->Fill(eta2);

    // Select event only if the sum of the energy of its recHits
    // is close enough to the gen energy
    int count = 0;
    int size = 0;
    float energy = 0.;
    float energy_tmp = 0.;
    for (const auto& it_caloPart : caloParticles) {
      count++;
      const SimClusterRefVector& simClusterRefVector = it_caloPart.simClusters();
      size += simClusterRefVector.size();
      for (const auto& it_sc : simClusterRefVector) {
        const SimCluster& simCluster = (*(it_sc));
        const std::vector<std::pair<uint32_t, float>>& hits_and_fractions = simCluster.hits_and_fractions();
        for (const auto& it_haf : hits_and_fractions) {
          if (hitmap.count(it_haf.first))
            energy += hitmap.at(it_haf.first)->energy() * it_haf.second;
        }  //hits and fractions
      }    // simcluster
      if (count == 1) {
        energy1_->Fill(energy);
        energy_tmp = energy;
      } else {
        energy2_->Fill(energy - energy_tmp);
      }
    }  // caloParticle
    energytot_->Fill(energy);
    if (filterOnEnergyAndCaloP_ && (energy < 2. * 0.8 * 80 or size != 2))
      return;

    deltaEtaPhi_->Fill(eta1 - eta2, phi1 - phi2);

    for (const auto& it_caloPart : caloParticles) {
      const SimClusterRefVector& simClusterRefVector = it_caloPart.simClusters();
      IfLogTrace(debug_ > 0, "HGCalShowerSeparation") << ">>> " << simClusterRefVector.size() << std::endl;
      for (const auto& it_sc : simClusterRefVector) {
        const SimCluster& simCluster = (*(it_sc));
        if (simCluster.energy() < 80 * 0.8)
          continue;
        scEnergy_->Fill(simCluster.energy());
        IfLogTrace(debug_ > 1, "HGCalShowerSeparation")
            << ">>> SC.energy(): " << simCluster.energy() << " SC.simEnergy(): " << simCluster.simEnergy() << std::endl;
        const std::vector<std::pair<uint32_t, float>>& hits_and_fractions = simCluster.hits_and_fractions();

        for (const auto& it_haf : hits_and_fractions) {
          if (!hitmap.count(it_haf.first))
            continue;
          unsigned int hitlayer = recHitTools_.getLayerWithOffset(it_haf.first);
          auto global = recHitTools_.getPosition(it_haf.first);
          float globalx = global.x();
          float globaly = global.y();
          float globalz = global.z();
          if (globalz == 0)
            continue;
          auto rho1 = globalz * tan(theta1);
          auto rho2 = globalz * tan(theta2);
          auto x1 = rho1 * cos(phi1);
          auto y1 = rho1 * sin(phi1);
          auto x2 = rho2 * cos(phi2);
          auto y2 = rho2 * sin(phi2);
          auto half_point_x = (x1 + x2) / 2.;
          auto half_point_y = (y1 + y2) / 2.;
          auto half_point = sqrt((x1 - half_point_x) * (x1 - half_point_x) + (y1 - half_point_y) * (y1 - half_point_y));
          auto d_len = sqrt((x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1));
          auto dn_x = (x2 - x1) / d_len;
          auto dn_y = (y2 - y1) / d_len;
          auto distance = (globalx - x1) * dn_x + (globaly - y1) * dn_y;
          distance -= half_point;
          auto idealDistance = sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
          if (hitmap.count(it_haf.first)) {
            profileOnLayer_[hitlayer]->Fill(10. * (globalx - half_point_x),
                                            10. * (globaly - half_point_y),
                                            hitmap.at(it_haf.first)->energy() * it_haf.second);
            profileOnLayer_[55]->Fill(10. * (globalx - half_point_x),
                                      10. * (globaly - half_point_y),
                                      hitmap.at(it_haf.first)->energy() * it_haf.second);
            globalProfileOnLayer_[hitlayer]->Fill(globalx, globaly, hitmap.at(it_haf.first)->energy() * it_haf.second);
            globalProfileOnLayer_[55]->Fill(globalx, globaly, hitmap.at(it_haf.first)->energy() * it_haf.second);
            layerEnergy_->Fill(hitlayer, hitmap.at(it_haf.first)->energy());
            layerDistance_->Fill(hitlayer, std::abs(10. * distance), hitmap.at(it_haf.first)->energy() * it_haf.second);
            etaPhi_->Fill(global.eta(), global.phi());
            distanceOnLayer_[hitlayer]->Fill(10. * distance);                  //,
            idealDistanceOnLayer_[hitlayer]->Fill(10. * idealDistance);        //,
            idealDeltaXY_[hitlayer]->Fill(10. * (x1 - x2), 10. * (y1 - y2));   //,
            centers_[hitlayer]->Fill(10. * half_point_x, 10. * half_point_y);  //,
            IfLogTrace(debug_ > 0, "HGCalShowerSeparation")
                << ">>> " << distance << " " << hitlayer << " " << hitmap.at(it_haf.first)->energy() * it_haf.second
                << std::endl;
            showerProfile_->Fill(10. * distance, hitlayer, hitmap.at(it_haf.first)->energy() * it_haf.second);
          }
        }  // end simHit
      }    // end simCluster
    }      // end caloparticle
  }
}

bookHistograms()

void HGCalShowerSeparation::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &    )

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 88 of file HGCalShowerSeparation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), centers_, deltaEtaPhi_, distanceOnLayer_, energy1_, energy2_, energytot_, eta1_, eta2_, etaPhi_, globalProfileOnLayer_, mps_fire::i, idealDeltaXY_, idealDistanceOnLayer_, layerDistance_, layerEnergy_, layers_, profileOnLayer_, scEnergy_, dqm::implementation::NavigatorBase::setCurrentFolder(), showerProfile_, AlCaHLTBitMon_QueryRunRegistry::string, and cond::impl::to_string().

                                                                                                             {
  ibooker.cd();
  ibooker.setCurrentFolder("HGCalShowerSeparation");
  scEnergy_ = ibooker.book1D("SCEnergy", "SCEnergy", 240, 0., 120.);
  eta1_ = ibooker.book1D("eta1", "eta1", 80, 0., 4.);
  eta2_ = ibooker.book1D("eta2", "eta2", 80, 0., 4.);
  energy1_ = ibooker.book1D("energy1", "energy1", 240, 0., 120.);
  energy2_ = ibooker.book1D("energy2", "energy2", 240, 0., 120.);
  energytot_ = ibooker.book1D("energytot", "energytot", 200, 100., 200.);
  showerProfile_ = ibooker.book2D("ShowerProfile", "ShowerProfile", 800, -400., 400., layers_, 0., (float)layers_);
  layerEnergy_ = ibooker.book2D("LayerEnergy", "LayerEnergy", 60, 0., 60., 50, 0., 0.1);
  layerDistance_ = ibooker.book2D("LayerDistance", "LayerDistance", 60, 0., 60., 400, -400., 400.);
  etaPhi_ = ibooker.book2D("EtaPhi", "EtaPhi", 800, -4., 4., 800, -4., 4.);
  deltaEtaPhi_ = ibooker.book2D("DeltaEtaPhi", "DeltaEtaPhi", 100, -0.5, 0.5, 100, -0.5, 0.5);
  for (int i = 0; i < layers_; ++i) {
    profileOnLayer_.push_back(ibooker.book2D(std::string("ProfileOnLayer_") + std::to_string(i),
                                             std::string("ProfileOnLayer_") + std::to_string(i),
                                             120,
                                             -600.,
                                             600.,
                                             120,
                                             -600.,
                                             600.));
    globalProfileOnLayer_.push_back(ibooker.book2D(std::string("GlobalProfileOnLayer_") + std::to_string(i),
                                                   std::string("GlobalProfileOnLayer_") + std::to_string(i),
                                                   320,
                                                   -160.,
                                                   160.,
                                                   320,
                                                   -160.,
                                                   160.));
    distanceOnLayer_.push_back(ibooker.book1D(std::string("DistanceOnLayer_") + std::to_string(i),
                                              std::string("DistanceOnLayer_") + std::to_string(i),
                                              120,
                                              -600.,
                                              600.));
    idealDistanceOnLayer_.push_back(ibooker.book1D(std::string("IdealDistanceOnLayer_") + std::to_string(i),
                                                   std::string("IdealDistanceOnLayer_") + std::to_string(i),
                                                   120,
                                                   -600.,
                                                   600.));
    idealDeltaXY_.push_back(ibooker.book2D(std::string("IdealDeltaXY_") + std::to_string(i),
                                           std::string("IdealDeltaXY_") + std::to_string(i),
                                           800,
                                           -400.,
                                           400.,
                                           800,
                                           -400.,
                                           400.));
    centers_.push_back(ibooker.book2D(std::string("Centers_") + std::to_string(i),
                                      std::string("Centers_") + std::to_string(i),
                                      320,
                                      -1600.,
                                      1600.,
                                      320,
                                      -1600.,
                                      1600.));
  }
}

fillDescriptions()

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

static

Definition at line 267 of file HGCalShowerSeparation.cc.

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

                                                                                       {
  edm::ParameterSetDescription desc;
  desc.add<int>("debug", 1);
  desc.add<bool>("filterOnEnergyAndCaloP", false);
  desc.add<edm::InputTag>("caloParticles", edm::InputTag("mix", "MergedCaloTruth"));
  desc.add<edm::InputTag>("hitMapTag", edm::InputTag("hgcalRecHitMapProducer"));
  descriptions.add("hgcalShowerSeparationDefault", desc);
}

fillWithRecHits()

void HGCalShowerSeparation::fillWithRecHits ( std::unordered_map< DetId, const HGCRecHit *> &  , DetId  , unsigned  int, float  , int &  , float &    )

private

Member Data Documentation

caloParticles_

edm::EDGetTokenT<std::vector<CaloParticle> > HGCalShowerSeparation::caloParticles_

private

Definition at line 50 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and HGCalShowerSeparation().

centers_

std::vector<MonitorElement*> HGCalShowerSeparation::centers_

private

Definition at line 73 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

debug_

int HGCalShowerSeparation::debug_

private

Definition at line 53 of file HGCalShowerSeparation.cc.

Referenced by analyze().

deltaEtaPhi_

MonitorElement* HGCalShowerSeparation::deltaEtaPhi_

private

Definition at line 67 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

distanceOnLayer_

std::vector<MonitorElement*> HGCalShowerSeparation::distanceOnLayer_

private

Definition at line 70 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

energy1_

MonitorElement* HGCalShowerSeparation::energy1_

private

Definition at line 59 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

energy2_

MonitorElement* HGCalShowerSeparation::energy2_

private

Definition at line 60 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

energytot_

MonitorElement* HGCalShowerSeparation::energytot_

private

Definition at line 61 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

eta1_

MonitorElement* HGCalShowerSeparation::eta1_

private

Definition at line 57 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

eta2_

MonitorElement* HGCalShowerSeparation::eta2_

private

Definition at line 58 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

etaPhi_

MonitorElement* HGCalShowerSeparation::etaPhi_

private

Definition at line 66 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

filterOnEnergyAndCaloP_

bool HGCalShowerSeparation::filterOnEnergyAndCaloP_

private

Definition at line 54 of file HGCalShowerSeparation.cc.

Referenced by analyze().

globalProfileOnLayer_

std::vector<MonitorElement*> HGCalShowerSeparation::globalProfileOnLayer_

private

Definition at line 69 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

hitMap_

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

private

Definition at line 49 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and HGCalShowerSeparation().

idealDeltaXY_

std::vector<MonitorElement*> HGCalShowerSeparation::idealDeltaXY_

private

Definition at line 72 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

idealDistanceOnLayer_

std::vector<MonitorElement*> HGCalShowerSeparation::idealDistanceOnLayer_

private

Definition at line 71 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

layerDistance_

MonitorElement* HGCalShowerSeparation::layerDistance_

private

Definition at line 65 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

layerEnergy_

MonitorElement* HGCalShowerSeparation::layerEnergy_

private

Definition at line 64 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

layers_

constexpr int HGCalShowerSeparation::layers_ = 52

staticprivate

Definition at line 75 of file HGCalShowerSeparation.cc.

Referenced by bookHistograms().

profileOnLayer_

std::vector<MonitorElement*> HGCalShowerSeparation::profileOnLayer_

private

Definition at line 68 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

recHitTools_

hgcal::RecHitTools HGCalShowerSeparation::recHitTools_

private

Definition at line 55 of file HGCalShowerSeparation.cc.

Referenced by analyze().

scEnergy_

MonitorElement* HGCalShowerSeparation::scEnergy_

private

Definition at line 62 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

showerProfile_

MonitorElement* HGCalShowerSeparation::showerProfile_

private

Definition at line 63 of file HGCalShowerSeparation.cc.

Referenced by analyze(), and bookHistograms().

tok_geom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HGCalShowerSeparation::tok_geom_

private

Definition at line 51 of file HGCalShowerSeparation.cc.

Referenced by analyze().