Inheritance diagram for util::LArPropTest:

Public Member Functions
	LArPropTest (fhicl::ParameterSet const &pset)

Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDAnalyzer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept

	Analyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	Analyzer (Table< Config > const &config)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	beginJob () override

void	analyze (const art::Event &evt) override

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 23 of file LArPropTest_module.cc.

Constructor & Destructor Documentation

util::LArPropTest::LArPropTest ( fhicl::ParameterSet const & pset )

explicit

Definition at line 34 of file LArPropTest_module.cc.

34 : EDAnalyzer(pset) {}

art::EDAnalyzer::EDAnalyzer

EDAnalyzer(fhicl::ParameterSet const &pset)

Definition: EDAnalyzer.h:25

Member Function Documentation

void util::LArPropTest::analyze ( const art::Event & evt )

overrideprivate

Definition at line 130 of file LArPropTest_module.cc.

131 {}

void util::LArPropTest::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 37 of file LArPropTest_module.cc.

   {
     // Make sure assert is enabled.
 
     bool assert_flag = false;
     assert((assert_flag = true, assert_flag));
     if (!assert_flag) {
       std::cerr << "Assert is disabled" << std::endl;
       abort();
     }
 
     // Get services.
 
     detinfo::LArProperties const* larprop = lar::providerFrom<detinfo::LArPropertiesService>();
     auto const detprop =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataForJob();
 
     // Test (default) accessors.
 
     std::cout << "Density = " << detprop.Density() << " g/cm^3" << std::endl;
     std::cout << "Drift velocity = " << detprop.DriftVelocity() << " cm/usec" << std::endl;
     std::cout << "Efield = " << detprop.Efield() << " kV/cm" << std::endl;
     std::cout << "Temperature = " << detprop.Temperature() << " Kelvin" << std::endl;
     std::cout << "Electron lifetime = " << detprop.ElectronLifetime() << " usec" << std::endl;
     std::cout << "Radiation Length = " << larprop->RadiationLength() << " g/cm^2" << std::endl;
     std::cout << "Radiation Length = " << larprop->RadiationLength() / detprop.Density() << " cm"
               << std::endl;
 
     // Make sure default values are acting correctly.
 
     assert(detprop.Density() == detprop.Density(detprop.Temperature()));
     assert(detprop.Density() != detprop.Density(detprop.Temperature() + 0.1));
     assert(detprop.DriftVelocity() == detprop.DriftVelocity(detprop.Efield()));
     assert(detprop.DriftVelocity() ==
            detprop.DriftVelocity(detprop.Efield(), detprop.Temperature()));
 
     // Drift velocity vs. electric field.
 
     std::cout << "\nDrift Velocity vs. Electric Field.\n"
               << "      E (kV/cm)      v (cm/us)" << std::fixed << std::endl;
     for (int i = 0; i < 3; ++i) {
       double e = 0.5;
       if (i == 1)
         e = 0.666667;
       else if (i == 2)
         e = 0.8;
       double v = detprop.DriftVelocity(e);
       std::cout << std::setprecision(3) << std::setw(15) << e << std::setprecision(4)
                 << std::setw(15) << v << std::endl;
     }
 
     // dE/dx.
 
     std::cout
       << "\nCompare "
          "http://pdg.lbl.gov/2011/AtomicNuclearProperties/MUON_ELOSS_TABLES/muonloss_289.dat\n"
       << std::endl;
 
     double mass = 0.10565839; // Muon.
     double fact[16] = {
       1.0, 1.2, 1.4, 1.7, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 7.0, 8.0, 9.0};
 
     std::ios_base::fmtflags flags = std::cout.flags();
     std::cout << "     T         p     Ionization dE/dx|_R\n"
               << "   [MeV]    [MeV/c]  ---[MeV cm^2/g]----\n"
               << std::scientific << std::setprecision(3);
 
     // Loop over kinetic energy.
 
     for (double tbase = 1.; tbase <= 1.e9; tbase *= 10.) {
       for (int i = 0; i < 16; ++i) {
 
         // T in MeV.
 
         double t = tbase * fact[i];
         if (t > 1.e9) break;
 
         // Calculate momentum in GeV/C
 
         double p = std::sqrt(1.e-6 * t * t + 2.e-3 * t * mass);
 
         // Calculate restricted and unrestricted dE/dx.
 
         double dedxr = detprop.Eloss(p, mass, 0.05) / detprop.Density(); // Restricted.
         double dedx = detprop.Eloss(p, mass, 0.) / detprop.Density();    // Unrestricted.
         std::cout << std::setw(10) << t << std::setw(10) << 1000. * p << std::setw(10) << dedx
                   << std::setw(10) << dedxr << std::endl;
       }
     }
     std::cout.flags(flags);
   }

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

lardata/test/DetectorInfoServices/LArPropTest_module.cc

Public Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation