#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 (uint32_t &id)

void	recoHFShowerLibrary (const FSimTrack &myTrack)

	~FastHFShowerLibrary ()

Private Attributes
bool	applyFidCut

const edm::ParameterSet	fast

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 42 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 fast, and hitMap.

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

53 { return hitMap; };

FastHFShowerLibrary::hitMap

std::map< CaloHitID, float > hitMap

Definition: FastHFShowerLibrary.h:63

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

Definition at line 48 of file FastHFShowerLibrary.cc.

References fast, edm::EventSetup::get(), hcalConstants, hfshower, g4SimHits_cfi::HFShowerLibrary, initializeOnce, name, numberingFromDDD, 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 = (HcalDDDSimConstants*)(&(*hdc));
 
   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 ( uint32_t & id )

Definition at line 131 of file FastHFShowerLibrary.cc.

References HcalDetId::depth(), stringResolutionProvider_cfi::eta, hcalConstants, HcalForward, HcalDetId::ieta(), HcalDetId::iphi(), HcalDDDSimConstants::maxHFDepth(), phi, DetId::rawId(), and HcalDetId::subdet().

Referenced by ~FastHFShowerLibrary().

                                                   {
 
   HcalDetId hid(id);
   if (hid.subdet() == HcalForward) {
     int eta = hid.ieta();
     int phi = hid.iphi();
     if (hcalConstants->maxHFDepth(eta,phi) > 2) {
       if (hid.depth() <= 2) {
     int dep = (G4UniformRand() > 0.5) ? (2+hid.depth()) : hid.depth();
     id = HcalDetId(HcalForward,eta,phi,dep).rawId();
       }
     }
   }
 }

void FastHFShowerLibrary::recoHFShowerLibrary ( const FSimTrack & myTrack )

Definition at line 75 of file FastHFShowerLibrary.cc.

References applyFidCut, particleFlowClusterECALTimeSelected_cfi::depth, HcalNumberingScheme::getUnitID(), hfshower, hitMap, hfClusterShapes_cfi::hits, mps_fire::i, FSimTrack::id(), numberingFromDDD, numberingScheme, convertSQLiteXML::ok, FSimTrack::onVFcal(), HFFibreFiducial::PMTNumber(), ntuplemaker::time, tmp, CoreSimTrack::type(), RawParticle::vertex(), FSimTrack::vfcalEntrance(), RawParticle::X(), RawParticle::Y(), and RawParticle::Z().

Referenced by CalorimetryManager::HDShowerSimulation(), CalorimetryManager::reconstruct(), 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, pos, depth, lay);
       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
 
 }