AuxDetSimChannel.cxx
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 ///
3 /// \file lardataobj/Simulation/AuxDetSimChannel.cxx
4 ///
5 /// \author miceli@fnal.gov
6 ///
7 ////////////////////////////////////////////////////////////////////////
8 
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>=0? ide.trackID+offset : 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  auto tid = std::abs(ide.trackID)+offset;
104 
105  if( tid < range_trackID.first )
106  range_trackID.first = tid;
107  if( tid > range_trackID.second )
108  range_trackID.second = tid;
109  }
110 
111  return range_trackID;
112  }
113 
114  //----------------------------------------------------------------------------
115  bool AuxDetSimChannel::operator< (const sim::AuxDetSimChannel& other) const
116  {
117  if(fAuxDetID < other.AuxDetID() ) return true;
118 
119  return fAuxDetSensitiveID < other.AuxDetSensitiveID();
120  }
121 
122  //----------------------------------------------------------------------------
123  bool AuxDetSimChannel::operator== (const sim::AuxDetSimChannel& other) const
124  {
125  return (fAuxDetID == other.AuxDetID() && fAuxDetSensitiveID == other.AuxDetSensitiveID());
126  }
127 
128 
129 }//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: SignalShapingService.h:21
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: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:115
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:42
sim::AuxDetSimChannel
Collection of particles crossing one auxiliary detector cell.
Definition: AuxDetSimChannel.h:60
abs
T abs(T value)
Definition: sparse_vector_test.cc:30
sim::AuxDetSimChannel::AuxDetID
uint32_t AuxDetID() const
Definition: AuxDetSimChannel.h:110
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:112
sim::AuxDetIDE::exitZ
float exitZ
Exit position Z of particle.
Definition: AuxDetSimChannel.h:43
max
static int max(int a, int b)
Definition: fortranscanner.cpp:60366
sim
Code to link reconstructed objects back to the MC truth information.
Definition: PedestalAdditionService.h:17
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
cet::sqlite::min
T min(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:55
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:111
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:123
Generated by   doxygen 1.8.11