#include <DetectorPropertiesProtoDUNEsp.h>

Inheritance diagram for spdp::DetectorPropertiesProtoDUNEsp:

Classes
struct	Configuration_t
	Structure for configuration parameters. More...

struct	SternheimerParameters_t
	Parameters for Sternheimer density effect corrections. More...

Public Types
using	providers_type = lar::ProviderPack< geo::GeometryCore, detinfo::LArProperties >
	List of service providers we depend on. More...

Public Member Functions
	DetectorPropertiesProtoDUNEsp ()

	DetectorPropertiesProtoDUNEsp (fhicl::ParameterSet const &pset, const geo::GeometryCore geo, const detinfo::LArProperties lp, std::set< std::string > const &ignore_params={})

	DetectorPropertiesProtoDUNEsp (fhicl::ParameterSet const &pset, providers_type providers, std::set< std::string > const &ignore_params={})
	Constructs the provider and sets up the dependencies. More...

	DetectorPropertiesProtoDUNEsp (DetectorPropertiesProtoDUNEsp const &)=delete

virtual	~DetectorPropertiesProtoDUNEsp ()=default

void	ValidateAndConfigure (fhicl::ParameterSet const &p, std::set< std::string > const &ignore_params={})
	Configures the provider, first validating the configuration. More...

void	Configure (Configuration_t const &config)
	Extracts the relevant configuration from the specified object. More...

Configuration_t	ValidateConfiguration (fhicl::ParameterSet const &p, std::set< std::string > const &ignore_params={})
	Validates the specified configuration. More...

bool	UpdateHV (std::string filename)

bool	UpdateReadoutWindowSize (detinfo::DetectorClocksData const &clockData, std::string filename)

bool	UpdateTemp (int run)

detinfo::DetectorPropertiesData	DataFor (const detinfo::DetectorClocksData &clockData) const override

void	Setup (providers_type providers)
	Sets all the providers at once. More...

void	SetGeometry (const geo::GeometryCore *g)

void	SetLArProperties (const detinfo::LArProperties *lp)

void	SetNumberTimeSamples (unsigned int nsamp)

virtual double	Efield (unsigned int planegap=0) const override
	kV/cm More...

virtual double	DriftVelocity (double efield=0., double temperature=0.) const override
	cm/us More...

virtual double	BirksCorrection (double dQdX) const override
	dQ/dX in electrons/cm, returns dE/dX in MeV/cm. More...

virtual double	BirksCorrection (double dQdX, double EField) const override

virtual double	ModBoxCorrection (double dQdX) const override

virtual double	ModBoxCorrection (double dQdX, double EField) const override

virtual double	ElectronLifetime () const override
	Returns the attenuation constant for ionization electrons. More...

virtual double	Density (double temperature) const override
	Returns argon density at a given temperature. More...

virtual double	Density () const override
	Returns argon density at the temperature from Temperature() More...

virtual double	Temperature () const override
	In kelvin. More...

virtual double	Eloss (double mom, double mass, double tcut) const override
	Restricted mean energy loss (dE/dx) More...

virtual double	ElossVar (double mom, double mass) const override
	Energy loss fluctuation ( $\sigma_{E}^2 / x$ ) More...

virtual double	ElectronsToADC () const override

virtual unsigned int	NumberTimeSamples () const override

virtual unsigned int	ReadOutWindowSize () const override

virtual double	TimeOffsetU () const override

virtual double	TimeOffsetV () const override

virtual double	TimeOffsetZ () const override

virtual double	TimeOffsetY () const override

virtual bool	SimpleBoundary () const override

void	CheckIfConfigured () const

Public Member Functions inherited from detinfo::DetectorProperties
	DetectorProperties ()=default

	DetectorProperties (const DetectorProperties &)=delete

	DetectorProperties (DetectorProperties &&)=delete

DetectorProperties &	operator= (const DetectorProperties &)=delete

DetectorProperties &	operator= (DetectorProperties &&)=delete

virtual	~DetectorProperties ()=default

virtual DetectorPropertiesData	DataFor (DetectorClocksData const &clockData) const =0

Protected Member Functions
detinfo::DetectorPropertiesData	CalculateXTicksParams (detinfo::DetectorClocksData const &clockData) const

std::string	CheckTimeOffsetConfigurationAfterSetup () const
	Checks the configuration of time offsets. More...

void	CheckConfigurationAfterSetup () const

Protected Attributes
const detinfo::LArProperties *	fLP

const geo::GeometryCore *	fGeo

bool	fGetHVDriftfromMetaData

bool	fGetReadOutWindowSizefromMetaData

bool	fUseRunDependentTemperature

double	fHV_cath

std::vector< double >	fEfield
	kV/cm (per inter-plane volume) More...

double	fElectronlifetime
	microseconds More...

double	fTemperature
	kelvin More...

double	fSamplingRate
	in ns More...

double	fElectronsToADC
	conversion factor for # of ionization electrons to 1 ADC count More...

unsigned int	fNumberTimeSamples
	number of clock ticks per event More...

unsigned int	fReadOutWindowSize
	number of clock ticks per readout window More...

double	fTimeOffsetU
	time offset to convert spacepoint coordinates to hit times on view U More...

double	fTimeOffsetV
	time offset to convert spacepoint coordinates to hit times on view V More...

double	fTimeOffsetZ
	time offset to convert spacepoint coordinates to hit times on view Z More...

double	fTimeOffsetY
	time offset to convert spacepoint coordinates to hit times on view Y More...

double	fTimeOffsetX
	time offset to convert spacepoint coordinates to hit times on view X More...

