Pi0CalibXMLwriter.cc

Go to the documentation of this file.

#include "Calibration/Tools/interface/Pi0CalibXMLwriter.h"
#include <string>
#include <cstdio>

Pi0CalibXMLwriter::Pi0CalibXMLwriter(EcalSubdetector subdet) : subdet_(subdet)
{
    
    char filename[128];
    if (subdet_==EcalEndcap) {
      sprintf(filename,"endcap_iniCalib.xml");
    } else {
      sprintf(filename,"barrel_iniCalib.xml");
    }
    FILENAME = fopen(filename,"w");
    fprintf(FILENAME,"<?xml version=\"1.0\" ?>\n");
    fprintf(FILENAME,"<CalibrationConstants>\n");
    if (subdet==EcalEndcap) {
      fprintf(FILENAME,"<EcalEndcap>\n");
    } else {
      fprintf(FILENAME,"<EcalBarrel>\n");
    }

}

Pi0CalibXMLwriter::Pi0CalibXMLwriter(EcalSubdetector subdet, int loop) : subdet_(subdet), loop_(loop)
{
    
    char filename[128];
    if (subdet_==EcalEndcap) {
      sprintf(filename,"endcap_calib_loop_%d.xml",loop);
    } else {
      sprintf(filename,"barrel_calib_loop_%d.xml",loop);
    }
    FILENAME = fopen(filename,"w");
    fprintf(FILENAME,"<?xml version=\"1.0\" ?>\n");
    fprintf(FILENAME,"<CalibrationConstants>\n");
    if (subdet==EcalEndcap) {
      fprintf(FILENAME,"<EcalEndcap>\n");
    } else {
      fprintf(FILENAME,"<EcalBarrel>\n");
    }

}

Pi0CalibXMLwriter::~Pi0CalibXMLwriter()
{
    if (subdet_==EcalEndcap) {
      fprintf(FILENAME,"<EcalEndcap>\n");
    } else {
      fprintf(FILENAME,"<EcalBarrel>\n");
    }
    fprintf(FILENAME,"</CalibrationConstants>\n");
    fclose(FILENAME);
}

void Pi0CalibXMLwriter::writeLine(EBDetId const & det, float calib)
{
int eta=det.ieta();
int phi=det.iphi();
fprintf(FILENAME,"<Cell eta_index=\"%d\" phi_index=\"%d\" scale_factor=\"%f\"/>\n",eta,phi,calib);
}


void Pi0CalibXMLwriter::writeLine(EEDetId const & det, float calib)
{
int x=det.ix();
int y=det.iy();
int z=det.zside()>0 ? 1 : 0;
fprintf(FILENAME,"<Cell x_index=\"%d\" y_index=\"%d\" z_index=\"%d\" scale_factor=\"%f\"/>\n",x,y,z,calib);
}