Public Member Functions | Private Member Functions | Private Attributes | List of all members
opdet::DUNE35tonOpDetResponse Class Reference

#include <DUNE35tonOpDetResponse.h>

Inheritance diagram for opdet::DUNE35tonOpDetResponse:
opdet::OpDetResponseInterface

Public Member Functions

 DUNE35tonOpDetResponse (fhicl::ParameterSet const &pset, art::ActivityRegistry &reg)
 
 ~DUNE35tonOpDetResponse () throw ()
 
- Public Member Functions inherited from opdet::OpDetResponseInterface
virtual ~OpDetResponseInterface ()=default
 
virtual void reconfigure (fhicl::ParameterSet const &p)
 
virtual int NOpChannels () const
 
virtual int readoutToGeoChannel (int readoutChannel) const
 
virtual bool detected (int OpChannel, const sim::OnePhoton &Phot, int &newOpChannel) const
 
virtual bool detected (int OpChannel, const sim::OnePhoton &Phot) const
 
virtual bool detectedLite (int OpChannel, int &newOpChannel) const
 
virtual bool detectedLite (int OpChannel) const
 
virtual float wavelength (double energy) const
 

Private Member Functions

virtual void doReconfigure (fhicl::ParameterSet const &p)
 
virtual int doNOpChannels () const
 
virtual bool doDetected (int OpDet, const sim::OnePhoton &Phot, int &newOpChannel) const
 
virtual bool doDetectedLite (int OpDet, int &newOpChannel) const
 

Private Attributes

float fQE
 
float fWavelengthCutLow
 
float fWavelengthCutHigh
 
bool fLightGuideAttenuation
 
std::string fChannelConversion
 
bool fFullSimChannelConvert
 
bool fFastSimChannelConvert
 
int fLongAxis
 

Detailed Description

Definition at line 21 of file DUNE35tonOpDetResponse.h.

Constructor & Destructor Documentation

opdet::DUNE35tonOpDetResponse::DUNE35tonOpDetResponse ( fhicl::ParameterSet const &  pset,
art::ActivityRegistry &  reg 
)

Definition at line 24 of file DUNE35tonOpDetResponse_service.cc.

26  {
27  this->doReconfigure(pset);
28  }
opdet::DUNE35tonOpDetResponse::doReconfigure
virtual void doReconfigure(fhicl::ParameterSet const &p)
Definition: DUNE35tonOpDetResponse_service.cc:36
opdet::DUNE35tonOpDetResponse::~DUNE35tonOpDetResponse ( )
throw (
)

Definition at line 31 of file DUNE35tonOpDetResponse_service.cc.

32  { }

Member Function Documentation

bool opdet::DUNE35tonOpDetResponse::doDetected ( int  OpDet,
const sim::OnePhoton &  Phot,
int &  newOpChannel 
) const
privatevirtual

Implements opdet::OpDetResponseInterface.

Definition at line 91 of file DUNE35tonOpDetResponse_service.cc.

