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(){};
        

        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;
        }
    };

}

#endif