double	fHasTimeOffsetU = false
	whether time offset was configured for view U More...

double	fHasTimeOffsetV = false
	whether time offset was configured for view V More...

double	fHasTimeOffsetZ = false
	whether time offset was configured for view Z More...

double	fHasTimeOffsetY = false
	whether time offset was configured for view Y More...

double	fHasTimeOffsetX = false
	whether time offset was configured for view X More...

SternheimerParameters_t	fSternheimerParameters
	Sternheimer parameters. More...

bool	fSimpleBoundary

Detailed Description

Definition at line 59 of file DetectorPropertiesProtoDUNEsp.h.

Member Typedef Documentation

using spdp::DetectorPropertiesProtoDUNEsp::providers_type = lar::ProviderPack< geo::GeometryCore, detinfo::LArProperties >

List of service providers we depend on.

Definition at line 65 of file DetectorPropertiesProtoDUNEsp.h.

Constructor & Destructor Documentation

spdp::DetectorPropertiesProtoDUNEsp::DetectorPropertiesProtoDUNEsp ( )

Definition at line 56 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                :
     fLP(0), fGeo(0)
   {
 
   }

spdp::DetectorPropertiesProtoDUNEsp::DetectorPropertiesProtoDUNEsp	(	fhicl::ParameterSet const &	pset,
		const geo::GeometryCore *	geo,
		const detinfo::LArProperties *	lp,
		std::set< std::string > const &	ignore_params = `{}`
	)

Definition at line 63 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                 :
     fLP(lp), fGeo(geo)
   {
     {
       mf::LogInfo debug("setupProvider<DetectorPropertiesStandard>");
 
       debug << "Asked to ignore " << ignore_params.size() << " keys:";
       for (auto const& key: ignore_params) debug << " '" << key << "'";
     }
 
     ValidateAndConfigure(pset, ignore_params);
 
   }

spdp::DetectorPropertiesProtoDUNEsp::DetectorPropertiesProtoDUNEsp	(	fhicl::ParameterSet const &	pset,
		providers_type	providers,
		std::set< std::string > const &	ignore_params = `{}`
	)

Constructs the provider and sets up the dependencies.

Parameters

pset	FHiCL parameter set for provider configuration
providers	pack of providers `DetectorPropertiesStandard` depends on
ignore_params	unknown configuration keys in `pset` to be tolerated

See also: Setup()

Definition at line 82 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                 :
     DetectorPropertiesProtoDUNEsp(pset,
                                   providers.get<geo::GeometryCore>(),
                                   providers.get<detinfo::LArProperties>(),
                                   ignore_params
                                   )
   {}

spdp::DetectorPropertiesProtoDUNEsp::DetectorPropertiesProtoDUNEsp ( DetectorPropertiesProtoDUNEsp const & )

delete

virtual spdp::DetectorPropertiesProtoDUNEsp::~DetectorPropertiesProtoDUNEsp ( )

virtualdefault

Member Function Documentation

double spdp::DetectorPropertiesProtoDUNEsp::BirksCorrection ( double dQdX ) const

overridevirtual

dQ/dX in electrons/cm, returns dE/dX in MeV/cm.

Implements detinfo::DetectorProperties.

Definition at line 482 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     return BirksCorrection(dQdx, Efield());
   }

double spdp::DetectorPropertiesProtoDUNEsp::BirksCorrection	(	double	dQdX,
		double	EField
	)		const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 486 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     // Correction for charge quenching using parameterization from
     // S.Amoruso et al., NIM A 523 (2004) 275
 
     double  A3t    = util::kRecombA;
     double  K3t    = util::kRecombk;                     // in KV/cm*(g/cm^2)/MeV
     double  rho    = Density();                    // LAr density in g/cm^3
     double Wion    = 1000./util::kGeVToElectrons;        // 23.6 eV = 1e, Wion in MeV/e
     K3t           /= rho;                                // KV/MeV
     double dEdx    = dQdx/(A3t/Wion-K3t/E_field*dQdx);    //MeV/cm
 
     return dEdx;
   }

detinfo::DetectorPropertiesData spdp::DetectorPropertiesProtoDUNEsp::CalculateXTicksParams ( detinfo::DetectorClocksData const & clockData ) const

protected

