Skip to content

Commit

Permalink
Add Et distributions for BC & SC energy ECAL DQM plots
Browse files Browse the repository at this point in the history
  • Loading branch information
Kyungmin Park committed Nov 26, 2024
1 parent e601d0d commit 43c39f5
Show file tree
Hide file tree
Showing 2 changed files with 91 additions and 2 deletions.
75 changes: 75 additions & 0 deletions DQM/EcalMonitorTasks/python/ClusterTask_cfi.py
Original file line number Diff line number Diff line change
Expand Up @@ -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'),
Expand Down Expand Up @@ -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'),
Expand Down Expand Up @@ -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'),
Expand All @@ -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'),
Expand Down Expand Up @@ -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'),
Expand All @@ -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'),
Expand Down
18 changes: 16 additions & 2 deletions DQM/EcalMonitorTasks/src/ClusterTask.cc
Original file line number Diff line number Diff line change
Expand Up @@ -209,7 +209,9 @@ namespace ecaldqm {

void ClusterTask::runOnBasicClusters(edm::View<reco::CaloCluster> 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"));
Expand Down Expand Up @@ -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);
Expand Down Expand Up @@ -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"));
Expand Down Expand Up @@ -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(
Expand All @@ -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());
Expand Down

0 comments on commit 43c39f5

Please sign in to comment.