HGTD_SmearedDigitizationTool Class Reference

#include <HGTD_SmearedDigitizationTool.h>

Inheritance diagram for HGTD_SmearedDigitizationTool:

Collaboration diagram for HGTD_SmearedDigitizationTool:

Public Types
using	Cluster_t = HGTD_Cluster
using	ClusterCollection_t = HGTD_ClusterCollection
using	ClusterContainer_t = HGTD_ClusterContainer
using	HGTD_DetElement_RIO_Map_t

Public Member Functions
	HGTD_SmearedDigitizationTool (const std::string &type, const std::string &name, const IInterface *parent)
	~HGTD_SmearedDigitizationTool ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual StatusCode	processAllSubEvents (const EventContext &ctx) override

Private Member Functions
StatusCode	retrieveTruth (PRD_MultiTruthCollection *&prd_truth_coll)
TimedHitCollection< SiHit >	setupTimedHitCollection ()
StatusCode	digitize (const EventContext &ctx, TimedHitCollection< SiHit > &timed_hit_collection, HGTD_DetElement_RIO_Map_t &det_element_rio_map)
StatusCode	fillMultiTruthCollection (PRD_MultiTruthCollection *, Cluster_t *, const TimedHitPtr< SiHit > &, const EventContext &)
StatusCode	fillClusterContainer (HGTD_DetElement_RIO_Map_t &det_element_rio_map, ClusterContainer_t &cluster_container)
float	smearPosition (float pos, float sig, float boundary, CLHEP::HepRandomEngine *rndmEngine)
float	smearMeanTime (float time, float time_res, CLHEP::HepRandomEngine *rndmEngine)

Private Attributes
const HGTD_DetectorManager *	m_hgtd_det_manager {nullptr}
const HGTD_ID *	m_hgtd_idhelper {nullptr}
std::string	m_prd_truth_coll_name {"PRD_MultiTruthPixel"}
ServiceHandle< PileUpMergeSvc >	m_merge_svc {this, "MergeSvc", "PileUpMergeSvc", "Merge service used in HGTD fast digitization"}
ServiceHandle< IAthRNGSvc >	m_rndm_svc {this, "RndmSvc", "AthRNGSvc", "Random number service used in HGTD fast digitization"}
	Random number service.
StringProperty	m_cluster_name {this, "HGTD_ClustersContainerName", "HGTD_Cluster", "Name of the HGTD cluster container"}
StringProperty	m_si_hit_collection_name {this, "SiHitCollectionName", "HGTD_Hits", "Name of the Si Hit collection"}
BooleanProperty	m_smear_intersection_position {this, "SmearIntersectionPosition", true, ""}
BooleanProperty	m_smear_mean_time {this, "SmearMeanTime", true, ""}
BooleanProperty	m_write_tree {this, "WriteTree", false, ""}
FloatProperty	m_pitch_x {this, "pitch_X", 1.3, ""}
FloatProperty	m_pitch_y {this, "pitch_Y", 1.3, ""}
float	m_time_res {0.035}
ServiceHandle< ITHistSvc >	m_hist_svc {this, "THistSvc", "THistSvc"}
std::unique_ptr< TFile >	m_output_file
std::unique_ptr< TTree >	m_tree
float	m_x_hit {}
float	m_y_hit {}
float	m_x_hit_smeared {}
float	m_y_hit_smeared {}
float	m_x_entry_hit {}
float	m_y_entry_hit {}
float	m_z_entry_hit {}
float	m_x_exit_hit {}
float	m_y_exit_hit {}
float	m_z_exit_hit {}
float	m_hit_time {}
float	m_hit_time_smeared {}
float	m_x_cluster_global {}
float	m_y_cluster_global {}
float	m_z_cluster_global {}
float	m_err_x_hit {}
float	m_err_y_hit {}

structors and AlgTool implementation
virtual bool	toProcess (int bunchXing) const override
	the method this base class helps implementing
virtual StatusCode	processBunchXing (int, SubEventIterator, SubEventIterator) override
virtual bool	filterPassed () const override
	dummy implementation of passing filter
virtual void	resetFilter () override
	dummy implementation of filter reset
Gaudi::Property< int >	m_firstXing
Gaudi::Property< int >	m_lastXing
Gaudi::Property< int >	m_vetoPileUpTruthLinks
bool	m_filterPassed {true}

Detailed Description

Definition at line 44 of file HGTD_SmearedDigitizationTool.h.

Member Typedef Documentation

Cluster_t

