InvMatrixUtils.cc File Reference

#include "Calibration/Tools/interface/InvMatrixUtils.h"
#include "Calibration/Tools/interface/InvMatrixCommonDefs.h"
#include "TStyle.h"
#include "TROOT.h"
#include "CLHEP/Geometry/Point3D.h"
#include "Calibration/Tools/interface/matrixSaver.h"
#include <iostream>
#include <string>
#include <fstream>
#include <sstream>
#include <vector>
#include <cassert>

Go to the source code of this file.

Functions
double	effectiveSigma (TH1F &histogram, int vSteps)
double	etaCorrE1E25 (int eta)
double	etaCorrE1E49 (int eta)
double	etaCorrE1E9 (int eta)
int	etaFromXtal (const int &xtal)
int	extract (std::vector< int > *output, const std::string &dati)
std::pair< int, int >	findSupport (TH1F &histogram, double thres)
double	get3x3 (const Float_t energy[7][7])
double	get5x5 (const Float_t energy[7][7])
TCanvas *	getGlobalCanvas (std::string name)
TFile *	getGlobalTFile (std::string name)
CLHEP::HepMatrix *	getSavedMatrix (const std::string &name)
int	ietaFromXtal (const int &xtal)
int	iphiFromXtal (const int &xtal)
void	mtrTransfer (double output[85][20], CLHEP::HepMatrix *input, double Default)
int	phiFromXtal (const int &xtal)
int	readCMSSWcoeff (CLHEP::HepMatrix &calibcoeff, const std::string &inputFileName, double defaultVal)
int	readCMSSWcoeffForComparison (CLHEP::HepMatrix &calibcoeff, const std::string &inputFileName)
int	saveGlobalTFile (std::string name)
void	setStyle ()
TH1D *	smartError (TH1D *strip)
TH1D *	smartGausProfile (TH2F *strip, double width)
TH1D *	smartProfile (TH2F *strip, double width)
HepGeom::Point3D< Float_t >	TBposition (const Float_t amplit[7][7], const Float_t beamEne, const Float_t w0, const Float_t x0, const Float_t a0, const Float_t sideX, const Float_t sideY)
int	translateCoeff (const CLHEP::HepMatrix &calibcoeff, const CLHEP::HepMatrix &sigmaMatrix, const CLHEP::HepMatrix &statisticMatrix, std::string SMnumber, double calibThres, std::string fileName, std::string genTag, std::string method, std::string version, std::string type)
int	writeCalibTxt (const CLHEP::HepMatrix &AmplitudeMatrix, const CLHEP::HepMatrix &SigmaMatrix, const CLHEP::HepMatrix &StatisticMatrix, std::string fileName)
int	writeCMSSWCoeff (const CLHEP::HepMatrix &amplMatrix, double calibThres, float ERef, const CLHEP::HepMatrix &sigmaMatrix, const CLHEP::HepMatrix &statisticMatrix, std::string fileName, std::string genTag, std::string method, std::string version, std::string type)
int	writeCMSSWCoeff (const CLHEP::HepMatrix &amplMatrix, double calibThres, int etaRef, int phiRef, const CLHEP::HepMatrix &sigmaMatrix, const CLHEP::HepMatrix &statisticMatrix, std::string fileName, std::string genTag, std::string method, std::string version, std::string type)
int	xtalFromEtaPhi (const int &myEta, const int &myPhi)
int	xtalFromiEtaiPhi (const int &iEta, const int &iPhi)

Function Documentation

effectiveSigma()

double effectiveSigma	(	TH1F &	histogram,
		int	vSteps = 100
	)

find the effective sigma as the half width of the sub-distribution containing 68.3% of the total distribution

Definition at line 619 of file InvMatrixUtils.cc.

                                                   {
  double totInt = histogram.Integral();
  int maxBin = histogram.GetMaximumBin();
  int maxBinVal = int(histogram.GetBinContent(maxBin));
  int totBins = histogram.GetNbinsX();
  double area = totInt;
  double threshold = 0;
  double vStep = maxBinVal / vSteps;
  int leftBin = 1;
  int rightBin = totBins - 1;
  //loop over the vertical range
  while (area / totInt > 0.683) {
    threshold += vStep;
    // loop toward the left
    for (int back = maxBin; back > 0; --back) {
      if (histogram.GetBinContent(back) < threshold) {
        leftBin = back;
        break;
      }
    }  // loop toward the left
 
    // loop toward the right
    for (int fwd = maxBin; fwd < totBins; ++fwd) {
      if (histogram.GetBinContent(fwd) < threshold) {
        rightBin = fwd;
        break;
      }
    }  // loop toward the right
    area = histogram.Integral(leftBin, rightBin);
  }  //loop over the vertical range
 
  histogram.GetXaxis()->SetRange(leftBin, rightBin);
  // double sigmaEff = histogram.GetRMS () ;
  double halfWidthRange = 0.5 * (histogram.GetBinCenter(rightBin) - histogram.GetBinCenter(leftBin));
  return halfWidthRange;
}