Definition at line 533 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     CheckIfConfigured();
 
     double samplingRate   = sampling_rate(clockData);
     double efield         = Efield();
     double temperature    = Temperature();
     double driftVelocity  = DriftVelocity(efield, temperature);
 
     double const x_ticks_coefficient = 0.001 * driftVelocity * samplingRate;
     double triggerOffset  = trigger_offset(clockData);
     std::vector<std::vector<std::vector<double>>> x_ticks_offsets(fGeo->Ncryostats());
     std::vector<std::vector<double>> drift_direction(fGeo->Ncryostats());
     for(size_t cstat = 0; cstat < fGeo->Ncryostats(); ++cstat){
       x_ticks_offsets[cstat].resize(fGeo->Cryostat(cstat).NTPC());
       drift_direction[cstat].resize(fGeo->Cryostat(cstat).NTPC());
       for(size_t tpc = 0; tpc < fGeo->Cryostat(cstat).NTPC(); ++tpc) {
         const geo::TPCGeo& tpcgeom = fGeo->Cryostat(cstat).TPC(tpc);
         const double dir((tpcgeom.DriftDirection() == geo::kNegX) ? +1.0 :-1.0);
         drift_direction[cstat][tpc] = dir;
         int nplane = tpcgeom.Nplanes();
         x_ticks_offsets[cstat][tpc].resize(nplane, 0.);
 
         for(int plane = 0; plane < nplane; ++plane) {
           const geo::PlaneGeo& pgeom = tpcgeom.Plane(plane);
           const double* xyz = tpcgeom.PlaneLocation(0);
           x_ticks_offsets[cstat][tpc][plane] = -xyz[0]/(dir * x_ticks_coefficient) + triggerOffset;
 
           // Add view dependent offset
           // FIXME the offset should be plane-dependent
           geo::View_t view = pgeom.View();
           switch (view) {
           case geo::kU:
             x_ticks_offsets[cstat][tpc][plane] += fTimeOffsetU;
             break;
           case geo::kV:
             x_ticks_offsets[cstat][tpc][plane] += fTimeOffsetV;
             break;
           case geo::kZ:
             x_ticks_offsets[cstat][tpc][plane] += fTimeOffsetZ;
             break;
           case geo::kY:
             x_ticks_offsets[cstat][tpc][plane] += fTimeOffsetY;
             break;
           case geo::kX:
             x_ticks_offsets[cstat][tpc][plane] += fTimeOffsetX;
             break;
           default:
             throw cet::exception(__FUNCTION__) << "Bad view = " << view << "\n" ;
           } // switch
         }
       }
     }
     return detinfo::DetectorPropertiesData{*this, x_ticks_coefficient, move(x_ticks_offsets), move(drift_direction)};
   }

void spdp::DetectorPropertiesProtoDUNEsp::CheckConfigurationAfterSetup ( ) const

protected

Checks that provider configuration is complete, using setup information.

Definition at line 633 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                          {
 
     std::string errors;
 
     errors += CheckTimeOffsetConfigurationAfterSetup();
 
     if (!errors.empty()) {
       throw cet::exception("DetectorPropertiesStandard")
         << "Detected configuration errors: \n" << errors;
     }
 
   } // DetectorPropertiesStandard::CheckConfigurationAfterSetup()

void spdp::DetectorPropertiesProtoDUNEsp::CheckIfConfigured ( ) const

Verifies that the provider is in a fully configured status

Exceptions

cet::exception (category DetectorPropertiesStandard) if not ok

Definition at line 522 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     if (!fGeo) throw cet::exception(__FUNCTION__) << "Geometry is uninitialized!";
     if (!fLP) throw cet::exception(__FUNCTION__) << "LArPropertiesStandard is uninitialized!";
   }

std::string spdp::DetectorPropertiesProtoDUNEsp::CheckTimeOffsetConfigurationAfterSetup ( ) const

protected

Checks the configuration of time offsets.

Definition at line 591 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
 
     std::ostringstream errors;
     auto const views = fGeo->Views();
 
     if ((views.count(geo::kU) != 0) != fHasTimeOffsetU) {
       if (fHasTimeOffsetU)
         errors << "TimeOffsetU has been specified, but no U view is present.\n";
       else
         errors << "TimeOffsetU missing for view U.\n";
     }
     if ((views.count(geo::kV) != 0) != fHasTimeOffsetV) {
       if (fHasTimeOffsetV)
         errors << "TimeOffsetV has been specified, but no V view is present.\n";
       else
         errors << "TimeOffsetV missing for view Z.\n";
     }
     if ((views.count(geo::kZ) != 0) != fHasTimeOffsetZ) {
       if (fHasTimeOffsetZ)
         errors << "TimeOffsetZ has been specified, but no Z view is present.\n";
       else
         errors << "TimeOffsetZ missing for view Z.\n";
     }
     if ((views.count(geo::kY) != 0) != fHasTimeOffsetY) {
       if (fHasTimeOffsetY)
         errors << "TimeOffsetY has been specified, but no Y view is present.\n";
       else
         errors << "TimeOffsetY missing for view Y.\n";
     }
     if ((views.count(geo::kX) != 0) != fHasTimeOffsetX) {
       if (fHasTimeOffsetX)
         errors << "TimeOffsetX has been specified, but no X view is present.\n";
       else
         errors << "TimeOffsetX missing for view X.\n";
     }
 
     return errors.str();
 
   } // DetectorPropertiesStandard::CheckTimeOffsetConfigurationAfterSetup()

void spdp::DetectorPropertiesProtoDUNEsp::Configure ( Configuration_t const & config )

Extracts the relevant configuration from the specified object.

Definition at line 229 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                              {
 
     fEfield                     = config.Efield();
     fGetHVDriftfromMetaData    = config.fGetHVDriftfromMetaData();
     fGetHVDriftfromMetaData    = false;  // DLA 2021-11-11  Redmine 26419
     fGetReadOutWindowSizefromMetaData = config.fGetReadOutWindowSizefromMetaData();
     fUseRunDependentTemperature = config.fUseRunDependentTemperature();
     fElectronlifetime           = config.Electronlifetime();
     fTemperature                = config.Temperature();
     fElectronsToADC             = config.ElectronsToADC();
     fNumberTimeSamples          = config.NumberTimeSamples();
     fReadOutWindowSize          = config.ReadOutWindowSize();
     fHasTimeOffsetU = config.TimeOffsetU(fTimeOffsetU);
     fHasTimeOffsetV = config.TimeOffsetV(fTimeOffsetV);
     fHasTimeOffsetZ = config.TimeOffsetZ(fTimeOffsetZ);
     fHasTimeOffsetY = config.TimeOffsetY(fTimeOffsetY);
     fHasTimeOffsetX = config.TimeOffsetX(fTimeOffsetX);
 
     fSternheimerParameters.a    = config.SternheimerA();
     fSternheimerParameters.k    = config.SternheimerK();
     fSternheimerParameters.x0   = config.SternheimerX0();
     fSternheimerParameters.x1   = config.SternheimerX1();
     fSternheimerParameters.cbar = config.SternheimerCbar();
     fSimpleBoundary = config.SimpleBoundary();
 
   } // DetectorPropertiesStandard::Configure()

