util::DetectorPropertiesServiceArgoNeuT Class Reference

#include "DetectorPropertiesServiceArgoNeuT.h"

Inheritance diagram for util::DetectorPropertiesServiceArgoNeuT:

Public Types
using	provider_type = DetectorPropertiesServiceArgoNeuT
	type of service provider, that is this very same More...

Public Member Functions
	DetectorPropertiesServiceArgoNeuT (fhicl::ParameterSet const &pset, art::ActivityRegistry &reg)
virtual void	reconfigure (fhicl::ParameterSet const &pset) override
virtual const detinfo::DetectorProperties *	provider () const override
	Returns our service provider, that is this very same class. More...
virtual double	Efield (unsigned int planegap=0) const override
	< kV/cm More...
virtual double	DriftVelocity (double efield=0., double temperature=0.) const override
virtual double	BirksCorrection (double dQdX) const override
	dQ/dX in electrons/cm, returns dE/dX in MeV/cm. More...
virtual double	ModBoxCorrection (double dQdX) const override
virtual double	ElectronLifetime () const override
	Returns the attenuation constant for ionization electrons. More...
virtual double	Temperature () const override
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	Eloss (double mom, double mass, double tcut) const override
	Restricted mean energy loss ( ) More...
virtual double	ElossVar (double mom, double mass) const override
	Energy loss fluctuation ( ) More...
virtual double	SamplingRate () const override
	Returns the period of the TPC readout electronics clock. More...
virtual double	ElectronsToADC () const override
virtual unsigned int	NumberTimeSamples () const override
virtual unsigned int	ReadOutWindowSize () const override
virtual int	TriggerOffset () const override
virtual double	TimeOffsetU () const override
virtual double	TimeOffsetV () const override
virtual double	TimeOffsetZ () const override
virtual double	TimeOffsetY () const override
virtual double	ConvertXToTicks (double X, int p, int t, int c) const override
virtual double	ConvertXToTicks (double X, geo::PlaneID const &planeid) const override
virtual double	ConvertTicksToX (double ticks, int p, int t, int c) const override
virtual double	ConvertTicksToX (double ticks, geo::PlaneID const &planeid) const override
virtual double	GetXTicksOffset (int p, int t, int c) const override
virtual double	GetXTicksOffset (geo::PlaneID const &planeid) const override
virtual double	GetXTicksCoefficient (int t, int c) const override
virtual double	GetXTicksCoefficient (geo::TPCID const &tpcid) const override
virtual double	GetXTicksCoefficient () const override
virtual double	ConvertTDCToTicks (double tdc) const override
virtual double	ConvertTicksToTDC (double ticks) const override
virtual bool	SimpleBoundary () const override

Private Member Functions
void	checkDBstatus () const
void	preProcessEvent (const art::Event &evt)
void	postOpenFile (std::string const &filename)
void	CalculateXTicksParams () const

Static Private Member Functions
static bool	isDetectorPropertiesServiceArgoNeuT (const fhicl::ParameterSet &ps)

Private Attributes
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 offsets to convert spacepoint More...
double	fTimeOffsetV
	coordinates to hit times on each More...
double	fTimeOffsetZ
	view More...
bool	fInheritNumberTimeSamples
	Flag saying whether to inherit NumberTimeSamples. More...
bool	fXTicksParamsLoaded
	calculations More...
double	fXTicksCoefficient
	Parameters for x<–>ticks. More...
std::vector< std::vector< std::vector< double > > >	fXTicksOffsets
std::vector< std::vector< double > >	fDriftDirection
fhicl::ParameterSet	fPS
	Original parameter set. More...
bool	fAlreadyReadFromDB
	tests whether the values have alread been picked up from the Database More...
detinfo::ElecClock	fTPCClock
	TPC electronics clock. More...
util::LArPropertiesServiceArgoNeuT const *	fLP
	Pointer to the specific LArPropertiesServiceArgoNeuT service (provider) More...
bool	fSimpleBoundary

Detailed Description

Definition at line 33 of file DetectorPropertiesServiceArgoNeuT.h.

Member Typedef Documentation

using util::DetectorPropertiesServiceArgoNeuT::provider_type = DetectorPropertiesServiceArgoNeuT

type of service provider, that is this very same

Definition at line 45 of file DetectorPropertiesServiceArgoNeuT.h.

Constructor & Destructor Documentation

util::DetectorPropertiesServiceArgoNeuT::DetectorPropertiesServiceArgoNeuT	(	fhicl::ParameterSet const &	pset,
		art::ActivityRegistry &	reg
	)