References custom_jme_cff::area, createfilelist::int, cms::cuda::allocator::maxBin, and remoteMonitoring_LED_IterMethod_cfg::threshold.

etaCorrE1E25()

double etaCorrE1E25

(

int

eta

)

correction for eta containment for 5*5 cluster

Definition at line 697 of file InvMatrixUtils.cc.

                             {
  double p0 = 0.807883;
  double p1 = 0.000182551;
  double p2 = -5.76961e-06;
  double p3 = 7.41903e-08;
  double p4 = -2.25384e-10;
 
  double corr;
  if (eta < 6)
    corr = p0;
  else
    corr = p0 + p1 * eta + p2 * eta * eta + p3 * eta * eta * eta + p4 * eta * eta * eta * eta;
  return corr / p0;
}

References alignCSCRings::corr, PVValHelper::eta, p1, p2, p3, and p4.

etaCorrE1E49()

double etaCorrE1E49

(

int

eta

)

correction for eta containment for 7*7 cluster

Definition at line 713 of file InvMatrixUtils.cc.

                             {
  double p0 = 0.799895;
  double p1 = 0.000235487;
  double p2 = -8.26496e-06;
  double p3 = 1.21564e-07;
  double p4 = -4.83286e-10;
 
  double corr;
  if (eta < 8)
    corr = p0;
  else
    corr = p0 + p1 * eta + p2 * eta * eta + p3 * eta * eta * eta + p4 * eta * eta * eta * eta;
  return corr / p0;
}

References alignCSCRings::corr, PVValHelper::eta, p1, p2, p3, and p4.

etaCorrE1E9()

double etaCorrE1E9

(

int

eta

)

correction for eta containment for 3*3 cluster

Definition at line 729 of file InvMatrixUtils.cc.

                            {
  if (eta < 4)
    return 1.0;
  // grazie Paolo
  double p0 = 0.834629;
  double p1 = 0.00015254;
  double p2 = -4.91784e-06;
  double p3 = 6.54652e-08;
  double p4 = -2.4894e-10;
 
  double corr;
  if (eta < 6)
    corr = p0;
  else
    corr = p0 + p1 * eta + p2 * eta * eta + p3 * eta * eta * eta + p4 * eta * eta * eta * eta;
  return corr / p0;
}

References alignCSCRings::corr, PVValHelper::eta, p1, p2, p3, and p4.

etaFromXtal()

int etaFromXtal

(

const int &

xtal

)

get the eta coord [0,84]

Definition at line 229 of file InvMatrixUtils.cc.

                                 {
  //  return floor (static_cast<double> ((xtal-1) / 20)) ;
  return int(floor((xtal - 1) / 20));
}

References createfilelist::int.

Referenced by ietaFromXtal(), phiFromXtal(), readCMSSWcoeff(), and readCMSSWcoeffForComparison().

extract()

int extract	(	std::vector< int > *	output,
		const std::string &	dati
	)

to read a file containing unserted integers while avoiding comment lines

Definition at line 251 of file InvMatrixUtils.cc.

                                                           {
  std::ifstream _dati(dati.c_str());
  // loop over the file
  while (!_dati.eof()) {
    // get the line
    std::string dataline;
    do {
      getline(_dati, dataline, '\n');
    } while (*dataline.begin() == '#');
    std::stringstream linea(dataline);
    // loop over the line
    while (!linea.eof()) {
      int buffer = -1;
      linea >> buffer;
      if (buffer != -1)
        output->push_back(buffer);
    }  // loop over the line
  }    // loop over the file
  return output->size();
}

References edmScanValgrind::buffer, convertSQLitetoXML_cfg::output, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by clone(), and SummaryPlotFactory< CommissioningAnalysis * >::init().

findSupport()

std::pair<int, int> findSupport	(	TH1F &	histogram,
		double	thres = 0.
	)

find the support of the histogram above a threshold return the min and max bins

