FieldReaders: read offset from configs again (09fddf63) · Commits · XRFLab109 / Allpix Squared

src/modules/DopingProfileReader/DopingProfileReaderModule.cpp

+12 −0

Original line number	Diff line number	Diff line
		@@ -59,11 +59,23 @@ void DopingProfileReaderModule::initialize() {
		field_scale = {{scales.x(), scales.y()}};
		}

		// Get the field offset in fractions of the field size, default is 0.0x0.0, i.e. no offset
		auto offset = config_.get<ROOT::Math::XYVector>("field_offset", {0.0, 0.0});
		if(offset.x() > 1.0 \|\| offset.y() > 1.0) {
		throw InvalidValueError(
		config_, "field_offset", "shifting doping profile by more than one pixel (offset > 1.0) is not allowed");
		}
		if(offset.x() < 0.0 \|\| offset.y() < 0.0) {
		throw InvalidValueError(config_, "field_offset", "offsets for the doping profile have to be positive");
		}
		LOG(DEBUG) << "Doping profile has offset of " << offset << " fractions of the field size";

		detector_->setDopingProfileGrid(field_data.getData(),
		field_data.getDimensions(),
		field_data.getSize(),
		field_mapping,
		field_scale,
		{{offset.x(), offset.y()}},
		thickness_domain);

		} else if(field_model == DopingProfile::CONSTANT) {

+12 −0

Original line number	Diff line number	Diff line
		@@ -83,11 +83,23 @@ void ElectricFieldReaderModule::initialize() {
		field_scale = {{scales.x(), scales.y()}};
		}

		// Get the field offset in fractions of the field size, default is 0.0x0.0, i.e. no offset
		auto offset = config_.get<ROOT::Math::XYVector>("field_offset", {0.0, 0.0});
		if(offset.x() > 1.0 \|\| offset.y() > 1.0) {
		throw InvalidValueError(
		config_, "field_offset", "shifting electric field by more than one pixel (offset > 1.0) is not allowed");
		}
		if(offset.x() < 0.0 \|\| offset.y() < 0.0) {
		throw InvalidValueError(config_, "field_offset", "offsets for the electric field have to be positive");
		}
		LOG(DEBUG) << "Electric field has offset of " << offset << " fractions of the field size";

		detector_->setElectricFieldGrid(field_data.getData(),
		field_data.getDimensions(),
		field_data.getSize(),
		field_mapping,
		field_scale,
		{{offset.x(), offset.y()}},
		thickness_domain);
		} else if(field_model == ElectricField::CONSTANT) {
		LOG(TRACE) << "Adding constant electric field";

+13 −0

Original line number	Diff line number	Diff line
		@@ -63,12 +63,25 @@ void WeightingPotentialReaderModule::initialize() {
		field_scale = {{scales.x(), scales.y()}};
		}

		// Get the field offset in fractions of the field size, default is 0.0x0.0, i.e. no offset
		auto offset = config_.get<ROOT::Math::XYVector>("field_offset", {0.0, 0.0});
		if(offset.x() > 1.0 \|\| offset.y() > 1.0) {
		throw InvalidValueError(config_,
		"field_offset",
		"shifting weighting potential by more than one pixel (offset > 1.0) is not allowed");
		}
		if(offset.x() < 0.0 \|\| offset.y() < 0.0) {
		throw InvalidValueError(config_, "field_offset", "offsets for the weighting potential have to be positive");
		}
		LOG(DEBUG) << "Weighting potential has offset of " << offset << " fractions of the field size";

		// Set the field grid, provide scale factors as fraction of the pixel pitch for correct scaling:
		detector_->setWeightingPotentialGrid(field_data.getData(),
		field_data.getDimensions(),
		field_data.getSize(),
		field_mapping,
		field_scale,
		{{offset.x(), offset.y()}},
		thickness_domain);
		} else if(field_model == WeightingPotential::PAD) {
		LOG(TRACE) << "Adding weighting potential from pad in plane condenser";