#include "DQM/TrackerRemapper/interface/Phase1PixelROCMaps.h"
#include <bitset>
#include <cstdint>
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <numeric>
#include <sstream>
#include <string>
#include <vector>
#include "TCanvas.h"
#include "TStyle.h"

Enumerations
enum	Region { Region::Barrel = 0, Region::Forward = 1, Region::Full = 2 }

Functions
int	main (int argc, char *argv[])

Region	parseRegion (const std::string &regionStr)

std::string	regionToString (Region region)

void	showHelp (const std::string &scriptName)

Enumeration Type Documentation

◆ Region

enum Region

strong

Enumerator
Barrel
Forward
Full

Definition at line 16 of file printPixelROCsMap.cc.

                   {
   Barrel = 0,   // Assume 0 for barrel
   Forward = 1,  // 1 for forward
   Full = 2      // 2 for full
 };

Function Documentation

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Sticking ADC counts in neughbouring TS for HB:

Sticking ADC counts in neughbouring TS for HE:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Definition at line 58 of file printPixelROCsMap.cc.

References dir2webdir::argc, GCPpyPlots::argv, Barrel, svgfig::canvas(), DMR_cfg::cerr, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), gather_cfg::cout, ALCARECOPPSCalTrackBasedSel_cff::detid, Phase1PixelROCMaps::drawBarrelMaps(), Phase1PixelROCMaps::drawForwardMaps(), Phase1PixelROCMaps::drawMaps(), geometryDiff::file, MillePedeFileConverter_cfg::fileName, Phase1PixelROCMaps::fillSelectedRocs(), Phase1PixelROCMaps::fillWholeModule(), Forward, Full, mps_fire::i, makeListRunsInFiles::inputFile, mps_splice::line, parseRegion(), nano_mu_digi_cff::region, regionToString(), PixelMapPlotter::rocs, showHelp(), AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

                                  {
   static constexpr std::array<int, 3> k_height = {{1200, 600, 1600}};
 
   std::string inputFile;
   std::string inputROCsFile;
   Region region = Region::Full;  // Default value: Full region
   std::vector<std::pair<uint32_t, float>> detidValues;
   std::vector<std::pair<uint32_t, std::bitset<16>>> detidRocs;
 
   // If no arguments are passed or --h/--help is passed, show the help message
   if (argc == 1) {
     showHelp(argv[0]);
     return 0;
   }
 
   // Parse command line arguments
   for (int i = 1; i < argc; ++i) {
     std::string arg = argv[i];
 
     if (arg == "--h" || arg == "--help") {
       showHelp(argv[0]);
       return 0;  // Exit after displaying help
     } else if (arg == "--input-file" && i + 1 < argc) {
       inputFile = argv[++i];
     } else if (arg == "--input-ROCs" && i + 1 < argc) {
       inputROCsFile = argv[++i];
     } else if (arg == "--region" && i + 1 < argc) {
       std::string regionArg = argv[++i];
       try {
         region = parseRegion(regionArg);  // Parse region from string
       } catch (const std::invalid_argument&) {
         std::cerr << "Invalid value for --region: " << regionArg << "\n";
         showHelp(argv[0]);
         return 1;
       }
     } else {
       // Assume it's a DetId, convert to uint32_t
       try {
         uint32_t detid = std::stoul(arg);
         detidValues.emplace_back(detid, 1.0);  // Default value is 1.0
       } catch (const std::invalid_argument&) {
         std::cerr << "Invalid argument: " << arg << "\n";
         showHelp(argv[0]);
         return 1;
       }
     }
   }
 
   // If --input-file is provided, read from file
   if (!inputFile.empty()) {
     std::ifstream file(inputFile);
     if (!file) {
       std::cerr << "Error: Unable to open input file " << inputFile << std::endl;
       return 1;
     }
 
     std::string line;
     while (std::getline(file, line)) {
       std::istringstream iss(line);
       uint32_t detid;
       float value = 1.0;  // Default value
 
       iss >> detid;
       if (iss >> value) {  // If a second column exists, read it as value
         detidValues.emplace_back(detid, value);
       } else {
         detidValues.emplace_back(detid, 1.0);
       }
     }
   }
 
   // If --input-ROCs is provided, read from file
   if (!inputROCsFile.empty()) {
     std::ifstream file(inputROCsFile);
     if (!file) {
       std::cerr << "Error: Unable to open input ROCs file " << inputROCsFile << std::endl;
       return 1;
     }
 
     std::string line;
     while (std::getline(file, line)) {
       std::istringstream iss(line);
       uint32_t detid;
       std::string rocBits;
       iss >> detid >> rocBits;
 
       if (rocBits.length() == 16) {
         std::bitset<16> rocs(rocBits);
         detidRocs.emplace_back(detid, rocs);
       } else {
         std::cerr << "Error: Invalid ROC bits string for detid " << detid << std::endl;
         return 1;
       }
     }
   }
 
   // Create the map and fill it
   Phase1PixelROCMaps theMap("");
 
   // Fill with detidValues if --input-file or command line DetIds are used
   for (const auto& [detid, value] : detidValues) {
     theMap.fillWholeModule(detid, value);
   }
 
   // Fill with detidRocs if --input-ROCs is used
   for (const auto& [detid, rocs] : detidRocs) {
     theMap.fillSelectedRocs(detid, rocs, 1.0);  // Default value 1.0
   }
 
   // Construct the full label string
   std::string title = "Marked Pixel ROCs - " + regionToString(region);
 
   // Draw and save the map
   TCanvas canvas("Summary", "Summary", 1200, k_height[static_cast<int>(region)]);
   if (region == Region::Full) {
     theMap.drawMaps(canvas, title.c_str());
   } else if (region == Region::Barrel) {
     theMap.drawBarrelMaps(canvas, title.c_str());
   } else if (region == Region::Forward) {
     theMap.drawForwardMaps(canvas, title.c_str());
   }
 
   // Construct the filename string based on the region
   std::string fileName = "Phase1PixelROCMap_" + regionToString(region) + ".png";
   // Save the canvas with the constructed filename
   canvas.SaveAs(fileName.c_str());
 
   std::cout << "Filled Phase1 Pixel ROC map with " << detidValues.size() + detidRocs.size() << " detids." << std::endl;
 
   return 0;
 }

