HFShowerPMT Class Reference

#include <HFShowerPMT.h>

Public Member Functions
double	getHits (const G4Step *aStep)
double	getRadius ()
	HFShowerPMT (const std::string &name, const HcalDDDSimConstants *hcons, const HcalSimulationParameters *hps, edm::ParameterSet const &p)
virtual	~HFShowerPMT ()

Private Attributes
std::unique_ptr< HFCherenkov >	cherenkov_
const HcalDDDSimConstants *	hcalConstant_
const HcalSimulationParameters *	hcalsimpar_
int	indexF
int	indexR
double	pePerGeV
std::vector< int >	pmtFib1
std::vector< int >	pmtFib2
std::vector< int >	pmtR1
std::vector< int >	pmtR2
std::vector< double >	rTable

Detailed Description

Definition at line 19 of file HFShowerPMT.h.

Constructor & Destructor Documentation

HFShowerPMT()

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

Definition at line 19 of file HFShowerPMT.cc.

References cherenkov_, edm::ParameterSet::getParameter(), HcalDDDSimConstants::getRTableHF(), hcalConstant_, hcalsimpar_, cuy::ii, AlCaHLTBitMon_ParallelJobs::p, pePerGeV, pmtFib1, pmtFib2, HcalSimulationParameters::pmtFiberLeft_, HcalSimulationParameters::pmtFiberRight_, HcalSimulationParameters::pmtLeft_, pmtR1, pmtR2, HcalSimulationParameters::pmtRight_, and rTable.

    : hcalConstant_(hcons), hcalsimpar_(hps) {
  edm::ParameterSet m_HF = p.getParameter<edm::ParameterSet>("HFShowerPMT");
  pePerGeV = m_HF.getParameter<double>("PEPerGeVPMT");

  //Special Geometry parameters
  pmtR1 = hcalsimpar_->pmtRight_;
  pmtFib1 = hcalsimpar_->pmtFiberRight_;
  pmtR2 = hcalsimpar_->pmtLeft_;
  pmtFib2 = hcalsimpar_->pmtFiberLeft_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerPMT: gets the Index matches for " << pmtR1.size() << " PMTs";
  for (unsigned int ii = 0; ii < pmtR1.size(); ii++) {
    edm::LogVerbatim("HFShower") << "HFShowerPMT: rIndexR[" << ii << "] = " << pmtR1[ii] << " fibreR[" << ii
                                 << "] = " << pmtFib1[ii] << " rIndexL[" << ii << "] = " << pmtR2[ii] << " fibreL["
                                 << ii << "] = " << pmtFib2[ii];
  }
#endif
  cherenkov_ = std::make_unique<HFCherenkov>(m_HF);

  // Special Geometry parameters
  rTable = hcalConstant_->getRTableHF();
#ifdef EDM_ML_DEBUG
  std::stringstream sss;
  for (unsigned int ig = 0; ig < rTable.size(); ++ig) {
    if (ig / 10 * 10 == ig) {
      sss << "\n";
    }
    sss << "  " << rTable[ig] / CLHEP::cm;
  }
  edm::LogVerbatim("HFShowerPMT") << "HFShowerPMT: " << rTable.size() << " rTable(cm):" << sss.str();
#endif
}

~HFShowerPMT()

HFShowerPMT::~HFShowerPMT ( )

virtual

Definition at line 56 of file HFShowerPMT.cc.

{}

Member Function Documentation

getHits()

double HFShowerPMT::getHits

(

const G4Step *

aStep

)

Definition at line 58 of file HFShowerPMT.cc.