using HGTD_SmearedDigitizationTool::Cluster_t = HGTD_Cluster

Definition at line 46 of file HGTD_SmearedDigitizationTool.h.

ClusterCollection_t

using HGTD_SmearedDigitizationTool::ClusterCollection_t = HGTD_ClusterCollection

Definition at line 47 of file HGTD_SmearedDigitizationTool.h.

ClusterContainer_t

using HGTD_SmearedDigitizationTool::ClusterContainer_t = HGTD_ClusterContainer

Definition at line 48 of file HGTD_SmearedDigitizationTool.h.

HGTD_DetElement_RIO_Map_t

using HGTD_SmearedDigitizationTool::HGTD_DetElement_RIO_Map_t

Initial value:

 
std::multimap<IdentifierHash, const Cluster_t*>

Definition at line 50 of file HGTD_SmearedDigitizationTool.h.

Constructor & Destructor Documentation

HGTD_SmearedDigitizationTool()

HGTD_SmearedDigitizationTool::HGTD_SmearedDigitizationTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 26 of file HGTD_SmearedDigitizationTool.cxx.

                                                                          :
  PileUpToolBase(type, name, parent) {
}

~HGTD_SmearedDigitizationTool()

HGTD_SmearedDigitizationTool::~HGTD_SmearedDigitizationTool ( )

default

Member Function Documentation

digitize()

StatusCode HGTD_SmearedDigitizationTool::digitize ( const EventContext & ctx, TimedHitCollection< SiHit > & timed_hit_collection, HGTD_DetElement_RIO_Map_t & det_element_rio_map )

private

Definition at line 161 of file HGTD_SmearedDigitizationTool.cxx.