detinfo::DetectorPropertiesData spdp::DetectorPropertiesProtoDUNEsp::DataFor ( const detinfo::DetectorClocksData & clockData ) const

override

Definition at line 223 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     return CalculateXTicksParams(clockData);
   }

double spdp::DetectorPropertiesProtoDUNEsp::Density ( double temperature ) const

overridevirtual

Returns argon density at a given temperature.

Parameters

temperature the temperature in kelvin

Returns: argon density in g/cm^3

Density is nearly a linear function of temperature. See the NIST tables for details Slope is between -6.2 and -6.1, intercept is 1928 kg/m^3. This parameterization will be good to better than 0.5%.g/cm^3

Implements detinfo::DetectorProperties.

Definition at line 304 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     // Default temperature use internal value.
     if(temperature == 0.)
       temperature = Temperature();
 
     double density = -0.00615*temperature + 1.928;
 
     return density;
   } // DetectorPropertiesStandard::Density()

virtual double spdp::DetectorPropertiesProtoDUNEsp::Density ( ) const

inlineoverridevirtual

Returns argon density at the temperature from Temperature()

Reimplemented from detinfo::DetectorProperties.

Definition at line 261 of file DetectorPropertiesProtoDUNEsp.h.

261 { return Density(Temperature()); }

spdp::DetectorPropertiesProtoDUNEsp::Temperature

virtual double Temperature() const override

In kelvin.

Definition: DetectorPropertiesProtoDUNEsp.h:264

spdp::DetectorPropertiesProtoDUNEsp::Density

virtual double Density() const override

Returns argon density at the temperature from Temperature()

Definition: DetectorPropertiesProtoDUNEsp.h:261

double spdp::DetectorPropertiesProtoDUNEsp::DriftVelocity	(	double	efield = `0.`,
		double	temperature = `0.`
	)		const

overridevirtual

cm/us

Implements detinfo::DetectorProperties.

Definition at line 402 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                                              {
     // Drift Velocity as a function of Electric Field and LAr Temperature
     // from : W. Walkowiak, NIM A 449 (2000) 288-294
     //
     // Efield should have units of kV/cm
     // Temperature should have units of Kelvin
     // Default Efield, use internal value.
     if(efield == 0.)
       efield = Efield();
     //
     if(efield > 4.0)
       mf::LogWarning("DetectorPropertiesStandard") << "DriftVelocity Warning! : E-field value of "
                                                    << efield
                                                    << " kV/cm is outside of range covered by drift"
                                                    << " velocity parameterization. Returned value"
                                                    << " may not be correct";
     // Default temperature use internal value.
     if(temperature == 0.)
       temperature = Temperature();
     if(temperature < 87.0 || temperature > 94.0)
       mf::LogWarning("DetectorPropertiesStandard") << "DriftVelocity Warning! : Temperature value of "
                                                    << temperature
                                                    << " K is outside of range covered by drift velocity"
                                                    << " parameterization. Returned value may not be"
                                                    << " correct";
     double tshift = -87.203+temperature;
     double xFit = 0.0938163-0.0052563*tshift-0.0001470*tshift*tshift;
     double uFit = 5.18406+0.01448*tshift-0.003497*tshift*tshift-0.000516*tshift*tshift*tshift;
     double vd;
     // Icarus Parameter Set, use as default
     double  P1 = -0.04640; // K^-1
     double  P2 = 0.01712;  // K^-1
     double  P3 = 1.88125;   // (kV/cm)^-1
     double  P4 =  0.99408;    // kV/cm
     double  P5 =  0.01172;   // (kV/cm)^-P6
     double  P6 =  4.20214;
     double  T0 =  105.749;  // K
     // Walkowiak Parameter Set
     double    P1W = -0.01481; // K^-1
     double  P2W = -0.0075;  // K^-1
     double   P3W =  0.141;   // (kV/cm)^-1
     double   P4W =  12.4;    // kV/cm
     double   P5W =  1.627;   // (kV/cm)^-P6
     double   P6W =  0.317;
     double   T0W =  90.371;  // K
     // From Craig Thorne . . . currently not documented
     // smooth transition from linear at small fields to
     //     icarus fit at most fields to Walkowiak at very high fields
     if (efield < xFit) vd=efield*uFit;
     else if (efield<0.619) {
       vd = ((P1*(temperature-T0)+1)
             *(P3*efield*std::log(1+P4/efield) + P5*std::pow(efield,P6))
             +P2*(temperature-T0));
     }
     else if (efield<0.699) {
       vd = 12.5*(efield-0.619)*((P1W*(temperature-T0W)+1)
                                 *(P3W*efield*std::log(1+P4W/efield) + P5W*std::pow(efield,P6W))
                                 +P2W*(temperature-T0W))+
         12.5*(0.699-efield)*((P1*(temperature-T0)+1)
                              *(P3*efield*std::log(1+P4/efield) + P5*std::pow(efield,P6))
                              +P2*(temperature-T0));
     }
     else {
       vd = ((P1W*(temperature-T0W)+1)
             *(P3W*efield*std::log(1+P4W/efield) + P5W*std::pow(efield,P6W))
             +P2W*(temperature-T0W));
     }
     vd /= 10.;
     return vd; // in cm/us
   }

