HFShowerLibrary Class Reference

#include <HFShowerLibrary.h>

Classes
struct	BranchCache
struct	BranchReader
struct	FileParams
struct	Hit
struct	Params
struct	VersionInfo

Public Member Functions
std::vector< Hit >	fillHits (const G4ThreeVector &p, const G4ThreeVector &v, int parCode, double parEnergy, bool &ok, double weight, double time, bool onlyLong=false)
std::vector< Hit >	getHits (const G4Step *aStep, bool &ok, double weight, bool onlyLong=false)
	HFShowerLibrary (const std::string &name, const HcalDDDSimConstants *hcons, const HcalSimulationParameters *hps, edm::ParameterSet const &p)
	HFShowerLibrary (Params const &, FileParams const &, HFFibre::Params)
	~HFShowerLibrary ()

Private Types
enum	FileFormat { FileFormat::kOld, FileFormat::kNew, FileFormat::kNewV3 }

Private Member Functions
HFShowerPhotonCollection	extrapolate (int, double)
HFShowerPhotonCollection	getRecord (int, int) const
	HFShowerLibrary (const HcalDDDSimConstants *hcons, const HcalSimulationParameters *hps, edm::ParameterSet const &hfShower, edm::ParameterSet const &hfShowerLibrary)
HFShowerPhotonCollection	interpolate (int, double)
VersionInfo	loadEventInfo (TBranch *, int fileVersion)
bool	rInside (double r) const
void	storePhoton (HFShowerPhoton const &iPhoton, HFShowerPhotonCollection &iPhotons) const

Private Attributes
bool	applyFidCut_
double	backProb_
double	dphi_
BranchReader	emBranch_
bool	equalizeTimeShift_
int	evtPerBin_
HFFibre	fibre_
std::vector< double >	gpar_
BranchReader	hadBranch_
std::unique_ptr< TFile >	hf_
int	nMomBin_
std::vector< double >	pmom_
double	probMax_
double	rMax_
double	rMin_
int	totEvents_
bool	verbose_

Detailed Description

Definition at line 28 of file HFShowerLibrary.h.

Member Enumeration Documentation

FileFormat

enum HFShowerLibrary::FileFormat

strongprivate

Enumerator
kOld
kNew
kNewV3

Definition at line 90 of file HFShowerLibrary.h.

{ kOld, kNew, kNewV3 };

Constructor & Destructor Documentation

HFShowerLibrary() [1/3]

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

Definition at line 77 of file HFShowerLibrary.cc.

    : HFShowerLibrary(
          hcons,
          hps,
          p.getParameter<edm::ParameterSet>("HFShower").getParameter<edm::ParameterSet>("HFShowerBlock"),
          p.getParameter<edm::ParameterSet>("HFShowerLibrary").getParameter<edm::ParameterSet>("HFLibraryFileBlock")) {}

~HFShowerLibrary()

HFShowerLibrary::~HFShowerLibrary

(

)

Definition at line 197 of file HFShowerLibrary.cc.

References hf_.

                                  {
  if (hf_)
    hf_->Close();
}

HFShowerLibrary() [2/3]

HFShowerLibrary::HFShowerLibrary	(	Params const &	iParams,
		FileParams const &	iFileParams,
		HFFibre::Params	iFibreParams
	)

Definition at line 87 of file HFShowerLibrary.cc.

