LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
AuxDetSimChannel.cxx
Go to the documentation of this file.
1 
9 // our header
10 #include "lardataobj/Simulation/AuxDetSimChannel.h"
11 
12 // C/C++ standard library
13 #include <limits> // std::numeric_limits<>
14 #include <stdexcept>
15 
16 // LArSoft headers
17 #include "larcoreobj/SimpleTypesAndConstants/PhysicalConstants.h" // util::kBogusX
18 
19 namespace sim {
20 
21  // Default constructor
22  //-------------------------------------------------
23  AuxDetIDE::AuxDetIDE()
24  : trackID(util::kBogusI)
25  , energyDeposited(util::kBogusF)
26  , entryX(util::kBogusF)
27  , entryY(util::kBogusF)
28  , entryZ(util::kBogusF)
29  , entryT(util::kBogusF)
30  , exitX(util::kBogusF)
31  , exitY(util::kBogusF)
32  , exitZ(util::kBogusF)
33  , exitT(util::kBogusF)
34  , exitMomentumX(util::kBogusF)
35  , exitMomentumY(util::kBogusF)
36  , exitMomentumZ(util::kBogusF)
37  {}
38 
39  // Copy with offset constructor
40  //-------------------------------------------------
41  AuxDetIDE::AuxDetIDE(AuxDetIDE const& ide, int offset)
42  : trackID(ide.trackID >= 0 ? ide.trackID + offset : ide.trackID - offset)
43  , energyDeposited(ide.energyDeposited)
44  , entryX(ide.entryX)
45  , entryY(ide.entryY)
46  , entryZ(ide.entryZ)
47  , entryT(ide.entryT)
48  , exitX(ide.exitX)
49  , exitY(ide.exitY)
50  , exitZ(ide.exitZ)
51  , exitT(ide.exitT)
52  , exitMomentumX(ide.exitMomentumX)
53  , exitMomentumY(ide.exitMomentumY)
54  , exitMomentumZ(ide.exitMomentumZ)
55  {}
56 
57  //----------------------------------------------------------------------------
58  AuxDetSimChannel::AuxDetSimChannel()
59  : fAuxDetID(std::numeric_limits<uint32_t>::max())
60  , fAuxDetSensitiveID(std::numeric_limits<uint32_t>::max())
61  {}
62 
63  //----------------------------------------------------------------------------
64  AuxDetSimChannel::AuxDetSimChannel(uint32_t inputAuxDetID, uint32_t inputAuxDetSensitiveID)
65  : fAuxDetID(inputAuxDetID), fAuxDetSensitiveID(inputAuxDetSensitiveID)
66  {}
67 
68  //----------------------------------------------------------------------------
69  AuxDetSimChannel::AuxDetSimChannel(uint32_t inputAuxDetID,
70  const std::vector<sim::AuxDetIDE>& inputAuxDetIDEs,
71  uint32_t inputAuxDetSensitiveID)
72  : fAuxDetID(inputAuxDetID)
73  , fAuxDetSensitiveID(inputAuxDetSensitiveID)
74  , fAuxDetIDEs(inputAuxDetIDEs)
75  {}
76 
77  //----------------------------------------------------------------------------
78  AuxDetSimChannel::AuxDetSimChannel(uint32_t inputAuxDetID,
79  std::vector<sim::AuxDetIDE>&& inputAuxDetIDEs,
80  uint32_t inputAuxDetSensitiveID)
81  : fAuxDetID(inputAuxDetID)
82  , fAuxDetSensitiveID(inputAuxDetSensitiveID)
83  , fAuxDetIDEs(inputAuxDetIDEs)
84  {}
85 
86  //----------------------------------------------------------------------------
87  std::pair<int, int> AuxDetSimChannel::MergeAuxDetSimChannel(const AuxDetSimChannel& chan,
88  int offset)
89  {
90  if (this->fAuxDetID != chan.AuxDetID() && this->fAuxDetSensitiveID != chan.AuxDetSensitiveID())
91  throw std::runtime_error("ERROR AuxDetSimChannel Merge: Trying to merge different channels!");
92 
93  std::pair<int, int> range_trackID(std::numeric_limits<int>::max(),
94  std::numeric_limits<int>::min());
95 
96  for (auto const& ide : AuxDetIDEs()) {
97  this->fAuxDetIDEs.emplace_back(ide, offset);
98 
99  auto tid = std::abs(ide.trackID) + offset;
100 
101  if (tid < range_trackID.first) range_trackID.first = tid;
102  if (tid > range_trackID.second) range_trackID.second = tid;
103  }
104 
105  return range_trackID;
106  }
107 
108  //----------------------------------------------------------------------------
109  bool AuxDetSimChannel::operator<(const sim::AuxDetSimChannel& other) const
110  {
111  if (fAuxDetID < other.AuxDetID()) return true;
112 
113  return fAuxDetSensitiveID < other.AuxDetSensitiveID();
114  }
115 
116  //----------------------------------------------------------------------------
117  bool AuxDetSimChannel::operator==(const sim::AuxDetSimChannel& other) const
118  {
119  return (fAuxDetID == other.AuxDetID() && fAuxDetSensitiveID == other.AuxDetSensitiveID());
120  }
121 
122 } //namespace sim
sim::AuxDetSimChannel::fAuxDetSensitiveID
uint32_t fAuxDetSensitiveID
integer used to retrieve AuxDetSensitiveGeo object
Definition: AuxDetSimChannel.h:68
sim::AuxDetSimChannel::fAuxDetIDEs
std::vector< sim::AuxDetIDE > fAuxDetIDEs
one sim::AuxDetIDE for each G4 track id
Definition: AuxDetSimChannel.h:69
util
Namespace for general, non-LArSoft-specific utilities.
Definition: PIDAAlg.h:26
sim::AuxDetSimChannel::MergeAuxDetSimChannel
std::pair< int, int > MergeAuxDetSimChannel(const AuxDetSimChannel &, int)
Definition: AuxDetSimChannel.cxx:87
sim::AuxDetIDE::trackID
int trackID
Geant4 supplied track ID.
Definition: AuxDetSimChannel.h:35
sim::AuxDetIDE::exitMomentumZ
float exitMomentumZ
Exit Z-Momentum of particle.
Definition: AuxDetSimChannel.h:47
sim::AuxDetSimChannel::operator<
bool operator<(const AuxDetSimChannel &other) const
Definition: AuxDetSimChannel.cxx:109
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
sim::AuxDetSimChannel::AuxDetSimChannel
AuxDetSimChannel()
Default constructor (invalid, empty data)
Definition: AuxDetSimChannel.cxx:58
std
STL namespace.
util::kBogusI
constexpr int kBogusI
obviously bogus integer value
Definition: PhysicalConstants.h:72
sim::AuxDetIDE::exitY
float exitY
Exit position Y of particle.
Definition: AuxDetSimChannel.h:42
sim::AuxDetSimChannel
Collection of particles crossing one auxiliary detector cell.
Definition: AuxDetSimChannel.h:60
AuxDetSimChannel.h
object containing MC truth information necessary for making RawDigits and doing back tracking ...
sim::AuxDetSimChannel::AuxDetID
uint32_t AuxDetID() const
Definition: AuxDetSimChannel.h:111
sim::AuxDetIDE::entryT
float entryT
Entry time of particle.
Definition: AuxDetSimChannel.h:40
sim::AuxDetSimChannel::fAuxDetID
uint32_t fAuxDetID
geo->AuxDet(auxDetID), integer used to retrieve AuxDetGeo objec
Definition: AuxDetSimChannel.h:67
sim::AuxDetIDE::exitMomentumX
float exitMomentumX
Exit X-Momentum of particle.
Definition: AuxDetSimChannel.h:45
sim::AuxDetIDE::exitT
float exitT
Exit time of particle.
Definition: AuxDetSimChannel.h:44
sim::AuxDetSimChannel::AuxDetIDEs
std::vector< sim::AuxDetIDE > const & AuxDetIDEs() const
Definition: AuxDetSimChannel.h:119
sim::AuxDetIDE::exitZ
float exitZ
Exit position Z of particle.
Definition: AuxDetSimChannel.h:43
sim
Monte Carlo Simulation.
Definition: FlashHypothesisCreator.h:21
sim::AuxDetIDE::entryZ
float entryZ
Entry position Z of particle.
Definition: AuxDetSimChannel.h:39
sim::AuxDetIDE::exitX
float exitX
Exit position X of particle.
Definition: AuxDetSimChannel.h:41
sim::AuxDetIDE::energyDeposited
float energyDeposited
total energy deposited for this track ID and time
Definition: AuxDetSimChannel.h:36
util::kBogusF
constexpr float kBogusF
obviously bogus float value
Definition: PhysicalConstants.h:73
sim::AuxDetIDE::entryX
float entryX
Entry position X of particle.
Definition: AuxDetSimChannel.h:37
sim::AuxDetIDE::entryY
float entryY
Entry position Y of particle.
Definition: AuxDetSimChannel.h:38
sim::AuxDetSimChannel::AuxDetSensitiveID
uint32_t AuxDetSensitiveID() const
Definition: AuxDetSimChannel.h:115
sim::AuxDetIDE
MC truth information to make RawDigits and do back tracking.
Definition: AuxDetSimChannel.h:27
sim::AuxDetIDE::exitMomentumY
float exitMomentumY
Exit Y-Momentum of particle.
Definition: AuxDetSimChannel.h:46
PhysicalConstants.h
Collection of Physical constants used in LArSoft.
sim::AuxDetSimChannel::operator==
bool operator==(const AuxDetSimChannel &other) const
Definition: AuxDetSimChannel.cxx:117