Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

MTD validation: update cluster histograms #37056

Merged
merged 6 commits into from Feb 27, 2022
Merged

MTD validation: update cluster histograms #37056

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
e6b8fa3
Add BTL cluster error monitoring
fabiocos Feb 3, 2022
c1d7906
Update BTL cluster monitoring
fabiocos Jan 10, 2022
df84240
Add position plots for ETL, and fixes
fabiocos Feb 4, 2022
f641616
Move space resolution plots to local coordinates and under special flag
fabiocos Feb 8, 2022
2c51320
Keep the minimal histograms for unmatched clusters
fabiocos Feb 24, 2022
5ce5f7c
Move unmatched clauster energy to log10
fabiocos Feb 25, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
111 changes: 99 additions & 12 deletions Validation/MtdValidation/plugins/BtlLocalRecoValidation.cc
View file
Open in desktop
Expand Up @@ -26,6 +26,7 @@
#include "DataFormats/ForwardDetId/interface/BTLDetId.h"
#include "DataFormats/FTLRecHit/interface/FTLRecHitCollections.h"
#include "DataFormats/FTLRecHit/interface/FTLClusterCollections.h"
#include "DataFormats/TrackerRecHit2D/interface/MTDTrackingRecHit.h"

#include "SimDataFormats/CrossingFrame/interface/CrossingFrame.h"
#include "SimDataFormats/CrossingFrame/interface/MixCollection.h"
Expand Down Expand Up @@ -61,6 +62,8 @@ class BtlLocalRecoValidation : public DQMEDAnalyzer {

void analyze(const edm::Event&, const edm::EventSetup&) override;

bool isSameCluster(const FTLCluster&, const FTLCluster&);

// ------------ member data ------------

const std::string folder_;
Expand All @@ -73,6 +76,7 @@ class BtlLocalRecoValidation : public DQMEDAnalyzer {
edm::EDGetTokenT<FTLUncalibratedRecHitCollection> btlUncalibRecHitsToken_;
edm::EDGetTokenT<CrossingFrame<PSimHit> > btlSimHitsToken_;
edm::EDGetTokenT<FTLClusterCollection> btlRecCluToken_;
edm::EDGetTokenT<MTDTrackingDetSetVector> mtdTrackingHitToken_;

edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> mtdgeoToken_;
edm::ESGetToken<MTDTopology, MTDTopologyRcd> mtdtopoToken_;
Expand Down Expand Up @@ -126,18 +130,29 @@ class BtlLocalRecoValidation : public DQMEDAnalyzer {
MonitorElement* meCluEta_;
MonitorElement* meCluHits_;
MonitorElement* meCluZvsPhi_;
MonitorElement* meCluEnergyvsEta_;
MonitorElement* meCluHitsvsEta_;

MonitorElement* meCluTimeRes_;
MonitorElement* meCluEnergyRes_;
MonitorElement* meCluTResvsE_;
MonitorElement* meCluTResvsEta_;
MonitorElement* meCluTPullvsE_;
MonitorElement* meCluTPullvsEta_;
MonitorElement* meCluRhoRes_;
MonitorElement* meCluPhiRes_;
MonitorElement* meCluXRes_;
MonitorElement* meCluYRes_;
MonitorElement* meCluLocalXRes_;
MonitorElement* meCluLocalYRes_;
MonitorElement* meCluZRes_;
MonitorElement* meCluLocalXPull_;
MonitorElement* meCluLocalYPull_;
MonitorElement* meCluZPull_;
MonitorElement* meCluYXLocal_;
MonitorElement* meCluYXLocalSim_;
MonitorElement* meCluXLocalErr_;
MonitorElement* meCluYLocalErr_;

MonitorElement* meUnmatchedCluEnergy_;

// --- UncalibratedRecHits histograms

Expand All @@ -158,6 +173,11 @@ class BtlLocalRecoValidation : public DQMEDAnalyzer {
MonitorElement* meTimeResEtavsQ_[nBinsEta_][nBinsEtaQ_];
};

bool BtlLocalRecoValidation::isSameCluster(const FTLCluster& clu1, const FTLCluster& clu2) {
return clu1.id() == clu2.id() && clu1.size() == clu2.size() && clu1.x() == clu2.x() && clu1.y() == clu2.y() &&
clu1.time() == clu2.time();
}

// ------------ constructor and destructor --------------
BtlLocalRecoValidation::BtlLocalRecoValidation(const edm::ParameterSet& iConfig)
: folder_(iConfig.getParameter<std::string>("folder")),
Expand All @@ -171,6 +191,7 @@ BtlLocalRecoValidation::BtlLocalRecoValidation(const edm::ParameterSet& iConfig)
consumes<FTLUncalibratedRecHitCollection>(iConfig.getParameter<edm::InputTag>("uncalibRecHitsTag"));
btlSimHitsToken_ = consumes<CrossingFrame<PSimHit> >(iConfig.getParameter<edm::InputTag>("simHitsTag"));
btlRecCluToken_ = consumes<FTLClusterCollection>(iConfig.getParameter<edm::InputTag>("recCluTag"));
mtdTrackingHitToken_ = consumes<MTDTrackingDetSetVector>(iConfig.getParameter<edm::InputTag>("trkHitTag"));

mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
Expand All @@ -193,6 +214,7 @@ void BtlLocalRecoValidation::analyze(const edm::Event& iEvent, const edm::EventS
auto btlRecHitsHandle = makeValid(iEvent.getHandle(btlRecHitsToken_));
auto btlSimHitsHandle = makeValid(iEvent.getHandle(btlSimHitsToken_));
auto btlRecCluHandle = makeValid(iEvent.getHandle(btlRecCluToken_));
auto mtdTrkHitHandle = makeValid(iEvent.getHandle(mtdTrackingHitToken_));
MixCollection<PSimHit> btlSimHits(btlSimHitsHandle.product());

// --- Loop over the BTL SIM hits
Expand Down Expand Up @@ -373,6 +395,22 @@ void BtlLocalRecoValidation::analyze(const edm::Event& iEvent, const edm::EventS

} // ihit loop

// Find the MTDTrackingRecHit corresponding to the cluster
MTDTrackingRecHit* comp(nullptr);
bool matchClu = false;
const auto& trkHits = (*mtdTrkHitHandle)[detIdObject];
for (const auto& trkHit : trkHits) {
if (isSameCluster(trkHit.mtdCluster(), cluster)) {
comp = trkHit.clone();
matchClu = true;
break;
}
}
if (!matchClu) {
edm::LogWarning("BtlLocalRecoValidation")
<< "No valid TrackingRecHit corresponding to cluster, detId = " << detIdObject.rawId();
}

// --- Fill the cluster resolution histograms
if (cluTimeSIM > 0. && cluEneSIM > 0.) {
cluTimeSIM /= cluEneSIM;
Expand All @@ -391,23 +429,40 @@ void BtlLocalRecoValidation::analyze(const edm::Event& iEvent, const edm::EventS
meCluRhoRes_->Fill(rho_res);
meCluPhiRes_->Fill(phi_res);

if (LocalPosDebug_) {
float x_res = global_point.x() - cluGlobalPosSIM.x();
float y_res = global_point.y() - cluGlobalPosSIM.y();
float z_res = global_point.z() - cluGlobalPosSIM.z();
float xlocal_res = local_point.x() - cluLocalPosSIM.x();
float ylocal_res = local_point.y() - cluLocalPosSIM.y();
float z_res = global_point.z() - cluGlobalPosSIM.z();

meCluZRes_->Fill(z_res);

meCluXRes_->Fill(x_res);
meCluYRes_->Fill(y_res);
meCluZRes_->Fill(z_res);
if (LocalPosDebug_) {
if (matchClu && comp != nullptr) {
meCluLocalXRes_->Fill(xlocal_res);
meCluLocalYRes_->Fill(ylocal_res);
meCluLocalXPull_->Fill(xlocal_res / std::sqrt(comp->localPositionError().xx()));
meCluLocalYPull_->Fill(ylocal_res / std::sqrt(comp->localPositionError().yy()));
meCluZPull_->Fill(z_res / std::sqrt(comp->globalPositionError().czz()));
meCluXLocalErr_->Fill(std::sqrt(comp->localPositionError().xx()));
meCluYLocalErr_->Fill(std::sqrt(comp->localPositionError().yy()));
}

meCluYXLocal_->Fill(local_point.x(), local_point.y());
meCluYXLocalSim_->Fill(cluLocalPosSIM.x(), cluLocalPosSIM.y());
}

meCluEnergyvsEta_->Fill(std::abs(cluGlobalPosSIM.eta()), cluster.energy());
meCluHitsvsEta_->Fill(std::abs(cluGlobalPosSIM.eta()), cluster.size());

meCluTResvsEta_->Fill(std::abs(cluGlobalPosSIM.eta()), time_res);
meCluTResvsE_->Fill(cluEneSIM, time_res);

meCluTPullvsEta_->Fill(std::abs(cluGlobalPosSIM.eta()), time_res / cluster.timeError());
meCluTPullvsE_->Fill(cluEneSIM, time_res / cluster.timeError());

} // if ( cluTimeSIM > 0. && cluEneSIM > 0. )
else {
meUnmatchedCluEnergy_->Fill(std::log10(cluster.energy()));
}

} // cluster loop

Expand Down Expand Up @@ -586,10 +641,30 @@ void BtlLocalRecoValidation::bookHistograms(DQMStore::IBooker& ibook,
meCluHits_ = ibook.book1D("BtlCluHitNumber", "BTL hits per cluster; Cluster size", 10, 0, 10);
meCluZvsPhi_ = ibook.book2D(
"BtlOccupancy", "BTL cluster Z vs #phi;Z_{RECO} [cm]; #phi_{RECO} [rad]", 144, -260., 260., 50, -3.2, 3.2);
meCluEnergyvsEta_ = ibook.bookProfile(
"BtlCluEnergyVsEta", "BTL cluster energy vs #eta; |#eta_{RECO}|; E_{RECO} [cm]", 30, 0., 1.55, 0., 20., "S");
meCluHitsvsEta_ = ibook.bookProfile(
"BtlCluHitsVsEta", "BTL hits per cluster vs #eta; |#eta_{RECO}|;Cluster size", 30, 0., 1.55, 0., 10., "S");

meCluTimeRes_ = ibook.book1D("BtlCluTimeRes", "BTL cluster time resolution;T_{RECO}-T_{SIM} [ns]", 100, -0.5, 0.5);
meCluEnergyRes_ =
ibook.book1D("BtlCluEnergyRes", "BTL cluster energy resolution;E_{RECO}-E_{SIM} [MeV]", 100, -0.5, 0.5);
meCluTResvsE_ = ibook.bookProfile("BtlCluTResvsE",
"BTL cluster time resolution vs E;E_{SIM} [MeV];(T_{RECO}-T_{SIM}) [ns]",
20,
0.,
20.,
-0.5,
0.5,
"S");
meCluTResvsEta_ = ibook.bookProfile("BtlCluTResvsEta",
"BTL cluster time resolution vs #eta;|#eta_{RECO}|;(T_{RECO}-T_{SIM}) [ns]",
30,
0,
1.55,
-0.5,
0.5,
"S");
meCluTPullvsE_ = ibook.bookProfile("BtlCluTPullvsE",
"BTL cluster time pull vs E;E_{SIM} [MeV];(T_{RECO}-T_{SIM})/#sigma_{T_{RECO}}",
20,
Expand All @@ -611,10 +686,17 @@ void BtlLocalRecoValidation::bookHistograms(DQMStore::IBooker& ibook,
ibook.book1D("BtlCluRhoRes", "BTL cluster #rho resolution;#rho_{RECO}-#rho_{SIM} [cm]", 100, -0.5, 0.5);
meCluPhiRes_ =
ibook.book1D("BtlCluPhiRes", "BTL cluster #phi resolution;#phi_{RECO}-#phi_{SIM} [rad]", 100, -0.03, 0.03);
meCluZRes_ = ibook.book1D("BtlCluZRes", "BTL cluster Z resolution;Z_{RECO}-Z_{SIM} [cm]", 100, -0.2, 0.2);
if (LocalPosDebug_) {
meCluXRes_ = ibook.book1D("BtlCluXRes", "BTL cluster X resolution;X_{RECO}-X_{SIM} [cm]", 100, -3.1, 3.1);
meCluYRes_ = ibook.book1D("BtlCluYRes", "BTL cluster Y resolution;Y_{RECO}-Y_{SIM} [cm]", 100, -3.1, 3.1);
meCluZRes_ = ibook.book1D("BtlCluZRes", "BTL cluster Z resolution;Z_{RECO}-Z_{SIM} [cm]", 100, -0.2, 0.2);
meCluLocalXRes_ =
ibook.book1D("BtlCluLocalXRes", "BTL cluster local X resolution;X_{RECO}-X_{SIM} [cm]", 100, -3.1, 3.1);
meCluLocalYRes_ =
ibook.book1D("BtlCluLocalYRes", "BTL cluster local Y resolution;Y_{RECO}-Y_{SIM} [cm]", 100, -3.1, 3.1);
meCluLocalXPull_ =
ibook.book1D("BtlCluLocalXPull", "BTL cluster local X pull;X_{RECO}-X_{SIM}/sigmaX_[RECO]", 100, -5., 5.);
meCluLocalYPull_ =
ibook.book1D("BtlCluLocalYPull", "BTL cluster local Y pull;Y_{RECO}-Y_{SIM}/sigmaY_[RECO]", 100, -5., 5.);
meCluZPull_ = ibook.book1D("BtlCluZPull", "BTL cluster Z pull;Z_{RECO}-Z_{SIM}/sigmaZ_[RECO]", 100, -5., 5.);
meCluYXLocal_ = ibook.book2D("BtlCluYXLocal",
"BTL cluster local Y vs X;X^{local}_{RECO} [cm];Y^{local}_{RECO} [cm]",
200,
Expand All @@ -631,7 +713,11 @@ void BtlLocalRecoValidation::bookHistograms(DQMStore::IBooker& ibook,
200,
-2.8,
2.8);
meCluXLocalErr_ = ibook.book1D("BtlCluXLocalErr", "BTL cluster X local error;sigmaX_{RECO,loc} [cm]", 30, 0., 3.);
meCluYLocalErr_ = ibook.book1D("BtlCluYLocalErr", "BTL cluster Y local error;sigmaY_{RECO,loc} [cm]", 30, 0., 0.9);
}
meUnmatchedCluEnergy_ =
ibook.book1D("BtlUnmatchedCluEnergy", "BTL unmatched cluster log10(energy);log10(E_{RECO} [MeV])", 5, -3, 2);

// --- UncalibratedRecHits histograms

Expand Down Expand Up @@ -675,6 +761,7 @@ void BtlLocalRecoValidation::fillDescriptions(edm::ConfigurationDescriptions& de
desc.add<edm::InputTag>("uncalibRecHitsTag", edm::InputTag("mtdUncalibratedRecHits", "FTLBarrel"));
desc.add<edm::InputTag>("simHitsTag", edm::InputTag("mix", "g4SimHitsFastTimerHitsBarrel"));
desc.add<edm::InputTag>("recCluTag", edm::InputTag("mtdClusters", "FTLBarrel"));
desc.add<edm::InputTag>("trkHitTag", edm::InputTag("mtdTrackingRecHits"));
desc.add<double>("HitMinimumEnergy", 1.); // [MeV]
desc.add<bool>("LocalPositionDebug", false);
desc.add<bool>("UncalibRecHitsPlots", false);
Expand Down
66 changes: 64 additions & 2 deletions Validation/MtdValidation/plugins/EtlLocalRecoValidation.cc
View file
Open in desktop
Expand Up @@ -26,6 +26,7 @@
#include "DataFormats/ForwardDetId/interface/ETLDetId.h"
#include "DataFormats/FTLRecHit/interface/FTLRecHitCollections.h"
#include "DataFormats/FTLRecHit/interface/FTLClusterCollections.h"
#include "DataFormats/TrackerRecHit2D/interface/MTDTrackingRecHit.h"

#include "SimDataFormats/CrossingFrame/interface/CrossingFrame.h"
#include "SimDataFormats/CrossingFrame/interface/MixCollection.h"
Expand Down Expand Up @@ -58,6 +59,8 @@ class EtlLocalRecoValidation : public DQMEDAnalyzer {

void analyze(const edm::Event&, const edm::EventSetup&) override;

bool isSameCluster(const FTLCluster&, const FTLCluster&);

// ------------ member data ------------

const std::string folder_;
Expand All @@ -71,6 +74,7 @@ class EtlLocalRecoValidation : public DQMEDAnalyzer {
edm::EDGetTokenT<FTLUncalibratedRecHitCollection> etlUncalibRecHitsToken_;
edm::EDGetTokenT<CrossingFrame<PSimHit> > etlSimHitsToken_;
edm::EDGetTokenT<FTLClusterCollection> etlRecCluToken_;
edm::EDGetTokenT<MTDTrackingDetSetVector> mtdTrackingHitToken_;

edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> mtdgeoToken_;
edm::ESGetToken<MTDTopology, MTDTopologyRcd> mtdtopoToken_;
Expand Down Expand Up @@ -120,8 +124,14 @@ class EtlLocalRecoValidation : public DQMEDAnalyzer {
MonitorElement* meCluXRes_[2];
MonitorElement* meCluYRes_[2];
MonitorElement* meCluZRes_[2];
MonitorElement* meCluXPull_[2];
MonitorElement* meCluYPull_[2];
MonitorElement* meCluYXLocal_[2];
MonitorElement* meCluYXLocalSim_[2];
MonitorElement* meCluXLocalErr_[2];
MonitorElement* meCluYLocalErr_[2];

MonitorElement* meUnmatchedCluEnergy_[2];

// --- UncalibratedRecHits histograms

Expand All @@ -137,6 +147,11 @@ class EtlLocalRecoValidation : public DQMEDAnalyzer {
MonitorElement* meTimeResEta_[2][nBinsEta_];
};

bool EtlLocalRecoValidation::isSameCluster(const FTLCluster& clu1, const FTLCluster& clu2) {
return clu1.id() == clu2.id() && clu1.size() == clu2.size() && clu1.x() == clu2.x() && clu1.y() == clu2.y() &&
clu1.time() == clu2.time();
}

// ------------ constructor and destructor --------------
EtlLocalRecoValidation::EtlLocalRecoValidation(const edm::ParameterSet& iConfig)
: folder_(iConfig.getParameter<std::string>("folder")),
Expand All @@ -151,6 +166,7 @@ EtlLocalRecoValidation::EtlLocalRecoValidation(const edm::ParameterSet& iConfig)
consumes<FTLUncalibratedRecHitCollection>(iConfig.getParameter<edm::InputTag>("uncalibRecHitsTag"));
etlSimHitsToken_ = consumes<CrossingFrame<PSimHit> >(iConfig.getParameter<edm::InputTag>("simHitsTag"));
etlRecCluToken_ = consumes<FTLClusterCollection>(iConfig.getParameter<edm::InputTag>("recCluTag"));
mtdTrackingHitToken_ = consumes<MTDTrackingDetSetVector>(iConfig.getParameter<edm::InputTag>("trkHitTag"));

mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
Expand Down Expand Up @@ -182,6 +198,7 @@ void EtlLocalRecoValidation::analyze(const edm::Event& iEvent, const edm::EventS
auto etlRecHitsHandle = makeValid(iEvent.getHandle(etlRecHitsToken_));
auto etlSimHitsHandle = makeValid(iEvent.getHandle(etlSimHitsToken_));
auto etlRecCluHandle = makeValid(iEvent.getHandle(etlRecCluToken_));
auto mtdTrkHitHandle = makeValid(iEvent.getHandle(mtdTrackingHitToken_));
MixCollection<PSimHit> etlSimHits(etlSimHitsHandle.product());

// --- Loop over the ETL SIM hits
Expand Down Expand Up @@ -437,10 +454,25 @@ void EtlLocalRecoValidation::analyze(const edm::Event& iEvent, const edm::EventS

} // ihit loop

// Find the MTDTrackingRecHit corresponding to the cluster
MTDTrackingRecHit* comp(nullptr);
bool matchClu = false;
const auto& trkHits = (*mtdTrkHitHandle)[detIdObject];
for (const auto& trkHit : trkHits) {
if (isSameCluster(trkHit.mtdCluster(), cluster)) {
comp = trkHit.clone();
matchClu = true;
break;
}
}
if (!matchClu) {
edm::LogWarning("BtlLocalRecoValidation")
<< "No valid TrackingRecHit corresponding to cluster, detId = " << detIdObject.rawId();
}

// --- Fill the cluster resolution histograms
int iside = (cluId.zside() == -1 ? 0 : 1);
if (cluTimeSIM > 0. && cluEneSIM > 0.) {
int iside = (cluId.zside() == -1 ? 0 : 1);

cluTimeSIM /= cluEneSIM;

Local3DPoint cluLocalPosSIM(cluLocXSIM / cluEneSIM, cluLocYSIM / cluEneSIM, cluLocZSIM / cluEneSIM);
Expand All @@ -462,11 +494,20 @@ void EtlLocalRecoValidation::analyze(const edm::Event& iEvent, const edm::EventS
meCluTPullvsE_[iside]->Fill(cluEneSIM, time_res / cluster.timeError());

if (LocalPosDebug_) {
if (matchClu && comp != nullptr) {
meCluXPull_[iside]->Fill(x_res / std::sqrt(comp->globalPositionError().cxx()));
meCluYPull_[iside]->Fill(y_res / std::sqrt(comp->globalPositionError().cyy()));
meCluXLocalErr_[iside]->Fill(std::sqrt(comp->localPositionError().xx()));
meCluYLocalErr_[iside]->Fill(std::sqrt(comp->localPositionError().yy()));
}
meCluYXLocal_[iside]->Fill(local_point.x(), local_point.y());
meCluYXLocalSim_[iside]->Fill(cluLocalPosSIM.x(), cluLocalPosSIM.y());
}

} // if ( cluTimeSIM > 0. && cluEneSIM > 0. )
else {
meUnmatchedCluEnergy_[iside]->Fill(std::log10(cluster.energy()));
}

} // cluster loop

Expand Down Expand Up @@ -962,6 +1003,14 @@ void EtlLocalRecoValidation::bookHistograms(DQMStore::IBooker& ibook,
meCluZRes_[1] =
ibook.book1D("EtlCluZResZpos", "ETL cluster Z resolution (+Z);Z_{RECO}-Z_{SIM} [cm]", 100, -0.003, 0.003);
if (LocalPosDebug_) {
meCluXPull_[0] =
ibook.book1D("EtlCluXPullZneg", "ETL cluster X pull (-Z);X_{RECO}-X_{SIM}/sigmaX_[RECO] [cm]", 100, -5., 5.);
meCluXPull_[1] =
ibook.book1D("EtlCluXPullZpos", "ETL cluster X pull (+Z);X_{RECO}-X_{SIM}/sigmaX_[RECO] [cm]", 100, -5., 5.);
meCluYPull_[0] =
ibook.book1D("EtlCluYPullZneg", "ETL cluster Y pull (-Z);Y_{RECO}-Y_{SIM}/sigmaY_[RECO] [cm]", 100, -5., 5.);
meCluYPull_[1] =
ibook.book1D("EtlCluYPullZpos", "ETL cluster Y pull (+Z);Y_{RECO}-Y_{SIM}/sigmaY_[RECO] [cm]", 100, -5., 5.);
meCluYXLocal_[0] = ibook.book2D("EtlCluYXLocalZneg",
"ETL cluster local Y vs X (-Z);X^{local}_{RECO} [cm];Y^{local}_{RECO} [cm]",
100,
Expand Down Expand Up @@ -994,7 +1043,19 @@ void EtlLocalRecoValidation::bookHistograms(DQMStore::IBooker& ibook,
200,
-1.1,
1.1);
meCluXLocalErr_[0] =
ibook.book1D("EtlCluXLocalErrNeg", "ETL cluster X local error (-Z);sigmaX_{RECO,loc} [cm]", 50, 0., 0.2);
meCluXLocalErr_[1] =
ibook.book1D("EtlCluXLocalErrPos", "ETL cluster X local error (+Z);sigmaX_{RECO,loc} [cm]", 50, 0., 0.2);
meCluYLocalErr_[0] =
ibook.book1D("EtlCluYLocalErrNeg", "ETL cluster Y local error (-Z);sigmaY_{RECO,loc} [cm]", 50., 0., 0.2);
meCluYLocalErr_[1] =
ibook.book1D("EtlCluYLocalErrPos", "ETL cluster Y local error (+Z);sigmaY_{RECO,loc} [cm]", 50, 0., 0.2);
}
meUnmatchedCluEnergy_[0] = ibook.book1D(
"EtlUnmatchedCluEnergyNeg", "ETL unmatched cluster log10(energy) (-Z);log10(E_{RECO} [MeV])", 5, -3, 2);
meUnmatchedCluEnergy_[1] = ibook.book1D(
"EtlUnmatchedCluEnergyPos", "ETL unmatched cluster log10(energy) (+Z);log10(E_{RECO} [MeV])", 5, -3, 2);

// --- UncalibratedRecHits histograms

Expand Down Expand Up @@ -1029,6 +1090,7 @@ void EtlLocalRecoValidation::fillDescriptions(edm::ConfigurationDescriptions& de
desc.add<edm::InputTag>("uncalibRecHitsTag", edm::InputTag("mtdUncalibratedRecHits", "FTLEndcap"));
desc.add<edm::InputTag>("simHitsTag", edm::InputTag("mix", "g4SimHitsFastTimerHitsEndcap"));
desc.add<edm::InputTag>("recCluTag", edm::InputTag("mtdClusters", "FTLEndcap"));
desc.add<edm::InputTag>("trkHitTag", edm::InputTag("mtdTrackingRecHits"));
desc.add<double>("hitMinimumEnergy1Dis", 1.); // [MeV]
desc.add<double>("hitMinimumEnergy2Dis", 0.001); // [MeV]
desc.add<bool>("LocalPositionDebug", false);
Expand Down