Add histograms and alter descriptions (264a5f90) · Commits · XRFLab109 / Allpix Squared

src/modules/GenericPropagation/GenericPropagationModule.cpp

+18 −6

Original line number	Diff line number	Diff line
		@@ -548,7 +548,11 @@ void GenericPropagationModule::initialize() {
		1);

		trapped_histo_ = CreateHistogram<TH1D>(
		"trapping_histo", "Fraction of trapped charge carriers;trapping [N / N_{total}] ;number of events", 100, 0, 1);
		"trapping_histo",
		"Fraction of trapped charge carriers at final state;trapping [N / N_{total}] ;number of events",
		100,
		0,
		1);

		recombination_time_histo_ =
		CreateHistogram<TH1D>("recombination_time_histo",
		@@ -556,8 +560,15 @@ void GenericPropagationModule::initialize() {
		static_cast<int>(Units::convert(integration_time_, "ns") * 5),
		0,
		static_cast<double>(Units::convert(integration_time_, "ns")));
		trapping_time_histo_ = CreateHistogram<TH1D>("trapping_time_histo",
		"Time until trapping of charge carriers;time [ns];charge carriers",
		trapping_time_histo_ =
		CreateHistogram<TH1D>("trapping_time_histo",
		"Absolute time until trapping of charge carriers;time [ns];charge carriers",
		static_cast<int>(Units::convert(integration_time_, "ns") * 5),
		0,
		static_cast<double>(Units::convert(integration_time_, "ns")));
		detrapping_time_histo_ =
		CreateHistogram<TH1D>("detrapping_time_histo",
		"Time from trapping until detrapping of charge carriers;time [ns];charge carriers",
		static_cast<int>(Units::convert(integration_time_, "ns") * 5),
		0,
		static_cast<double>(Units::convert(integration_time_, "ns")));
		@@ -929,6 +940,7 @@ void GenericPropagationModule::finalize() {
		trapped_histo_->Write();
		recombination_time_histo_->Write();
		trapping_time_histo_->Write();
		detrapping_time_histo_->Write();
		}

		long double average_time = static_cast<long double>(total_time_picoseconds_) / 1e3 /

src/modules/GenericPropagation/GenericPropagationModule.hpp

+1 −0

Original line number	Diff line number	Diff line
		@@ -134,6 +134,7 @@ namespace allpix {
		Histogram<TH1D> trapped_histo_;
		Histogram<TH1D> recombination_time_histo_;
		Histogram<TH1D> trapping_time_histo_;
		Histogram<TH1D> detrapping_time_histo_;
		std::mutex stats_mutex_;
		};

src/modules/TransientPropagation/TransientPropagationModule.cpp

+18 −0

Original line number	Diff line number	Diff line
		@@ -165,8 +165,26 @@ void TransientPropagationModule::initialize() {
		100,
		0,
		1);
		recombination_time_histo_ =
		CreateHistogram<TH1D>("recombination_time_histo",
		"Time until recombination of charge carriers;time [ns];charge carriers",
		static_cast<int>(Units::convert(integration_time_, "ns") * 5),
		0,
		static_cast<double>(Units::convert(integration_time_, "ns")));
		trapped_histo_ = CreateHistogram<TH1D>(
		"trapping_histo", "Fraction of trapped charge carriers;trapping [N / N_{total}] ;number of events", 100, 0, 1);
		trapping_time_histo_ =
		CreateHistogram<TH1D>("trapping_time_histo",
		"Absolute time until trapping of charge carriers;time [ns];charge carriers",
		static_cast<int>(Units::convert(integration_time_, "ns") * 5),
		0,
		static_cast<double>(Units::convert(integration_time_, "ns")));
		detrapping_time_histo_ =
		CreateHistogram<TH1D>("detrapping_time_histo",
		"Time from trapping until detrapping of charge carriers;time [ns];charge carriers",
		static_cast<int>(Units::convert(integration_time_, "ns") * 5),
		0,
		static_cast<double>(Units::convert(integration_time_, "ns")));
		}
		}

src/modules/TransientPropagation/TransientPropagationModule.hpp

+3 −0

Original line number	Diff line number	Diff line
		@@ -122,5 +122,8 @@ namespace allpix {
		Histogram<TH1D> drift_time_histo_;
		Histogram<TH1D> recombine_histo_;
		Histogram<TH1D> trapped_histo_;
		Histogram<TH1D> recombination_time_histo_;
		Histogram<TH1D> trapping_time_histo_;
		Histogram<TH1D> detrapping_time_histo_;
		};
		} // namespace allpix