◆ parseRegion()

Region parseRegion ( const std::string & regionStr )

Definition at line 46 of file printPixelROCsMap.cc.

References Barrel, Forward, and Full.

Referenced by main().

                                                {
   if (regionStr == "barrel") {
     return Region::Barrel;
   } else if (regionStr == "forward") {
     return Region::Forward;
   } else if (regionStr == "full") {
     return Region::Full;
   } else {
     throw std::invalid_argument("Invalid region value");
   }
 }

◆ regionToString()

std::string regionToString ( Region region )

Definition at line 32 of file printPixelROCsMap.cc.

References Barrel, Forward, Full, and nano_mu_digi_cff::region.

Referenced by main().

                                         {
   switch (region) {
     case Region::Barrel:
       return "barrel";
     case Region::Forward:
       return "forward";
     case Region::Full:
       return "full";
     default:
       return "unknown";
   }
 }

◆ showHelp()

void showHelp ( const std::string & scriptName )

Definition at line 22 of file printPixelROCsMap.cc.

References gather_cfg::cout, and Herwig7_Dummy_ReadLHE_GEN_SIM::scriptName.

Referenced by main().

                                            {
   std::cout << "Usage: " << scriptName << " [options] <detid>\n"
             << "  --input-file <filename>        Specify the input file\n"
             << "  --input-ROCs <filename>        Specify the input ROCs file\n"
             << "  --region <barrel|forward|full> Specify the region (default: full)\n"
             << "  --h or --help                  Show this help message\n"
             << "  <detid>                        Provide DetId (list of DetIds)\n";
 }

Enumerations

Functions

Enumeration Type Documentation

◆ Region

Function Documentation

◆ main()

◆ parseRegion()

◆ regionToString()

◆ showHelp()