FCS_StepInfoSD Class Reference

Common sensitive detector class for LAr systems. More...

#include <FCS_StepInfoSD.h>

Inheritance diagram for FCS_StepInfoSD:

Collaboration diagram for FCS_StepInfoSD:

Public Member Functions
	FCS_StepInfoSD (G4String a_name, const FCS_Param::Config &config)
	Constructor. More...
virtual G4bool	ProcessHits (G4Step *a_step, G4TouchableHistory *) override
	Main processing method. More...
void	EndOfAthenaEvent (ISF_FCS_Parametrization::FCS_StepInfoCollection *hitContnainer)
	End of athena event processing. More...
void	setupHelpers (const LArEM_ID *EM, const LArFCAL_ID *FCAL, const LArHEC_ID *HEC, const TileID *tile)
	Sets the ID helper pointers. More...

Protected Member Functions
void	getCaloDDManager ()
	Keep a map instead of trying to keep the full vector. More...
void	update_map (const CLHEP::Hep3Vector &l_vec, const Identifier &l_identifier, double l_energy, double l_time, bool l_valid, int l_detector)

Protected Attributes
FCS_Param::Config	m_config
const LArEM_ID *	m_larEmID {nullptr}
	Pointers to the identifier helpers. More...
const LArFCAL_ID *	m_larFcalID {nullptr}
const LArHEC_ID *	m_larHecID {nullptr}
const TileID *	m_tileID {nullptr}
CxxUtils::CachedPointer< const CaloDetDescrManager >	m_calo_dd_man
std::map< Identifier, std::vector< ISF_FCS_Parametrization::FCS_StepInfo * > * >	m_hit_map

Private Member Functions
double	getMaxTime (const CaloCell_ID::CaloSample &layer) const

Detailed Description

Common sensitive detector class for LAr systems.

This SD implementation saves the standard LArHits. See LArG4CalibSD for an SD that handles calibration hits.

Definition at line 71 of file FCS_StepInfoSD.h.

Constructor & Destructor Documentation

FCS_StepInfoSD()

FCS_StepInfoSD::FCS_StepInfoSD	(	G4String	a_name,
		const FCS_Param::Config &	config
	)

Constructor.

Definition at line 21 of file FCS_StepInfoSD.cxx.

    : G4VSensitiveDetector(std::move(a_name)),
      m_config(config),
      m_calo_dd_man(nullptr) {}

Member Function Documentation

EndOfAthenaEvent()

void FCS_StepInfoSD::EndOfAthenaEvent

(

ISF_FCS_Parametrization::FCS_StepInfoCollection *

hitContnainer

)

End of athena event processing.

Definition at line 146 of file FCS_StepInfoSD.cxx.

                                                                 {
  // Unpack map into vector
  for (auto it : m_hit_map) {
    for (auto* a_s : *it.second) {
      // Giving away ownership of the objects!
      hitContainer->push_back(a_s);
    }
    it.second->clear();
    delete it.second;
  }  // Vector of IDs in the map
  m_hit_map.clear();
  if (m_config.verboseLevel > 4) {
    G4cout << this->GetName()
           << " DEBUG EndOfAthenaEvent: After initial cleanup, N="
           << hitContainer->size() << G4endl;
  }
  return;
}

getCaloDDManager()

void FCS_StepInfoSD::getCaloDDManager ( )

protected

Keep a map instead of trying to keep the full vector.

At the end of the event we'll push the map back into the FCS_StepInfoCollection in StoreGate.

Definition at line 49 of file FCS_StepInfoSD.cxx.

                                      {
  SG::ReadCondHandleKey<CaloDetDescrManager> caloMgrKey{"CaloDetDescrManager"};
  if (caloMgrKey.initialize().isFailure()) {
    G4ExceptionDescription description;
    description << "Failed to get CaloDetDescrManager!";
    G4Exception("FCS_StepInfoSD", "FCSBadCall", FatalException, description);
    abort();
  }
  SG::ReadCondHandle<CaloDetDescrManager> caloMgr(
      caloMgrKey, Gaudi::Hive::currentContext());
  m_calo_dd_man.set(*caloMgr);
}

getMaxTime()

double FCS_StepInfoSD::getMaxTime ( const CaloCell_ID::CaloSample &  layer ) const

inlineprivate

NB The result of this function should actually be constant for each SD

Definition at line 34 of file FCS_StepInfoSD.cxx.

                                              {
  if (layer >= CaloCell_ID::PreSamplerB && layer <= CaloCell_ID::EME3) {
    return m_config.m_maxTimeLAr;
  }
  if (layer >= CaloCell_ID::HEC0 && layer <= CaloCell_ID::HEC3) {
    return m_config.m_maxTimeHEC;
  }
  if (layer >= CaloCell_ID::FCAL0 && layer <= CaloCell_ID::FCAL2) {
    return m_config.m_maxTimeFCAL;
  }
  return m_config.m_maxTime;
}

ProcessHits()

G4bool FCS_StepInfoSD::ProcessHits ( G4Step *  a_step, G4TouchableHistory *    )

overridevirtual

Main processing method.

Reimplemented in LArFCS_StepInfoSD, and TileFCS_StepInfoSD.

Definition at line 26 of file FCS_StepInfoSD.cxx.

                                                               {
  G4ExceptionDescription description;
  description << "ProcessHits: Base class method should not be called!!!";
  G4Exception("FCS_StepInfoSD", "FCSBadCall", FatalException, description);
  abort();
  return false;
}

