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

#include <ProtoDUNEOpDetResponse.h>

Inheritance diagram for opdet::ProtoDUNEOpDetResponse:
opdet::OpDetResponseInterface

Public Member Functions

 ProtoDUNEOpDetResponse (fhicl::ParameterSet const &pset, art::ActivityRegistry &reg)
 
 ~ProtoDUNEOpDetResponse () throw ()
 
bool detectedLiteWithChannel (int OpDet, int &newOpChannel, int &hardwareChannel) const
 
- 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
 
bool doDetectedLiteWithChannel (int OpDet, int &newOpChannel, int &hardwareChannel) const
 

Private Attributes

float fQE
 
float fQEArapucaBeam
 
float fQEArapucaNonBeam
 
float fWavelengthCutLow
 
float fWavelengthCutHigh
 
bool fLightGuideAttenuation
 
double lambdaShort
 
double lambdaLong
 
double fracShort
 
double fracLong
 
std::string fWireTransmissionFile
 
std::string fChannelConversion
 
bool fFullSimChannelConvert
 
bool fFastSimChannelConvert
 
int fLongAxis
 

Detailed Description

Definition at line 27 of file ProtoDUNEOpDetResponse.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file ProtoDUNEOpDetResponse_service.cc.

25  {
26  this->doReconfigure(pset);
27  }
opdet::ProtoDUNEOpDetResponse::doReconfigure
virtual void doReconfigure(fhicl::ParameterSet const &p)
Definition: ProtoDUNEOpDetResponse_service.cc:35
opdet::ProtoDUNEOpDetResponse::~ProtoDUNEOpDetResponse ( )
throw (
)

Definition at line 30 of file ProtoDUNEOpDetResponse_service.cc.

31  { }

Member Function Documentation

bool opdet::ProtoDUNEOpDetResponse::detectedLiteWithChannel ( int  OpDet,
int &  newOpChannel,
int &  hardwareChannel 
) const
inline

Definition at line 36 of file ProtoDUNEOpDetResponse.h.

37  {
38  return doDetectedLiteWithChannel( OpDet, newOpChannel, hardwareChannel);
39  }
opdet::ProtoDUNEOpDetResponse::doDetectedLiteWithChannel
bool doDetectedLiteWithChannel(int OpDet, int &newOpChannel, int &hardwareChannel) const
Definition: ProtoDUNEOpDetResponse_service.cc:202
bool opdet::ProtoDUNEOpDetResponse::doDetected ( int  OpDet,
const sim::OnePhoton &  Phot,
int &  newOpChannel 
) const
privatevirtual

Implements opdet::OpDetResponseInterface.

Definition at line 99 of file ProtoDUNEOpDetResponse_service.cc.