Definition at line 658 of file InvMatrixUtils.cc.

                                                             {
  int totBins = histogram.GetNbinsX();
  if (thres >= histogram.GetMaximum())
    return std::pair<int, int>(0, totBins);
 
  int leftBin = totBins - 1;
  // search from left for the minimum
  for (int bin = 1; bin < totBins; ++bin) {
    if (histogram.GetBinContent(bin) > thres) {
      leftBin = bin;
      break;
    }
  }  // search from left for the minimum
  int rightBin = 1;
  // search from right for the maximum
  for (int bin = totBins - 1; bin > 0; --bin) {
    if (histogram.GetBinContent(bin) > thres) {
      rightBin = bin;
      break;
    }
  }  // search from right for the maximum
  return std::pair<int, int>(leftBin, rightBin);
}

References newFWLiteAna::bin.

get3x3()

double get3x3

(

const Float_t

energy[7][7]

)

get the energy in the 3x3 from the 7x7 array around the most energetic crystal

Definition at line 165 of file InvMatrixUtils.cc.

                                          {
  double total = 0.;
 
  for (int eta = 2; eta < 5; ++eta)
    for (int phi = 2; phi < 5; ++phi)
      total += energy[eta][phi];
 
  return total;
}

References HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, PVValHelper::phi, and dqmMemoryStats::total.

Referenced by TBposition().

get5x5()

double get5x5

(

const Float_t

energy[7][7]

)

get the energy in the 5x5 from the 7x7 array around the most energetic crystal

Definition at line 151 of file InvMatrixUtils.cc.

                                          {
  double total = 0.;
 
  for (int eta = 1; eta < 6; ++eta)
    for (int phi = 1; phi < 6; ++phi)
      total += energy[eta][phi];
 
  return total;
}

References HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, PVValHelper::phi, and dqmMemoryStats::total.

getGlobalCanvas()

TCanvas* getGlobalCanvas

(

std::string

name = "Inv MatrixCanvas"

)

search for an existing canvas with the name and returns the poiter to it

Definition at line 39 of file InvMatrixUtils.cc.

                                         {
  setStyle();
  TCanvas *globalCanvas = static_cast<TCanvas *>(gROOT->FindObject(name.c_str()));
  if (globalCanvas) {
    globalCanvas->Clear();
    globalCanvas->UseCurrentStyle();
    globalCanvas->SetWindowSize(700, 600);
 
  } else {
    globalCanvas = new TCanvas(name.c_str(), name.c_str(), 700, 600);
  }
  return globalCanvas;
}

References Skims_PA_cff::name, and setStyle().

getGlobalTFile()

TFile* getGlobalTFile

(

std::string

name = "Inv MatrixTFile.root"

)

search for an existing TFile with the name and returns the poiter to it

Definition at line 55 of file InvMatrixUtils.cc.

                                      {
  //    std::cout << "writing " << name << std::endl ;
  //    setStyle () ;
  TFile *globalTFile = (TFile *)gROOT->FindObject(name.c_str());
  if (!globalTFile) {
    //        std::cout << "does not exist. creating it " << std::endl;
    globalTFile = new TFile(name.c_str(), "RECREATE");
  }
 
  return globalTFile;
}

References Skims_PA_cff::name.

getSavedMatrix()

CLHEP::HepMatrix* getSavedMatrix

(

const std::string &

name

)

search for an existing calib matrix saved with the name and returns the poiter to it, the deletion is responsiblity of the user

Definition at line 81 of file InvMatrixUtils.cc.

                                                      {
  matrixSaver reader;
  CLHEP::HepMatrix *savedMatrix;
  if (reader.touch(name)) {
    savedMatrix = static_cast<CLHEP::HepMatrix *>(reader.getMatrix(name));
  } else {
    savedMatrix = new CLHEP::HepMatrix(SCMaxEta, SCMaxPhi, 0);
  }
 
  return savedMatrix;
}

References Skims_PA_cff::name, DQM::reader, SCMaxEta, and SCMaxPhi.

ietaFromXtal()

int ietaFromXtal

(

const int &

xtal

)

get the eta coord [1,85]

Definition at line 243 of file InvMatrixUtils.cc.

{ return etaFromXtal(xtal) + 1; }

References etaFromXtal().

iphiFromXtal()

int iphiFromXtal

(

const int &

xtal

)

get the phi coord [1,20]

Definition at line 247 of file InvMatrixUtils.cc.

{ return phiFromXtal(xtal) + 1; }

References phiFromXtal().

mtrTransfer()

void mtrTransfer	(	double	output[85][20],
		CLHEP::HepMatrix *	input,
		double	Default
	)

transfers a CLHEP matrix into a double array with the size of a supermodule

