#include <FastHFShowerLibrary.h>

Public Member Functions
	FastHFShowerLibrary (edm::ParameterSet const &p)

const std::map< CaloHitID, float > &	getHitsMap ()

const void	initHFShowerLibrary (const edm::EventSetup &)

void	modifyDepth (HcalNumberingFromDDD::HcalID &id)

void	recoHFShowerLibrary (const FSimTrack &myTrack)

void	SetRandom (const RandomEngineAndDistribution *)

	~FastHFShowerLibrary ()

Private Attributes
bool	applyFidCut

const edm::ParameterSet	fast

const HcalDDDSimConstants *	hcalConstants

std::unique_ptr< HFShowerLibrary >	hfshower

std::map< CaloHitID, float >	hitMap

std::string	name

std::unique_ptr< HcalNumberingFromDDD >	numberingFromDDD

HcalNumberingScheme	numberingScheme

Detailed Description

Definition at line 38 of file FastHFShowerLibrary.h.

Constructor & Destructor Documentation

◆ FastHFShowerLibrary()

FastHFShowerLibrary::FastHFShowerLibrary ( edm::ParameterSet const & p )

Definition at line 42 of file FastHFShowerLibrary.cc.

                                                                  : fast(p) {
   edm::ParameterSet m_HS = p.getParameter<edm::ParameterSet>("HFShowerLibrary");
   applyFidCut = m_HS.getParameter<bool>("ApplyFiducialCut");
 }

References applyFidCut, edm::ParameterSet::getParameter(), and AlCaHLTBitMon_ParallelJobs::p.

◆ ~FastHFShowerLibrary()

FastHFShowerLibrary::~FastHFShowerLibrary ( )

inline

Definition at line 42 of file FastHFShowerLibrary.h.

42 { ; }

Member Function Documentation

◆ getHitsMap()

const std::map<CaloHitID, float>& FastHFShowerLibrary::getHitsMap ( )

inline

Definition at line 48 of file FastHFShowerLibrary.h.

48 { return hitMap; };

References hitMap.

Referenced by CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::reconstructTrack().

◆ initHFShowerLibrary()

const void FastHFShowerLibrary::initHFShowerLibrary ( const edm::EventSetup & iSetup )

Definition at line 47 of file FastHFShowerLibrary.cc.

                                                                                {
   edm::LogInfo("FastCalorimetry") << "initHFShowerLibrary::initialization";
  
   edm::ESHandle<HcalDDDSimConstants> hdc;
   iSetup.get<HcalSimNumberingRecord>().get(hdc);
   hcalConstants = hdc.product();
  
   edm::ESHandle<HcalDDDSimulationConstants> hdsc;
   iSetup.get<HcalSimNumberingRecord>().get(hdsc);
   const HcalDDDSimulationConstants* hsps = hdsc.product();
  
   std::string name = "HcalHits";
   numberingFromDDD.reset(new HcalNumberingFromDDD(hcalConstants));
   hfshower.reset(new HFShowerLibrary(name, hcalConstants, hsps->hcalsimpar(), fast));
  
   //only one thread can be allowed to setup the G4 physics table.
   std::call_once(initializeOnce, []() {
     // Geant4 particles
     G4DecayPhysics decays;
     decays.ConstructParticle();
     G4ParticleTable* partTable = G4ParticleTable::GetParticleTable();
     partTable->SetReadiness();
   });
   //G4ParticleTable* partTable = G4ParticleTable::GetParticleTable();
   //hfshower->initRun(partTable, hcalConstants);  // init particle code
 }

