Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
#include <iostream>
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "CondTools/Ecal/interface/EcalIntercalibConstantsXMLTranslator.h"
#include "CondTools/Ecal/interface/EcalCondHeader.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"

#include <string>
#include <sstream>
#include <fstream>

using namespace std;

void usage() {
  cout << endl;
  cout << "array2xml [arrayfile] [xmlfile]" << endl;
  cout << "Read coefficients from straight array [denseindex]" << endl;
  cout << "and write in xml format" << endl;
}

int main(int argc, char* argv[]) {
  if (argc != 3) {
    usage();
    exit(0);
  }

  string arrayfilename(argv[1]);
  string xmlfilename(argv[2]);
  fstream arrayfile(arrayfilename.c_str(), ios::in);

  EcalIntercalibConstants rcd;

  float c = 0;
  int idx = 0;

  while (arrayfile >> c) {
    uint32_t id = EBDetId::unhashIndex(idx);
    rcd[id] = c;
    ++idx;
  }
  cout << idx << endl;

  for (int cellid = 0; cellid < EEDetId::kSizeForDenseIndexing; ++cellid) {  // loop on EB cells

    if (EEDetId::validHashIndex(cellid)) {
      uint32_t rawid = EEDetId::unhashIndex(cellid);

      rcd[rawid] = 0.0;

    }  // if
  }

  // write new format
  EcalCondHeader h;
  EcalIntercalibConstantsXMLTranslator::writeXML(xmlfilename, h, rcd);
}