double spdp::DetectorPropertiesProtoDUNEsp::Efield ( unsigned int planegap = 0 ) const

overridevirtual

kV/cm

Implements detinfo::DetectorProperties.

Definition at line 293 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     if(planegap >= fEfield.size())
       throw cet::exception("DetectorPropertiesStandard") << "requesting Electric field in a plane gap that is not defined\n";
 
 
     return fEfield[planegap];
   }

virtual double spdp::DetectorPropertiesProtoDUNEsp::ElectronLifetime ( ) const

inlineoverridevirtual

Returns the attenuation constant for ionization electrons.

Returns: the attenuation constant [µs]

The returned constant τ can be used to know the attenuation the ionization charge undergoes after drifting for a certain time t: $A(t) = e^{-t/\tau}$ (t is measured in microseconds).

This is a uniform, constant value for the detector.

Implements detinfo::DetectorProperties.

Definition at line 245 of file DetectorPropertiesProtoDUNEsp.h.

245 { return fElectronlifetime; } //< microseconds

spdp::DetectorPropertiesProtoDUNEsp::fElectronlifetime

double fElectronlifetime

microseconds

Definition: DetectorPropertiesProtoDUNEsp.h:333

virtual double spdp::DetectorPropertiesProtoDUNEsp::ElectronsToADC ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 291 of file DetectorPropertiesProtoDUNEsp.h.

291 { return fElectronsToADC; }

spdp::DetectorPropertiesProtoDUNEsp::fElectronsToADC

double fElectronsToADC

conversion factor for # of ionization electrons to 1 ADC count

Definition: DetectorPropertiesProtoDUNEsp.h:336

double spdp::DetectorPropertiesProtoDUNEsp::Eloss	(	double	mom,
		double	mass,
		double	tcut
	)		const

overridevirtual

Restricted mean energy loss (dE/dx)

Parameters

mom	momentum of incident particle [GeV/c]
mass	mass of incident particle [GeV/c^2]
tcut	maximum kinetic energy of delta rays [MeV]; 0 for unlimited

Returns: the restricted mean energy loss (dE/dx) in units of MeV/cm

Returned value is always positive. For unrestricted mean energy loss, set tcut = 0 (special case), or tcut large.

Based on Bethe-Bloch formula as contained in particle data book. Material parameters are from the configuration.

Implements detinfo::DetectorProperties.

Definition at line 333 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     // Some constants.
 
     double K = 0.307075;     // 4 pi N_A r_e^2 m_e c^2 (MeV cm^2/mol).
     double me = 0.510998918; // Electron mass (MeV/c^2).
 
     // Calculate kinematic quantities.
 
     double bg = mom / mass;           // beta*gamma.
     double gamma = sqrt(1. + bg*bg);  // gamma.
     double beta = bg / gamma;         // beta (velocity).
     double mer = 0.001 * me / mass;   // electron mass / mass of incident particle.
     double tmax = 2.*me* bg*bg / (1. + 2.*gamma*mer + mer*mer);  // Maximum delta ray energy (MeV).
 
     // Make sure tcut does not exceed tmax.
 
     if(tcut == 0. || tcut > tmax)
       tcut = tmax;
 
     // Calculate density effect correction (delta).
 
     double x = std::log10(bg);
     double delta = 0.;
     if(x >= fSternheimerParameters.x0) {
       delta = 2. * std::log(10.) * x - fSternheimerParameters.cbar;
       if(x < fSternheimerParameters.x1)
         delta += fSternheimerParameters.a * std::pow(fSternheimerParameters.x1 - x, fSternheimerParameters.k);
     }
 
     // Calculate stopping number.
 
     double B = 0.5 * std::log(2.*me*bg*bg*tcut / (1.e-12 * sqr(fLP->ExcitationEnergy())))
       - 0.5 * beta*beta * (1. + tcut / tmax) - 0.5 * delta;
 
     // Don't let the stopping number become negative.
 
     if(B < 1.)
       B = 1.;
 
     // Calculate dE/dx.
 
     double dedx = Density() * K*fLP->AtomicNumber()*B / (fLP->AtomicMass() * beta*beta);
 
     // Done.
 
     return dedx;
   } // DetectorPropertiesStandard::Eloss()

double spdp::DetectorPropertiesProtoDUNEsp::ElossVar	(	double	mom,
		double	mass
	)		const

overridevirtual

Energy loss fluctuation ( $\sigma_{E}^2 / x$ )

Parameters

mom	momentum of incident particle in [GeV/c]
mass	mass of incident particle [GeV/c^2]

Returns: energy loss fluctuation in MeV^2/cm

Based on Bichsel formula referred to but not given in PDG.

Implements detinfo::DetectorProperties.

Definition at line 383 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     // Some constants.
 
     double K = 0.307075;     // 4 pi N_A r_e^2 m_e c^2 (MeV cm^2/mol).
     double me = 0.510998918; // Electron mass (MeV/c^2).
 
     // Calculate kinematic quantities.
 
     double bg = mom / mass;          // beta*gamma.
     double gamma2 = 1. + bg*bg;      // gamma^2.
     double beta2 = bg*bg / gamma2;   // beta^2.
 
     // Calculate final result.
 
     double result = gamma2 * (1. - 0.5 * beta2) * me * (fLP->AtomicNumber() / fLP->AtomicMass()) * K * Density();
     return result;
   } // DetectorPropertiesStandard::ElossVar()

double spdp::DetectorPropertiesProtoDUNEsp::ModBoxCorrection ( double dQdX ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 503 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     return ModBoxCorrection(dQdx, Efield());
   }