References cms::cuda::assert(), backProb_, HFShowerLibrary::FileParams::branchEvInfo_, HFShowerLibrary::FileParams::cacheBranches_, visDQMUpload::context, dphi_, emBranch_, HFShowerLibrary::FileParams::emBranchName_, equalizeTimeShift_, edmPickEvents::event, evtPerBin_, Exception, DQMFileSaver_cfi::fileFormat, HFShowerLibrary::FileParams::fileName_, HFShowerLibrary::FileParams::fileVersion_, hadBranch_, HFShowerLibrary::FileParams::hadBranchName_, hf_, mps_fire::i, info(), kNew, kNewV3, kOld, HFShowerLibrary::VersionInfo::libVers_, HFShowerLibrary::VersionInfo::listVersion_, loadEventInfo(), nMomBin_, hltrates_dqm_sourceclient-live_cfg::offset, pmom_, probMax_, rMax_, rMin_, AlCaHLTBitMon_QueryRunRegistry::string, totEvents_, and verbose_.

    : fibre_(iFibreParams),
      hf_(),
      emBranch_(),
      hadBranch_(),
      verbose_(iParams.verbose_),
      applyFidCut_(iParams.applyFidCut_),
      equalizeTimeShift_(iParams.equalizeTimeShift_),
      probMax_(iParams.probMax_),
      backProb_(iParams.backProb_),
      dphi_(iParams.dphi_),
      rMin_(iFibreParams.rTableHF_.front()),
      rMax_(iFibreParams.rTableHF_.back()),
      gpar_(iFibreParams.gParHF_) {
  std::string pTreeName = iFileParams.fileName_;

  const char* nTree = pTreeName.c_str();
  {
    //It is not safe to open a Tfile in one thread and close in another without adding the following:
    TDirectory::TContext context;
    hf_ = std::unique_ptr<TFile>(TFile::Open(nTree));
  }
  if (!hf_->IsOpen()) {
    edm::LogError("HFShower") << "HFShowerLibrary: opening " << nTree << " failed";
    throw cms::Exception("Unknown", "HFShowerLibrary") << "Opening of " << pTreeName << " fails\n";
  } else {
    edm::LogVerbatim("HFShower") << "HFShowerLibrary: opening " << nTree << " successfully";
  }

  auto fileFormat = FileFormat::kOld;
  const int fileVersion = iFileParams.fileVersion_;

  auto newForm = iFileParams.branchEvInfo_.empty();
  TTree* event(nullptr);
  if (newForm) {
    fileFormat = FileFormat::kNew;
    event = (TTree*)hf_->Get("HFSimHits");
  } else {
    event = (TTree*)hf_->Get("Events");
  }
  VersionInfo versionInfo;
  if (event) {
    TBranch* evtInfo(nullptr);
    if (!newForm) {
      std::string info = iFileParams.branchEvInfo_;
      evtInfo = event->GetBranch(info.c_str());
    }
    if (evtInfo || newForm) {
      versionInfo = loadEventInfo(evtInfo, fileVersion);
    } else {
      edm::LogError("HFShower") << "HFShowerLibrary: HFShowerLibrayEventInfo"
                                << " Branch does not exist in Event";
      throw cms::Exception("Unknown", "HFShowerLibrary") << "Event information absent\n";
    }
  } else {
    edm::LogError("HFShower") << "HFShowerLibrary: Events Tree does not "
                              << "exist";
    throw cms::Exception("Unknown", "HFShowerLibrary") << "Events tree absent\n";
  }

  edm::LogVerbatim("HFShower").log([&](auto& logger) {
    logger << "HFShowerLibrary: Library " << versionInfo.libVers_ << " ListVersion " << versionInfo.listVersion_
           << " File version " << fileVersion << " Events Total " << totEvents_ << " and " << evtPerBin_
           << " per bin\n";
    logger << "HFShowerLibrary: Energies (GeV) with " << nMomBin_ << " bins\n";
    for (int i = 0; i < nMomBin_; ++i) {
      if (i / 10 * 10 == i && i > 0) {
        logger << "\n";
      }
      logger << "  " << pmom_[i] / CLHEP::GeV;
    }
  });

  std::string nameBr = iFileParams.emBranchName_;
  auto emBranch = event->GetBranch(nameBr.c_str());
  if (verbose_)
    emBranch->Print();
  nameBr = iFileParams.hadBranchName_;
  auto hadBranch = event->GetBranch(nameBr.c_str());
  if (verbose_)
    hadBranch->Print();

  if (emBranch->GetClassName() == std::string("vector<float>")) {
    assert(fileFormat == FileFormat::kNew);
    fileFormat = FileFormat::kNewV3;
  }
  emBranch_ = BranchReader(emBranch, fileFormat, 0, iFileParams.cacheBranches_ ? totEvents_ : 0);
  size_t offset = 0;
  if ((fileFormat == FileFormat::kNewV3 && fileVersion < 3) || (fileFormat == FileFormat::kNew && fileVersion < 2)) {
    //NOTE: for this format, the hadBranch is all empty up to
    // totEvents_ (which is more like 1/2*GenEntries())
    offset = totEvents_;
  }
  hadBranch_ = BranchReader(hadBranch, fileFormat, offset, iFileParams.cacheBranches_ ? totEvents_ : 0);

  edm::LogVerbatim("HFShower") << " HFShowerLibrary:Branch " << iFileParams.emBranchName_ << " has "
                               << emBranch->GetEntries() << " entries and Branch " << iFileParams.hadBranchName_
                               << " has " << hadBranch->GetEntries()
                               << " entries\n HFShowerLibrary::No packing information - Assume x, y, z are not in "
                                  "packed form\n Maximum probability cut off "
                               << probMax_ << "  Back propagation of light probability " << backProb_
                               << " Flag for equalizing Time Shift for different eta " << equalizeTimeShift_;

  edm::LogVerbatim("HFShower") << "HFShowerLibrary: rMIN " << rMin_ / CLHEP::cm << " cm and rMax " << rMax_ / CLHEP::cm
                               << " (Half) Phi Width of wedge " << dphi_ / CLHEP::deg;
  if (iFileParams.cacheBranches_) {
    hf_.reset();
  }
}