100  {
101 
102  // Find the Optical Detector using the geometry service
103  art::ServiceHandle<geo::Geometry> geom;
104 
105  if (fFullSimChannelConvert){
106  // Override default number of channels for Fiber and Plank
107  float NOpHardwareChannels = geom->NOpHardwareChannels(OpDet);
108  int hardwareChannel = (int) ( CLHEP::RandFlat::shoot(1.0) * NOpHardwareChannels );
109  newOpChannel = geom->OpChannel(OpDet, hardwareChannel);
110  }else{
111  newOpChannel = OpDet;
112  }
113 
114  // Get the length of the photon detector
115  const TGeoNode* node = geom->OpDetGeoFromOpDet(OpDet).Node();
116  TGeoBBox *box = (TGeoBBox*)node->GetVolume()->GetShape();
117  double opdetLength = 0;
118  double sipmDistance = 0;
119  if (fLongAxis == 0) {
120  opdetLength = box->GetDX();
121  sipmDistance = opdetLength - Phot.FinalLocalPosition.x();
122  }
123  else if (fLongAxis == 1) {
124  opdetLength = box->GetDY();
125  sipmDistance = opdetLength - Phot.FinalLocalPosition.y();
126  }
127  else if (fLongAxis == 2) {
128  opdetLength = box->GetDZ();
129  sipmDistance = opdetLength - Phot.FinalLocalPosition.z();
130  }
131  else {
132  mf::LogError("ProtoDUNEOpDetResponse") << "Unknown axis, fLongAxis = " << fLongAxis;
133  std::abort();
134  }
135 
136  if(opdetLength > 50){ //i.e. if this is a paddle (this is half length in units of cm)
137  // Check QE
138  if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
139 
140  double wavel = wavelength(Phot.Energy);
141  // Check wavelength acceptance
142  if (wavel < fWavelengthCutLow) return false;
143  if (wavel > fWavelengthCutHigh) return false;
144 
145  if (fLightGuideAttenuation) {
146  // Throw away some photons based on attenuation
147  double AttenuationProb = fracShort*exp(-sipmDistance/lambdaShort) + fracLong*exp(-sipmDistance/lambdaLong);
148 
149  if ( CLHEP::RandFlat::shoot(1.0) > AttenuationProb ) return false;
150  }
151  }else{ //if it is Arapuca
152  if(Phot.FinalLocalPosition.x()<0){ //beam side
153  // Check QE
154  if ( CLHEP::RandFlat::shoot(1.0) > fQEArapucaBeam ) return false;
155  }else{ //non-beam side
156  // Check QE
157  if ( CLHEP::RandFlat::shoot(1.0) > fQEArapucaNonBeam ) return false;
158  }
159  }
160  return true;
161  }
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
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
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::ProtoDUNEOpDetResponse::doDetectedLite ( int  OpDet,
int &  newOpChannel 
) const
privatevirtual

Implements opdet::OpDetResponseInterface.

Definition at line 164 of file ProtoDUNEOpDetResponse_service.cc.

165  {
166  // Find the Optical Detector using the geometry service
167  art::ServiceHandle<geo::Geometry> geom;
168 
169  if (fFastSimChannelConvert){
170  // Here OpDet must be opdet since we are introducing
171  // channel mapping here.
172  float NOpHardwareChannels = geom->NOpHardwareChannels(OpDet);
173  int hardwareChannel = (int) ( CLHEP::RandFlat::shoot(1.0) * NOpHardwareChannels );
174  newOpChannel = geom->OpChannel(OpDet, hardwareChannel);
175  }else{
176  newOpChannel = OpDet;
177  }
178 
179  double opdetLength = 0;
180  const TGeoNode* node = geom->OpDetGeoFromOpDet(OpDet).Node();
181  TGeoBBox *box = (TGeoBBox*)node->GetVolume()->GetShape();
182  opdetLength = box->GetDX();
183  if(opdetLength > 50){ //i.e. if this is a paddle (this is half length in units of cm)
184  // Check QE
185  if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
186  }else{ //if it is Arapuca
187  double xyz[3];
188  geom->OpDetGeoFromOpChannel(OpDet).GetCenter(xyz);
189  if(xyz[0]<0){ //beam side
190  // Check QE
191  if ( CLHEP::RandFlat::shoot(1.0) > fQEArapucaBeam ) return false;
192  }else{ //non-beam side
193  // Check QE
194  if ( CLHEP::RandFlat::shoot(1.0) > fQEArapucaNonBeam ) return false;
195  }
196  }
197 
198  return true;
199  }
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
geo::GeometryCore::NOpHardwareChannels
unsigned int NOpHardwareChannels(int opDet) const
Definition: GeometryCore.cxx:265
geo::OpDetGeo::GetCenter
void GetCenter(double *xyz, double localz=0.0) const
Definition: OpDetGeo.cxx:40
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
geo::GeometryCore::OpDetGeoFromOpChannel
OpDetGeo const & OpDetGeoFromOpChannel(unsigned int OpChannel) const
Returns the geo::OpDetGeo object for the given channel number.
Definition: GeometryCore.cxx:1824
bool opdet::ProtoDUNEOpDetResponse::doDetectedLiteWithChannel ( int  OpDet,
int &  newOpChannel,
int &  hardwareChannel 
) const
private