92  {
93 
94  // Find the Optical Detector using the geometry service
95  art::ServiceHandle<geo::Geometry> geom;
96  const TGeoNode* node = geom->OpDetGeoFromOpDet(OpDet).Node();
97 
98  // Identify the photon detector type
99  int pdtype;
100  std::string detname = node->GetName();
101  boost::to_lower(detname);
102  if (detname.find("bar") != std::string::npos ) pdtype = 0;
103  else if (detname.find("fiber") != std::string::npos ) pdtype = 1;
104  else if (detname.find("plank") != std::string::npos ) pdtype = 2;
105  else pdtype = -1;
106 
107 
108  if (fFullSimChannelConvert){
109  // Override default number of channels for Fiber and Plank
110  float NOpHardwareChannels = geom->NOpHardwareChannels(OpDet);
111  if (pdtype == 1) NOpHardwareChannels = 3;
112  if (pdtype == 2) NOpHardwareChannels = 2;
113 
114  int hardwareChannel = (int) ( CLHEP::RandFlat::shoot(1.0) * NOpHardwareChannels );
115  newOpChannel = geom->OpChannel(OpDet, hardwareChannel);
116  }
117  else{
118  newOpChannel = OpDet;
119  }
120 
121  // Check QE
122  if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
123 
124  double wavel = wavelength(Phot.Energy);
125  // Check wavelength acceptance
126  if (wavel < fWavelengthCutLow) return false;
127  if (wavel > fWavelengthCutHigh) return false;
128 
129  if (fLightGuideAttenuation) {
130  // Get the length of the photon detector
131  TGeoBBox *box = (TGeoBBox*)node->GetVolume()->GetShape();
132  double opdetHalfLength = 0;
133  double sipmDistance = 0;
134 
135  if (fLongAxis == 0) {
136  opdetHalfLength = box->GetDX();
137  sipmDistance = opdetHalfLength - Phot.FinalLocalPosition.x();
138  }
139  else if (fLongAxis == 1) {
140  opdetHalfLength = box->GetDY();
141  sipmDistance = opdetHalfLength - Phot.FinalLocalPosition.y();
142  }
143  else if (fLongAxis == 2) {
144  opdetHalfLength = box->GetDZ();
145  sipmDistance = opdetHalfLength - Phot.FinalLocalPosition.z();
146  }
147  else {
148  mf::LogError("DUNE35tonOpDetResponse") << "Unknown axis, fLongAxis = " << fLongAxis;
149  assert(false);
150  }
151 
152 
153 
154  if (pdtype == 0) {
155  double normalize = 0.6961; // Normalize mean performance to be the same in all PDs
156  double lambdaShort = 5.56; // cm
157  double lambdaLong = 44.13; // cm
158  double fracShort = 0.40 * normalize;
159  double fracLong = 0.60 * normalize;
160 
161  double AttenuationProb = fracShort*exp(-sipmDistance/lambdaShort) + fracLong*exp(-sipmDistance/lambdaLong);
162 
163  //std::cout <<"GRAPH: " << OpDet << "," << Phot.FinalLocalPosition.y() << "," << sipmDistance << "," << AttenuationProb << std::endl;
164 
165  //mf::LogVerbatim("DUNE35tonOpDetResponse") << "OpDet: " << OpDet << " is a " << pdtype
166  // << " with length " << opdetHalfLength << " in detector "
167  // << box->GetDX() << " x " << box->GetDY() << " x " << box->GetDZ()
168  // << " named " << detname;
169  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Local Position = (" << Phot.FinalLocalPosition.x()
170  // << ", " << Phot.FinalLocalPosition.y() << ", " << Phot.FinalLocalPosition.z() << ")";
171  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Distance to SiPM = " << sipmDistance << " along axis " << fLongAxis;
172  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Attenuation Probability = " << AttenuationProb;
173 
174  // Throw away some photons based on attenuation
175  if ( CLHEP::RandFlat::shoot(1.0) > AttenuationProb ) return false;
176  }
177  else if (pdtype == 1) {
178  double normalize = 1.0; // Normalize mean performance to be the same in all PDs
179  double lambda = 14.6; // cm
180 
181  // Throw away some photons based on attenuation
182  double AttenuationProb = normalize*exp(-sipmDistance/lambda);
183 
184  //std::cout <<"GRAPH: " << OpDet << "," << Phot.FinalLocalPosition.y() << "," << sipmDistance << "," << AttenuationProb << std::endl;
185 
186  //mf::LogVerbatim("DUNE35tonOpDetResponse") << "OpDet: " << OpDet << " is a " << pdtype
187  // << " with length " << opdetHalfLength << " in detector "
188  // << box->GetDX() << " x " << box->GetDY() << " x " << box->GetDZ()
189  // << " named " << detname;
190  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Local Position = (" << Phot.FinalLocalPosition.x()
191  // << ", " << Phot.FinalLocalPosition.y() << ", " << Phot.FinalLocalPosition.z() << ")";
192  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Distance to SiPM = " << sipmDistance << " along axis " << fLongAxis;
193  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Attenuation Probability = " << AttenuationProb;
194 
195  // Throw away some photons based on attenuation
196  if ( CLHEP::RandFlat::shoot(1.0) > AttenuationProb ) return false;
197  }
198  else if (pdtype == 2) {
199  double normalize = 0.4498; // Normalize mean performance to be the same in all PDs
200  double lambda = 48.4; // cm
201  double altDistance = 2*opdetHalfLength - sipmDistance;
202  double frac = 0.5 * normalize;
203 
204  double AttenuationProb = frac*exp(-sipmDistance/lambda) + frac*exp(-altDistance/lambda);
205 
206  //std::cout <<"GRAPH: " << OpDet << "," << Phot.FinalLocalPosition.y() << "," << sipmDistance << "," << AttenuationProb << std::endl;
207 
208  //mf::LogVerbatim("DUNE35tonOpDetResponse") << "OpDet: " << OpDet << " is a " << pdtype
209  // << " with length " << opdetHalfLength << " in detector "
210  // << box->GetDX() << " x " << box->GetDY() << " x " << box->GetDZ()
211  // << " named " << detname;
212  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Local Position = (" << Phot.FinalLocalPosition.x()
213  // << ", " << Phot.FinalLocalPosition.y() << ", " << Phot.FinalLocalPosition.z() << ")";
214  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Distance to SiPM = " << sipmDistance << " along axis " << fLongAxis;
215  //mf::LogVerbatim("DUNE35tonOpDetResponse") << " Attenuation Probability = " << AttenuationProb;
216 
217  // Throw away some photons based on attenuation
218  if ( CLHEP::RandFlat::shoot(1.0) > AttenuationProb ) return false;
219  }
220  else {
221  mf::LogWarning("DUNE35tonOpDetResponse") << "OpDet: " << OpDet << " is an unknown PD type named: " << detname
222  << ". Assuming no attenuation.";
223  }
224 
225  }
226 
227  return true;
228  }
geo::OpDetGeo::Node
const TGeoNode * Node() const
Returns the ROOT object describing the detector geometry.
Definition: OpDetGeo.h:146
generate_datataset.normalize
def normalize(value)
Definition: generate_datataset.py:81
geo::GeometryCore::OpDetGeoFromOpDet
OpDetGeo const & OpDetGeoFromOpDet(unsigned int OpDet) const
Returns the geo::OpDetGeo object for the given detector number.
Definition: GeometryCore.cxx:1830
string
std::string string
Definition: nybbler.cc:12
geo::GeometryCore::NOpHardwareChannels
unsigned int NOpHardwareChannels(int opDet) const
Definition: GeometryCore.cxx:265
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:211
keras_to_tensorflow.int
int
Definition: keras_to_tensorflow.py:69
geo::GeometryCore::OpChannel
unsigned int OpChannel(int detNum, int hardwareChannel) const
Convert detector number and hardware channel to unique channel.
Definition: GeometryCore.cxx:271
NOpHardwareChannels
Index NOpHardwareChannels(Index opDet)
Definition: sspmapmaker_v1.c:195
sim::OnePhoton::FinalLocalPosition
geo::OpticalPoint_t FinalLocalPosition
Where photon enters the optical detector in local coordinates [cm].
Definition: SimPhotons.h:75
mf::LogWarning
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
Definition: MessageLogger.h:210
sim::OnePhoton::Energy
float Energy
Scintillation photon energy [GeV].
Definition: SimPhotons.h:82
opdet::OpDetResponseInterface::wavelength
virtual float wavelength(double energy) const
Definition: OpDetResponseInterface.h:126
bool opdet::DUNE35tonOpDetResponse::doDetectedLite ( int  OpDet,
int &  newOpChannel 
) const
privatevirtual

Implements opdet::OpDetResponseInterface.

Definition at line 231 of file DUNE35tonOpDetResponse_service.cc.

232  {
233  if (fFastSimChannelConvert){
234 
235  // Find the Optical Detector using the geometry service
236  art::ServiceHandle<geo::Geometry> geom;
237  // Here OpDet must be opdet since we are introducing
238  // channel mapping here.
239  const TGeoNode* node = geom->OpDetGeoFromOpDet(OpDet).Node();
240 
241 
242  // Identify the photon detector type
243  int pdtype;
244  std::string detname = node->GetName();
245  boost::to_lower(detname);
246  if (detname.find("bar") != std::string::npos ) pdtype = 0;
247  else if (detname.find("fiber") != std::string::npos ) pdtype = 1;
248  else if (detname.find("plank") != std::string::npos ) pdtype = 2;
249  else pdtype = -1;
250 
251  // Override default number of channels for Fiber and Plank
252  float NOpHardwareChannels = geom->NOpHardwareChannels(OpDet);
253  if (pdtype == 1) NOpHardwareChannels = 8;
254  if (pdtype == 2) NOpHardwareChannels = 2;
255 
256  int hardwareChannel = (int) ( CLHEP::RandFlat::shoot(1.0) * NOpHardwareChannels );
257  newOpChannel = geom->OpChannel(OpDet, hardwareChannel);
258  }
259  else{
260  newOpChannel = OpDet;
261  }
262 
263  // Check QE
264  if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
265 
266 
267  return true;
268  }
geo::OpDetGeo::Node
const TGeoNode * Node() const
Returns the ROOT object describing the detector geometry.
Definition: OpDetGeo.h:146
geo::GeometryCore::OpDetGeoFromOpDet
OpDetGeo const & OpDetGeoFromOpDet(unsigned int OpDet) const
Returns the geo::OpDetGeo object for the given detector number.
Definition: GeometryCore.cxx:1830
string
std::string string
Definition: nybbler.cc:12
geo::GeometryCore::NOpHardwareChannels
unsigned int NOpHardwareChannels(int opDet) const
Definition: GeometryCore.cxx:265
keras_to_tensorflow.int
int
Definition: keras_to_tensorflow.py:69
geo::GeometryCore::OpChannel
unsigned int OpChannel(int detNum, int hardwareChannel) const
Convert detector number and hardware channel to unique channel.
Definition: GeometryCore.cxx:271
NOpHardwareChannels
Index NOpHardwareChannels(Index opDet)
Definition: sspmapmaker_v1.c:195
int opdet::DUNE35tonOpDetResponse::doNOpChannels ( ) const
privatevirtual

Reimplemented from opdet::OpDetResponseInterface.

Definition at line 79 of file DUNE35tonOpDetResponse_service.cc.

80  {
81  art::ServiceHandle<geo::Geometry> geom;
82  if (fFastSimChannelConvert || fFullSimChannelConvert)
83  return geom->NOpChannels();
84  else
85  return geom->NOpDets();
86 
87  }
geo::GeometryCore::NOpChannels
unsigned int NOpChannels() const
Number of electronics channels for all the optical detectors.
Definition: GeometryCore.cxx:253
geo::GeometryCore::NOpDets
unsigned int NOpDets() const
Number of OpDets in the whole detector.
Definition: GeometryCore.cxx:244
void opdet::DUNE35tonOpDetResponse::doReconfigure ( fhicl::ParameterSet const &  p)
privatevirtual

Implements opdet::OpDetResponseInterface.

Definition at line 36 of file DUNE35tonOpDetResponse_service.cc.

37  {
38  double tempfQE = pset.get<double>("QuantumEfficiency");
39  fWavelengthCutLow = pset.get<double>("WavelengthCutLow");
40  fWavelengthCutHigh = pset.get<double>("WavelengthCutHigh");
41  fLightGuideAttenuation = pset.get<bool>("LightGuideAttenuation");
42  fChannelConversion = pset.get<std::string>("ChannelConversion");
43  std::string tmpAxis = pset.get<std::string>("LongAxis");
44 
45  boost::algorithm::to_lower(tmpAxis);
46 
47  if (tmpAxis == "x") fLongAxis = 0;
48  if (tmpAxis == "y") fLongAxis = 1;
49  if (tmpAxis == "z") fLongAxis = 2;
50 
51  // Only allow channel conversion once - so it must be set to happen
52  // either during full simulation (library generation) or during
53  // fast simulation (library use).
54 
55  boost::algorithm::to_lower(fChannelConversion);
56 
57  fFullSimChannelConvert = false;
58  fFastSimChannelConvert = false;
59 
60  if (fChannelConversion == "full") fFullSimChannelConvert = true;
61  if (fChannelConversion == "fast") fFastSimChannelConvert = true;
62 
63  // Correct out the prescaling applied during simulation
64  auto const *LarProp = lar::providerFrom<detinfo::LArPropertiesService>();
65  fQE = tempfQE / LarProp->ScintPreScale();
66 
67  if (fQE > 1.0001 ) {
68  mf::LogError("DUNE35tonOpDetResponse_service") << "Quantum efficiency set in OpDetResponse_service, " << tempfQE
69  << " is too large. It is larger than the prescaling applied during simulation, "
70  << LarProp->ScintPreScale()
71  << ". Final QE must be equalt to or smaller than the QE applied at simulation time.";
72  assert(false);
73  }
74 
75  }
string
std::string string
Definition: nybbler.cc:12
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:211

Member Data Documentation

std::string opdet::DUNE35tonOpDetResponse::fChannelConversion
private

Definition at line 44 of file DUNE35tonOpDetResponse.h.

bool opdet::DUNE35tonOpDetResponse::fFastSimChannelConvert
private

Definition at line 46 of file DUNE35tonOpDetResponse.h.

bool opdet::DUNE35tonOpDetResponse::fFullSimChannelConvert
private

Definition at line 45 of file DUNE35tonOpDetResponse.h.

bool opdet::DUNE35tonOpDetResponse::fLightGuideAttenuation
private

Definition at line 42 of file DUNE35tonOpDetResponse.h.

int opdet::DUNE35tonOpDetResponse::fLongAxis
private

Definition at line 48 of file DUNE35tonOpDetResponse.h.

float opdet::DUNE35tonOpDetResponse::fQE
private

Definition at line 37 of file DUNE35tonOpDetResponse.h.

float opdet::DUNE35tonOpDetResponse::fWavelengthCutHigh
private

Definition at line 40 of file DUNE35tonOpDetResponse.h.

float opdet::DUNE35tonOpDetResponse::fWavelengthCutLow
private

Definition at line 39 of file DUNE35tonOpDetResponse.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11