HFDarkening Class Reference

#include <HFDarkening.h>

Public Member Functions
double	degradation (double mrad)
double	dose (unsigned int layer, double radius)
	HFDarkening (const edm::ParameterSet &pset)
double	int_lumi (double intlumi)
	~HFDarkening ()

Static Public Attributes
static const unsigned int	lowZLimit = 1115
static const unsigned int	numberOfRLayers = 13
static const unsigned int	numberOfZLayers = 33
static const unsigned int	upperZLimit = 1280

Private Attributes
double	HFDoseLayerDarkeningPars [numberOfZLayers][numberOfRLayers]

Static Private Attributes
static const unsigned int	_numberOfRLayers = numberOfRLayers
static const unsigned int	_numberOfZLayers = numberOfZLayers

Detailed Description

Definition at line 15 of file HFDarkening.h.

Constructor & Destructor Documentation

HFDarkening::HFDarkening

(

const edm::ParameterSet &

pset

)

Definition at line 17 of file HFDarkening.cc.

References edm::ParameterSet::getParameter(), mps_fire::i, dttmaxenums::R, and DOFs::Z.

                                                    {
  //HF area of consideration is 1115 cm from interaction point to 1280cm in z-axis
  //Radius (cm) - 13 cm from Beam pipe to 130cm (the top of HF active area)
  //Dose in MRad
  
  vecOfDoubles HFDosePars = pset.getParameter<vecOfDoubles>("doseLayerDepth");
  int i = 0;
  for (int Z = 0; Z != _numberOfZLayers; ++Z) {
    for (int R = 0; R != _numberOfRLayers; ++R) {
      HFDoseLayerDarkeningPars[Z][R] = HFDosePars[i];
      ++i;
    }
  }
}

HFDarkening::~HFDarkening

(

)

Definition at line 32 of file HFDarkening.cc.

{ }

Member Function Documentation

double HFDarkening::degradation

(

double

mrad

)

Definition at line 62 of file HFDarkening.cc.

References JetChargeProducer_cfi::exp, and funct::pow().

                                           {
  return (exp(-1.44*pow(mrad/100,0.44)*0.2/4.343));
}

double HFDarkening::dose	(	unsigned int	layer,
		double	radius
	)

Definition at line 34 of file HFDarkening.cc.

References TCMET_cfi::radius.

                                                          {
  // Radii are 13-17, 17-20, 20-24, 24-29, 29-34, 34-41, 41-48, 48-58, 58-69, 69-82, 82-98, 98-116, 116-130
  // These radii are specific to the geometry of the dose map, which closely matches HF Tower Geometry,
  // but not exactly.
  if (layer > (_numberOfZLayers-1))
  {
    return 0.;
  }
  
  int radius = 0;
  if (Radius <= 17.0) radius = 0;
  else if (Radius <= 20.0) radius = 1;
  else if (Radius <= 24.0) radius = 2;
  else if (Radius <= 29.0) radius = 3;
  else if (Radius <= 34.0) radius = 4;
  else if (Radius <= 41.0) radius = 5;
  else if (Radius <= 48.0) radius = 6;
  else if (Radius <= 58.0) radius = 7;
  else if (Radius <= 69.0) radius = 8;
  else if (Radius <= 82.0) radius = 9;
  else if (Radius <= 98.0) radius = 10;
  else if (Radius <= 116.0) radius = 11;
  else if (Radius <= 130.0) radius = 12;
  else return 0.;
  
  return HFDoseLayerDarkeningPars[layer][radius];
}

double HFDarkening::int_lumi

(

double

intlumi

)

Definition at line 66 of file HFDarkening.cc.

                                           {
  return (intlumi/500.);
}

Member Data Documentation

const unsigned int HFDarkening::_numberOfRLayers = numberOfRLayers

staticprivate

Definition at line 35 of file HFDarkening.h.

const unsigned int HFDarkening::_numberOfZLayers = numberOfZLayers

staticprivate

Definition at line 34 of file HFDarkening.h.

double HFDarkening::HFDoseLayerDarkeningPars[numberOfZLayers][numberOfRLayers]

private

Definition at line 33 of file HFDarkening.h.

const unsigned int HFDarkening::lowZLimit = 1115

static

Definition at line 29 of file HFDarkening.h.

Referenced by HCalSD::getFromLibrary().

const unsigned int HFDarkening::numberOfRLayers = 13

static

Definition at line 27 of file HFDarkening.h.

const unsigned int HFDarkening::numberOfZLayers = 33

static

Definition at line 26 of file HFDarkening.h.

Referenced by HCalSD::getFromLibrary().

const unsigned int HFDarkening::upperZLimit = 1280

static

Definition at line 30 of file HFDarkening.h.

Referenced by HCalSD::getFromLibrary().