Definition at line 684 of file InvMatrixUtils.cc.

                                                                                           {
  for (int eta = 0; eta < SCMaxEta; ++eta)
    for (int phi = 0; phi < SCMaxPhi; ++phi) {
      if ((*input)[eta][phi])
        output[eta][phi] = (*input)[eta][phi];
      else
        output[eta][phi] = Default;
    }
  return;
}

References Default, PVValHelper::eta, input, convertSQLitetoXML_cfg::output, PVValHelper::phi, SCMaxEta, and SCMaxPhi.

phiFromXtal()

int phiFromXtal

(

const int &

xtal

)

get the phi coord [0,19]

Definition at line 236 of file InvMatrixUtils.cc.

                                 {
  int phi = (xtal - 1) - 20 * etaFromXtal(xtal);
  return (20 - phi - 1);
}

References etaFromXtal(), and PVValHelper::phi.

Referenced by iphiFromXtal(), readCMSSWcoeff(), readCMSSWcoeffForComparison(), and xtalFromEtaPhi().

readCMSSWcoeff()

int readCMSSWcoeff	(	CLHEP::HepMatrix &	calibcoeff,
		const std::string &	inputFileName,
		double	defaultVal = 1.
	)

translates the calib coefficients format, from the CMSSW one to the TB06Studies one

Definition at line 443 of file InvMatrixUtils.cc.

                                                                                                  {
  std::ifstream CMSSWfile;
  CMSSWfile.open(inputFileName.c_str());
  std::string buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  while (!CMSSWfile.eof()) {
    int xtalnum;
    CMSSWfile >> xtalnum;
    double coeff;
    CMSSWfile >> coeff;
    double buffer;
    CMSSWfile >> buffer;
    int good;
    CMSSWfile >> good;
    CMSSWfile >> good;
    if (!good)
      coeff = defaultVal;  //FIXME 0 o 1?
    calibcoeff[etaFromXtal(xtalnum)][phiFromXtal(xtalnum)] = coeff;
  }
  return 0;
}

References edmScanValgrind::buffer, etaFromXtal(), InefficientDoubleROC::inputFileName, phiFromXtal(), and AlCaHLTBitMon_QueryRunRegistry::string.

readCMSSWcoeffForComparison()

int readCMSSWcoeffForComparison	(	CLHEP::HepMatrix &	calibcoeff,
		const std::string &	inputFileName
	)

translates the calib coefficients format, from the CMSSW one to the TB06Studies one

Definition at line 472 of file InvMatrixUtils.cc.

                                                                                            {
  std::ifstream CMSSWfile;
  CMSSWfile.open(inputFileName.c_str());
  std::string buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  CMSSWfile >> buffer;
  while (!CMSSWfile.eof()) {
    int xtalnum;
    CMSSWfile >> xtalnum;
    double coeff;
    CMSSWfile >> coeff;
    double buffer;
    CMSSWfile >> buffer;
    int good;
    CMSSWfile >> good;
    CMSSWfile >> good;
    if (!good)
      coeff = 0.;  //FIXME 0 o 1?
    calibcoeff[etaFromXtal(xtalnum)][phiFromXtal(xtalnum)] = coeff;
  }
  return 0;
}

References edmScanValgrind::buffer, etaFromXtal(), InefficientDoubleROC::inputFileName, phiFromXtal(), and AlCaHLTBitMon_QueryRunRegistry::string.

saveGlobalTFile()

int saveGlobalTFile

(

std::string

name = "Inv MatrixFile.root"

)

search for an existing TFile with the name and saves it to disk with his name

Definition at line 69 of file InvMatrixUtils.cc.

                                    {
  TFile *globalTFile = static_cast<TFile *>(gROOT->FindObject(name.c_str()));
  if (!globalTFile)
    return 1;
  globalTFile->Write();
  globalTFile->Close();
  delete globalTFile;
  return 0;
}

References Skims_PA_cff::name.

setStyle()

void setStyle

(

)

set the style for the printout

Definition at line 20 of file InvMatrixUtils.cc.

                {
  gROOT->SetStyle("Plain");
  gStyle->SetTextSize(0.5);
  //gStyle->SetOptStat (1111111) ;
  gStyle->SetOptStat(0);
  //gStyle->SetOptFit (1111) ;
  gStyle->SetOptFit(0);
  gStyle->SetTitleBorderSize(0);
  gStyle->SetTitleX(0.08);
  gStyle->SetTitleY(0.97);
  gStyle->SetPalette(1, nullptr);
  gStyle->SetStatColor(0);
  gStyle->SetFrameFillStyle(0);
  gStyle->SetFrameFillColor(0);
  return;
}

