Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:07:59

0001 #include "CommonTools/UtilAlgos/interface/TFileService.h"
0002 #include "DataFormats/EcalDetId/interface/ESDetId.h"
0003 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
0004 #include "FWCore/Framework/interface/EventSetup.h"
0005 #include "FWCore/Framework/interface/MakerMacros.h"
0006 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0007 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
0008 #include "FWCore/ServiceRegistry/interface/Service.h"
0009 #include "Geometry/Records/interface/CaloGeometryRecord.h"
0010 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
0011 #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
0012 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
0013 #include "Geometry/CaloGeometry/interface/EZArrayFL.h"
0014 
0015 #include "TH2.h"
0016 
0017 #include <iomanip>
0018 #include <iostream>
0019 #include <sstream>  // for ostringstream
0020 
0021 typedef EZArrayFL<GlobalPoint> CornersVec;
0022 
0023 class EcalPreshowerCellParameterDump : public edm::one::EDAnalyzer<edm::one::SharedResources> {
0024 public:
0025   explicit EcalPreshowerCellParameterDump(const edm::ParameterSet&);
0026   ~EcalPreshowerCellParameterDump() override {}
0027 
0028   static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
0029 
0030   void beginJob() override {}
0031   void analyze(edm::Event const& iEvent, edm::EventSetup const&) override;
0032   void endJob() override {}
0033 
0034 private:
0035   const bool debug_;
0036   const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> tok_geom_;
0037   std::vector<TH2D*> hist_;
0038 };
0039 
0040 EcalPreshowerCellParameterDump::EcalPreshowerCellParameterDump(const edm::ParameterSet& ps)
0041     : debug_(ps.getUntrackedParameter<bool>("debug", false)),
0042       tok_geom_(esConsumes<CaloGeometry, CaloGeometryRecord>()) {
0043   usesResource(TFileService::kSharedResource);
0044 
0045   if (debug_) {
0046     edm::Service<TFileService> fs;
0047     for (short iz = 0; iz < 2; ++iz) {
0048       short zside = 2 * iz - 1;
0049       for (short lay = 1; lay <= 2; ++lay) {
0050         std::ostringstream name, title;
0051         name << "hist" << iz << lay;
0052         title << "y vs. x (zside = " << zside << ",layer = " << lay << ")";
0053         hist_.emplace_back(
0054             fs->make<TH2D>(name.str().c_str(), title.str().c_str(), 5000, -125.0, 125.0, 5000, -125.0, 125.0));
0055       }
0056     }
0057   }
0058 }
0059 
0060 void EcalPreshowerCellParameterDump::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
0061   edm::ParameterSetDescription desc;
0062   desc.addUntracked<bool>("debug", false);
0063   descriptions.add("ecalPreshowerCellParameterDump", desc);
0064 }
0065 
0066 void EcalPreshowerCellParameterDump::analyze(const edm::Event& /*iEvent*/, const edm::EventSetup& iSetup) {
0067   const CaloGeometry* geo = &iSetup.getData(tok_geom_);
0068   const CaloSubdetectorGeometry* ecalGeom =
0069       static_cast<const CaloSubdetectorGeometry*>(geo->getSubdetectorGeometry(DetId::Ecal, EcalPreshower));
0070 
0071   std::cout << "\n\nStudy Detector = Ecal SubDetector = ES"
0072             << "\n======================================\n\n";
0073   const std::vector<DetId>& ids = ecalGeom->getValidDetIds(DetId::Ecal, EcalPreshower);
0074   int nall(0);
0075   for (auto id : ids) {
0076     ++nall;
0077     auto geom = ecalGeom->getGeometry(id);
0078     ESDetId esid(id);
0079 
0080     if (debug_) {
0081       std::cout << nall << " " << esid.rawId() << " " << std::setprecision(6) << geom->getPosition() << std::endl;
0082       unsigned int hid = ((esid.zside() + 1) + esid.plane() - 1);
0083       if (hid < hist_.size())
0084         hist_[hid]->Fill(geom->getPosition().x(), geom->getPosition().y());
0085     } else {
0086       std::cout << "Cell[" << nall << "] " << esid << " geom->getPosition " << std::setprecision(4)
0087                 << geom->getPosition() << " BackPoint " << geom->getBackPoint() << " [rho,eta:etaSpan,phi:phiSpan] ("
0088                 << geom->rhoPos() << ", " << geom->etaPos() << ":" << geom->etaSpan() << ", " << geom->phiPos() << ":"
0089                 << geom->phiSpan() << ")";
0090 
0091       const CaloCellGeometry::CornersVec& corners(geom->getCorners());
0092 
0093       for (unsigned int ci(0); ci != corners.size(); ci++) {
0094         std::cout << " Corner: " << ci << "  Location" << corners[ci] << " ; ";
0095       }
0096 
0097       std::cout << std::endl;
0098     }
0099   }
0100   std::cout << "\n\nDumps a total of : " << nall << " cells of the detector\n" << std::endl;
0101 }
0102 
0103 DEFINE_FWK_MODULE(EcalPreshowerCellParameterDump);