Definition at line 27 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References art::errors::Configuration, fLP, postOpenFile(), preProcessEvent(), reconfigure(), art::ActivityRegistry::sPostOpenFile, and art::ActivityRegistry::sPreProcessEvent.

    : fNumberTimeSamples(pset.get< unsigned int >("NumberTimeSamples"))
  {
    fLP = dynamic_cast<util::LArPropertiesServiceArgoNeuT const*>
      (lar::providerFrom<detinfo::LArPropertiesService>());
    if (!fLP) {
      // this legacy service works only coupled to the corresponding
      // LArProperties legacy service:
      throw art::Exception(art::errors::Configuration)
        << "DetectorPropertiesServiceArgoNeuT service requires"
        " LArPropertiesServiceArgoNeuT";
    }
    
    this->reconfigure(pset);

    // Register for callbacks.

    reg.sPostOpenFile.watch    (this, &DetectorPropertiesServiceArgoNeuT::postOpenFile);
    reg.sPreProcessEvent.watch (this, &DetectorPropertiesServiceArgoNeuT::preProcessEvent);
  }

Member Function Documentation

virtual double util::DetectorPropertiesServiceArgoNeuT::BirksCorrection ( double  dQdX ) const

inlineoverridevirtual

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

Implements detinfo::DetectorProperties.

Definition at line 62 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::BirksCorrection(), and fLP.

        { return fLP->BirksCorrection(dQdX); }

void util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams ( ) const

private

< special case for ArgoNeuT

Definition at line 170 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References geo::GeometryCore::Cryostat(), dir, geo::TPCGeo::DriftDirection(), DriftVelocity(), Efield(), fDriftDirection, fLP, fTimeOffsetU, fTimeOffsetV, fTimeOffsetZ, fXTicksCoefficient, fXTicksOffsets, fXTicksParamsLoaded, geo::kNegX, geo::kU, geo::kV, geo::kZ, geo::GeometryCore::Ncryostats(), geo::TPCGeo::Nplanes(), geo::CryostatGeo::NTPC(), geo::TPCGeo::Plane(), geo::TPCGeo::PlaneLocation(), geo::TPCGeo::PlanePitch(), SamplingRate(), util::LArPropertiesServiceArgoNeuT::Temperature(), geo::CryostatGeo::TPC(), TriggerOffset(), and geo::PlaneGeo::View().

Referenced by ConvertTicksToX(), ConvertXToTicks(), GetXTicksCoefficient(), GetXTicksOffset(), and SimpleBoundary().

  {
    art::ServiceHandle<geo::Geometry>            geo;

    double samplingRate   = SamplingRate();
    double efield         = Efield();
    double temperature    = fLP->Temperature();
    double driftVelocity  = DriftVelocity(efield, temperature);
    
    fXTicksCoefficient    = 0.001 * driftVelocity * samplingRate;

    double triggerOffset  = TriggerOffset();

    fXTicksOffsets.clear();
    fXTicksOffsets.resize(geo->Ncryostats());

    fDriftDirection.clear();
    fDriftDirection.resize(geo->Ncryostats());

    for(size_t cstat = 0; cstat < geo->Ncryostats(); ++cstat){
      fXTicksOffsets[cstat].resize(geo->Cryostat(cstat).NTPC());
      fDriftDirection[cstat].resize(geo->Cryostat(cstat).NTPC());

      for(size_t tpc = 0; tpc < geo->Cryostat(cstat).NTPC(); ++tpc) {
        const geo::TPCGeo& tpcgeom = geo->Cryostat(cstat).TPC(tpc);

        const double dir((tpcgeom.DriftDirection() == geo::kNegX) ? +1.0 :-1.0);
        fDriftDirection[cstat][tpc] = dir;

        int nplane = tpcgeom.Nplanes();
        fXTicksOffsets[cstat][tpc].resize(nplane, 0.);
        for(int plane = 0; plane < nplane; ++plane) {
          const geo::PlaneGeo& pgeom = tpcgeom.Plane(plane);
          
          
          // Get field in gap between planes
          double efieldgap[3];
          double driftVelocitygap[3];
          double fXTicksCoefficientgap[3];
          for (int igap = 0; igap<3; ++igap){
            efieldgap[igap] = Efield(igap);
            driftVelocitygap[igap] = DriftVelocity(efieldgap[igap], temperature);
            fXTicksCoefficientgap[igap] = 0.001 * driftVelocitygap[igap] * samplingRate;
          }
          
          // Calculate geometric time offset.
          // only works if xyz[0]<=0
          const double* xyz = tpcgeom.PlaneLocation(0);
          
          fXTicksOffsets[cstat][tpc][plane] = -xyz[0]/(dir * fXTicksCoefficient) + triggerOffset;

          if (nplane==3){
            /*
         |    ---------- plane = 2 (collection)
         |                      Coeff[2]
         |    ---------- plane = 1 (2nd induction)
         |                      Coeff[1]
         |    ---------- plane = 0 (1st induction) x = xyz[0]
         |                      Coeff[0]
         |    ---------- x = 0
         V     For plane = 0, t offset is -xyz[0]/Coeff[0]
         x   */
            for (int ip = 0; ip < plane; ++ip){
              fXTicksOffsets[cstat][tpc][plane] += tpcgeom.PlanePitch(ip,ip+1)/fXTicksCoefficientgap[ip+1];
            }
          }       
          else if (nplane==2){ 
            /*
         |    ---------- plane = 1 (collection)
         |                      Coeff[2]
         |    ---------- plane = 0 (2nd induction) x = xyz[0]
         |    ---------- x = 0, Coeff[1]
         V    ---------- first induction plane
         x                      Coeff[0]
For plane = 0, t offset is pitch/Coeff[1] - (pitch+xyz[0])/Coeff[0]
                         = -xyz[0]/Coeff[0] - pitch*(1/Coeff[0]-1/Coeff[1])
            */
            for (int ip = 0; ip < plane; ++ip){
              fXTicksOffsets[cstat][tpc][plane] += tpcgeom.PlanePitch(ip,ip+1)/fXTicksCoefficientgap[ip+2];
            }
            fXTicksOffsets[cstat][tpc][plane] -= tpcgeom.PlanePitch()*(1/fXTicksCoefficient-1/fXTicksCoefficientgap[1]);
          }
          
          // Add view dependent offset
          geo::View_t view = pgeom.View();
          if(view == geo::kU)
            fXTicksOffsets[cstat][tpc][plane] += fTimeOffsetU;
          else if(view == geo::kV)
            fXTicksOffsets[cstat][tpc][plane] += fTimeOffsetV;
          else if(view == geo::kZ)
            fXTicksOffsets[cstat][tpc][plane] += fTimeOffsetZ;
          else
            throw cet::exception("DetectorPropertiesServiceArgoNeuT") << "Bad view = "
                                                       << view << "\n" ;
        }       
      }
    }

    fXTicksParamsLoaded=true;
  }