setupHelpers()

void FCS_StepInfoSD::setupHelpers ( const LArEM_ID *  EM, const LArFCAL_ID *  FCAL, const LArHEC_ID *  HEC, const TileID *  tile  )

inline

Sets the ID helper pointers.

Definition at line 84 of file FCS_StepInfoSD.h.

                                                              {
    m_larEmID = EM;
    m_larFcalID = FCAL;
    m_larHecID = HEC;
    m_tileID = tile;
  }

update_map()

void FCS_StepInfoSD::update_map ( const CLHEP::Hep3Vector &  l_vec, const Identifier &  l_identifier, double  l_energy, double  l_time, bool  l_valid, int  l_detector  )

protected

Definition at line 62 of file FCS_StepInfoSD.cxx.

{
  // NB l_identifier refers to:
  // - the cell identifier for LAr
  // - the PMT identifier for Tile
 
  // Drop any hits that don't have a good identifier attached
  if (!m_calo_dd_man.get()->get_element(l_identifier)) {
    if (m_config.verboseLevel > 4) {
      G4cout << this->GetName() << " DEBUG update_map: bad identifier: "
             << l_identifier.getString() << " skipping this hit." << G4endl;
    }
    return;
  }
 
  auto map_item = m_hit_map.find(l_identifier);
  if (map_item == m_hit_map.end()) {
    m_hit_map[l_identifier] =
        new std::vector<ISF_FCS_Parametrization::FCS_StepInfo*>;
    m_hit_map[l_identifier]->reserve(200);
    m_hit_map[l_identifier]->push_back(
        new ISF_FCS_Parametrization::FCS_StepInfo(l_vec, l_identifier, l_energy,
                                                  l_time, l_valid, l_detector));
  } else {
 
    // Get the appropriate merging limits
    const CaloCell_ID::CaloSample& layer =
        m_calo_dd_man.get()->get_element(l_identifier)->getSampling();
 
    double timeWindow = m_config.m_maxTime;
    const double distWinLong = m_config.m_maxRadiusLongitudinal.at(layer);
    const double distWinLat = m_config.m_maxRadiusLateral.at(layer);
 
    const double tsame(this->getMaxTime(layer));
    bool match = false;
    for (auto* map_it : *map_item->second) {
      // Time check ... both a global flag and a check on the layer
      const double delta_t = std::fabs(map_it->time() - l_time);
      if (delta_t >= tsame) {
        continue;
      }
      if (delta_t >= timeWindow) {
        continue;
      }
 
      // Distance check
      const CLHEP::Hep3Vector& currentPosition = map_it->position();
      const double currentPosition_mag = currentPosition.mag();
      const double proj_longitudinal =
          currentPosition.dot(l_vec) / currentPosition_mag;
      const double delta_longitudinal = currentPosition_mag - proj_longitudinal;
      if (std::fabs(delta_longitudinal) >= distWinLong) {
        continue;
      }
 
      // Lateral distance check
      double delta_lateral_2 = l_vec.mag2() - proj_longitudinal * proj_longitudinal;
      if (delta_lateral_2 < 0) {
        delta_lateral_2 = 0;  // Avoid negative square root
      }
      const double delta_lateral =
          std::sqrt(delta_lateral_2);
      if (delta_lateral >= distWinLat) {
        continue;
      }
 
      // Found a match.  Make a temporary that will be deleted!
      const ISF_FCS_Parametrization::FCS_StepInfo my_info(
          l_vec, l_identifier, l_energy, l_time, l_valid, l_detector);
      *map_it += my_info;
      match = true;
      break;
    }  // End of search for match in time and space
    if (!match) {
      map_item->second->push_back(new ISF_FCS_Parametrization::FCS_StepInfo(
          l_vec, l_identifier, l_energy, l_time, l_valid, l_detector));
    }  // Didn't match
  }  // ID already in the map
  return;
}  // That's it for updating the map!

Member Data Documentation

m_calo_dd_man

CxxUtils::CachedPointer<const CaloDetDescrManager> FCS_StepInfoSD::m_calo_dd_man

protected

Definition at line 106 of file FCS_StepInfoSD.h.

m_config

FCS_Param::Config FCS_StepInfoSD::m_config

protected

Definition at line 100 of file FCS_StepInfoSD.h.

m_hit_map

std::map<Identifier, std::vector<ISF_FCS_Parametrization::FCS_StepInfo*>*> FCS_StepInfoSD::m_hit_map

protected

Definition at line 108 of file FCS_StepInfoSD.h.

m_larEmID

const LArEM_ID* FCS_StepInfoSD::m_larEmID {nullptr}

protected

Pointers to the identifier helpers.

Definition at line 102 of file FCS_StepInfoSD.h.

m_larFcalID

const LArFCAL_ID* FCS_StepInfoSD::m_larFcalID {nullptr}

protected

Definition at line 103 of file FCS_StepInfoSD.h.

m_larHecID

const LArHEC_ID* FCS_StepInfoSD::m_larHecID {nullptr}

protected

Definition at line 104 of file FCS_StepInfoSD.h.

m_tileID

const TileID* FCS_StepInfoSD::m_tileID {nullptr}

protected

Definition at line 105 of file FCS_StepInfoSD.h.

---

The documentation for this class was generated from the following files:

• FCS_StepInfoSD.h
• FCS_StepInfoSD.cxx