{
  ATH_MSG_DEBUG("--- HGTD_SmearedDigitizationTool: in digitize() ---");
 
  // Set the RNG to use for this event.
  ATHRNG::RNGWrapper* rngWrapper = m_rndm_svc->getEngine(this);
  rngWrapper->setSeed(name(), ctx);
  CLHEP::HepRandomEngine* rndmEngine = rngWrapper->getEngine(ctx);
 
  PRD_MultiTruthCollection* prd_truth_coll = nullptr;
  ATH_CHECK(retrieveTruth(prd_truth_coll));
 
  // i, e are iterator and end of one detector element
  TimedHitCollection<SiHit>::const_iterator i, e;
 
  while (timed_hit_collection.nextDetectorElement(i, e) and i != e) {
    const TimedHitPtr<SiHit>& hit(*i++);
 
    // FIXME should I cut on the deposited energy??
    // n electrons = hit->energyLoss() / electron ionisation energy
    // n * e * gain  > 4fCs
 
    int endcap = hit->getBarrelEndcap();
    int layer = hit->getLayerDisk();
    int phi_module = hit->getPhiModule();
    int eta_module = hit->getEtaModule();
 
    ATH_MSG_DEBUG("HGTD DetectorElement --> endcap "
                  << endcap << ", layer_disk " << layer << ", phi_module "
                  << phi_module << ", eta_module " << eta_module);
 
    const InDetDD::SolidStateDetectorElementBase* curr_det_element =
        m_hgtd_det_manager->getDetectorElement(endcap, layer, phi_module,
                                               eta_module);
 
    if (!curr_det_element) {
      ATH_MSG_ERROR("could not get detector element for SolidStateDetector");
      continue;
    }
 
    IdentifierHash wafer_id =
        m_hgtd_idhelper->wafer_hash(curr_det_element->identify());
 
    HepGeom::Point3D<double> hit_loc_start_pos = hit->localStartPosition();
    HepGeom::Point3D<double> hit_loc_end_pos = hit->localEndPosition();
 
    HepGeom::Point3D<double> hit_glob_start_pos =
        curr_det_element->hitLocalToLocal3D(hit_loc_start_pos);
    HepGeom::Point3D<double> hit_glob_end_pos =
        curr_det_element->hitLocalToLocal3D(hit_loc_end_pos);
 
    double globa_entry_x = hit_glob_start_pos.x();
    double globa_entry_y = hit_glob_start_pos.y();
    m_x_entry_hit = globa_entry_x;
    m_y_entry_hit = globa_entry_y;
    m_z_entry_hit = hit_glob_start_pos.z();
 
    double globa_exit_x = hit_glob_end_pos.x();
    double globa_exit_y = hit_glob_end_pos.y();
    m_x_exit_hit = globa_exit_x;
    m_y_exit_hit = globa_exit_y;
    m_z_exit_hit = hit_glob_end_pos.z();
 
    double dist_x = std::abs(std::abs(globa_exit_x) - std::abs(globa_entry_x));
    double dist_y = std::abs(std::abs(globa_exit_y) - std::abs(globa_entry_y));
 
    Amg::Vector2D local_entry(globa_entry_x, globa_entry_y);
    Amg::Vector2D local_exit(globa_exit_x, globa_exit_y);
 
    // get the identifier of the entry and the exit
    Identifier entry_id = curr_det_element->identifierOfPosition(local_entry);
    Identifier exit_id = curr_det_element->identifierOfPosition(local_exit);
 
    // now get the cellIds and check whether they're valid
    InDetDD::SiCellId entry_cell_id =
        curr_det_element->cellIdFromIdentifier(entry_id);
    InDetDD::SiCellId exit_cell_id =
        curr_det_element->cellIdFromIdentifier(exit_id);
 
    ATH_MSG_DEBUG("--- HGTD_SmearedDigitizationTool: entryId "
                  << entry_id << " --- exitId " << exit_id);
    ATH_MSG_DEBUG("--- HGTD_SmearedDigitizationTool: entryCellId "
                  << entry_cell_id << " --- exitCellId " << exit_cell_id);
 
    if (not entry_cell_id.isValid() or not exit_cell_id.isValid()) {
      continue;
    }
 
    // the intersecetion id and cellId of it. yes, this is correct
    double intersection_x = 0.5 * (globa_entry_x + globa_exit_x);
    double intersection_y = 0.5 * (globa_entry_y + globa_exit_y);
    m_x_hit = intersection_x;
    m_y_hit = intersection_y;
 
    // calculate how many pixels the particle crosses in x and y direction
    double times_x = floor(dist_x / m_pitch_x);
    double times_y = floor(dist_y / m_pitch_y);
 
    double sigma_x = m_pitch_x / std::sqrt(12);
    double sigma_y = m_pitch_y / std::sqrt(12);
 
    int element_x = times_x + 1;
    int element_y = times_y + 1;
 
    // Smear pixel in x and y direction
    // FIXME make sure the length width dimensions are correct!!!
    if (m_smear_intersection_position) {
      intersection_x = smearPosition(intersection_x, sigma_x,
                                     curr_det_element->width() / 2.,
                                     rndmEngine);
      intersection_y = smearPosition(intersection_y, sigma_y,
                                     curr_det_element->length() / 2.,
                                     rndmEngine);
    }
    m_x_hit_smeared = intersection_x;
    m_y_hit_smeared = intersection_y;
 
    Amg::Vector2D intersection(intersection_x, intersection_y);
 
    Identifier intersection_id =
        curr_det_element->identifierOfPosition(intersection);
 
    // the pixel positions and other elements for the geometrical clustering
    std::vector<Identifier> rdo_list = {intersection_id};
 
    InDetDD::SiCellId current_cell_id =
        curr_det_element->cellIdFromIdentifier(intersection_id);
 
    if (!current_cell_id.isValid()) {
      continue;
    }
 
    double length_x = element_x * m_pitch_x;
    double length_y = element_y * m_pitch_y;
 
    InDet::SiWidth si_width(Amg::Vector2D(element_x, element_y),
                            Amg::Vector2D(length_x, length_y));
 
    AmgSymMatrix(2) covariance;
    covariance.setIdentity();
    covariance(Trk::locX, Trk::locX) = sigma_x * sigma_x;
    covariance(Trk::locY, Trk::locY) = sigma_y * sigma_y;
    Amg::MatrixX cluster_err = Amg::MatrixX(covariance);
 
    float hit_time = hit->meanTime();
    m_hit_time = hit_time;
    if (m_smear_mean_time) {
      hit_time = smearMeanTime(hit_time, m_time_res, rndmEngine);
    }
    m_hit_time_smeared = hit_time;
 
    std::vector<int> tot_vec = {0}; // dummy, not used in reco currently
 
    Cluster_t* cluster = new Cluster_t(intersection_id, intersection, std::move(rdo_list),
                                       si_width, curr_det_element, std::move(cluster_err),
                                       hit_time, m_time_res, std::move(tot_vec));
 
    m_x_cluster_global = (cluster->globalPosition()).x();
    m_y_cluster_global = (cluster->globalPosition()).y();
    m_z_cluster_global = (cluster->globalPosition()).z();
    m_err_x_hit = Amg::error(cluster->localCovariance(), Trk::locX);
    m_err_y_hit = Amg::error(cluster->localCovariance(), Trk::locY);
 
    if (m_write_tree) {
      m_tree->Fill();
    }
 
    det_element_rio_map.insert(
        std::pair<IdentifierHash, const Cluster_t*>(wafer_id, cluster));
 
    ATH_CHECK(fillMultiTruthCollection(prd_truth_coll, cluster, hit, ctx));
 
  } // END while
  return StatusCode::SUCCESS;
}

