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