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