CaloVectors.h

Go to the documentation of this file.

#include "RecoEcal/EgammaClusterAlgos/interface/EgammaSCEnergyCorrectionAlgo.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"

math::XYZTLorentzVector DetectorVector(const reco::SuperClusterRef& sc)
{
    double pt = sc->energy()/cosh(sc->eta());
    math::XYZTLorentzVector detVec(pt*cos(sc->phi()), pt*sin(sc->phi()), pt*sinh(sc->eta()), sc->energy());
    return detVec;
}
math::XYZTLorentzVector DetectorVector(const reco::SuperCluster sc)
{
    double pt = sc.energy()/cosh(sc.eta());
    math::XYZTLorentzVector detVec(pt*cos(sc.phi()), pt*sin(sc.phi()), pt*sinh(sc.eta()), sc.energy());
    return detVec;
}
math::XYZTLorentzVector DetectorVector(const GlobalPoint& pos, const math::XYZPoint& vertex, double energy)
{
    math::XYZPoint hitPos(pos.x(), pos.y(), pos.z());
    math::XYZVector Vec = hitPos - vertex;
    double eta = Vec.Eta();
    double phi = Vec.Phi();
    double pt = energy/cosh(eta);
    math::XYZTLorentzVector detVec(pt*cos(phi), pt*sin(phi), pt*sinh(eta), energy);
    return detVec;
}
math::XYZTLorentzVector PhysicsVector(const math::XYZPoint& vertex, const reco::SuperCluster& sc)
{
    math::XYZVector Vec = sc.position() - vertex;
    double eta = Vec.Eta();
    double phi = Vec.Phi();
    double pt = sc.energy()/cosh(eta);
    math::XYZTLorentzVector probe(pt*cos(phi), pt*sin(phi), pt*sinh(eta), sc.energy());
    return probe;
}
math::XYZTLorentzVector PhysicsVectorRaw(const math::XYZPoint& vertex, const reco::SuperCluster& sc)
{
    math::XYZVector Vec = sc.position() - vertex;
    double eta = Vec.Eta();
    double phi = Vec.Phi();
    double pt = sc.rawEnergy()/cosh(eta);
    math::XYZTLorentzVector probe(pt*cos(phi), pt*sin(phi), pt*sinh(eta), sc.rawEnergy());
    return probe;
}