void util::DetectorPropertiesServiceArgoNeuT::checkDBstatus ( ) const

private

Definition at line 103 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References fAlreadyReadFromDB.

Referenced by SimpleBoundary().

{
  bool fToughErrorTreatment= art::ServiceHandle<util::DatabaseUtil>()->ToughErrorTreatment();
  bool fShouldConnect =  art::ServiceHandle<util::DatabaseUtil>()->ShouldConnect();
  //Have not read from DB, should read and requested tough treatment
    if(!fAlreadyReadFromDB && fToughErrorTreatment && fShouldConnect )
      throw cet::exception("DetectorPropertiesServiceArgoNeuT") << " Extracting values from DetectorPropertiesServiceArgoNeuT before they "
              << " have been read in from database. \n "
              << "Set ToughErrorTreatment or ShouldConnect "
              << " to false in databaseutil.fcl if you want "
              << " to avoid this. \n";
   //Have not read from DB, should read and requested soft treatment
    else if(!fAlreadyReadFromDB && !fToughErrorTreatment && fShouldConnect )
      mf::LogWarning("DetectorPropertiesServiceArgoNeuT") <<  "!!! Extracting values from DetectorPropertiesServiceArgoNeuT before they "
              << " have been read in from the database. \n "
              << " You may not be using the correct values of "
              << " T0!"
              << " You should not be initializing"
              << " Database originating values in BeginJob()s or constructors."
              << " You have been warned !!! \n ";

    //In other cases, either already read from DB, or should not connect so it doesn't matter
}