Definition at line 202 of file ProtoDUNEOpDetResponse_service.cc.

203  {
204  // Find the Optical Detector using the geometry service
205  art::ServiceHandle<geo::Geometry> geom;
206 
207  if (fFastSimChannelConvert){
208  // Here OpDet must be opdet since we are introducing
209  // channel mapping here.
210  float NOpHardwareChannels = geom->NOpHardwareChannels(OpDet);
211  hardwareChannel = (int) ( CLHEP::RandFlat::shoot(1.0) * NOpHardwareChannels );
212  newOpChannel = geom->OpChannel(OpDet, hardwareChannel);
213  }else{
214  newOpChannel = OpDet;
215  }
216 
217  double opdetLength = 0;
218  const TGeoNode* node = geom->OpDetGeoFromOpDet(OpDet).Node();
219  TGeoBBox *box = (TGeoBBox*)node->GetVolume()->GetShape();
220  opdetLength = box->GetDX();
221  if(opdetLength > 50){ //i.e. if this is a paddle
222  // Check QE
223  if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
224  }else{ //if it is Arapuca
225  double xyz[3];
226  geom->OpDetGeoFromOpChannel(OpDet).GetCenter(xyz);
227  if(xyz[0]<0){ //beam side
228  // Check QE
229  if ( CLHEP::RandFlat::shoot(1.0) > fQEArapucaBeam ) return false;
230  }else{ //non-beam side
231  // Check QE
232  if ( CLHEP::RandFlat::shoot(1.0) > fQEArapucaNonBeam ) return false;
233  }
234  }
235 
236  return true;
237  }
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
geo::GeometryCore::NOpHardwareChannels
unsigned int NOpHardwareChannels(int opDet) const
Definition: GeometryCore.cxx:265
geo::OpDetGeo::GetCenter
void GetCenter(double *xyz, double localz=0.0) const
Definition: OpDetGeo.cxx:40
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
geo::GeometryCore::OpDetGeoFromOpChannel
OpDetGeo const & OpDetGeoFromOpChannel(unsigned int OpChannel) const
Returns the geo::OpDetGeo object for the given channel number.
Definition: GeometryCore.cxx:1824
int opdet::ProtoDUNEOpDetResponse::doNOpChannels ( ) const
privatevirtual

Reimplemented from opdet::OpDetResponseInterface.

Definition at line 87 of file ProtoDUNEOpDetResponse_service.cc.

88  {
89  art::ServiceHandle<geo::Geometry> geom;
90  if (fFastSimChannelConvert || fFullSimChannelConvert)
91  return geom->NOpChannels();
92  else
93  return geom->NOpDets();
94 
95  }
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::ProtoDUNEOpDetResponse::doReconfigure ( fhicl::ParameterSet const &  p)
privatevirtual

Implements opdet::OpDetResponseInterface.

Definition at line 35 of file ProtoDUNEOpDetResponse_service.cc.