Referenced by getGlobalCanvas().

smartError()

TH1D* smartError

(

TH1D *

strip

)

Definition at line 598 of file InvMatrixUtils.cc.

                              {
  double xmin = strip->GetXaxis()->GetXmin();
  double xmax = strip->GetXaxis()->GetXmax();
  int stripsBins = strip->GetNbinsX();
 
  std::string name = strip->GetName();
  name += "_error";
  TH1D *error = new TH1D(name.c_str(), strip->GetTitle(), stripsBins, xmin, xmax);
 
  int binmin = 1;
  int ngroup = 1;  // bins per step
  int binmax = stripsBins;
  for (int bin = binmin; bin <= binmax; bin += ngroup) {
    double dummyError = strip->GetBinError(bin);
    error->SetBinContent(bin, dummyError);
  }
  return error;
}

References newFWLiteAna::bin, relativeConstraints::error, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, digitizers_cfi::strip, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

smartGausProfile()

TH1D* smartGausProfile	(	TH2F *	strip,
		double	width
	)

smart profiling by fixing gaussian parameters and range from a first averaging

Definition at line 546 of file InvMatrixUtils.cc.

                                                  {
  TProfile *stripProfile = strip->ProfileX();
 
  // (from FitSlices of TH2.h)
 
  double xmin = stripProfile->GetXaxis()->GetXmin();
  double xmax = stripProfile->GetXaxis()->GetXmax();
  int profileBins = stripProfile->GetNbinsX();
 
  std::string name = strip->GetName();
  name += "_smartGaus";
  TH1D *prof = new TH1D(name.c_str(), strip->GetTitle(), profileBins, xmin, xmax);
 
  int cut = 0;  // minimum number of entries per fitted bin
  int nbins = strip->GetXaxis()->GetNbins();
  int binmin = 1;
  int ngroup = 1;  // bins per step
  int binmax = nbins;
 
  // loop over the strip bins
  for (int bin = binmin; bin <= binmax; bin += ngroup) {
    TH1D *hpy = strip->ProjectionY("_temp", bin, bin + ngroup - 1, "e");
    if (hpy == nullptr)
      continue;
    int nentries = Int_t(hpy->GetEntries());
    if (nentries == 0 || nentries < cut) {
      delete hpy;
      continue;
    }
 
    Int_t biny = bin + ngroup / 2;
 
    TF1 *gaussian =
        new TF1("gaussian", "gaus", hpy->GetMean() - width * hpy->GetRMS(), hpy->GetMean() + width * hpy->GetRMS());
    gaussian->SetParameter(1, hpy->GetMean());
    gaussian->SetParameter(2, hpy->GetRMS());
    hpy->Fit("gaussian", "RQL");
 
    hpy->GetXaxis()->SetRangeUser(hpy->GetMean() - width * hpy->GetRMS(), hpy->GetMean() + width * hpy->GetRMS());
    prof->Fill(strip->GetXaxis()->GetBinCenter(biny), gaussian->GetParameter(1));
    prof->SetBinError(biny, gaussian->GetParameter(2));
 
    delete gaussian;
    delete hpy;
  }  // loop over the bins
 
  delete stripProfile;
  return prof;
}

References newFWLiteAna::bin, TkAlMuonSelectors_cfi::cut, Skims_PA_cff::name, LaserClient_cfi::nbins, AlCaHLTBitMon_QueryRunRegistry::string, digitizers_cfi::strip, tools::TF1, ApeEstimator_cff::width, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

smartProfile()

TH1D* smartProfile	(	TH2F *	strip,
		double	width
	)

smart profiling by double averaging

Definition at line 501 of file InvMatrixUtils.cc.

                                              {
  TProfile *stripProfile = strip->ProfileX();
 
  // (from FitSlices of TH2.h)
 
  double xmin = stripProfile->GetXaxis()->GetXmin();
  double xmax = stripProfile->GetXaxis()->GetXmax();
  int profileBins = stripProfile->GetNbinsX();
 
  std::string name = strip->GetName();
  name += "_smart";
  TH1D *prof = new TH1D(name.c_str(), strip->GetTitle(), profileBins, xmin, xmax);
 
  int cut = 0;  // minimum number of entries per fitted bin
  int nbins = strip->GetXaxis()->GetNbins();
  int binmin = 1;
  int ngroup = 1;  // bins per step
  int binmax = nbins;
 
  // loop over the strip bins
  for (int bin = binmin; bin <= binmax; bin += ngroup) {
    TH1D *hpy = strip->ProjectionY("_temp", bin, bin + ngroup - 1, "e");
    if (hpy == nullptr)
      continue;
    int nentries = Int_t(hpy->GetEntries());
    if (nentries == 0 || nentries < cut) {
      delete hpy;
      continue;
    }
 
    Int_t biny = bin + ngroup / 2;
 
    hpy->GetXaxis()->SetRangeUser(hpy->GetMean() - width * hpy->GetRMS(), hpy->GetMean() + width * hpy->GetRMS());
    prof->Fill(strip->GetXaxis()->GetBinCenter(biny), hpy->GetMean());
    prof->SetBinError(biny, hpy->GetRMS());
 
    delete hpy;
  }  // loop over the bins
 
  delete stripProfile;
  return prof;
}