double spdp::DetectorPropertiesProtoDUNEsp::ModBoxCorrection	(	double	dQdX,
		double	EField
	)		const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 507 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
     // Modified Box model correction has better behavior than the Birks
     // correction at high values of dQ/dx.
     double  rho    = Density();                    // LAr density in g/cm^3
     double Wion    = 1000./util::kGeVToElectrons;        // 23.6 eV = 1e, Wion in MeV/e
     double Beta    = util::kModBoxB / (rho * E_field);
     double Alpha   = util::kModBoxA;
     double dEdx = (exp(Beta * Wion * dQdx ) - Alpha) / Beta;
 
     return dEdx;
 
   }

virtual unsigned int spdp::DetectorPropertiesProtoDUNEsp::NumberTimeSamples ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 292 of file DetectorPropertiesProtoDUNEsp.h.

292 { return fNumberTimeSamples; }

spdp::DetectorPropertiesProtoDUNEsp::fNumberTimeSamples

unsigned int fNumberTimeSamples

number of clock ticks per event

Definition: DetectorPropertiesProtoDUNEsp.h:337

virtual unsigned int spdp::DetectorPropertiesProtoDUNEsp::ReadOutWindowSize ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 293 of file DetectorPropertiesProtoDUNEsp.h.

293 { return fReadOutWindowSize; }

spdp::DetectorPropertiesProtoDUNEsp::fReadOutWindowSize

unsigned int fReadOutWindowSize

number of clock ticks per readout window

Definition: DetectorPropertiesProtoDUNEsp.h:338

void spdp::DetectorPropertiesProtoDUNEsp::SetGeometry ( const geo::GeometryCore * g )

inline

Definition at line 233 of file DetectorPropertiesProtoDUNEsp.h.

233 { fGeo = g; }

genie::units::g

static constexpr double g

Definition: Units.h:144

spdp::DetectorPropertiesProtoDUNEsp::fGeo

const geo::GeometryCore * fGeo

Definition: DetectorPropertiesProtoDUNEsp.h:324

void spdp::DetectorPropertiesProtoDUNEsp::SetLArProperties ( const detinfo::LArProperties * lp )

inline

Definition at line 234 of file DetectorPropertiesProtoDUNEsp.h.

234 { fLP = lp; }

spdp::DetectorPropertiesProtoDUNEsp::fLP

const detinfo::LArProperties * fLP

Definition: DetectorPropertiesProtoDUNEsp.h:323

void spdp::DetectorPropertiesProtoDUNEsp::SetNumberTimeSamples ( unsigned int nsamp )

inline

Definition at line 235 of file DetectorPropertiesProtoDUNEsp.h.

235 { fNumberTimeSamples=nsamp;}

spdp::DetectorPropertiesProtoDUNEsp::fNumberTimeSamples

unsigned int fNumberTimeSamples

number of clock ticks per event

Definition: DetectorPropertiesProtoDUNEsp.h:337

void spdp::DetectorPropertiesProtoDUNEsp::Setup ( providers_type providers )

Sets all the providers at once.

Parameters

providers the pack of service providers we depend on

Example:

lar::DetectorPropertiesStandard::providers_type providers;
providers.set(lar::providerFrom<geo::Geometry>());
providers.set(lar::providerFrom<detinfo::LArPropertiesService>());
providers.set(lar::providerFrom<detinfo::DetectorClocksService>());
detprop->Setup(providers);

Definition at line 282 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                     {
 
     SetGeometry(providers.get<geo::GeometryCore>());
     SetLArProperties(providers.get<detinfo::LArProperties>());
 
     CheckConfigurationAfterSetup();
 
   } // DetectorPropertiesStandard::Setup()

virtual bool spdp::DetectorPropertiesProtoDUNEsp::SimpleBoundary ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 301 of file DetectorPropertiesProtoDUNEsp.h.

301 { return fSimpleBoundary; }

spdp::DetectorPropertiesProtoDUNEsp::fSimpleBoundary

bool fSimpleBoundary

Definition: DetectorPropertiesProtoDUNEsp.h:352

virtual double spdp::DetectorPropertiesProtoDUNEsp::Temperature ( ) const

inlineoverridevirtual

In kelvin.

Implements detinfo::DetectorProperties.

Definition at line 264 of file DetectorPropertiesProtoDUNEsp.h.

264 { return fTemperature; }

spdp::DetectorPropertiesProtoDUNEsp::fTemperature

double fTemperature

kelvin

Definition: DetectorPropertiesProtoDUNEsp.h:334

virtual double spdp::DetectorPropertiesProtoDUNEsp::TimeOffsetU ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 294 of file DetectorPropertiesProtoDUNEsp.h.

294 { return fTimeOffsetU; };

spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetU

double fTimeOffsetU

time offset to convert spacepoint coordinates to hit times on view U

Definition: DetectorPropertiesProtoDUNEsp.h:339

virtual double spdp::DetectorPropertiesProtoDUNEsp::TimeOffsetV ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 295 of file DetectorPropertiesProtoDUNEsp.h.

295 { return fTimeOffsetV; };

spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetV

double fTimeOffsetV

time offset to convert spacepoint coordinates to hit times on view V

Definition: DetectorPropertiesProtoDUNEsp.h:340

virtual double spdp::DetectorPropertiesProtoDUNEsp::TimeOffsetY ( ) const

inlineoverridevirtual

Reimplemented from detinfo::DetectorProperties.

Definition at line 297 of file DetectorPropertiesProtoDUNEsp.h.

