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



}