HFShowerLibrary() [3/3]

HFShowerLibrary::HFShowerLibrary ( const HcalDDDSimConstants *  hcons, const HcalSimulationParameters *  hps, edm::ParameterSet const &  hfShower, edm::ParameterSet const &  hfShowerLibrary  )

private

Definition at line 69 of file HFShowerLibrary.cc.

    : HFShowerLibrary(paramsFrom(hfShower, hfShowerLibrary, hcons->getPhiTableHF().front()),
                      fileParamsFrom(hfShowerLibrary),
                      HFFibre::Params(hfShower.getParameter<double>("CFibre"), hcons, hps)) {}

Member Function Documentation

extrapolate()

HFShowerPhotonCollection HFShowerLibrary::extrapolate ( int  type, double  pin  )

private

Definition at line 656 of file HFShowerLibrary.cc.

References evtPerBin_, getRecord(), createfilelist::int, dqmiolumiharvest::j, nMomBin_, l1ctLayer1_cff::nPhoton, displacedMuons_cfi::photon, BPHMonitor_cfi::photons, pmom_, alignCSCRings::r, storePhoton(), totEvents_, and w().

Referenced by fillHits().

                                                                          {
  int nrec = int(pin / pmom_[nMomBin_ - 1]);
  double w = (pin - pmom_[nMomBin_ - 1] * nrec) / pmom_[nMomBin_ - 1];
  nrec++;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerLibrary:: Extrapolate for Energy " << pin / CLHEP::GeV << " GeV with "
                               << nMomBin_ << " momentum bins and " << evtPerBin_ << " entries/bin -- "
                               << "total " << totEvents_ << " using " << nrec << " records";
#endif
  std::vector<int> irc(nrec);

  for (int ir = 0; ir < nrec; ir++) {
    double r = G4UniformRand();
    irc[ir] = int(evtPerBin_ * 0.5 * r) + (nMomBin_ - 1) * evtPerBin_ + 1;
    if (irc[ir] < 1) {
      edm::LogWarning("HFShower") << "HFShowerLibrary:: Illegal irc[" << ir << "] = " << irc[ir] << " now set to 1";
      irc[ir] = 1;
    } else if (irc[ir] > totEvents_) {
      edm::LogWarning("HFShower") << "HFShowerLibrary:: Illegal irc[" << ir << "] = " << irc[ir] << " now set to "
                                  << totEvents_;
      irc[ir] = totEvents_;
#ifdef EDM_ML_DEBUG
    } else {
      edm::LogVerbatim("HFShower") << "HFShowerLibrary::Extrapolation use irc[" << ir << "] = " << irc[ir];
#endif
    }
  }

  HFShowerPhotonCollection pe;
  std::size_t npold = 0;
  for (int ir = 0; ir < nrec; ir++) {
    if (irc[ir] > 0) {
      auto const photons = getRecord(type, irc[ir]);
      int nPhoton = photons.size();
      npold += nPhoton;
      pe.reserve(pe.size() + nPhoton);
      for (auto const& photon : photons) {
        double r = G4UniformRand();
        if (ir != nrec - 1 || r < w) {
          storePhoton(photon, pe);
        }
      }
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HFShower") << "HFShowerLibrary: Record [" << ir << "] = " << irc[ir] << " npold = " << npold;
#endif
    }
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerLibrary:: uses " << npold << " photons";
#endif

  if (pe.size() > npold || npold == 0)
    edm::LogWarning("HFShower") << "HFShowerLibrary: Extrapolation Warning == " << nrec << " records " << irc[0] << ", "
                                << irc[1] << ", ... gives a buffer of " << npold << " photons and fills " << pe.size()
                                << " *****";
#ifdef EDM_ML_DEBUG
  else
    edm::LogVerbatim("HFShower") << "HFShowerLibrary: Extrapolation == " << nrec << " records " << irc[0] << ", "
                                 << irc[1] << ", ... gives a buffer of " << npold << " photons and fills " << pe.size()
                                 << " PE";
  for (std::size_t j = 0; j < pe.size(); j++)
    edm::LogVerbatim("HFShower") << "Photon " << j << " " << pe[j];
#endif
  return pe;
}

fillHits()

std::vector< HFShowerLibrary::Hit > HFShowerLibrary::fillHits	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v,
		int	parCode,
		double	parEnergy,
		bool &	ok,
		double	weight,
		double	time,
		bool	onlyLong = false
	)

