Commit 0c369f30 authored by Lennart Huth's avatar Lennart Huth
Browse files

volume scattering length to radiation length

parent efd7c2ad
......@@ -19,10 +19,8 @@ Clusters from the first plane in Z (named the seed plane) are related to cluster
* `timestamp_from`: Defines the detector which provides the track timestamp. This detector also needs to be set as `required_detector`. If empty, the average timestamp of all clusters on the track will be used. Empty by default.
* `track_model`: Select the track model used for reconstruction. A simple line fit ignoring scattering (`straightline`) and a General-Broken-Lines (`gbl`) are currently supported. Defaults to `straightline`.
* `momentum`: Set the beam momentum. Defaults to 5 GeV
* `volume_scattering`: Select if volume scattering will be taken into account
- defaults to false
* `volume_scattering_length`: Define the scattering length of the volume
around the telescope. Defaults to `304.2 m`
* `volume_scattering`: Select if volume scattering will be taken into account - defaults to false
* `volume_radiation_length`: Define the radiation length of the volume around the telescope. Defaults to dry air with a radiation length of`304.2 m`
### Plots produced
......
......@@ -34,7 +34,7 @@ Tracking4D::Tracking4D(Configuration config, std::vector<std::shared_ptr<Detecto
timestampFrom = m_config.get<std::string>("timestamp_from", {});
trackModel = m_config.get<std::string>("track_model", "straightline");
momentum = m_config.get<double>("momentum", Units::get<double>(5, "GeV"));
volumeScatteringLength = m_config.get<double>("volume_scattering_length", Units::get<double>(304.2, "m"));
volumeRadiationLength = m_config.get<double>("volume_radiation_length", Units::get<double>(304.2, "m"));
useVolumeScatterer = m_config.get<bool>("volume_scattering", false);
// print a warning if volumeScatterer are used as this causes fit failures
......@@ -194,7 +194,7 @@ StatusCode Tracking4D::run(std::shared_ptr<Clipboard> clipboard) {
// Add the cluster to the track
track->addCluster(cluster);
track->setTimestamp(cluster->timestamp());
track->setVolumeScatter(volumeScatteringLength);
track->setVolumeScatter(volumeRadiationLength);
track->useVolumeScatterer(useVolumeScatterer);
refTrack->addCluster(cluster);
refTrack->setTimestamp(cluster->timestamp());
......
......@@ -47,7 +47,7 @@ namespace corryvreckan {
std::map<std::string, TH1F*> pullY;
// Cuts for tracking
double momentum{};
double volumeScatteringLength{};
double volumeRadiationLength{};
double time_cut_reference_;
size_t minHitsOnTrack;
bool excludeDUT;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment