HFFibre Class Reference

#include <HFFibre.h>

Public Member Functions
double	attLength (double lambda)
	HFFibre (const std::string &name, const HcalDDDSimConstants *hcons, const HcalSimulationParameters *hps, edm::ParameterSet const &p)
double	tShift (const G4ThreeVector &point, int depth, int fromEndAbs=0)
double	zShift (const G4ThreeVector &point, int depth, int fromEndAbs=0)
	~HFFibre ()=default

Private Attributes
std::vector< double >	attL
double	cFibre
std::vector< double >	gpar
const HcalDDDSimConstants *	hcalConstant_
const HcalSimulationParameters *	hcalsimpar_
double	lambLim [2]
std::vector< double >	longFL
int	nBinAtt
int	nBinR
std::vector< double >	radius
std::vector< double >	shortFL

Detailed Description

Definition at line 18 of file HFFibre.h.

Constructor & Destructor Documentation

HFFibre::HFFibre	(	const std::string &	name,
		const HcalDDDSimConstants *	hcons,
		const HcalSimulationParameters *	hps,
		edm::ParameterSet const &	p
	)

Definition at line 16 of file HFFibre.cc.

References HcalSimulationParameters::attenuationLength_, attL, cFibre, HcalDDDSimConstants::getGparHF(), edm::ParameterSet::getParameter(), HcalDDDSimConstants::getRTableHF(), gpar, hcalConstant_, hcalsimpar_, mps_fire::i, HcalSimulationParameters::lambdaLimits_, lambLim, HcalSimulationParameters::longFiberLength_, longFL, visualization-live-secondInstance_cfg::m, nBinAtt, nBinR, radius, HcalSimulationParameters::shortFiberLength_, and shortFL.

    : hcalConstant_(hcons), hcalsimpar_(hps) {
  edm::ParameterSet m_HF =
      (p.getParameter<edm::ParameterSet>("HFShower")).getParameter<edm::ParameterSet>("HFShowerBlock");
  cFibre = c_light * (m_HF.getParameter<double>("CFibre"));
  edm::LogVerbatim("HFShower") << "HFFibre:: Speed of light in fibre " << cFibre / (CLHEP::m / CLHEP::ns) << " m/ns";
  // Attenuation length
  attL = hcalsimpar_->attenuationLength_;
  nBinAtt = static_cast<int>(attL.size());
  std::stringstream ss1;
  for (int it = 0; it < nBinAtt; it++) {
    if (it / 10 * 10 == it) {
      ss1 << "\n";
    }
    ss1 << "  " << attL[it] * CLHEP::cm;
  }
  edm::LogVerbatim("HFShower") << "HFFibre: " << nBinAtt << " attL(1/cm): " << ss1.str();
  // Limits on Lambda
  std::vector<int> nvec = hcalsimpar_->lambdaLimits_;
  lambLim[0] = nvec[0];
  lambLim[1] = nvec[1];
  edm::LogVerbatim("HFShower") << "HFFibre: Limits on lambda " << lambLim[0] << " and " << lambLim[1];
  // Fibre Lengths
  longFL = hcalsimpar_->longFiberLength_;
  std::stringstream ss2;
  for (unsigned int it = 0; it < longFL.size(); it++) {
    if (it / 10 * 10 == it) {
      ss2 << "\n";
    }
    ss2 << "  " << longFL[it] / CLHEP::cm;
  }
  edm::LogVerbatim("HFShower") << "HFFibre: " << longFL.size() << " Long Fibre Length(cm):" << ss2.str();
  shortFL = hcalsimpar_->shortFiberLength_;
  std::stringstream ss3;
  for (unsigned int it = 0; it < shortFL.size(); it++) {
    if (it / 10 * 10 == it) {
      ss3 << "\n";
    }
    ss3 << "  " << shortFL[it] / CLHEP::cm;
  }
  edm::LogVerbatim("HFShower") << "HFFibre: " << shortFL.size() << " Short Fibre Length(cm):" << ss3.str();

  // Now geometry parameters
  gpar = hcalConstant_->getGparHF();
  radius = hcalConstant_->getRTableHF();

  nBinR = static_cast<int>(radius.size());
  std::stringstream sss;
  for (int i = 0; i < nBinR; ++i) {
    if (i / 10 * 10 == i) {
      sss << "\n";
    }
    sss << "  " << radius[i] / CLHEP::cm;
  }
  edm::LogVerbatim("HFShower") << "HFFibre: " << radius.size() << " rTable(cm):" << sss.str();
}