Definition at line 243 of file HFShowerLibrary.cc.

References funct::abs(), applyFidCut_, HFFibre::attLength(), backProb_, funct::cos(), hcalRecHitTable_cff::depth, HFShowerLibrary::Hit::depth, dphi_, equalizeTimeShift_, JetChargeProducer_cfi::exp, extrapolate(), fibre_, gpar_, mps_fire::i, createfilelist::int, interpolate(), G4TrackToParticleID::isGammaElectronPositron(), G4TrackToParticleID::isStableHadron(), nMomBin_, convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, displacedMuons_cfi::photon, pmom_, HFShowerLibrary::Hit::position, probMax_, alignCSCRings::r, diffTwoXMLs::r1, diffTwoXMLs::r2, rInside(), funct::sin(), remoteMonitoring_LASER_era2018_cfg::threshold, HFShowerLibrary::Hit::time, funct::true, HFFibre::tShift(), geometryCSVtoXML::xx, geometryCSVtoXML::yy, hit::z, HFFibre::zShift(), gpuVertexFinder::zv, and geometryCSVtoXML::zz.

Referenced by getHits().

                                                                           {
  std::vector<HFShowerLibrary::Hit> hit;
  ok = false;
  bool isEM = G4TrackToParticleID::isGammaElectronPositron(parCode);
  // shower is built only for gamma, e+- and stable hadrons
  if (!isEM && !G4TrackToParticleID::isStableHadron(parCode)) {
    return hit;
  }
  ok = true;

  // remove low-energy component
  const double threshold = 50 * MeV;
  if (pin < threshold) {
    return hit;
  }

  double pz = momDir.z();
  double zint = hitPoint.z();

  // if particle moves from interaction point or "backwards (halo)
  bool backward = (pz * zint < 0.) ? true : false;

  double sphi = sin(momDir.phi());
  double cphi = cos(momDir.phi());
  double ctheta = cos(momDir.theta());
  double stheta = sin(momDir.theta());

  HFShowerPhotonCollection pe;
  if (isEM) {
    if (pin < pmom_[nMomBin_ - 1]) {
      pe = interpolate(0, pin);
    } else {
      pe = extrapolate(0, pin);
    }
  } else {
    if (pin < pmom_[nMomBin_ - 1]) {
      pe = interpolate(1, pin);
    } else {
      pe = extrapolate(1, pin);
    }
  }

  std::size_t nHit = 0;
  HFShowerLibrary::Hit oneHit;
#ifdef EDM_ML_DEBUG
  int i = 0;
#endif
  for (auto const& photon : pe) {
    double zv = std::abs(photon.z());  // abs local z
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HFShower") << "HFShowerLibrary: Hit " << i++ << " " << photon << " zv " << zv;
#endif
    if (zv <= gpar_[1] && photon.lambda() > 0 && (photon.z() >= 0 || (zv > gpar_[0] && (!onlyLong)))) {
      int depth = 1;
      if (onlyLong) {
      } else if (!backward) {  // fully valid only for "front" particles
        if (photon.z() < 0)
          depth = 2;                 // with "front"-simulated shower lib.
      } else {                       // for "backward" particles - almost equal
        double r = G4UniformRand();  // share between L and S fibers
        if (r > 0.5)
          depth = 2;
      }

      // Updated coordinate transformation from local
      //  back to global using two Euler angles: phi and theta
      double pex = photon.x();
      double pey = photon.y();

      double xx = pex * ctheta * cphi - pey * sphi + zv * stheta * cphi;
      double yy = pex * ctheta * sphi + pey * cphi + zv * stheta * sphi;
      double zz = -pex * stheta + zv * ctheta;

      G4ThreeVector pos = hitPoint + G4ThreeVector(xx, yy, zz);
      zv = std::abs(pos.z()) - gpar_[4] - 0.5 * gpar_[1];
      G4ThreeVector lpos = G4ThreeVector(pos.x(), pos.y(), zv);

      zv = fibre_.zShift(lpos, depth, 0);  // distance to PMT !

      double r = pos.perp();
      double p = fibre_.attLength(photon.lambda());
      double fi = pos.phi();
      if (fi < 0)
        fi += CLHEP::twopi;
      int isect = int(fi / dphi_) + 1;
      isect = (isect + 1) / 2;
      double dfi = ((isect * 2 - 1) * dphi_ - fi);
      if (dfi < 0)
        dfi = -dfi;
      double dfir = r * sin(dfi);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HFShower") << "HFShowerLibrary: Position shift " << xx << ", " << yy << ", " << zz << ": "
                                   << pos << " R " << r << " Phi " << fi << " Section " << isect << " R*Dfi " << dfir
                                   << " Dist " << zv;
#endif
      zz = std::abs(pos.z());
      double r1 = G4UniformRand();
      double r2 = G4UniformRand();
      double r3 = backward ? G4UniformRand() : -9999.;
      if (!applyFidCut_)
        dfir += gpar_[5];

#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HFShower") << "HFShowerLibrary: rLimits " << rInside(r) << " attenuation " << r1 << ":"
                                   << exp(-p * zv) << " r2 " << r2 << " r3 " << r3 << " rDfi " << gpar_[5] << " zz "
                                   << zz << " zLim " << gpar_[4] << ":" << gpar_[4] + gpar_[1] << "\n"
                                   << "  rInside(r) :" << rInside(r) << "  r1 <= exp(-p*zv) :" << (r1 <= exp(-p * zv))
                                   << "  r2 <= probMax :" << (r2 <= probMax_ * weight)
                                   << "  r3 <= backProb :" << (r3 <= backProb_)
                                   << "  dfir > gpar[5] :" << (dfir > gpar_[5])
                                   << "  zz >= gpar[4] :" << (zz >= gpar_[4])
                                   << "  zz <= gpar[4]+gpar[1] :" << (zz <= gpar_[4] + gpar_[1]);
#endif
      if (rInside(r) && r1 <= exp(-p * zv) && r2 <= probMax_ * weight && dfir > gpar_[5] && zz >= gpar_[4] &&
          zz <= gpar_[4] + gpar_[1] && r3 <= backProb_ && (depth != 2 || zz >= gpar_[4] + gpar_[0])) {
        double tdiff = (equalizeTimeShift_) ? (fibre_.tShift(lpos, depth, -1)) : (fibre_.tShift(lpos, depth, 1));
        oneHit.position = pos;
        oneHit.depth = depth;
        oneHit.time = (tSlice + (photon.t()) + tdiff);
        hit.push_back(oneHit);

#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HFShower") << "HFShowerLibrary: Final Hit " << nHit << " position " << (hit[nHit].position)
                                     << " Depth " << (hit[nHit].depth) << " Time " << tSlice << ":" << photon.t() << ":"
                                     << tdiff << ":" << (hit[nHit].time);
#endif
        ++nHit;
      }
#ifdef EDM_ML_DEBUG
      else
        edm::LogVerbatim("HFShower") << "HFShowerLibrary: REJECTED !!!";
#endif
      if (onlyLong && zz >= gpar_[4] + gpar_[0] && zz <= gpar_[4] + gpar_[1]) {
        r1 = G4UniformRand();
        r2 = G4UniformRand();
        if (rInside(r) && r1 <= exp(-p * zv) && r2 <= probMax_ && dfir > gpar_[5]) {
          double tdiff = (equalizeTimeShift_) ? (fibre_.tShift(lpos, 2, -1)) : (fibre_.tShift(lpos, 2, 1));
          oneHit.position = pos;
          oneHit.depth = 2;
          oneHit.time = (tSlice + (photon.t()) + tdiff);
          hit.push_back(oneHit);
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("HFShower") << "HFShowerLibrary: Final Hit " << nHit << " position " << (hit[nHit].position)
                                       << " Depth " << (hit[nHit].depth) << " Time " << tSlice << ":" << photon.t()
                                       << ":" << tdiff << ":" << (hit[nHit].time);
#endif
          ++nHit;
        }
      }
    }
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerLibrary: Total Hits " << nHit << " out of " << pe.size() << " PE";
#endif
  if (nHit > pe.size() && !onlyLong) {
    edm::LogWarning("HFShower") << "HFShowerLibrary: Hit buffer " << pe.size() << " smaller than " << nHit << " Hits";
  }
  return hit;
}

getHits()

std::vector< HFShowerLibrary::Hit > HFShowerLibrary::getHits	(	const G4Step *	aStep,
		bool &	ok,
		double	weight,
		bool	onlyLong = false
	)

Definition at line 202 of file HFShowerLibrary.cc.

References funct::cos(), fillHits(), gpar_, funct::sin(), and HLT_2023v12_cff::track.

                                                                          {
  auto const preStepPoint = aStep->GetPreStepPoint();
  auto const postStepPoint = aStep->GetPostStepPoint();
  auto const track = aStep->GetTrack();
  // Get Z-direction
  auto const aParticle = track->GetDynamicParticle();
  const G4ThreeVector& momDir = aParticle->GetMomentumDirection();
  const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
  int parCode = track->GetDefinition()->GetPDGEncoding();

  // VI: for ions use internally pdg code of alpha in order to keep
  // consistency with previous simulation
  if (track->GetDefinition()->IsGeneralIon()) {
    parCode = 1000020040;
  }

#ifdef EDM_ML_DEBUG
  G4String partType = track->GetDefinition()->GetParticleName();
  const G4ThreeVector localPos = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
  double zoff = localPos.z() + 0.5 * gpar_[1];

  edm::LogVerbatim("HFShower") << "HFShowerLibrary::getHits " << partType << " of energy "
                               << track->GetKineticEnergy() / CLHEP::GeV << " GeV weight= " << weight
                               << " onlyLong: " << onlyLong << "  dir.orts " << momDir.x() << ", " << momDir.y() << ", "
                               << momDir.z() << "  Pos x,y,z = " << hitPoint.x() << "," << hitPoint.y() << ","
                               << hitPoint.z() << " (" << zoff << ")   sphi,cphi,stheta,ctheta  = " << sin(momDir.phi())
                               << "," << cos(momDir.phi()) << ", " << sin(momDir.theta()) << "," << cos(momDir.theta());
#endif

  double tSlice = (postStepPoint->GetGlobalTime()) / CLHEP::nanosecond;

  // use kinetic energy for protons and ions
  double pin = (track->GetDefinition()->GetBaryonNumber() > 0) ? preStepPoint->GetKineticEnergy()
                                                               : preStepPoint->GetTotalEnergy();

  return fillHits(hitPoint, momDir, parCode, pin, isKilled, weight, tSlice, onlyLong);
}

getRecord()

HFShowerPhotonCollection HFShowerLibrary::getRecord ( int  type, int  record  ) const

private

Definition at line 527 of file HFShowerLibrary.cc.

References emBranch_, HFShowerLibrary::BranchReader::getRecord(), hadBranch_, dqmiolumiharvest::j, l1ctLayer1_cff::nPhoton, and AlCaHarvesting_cff::record.

Referenced by extrapolate(), and interpolate().

                                                                              {
  HFShowerPhotonCollection photo;
  if (type > 0) {
    photo = hadBranch_.getRecord(record);
  } else {
    photo = emBranch_.getRecord(record);
  }
#ifdef EDM_ML_DEBUG
  int nPhoton = photo.size();
  edm::LogVerbatim("HFShower") << "HFShowerLibrary::getRecord: Record " << record << " of type " << type << " with "
                               << nPhoton << " photons";
  for (int j = 0; j < nPhoton; j++)
    edm::LogVerbatim("HFShower") << "Photon " << j << " " << photo[j];
#endif
  return photo;
}

interpolate()

HFShowerPhotonCollection HFShowerLibrary::interpolate ( int  type, double  pin  )

private

Definition at line 575 of file HFShowerLibrary.cc.

References evtPerBin_, getRecord(), createfilelist::int, dqmiolumiharvest::j, nMomBin_, l1ctLayer1_cff::nPhoton, displacedMuons_cfi::photon, BPHMonitor_cfi::photons, pmom_, alignCSCRings::r, storePhoton(), totEvents_, and w().

Referenced by fillHits().

                                                                          {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerLibrary:: Interpolate for Energy " << pin / CLHEP::GeV << " GeV with "
                               << nMomBin_ << " momentum bins and " << evtPerBin_ << " entries/bin -- total "
                               << totEvents_;
#endif
  int irc[2] = {0, 0};
  double w = 0.;
  double r = G4UniformRand();

  if (pin < pmom_[0]) {
    w = pin / pmom_[0];
    irc[1] = int(evtPerBin_ * r) + 1;
    irc[0] = 0;
  } else {
    for (int j = 0; j < nMomBin_ - 1; j++) {
      if (pin >= pmom_[j] && pin < pmom_[j + 1]) {
        w = (pin - pmom_[j]) / (pmom_[j + 1] - pmom_[j]);
        if (j == nMomBin_ - 2) {
          irc[1] = int(evtPerBin_ * 0.5 * r);
        } else {
          irc[1] = int(evtPerBin_ * r);
        }
        irc[1] += (j + 1) * evtPerBin_ + 1;
        r = G4UniformRand();
        irc[0] = int(evtPerBin_ * r) + 1 + j * evtPerBin_;
        if (irc[0] < 0) {
          edm::LogWarning("HFShower") << "HFShowerLibrary:: Illegal irc[0] = " << irc[0] << " now set to 0";
          irc[0] = 0;
        } else if (irc[0] > totEvents_) {
          edm::LogWarning("HFShower") << "HFShowerLibrary:: Illegal irc[0] = " << irc[0] << " now set to "
                                      << totEvents_;
          irc[0] = totEvents_;
        }
      }
    }
  }
  if (irc[1] < 1) {
    edm::LogWarning("HFShower") << "HFShowerLibrary:: Illegal irc[1] = " << irc[1] << " now set to 1";
    irc[1] = 1;
  } else if (irc[1] > totEvents_) {
    edm::LogWarning("HFShower") << "HFShowerLibrary:: Illegal irc[1] = " << irc[1] << " now set to " << totEvents_;
    irc[1] = totEvents_;
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerLibrary:: Select records " << irc[0] << " and " << irc[1] << " with weights "
                               << 1 - w << " and " << w;
#endif
  HFShowerPhotonCollection pe;

  std::size_t npold = 0;
  for (int ir = 0; ir < 2; ir++) {
    if (irc[ir] > 0) {
      auto photons = getRecord(type, irc[ir]);
      int nPhoton = photons.size();
      npold += nPhoton;
      pe.reserve(pe.size() + nPhoton);
      for (auto const& photon : photons) {
        r = G4UniformRand();
        if ((ir == 0 && r > w) || (ir > 0 && r < w)) {
          storePhoton(photon, pe);
        }
      }
    }
  }

  if ((pe.size() > npold || (npold == 0 && irc[0] > 0)) && !(pe.empty() && npold == 0))
    edm::LogWarning("HFShower") << "HFShowerLibrary: Interpolation Warning =="
                                << " records " << irc[0] << " and " << irc[1] << " gives a buffer of " << npold
                                << " photons and fills " << pe.size() << " *****";
#ifdef EDM_ML_DEBUG
  else
    edm::LogVerbatim("HFShower") << "HFShowerLibrary: Interpolation == records " << irc[0] << " and " << irc[1]
                                 << " gives a buffer of " << npold << " photons and fills " << pe.size() << " PE";
  for (std::size_t j = 0; j < pe.size(); j++)
    edm::LogVerbatim("HFShower") << "Photon " << j << " " << pe[j];
#endif
  return pe;
}

loadEventInfo()

HFShowerLibrary::VersionInfo HFShowerLibrary::loadEventInfo ( TBranch *  branch, int  fileVersion  )

private

Definition at line 544 of file HFShowerLibrary.cc.

References MicroEventContent_cff::branch, evtPerBin_, mps_fire::i, HFShowerLibrary::VersionInfo::libVers_, HFShowerLibrary::VersionInfo::listVersion_, nMomBin_, pmom_, and totEvents_.

Referenced by HFShowerLibrary().

                                                                                          {
  VersionInfo versionInfo;
  if (branch) {
    std::vector<HFShowerLibraryEventInfo> eventInfoCollection;
    branch->SetAddress(&eventInfoCollection);
    branch->GetEntry(0);
    edm::LogVerbatim("HFShower") << "HFShowerLibrary::loadEventInfo loads EventInfo Collection of size "
                                 << eventInfoCollection.size() << " records";

    totEvents_ = eventInfoCollection[0].totalEvents();
    nMomBin_ = eventInfoCollection[0].numberOfBins();
    evtPerBin_ = eventInfoCollection[0].eventsPerBin();
    versionInfo.libVers_ = eventInfoCollection[0].showerLibraryVersion();
    versionInfo.listVersion_ = eventInfoCollection[0].physListVersion();
    pmom_ = eventInfoCollection[0].energyBins();
  } else {
    edm::LogVerbatim("HFShower") << "HFShowerLibrary::loadEventInfo loads EventInfo from hardwired"
                                 << " numbers";

    nMomBin_ = 16;
    evtPerBin_ = (fileVersion == 0) ? 5000 : 10000;
    totEvents_ = nMomBin_ * evtPerBin_;
    versionInfo.libVers_ = (fileVersion == 0) ? 1.1 : 1.2;
    versionInfo.listVersion_ = 3.6;
    pmom_ = {2, 3, 5, 7, 10, 15, 20, 30, 50, 75, 100, 150, 250, 350, 500, 1000};
  }
  for (int i = 0; i < nMomBin_; i++)
    pmom_[i] *= CLHEP::GeV;
  return versionInfo;
}

rInside()

bool HFShowerLibrary::rInside ( double  r ) const

private

Definition at line 411 of file HFShowerLibrary.cc.

References alignCSCRings::r, rMax_, and rMin_.

Referenced by fillHits().

{ return (r >= rMin_ && r <= rMax_); }

storePhoton()

void HFShowerLibrary::storePhoton ( HFShowerPhoton const &  iPhoton, HFShowerPhotonCollection &  iPhotons  ) const

private

Definition at line 722 of file HFShowerLibrary.cc.

Referenced by extrapolate(), and interpolate().

                                                                                                         {
  iPhotons.push_back(iPhoton);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerLibrary: storePhoton " << iPhoton << " npe " << iPhotons.size() << " "
                               << iPhotons.back();
#endif
}

Member Data Documentation

applyFidCut_

bool HFShowerLibrary::applyFidCut_

private

Definition at line 139 of file HFShowerLibrary.h.

Referenced by fillHits().

backProb_

double HFShowerLibrary::backProb_

private

Definition at line 144 of file HFShowerLibrary.h.

Referenced by fillHits(), and HFShowerLibrary().

dphi_

double HFShowerLibrary::dphi_

private

Definition at line 145 of file HFShowerLibrary.h.

Referenced by fillHits(), and HFShowerLibrary().

emBranch_

BranchReader HFShowerLibrary::emBranch_

private

Definition at line 136 of file HFShowerLibrary.h.

Referenced by getRecord(), and HFShowerLibrary().

equalizeTimeShift_

bool HFShowerLibrary::equalizeTimeShift_

private

Definition at line 143 of file HFShowerLibrary.h.

Referenced by fillHits(), and HFShowerLibrary().

evtPerBin_

int HFShowerLibrary::evtPerBin_

private

Definition at line 140 of file HFShowerLibrary.h.

Referenced by extrapolate(), HFShowerLibrary(), interpolate(), and loadEventInfo().

fibre_

HFFibre HFShowerLibrary::fibre_

private

Definition at line 133 of file HFShowerLibrary.h.

Referenced by fillHits().

gpar_

std::vector<double> HFShowerLibrary::gpar_

private

Definition at line 146 of file HFShowerLibrary.h.

Referenced by fillHits(), and getHits().

hadBranch_

BranchReader HFShowerLibrary::hadBranch_

private

Definition at line 137 of file HFShowerLibrary.h.

Referenced by getRecord(), and HFShowerLibrary().

hf_

std::unique_ptr<TFile> HFShowerLibrary::hf_

private

Definition at line 134 of file HFShowerLibrary.h.

Referenced by HFShowerLibrary(), and ~HFShowerLibrary().

nMomBin_

int HFShowerLibrary::nMomBin_

private

Definition at line 140 of file HFShowerLibrary.h.

Referenced by extrapolate(), fillHits(), HFShowerLibrary(), interpolate(), and loadEventInfo().

pmom_

std::vector<double> HFShowerLibrary::pmom_

private

Definition at line 141 of file HFShowerLibrary.h.

Referenced by extrapolate(), fillHits(), HFShowerLibrary(), interpolate(), and loadEventInfo().

probMax_

double HFShowerLibrary::probMax_

private

Definition at line 144 of file HFShowerLibrary.h.

Referenced by fillHits(), and HFShowerLibrary().

rMax_

double HFShowerLibrary::rMax_

private

Definition at line 145 of file HFShowerLibrary.h.

Referenced by HFShowerLibrary(), and rInside().

rMin_

double HFShowerLibrary::rMin_

private

Definition at line 145 of file HFShowerLibrary.h.

Referenced by HFShowerLibrary(), and rInside().

totEvents_

int HFShowerLibrary::totEvents_

private

Definition at line 140 of file HFShowerLibrary.h.

Referenced by extrapolate(), HFShowerLibrary(), interpolate(), and loadEventInfo().

verbose_

bool HFShowerLibrary::verbose_

private

Definition at line 139 of file HFShowerLibrary.h.

Referenced by HFShowerLibrary().