297 { return fTimeOffsetY; };

spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetY

double fTimeOffsetY

time offset to convert spacepoint coordinates to hit times on view Y

Definition: DetectorPropertiesProtoDUNEsp.h:342

virtual double spdp::DetectorPropertiesProtoDUNEsp::TimeOffsetZ ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 296 of file DetectorPropertiesProtoDUNEsp.h.

296 { return fTimeOffsetZ; };

spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetZ

double fTimeOffsetZ

time offset to convert spacepoint coordinates to hit times on view Z

Definition: DetectorPropertiesProtoDUNEsp.h:341

bool spdp::DetectorPropertiesProtoDUNEsp::UpdateHV ( std::string filename )

Definition at line 146 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
 
 
 
     bool retVal=false;
 
     if(fGetHVDriftfromMetaData){
       retVal = true;
 
       art::ServiceHandle<ifdh_ns::IFDH> ifdh;
       auto metadata=ifdh->getMetadata(filename);
 
       int run = 0;
       std::string run_str="Runs: ";
       if(metadata.find(run_str)!=std::string::npos){
         int n1 = metadata.find(run_str);
         n1 += run_str.length();
         int n2 = metadata.find(".", n1);
         run = std::stoi(metadata.substr(n1, n2-n1));
         std::cout<<"Run number from metadata: "<<run<<std::endl;
       }
 
       std::string hv_str="detector.hv_value: ";
       if(metadata.find(hv_str)!=std::string::npos){
         int n1 = metadata.find(hv_str);
         n1 += hv_str.length();
         int n2 = metadata.find("\n", n1);
 
         fHV_cath=std::stod(metadata.substr(n1, n2-n1));
         std::cout<<"Using HV on cathode as: "<<fHV_cath<<"KV,  Value retreived from samweb MetaData"<<std::endl;
       }
       else{
         std::cout<<"HV Metadata entry not found"<<std::endl;
         std::cout<<"Run has no nomninal HV value entered in Metadata, default to 180KV"<<std::endl;
         fHV_cath=180;
       }
 
 
 
       double Gplane_bias=(-0.665); //
       double Uplane_bias=(-0.370);
       double Vplane_bias=0;
       double Xplane_bias=(0.820);
 
 
       //double Gplane_bias=(-0.665*(fHV_cath/180)); // this is scaling, wasnt actually done.
       //double Uplane_bias=(-0.370*(fHV_cath/180));
       //double Vplane_bias=0;
       //double Xplane_bias=0.820*(fHV_cath/180);
       fEfield={fHV_cath/360,std::fabs(Gplane_bias-Uplane_bias)/0.47625,std::fabs(Uplane_bias-Vplane_bias)/0.47625,std::fabs(Vplane_bias-Xplane_bias)/0.47625};
       //fEfield[0]=fHV_cath/360;
       //fEfield[1]=std::fabs(Gplane_bias-Uplane_bias)/0.475;
       //fEfield[2]=std::fabs(Uplane_bias-Vplane_bias)/0.475;
       //fEfield[3]=std::fabs(Vplane_bias-Xplane_bias)/0.475;
       if (std::abs(fHV_cath-180) < 1e-6){
         // Use the corrected E field from Flavio, Francesco and Stefania
         // Ref: https://indico.fnal.gov/event/20939/contribution/1/material/slides/0.pdf
         if (run < 6725){ //first run in Feb 8, 2019
           fEfield[0] = 0.4867;
         }
         else{
           fEfield[0] = 0.4995;
         }
       }
       std::cout<<"Calculated E field in 4 plane gaps as: "<<fEfield[0]<<","<<fEfield[1]<<","<<fEfield[2]<<","<<fEfield[3]<<std::endl;
     }//End GetHVDriftfrom MetaData if
 
 
     return retVal;
   }

bool spdp::DetectorPropertiesProtoDUNEsp::UpdateReadoutWindowSize	(	detinfo::DetectorClocksData const &	clockData,
		std::string	filename
	)

Definition at line 112 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                                                  {
 
     bool retVal=false;
 
     if(fGetReadOutWindowSizefromMetaData){
       art::ServiceHandle<ifdh_ns::IFDH> ifdh;
       auto metadata=ifdh->getMetadata(filename);
       retVal = true;
       std::string window_str="DUNE_data.readout_window: ";
       if(metadata.find(window_str)!=std::string::npos){
         int n1 = metadata.find(window_str);
         n1 += window_str.length();
         int n2 = metadata.find("\n", n1);
         double window=std::stod(metadata.substr(n1, n2-n1)); //milliseconds
         double ticks=window*1000/(clockData.TPCClock().TickPeriod()); //sampling rate 2Mhz
         fNumberTimeSamples=ticks;
         fReadOutWindowSize=ticks;
         std::cout<<"Setting ReadOutWindowSize and NumberTimeSamples as : "<<fReadOutWindowSize<<" (Ticks), Value retrieved from Metadata"<<std::endl;
       }
       else{
         std::cout<<"Run has no specific readout window size entry in Metadata, default to fcl configured value : "<<fReadOutWindowSize<<" (Ticks)"<<std::endl;
       }
     }
 
     return retVal;
 
   }

bool spdp::DetectorPropertiesProtoDUNEsp::UpdateTemp ( int run )

Definition at line 93 of file DetectorPropertiesProtoDUNEsp.cxx.

   {
 
 
     if(fUseRunDependentTemperature){ // updata the temperature based on run number for data
 
       if ((run > 5903) & (run < 6930)){ //first runs after Nov 17 2019 and March 1 2019, where tempreture average was lower
         fTemperature = 87.36;
       }
       else if (run >= 6930)
       {
         fTemperature = 87.65;
       }
     }
 
 
     return true;
   }

