print_TrackerMap.cc

Go to the documentation of this file.

#include <iostream>
#include <string>
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CommonTools/TrackerMap/interface/TrackerMap.h"
 
void printTrackerMap(const std::string filename,
                     const std::string title,
                     const std::string outfile,
                     const int size,
                     const std::string logscale,
                     std::string withpixel,
                     const float min,
                     const float max);
 
int main(int argc, char* argv[]) {
  if (argc >= 4) {
    int size = 4500;
    float min = 10.;
    float max = -1.;
    std::string logscale = "False";
    std::string withpixel = "False";
 
    char* filename = argv[1];
    char* title = argv[2];
    char* outfile = argv[3];
    std::cout << "ready to use file " << filename << " to create the map " << title << " and save it in " << outfile
              << std::endl;
 
    if (argc >= 5) {
      size = atoi(argv[4]);
    }
    if (argc >= 6) {
      logscale = argv[5];
      withpixel = argv[6];
    }
    if (argc >= 9) {
      min = atof(argv[7]);
      max = atof(argv[8]);
    }
 
    printTrackerMap(filename, title, outfile, size, logscale, withpixel, min, max);
 
  } else {
    std::cout << "Wrong number of arguments: " << argc << std::endl;
    return -1;
  }
 
  return 0;
}
 
void printTrackerMap(const std::string filename,
                     const std::string title,
                     const std::string outfile,
                     const int size,
                     const std::string logscale,
                     std::string withpixel,
                     const float min,
                     const float max) {
  edm::ParameterSet pset;
 
  if (logscale == "True")
    pset.addUntrackedParameter<bool>("logScale", true);
  if (logscale == "False")
    pset.addUntrackedParameter<bool>("logScale", false);
 
  TrackerMap themap(pset);
  themap.setTitle(title);  // title as input parameter
  double ratio = 2400. / 4500.;
  if (withpixel == "True") {
    themap.addPixel(true);
    themap.onlyPixel(false);
    ratio = 8. / 19.;
  }
  if (withpixel == "Only") {
    themap.addPixel(false);
    themap.onlyPixel(true);
    ratio = 16. / 9.;
  }
  if (withpixel == "False") {
    themap.addPixel(false);
    themap.onlyPixel(false);
    ratio = 8. / 15.;
  }
 
  std::ifstream input(filename);
 
  unsigned int detid;
  float val;
 
  while (input >> detid >> val) {
    themap.fill_current_val(detid, val);
  }
 
  std::cout << "preparing a map with size " << size << "x" << int(size * ratio) << std::endl;
 
  themap.save(true, min, max, outfile, size, int(size * ratio));  // name and size as input parameters
}