fillClusterContainer()

StatusCode HGTD_SmearedDigitizationTool::fillClusterContainer ( HGTD_DetElement_RIO_Map_t & det_element_rio_map, ClusterContainer_t & cluster_container )

private

Definition at line 375 of file HGTD_SmearedDigitizationTool.cxx.

                                                                                                     {
  ATH_MSG_DEBUG(
      "--- HGTD_SmearedDigitizationTool: in fillClusterContainer() ---");
  using RIO_map_t = HGTD_DetElement_RIO_Map_t;
 
  RIO_map_t::iterator i = det_element_rio_map.begin();
  RIO_map_t::iterator e = det_element_rio_map.end();
 
  for (; i != e; i = det_element_rio_map.upper_bound(i->first)) {
 
    std::pair<RIO_map_t::iterator, RIO_map_t::iterator> range;
    range = det_element_rio_map.equal_range(i->first);
 
    RIO_map_t::iterator first_det_elem = range.first;
 
    IdentifierHash wafer_id = first_det_elem->first;
 
    const InDetDD::SolidStateDetectorElementBase* det_element =
        m_hgtd_det_manager->getDetectorElement(wafer_id);
 
    ClusterCollection_t* cluster_coll = new ClusterCollection_t(wafer_id);
    cluster_coll->setIdentifier(det_element->identify());
 
    for (RIO_map_t::iterator iter = range.first; iter != range.second; ++iter) {
 
      Cluster_t* cluster = const_cast<Cluster_t*>((*iter).second);
      cluster->setHashAndIndex(cluster_coll->identifyHash(),
                               cluster_coll->size());
      cluster_coll->push_back(cluster);
    }
 
    if (cluster_container.addCollection(cluster_coll, wafer_id)
            .isFailure()) {
      ATH_MSG_WARNING("Could not add collection to Identifiable container!");
    }
  } // end for
 
  return StatusCode::SUCCESS;
}

fillMultiTruthCollection()

StatusCode HGTD_SmearedDigitizationTool::fillMultiTruthCollection ( PRD_MultiTruthCollection * map, Cluster_t * cluster, const TimedHitPtr< SiHit > & hit, const EventContext & ctx )

private

Definition at line 359 of file HGTD_SmearedDigitizationTool.cxx.

                                                                                                                                                                           {
 
  const HepMcParticleLink trk_link = HepMcParticleLink::getRedirectedLink(hit->particleLink(), hit.eventId(), ctx); // This link should now correctly resolve to the TruthEvent McEventCollection in the main StoreGateSvc.
 
  ATH_MSG_DEBUG("Truth map filling with cluster "
                << *cluster << " and link = " << trk_link);
 
  if (trk_link.isValid()) {
    if (!HepMC::ignoreTruthLink(trk_link, m_vetoPileUpTruthLinks)) {
      map->insert(std::make_pair(cluster->identify(), trk_link));
    }
  }
 
  return StatusCode::SUCCESS;
}

filterPassed()

virtual bool PileUpToolBase::filterPassed ( ) const

inlineoverridevirtualinherited

dummy implementation of passing filter

Definition at line 49 of file PileUpToolBase.h.

{ return m_filterPassed; }

finalize()

StatusCode HGTD_SmearedDigitizationTool::finalize ( )

overridevirtual