References bJpsiMuMuTrigSettings_cff::decays, fast, edm::EventSetup::get(), get, hcalConstants, HcalDDDSimulationConstants::hcalsimpar(), hfshower, Calorimetry_cff::HFShowerLibrary, initializeOnce, name, numberingFromDDD, edm::ESHandle< T >::product(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by FamosManager::setupGeometryAndField().

◆ modifyDepth()

void FastHFShowerLibrary::modifyDepth ( HcalNumberingFromDDD::HcalID & id )

Definition at line 136 of file FastHFShowerLibrary.cc.

                                                                     {
   if (id.subdet == HcalForward) {
     int ieta = (id.zside == 0) ? -id.etaR : id.etaR;
     if (hcalConstants->maxHFDepth(ieta, id.phis) > 2) {
       if (id.depth <= 2) {
         if (G4UniformRand() > 0.5)
           id.depth += 2;
       }
     }
   }
 }

References LEDCalibrationChannels::depth, hcalConstants, HcalForward, LEDCalibrationChannels::ieta, and HcalDDDSimConstants::maxHFDepth().

Referenced by recoHFShowerLibrary().

◆ recoHFShowerLibrary()

void FastHFShowerLibrary::recoHFShowerLibrary ( const FSimTrack & myTrack )

Definition at line 81 of file FastHFShowerLibrary.cc.

                                                                       {
 #ifdef DebugLog
   edm::LogInfo("FastCalorimetry") << "FastHFShowerLibrary: recoHFShowerLibrary ";
 #endif
  
   if (!myTrack.onVFcal()) {
 #ifdef DebugLog
     edm::LogInfo("FastCalorimetry") << "FastHFShowerLibrary: we should not be here ";
 #endif
   }
  
   hitMap.clear();
   double eGen = 1000. * myTrack.vfcalEntrance().e();  // energy in [MeV]
   double delZv = (myTrack.vfcalEntrance().vertex().Z() > 0.0) ? 50.0 : -50.0;
   G4ThreeVector vertex(10. * myTrack.vfcalEntrance().vertex().X(),
                        10. * myTrack.vfcalEntrance().vertex().Y(),
                        10. * myTrack.vfcalEntrance().vertex().Z() + delZv);  // in [mm]
  
   G4ThreeVector direction(
       myTrack.vfcalEntrance().Vect().X(), myTrack.vfcalEntrance().Vect().Y(), myTrack.vfcalEntrance().Vect().Z());
  
   bool ok;
   double weight = 1.0;  // rad. damage
   int parCode = myTrack.type();
   double tSlice = 0.1 * vertex.mag() / 29.98;
  
   std::vector<HFShowerLibrary::Hit> hits =
       hfshower->fillHits(vertex, direction, parCode, eGen, ok, weight, tSlice, false);
  
   for (unsigned int i = 0; i < hits.size(); ++i) {
     G4ThreeVector pos = hits[i].position;
     int depth = hits[i].depth;
     double time = hits[i].time;
     if (!applyFidCut || (HFFibreFiducial::PMTNumber(pos) > 0)) {
       //    if (!applyFidCut || (applyFidCut && HFFibreFiducial::PMTNumber(pos)>0)) {
       int det = 5;
       int lay = 1;
       uint32_t id = 0;
       HcalNumberingFromDDD::HcalID tmp =
           numberingFromDDD->unitID(det, math::XYZVectorD(pos.x(), pos.y(), pos.z()), depth, lay);
       modifyDepth(tmp);
       id = numberingScheme.getUnitID(tmp);
  
       CaloHitID current_id(id, time, myTrack.id());
       std::map<CaloHitID, float>::iterator cellitr;
       cellitr = hitMap.find(current_id);
       if (cellitr == hitMap.end()) {
         hitMap.insert(std::pair<CaloHitID, float>(current_id, 1.0));
       } else {
         cellitr->second += 1.0;
       }
     }  // end of isItinFidVolume check
   }    // end loop over hits
 }

References applyFidCut, LEDCalibrationChannels::depth, RawParticle::e(), HcalNumberingScheme::getUnitID(), hfshower, hitMap, hfClusterShapes_cfi::hits, mps_fire::i, FSimTrack::id(), modifyDepth(), numberingFromDDD, numberingScheme, convertSQLiteXML::ok, FSimTrack::onVFcal(), HFFibreFiducial::PMTNumber(), ntuplemaker::time, createJobs::tmp, CoreSimTrack::type(), RawParticle::Vect(), bphysicsOniaDQM_cfi::vertex, RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::reconstructTrack().

◆ SetRandom()

void FastHFShowerLibrary::SetRandom ( const RandomEngineAndDistribution * rnd )

Definition at line 74 of file FastHFShowerLibrary.cc.

                                                                           {
   // define Geant4 engine per thread
   G4Random::setTheEngine(&(rnd->theEngine()));
   LogDebug("FastHFShowerLibrary::recoHFShowerLibrary")
       << "Begin of event " << G4UniformRand() << "  " << rnd->theEngine().name() << "  " << rnd->theEngine();
 }