Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2022-05-26 22:39:55

0001 // -*- C++ -*-
0002 //
0003 // Package:    Validation/MtdValidation
0004 // Class:      BtlDigiHitsValidation
0005 //
0006 /**\class BtlDigiHitsValidation BtlDigiHitsValidation.cc Validation/MtdValidation/plugins/BtlDigiHitsValidation.cc
0007 
0008  Description: BTL DIGI hits validation
0009 
0010  Implementation:
0011      [Notes on implementation]
0012 */
0013 
0014 #include <string>
0015 
0016 #include "FWCore/Framework/interface/Frameworkfwd.h"
0017 #include "FWCore/Framework/interface/Event.h"
0018 #include "FWCore/Framework/interface/MakerMacros.h"
0019 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0020 
0021 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
0022 #include "DQMServices/Core/interface/DQMStore.h"
0023 
0024 #include "DataFormats/Common/interface/ValidHandle.h"
0025 #include "DataFormats/ForwardDetId/interface/BTLDetId.h"
0026 #include "DataFormats/FTLDigi/interface/FTLDigiCollections.h"
0027 
0028 #include "Geometry/Records/interface/MTDDigiGeometryRecord.h"
0029 #include "Geometry/Records/interface/MTDTopologyRcd.h"
0030 #include "Geometry/MTDGeometryBuilder/interface/MTDGeometry.h"
0031 #include "Geometry/MTDNumberingBuilder/interface/MTDTopology.h"
0032 
0033 #include "Geometry/MTDGeometryBuilder/interface/ProxyMTDTopology.h"
0034 #include "Geometry/MTDGeometryBuilder/interface/RectangularMTDTopology.h"
0035 
0036 #include "Geometry/MTDCommonData/interface/MTDTopologyMode.h"
0037 
0038 class BtlDigiHitsValidation : public DQMEDAnalyzer {
0039 public:
0040   explicit BtlDigiHitsValidation(const edm::ParameterSet&);
0041   ~BtlDigiHitsValidation() override;
0042 
0043   static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
0044 
0045 private:
0046   void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
0047 
0048   void analyze(const edm::Event&, const edm::EventSetup&) override;
0049 
0050   // ------------ member data ------------
0051 
0052   const std::string folder_;
0053   const bool optionalPlots_;
0054 
0055   edm::EDGetTokenT<BTLDigiCollection> btlDigiHitsToken_;
0056 
0057   edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> mtdgeoToken_;
0058   edm::ESGetToken<MTDTopology, MTDTopologyRcd> mtdtopoToken_;
0059 
0060   // --- histograms declaration
0061 
0062   MonitorElement* meNhits_[2];
0063 
0064   MonitorElement* meHitCharge_[2];
0065   MonitorElement* meHitTime_[2];
0066 
0067   MonitorElement* meOccupancy_[2];
0068 
0069   //local position monitoring
0070   MonitorElement* meLocalOccupancy_[2];
0071   MonitorElement* meHitXlocal_[2];
0072   MonitorElement* meHitYlocal_[2];
0073   MonitorElement* meHitZlocal_[2];
0074 
0075   MonitorElement* meHitX_[2];
0076   MonitorElement* meHitY_[2];
0077   MonitorElement* meHitZ_[2];
0078   MonitorElement* meHitPhi_[2];
0079   MonitorElement* meHitEta_[2];
0080 
0081   MonitorElement* meHitTvsQ_[2];
0082   MonitorElement* meHitQvsPhi_[2];
0083   MonitorElement* meHitQvsEta_[2];
0084   MonitorElement* meHitQvsZ_[2];
0085   MonitorElement* meHitTvsPhi_[2];
0086   MonitorElement* meHitTvsEta_[2];
0087   MonitorElement* meHitTvsZ_[2];
0088 };
0089 
0090 // ------------ constructor and destructor --------------
0091 BtlDigiHitsValidation::BtlDigiHitsValidation(const edm::ParameterSet& iConfig)
0092     : folder_(iConfig.getParameter<std::string>("folder")),
0093       optionalPlots_(iConfig.getParameter<bool>("optionalPlots")) {
0094   btlDigiHitsToken_ = consumes<BTLDigiCollection>(iConfig.getParameter<edm::InputTag>("inputTag"));
0095   mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
0096   mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
0097 }
0098 
0099 BtlDigiHitsValidation::~BtlDigiHitsValidation() {}
0100 
0101 // ------------ method called for each event  ------------
0102 void BtlDigiHitsValidation::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
0103   using namespace edm;
0104 
0105   auto geometryHandle = iSetup.getTransientHandle(mtdgeoToken_);
0106   const MTDGeometry* geom = geometryHandle.product();
0107 
0108   auto topologyHandle = iSetup.getTransientHandle(mtdtopoToken_);
0109   const MTDTopology* topology = topologyHandle.product();
0110 
0111   auto btlDigiHitsHandle = makeValid(iEvent.getHandle(btlDigiHitsToken_));
0112 
0113   // --- Loop over the BLT DIGI hits
0114 
0115   unsigned int n_digi_btl[2] = {0, 0};
0116 
0117   for (const auto& dataFrame : *btlDigiHitsHandle) {
0118     BTLDetId detId = dataFrame.id();
0119     DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology->getMTDTopologyMode()));
0120     const MTDGeomDet* thedet = geom->idToDet(geoId);
0121     if (thedet == nullptr)
0122       throw cms::Exception("BtlDigiHitsValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
0123                                                     << detId.rawId() << ") is invalid!" << std::dec << std::endl;
0124     const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
0125     const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
0126 
0127     Local3DPoint local_point(0., 0., 0.);
0128     local_point = topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
0129     const auto& global_point = thedet->toGlobal(local_point);
0130 
0131     const auto& sample_L = dataFrame.sample(0);
0132     const auto& sample_R = dataFrame.sample(1);
0133 
0134     uint32_t adc[2] = {sample_L.data(), sample_R.data()};
0135     uint32_t tdc[2] = {sample_L.toa(), sample_R.toa()};
0136 
0137     for (int iside = 0; iside < 2; ++iside) {
0138       if (adc[iside] == 0)
0139         continue;
0140 
0141       meHitCharge_[iside]->Fill(adc[iside]);
0142       meHitTime_[iside]->Fill(tdc[iside]);
0143 
0144       meOccupancy_[iside]->Fill(global_point.z(), global_point.phi());
0145 
0146       if (optionalPlots_) {
0147         meLocalOccupancy_[iside]->Fill(local_point.x(), local_point.y());
0148         meHitXlocal_[iside]->Fill(local_point.x());
0149         meHitYlocal_[iside]->Fill(local_point.y());
0150         meHitZlocal_[iside]->Fill(local_point.z());
0151       }
0152 
0153       meHitX_[iside]->Fill(global_point.x());
0154       meHitY_[iside]->Fill(global_point.y());
0155       meHitZ_[iside]->Fill(global_point.z());
0156       meHitPhi_[iside]->Fill(global_point.phi());
0157       meHitEta_[iside]->Fill(global_point.eta());
0158 
0159       meHitTvsQ_[iside]->Fill(adc[iside], tdc[iside]);
0160       meHitQvsPhi_[iside]->Fill(global_point.phi(), adc[iside]);
0161       meHitQvsEta_[iside]->Fill(global_point.eta(), adc[iside]);
0162       meHitQvsZ_[iside]->Fill(global_point.z(), adc[iside]);
0163       meHitTvsPhi_[iside]->Fill(global_point.phi(), tdc[iside]);
0164       meHitTvsEta_[iside]->Fill(global_point.eta(), tdc[iside]);
0165       meHitTvsZ_[iside]->Fill(global_point.z(), tdc[iside]);
0166 
0167       n_digi_btl[iside]++;
0168 
0169     }  // iside loop
0170 
0171   }  // dataFrame loop
0172 
0173   if (n_digi_btl[0] > 0)
0174     meNhits_[0]->Fill(log10(n_digi_btl[0]));
0175   if (n_digi_btl[1] > 0)
0176     meNhits_[1]->Fill(log10(n_digi_btl[1]));
0177 }
0178 
0179 // ------------ method for histogram booking ------------
0180 void BtlDigiHitsValidation::bookHistograms(DQMStore::IBooker& ibook,
0181                                            edm::Run const& run,
0182                                            edm::EventSetup const& iSetup) {
0183   ibook.setCurrentFolder(folder_);
0184   // --- histograms booking
0185 
0186   meNhits_[0] = ibook.book1D("BtlNhitsL", "Number of BTL DIGI hits (L);log_{10}(N_{DIGI})", 100, 0., 5.25);
0187   meNhits_[1] = ibook.book1D("BtlNhitsR", "Number of BTL DIGI hits (R);log_{10}(N_{DIGI})", 100, 0., 5.25);
0188 
0189   meHitCharge_[0] = ibook.book1D("BtlHitChargeL", "BTL DIGI hits charge (L);Q_{DIGI} [ADC counts]", 100, 0., 1024.);
0190   meHitCharge_[1] = ibook.book1D("BtlHitChargeR", "BTL DIGI hits charge (R);Q_{DIGI} [ADC counts]", 100, 0., 1024.);
0191   meHitTime_[0] = ibook.book1D("BtlHitTimeL", "BTL DIGI hits ToA (L);ToA_{DIGI} [TDC counts]", 100, 0., 1024.);
0192   meHitTime_[1] = ibook.book1D("BtlHitTimeR", "BTL DIGI hits ToA (R);ToA_{DIGI} [TDC counts]", 100, 0., 1024.);
0193   meOccupancy_[0] = ibook.book2D("BtlOccupancyL",
0194                                  "BTL DIGI hits occupancy (L);Z_{DIGI} [cm]; #phi_{DIGI} [rad]",
0195                                  65,
0196                                  -260.,
0197                                  260.,
0198                                  126,
0199                                  -3.15,
0200                                  3.15);
0201   meOccupancy_[1] = ibook.book2D("BtlOccupancyR",
0202                                  "BTL DIGI hits occupancy (R);Z_{DIGI} [cm]; #phi_{DIGI} [rad]",
0203                                  65,
0204                                  -260.,
0205                                  260.,
0206                                  126,
0207                                  -3.15,
0208                                  3.15);
0209   if (optionalPlots_) {
0210     meLocalOccupancy_[0] = ibook.book2D("BtlLocalOccupancyL",
0211                                         "BTL DIGI hits local occupancy (L);X_{DIGI} [cm]; Y_{DIGI} [cm]",
0212                                         100,
0213                                         -10.,
0214                                         10,
0215                                         60,
0216                                         -3.,
0217                                         3.);
0218     meLocalOccupancy_[1] = ibook.book2D(
0219         "BtlLocalOccupancyR", "BTL DIGI hits occupancy (R);X_{DIGI} [cm]; Y_{DIGI} [cm]", 100, -10., 10., 60, -3., 3.);
0220     meHitXlocal_[0] = ibook.book1D("BtlHitXlocalL", "BTL DIGI local X (L);X_{DIGI}^{LOC} [cm]", 100, -10., 10.);
0221     meHitXlocal_[1] = ibook.book1D("BtlHitXlocalR", "BTL DIGI local X (R);X_{DIGI}^{LOC} [cm]", 100, -10., 10.);
0222     meHitYlocal_[0] = ibook.book1D("BtlHitYlocalL", "BTL DIGI local Y (L);Y_{DIGI}^{LOC} [cm]", 60, -3., 3.);
0223     meHitYlocal_[1] = ibook.book1D("BtlHitYlocalR", "BTL DIGI local Y (R);Y_{DIGI}^{LOC} [cm]", 60, -3., 3.);
0224     meHitZlocal_[0] = ibook.book1D("BtlHitZlocalL", "BTL DIGI local z (L);z_{DIGI}^{LOC} [cm]", 10, -1, 1);
0225     meHitZlocal_[1] = ibook.book1D("BtlHitZlocalR", "BTL DIGI local z (R);z_{DIGI}^{LOC} [cm]", 10, -1, 1);
0226   }
0227 
0228   meHitX_[0] = ibook.book1D("BtlHitXL", "BTL DIGI hits X (L);X_{DIGI} [cm]", 60, -120., 120.);
0229   meHitX_[1] = ibook.book1D("BtlHitXR", "BTL DIGI hits X (R);X_{DIGI} [cm]", 60, -120., 120.);
0230   meHitY_[0] = ibook.book1D("BtlHitYL", "BTL DIGI hits Y (L);Y_{DIGI} [cm]", 60, -120., 120.);
0231   meHitY_[1] = ibook.book1D("BtlHitYR", "BTL DIGI hits Y (R);Y_{DIGI} [cm]", 60, -120., 120.);
0232   meHitZ_[0] = ibook.book1D("BtlHitZL", "BTL DIGI hits Z (L);Z_{DIGI} [cm]", 100, -260., 260.);
0233   meHitZ_[1] = ibook.book1D("BtlHitZR", "BTL DIGI hits Z (R);Z_{DIGI} [cm]", 100, -260., 260.);
0234   meHitPhi_[0] = ibook.book1D("BtlHitPhiL", "BTL DIGI hits #phi (L);#phi_{DIGI} [rad]", 126, -3.15, 3.15);
0235   meHitPhi_[1] = ibook.book1D("BtlHitPhiR", "BTL DIGI hits #phi (R);#phi_{DIGI} [rad]", 126, -3.15, 3.15);
0236   meHitEta_[0] = ibook.book1D("BtlHitEtaL", "BTL DIGI hits #eta (L);#eta_{DIGI}", 100, -1.55, 1.55);
0237   meHitEta_[1] = ibook.book1D("BtlHitEtaR", "BTL DIGI hits #eta (R);#eta_{DIGI}", 100, -1.55, 1.55);
0238 
0239   meHitTvsQ_[0] = ibook.bookProfile("BtlHitTvsQL",
0240                                     "BTL DIGI ToA vs charge (L);Q_{DIGI} [ADC counts];ToA_{DIGI} [TDC counts]",
0241                                     50,
0242                                     0.,
0243                                     1024.,
0244                                     0.,
0245                                     1024.);
0246   meHitTvsQ_[1] = ibook.bookProfile("BtlHitTvsQR",
0247                                     "BTL DIGI ToA vs charge (R);Q_{DIGI} [ADC counts];ToA_{DIGI} [TDC counts]",
0248                                     50,
0249                                     0.,
0250                                     1024.,
0251                                     0.,
0252                                     1024.);
0253   meHitQvsPhi_[0] = ibook.bookProfile("BtlHitQvsPhiL",
0254                                       "BTL DIGI charge vs #phi (L);#phi_{DIGI} [rad];Q_{DIGI} [ADC counts]",
0255                                       50,
0256                                       -3.15,
0257                                       3.15,
0258                                       0.,
0259                                       1024.);
0260   meHitQvsPhi_[1] = ibook.bookProfile("BtlHitQvsPhiR",
0261                                       "BTL DIGI charge vs #phi (R);#phi_{DIGI} [rad];Q_{DIGI} [ADC counts]",
0262                                       50,
0263                                       -3.15,
0264                                       3.15,
0265                                       0.,
0266                                       1024.);
0267   meHitQvsEta_[0] = ibook.bookProfile(
0268       "BtlHitQvsEtaL", "BTL DIGI charge vs #eta (L);#eta_{DIGI};Q_{DIGI} [ADC counts]", 50, -1.55, 1.55, 0., 1024.);
0269   meHitQvsEta_[1] = ibook.bookProfile(
0270       "BtlHitQvsEtaR", "BTL DIGI charge vs #eta (R);#eta_{DIGI};Q_{DIGI} [ADC counts]", 50, -1.55, 1.55, 0., 1024.);
0271   meHitQvsZ_[0] = ibook.bookProfile(
0272       "BtlHitQvsZL", "BTL DIGI charge vs Z (L);Z_{DIGI} [cm];Q_{DIGI} [ADC counts]", 50, -260., 260., 0., 1024.);
0273   meHitQvsZ_[1] = ibook.bookProfile(
0274       "BtlHitQvsZR", "BTL DIGI charge vs Z (R);Z_{DIGI} [cm];Q_{DIGI} [ADC counts]", 50, -260., 260., 0., 1024.);
0275   meHitTvsPhi_[0] = ibook.bookProfile(
0276       "BtlHitTvsPhiL", "BTL DIGI ToA vs #phi (L);#phi_{DIGI} [rad];ToA_{DIGI} [TDC counts]", 50, -3.15, 3.15, 0., 1024.);
0277   meHitTvsPhi_[1] = ibook.bookProfile(
0278       "BtlHitTvsPhiR", "BTL DIGI ToA vs #phi (R);#phi_{DIGI} [rad];ToA_{DIGI} [TDC counts]", 50, -3.15, 3.15, 0., 1024.);
0279   meHitTvsEta_[0] = ibook.bookProfile(
0280       "BtlHitTvsEtaL", "BTL DIGI ToA vs #eta (L);#eta_{DIGI};ToA_{DIGI} [TDC counts]", 50, -1.55, 1.55, 0., 1024.);
0281   meHitTvsEta_[1] = ibook.bookProfile(
0282       "BtlHitTvsEtaR", "BTL DIGI ToA vs #eta (R);#eta_{DIGI};ToA_{DIGI} [TDC counts]", 50, -1.55, 1.55, 0., 1024.);
0283   meHitTvsZ_[0] = ibook.bookProfile(
0284       "BtlHitTvsZL", "BTL DIGI ToA vs Z (L);Z_{DIGI} [cm];ToA_{DIGI} [TDC counts]", 50, -260., 260., 0., 1024.);
0285   meHitTvsZ_[1] = ibook.bookProfile(
0286       "BtlHitTvsZR", "BTL DIGI ToA vs Z (R);Z_{DIGI} [cm];ToA_{DIGI} [TDC counts]", 50, -260., 260., 0., 1024.);
0287 }
0288 
0289 // ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
0290 void BtlDigiHitsValidation::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
0291   edm::ParameterSetDescription desc;
0292 
0293   desc.add<std::string>("folder", "MTD/BTL/DigiHits");
0294   desc.add<edm::InputTag>("inputTag", edm::InputTag("mix", "FTLBarrel"));
0295   desc.add<bool>("optionalPlots", false);
0296 
0297   descriptions.add("btlDigiHitsDefaultValid", desc);
0298 }
0299 
0300 DEFINE_FWK_MODULE(BtlDigiHitsValidation);