36  {
37  double tempfQE = pset.get<double>("QuantumEfficiency");
38  double tempfQEAraBeam = pset.get<double>("QuantumEfficiencyArapucaBeam");
39  double tempfQEAraNonBeam = pset.get<double>("QuantumEfficiencyArapucaNonBeam");
40  fWavelengthCutLow = pset.get<double>("WavelengthCutLow");
41  fWavelengthCutHigh = pset.get<double>("WavelengthCutHigh");
42  fLightGuideAttenuation = pset.get<bool>("LightGuideAttenuation");
43  lambdaShort = pset.get<double>("LambdaShort");
44  lambdaLong = pset.get<double>("LambdaLong");
45  fracShort = pset.get<double>("FracShort");
46  fracLong = pset.get<double>("FracLong");
47  fChannelConversion = pset.get<std::string>("ChannelConversion");
48  std::string tmpAxis = pset.get<std::string>("LongAxis");
49 
50  boost::algorithm::to_lower(tmpAxis);
51 
52  if (tmpAxis == "x") fLongAxis = 0;
53  if (tmpAxis == "y") fLongAxis = 1;
54  if (tmpAxis == "z") fLongAxis = 2;
55 
56  // Only allow channel conversion once - so it must be set to happen
57  // either during full simulation (library generation) or during
58  // fast simulation (library use).
59 
60  boost::algorithm::to_lower(fChannelConversion);
61 
62  fFullSimChannelConvert = false;
63  fFastSimChannelConvert = false;
64 
65  if (fChannelConversion == "full") fFullSimChannelConvert = true;
66  if (fChannelConversion == "fast") fFastSimChannelConvert = true;
67 
68  // Correct out the prescaling applied during simulation
69  auto const *LarProp = lar::providerFrom<detinfo::LArPropertiesService>();
70  fQE = tempfQE / LarProp->ScintPreScale();
71  fQEArapucaBeam = tempfQEAraBeam / LarProp->ScintPreScale();
72  fQEArapucaNonBeam = tempfQEAraNonBeam / LarProp->ScintPreScale();
73 
74  if (fQE > 1.0001 || fQEArapucaBeam > 1.0001 || fQEArapucaNonBeam > 1.0001) {
75  mf::LogError("ProtoDUNEOpDetResponse_service") << "Quantum efficiency set in OpDetResponse_service, " << tempfQE << ", "
76  << tempfQEAraBeam << ", or " << tempfQEAraNonBeam
77  << " is too large. It is larger than the prescaling applied during simulation, "
78  << LarProp->ScintPreScale()
79  << ". Final QE must be equalt to or smaller than the QE applied at simulation time.";
80  std::abort();
81  }
82 
83  }
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::ProtoDUNEOpDetResponse::fChannelConversion
private

Definition at line 74 of file ProtoDUNEOpDetResponse.h.

bool opdet::ProtoDUNEOpDetResponse::fFastSimChannelConvert
private

Definition at line 76 of file ProtoDUNEOpDetResponse.h.

bool opdet::ProtoDUNEOpDetResponse::fFullSimChannelConvert
private

Definition at line 75 of file ProtoDUNEOpDetResponse.h.

bool opdet::ProtoDUNEOpDetResponse::fLightGuideAttenuation
private

Definition at line 60 of file ProtoDUNEOpDetResponse.h.

int opdet::ProtoDUNEOpDetResponse::fLongAxis
private

Definition at line 78 of file ProtoDUNEOpDetResponse.h.

float opdet::ProtoDUNEOpDetResponse::fQE
private

Definition at line 53 of file ProtoDUNEOpDetResponse.h.

float opdet::ProtoDUNEOpDetResponse::fQEArapucaBeam
private

Definition at line 54 of file ProtoDUNEOpDetResponse.h.

float opdet::ProtoDUNEOpDetResponse::fQEArapucaNonBeam
private

Definition at line 55 of file ProtoDUNEOpDetResponse.h.

double opdet::ProtoDUNEOpDetResponse::fracLong
private

Definition at line 64 of file ProtoDUNEOpDetResponse.h.

double opdet::ProtoDUNEOpDetResponse::fracShort
private

Definition at line 63 of file ProtoDUNEOpDetResponse.h.

float opdet::ProtoDUNEOpDetResponse::fWavelengthCutHigh
private

Definition at line 58 of file ProtoDUNEOpDetResponse.h.

float opdet::ProtoDUNEOpDetResponse::fWavelengthCutLow
private

Definition at line 57 of file ProtoDUNEOpDetResponse.h.

std::string opdet::ProtoDUNEOpDetResponse::fWireTransmissionFile
private

Definition at line 72 of file ProtoDUNEOpDetResponse.h.

double opdet::ProtoDUNEOpDetResponse::lambdaLong
private

Definition at line 62 of file ProtoDUNEOpDetResponse.h.

double opdet::ProtoDUNEOpDetResponse::lambdaShort
private

Definition at line 61 of file ProtoDUNEOpDetResponse.h.

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