Definition at line 75 of file HGTD_SmearedDigitizationTool.cxx.

                                                  {
 
  ATH_MSG_DEBUG("HGTD_SmearedDigitizationTool::finalize()");
 
  if (m_write_tree) {
    m_output_file->cd();
    m_tree->Write();
    m_output_file->Close();
  }
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode HGTD_SmearedDigitizationTool::initialize ( )

overridevirtual

Reimplemented from PileUpToolBase.

Definition at line 33 of file HGTD_SmearedDigitizationTool.cxx.

                                                    {
 
  ATH_MSG_DEBUG("HGTD_SmearedDigitizationTool::initialize()");
 
  // Get the HGTD Detector Manager
  ATH_CHECK(detStore()->retrieve(m_hgtd_det_manager, "HGTD"));
 
  // Get the HGTD ID helper
  ATH_CHECK(detStore()->retrieve(m_hgtd_idhelper, "HGTD_ID"));
 
  // locate the PileUpMergeSvc
  ATH_CHECK(m_merge_svc.retrieve());
 
  if (m_write_tree) {
    ATH_CHECK(m_hist_svc.retrieve());
    m_output_file = std::make_unique<TFile>("HGTD_SmearedDigiOutput.root", "RECREATE");
    m_tree = std::make_unique<TTree>("SmearedDigiTree", "SmearedDigiTree");
    m_tree->Branch("m_x_hit", &m_x_hit);
    m_tree->Branch("m_y_hit", &m_y_hit);
    m_tree->Branch("m_x_entry_hit", &m_x_entry_hit);
    m_tree->Branch("m_y_entry_hit", &m_y_entry_hit);
    m_tree->Branch("m_z_entry_hit", &m_z_entry_hit);
    m_tree->Branch("m_x_exit_hit", &m_x_exit_hit);
    m_tree->Branch("m_y_exit_hit", &m_y_exit_hit);
    m_tree->Branch("m_z_exit_hit", &m_z_exit_hit);
    m_tree->Branch("m_x_cluster_global", &m_x_cluster_global);
    m_tree->Branch("m_y_cluster_global", &m_y_cluster_global);
    m_tree->Branch("m_z_cluster_global", &m_z_cluster_global);
    m_tree->Branch("m_x_hit_smeared", &m_x_hit_smeared);
    m_tree->Branch("m_y_hit_smeared", &m_y_hit_smeared);
    m_tree->Branch("m_hit_time", &m_hit_time);
    m_tree->Branch("m_hit_time_smeared", &m_hit_time_smeared);
    m_tree->Branch("m_err_x_hit", &m_err_x_hit);
    m_tree->Branch("m_err_y_hit", &m_err_y_hit);
 
    ATH_CHECK(m_hist_svc->regTree("HGTD_ClusterSmearedDigi", m_tree.get()));
  }
 
  return StatusCode::SUCCESS;
}

processAllSubEvents()

StatusCode HGTD_SmearedDigitizationTool::processAllSubEvents ( const EventContext & ctx )

overridevirtual

Reimplemented from PileUpToolBase.

Definition at line 88 of file HGTD_SmearedDigitizationTool.cxx.

                                                                                    {
 
  ATH_MSG_DEBUG("[HGTD_SmearedDigitizationTool::processAllSubEvents]");
 
  auto cluster_container =
      new ClusterContainer_t(m_hgtd_idhelper->wafer_hash_max());
 
  cluster_container->cleanup();
 
  ATH_CHECK(evtStore()->record(cluster_container, m_cluster_name));
 
  auto timed_hit_collection = setupTimedHitCollection();
  HGTD_DetElement_RIO_Map_t det_element_rio_map;
  ATH_CHECK(digitize(ctx,
                     timed_hit_collection,
                     det_element_rio_map));
 
  // TODO: also needs an EventContext (is it the equivalent of createAndStoreRIOs in SiSmearedDigitizationTool)?
  ATH_CHECK(fillClusterContainer (det_element_rio_map,
                                  *cluster_container));
 
  return StatusCode::SUCCESS;
}

processBunchXing()

virtual StatusCode PileUpToolBase::processBunchXing ( int , SubEventIterator , SubEventIterator  )

inlineoverridevirtualinherited

Reimplemented in AFP_PileUpTool, ALFA_PileUpTool, BCM_DigitizationTool, CscDigitizationTool, HGTD_DigitizationTool, LArPileUpTool, LUCID_PileUpTool, MdtDigitizationTool, MergeCalibHitsTool, MergeGenericMuonSimHitCollTool, MergeHijingParsTool, MergeMcEventCollTool, MergeTrackRecordCollTool, MergeTruthJetsTool, MergeTruthParticlesTool, MM_DigitizationTool, MuonR4::MuonDigitizationTool, NewMergeMcEventCollTool, PixelDigitizationTool, PixelFastDigitizationTool, RpcDigitizationTool, SCT_FastDigitizationTool, SimpleMergeMcEventCollTool, SiSmearedDigitizationTool, sTgcDigitizationTool, TestPileUpTool, TgcDigitizationTool, TileHitVecToCntTool, TRTDigitizationTool, TRTFastDigitizationTool, and ZDC_PileUpTool.

Definition at line 43 of file PileUpToolBase.h.

                                                {
    return StatusCode::SUCCESS;
  }

resetFilter()

virtual void PileUpToolBase::resetFilter ( )

inlineoverridevirtualinherited

dummy implementation of filter reset

Reimplemented in MergeTruthJetsTool.

Definition at line 51 of file PileUpToolBase.h.

{ m_filterPassed=true; }

retrieveTruth()

StatusCode HGTD_SmearedDigitizationTool::retrieveTruth ( PRD_MultiTruthCollection *& prd_truth_coll )

private

Definition at line 113 of file HGTD_SmearedDigitizationTool.cxx.

                                                                                     {
 
  if ((evtStore()->contains<PRD_MultiTruthCollection>(m_prd_truth_coll_name))) {
    ATH_CHECK((evtStore()->retrieve(prd_truth_coll, m_prd_truth_coll_name)));
  } else {
    prd_truth_coll = new PRD_MultiTruthCollection;
    ATH_CHECK(evtStore()->record(prd_truth_coll, m_prd_truth_coll_name));
  }
 
  return StatusCode::SUCCESS;
}

setupTimedHitCollection()

TimedHitCollection< SiHit > HGTD_SmearedDigitizationTool::setupTimedHitCollection ( )

private

Definition at line 126 of file HGTD_SmearedDigitizationTool.cxx.

                                                      {
  // Define Hit Collection
  TimedHitCollection<SiHit> timed_hit_coll;
 
  //  get the container(s)
  typedef PileUpMergeSvc::TimedList<SiHitCollection>::type TimedHitCollList_t;
 
  // this is a list<pair<time_t, DataLink<SCTUncompressedHitCollection> > >
  TimedHitCollList_t hit_coll_list;
  unsigned int num_si_hits(0);
  if (!(m_merge_svc
            ->retrieveSubEvtsData(m_si_hit_collection_name.value(), hit_coll_list,
                                  num_si_hits)
            .isSuccess()) and
      hit_coll_list.empty()) {
    ATH_MSG_ERROR("Could not fill TimedHitCollList_t");
    return timed_hit_coll;
  }
 
  ATH_MSG_DEBUG(hit_coll_list.size() << " SiHitCollections with key "
                                     << m_si_hit_collection_name << " found");
 
  timed_hit_coll.reserve(num_si_hits);
 
  // loop on the hit collections
  for(auto hit_coll : hit_coll_list) {
    const SiHitCollection* si_collection(hit_coll.second);
    timed_hit_coll.insert(hit_coll.first, si_collection);
    ATH_MSG_DEBUG("SiHitCollection found with " << si_collection->size()
                                                << " hits");
  }
  return timed_hit_coll;
}

smearMeanTime()

float HGTD_SmearedDigitizationTool::smearMeanTime ( float time, float time_res, CLHEP::HepRandomEngine * rndmEngine )

private

Definition at line 355 of file HGTD_SmearedDigitizationTool.cxx.

                                                                                                               {
  return time + CLHEP::RandGauss::shoot(rndmEngine, 0., time_res);
}

smearPosition()

float HGTD_SmearedDigitizationTool::smearPosition ( float pos, float sig, float boundary, CLHEP::HepRandomEngine * rndmEngine )

private

Definition at line 339 of file HGTD_SmearedDigitizationTool.cxx.

                                                                                                     {
  ATH_MSG_DEBUG("[HGTD_SmearedDigitizationTool::smearPosition] pos: "
                << pos << " sig: " << sig << " boundary: " << boundary);
  float smeared_pos = pos;
  float smear_para = 0.;
  if (sig != 0.) {
    do {
      smear_para = CLHEP::RandGauss::shoot(rndmEngine, 0., sig);
 
    } while (std::abs(smeared_pos + smear_para) > boundary);
    smeared_pos += smear_para;
  }
  return smeared_pos;
}

toProcess()

virtual bool PileUpToolBase::toProcess ( int bunchXing ) const

inlineoverridevirtualinherited

the method this base class helps implementing

Reimplemented in MergeHijingParsTool, and MergeTrackRecordCollTool.

Definition at line 32 of file PileUpToolBase.h.

                                                       {
    //closed interval [m_firstXing,m_lastXing]
    return !((m_firstXing > bunchXing) || (bunchXing > m_lastXing));
  }

Member Data Documentation

m_cluster_name

StringProperty HGTD_SmearedDigitizationTool::m_cluster_name {this, "HGTD_ClustersContainerName", "HGTD_Cluster", "Name of the HGTD cluster container"}

private

Definition at line 93 of file HGTD_SmearedDigitizationTool.h.

{this, "HGTD_ClustersContainerName", "HGTD_Cluster", "Name of the HGTD cluster container"};

m_err_x_hit

float HGTD_SmearedDigitizationTool::m_err_x_hit {}

private

Definition at line 131 of file HGTD_SmearedDigitizationTool.h.

{};

m_err_y_hit

float HGTD_SmearedDigitizationTool::m_err_y_hit {}

private

Definition at line 132 of file HGTD_SmearedDigitizationTool.h.

{};

m_filterPassed

bool PileUpToolBase::m_filterPassed {true}

protectedinherited

Definition at line 60 of file PileUpToolBase.h.

{true}; 

m_firstXing

Gaudi::Property<int> PileUpToolBase::m_firstXing

protectedinherited

Initial value:

{this, "FirstXing", -999,
      "First bunch-crossing in which det is live"}

Definition at line 54 of file PileUpToolBase.h.

                                {this, "FirstXing", -999,
      "First bunch-crossing in which det is live"};

m_hgtd_det_manager

const HGTD_DetectorManager* HGTD_SmearedDigitizationTool::m_hgtd_det_manager {nullptr}

private

Definition at line 85 of file HGTD_SmearedDigitizationTool.h.

{nullptr};

m_hgtd_idhelper

const HGTD_ID* HGTD_SmearedDigitizationTool::m_hgtd_idhelper {nullptr}

private

Definition at line 86 of file HGTD_SmearedDigitizationTool.h.

{nullptr};

m_hist_svc

ServiceHandle<ITHistSvc> HGTD_SmearedDigitizationTool::m_hist_svc {this, "THistSvc", "THistSvc"}

private

Definition at line 108 of file HGTD_SmearedDigitizationTool.h.

{this, "THistSvc", "THistSvc"};

m_hit_time

float HGTD_SmearedDigitizationTool::m_hit_time {}

private

Definition at line 124 of file HGTD_SmearedDigitizationTool.h.

{};

m_hit_time_smeared

float HGTD_SmearedDigitizationTool::m_hit_time_smeared {}

private

Definition at line 125 of file HGTD_SmearedDigitizationTool.h.

{};

m_lastXing

Gaudi::Property<int> PileUpToolBase::m_lastXing

protectedinherited

Initial value:

{this, "LastXing", 999,
      "Last bunch-crossing in which det is live"}

Definition at line 56 of file PileUpToolBase.h.

                                {this, "LastXing", 999,
      "Last bunch-crossing in which det is live"};

m_merge_svc

ServiceHandle<PileUpMergeSvc> HGTD_SmearedDigitizationTool::m_merge_svc {this, "MergeSvc", "PileUpMergeSvc", "Merge service used in HGTD fast digitization"}

private

Definition at line 90 of file HGTD_SmearedDigitizationTool.h.

{this, "MergeSvc", "PileUpMergeSvc", "Merge service used in HGTD fast digitization"}; 

m_output_file

std::unique_ptr<TFile> HGTD_SmearedDigitizationTool::m_output_file

private

Definition at line 109 of file HGTD_SmearedDigitizationTool.h.

m_pitch_x

FloatProperty HGTD_SmearedDigitizationTool::m_pitch_x {this, "pitch_X", 1.3, ""}

private

Definition at line 102 of file HGTD_SmearedDigitizationTool.h.

{this, "pitch_X", 1.3, ""};

m_pitch_y

FloatProperty HGTD_SmearedDigitizationTool::m_pitch_y {this, "pitch_Y", 1.3, ""}

private

Definition at line 103 of file HGTD_SmearedDigitizationTool.h.

{this, "pitch_Y", 1.3, ""};

m_prd_truth_coll_name

std::string HGTD_SmearedDigitizationTool::m_prd_truth_coll_name {"PRD_MultiTruthPixel"}

private

Definition at line 88 of file HGTD_SmearedDigitizationTool.h.

{"PRD_MultiTruthPixel"};

m_rndm_svc

ServiceHandle<IAthRNGSvc> HGTD_SmearedDigitizationTool::m_rndm_svc {this, "RndmSvc", "AthRNGSvc", "Random number service used in HGTD fast digitization"}

private

Random number service.

Definition at line 91 of file HGTD_SmearedDigitizationTool.h.

{this, "RndmSvc", "AthRNGSvc", "Random number service used in HGTD fast digitization"}; 

m_si_hit_collection_name

StringProperty HGTD_SmearedDigitizationTool::m_si_hit_collection_name {this, "SiHitCollectionName", "HGTD_Hits", "Name of the Si Hit collection"}

private

Definition at line 95 of file HGTD_SmearedDigitizationTool.h.

{this, "SiHitCollectionName", "HGTD_Hits", "Name of the Si Hit collection"};

m_smear_intersection_position

BooleanProperty HGTD_SmearedDigitizationTool::m_smear_intersection_position {this, "SmearIntersectionPosition", true, ""}

private

Definition at line 97 of file HGTD_SmearedDigitizationTool.h.

{this, "SmearIntersectionPosition", true, ""};

m_smear_mean_time

BooleanProperty HGTD_SmearedDigitizationTool::m_smear_mean_time {this, "SmearMeanTime", true, ""}

private

Definition at line 98 of file HGTD_SmearedDigitizationTool.h.

{this, "SmearMeanTime", true, ""};

m_time_res

float HGTD_SmearedDigitizationTool::m_time_res {0.035}

private

Definition at line 104 of file HGTD_SmearedDigitizationTool.h.

{0.035}; // time resolution in ns

m_tree

std::unique_ptr<TTree> HGTD_SmearedDigitizationTool::m_tree

private

Definition at line 110 of file HGTD_SmearedDigitizationTool.h.

m_vetoPileUpTruthLinks

Gaudi::Property<int> PileUpToolBase::m_vetoPileUpTruthLinks

protectedinherited

Initial value:

{this, "VetoPileUpTruthLinks", true,
      "Ignore links to suppressed pile-up truth"}

Definition at line 58 of file PileUpToolBase.h.

                                           {this, "VetoPileUpTruthLinks", true,
      "Ignore links to suppressed pile-up truth"};

m_write_tree

BooleanProperty HGTD_SmearedDigitizationTool::m_write_tree {this, "WriteTree", false, ""}

private

Definition at line 99 of file HGTD_SmearedDigitizationTool.h.

{this, "WriteTree", false, ""};

m_x_cluster_global

float HGTD_SmearedDigitizationTool::m_x_cluster_global {}

private

Definition at line 127 of file HGTD_SmearedDigitizationTool.h.

{};

m_x_entry_hit

float HGTD_SmearedDigitizationTool::m_x_entry_hit {}

private

Definition at line 117 of file HGTD_SmearedDigitizationTool.h.

{};

m_x_exit_hit

float HGTD_SmearedDigitizationTool::m_x_exit_hit {}

private

Definition at line 120 of file HGTD_SmearedDigitizationTool.h.

{};

m_x_hit

float HGTD_SmearedDigitizationTool::m_x_hit {}

private

Definition at line 112 of file HGTD_SmearedDigitizationTool.h.

{};

m_x_hit_smeared

float HGTD_SmearedDigitizationTool::m_x_hit_smeared {}

private

Definition at line 114 of file HGTD_SmearedDigitizationTool.h.

{};

m_y_cluster_global

float HGTD_SmearedDigitizationTool::m_y_cluster_global {}

private

Definition at line 128 of file HGTD_SmearedDigitizationTool.h.

{};

m_y_entry_hit

float HGTD_SmearedDigitizationTool::m_y_entry_hit {}

private

Definition at line 118 of file HGTD_SmearedDigitizationTool.h.

{};

m_y_exit_hit

float HGTD_SmearedDigitizationTool::m_y_exit_hit {}

private

Definition at line 121 of file HGTD_SmearedDigitizationTool.h.

{};

m_y_hit

float HGTD_SmearedDigitizationTool::m_y_hit {}

private

Definition at line 113 of file HGTD_SmearedDigitizationTool.h.

{};

m_y_hit_smeared

float HGTD_SmearedDigitizationTool::m_y_hit_smeared {}

private

Definition at line 115 of file HGTD_SmearedDigitizationTool.h.

{};

m_z_cluster_global

float HGTD_SmearedDigitizationTool::m_z_cluster_global {}

private

Definition at line 129 of file HGTD_SmearedDigitizationTool.h.

{};

m_z_entry_hit

float HGTD_SmearedDigitizationTool::m_z_entry_hit {}

private

Definition at line 119 of file HGTD_SmearedDigitizationTool.h.

{};

m_z_exit_hit

float HGTD_SmearedDigitizationTool::m_z_exit_hit {}

private

Definition at line 122 of file HGTD_SmearedDigitizationTool.h.

{};

---

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

• HGTD_SmearedDigitizationTool.h
• HGTD_SmearedDigitizationTool.cxx