EcalTPGParamBuilder.cc File Reference

#include "EcalTPGParamBuilder.h"
#include "EcalTPGDBApp.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "CondFormats/EcalObjects/interface/EcalMGPAGainRatio.h"
#include "CalibCalorimetry/EcalLaserCorrection/interface/EcalLaserDbService.h"
#include "CalibCalorimetry/EcalLaserCorrection/interface/EcalLaserDbRecord.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalSimParameterMap.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <TF1.h>
#include <TH2F.h>
#include <TFile.h>
#include <TNtuple.h>
#include <ctime>
#include <iomanip>
#include <fstream>
#include <iostream>
#include <string>
#include <sstream>

Go to the source code of this file.

Functions
double	oneOverEtResolEt (double *x, double *par)

Function Documentation

oneOverEtResolEt()

double oneOverEtResolEt	(	double *	x,
		double *	par
	)

Definition at line 32 of file EcalTPGParamBuilder.cc.

References MillePedeFileConverter_cfg::e, mathSSE::sqrt(), and x.

Referenced by EcalTPGParamBuilder::computeLUT().

                                                {
  double Et = x[0];
  if (Et < 1e-6)
    return 1. / par[1];  // to avoid division by 0.
  double resolEt_overEt =
      sqrt((par[0] / sqrt(Et)) * (par[0] / sqrt(Et)) + (par[1] / Et) * (par[1] / Et) + par[2] * par[2]);
  return 1. / (Et * resolEt_overEt);
}