RawClusterUtility.h

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

// $Id: $

#ifndef CALOBASE_RAWCLUSTERUTILITY_H
#define CALOBASE_RAWCLUSTERUTILITY_H

#include "RawCluster.h"

#include <CLHEP/Vector/ThreeVector.h>

class RawClusterUtility
{
 public:
  virtual ~RawClusterUtility() {}

  //
  static inline float
  GetPolarAngle(const RawCluster& cluster, const CLHEP::Hep3Vector vertex) { return (cluster.get_position() - vertex).getTheta(); }
  //
  static inline float
  GetAzimuthAngle(const RawCluster& cluster, const CLHEP::Hep3Vector vertex) { return (cluster.get_position() - vertex).getPhi(); }
  //
  static inline float
  GetPseudorapidity(const RawCluster& cluster, const CLHEP::Hep3Vector vertex) { return (cluster.get_position() - vertex).pseudoRapidity(); }
  //
  static inline float
  GetET(const RawCluster& cluster, const CLHEP::Hep3Vector vertex)
  {
    const CLHEP::Hep3Vector cluster_vec(cluster.get_position() - vertex);
    const double mag = cluster_vec.mag();

    return mag > 0 ? cluster_vec.perp() / cluster_vec.mag() * cluster.get_energy() : 0;
  }
  //
  static inline float
  GetETCore(const RawCluster& cluster, const CLHEP::Hep3Vector vertex)
  {
    const CLHEP::Hep3Vector cluster_vec(cluster.get_position() - vertex);

    return cluster_vec.perp() / cluster_vec.mag() * cluster.get_ecore();
  }
  //
  static inline CLHEP::Hep3Vector GetEVec(const RawCluster& cluster, const CLHEP::Hep3Vector vertex) { return (cluster.get_position() - vertex).unit() * cluster.get_energy(); }
  static inline CLHEP::Hep3Vector GetECoreVec(const RawCluster& cluster, const CLHEP::Hep3Vector vertex) { return (cluster.get_position() - vertex).unit() * cluster.get_ecore(); }

 private:
  RawClusterUtility() {}
};

#endif