BarrelSimplifiedGeometry.h

Go to the documentation of this file.

#ifndef FASTSIM_BARRELSIMPLIFIEDGEOMETRY_H
#define FASTSIM_BARRELSIMPLIFIEDGEOMETRY_H
 
#include "FastSimulation/SimplifiedGeometryPropagator/interface/SimplifiedGeometry.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <TH1F.h>
 
// Author: L. Vanelderen, S. Kurz
// Date: 29 May 2017
 
namespace fastsim {
 
 
  class BarrelSimplifiedGeometry : public SimplifiedGeometry {
  public:
 
    BarrelSimplifiedGeometry(double radius) : SimplifiedGeometry(radius) {}
 
    BarrelSimplifiedGeometry(BarrelSimplifiedGeometry &&) = default;
 
    ~BarrelSimplifiedGeometry() override{};
 
 
    const double getRadius() const { return geomProperty_; }
 
 
    const double getThickness(const math::XYZTLorentzVector &position) const override {
      return thicknessHist_->GetBinContent(thicknessHist_->GetXaxis()->FindBin(fabs(position.Z())));
    }
 
 
    const double getThickness(const math::XYZTLorentzVector &position,
                              const math::XYZTLorentzVector &momentum) const override {
      // Do projection of norm(layer) on momentum vector
      // CosTheta = (momentum dot norm) / (length(momentum) / length(norm))
      return getThickness(position) /
             (fabs(momentum.X() * position.X() + momentum.Y() * position.Y()) /
              (momentum.P() * std::sqrt(position.X() * position.X() + position.Y() * position.Y())));
    }
 
 
    const double getMagneticFieldZ(const math::XYZTLorentzVector &position) const override {
      return magneticFieldHist_->GetBinContent(magneticFieldHist_->GetXaxis()->FindBin(fabs(position.z())));
    }
 
 
    bool isForward() const override { return false; }
  };
 
}  // namespace fastsim
 
#endif