References HLT_2024v13_cff::beta, cherenkov_, indexF, indexR, pePerGeV, BPHMonitor_cfi::photons, pmtFib1, pmtFib2, pmtR1, and pmtR2.

                                               {
  indexR = indexF = -1;

  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  const G4VTouchable* touch = preStepPoint->GetTouchable();
  int boxNo = touch->GetReplicaNumber(2);
  int pmtNo = touch->GetReplicaNumber(1);
  if (boxNo <= 1) {
    indexR = pmtR1[pmtNo - 1];
    indexF = pmtFib1[pmtNo - 1];
  } else {
    indexR = pmtR2[pmtNo - 1];
    indexF = pmtFib2[pmtNo - 1];
  }

#ifdef EDM_ML_DEBUG
  double edep = aStep->GetTotalEnergyDeposit();
  edm::LogVerbatim("HFShower") << "HFShowerPMT: Box " << boxNo << " PMT " << pmtNo << " Mapped Indices " << indexR
                               << ", " << indexF << " Edeposit " << edep / CLHEP::MeV << " MeV; PE "
                               << edep * pePerGeV / CLHEP::GeV;
#endif

  double photons = 0;
  if (indexR >= 0 && indexF > 0) {
    G4Track* aTrack = aStep->GetTrack();
    G4ParticleDefinition* particleDef = aTrack->GetDefinition();
    double stepl = aStep->GetStepLength();
    double beta = preStepPoint->GetBeta();
    G4ThreeVector pDir = aTrack->GetDynamicParticle()->GetMomentumDirection();
    G4ThreeVector localMom = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformAxis(pDir);
    photons = cherenkov_->computeNPEinPMT(particleDef, beta, localMom.x(), localMom.y(), localMom.z(), stepl);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HFShower") << "HFShowerPMT::getHits: for particle " << particleDef->GetParticleName() << " Step "
                                 << stepl << " Beta " << beta << " Direction " << pDir << " Local " << localMom
                                 << " p.e. " << photons;
#endif
  }
  return photons;
}

getRadius()

double HFShowerPMT::getRadius

(

)

Definition at line 98 of file HFShowerPMT.cc.

References indexF, indexR, alignCSCRings::r, and rTable.

                              {
  double r = 0.;
  if (indexR >= 0 && indexR + 1 < (int)(rTable.size()))
    r = 0.5 * (rTable[indexR] + rTable[indexR + 1]);
#ifdef EDM_ML_DEBUG
  else
    edm::LogVerbatim("HFShower") << "HFShowerPMT::getRadius: R " << indexR << " F " << indexF;
#endif
  if (indexF == 2)
    r = -r;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFShower") << "HFShowerPMT: Radius (" << indexR << "/" << indexF << ") " << r;
#endif
  return r;
}

Member Data Documentation

cherenkov_

std::unique_ptr<HFCherenkov> HFShowerPMT::cherenkov_

private

Definition at line 32 of file HFShowerPMT.h.

Referenced by getHits(), and HFShowerPMT().

hcalConstant_

const HcalDDDSimConstants* HFShowerPMT::hcalConstant_

private

Definition at line 30 of file HFShowerPMT.h.

Referenced by HFShowerPMT().

hcalsimpar_

const HcalSimulationParameters* HFShowerPMT::hcalsimpar_

private

Definition at line 31 of file HFShowerPMT.h.

Referenced by HFShowerPMT().

indexF

int HFShowerPMT::indexF

private

Definition at line 34 of file HFShowerPMT.h.

Referenced by getHits(), and getRadius().

indexR

int HFShowerPMT::indexR

private

Definition at line 34 of file HFShowerPMT.h.

Referenced by getHits(), and getRadius().

pePerGeV

double HFShowerPMT::pePerGeV

private

Definition at line 33 of file HFShowerPMT.h.

Referenced by getHits(), and HFShowerPMT().

pmtFib1

std::vector<int> HFShowerPMT::pmtFib1

private

Definition at line 36 of file HFShowerPMT.h.

Referenced by getHits(), and HFShowerPMT().

pmtFib2

std::vector<int> HFShowerPMT::pmtFib2

private

Definition at line 37 of file HFShowerPMT.h.

Referenced by getHits(), and HFShowerPMT().

pmtR1

std::vector<int> HFShowerPMT::pmtR1

private

Definition at line 36 of file HFShowerPMT.h.

Referenced by getHits(), and HFShowerPMT().

pmtR2

std::vector<int> HFShowerPMT::pmtR2

private

Definition at line 37 of file HFShowerPMT.h.

Referenced by getHits(), and HFShowerPMT().

rTable

std::vector<double> HFShowerPMT::rTable

private

Definition at line 35 of file HFShowerPMT.h.

Referenced by getRadius(), and HFShowerPMT().