MaterialInteraction.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file MaterialInteraction.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2021 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

#include "Acts/Definitions/Algebra.hpp"
#include "Acts/Detector/DetectorVolume.hpp"
#include "Acts/Geometry/GeometryIdentifier.hpp"
#include "Acts/Geometry/TrackingVolume.hpp"
#include "Acts/Material/MaterialSlab.hpp"

namespace Acts {

class Surface;

struct InteractionVolume {
  const TrackingVolume* trackingVolume = nullptr;
  const Experimental::DetectorVolume* detectorVolume = nullptr;

  InteractionVolume() = default;

  InteractionVolume(const TrackingVolume* tv) : trackingVolume(tv) {}

  InteractionVolume(const Experimental::DetectorVolume* dv)
      : detectorVolume(dv) {}

  GeometryIdentifier geometryId() const {
    if (trackingVolume != nullptr) {
      return trackingVolume->geometryId();
    } else if (detectorVolume != nullptr) {
      return detectorVolume->geometryId();
    } else {
      return GeometryIdentifier();
    }
  }

  bool empty() const {
    return trackingVolume == nullptr and detectorVolume == nullptr;
  }
};

struct MaterialInteraction {
  Vector3 position = Vector3(0., 0., 0);
  double time = 0.0;
  Vector3 direction = Vector3(0., 0., 0);
  double deltaP = 0.0;
  double sigmaPhi2 = 0.0;
  double sigmaTheta2 = 0.0;
  double sigmaQoP2 = 0.0;
  Vector3 intersection = Vector3(0., 0., 0);
  GeometryIdentifier intersectionID;
  const Surface* surface = nullptr;
  InteractionVolume volume{};
  bool updatedVolumeStep = false;
  double pathCorrection = 1.;
  MaterialSlab materialSlab;
};

struct RecordedMaterial {
  // The accumulated materialInX0
  double materialInX0 = 0.;
  double materialInL0 = 0.;
  std::vector<MaterialInteraction> materialInteractions;
};

using RecordedMaterialTrack =
    std::pair<std::pair<Acts::Vector3, Acts::Vector3>, RecordedMaterial>;

}  // namespace Acts