#include <FastHFShowerLibrary.h>

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

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

void const	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 40 of file FastHFShowerLibrary.h.

Constructor & Destructor Documentation

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

Definition at line 43 of file FastHFShowerLibrary.cc.

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

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

FastHFShowerLibrary::~FastHFShowerLibrary ( )

inline

Definition at line 46 of file FastHFShowerLibrary.h.

References initHFShowerLibrary(), modifyDepth(), and recoHFShowerLibrary().

46 {;}

Member Function Documentation

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

inline

Definition at line 53 of file FastHFShowerLibrary.h.

References hitMap, and SetRandom().

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

53 { return hitMap; };

FastHFShowerLibrary::hitMap

std::map< CaloHitID, float > hitMap

Definition: FastHFShowerLibrary.h:65

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

Definition at line 49 of file FastHFShowerLibrary.cc.

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

Referenced by FamosManager::setupGeometryAndField(), and ~FastHFShowerLibrary().

                                                                                {
 
   edm::LogInfo("FastCalorimetry") << "initHFShowerLibrary::initialization"; 
 
   edm::ESTransientHandle<DDCompactView> cpv;
   iSetup.get<IdealGeometryRecord>().get(cpv);
 
   edm::ESHandle<HcalDDDSimConstants>    hdc;
   iSetup.get<HcalSimNumberingRecord>().get(hdc);
   hcalConstants = hdc.product();
 
   std::string name = "HcalHits";
   numberingFromDDD.reset(new HcalNumberingFromDDD(hcalConstants));  
   hfshower.reset(new HFShowerLibrary(name,*cpv,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
 }

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

Definition at line 144 of file FastHFShowerLibrary.cc.

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

Referenced by recoHFShowerLibrary(), and ~FastHFShowerLibrary().

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

void FastHFShowerLibrary::recoHFShowerLibrary ( const FSimTrack & myTrack )

Definition at line 85 of file FastHFShowerLibrary.cc.

References applyFidCut, egammaForCoreTracking_cff::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, tmp, CoreSimTrack::type(), RawParticle::Vect(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

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

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

void FastHFShowerLibrary::SetRandom ( const RandomEngineAndDistribution * rnd )

Definition at line 76 of file FastHFShowerLibrary.cc.

References LogDebug, and RandomEngineAndDistribution::theEngine().

Referenced by getHitsMap(), CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::initialize().

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