|
44 | 44 | #include "FWCore/ParameterSet/interface/ParameterSet.h" |
45 | 45 | #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h" |
46 | 46 | #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h" |
47 | | -#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h" |
48 | 47 | #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h" |
49 | 48 | #include "MagneticField/Engine/interface/MagneticField.h" |
50 | 49 | #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h" |
@@ -248,13 +247,15 @@ void SiStripLorentzAnglePCLMonitor::analyze(const edm::Event& iEvent, const edm: |
248 | 247 | iHists_.h1_["track_eta"]->Fill(track->eta()); |
249 | 248 | iHists_.h1_["track_phi"]->Fill(track->phi()); |
250 | 249 | iHists_.h1_["track_validhits"]->Fill(track->numberOfValidHits()); |
251 | | - iHists_.h1_["track_chi2ndof"]->Fill((track->chi2() / track->ndof())); |
252 | | - iHists_.h2_["track_chi2xhits"]->Fill((track->chi2() / track->ndof()), track->numberOfValidHits()); |
| 250 | + |
| 251 | + const auto normChi2 = track->ndof() > 0 ? track->chi2() / track->ndof() : -1.; |
| 252 | + iHists_.h1_["track_chi2ndof"]->Fill(normChi2); |
| 253 | + iHists_.h2_["track_chi2xhits"]->Fill(normChi2, track->numberOfValidHits()); |
253 | 254 | iHists_.h2_["track_ptxhits"]->Fill(track->pt(), track->numberOfValidHits()); |
254 | 255 | iHists_.h2_["track_etaxhits"]->Fill(track->eta(), track->numberOfValidHits()); |
255 | | - iHists_.h2_["track_ptxchi2"]->Fill(track->pt(), (track->chi2() / track->ndof())); |
| 256 | + iHists_.h2_["track_ptxchi2"]->Fill(track->pt(), normChi2); |
256 | 257 | iHists_.h2_["track_ptxeta"]->Fill(track->pt(), track->eta()); |
257 | | - iHists_.h2_["track_etaxchi2"]->Fill(track->eta(), (track->chi2() / track->ndof())); |
| 258 | + iHists_.h2_["track_etaxchi2"]->Fill(track->eta(), normChi2); |
258 | 259 |
|
259 | 260 | edm::LogInfo("SiStripLorentzAnglePCLMonitor") |
260 | 261 | << " track pT()" << track->pt() << " track eta()" << track->eta() << std::endl; |
|
0 commit comments