HFFibre Class Reference

#include <HFFibre.h>

Classes
struct	Params

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

Private Attributes
std::vector< double >	attL_
double	cFibre_
std::vector< double >	gpar_
std::array< double, 2 >	lambLim_
std::vector< double >	longFL_
int	nBinAtt_
int	nBinR_
std::vector< double >	radius_
std::vector< double >	shortFL_

Detailed Description

Definition at line 19 of file HFFibre.h.

Constructor & Destructor Documentation

HFFibre() [1/2]

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

Definition at line 26 of file HFFibre.cc.

    : HFFibre(Params(p.getParameter<edm::ParameterSet>("HFShower")
                         .getParameter<edm::ParameterSet>("HFShowerBlock")
                         .getParameter<double>("CFibre"),
                     hcons,
                     hps)) {}

HFFibre() [2/2]

HFFibre::HFFibre

(

Params

iP

)

Definition at line 33 of file HFFibre.cc.

References attL_, cFibre_, gpar_, mps_fire::i, lambLim_, longFL_, visualization-live-secondInstance_cfg::m, nBinAtt_, nBinR_, radius_, and shortFL_.

    : cFibre_(c_light * iP.fractionOfSpeedOfLightInFibre_),
      gpar_(std::move(iP.gParHF_)),
      radius_(std::move(iP.rTableHF_)),
      shortFL_(std::move(iP.shortFibreLength_)),
      longFL_(std::move(iP.longFibreLength_)),
      attL_(std::move(iP.attenuationLength_)),
      lambLim_(iP.lambdaLimits_) {
  edm::LogVerbatim("HFShower") << "HFFibre:: Speed of light in fibre " << cFibre_ / (CLHEP::m / CLHEP::ns) << " m/ns";
  // Attenuation length
  nBinAtt_ = static_cast<int>(attL_.size());

  edm::LogVerbatim("HFShower").log([&](auto& logger) {
    logger << "HFFibre: " << nBinAtt_ << " attL(1/cm): ";
    for (int it = 0; it < nBinAtt_; it++) {
      if (it / 10 * 10 == it) {
        logger << "\n";
      }
      logger << "  " << attL_[it] * CLHEP::cm;
    }
  });
  edm::LogVerbatim("HFShower") << "HFFibre: Limits on lambda " << lambLim_[0] << " and " << lambLim_[1];
  // Fibre Lengths
  edm::LogVerbatim("HFShower").log([&](auto& logger) {
    logger << "HFFibre: " << longFL_.size() << " Long Fibre Length(cm):";
    for (unsigned int it = 0; it < longFL_.size(); it++) {
      if (it / 10 * 10 == it) {
        logger << "\n";
      }
      logger << "  " << longFL_[it] / CLHEP::cm;
    }
  });
  edm::LogVerbatim("HFShower").log([&](auto& logger) {
    logger << "HFFibre: " << shortFL_.size() << " Short Fibre Length(cm):";
    for (unsigned int it = 0; it < shortFL_.size(); it++) {
      if (it / 10 * 10 == it) {
        logger << "\n";
      }
      logger << "  " << shortFL_[it] / CLHEP::cm;
    }
  });
  // Now geometry parameters

  nBinR_ = static_cast<int>(radius_.size());
  edm::LogVerbatim("HFShower").log([&](auto& logger) {
    logger << "HFFibre: " << radius_.size() << " rTable(cm):";
    for (int i = 0; i < nBinR_; ++i) {
      if (i / 10 * 10 == i) {
        logger << "\n";
      }
      logger << "  " << radius_[i] / CLHEP::cm;
    }
  });

  edm::LogVerbatim("HFShower").log([&](auto& logger) {
    logger << "HFFibre: " << gpar_.size() << " gParHF:";
    for (std::size_t i = 0; i < gpar_.size(); ++i) {
      if (i / 10 * 10 == i) {
        logger << "\n";
      }
      logger << "  " << gpar_[i];
    }
  });
}

Member Function Documentation

attLength()

double HFFibre::attLength

(

double

lambda

)

const

Definition at line 98 of file HFFibre.cc.

References attL_, mps_fire::i, createfilelist::int, dqmiolumiharvest::j, lambLim_, and nBinAtt_.

Referenced by HFShowerLibrary::fillHits(), HFShower::getHits(), and HFShowerParam::HFShowerParam().

                                             {
  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;
}

tShift()

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

Definition at line 114 of file HFFibre.cc.

References cFibre_, hcalRecHitTable_cff::depth, point, hcalRecHitTable_cff::time, and zShift().

Referenced by HFShowerLibrary::fillHits(), HFShowerParam::getHits(), and HFShower::getHits().

                                                                                  {
  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;
}

zShift()

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

Definition at line 125 of file HFFibre.cc.

References hcalRecHitTable_cff::depth, gpar_, mps_fire::i, hcalRecHitTable_cff::ieta, longFL_, nBinR_, point, radius_, shortFL_, mathSSE::sqrt(), and geometryCSVtoXML::zz.

Referenced by HFShowerLibrary::fillHits(), HFShowerParam::getHits(), and tShift().

                                                                                  {  // point is z-local

  double zFibre = 0;
  int ieta = 0;
  double length = 250 * CLHEP::cm;
  double hR = std::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

attL_

std::vector<double> HFFibre::attL_

private

Definition at line 48 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

cFibre_

double HFFibre::cFibre_

private

Definition at line 45 of file HFFibre.h.

Referenced by HFFibre(), and tShift().

gpar_

std::vector<double> HFFibre::gpar_

private

Definition at line 46 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

lambLim_

std::array<double, 2> HFFibre::lambLim_

private

Definition at line 50 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

longFL_

std::vector<double> HFFibre::longFL_

private

Definition at line 47 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

nBinAtt_

int HFFibre::nBinAtt_

private

Definition at line 49 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

nBinR_

int HFFibre::nBinR_

private

Definition at line 49 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

radius_

std::vector<double> HFFibre::radius_

private

Definition at line 46 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

shortFL_

std::vector<double> HFFibre::shortFL_

private

Definition at line 47 of file HFFibre.h.

Referenced by HFFibre(), and zShift().