From 43c39f5dbdc2f5cba68781087dfa88908ca0b686 Mon Sep 17 00:00:00 2001 From: Kyungmin Park Date: Tue, 26 Nov 2024 23:12:14 +0100 Subject: [PATCH] Add Et distributions for BC & SC energy ECAL DQM plots --- .../python/ClusterTask_cfi.py | 75 +++++++++++++++++++ DQM/EcalMonitorTasks/src/ClusterTask.cc | 18 ++++- 2 files changed, 91 insertions(+), 2 deletions(-) diff --git a/DQM/EcalMonitorTasks/python/ClusterTask_cfi.py b/DQM/EcalMonitorTasks/python/ClusterTask_cfi.py index cdc3f7c68949c..1ea1561be6888 100644 --- a/DQM/EcalMonitorTasks/python/ClusterTask_cfi.py +++ b/DQM/EcalMonitorTasks/python/ClusterTask_cfi.py @@ -103,6 +103,19 @@ btype = cms.untracked.string('User'), description = cms.untracked.string('Basic cluster energy distribution.') ), + BCEt = cms.untracked.PSet( + path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT BC transverse energy'), + kind = cms.untracked.string('TH1F'), + otype = cms.untracked.string('Ecal2P'), + xaxis = cms.untracked.PSet( + high = cms.untracked.double(150.0), + nbins = cms.untracked.int32(50), + low = cms.untracked.double(0.0), + title = cms.untracked.string('energy (GeV)') + ), + btype = cms.untracked.string('User'), + description = cms.untracked.string('Basic cluster transverse energy distribution.') + ), BCSizeMap = cms.untracked.PSet( path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT BC size map%(suffix)s'), kind = cms.untracked.string('TProfile2D'), @@ -152,6 +165,16 @@ btype = cms.untracked.string('SuperCrystal'), description = cms.untracked.string('2D distribution of the mean energy of the basic clusters.') ), + BCEtMap = cms.untracked.PSet( + path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT BC transverse energy map%(suffix)s'), + kind = cms.untracked.string('TProfile2D'), + zaxis = cms.untracked.PSet( + title = cms.untracked.string('energy (GeV)') + ), + otype = cms.untracked.string('Ecal3P'), + btype = cms.untracked.string('SuperCrystal'), + description = cms.untracked.string('2D distribution of the mean transverse energy of the basic clusters.') + ), BCEtMapProjEta = cms.untracked.PSet( path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT BC ET projection eta%(suffix)s'), kind = cms.untracked.string('TProfile'), @@ -272,6 +295,19 @@ btype = cms.untracked.string('User'), description = cms.untracked.string('Super cluster energy distribution.') ), + SCEt = cms.untracked.PSet( + path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC transverse energy'), + kind = cms.untracked.string('TH1F'), + otype = cms.untracked.string('Ecal2P'), + xaxis = cms.untracked.PSet( + high = cms.untracked.double(150.0), + nbins = cms.untracked.int32(50), + low = cms.untracked.double(0.0), + title = cms.untracked.string('energy (GeV)') + ), + btype = cms.untracked.string('User'), + description = cms.untracked.string('Super cluster transverse energy distribution.') + ), SCRawE = cms.untracked.PSet( path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC raw energy'), kind = cms.untracked.string('TH1F'), @@ -285,6 +321,19 @@ btype = cms.untracked.string('User'), description = cms.untracked.string('Super cluster raw energy distribution.') ), + SCRawEt = cms.untracked.PSet( + path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC raw transverse energy'), + kind = cms.untracked.string('TH1F'), + otype = cms.untracked.string('Ecal2P'), + xaxis = cms.untracked.PSet( + high = cms.untracked.double(150.0), + nbins = cms.untracked.int32(50), + low = cms.untracked.double(0.0), + title = cms.untracked.string('energy (GeV)') + ), + btype = cms.untracked.string('User'), + description = cms.untracked.string('Super cluster raw transverse energy distribution.') + ), SCNcrystals = cms.untracked.PSet( path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC size (crystal)'), kind = cms.untracked.string('TH1F'), @@ -336,6 +385,19 @@ btype = cms.untracked.string('User'), description = cms.untracked.string('Energy distribution (raw energy) of the super clusters (low scale).') ), + SCEtLow = cms.untracked.PSet( + path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC transverse energy (low scale)'), + kind = cms.untracked.string('TH1F'), + otype = cms.untracked.string('Ecal2P'), + xaxis = cms.untracked.PSet( + high = cms.untracked.double(10.0), + nbins = cms.untracked.int32(50), + low = cms.untracked.double(0.0), + title = cms.untracked.string('energy (GeV)') + ), + btype = cms.untracked.string('User'), + description = cms.untracked.string('Transverse energy distribution (raw energy) of the super clusters (low scale).') + ), SCRawELow = cms.untracked.PSet( path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC raw energy (low scale)'), kind = cms.untracked.string('TH1F'), @@ -349,6 +411,19 @@ btype = cms.untracked.string('User'), description = cms.untracked.string('Energy distribution of the super clusters (low scale).') ), + SCRawEtLow = cms.untracked.PSet( + path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC raw transverse energy (low scale)'), + kind = cms.untracked.string('TH1F'), + otype = cms.untracked.string('Ecal2P'), + xaxis = cms.untracked.PSet( + high = cms.untracked.double(10.0), + nbins = cms.untracked.int32(50), + low = cms.untracked.double(0.0), + title = cms.untracked.string('energy (GeV)') + ), + btype = cms.untracked.string('User'), + description = cms.untracked.string('Transverse energy distribution of the super clusters (low scale).') + ), SCRawEHigh = cms.untracked.PSet( path = cms.untracked.string('%(subdet)s/%(prefix)sClusterTask/%(prefix)sCLT SC raw energy (high scale)'), kind = cms.untracked.string('TH1F'), diff --git a/DQM/EcalMonitorTasks/src/ClusterTask.cc b/DQM/EcalMonitorTasks/src/ClusterTask.cc index 7ac9405ff99bd..fdbd3846d6361 100644 --- a/DQM/EcalMonitorTasks/src/ClusterTask.cc +++ b/DQM/EcalMonitorTasks/src/ClusterTask.cc @@ -209,7 +209,9 @@ namespace ecaldqm { void ClusterTask::runOnBasicClusters(edm::View const& _bcs, Collections _collection) { MESet& meBCE(MEs_.at("BCE")); + MESet& meBCEt(MEs_.at("BCEt")); MESet& meBCEMap(MEs_.at("BCEMap")); + MESet& meBCEtMap(MEs_.at("BCEtMap")); MESet& meBCEMapProjEta(MEs_.at("BCEMapProjEta")); MESet& meBCEMapProjPhi(MEs_.at("BCEMapProjPhi")); MESet& meBCEtMapProjEta(MEs_.at("BCEtMapProjEta")); @@ -261,8 +263,10 @@ namespace ecaldqm { subdet = -EcalEndcap; meBCE.fill(getEcalDQMSetupObjects(), id, energy); + meBCEt.fill(getEcalDQMSetupObjects(), id, et); meBCEMap.fill(getEcalDQMSetupObjects(), id, energy); + meBCEtMap.fill(getEcalDQMSetupObjects(), id, et); meBCEMapProjEta.fill(getEcalDQMSetupObjects(), posEta, energy); meBCEMapProjPhi.fill(getEcalDQMSetupObjects(), subdet, posPhi, energy); meBCEtMapProjEta.fill(getEcalDQMSetupObjects(), posEta, et); @@ -361,9 +365,13 @@ namespace ecaldqm { EcalSubdetector subdet(isBarrel ? EcalBarrel : EcalEndcap); MESet& meSCE(MEs_.at("SCE")); + MESet& meSCEt(MEs_.at("SCEt")); MESet& meSCELow(MEs_.at("SCELow")); + MESet& meSCEtLow(MEs_.at("SCEtLow")); MESet& meSCRawE(MEs_.at("SCRawE")); + MESet& meSCRawEt(MEs_.at("SCRawEt")); MESet& meSCRawELow(MEs_.at("SCRawELow")); + MESet& meSCRawEtLow(MEs_.at("SCRawEtLow")); MESet& meSCRawEHigh(MEs_.at("SCRawEHigh")); MESet& meSCNBCs(MEs_.at("SCNBCs")); MESet& meSCNcrystals(MEs_.at("SCNcrystals")); @@ -395,9 +403,8 @@ namespace ecaldqm { for (reco::SuperClusterCollection::const_iterator scItr(_scs.begin()); scItr != _scs.end(); ++scItr) { DetId seedId(scItr->seed()->seed()); + math::XYZPoint const& position(scItr->position()); if (seedId.null()) { - math::XYZPoint const& position(scItr->position()); - GlobalPoint gp(position.x(), position.y(), position.z()); CaloSubdetectorGeometry const* subgeom( @@ -417,13 +424,20 @@ namespace ecaldqm { float energy(scItr->energy()); float rawEnergy(scItr->rawEnergy()); + float posEta(position.eta()); + float et(energy / std::cosh(posEta)); + float rawEt(rawEnergy / std::cosh(posEta)); float size(scItr->size()); meSCE.fill(getEcalDQMSetupObjects(), seedId, energy); + meSCEt.fill(getEcalDQMSetupObjects(), seedId, et); meSCELow.fill(getEcalDQMSetupObjects(), seedId, energy); + meSCEtLow.fill(getEcalDQMSetupObjects(), seedId, et); meSCRawE.fill(getEcalDQMSetupObjects(), seedId, rawEnergy); + meSCRawEt.fill(getEcalDQMSetupObjects(), seedId, rawEt); meSCRawELow.fill(getEcalDQMSetupObjects(), seedId, rawEnergy); + meSCRawEtLow.fill(getEcalDQMSetupObjects(), seedId, rawEt); meSCRawEHigh.fill(getEcalDQMSetupObjects(), seedId, rawEnergy); meSCNBCs.fill(getEcalDQMSetupObjects(), seedId, scItr->clustersSize());