void spdp::DetectorPropertiesProtoDUNEsp::ValidateAndConfigure	(	fhicl::ParameterSet const &	p,
		std::set< std::string > const &	ignore_params = `{}`
	)

Configures the provider, first validating the configuration.

Parameters

p	configuration parameter set
ignore_params	parameters to be ignored (optional)

This method will validate the parameter set (except for the parameters it's explicitly told to ignore) and extract the useful information out of it.

Definition at line 273 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                              {
     Configure(ValidateConfiguration(p, ignore_params));
   } // ValidateAndConfigure()

DetectorPropertiesProtoDUNEsp::Configuration_t spdp::DetectorPropertiesProtoDUNEsp::ValidateConfiguration	(	fhicl::ParameterSet const &	p,
		std::set< std::string > const &	ignore_params = `{}`
	)

Validates the specified configuration.

Parameters

p	configuration parameter set
ignore_params	parameters to be ignored (optional)

Returns: a parsed configuration object

See also: ValidateAndConfigure(), Configure()

This method will validate the parameter set (except for the parameters it's explicitly told to ignore) and it returns an object ready to be used with Configure().

Definition at line 258 of file DetectorPropertiesProtoDUNEsp.cxx.

                                                          {
     std::set<std::string> ignorable_keys = lar::IgnorableProviderConfigKeys();
     ignorable_keys.insert(ignore_params.begin(), ignore_params.end());
 
     // parses and validates the parameter set:
     fhicl::Table<Configuration_t> config_table { p, ignorable_keys };
 
     return std::move(config_table());
 
   } // DetectorPropertiesStandard::ValidateConfiguration()

Member Data Documentation

std::vector<double> spdp::DetectorPropertiesProtoDUNEsp::fEfield

protected

kV/cm (per inter-plane volume)

Definition at line 332 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fElectronlifetime

protected

microseconds

Definition at line 333 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fElectronsToADC

protected

conversion factor for # of ionization electrons to 1 ADC count

Definition at line 336 of file DetectorPropertiesProtoDUNEsp.h.

const geo::GeometryCore* spdp::DetectorPropertiesProtoDUNEsp::fGeo

protected

Definition at line 324 of file DetectorPropertiesProtoDUNEsp.h.

bool spdp::DetectorPropertiesProtoDUNEsp::fGetHVDriftfromMetaData

protected

Definition at line 328 of file DetectorPropertiesProtoDUNEsp.h.

bool spdp::DetectorPropertiesProtoDUNEsp::fGetReadOutWindowSizefromMetaData

protected

Definition at line 329 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fHasTimeOffsetU = false

protected

whether time offset was configured for view U

Definition at line 344 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fHasTimeOffsetV = false

protected

whether time offset was configured for view V

Definition at line 345 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fHasTimeOffsetX = false

protected

whether time offset was configured for view X

Definition at line 348 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fHasTimeOffsetY = false

protected

whether time offset was configured for view Y

Definition at line 347 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fHasTimeOffsetZ = false

protected

whether time offset was configured for view Z

Definition at line 346 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fHV_cath

protected

Definition at line 331 of file DetectorPropertiesProtoDUNEsp.h.

const detinfo::LArProperties* spdp::DetectorPropertiesProtoDUNEsp::fLP

protected

Definition at line 323 of file DetectorPropertiesProtoDUNEsp.h.

unsigned int spdp::DetectorPropertiesProtoDUNEsp::fNumberTimeSamples

protected

number of clock ticks per event

Definition at line 337 of file DetectorPropertiesProtoDUNEsp.h.

unsigned int spdp::DetectorPropertiesProtoDUNEsp::fReadOutWindowSize

protected

number of clock ticks per readout window

Definition at line 338 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fSamplingRate

protected

in ns

Definition at line 335 of file DetectorPropertiesProtoDUNEsp.h.

bool spdp::DetectorPropertiesProtoDUNEsp::fSimpleBoundary

protected

Definition at line 352 of file DetectorPropertiesProtoDUNEsp.h.

SternheimerParameters_t spdp::DetectorPropertiesProtoDUNEsp::fSternheimerParameters

protected

Sternheimer parameters.

Definition at line 350 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fTemperature

protected

kelvin

Definition at line 334 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetU

protected

time offset to convert spacepoint coordinates to hit times on view U

Definition at line 339 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetV

protected

time offset to convert spacepoint coordinates to hit times on view V

Definition at line 340 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetX

protected

time offset to convert spacepoint coordinates to hit times on view X

Definition at line 343 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetY

protected

time offset to convert spacepoint coordinates to hit times on view Y

Definition at line 342 of file DetectorPropertiesProtoDUNEsp.h.

double spdp::DetectorPropertiesProtoDUNEsp::fTimeOffsetZ

protected

time offset to convert spacepoint coordinates to hit times on view Z

Definition at line 341 of file DetectorPropertiesProtoDUNEsp.h.

bool spdp::DetectorPropertiesProtoDUNEsp::fUseRunDependentTemperature

protected

Definition at line 330 of file DetectorPropertiesProtoDUNEsp.h.

The documentation for this class was generated from the following files:

duneprototypes/duneprototypes/Protodune/singlephase/DetectorServices/Providers/DetectorPropertiesProtoDUNEsp.h
duneprototypes/duneprototypes/Protodune/singlephase/DetectorServices/Providers/DetectorPropertiesProtoDUNEsp.cxx

Classes

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation