Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2022-03-05 02:38:57

0001 
0002 /*----------------------------------------------------------------------
0003 
0004 Toy EDAnalyzer for testing purposes only.
0005 
0006 ----------------------------------------------------------------------*/
0007 
0008 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
0009 #include "FWCore/Framework/interface/ESHandle.h"
0010 #include "FWCore/Framework/interface/Event.h"
0011 #include "FWCore/Framework/interface/MakerMacros.h"
0012 #include <cmath>
0013 #include <fstream>
0014 #include <iostream>
0015 #include <map>
0016 #include <sstream>
0017 #include <stdexcept>
0018 #include <string>
0019 
0020 #include "FWCore/Framework/interface/EventSetup.h"
0021 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0022 
0023 #include "CondTools/DT/test/validate/DTROMapValidateDBRead.h"
0024 #include "CondFormats/DTObjects/interface/DTReadOutMapping.h"
0025 #include "CondFormats/DataRecord/interface/DTReadOutMappingRcd.h"
0026 
0027 DTROMapValidateDBRead::DTROMapValidateDBRead(edm::ParameterSet const& p)
0028     : dataFileName(p.getParameter<std::string>("chkFile")),
0029       elogFileName(p.getParameter<std::string>("logFile")),
0030       dtreadoutmappingToken_(esConsumes()) {}
0031 
0032 DTROMapValidateDBRead::DTROMapValidateDBRead(int i) : dtreadoutmappingToken_(esConsumes()) {}
0033 
0034 void DTROMapValidateDBRead::analyze(const edm::Event& e, const edm::EventSetup& context) {
0035   using namespace edm::eventsetup;
0036   // Context is not used.
0037   std::cout << " I AM IN RUN NUMBER " << e.id().run() << std::endl;
0038   std::cout << " ---EVENT NUMBER " << e.id().event() << std::endl;
0039   std::stringstream run_fn;
0040   //  run_fn << "run" << e.id().run() << dataFileName;
0041   //  std::ifstream chkFile( run_fn.str().c_str() );
0042   std::ifstream chkFile(dataFileName.c_str());
0043   std::ofstream logFile(elogFileName.c_str(), std::ios_base::app);
0044   auto ro = context.getHandle(dtreadoutmappingToken_);
0045   std::cout << ro->mapRobRos() << " " << ro->mapCellTdc() << std::endl;
0046   std::cout << std::distance(ro->begin(), ro->end()) << " data in the container" << std::endl;
0047   int whe;
0048   int sta;
0049   int sec;
0050   int qua;
0051   int lay;
0052   int cel;
0053   int ddu;
0054   int ros;
0055   int rob;
0056   int tdc;
0057   int cha;
0058   int ckwhe;
0059   int cksta;
0060   int cksec;
0061   int ckqua;
0062   int cklay;
0063   int ckcel;
0064   int ckddu;
0065   int ckros;
0066   int ckrob;
0067   int cktdc;
0068   int ckcha;
0069   int status;
0070   DTReadOutMapping::const_iterator iter = ro->begin();
0071   DTReadOutMapping::const_iterator iend = ro->end();
0072   while (iter != iend) {
0073     const DTReadOutGeometryLink& link = *iter;
0074     status = ro->readOutToGeometry(
0075         link.dduId, link.rosId, link.robId, link.tdcId, link.channelId, whe, sta, sec, qua, lay, cel);
0076     if (status)
0077       logFile << "ERROR while getting chan->cell map " << link.dduId << " " << link.rosId << " " << link.robId << " "
0078               << link.tdcId << " " << link.channelId << " , status = " << status << std::endl;
0079     if ((link.wheelId != whe) || (link.stationId != sta) || (link.sectorId != sec) || (link.slId != qua) ||
0080         (link.layerId != lay) || (link.cellId != cel))
0081       logFile << "MISMATCH WHEN READING chan->cell  " << link.dduId << " " << link.rosId << " " << link.robId << " "
0082               << link.tdcId << " " << link.channelId << " : " << link.wheelId << " " << link.stationId << " "
0083               << link.sectorId << " " << link.slId << " " << link.layerId << " " << link.cellId << " -> " << whe << " "
0084               << sta << " " << sec << " " << qua << " " << lay << " " << cel << std::endl;
0085     status = ro->geometryToReadOut(
0086         link.wheelId, link.stationId, link.sectorId, link.slId, link.layerId, link.cellId, ddu, ros, rob, tdc, cha);
0087     if (status)
0088       logFile << "ERROR while getting cell->chan map " << link.wheelId << " " << link.stationId << " " << link.sectorId
0089               << " " << link.slId << " " << link.layerId << " " << link.cellId << " , status = " << status << std::endl;
0090     if ((link.dduId != ddu) || (link.rosId != ros) || (link.robId != rob) || (link.tdcId != tdc) ||
0091         (link.channelId != cha))
0092       logFile << "MISMATCH WHEN READING cell->chan  " << link.wheelId << " " << link.stationId << " " << link.sectorId
0093               << " " << link.slId << " " << link.layerId << " " << link.cellId << " : " << link.dduId << " "
0094               << link.rosId << " " << link.robId << " " << link.tdcId << " " << link.channelId << " -> " << ddu << " "
0095               << ros << " " << rob << " " << tdc << " " << cha << std::endl;
0096     iter++;
0097   }
0098 
0099   while (chkFile >> ckddu >> ckros >> ckrob >> cktdc >> ckcha >> ckwhe >> cksta >> cksec >> ckqua >> cklay >> ckcel) {
0100     status = ro->readOutToGeometry(ckddu, ckros, ckrob, cktdc, ckcha, whe, sta, sec, qua, lay, cel);
0101     if ((ckwhe != whe) || (cksta != sta) || (cksec != sec) || (ckqua != qua) || (cklay != lay) || (ckcel != cel))
0102       logFile << "MISMATCH IN WRITING AND READING chan->cell map " << ckddu << " " << ckros << " " << ckrob << " "
0103               << cktdc << " " << ckcha << " : " << ckwhe << " " << cksta << " " << cksec << " " << ckqua << " " << cklay
0104               << " " << ckcel << " -> " << whe << " " << sta << " " << sec << " " << qua << " " << lay << " " << cel
0105               << std::endl;
0106     status = ro->geometryToReadOut(ckwhe, cksta, cksec, ckqua, cklay, ckcel, ddu, ros, rob, tdc, cha);
0107     if ((ckddu != ddu) || (ckros != ros) || (ckrob != rob) || (cktdc != tdc) || (ckcha != cha))
0108       logFile << "MISMATCH IN WRITING AND READING cell->chan map " << ckwhe << " " << cksta << " " << cksec << " "
0109               << ckqua << " " << cklay << " " << ckcel << " : " << ckddu << " " << ckros << " " << ckrob << " " << cktdc
0110               << " " << ckcha << " -> " << ddu << " " << ros << " " << rob << " " << tdc << " " << cha << std::endl;
0111   }
0112 }
0113 
0114 void DTROMapValidateDBRead::endJob() {
0115   std::ifstream logFile(elogFileName.c_str());
0116   char* line = new char[1000];
0117   int errors = 0;
0118   std::cout << "ReadOut Map validation result:" << std::endl;
0119   while (logFile.getline(line, 1000)) {
0120     std::cout << line << std::endl;
0121     errors++;
0122   }
0123   if (!errors) {
0124     std::cout << " ********************************* " << std::endl;
0125     std::cout << " ***                           *** " << std::endl;
0126     std::cout << " ***      NO ERRORS FOUND      *** " << std::endl;
0127     std::cout << " ***                           *** " << std::endl;
0128     std::cout << " ********************************* " << std::endl;
0129   }
0130   return;
0131 }
0132 
0133 DEFINE_FWK_MODULE(DTROMapValidateDBRead);