References newFWLiteAna::bin, TkAlMuonSelectors_cfi::cut, Skims_PA_cff::name, LaserClient_cfi::nbins, AlCaHLTBitMon_QueryRunRegistry::string, digitizers_cfi::strip, ApeEstimator_cff::width, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

TBposition()

HepGeom::Point3D<Float_t> TBposition	(	const Float_t	amplit[7][7],
		const Float_t	beamEne,
		const Float_t	w0,
		const Float_t	x0,
		const Float_t	a0,
		const Float_t	sideX,
		const Float_t	sideY
	)

return the impact position of the electron over ECAL

Definition at line 106 of file InvMatrixUtils.cc.

                                                          {
  // variables
  Float_t caloX = 0.;
  Float_t caloY = 0.;
  Float_t sumWeight = 0.;
  Float_t depth = x0 * (log(beamEne) + a0);  // shower depthh, in mm
  Float_t sin3 = 0.052335956;                // sin (3 degrees) , sin3 = sin(3.*3.141592654/180.)
 
  Float_t invE3x3 = 1. / get3x3(amplit);
 
  // loop over 3x3 crystals
  for (int eta = 2; eta <= 4; eta++) {
    for (int phi = 2; phi <= 4; phi++) {
      Float_t weight = log(amplit[eta][phi] * invE3x3) + w0;
      if (weight > 0) {
        caloX += (eta - 3) * sideX * weight;
        caloY -= (phi - 3) * sideY * weight;
        sumWeight += weight;
      }
    }
  }
 
  caloX /= sumWeight;
  caloY /= sumWeight;
 
  // correction for shower depthh
  caloX -= depth * sin3;
  caloY -= depth * sin3;
 
  // FIXME the z set to zero
  HepGeom::Point3D<Float_t> TBposition(caloX, caloY, 0);
 
  return TBposition;
}

References a0, LEDCalibrationChannels::depth, PVValHelper::eta, get3x3(), dqm-mbProfile::log, PVValHelper::phi, and mps_merge::weight.

translateCoeff()

int translateCoeff	(	const CLHEP::HepMatrix &	calibcoeff,
		const CLHEP::HepMatrix &	sigmaMatrix,
		const CLHEP::HepMatrix &	statisticMatrix,
		std::string	SMnumber = "1",
		double	calibThres = 0.01,
		std::string	fileName = "calibOutput.txt",
		std::string	genTag = "CAL_GENTAG",
		std::string	method = "CAL_METHOD",
		std::string	version = "CAL_VERSION",
		std::string	type = "CAL_TYPE"
	)

translates the calib coefficients format, from the TB06Studies one to the CMSSSW one

Definition at line 403 of file InvMatrixUtils.cc.

                                   {
  // open the file for output
  std::ofstream txt_outfile;
  txt_outfile.open(fileName.c_str());
  txt_outfile << SMnumber << "\n";  // super-module number
  txt_outfile << "-1\n";            // number of events
  txt_outfile << genTag << "\n";
  txt_outfile << method << "\n";
  txt_outfile << version << "\n";
  txt_outfile << type << "\n";
 
  // loop over crystals
  for (int eta = 0; eta < SCMaxEta; ++eta)
    for (int phi = 0; phi < SCMaxPhi; ++phi) {
      if (calibcoeff[eta][phi] < calibThres) {
        txt_outfile << xtalFromiEtaiPhi(eta + 1, phi + 1) << "\t" << 1 << "\t" << -1 << "\t" << -1 << "\t" << 0 << "\n";
        std::cout << "[translateCoefff][" << SMnumber << "]\t WARNING crystal " << xtalFromiEtaiPhi(eta + 1, phi + 1)
                  << " calib coeff below threshold: "
                  << "\t" << 1 << "\t" << -1 << "\t" << -1 << "\t" << 0 << "\n";
      } else
        txt_outfile << xtalFromiEtaiPhi(eta + 1, phi + 1) << "\t" << calibcoeff[eta][phi] << "\t"
                    << sigmaMatrix[eta][phi] << "\t" << statisticMatrix[eta][phi] << "\t" << 1 << "\n";
    }  // loop over crystals
 
  // save and close the file
  txt_outfile.close();
  return 0;
}