double util::DetectorPropertiesServiceArgoNeuT::ConvertTDCToTicks ( double  tdc ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 50 of file DetectorPropertiesServiceArgoNeuT_service.cc.

Referenced by GetXTicksCoefficient().

  {
    return lar::providerFrom<detinfo::DetectorClocksService>()->TPCTDC2Tick(tdc);
  }

double util::DetectorPropertiesServiceArgoNeuT::ConvertTicksToTDC ( double  ticks ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 56 of file DetectorPropertiesServiceArgoNeuT_service.cc.

Referenced by GetXTicksCoefficient().

  {
    return lar::providerFrom<detinfo::DetectorClocksService>()->TPCTick2TDC(ticks);
  }

double util::DetectorPropertiesServiceArgoNeuT::ConvertTicksToX ( double  ticks, int  p, int  t, int  c  ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 160 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References CalculateXTicksParams(), fDriftDirection, fXTicksCoefficient, fXTicksOffsets, and fXTicksParamsLoaded.

Referenced by ConvertTicksToX(), and ConvertXToTicks().

  {
    if(!fXTicksParamsLoaded) CalculateXTicksParams();
    return (ticks - fXTicksOffsets.at(c).at(t).at(p)) * fXTicksCoefficient * fDriftDirection.at(c).at(t);  
  }

virtual double util::DetectorPropertiesServiceArgoNeuT::ConvertTicksToX ( double  ticks, geo::PlaneID const &  planeid  ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 98 of file DetectorPropertiesServiceArgoNeuT.h.

References ConvertTicksToX(), geo::CryostatID::Cryostat, GetXTicksOffset(), geo::PlaneID::Plane, and geo::TPCID::TPC.

        { return ConvertTicksToX(ticks, planeid.Plane, planeid.TPC, planeid.Cryostat); }

double util::DetectorPropertiesServiceArgoNeuT::ConvertXToTicks ( double  X, int  p, int  t, int  c  ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 148 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References CalculateXTicksParams(), fDriftDirection, fXTicksCoefficient, fXTicksOffsets, and fXTicksParamsLoaded.

Referenced by ConvertXToTicks(), and TimeOffsetY().

  {
    if(!fXTicksParamsLoaded) CalculateXTicksParams();
    return (X / (fXTicksCoefficient * fDriftDirection.at(c).at(t)) +  fXTicksOffsets.at(c).at(t).at(p) );
  }

virtual double util::DetectorPropertiesServiceArgoNeuT::ConvertXToTicks ( double  X, geo::PlaneID const &  planeid  ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 95 of file DetectorPropertiesServiceArgoNeuT.h.

References ConvertTicksToX(), ConvertXToTicks(), geo::CryostatID::Cryostat, geo::PlaneID::Plane, and geo::TPCID::TPC.

        { return ConvertXToTicks(X, planeid.Plane, planeid.TPC, planeid.Cryostat); }

virtual double util::DetectorPropertiesServiceArgoNeuT::Density ( double  temperature ) const

inlineoverridevirtual

Returns argon density at a given temperature.

Parameters

temperature

the temperature in kelvin

Returns

argon density in g/cm^3

Implements detinfo::DetectorProperties.

Definition at line 74 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::Density(), and fLP.

        { return fLP->Density(temperature); }

virtual double util::DetectorPropertiesServiceArgoNeuT::Density ( ) const

inlineoverridevirtual

Returns argon density at the temperature from Temperature()

Reimplemented from detinfo::DetectorProperties.

Definition at line 76 of file DetectorPropertiesServiceArgoNeuT.h.

References Density(), and Temperature().

Referenced by Density().

{ return Density(Temperature()); }

virtual double util::DetectorPropertiesServiceArgoNeuT::DriftVelocity ( double  efield = 0., double  temperature = 0.  ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 58 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::DriftVelocity(), and fLP.

Referenced by CalculateXTicksParams().

        { return fLP->DriftVelocity(efield, temperature); }

virtual double util::DetectorPropertiesServiceArgoNeuT::Efield ( unsigned int  planegap = 0 ) const

inlineoverridevirtual

< kV/cm

Implements detinfo::DetectorProperties.

Definition at line 55 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::Efield(), and fLP.

Referenced by CalculateXTicksParams().

        { return fLP->Efield(planegap); }

virtual double util::DetectorPropertiesServiceArgoNeuT::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: (t is measured in microseconds).

This is a uniform, constant value for the detector.

Implements detinfo::DetectorProperties.

Definition at line 68 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::ElectronLifetime(), and fLP.

        { return fLP->ElectronLifetime(); }

virtual double util::DetectorPropertiesServiceArgoNeuT::ElectronsToADC ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 85 of file DetectorPropertiesServiceArgoNeuT.h.

References fElectronsToADC.

{ return fElectronsToADC; }

virtual double util::DetectorPropertiesServiceArgoNeuT::Eloss ( double  mom, double  mass, double  tcut  ) const

inlineoverridevirtual

Restricted mean energy loss ( )

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.

Implements detinfo::DetectorProperties.

Definition at line 78 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::Eloss(), and fLP.

        { return fLP->Eloss(mom, mass, tcut); }

virtual double util::DetectorPropertiesServiceArgoNeuT::ElossVar ( double  mom, double  mass  ) const

inlineoverridevirtual

Energy loss fluctuation ( )

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

Implements detinfo::DetectorProperties.

Definition at line 81 of file DetectorPropertiesServiceArgoNeuT.h.

References util::LArPropertiesServiceArgoNeuT::ElossVar(), and fLP.

        { return fLP->ElossVar(mom, mass); }

double util::DetectorPropertiesServiceArgoNeuT::GetXTicksCoefficient ( int  t, int  c  ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 274 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References CalculateXTicksParams(), fDriftDirection, fXTicksCoefficient, and fXTicksParamsLoaded.

  {
    if(!fXTicksParamsLoaded) CalculateXTicksParams();
    return fXTicksCoefficient * fDriftDirection.at(c).at(t);
  }

virtual double util::DetectorPropertiesServiceArgoNeuT::GetXTicksCoefficient ( geo::TPCID const &  tpcid ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 105 of file DetectorPropertiesServiceArgoNeuT.h.

References ConvertTDCToTicks(), ConvertTicksToTDC(), geo::CryostatID::Cryostat, GetXTicksCoefficient(), and geo::TPCID::TPC.

        { return GetXTicksCoefficient(tpcid.TPC, tpcid.Cryostat); }

double util::DetectorPropertiesServiceArgoNeuT::GetXTicksCoefficient ( ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 283 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References CalculateXTicksParams(), fXTicksCoefficient, and fXTicksParamsLoaded.

Referenced by GetXTicksCoefficient(), and GetXTicksOffset().

  {
    if(!fXTicksParamsLoaded) CalculateXTicksParams();
    return fXTicksCoefficient;
  }

double util::DetectorPropertiesServiceArgoNeuT::GetXTicksOffset ( int  p, int  t, int  c  ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 292 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References CalculateXTicksParams(), fXTicksOffsets, and fXTicksParamsLoaded.

Referenced by ConvertTicksToX(), and GetXTicksOffset().

  {
    if(!fXTicksParamsLoaded) CalculateXTicksParams();
    return fXTicksOffsets.at(c).at(t).at(p);    
  }

virtual double util::DetectorPropertiesServiceArgoNeuT::GetXTicksOffset ( geo::PlaneID const &  planeid ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 102 of file DetectorPropertiesServiceArgoNeuT.h.

References geo::CryostatID::Cryostat, GetXTicksCoefficient(), GetXTicksOffset(), geo::PlaneID::Plane, and geo::TPCID::TPC.

        { return GetXTicksOffset(planeid.Plane, planeid.TPC, planeid.Cryostat); }

bool util::DetectorPropertiesServiceArgoNeuT::isDetectorPropertiesServiceArgoNeuT ( const fhicl::ParameterSet &  ps )

staticprivate

Definition at line 404 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References d, DEFINE_ART_SERVICE_INTERFACE_IMPL, fhicl::ParameterSet::get_if_present(), and s.

Referenced by postOpenFile(), and SimpleBoundary().

  {
    // This method uses heuristics to determine whether the parameter
    // set passed as argument is a DetectorPropertiesServiceArgoNeuT configuration
    // parameter set.

    std::string s;
    double d;
    int i;
    unsigned int u;
    
    bool result = !ps.get_if_present("module_label", s);
    result = result && ps.get_if_present("TriggerOffset", i);
    result = result && ps.get_if_present("SamplingRate", d);
    result = result && ps.get_if_present("NumberTimeSamples", u);
    result = result && ps.get_if_present("ReadOutWindowSize", u);

    return result;
  }

virtual double util::DetectorPropertiesServiceArgoNeuT::ModBoxCorrection ( double  dQdX ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 65 of file DetectorPropertiesServiceArgoNeuT.h.

References fLP, and util::LArPropertiesServiceArgoNeuT::ModBoxCorrection().

        { return fLP->ModBoxCorrection(dQdX); }

virtual unsigned int util::DetectorPropertiesServiceArgoNeuT::NumberTimeSamples ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 86 of file DetectorPropertiesServiceArgoNeuT.h.

References fNumberTimeSamples.

{ return fNumberTimeSamples; }

void util::DetectorPropertiesServiceArgoNeuT::postOpenFile ( std::string const &  filename )

private

Definition at line 301 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References file, fInheritNumberTimeSamples, fNumberTimeSamples, fPS, fhicl::ParameterSet::get(), isDetectorPropertiesServiceArgoNeuT(), and fhicl::make_ParameterSet().

Referenced by DetectorPropertiesServiceArgoNeuT(), and SimpleBoundary().

  {
    // Use this method to figure out whether to inherit configuration
    // parameters from previous jobs.
    //
    // There is no way currently to correlate parameter sets saved in
    // sqlite RootFileDB with process history (from MetaData tree).
    // Therefore, we use the approach of scanning every historical
    // parameter set in RootFileDB, and finding all parameter sets
    // that appear to be DetectorPropertiesServiceArgoNeuT configurations.  If all
    // historical parameter sets are in agreement about the value of
    // an inherited parameter, then we accept the historical value,
    // print a message, and override the configuration parameter.  In
    // cases where the historical configurations are not in agreement
    // about the value of an inherited parameter, we ignore any
    // historical parameter values that are the same as the current
    // configured value of the parameter (that is, we resolve the
    // conflict in favor of parameters values that are different than
    // the current configuration).  If two or more historical values
    // differ from the current configuration, throw an exception.
    // Note that it is possible to give precendence to the current
    // configuration by disabling inheritance for that configuration
    // parameter.

    // Don't do anything if no parameters are supposed to be inherited.

    if(!fInheritNumberTimeSamples)
      return;

    // The only way to access art service metadata from the input file
    // is to open it as a separate TFile object.  Do that now.

    if(filename.size() != 0) {

      TFile* file = TFile::Open(filename.c_str(), "READ");
      if(file != 0 && !file->IsZombie() && file->IsOpen()) {

        // Open the sqlite datatabase.

        art::SQLite3Wrapper sqliteDB(file, "RootFileDB");

        // Loop over all stored ParameterSets.

        unsigned int iNumberTimeSamples = 0;  // Combined value of NumberTimeSamples.
        unsigned int nNumberTimeSamples = 0;  // Number of NumberTimeSamples parameters seen.

        sqlite3_stmt * stmt = 0;
        sqlite3_prepare_v2(sqliteDB, "SELECT PSetBlob from ParameterSets;", -1, &stmt, NULL);
        while (sqlite3_step(stmt) == SQLITE_ROW) {
          fhicl::ParameterSet ps;
          fhicl::make_ParameterSet(reinterpret_cast<char const *>(sqlite3_column_text(stmt, 0)), ps);
          // Is this a DetectorPropertiesServiceArgoNeuT parameter set?

          bool psok = isDetectorPropertiesServiceArgoNeuT(ps);
          if(psok) {

            // Check NumberTimeSamples

            if(fInheritNumberTimeSamples) {
              unsigned int newNumberTimeSamples = ps.get<unsigned int>("NumberTimeSamples");

              // Ignore parameter values that match the current configuration.

              if(newNumberTimeSamples != fPS.get<unsigned int>("NumberTimeSamples")) {
                if(nNumberTimeSamples == 0)
                  iNumberTimeSamples = newNumberTimeSamples;
                else if(newNumberTimeSamples != iNumberTimeSamples) {
                  throw cet::exception("DetectorPropertiesServiceArgoNeuT")
                    << "Historical values of NumberTimeSamples do not agree: "
                    << iNumberTimeSamples << " " << newNumberTimeSamples << "\n" ;
                }
                ++nNumberTimeSamples;
              }
            }
          }
        }

        // Done looping over parameter sets.
        // Now decide which parameters we will actually override.

        if(fInheritNumberTimeSamples && 
           nNumberTimeSamples != 0 && 
           iNumberTimeSamples != fNumberTimeSamples) {
          mf::LogInfo("DetectorPropertiesServiceArgoNeuT")
            << "Overriding configuration parameter NumberTimeSamples using historical value.\n"
            << "  Configured value:        " << fNumberTimeSamples << "\n"
            << "  Historical (used) value: " << iNumberTimeSamples << "\n";
          fNumberTimeSamples = iNumberTimeSamples;
        }
      }

      // Close file.
      if(file != 0) {
        if(file->IsOpen())
          file->Close();
        delete file;
      }
    }
  }

void util::DetectorPropertiesServiceArgoNeuT::preProcessEvent ( const art::Event &  evt )

private

Definition at line 96 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References fTPCClock.

Referenced by DetectorPropertiesServiceArgoNeuT(), and SimpleBoundary().

  {
    // Make sure TPC Clock is updated with DetectorClocksService (though in principle it shouldn't change)
    fTPCClock = lar::providerFrom<detinfo::DetectorClocksService>()->TPCClock();
  }

virtual const detinfo::DetectorProperties* util::DetectorPropertiesServiceArgoNeuT::provider ( ) const

inlineoverridevirtual

Returns our service provider, that is this very same class.

Implements detinfo::DetectorPropertiesService.

Definition at line 50 of file DetectorPropertiesServiceArgoNeuT.h.

        { return this; }

virtual unsigned int util::DetectorPropertiesServiceArgoNeuT::ReadOutWindowSize ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 87 of file DetectorPropertiesServiceArgoNeuT.h.

References fReadOutWindowSize, and TriggerOffset().

{ return fReadOutWindowSize; }

void util::DetectorPropertiesServiceArgoNeuT::reconfigure ( fhicl::ParameterSet const &  pset )

overridevirtual

Implements detinfo::DetectorPropertiesService.

Definition at line 62 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References d, fElectronsToADC, fInheritNumberTimeSamples, fNumberTimeSamples, fPS, fReadOutWindowSize, fSimpleBoundary, fTimeOffsetU, fTimeOffsetV, fTimeOffsetZ, fTPCClock, fXTicksParamsLoaded, fhicl::ParameterSet::get(), and fhicl::ParameterSet::get_if_present().

Referenced by DetectorPropertiesServiceArgoNeuT().

  {
    //fSamplingRate             = p.get< double        >("SamplingRate"     );
    double d;
    int    i;
    bool   b;
    if(p.get_if_present<double>("SamplingRate",d))
      throw cet::exception(__FUNCTION__) << "SamplingRate is a deprecated fcl parameter for DetectorPropertiesServiceArgoNeuT!";
    if(p.get_if_present<int>("TriggerOffset",i))
      throw cet::exception(__FUNCTION__) << "TriggerOffset is a deprecated fcl parameter for DetectorPropertiesServiceArgoNeuT!";
    if(p.get_if_present<bool>("InheritTriggerOffset",b))
      throw cet::exception(__FUNCTION__) << "InheritTriggerOffset is a deprecated fcl parameter for DetectorPropertiesServiceArgoNeuT!";

    fElectronsToADC           = p.get< double        >("ElectronsToADC"   );
    fNumberTimeSamples        = p.get< unsigned int >("NumberTimeSamples");
    fReadOutWindowSize        = p.get< unsigned int >("ReadOutWindowSize");
    fTimeOffsetU              = p.get< double        >("TimeOffsetU"      );
    fTimeOffsetV              = p.get< double        >("TimeOffsetV"      );
    fTimeOffsetZ              = p.get< double        >("TimeOffsetZ"      );
    fInheritNumberTimeSamples = p.get<bool           >("InheritNumberTimeSamples", false);
    
    fSimpleBoundary           = p.get<bool           >("SimpleBoundaryProcess", true);

    fXTicksParamsLoaded = false;

    fTPCClock = lar::providerFrom<detinfo::DetectorClocksService>()->TPCClock();

    // Save the parameter set.
    fPS = p;

    return;
  }

virtual double util::DetectorPropertiesServiceArgoNeuT::SamplingRate ( ) const

inlineoverridevirtual

Returns the period of the TPC readout electronics clock.

Returns

the period of the TPC readout electronics clock [µs]

See also

detinfo::DetectorClocks::TPCClock()

Implements detinfo::DetectorProperties.

Definition at line 84 of file DetectorPropertiesServiceArgoNeuT.h.

References fTPCClock, and detinfo::ElecClock::TickPeriod().

Referenced by CalculateXTicksParams().

{ return fTPCClock.TickPeriod() * 1.e3; }

virtual bool util::DetectorPropertiesServiceArgoNeuT::SimpleBoundary ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 114 of file DetectorPropertiesServiceArgoNeuT.h.

References CalculateXTicksParams(), checkDBstatus(), fSimpleBoundary, isDetectorPropertiesServiceArgoNeuT(), postOpenFile(), and preProcessEvent().

{ return fSimpleBoundary; }

virtual double util::DetectorPropertiesServiceArgoNeuT::Temperature ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 71 of file DetectorPropertiesServiceArgoNeuT.h.

References fLP, and util::LArPropertiesServiceArgoNeuT::Temperature().

Referenced by Density().

        { return fLP->Temperature(); }

virtual double util::DetectorPropertiesServiceArgoNeuT::TimeOffsetU ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 89 of file DetectorPropertiesServiceArgoNeuT.h.

References fTimeOffsetU.

{ return fTimeOffsetU; }

virtual double util::DetectorPropertiesServiceArgoNeuT::TimeOffsetV ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 90 of file DetectorPropertiesServiceArgoNeuT.h.

References fTimeOffsetV.

{ return fTimeOffsetV; }

virtual double util::DetectorPropertiesServiceArgoNeuT::TimeOffsetY ( ) const

inlineoverridevirtual

Reimplemented from detinfo::DetectorProperties.

Definition at line 92 of file DetectorPropertiesServiceArgoNeuT.h.

References ConvertXToTicks(), and X.

{ return 0; }

virtual double util::DetectorPropertiesServiceArgoNeuT::TimeOffsetZ ( ) const

inlineoverridevirtual

Implements detinfo::DetectorProperties.

Definition at line 91 of file DetectorPropertiesServiceArgoNeuT.h.

References fTimeOffsetZ.

{ return fTimeOffsetZ; }

int util::DetectorPropertiesServiceArgoNeuT::TriggerOffset ( ) const

overridevirtual

Implements detinfo::DetectorProperties.

Definition at line 128 of file DetectorPropertiesServiceArgoNeuT_service.cc.

References fTPCClock, and detinfo::ElecClock::Ticks().

Referenced by CalculateXTicksParams(), and ReadOutWindowSize().

{
  return fTPCClock.Ticks(lar::providerFrom<detinfo::DetectorClocksService>()->TriggerOffsetTPC() * -1.);
}

Member Data Documentation

bool util::DetectorPropertiesServiceArgoNeuT::fAlreadyReadFromDB

private

tests whether the values have alread been picked up from the Database

Definition at line 151 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by checkDBstatus().

std::vector<std::vector<double> > util::DetectorPropertiesServiceArgoNeuT::fDriftDirection

mutableprivate

Definition at line 147 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), ConvertTicksToX(), ConvertXToTicks(), and GetXTicksCoefficient().

double util::DetectorPropertiesServiceArgoNeuT::fElectronsToADC

private

conversion factor for # of ionization electrons to 1 ADC count

Definition at line 134 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by ElectronsToADC(), and reconfigure().

bool util::DetectorPropertiesServiceArgoNeuT::fInheritNumberTimeSamples

private

Flag saying whether to inherit NumberTimeSamples.

Definition at line 141 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by postOpenFile(), and reconfigure().

util::LArPropertiesServiceArgoNeuT const* util::DetectorPropertiesServiceArgoNeuT::fLP

private

Pointer to the specific LArPropertiesServiceArgoNeuT service (provider)

Definition at line 156 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by BirksCorrection(), CalculateXTicksParams(), Density(), DetectorPropertiesServiceArgoNeuT(), DriftVelocity(), Efield(), ElectronLifetime(), Eloss(), ElossVar(), ModBoxCorrection(), and Temperature().

unsigned int util::DetectorPropertiesServiceArgoNeuT::fNumberTimeSamples

private

number of clock ticks per event

Definition at line 135 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by NumberTimeSamples(), postOpenFile(), and reconfigure().

fhicl::ParameterSet util::DetectorPropertiesServiceArgoNeuT::fPS

private

Original parameter set.

Definition at line 149 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by postOpenFile(), and reconfigure().

unsigned int util::DetectorPropertiesServiceArgoNeuT::fReadOutWindowSize

private

number of clock ticks per readout window

Definition at line 136 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by ReadOutWindowSize(), and reconfigure().

bool util::DetectorPropertiesServiceArgoNeuT::fSimpleBoundary

private

Definition at line 158 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by reconfigure(), and SimpleBoundary().

double util::DetectorPropertiesServiceArgoNeuT::fTimeOffsetU

private

time offsets to convert spacepoint

Definition at line 137 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), reconfigure(), and TimeOffsetU().

double util::DetectorPropertiesServiceArgoNeuT::fTimeOffsetV

private

coordinates to hit times on each

Definition at line 138 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), reconfigure(), and TimeOffsetV().

double util::DetectorPropertiesServiceArgoNeuT::fTimeOffsetZ

private

view

Definition at line 139 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), reconfigure(), and TimeOffsetZ().

detinfo::ElecClock util::DetectorPropertiesServiceArgoNeuT::fTPCClock

private

TPC electronics clock.

Definition at line 153 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by preProcessEvent(), reconfigure(), SamplingRate(), and TriggerOffset().

double util::DetectorPropertiesServiceArgoNeuT::fXTicksCoefficient

mutableprivate

Parameters for x<–>ticks.

Definition at line 144 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), ConvertTicksToX(), ConvertXToTicks(), and GetXTicksCoefficient().

std::vector<std::vector<std::vector<double> > > util::DetectorPropertiesServiceArgoNeuT::fXTicksOffsets

mutableprivate

Definition at line 146 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), ConvertTicksToX(), ConvertXToTicks(), and GetXTicksOffset().

bool util::DetectorPropertiesServiceArgoNeuT::fXTicksParamsLoaded

mutableprivate

calculations

Definition at line 142 of file DetectorPropertiesServiceArgoNeuT.h.

Referenced by CalculateXTicksParams(), ConvertTicksToX(), ConvertXToTicks(), GetXTicksCoefficient(), GetXTicksOffset(), and reconfigure().

---

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

• lardata/v06_47_03/source/lardata/Utilities/DetectorPropertiesServiceArgoNeuT.h
• lardata/v06_47_03/source/lardata/Utilities/DetectorPropertiesServiceArgoNeuT_service.cc