ATLAS Offline Software
Functions
HepMC::StatusBased Namespace Reference

Functions

template<class T >
bool no_truth_link (const T &p)
 Method to establish if a particle (or barcode) corresponds to truth-suppressed pile-up TODO implement a status/id based version of is_truth_suppressed_pileup. More...
 
template<class T >
bool is_simulation_particle (const T &p)
 Helper function for SDO creation in PileUpTools TODO implement a status/id based version of ignoreTruthLink. More...
 
template<class T >
bool is_sim_secondary (const T &p)
 Method to establish if a particle is a new seondary created during the simulation based on the status value. More...
 
template<class T >
int generations (const T &p)
 Method to return how many interactions a particle has undergone during simulation based on the status value. More...
 
template<class T >
bool is_simulation_vertex (const T &v)
 Method to establish if the vertex was created during simulation from the status. More...
 
template<class T1 , class T2 >
bool is_same_generator_particle (const T1 &p1, const T2 &p2)
 Method to establish if two particles in the GenEvent actually represent the same generated particle. More...
 
template<class T1 , class T2 >
bool is_sim_descendant (const T1 &p1, const T2 &p2)
 Method to check if the first particle is a descendant of the second in the simulation, i.e. particle p1 was produced simulations particle p2. More...
 

Function Documentation

◆ generations()

template<class T >
int HepMC::StatusBased::generations ( const T &  p)
inline

Method to return how many interactions a particle has undergone during simulation based on the status value.

Definition at line 235 of file MagicNumbers.h.

235 { return (status(p)/SIM_STATUS_INCREMENT);}

◆ is_same_generator_particle()

template<class T1 , class T2 >
bool HepMC::StatusBased::is_same_generator_particle ( const T1 &  p1,
const T2 &  p2 
)
inline

Method to establish if two particles in the GenEvent actually represent the same generated particle.

Definition at line 242 of file MagicNumbers.h.

242  {
243  const int id1 = uniqueID(p1);
244  const int id2 = uniqueID(p2);
245  if (id1 == id2) { return true;} // simplest case
246  const int generations1 = generations(p1);
247  const int generations2 = generations(p2);
248  if (generations1 == generations2) { return false; } // if the id values don't match and the particles have the same generation number then they cannot be the same particle.
249  if constexpr (std::is_same_v<std::remove_const_t<remove_smart_pointer_t<std::remove_pointer_t<T1>>>, TrackRecord>) {
250  // No choice, but to get the history of one of the particles:
251  const int direction = (generations2 > generations1) ? -1 : 1;
252  std::deque<int> history2 = simulation_history( p2, direction );
253  if (std::find(history2.begin(),history2.end(), id1) == history2.end()) { return false; }
254  }
255  else {
256  // No choice, but to get the history of one of the particles:
257  const int direction = (generations1 > generations2) ? -1 : 1;
258  std::deque<int> history1 = simulation_history( p1, direction );
259  if (std::find(history1.begin(),history1.end(), id2) == history1.end()) { return false; }
260  }
261  return true;
262  }

◆ is_sim_descendant()

template<class T1 , class T2 >
bool HepMC::StatusBased::is_sim_descendant ( const T1 &  p1,
const T2 &  p2 
)
inline

Method to check if the first particle is a descendant of the second in the simulation, i.e. particle p1 was produced simulations particle p2.

Definition at line 266 of file MagicNumbers.h.

266  {
267  const int id1 = uniqueID(p1);
268  const int id2 = uniqueID(p2);
269  if (id1 == id2) { return true;} // simplest case
270  const int generations1 = generations(p1);
271  const int generations2 = generations(p2);
272  if (generations1 == generations2) { return false; } // if the id values don't match and the particles have the same generation number then they cannot be the same particle.
273  if constexpr (std::is_same_v<std::remove_const_t<remove_smart_pointer_t<std::remove_pointer_t<T1>>>, TrackRecord>) {
274  // No choice, but to get the descendents of p2
275  constexpr int descendents = 1;
276  std::deque<int> history2 = simulation_history( p2, descendents );
277  if (std::find(history2.begin(),history2.end(), id1) == history2.end()) { return false; }
278  }
279  else {
280  // No choice, but to get the history of p1
281  constexpr int ancestors = -1;
282  std::deque<int> history1 = simulation_history( p1, ancestors );
283  if (std::find(history1.begin(),history1.end(), id2) == history1.end()) { return false; }
284  }
285  return true;
286  }

◆ is_sim_secondary()

template<class T >
bool HepMC::StatusBased::is_sim_secondary ( const T &  p)
inline

Method to establish if a particle is a new seondary created during the simulation based on the status value.

Definition at line 232 of file MagicNumbers.h.

◆ is_simulation_particle()

template<class T >
bool HepMC::StatusBased::is_simulation_particle ( const T &  p)
inline

Helper function for SDO creation in PileUpTools TODO implement a status/id based version of ignoreTruthLink.

Method to establish if a particle was created during the simulation based on the status value

Definition at line 229 of file MagicNumbers.h.

229 { return (status(p)>SIM_STATUS_THRESHOLD);}

◆ is_simulation_vertex()

template<class T >
bool HepMC::StatusBased::is_simulation_vertex ( const T &  v)
inline

Method to establish if the vertex was created during simulation from the status.

Definition at line 238 of file MagicNumbers.h.

238 { return (status(v)>SIM_STATUS_THRESHOLD);}

◆ no_truth_link()

template<class T >
bool HepMC::StatusBased::no_truth_link ( const T &  p)
inline

Method to establish if a particle (or barcode) corresponds to truth-suppressed pile-up TODO implement a status/id based version of is_truth_suppressed_pileup.

Method to establish if a if the object is linked to something which was never saved to the HepMC Truth - for example particle was too low energy to be recorded

Definition at line 223 of file MagicNumbers.h.

223 { return (uniqueID(p) == 0);}
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
HepMC::SIM_STATUS_INCREMENT
constexpr int SIM_STATUS_INCREMENT
Constant defining the barcode threshold for regenerated particles, i.e. particles surviving an intera...
Definition: MagicNumbers.h:35
id2
HWIdentifier id2
Definition: LArRodBlockPhysicsV0.cxx:564
HepMC::simulation_history
std::deque< int > simulation_history(const T &p, int direction)
Function to calculate all the descendants(direction=1)/ancestors(direction=-1) of the particle.
Definition: MagicNumbers.h:185
HepMC::uniqueID
int uniqueID(const T &p)
Definition: MagicNumbers.h:113
HepMC::SIM_STATUS_THRESHOLD
constexpr int SIM_STATUS_THRESHOLD
Constant definiting the status threshold for simulated particles, eg. can be used to separate generat...
Definition: MagicNumbers.h:38
HepMC::status
int status(const GenVertexPtr &v1)
Definition: MagicNumbers.h:152
TrackRecord
Definition: TrackRecord.h:10
python.PyAthena.v
v
Definition: PyAthena.py:157
HepMC::remove_smart_pointer_t
typename remove_smart_pointer< T >::type remove_smart_pointer_t
Definition: MagicNumbers.h:74
HepMC::generations
int generations(const T &p)
Method to return how many interactions a particle has undergone during simulation (TODO migrate to be...
Definition: MagicNumbers.h:302