References gather_cfg::cout, PVValHelper::eta, MillePedeFileConverter_cfg::fileName, AlcaSiPixelAliHarvester0T_cff::method, PVValHelper::phi, SCMaxEta, SCMaxPhi, BeamSplash_cfg::version, and xtalFromiEtaiPhi().

writeCalibTxt()

int writeCalibTxt	(	const CLHEP::HepMatrix &	AmplitudeMatrix,
		const CLHEP::HepMatrix &	SigmaMatrix,
		const CLHEP::HepMatrix &	StatisticMatrix,
		std::string	fileName = "calibOutput.txt"
	)

to write the calibration constants file

Definition at line 276 of file InvMatrixUtils.cc.

                                      {
  // look for the reference crystal
  double reference = 0.;
  for (int eta = 0; eta < SCMaxEta; ++eta)
    for (int phi = 0; phi < SCMaxPhi; ++phi) {
      if (AmplitudeMatrix[eta][phi] && SigmaMatrix[eta][phi] < 100 /*FIXME sigmaCut*/) {
        reference = AmplitudeMatrix[eta][phi];
        std::cout << "[InvMatrixUtils][writeCalibTxt] reference crystal: "
                  << "(" << eta << "," << phi << ") -> " << reference << "\n";
        break;
      }
    }
  if (!reference) {
    std::cerr << "ERROR: no calibration coefficients found" << std::endl;
    return 1;
  }
 
  // open the file for output
  std::ofstream txt_outfile;
  txt_outfile.open(fileName.c_str());
  txt_outfile << "# xtal\tcoeff\tsigma\tevt\tisGood\n";
 
  // loop over the crystals
  for (int eta = 0; eta < SCMaxEta; ++eta)
    for (int phi = 0; phi < SCMaxPhi; ++phi) {
      int isGood = 1;
      if (AmplitudeMatrix[eta][phi] == 0)
        isGood = 0;
      if (SigmaMatrix[eta][phi] > 100 /*FIXME sigmaCut*/)
        isGood = 0;
      txt_outfile << xtalFromEtaPhi(eta, phi) << "\t" << AmplitudeMatrix[eta][phi] / reference << "\t"
                  << SigmaMatrix[eta][phi] << "\t" << StatisticMatrix[eta][phi] << "\t" << isGood << "\n";
    }
 
  // save and close the file
  txt_outfile.close();
  return 0;
}

References beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, PVValHelper::eta, MillePedeFileConverter_cfg::fileName, PVValHelper::phi, RecoTauValidation_cfi::reference, SCMaxEta, SCMaxPhi, and xtalFromEtaPhi().

writeCMSSWCoeff() [1/2]

int writeCMSSWCoeff	(	const CLHEP::HepMatrix &	amplMatrix,
		double	calibThres,
		float	ERef,
		const CLHEP::HepMatrix &	sigmaMatrix,
		const CLHEP::HepMatrix &	statisticMatrix,
		std::string	fileName = "calibOutput.txt",
		std::string	genTag = "CAL_GENTAG",
		std::string	method = "CAL_METHOD",
		std::string	version = "CAL_VERSION",
		std::string	type = "CAL_TYPE"
	)

to write the file fpr the CMSSW in the DB compliant format (using Energy as reference)

Definition at line 320 of file InvMatrixUtils.cc.

                                    {
  // open the file for output
  std::ofstream txt_outfile;
  txt_outfile.open(fileName.c_str());
  txt_outfile << "1\n";   // super-module number
  txt_outfile << "-1\n";  // number of events
  txt_outfile << genTag << "\n";
  txt_outfile << method << "\n";
  txt_outfile << version << "\n";
  txt_outfile << type << "\n";
 
  double reference = ERef;
 
  // loop over crystals
  for (int eta = 0; eta < SCMaxEta; ++eta)
    for (int phi = 0; phi < SCMaxPhi; ++phi) {
      if (amplMatrix[eta][phi] <= calibThres)
        txt_outfile << xtalFromiEtaiPhi(eta + 1, phi + 1) << "\t" << 1 << "\t" << -1 << "\t" << -1 << "\t" << 0 << "\n";
      else
        txt_outfile << xtalFromiEtaiPhi(eta + 1, phi + 1) << "\t" << reference / amplMatrix[eta][phi] << "\t"
                    << sigmaMatrix[eta][phi] << "\t" << statisticMatrix[eta][phi] << "\t" << 1 << "\n";
    }  // loop over crystals
 
  // save and close the file
  txt_outfile.close();
  return 0;
}