HFFibre::~HFFibre ( )

default

Member Function Documentation

double HFFibre::attLength

(

double

lambda

)

Definition at line 76 of file HFFibre.cc.

References attL, mps_fire::i, dqmiolumiharvest::j, lambLim, and nBinAtt.

                                       {
  int i = int(nBinAtt * (lambda - lambLim[0]) / (lambLim[1] - lambLim[0]));

  int j = i;
  if (i >= nBinAtt)
    j = nBinAtt - 1;
  else if (i < 0)
    j = 0;
  double att = attL[j];
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFFibre::attLength for Lambda " << lambda << " index " << i << " " << j
                               << " Att. Length " << att;
#endif
  return att;
}

double HFFibre::tShift	(	const G4ThreeVector &	point,
		int	depth,
		int	fromEndAbs = 0
	)

Definition at line 92 of file HFFibre.cc.

References cFibre, and zShift().

                                                                            {
  double zFibre = zShift(point, depth, fromEndAbs);
  double time = zFibre / cFibre;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFFibre::tShift for point " << point << " ( depth = " << depth
                               << ", traversed length = " << zFibre / CLHEP::cm << " cm) = " << time / CLHEP::ns
                               << " ns";
#endif
  return time;
}

double HFFibre::zShift	(	const G4ThreeVector &	point,
		int	depth,
		int	fromEndAbs = 0
	)

Definition at line 103 of file HFFibre.cc.

References gpar, mps_fire::i, longFL, nBinR, radius, shortFL, and mathSSE::sqrt().

Referenced by tShift().

                                                                            {  // point is z-local

  double zFibre = 0;
  int ieta = 0;
  double length = 250 * CLHEP::cm;
  double hR = sqrt((point.x()) * (point.x()) + (point.y()) * (point.y()));

  // Defines the Radius bin by radial subdivision
  if (fromEndAbs >= 0) {
    for (int i = nBinR - 1; i > 0; --i)
      if (hR < radius[i])
        ieta = nBinR - i - 1;
  }

  // Defines the full length of the fibre
  if (depth == 2) {
    if (static_cast<int>(shortFL.size()) > ieta)
      length = shortFL[ieta] + gpar[0];
  } else {
    if (static_cast<int>(longFL.size()) > ieta)
      length = longFL[ieta];
  }
  zFibre = length;  // from beginning of abs (full length)

  if (fromEndAbs > 0) {
    zFibre -= gpar[1];  // length from end of HF
  } else {
    double zz = 0.5 * gpar[1] + point.z();  // depth of point of photon emission (from beginning of HF)
    zFibre -= zz;                           // length of fiber from point of photon emission
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFFibre::zShift for point " << point << " (R = " << hR / CLHEP::cm
                               << " cm, Index = " << ieta << ", depth = " << depth
                               << ", Fibre Length = " << length / CLHEP::cm << " cm = " << zFibre / CLHEP::cm << " cm)";
#endif
  return zFibre;
}

Member Data Documentation

std::vector<double> HFFibre::attL

private

Definition at line 37 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

double HFFibre::cFibre

private

Definition at line 34 of file HFFibre.h.

Referenced by HFFibre(), and tShift().

std::vector<double> HFFibre::gpar

private

Definition at line 35 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

const HcalDDDSimConstants* HFFibre::hcalConstant_

private

Definition at line 32 of file HFFibre.h.

Referenced by HFFibre().

const HcalSimulationParameters* HFFibre::hcalsimpar_

private

Definition at line 33 of file HFFibre.h.

Referenced by HFFibre().

double HFFibre::lambLim[2]

private

Definition at line 39 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

std::vector<double> HFFibre::longFL

private

Definition at line 36 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

int HFFibre::nBinAtt

private

Definition at line 38 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

int HFFibre::nBinR

private

Definition at line 38 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

std::vector<double> HFFibre::radius

private

Definition at line 35 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

std::vector<double> HFFibre::shortFL

private

Definition at line 36 of file HFFibre.h.

Referenced by HFFibre(), and zShift().