#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	(	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, 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] / cm;
   }
   edm::LogVerbatim("HFShowerPMT") << "HFShowerPMT: " << rTable.size() << " rTable(cm):" << sss.str();
 #endif
 }

HFShowerPMT::~HFShowerPMT ( )

virtual

Definition at line 56 of file HFShowerPMT.cc.

56 {}

Member Function Documentation

double HFShowerPMT::getHits ( const G4Step * aStep )

Definition at line 58 of file HFShowerPMT.cc.

References HLT_FULL_cff::beta, cherenkov_, GeV, indexF, indexR, MeV, pePerGeV, interactiveExample::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 / MeV << " MeV; PE " << edep * pePerGeV / 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;
 }

double HFShowerPMT::getRadius ( )

Definition at line 97 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;
 }