References PVValHelper::eta, MillePedeFileConverter_cfg::fileName, AlcaSiPixelAliHarvester0T_cff::method, PVValHelper::phi, RecoTauValidation_cfi::reference, SCMaxEta, SCMaxPhi, BeamSplash_cfg::version, and xtalFromiEtaiPhi().

writeCMSSWCoeff() [2/2]

int writeCMSSWCoeff	(	const CLHEP::HepMatrix &	amplMatrix,
		double	calibThres,
		int	etaRef,
		int	phiRef,
		const CLHEP::HepMatrix &	sigmaMatrix,
		const CLHEP::HepMatrix &	statisticMatrix,
		std::string	fileName = "calibOutput.txt",
		std::string	genTag = "CAL_GENTAG",
		std::string	method = "CAL_METHOD",
		std::string	version = "CAL_VERSION",
		std::string	type = "CAL_TYPE"
	)

to write the file fpr the CMSSW in the DB compliant format (using Crystal as reference)

Definition at line 359 of file InvMatrixUtils.cc.

                                    {
  // open the file for output
  std::ofstream txt_outfile;
  txt_outfile.open(fileName.c_str());
  txt_outfile << "1\n";   // super-module number
  txt_outfile << "-1\n";  // number of events
  txt_outfile << genTag << "\n";
  txt_outfile << method << "\n";
  txt_outfile << version << "\n";
  txt_outfile << type << "\n";
 
  if (amplMatrix[etaRef - 1][phiRef - 1] == 0) {
    std::cerr << "The reference crystal: (" << etaRef << "," << phiRef << ") is out of range\n";
    return 1;
  }
  double reference = amplMatrix[etaRef - 1][phiRef - 1];
 
  // loop over crystals
  for (int eta = 0; eta < SCMaxEta; ++eta)
    for (int phi = 0; phi < SCMaxPhi; ++phi) {
      if (amplMatrix[eta][phi] <= calibThres)
        txt_outfile << xtalFromiEtaiPhi(eta + 1, phi + 1) << "\t" << 1 << "\t" << -1 << "\t" << -1 << "\t" << 0 << "\n";
      else
        txt_outfile << xtalFromiEtaiPhi(eta + 1, phi + 1) << "\t" << reference / amplMatrix[eta][phi] << "\t"
                    << sigmaMatrix[eta][phi] << "\t" << statisticMatrix[eta][phi] << "\t" << 1 << "\n";
    }  // loop over crystals
 
  // save and close the file
  txt_outfile.close();
  return 0;
}

References beam_dqm_sourceclient-live_cfg::cerr, PVValHelper::eta, MillePedeFileConverter_cfg::fileName, AlcaSiPixelAliHarvester0T_cff::method, PVValHelper::phi, RecoTauValidation_cfi::reference, SCMaxEta, SCMaxPhi, BeamSplash_cfg::version, and xtalFromiEtaiPhi().

xtalFromEtaPhi()

int xtalFromEtaPhi	(	const int &	myEta,
		const int &	myPhi
	)

to get the parameters from a congiguration file

Definition at line 204 of file InvMatrixUtils.cc.

                                                       {
  int xMin = 20 * myEta + 1;
  int xMax = 20 * (myEta + 1) + 1;
 
  int myCryst = 999999;
 
  for (int x = xMin; x < xMax; x++) {
    if (phiFromXtal(x) == myPhi)
      myCryst = x;
  }
  return myCryst;
}

References phiFromXtal(), vertices_cff::x, multiplicitycorr_cfi::xMax, and photonAnalyzer_cfi::xMin.

Referenced by writeCalibTxt().

xtalFromiEtaiPhi()

int xtalFromiEtaiPhi	(	const int &	iEta,
		const int &	iPhi
	)

to get the crystal number from iEta and iPhi iEta runs from 1 to 85 iPhi runs from 1 to 20

Definition at line 219 of file InvMatrixUtils.cc.

                                                       {
  assert(iEta >= 1);
  assert(iEta <= 85);
  assert(iPhi >= 1);
  assert(iPhi <= 20);
  return 20 * (iEta - 1) + 21 - iPhi;
}

References cms::cuda::assert(), and L1TowerCalibrationProducer_cfi::iEta.

Referenced by translateCoeff(), and writeCMSSWCoeff().