EventEtDensity.h

Go to the documentation of this file.

// -*- c++ -*-
 
/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef EVENTSHAPEEVENT_EVENTETDENSITY_H
#define EVENTSHAPEEVENT_EVENTETDENSITY_H
 
#include "fastjet/RangeDefinition.hh"
 
#include <cmath>
 
class EventEtDensity
{
public:

  EventEtDensity();
  EventEtDensity(double rho,double sigma=0.,double area=40.0*M_PI_2,
         double etaLow=-5.,double etaHigh=5.);
  EventEtDensity(double rho,double sigma,double area,
         const fastjet::RangeDefinition& range);
  // Use default copy constructor.
  ~EventEtDensity();

  enum DataIndex { RHOVALUE     = 2,
           RHOSIGMA     = 3,
           AREAVALUE    = 4,
           ETARANGELOW  = 0,
           ETARANGEHIGH = 1,
           MAXINDEX     = 5,
           INVALIDINDEX = -1 };

  double rhoValue;
  double rhoSigma;
  double areaValue;
  double etaRangeLow;
  double etaRangeHigh;
 
  /**/
  static size_t capacity();

  fastjet::RangeDefinition range() const;
  static fastjet::RangeDefinition range(const EventEtDensity& data);
  static bool copyData(const std::vector<double>& lparm,
               EventEtDensity&            data,
               size_t                     lindx=0);
  bool copyData(const std::vector<double>& lparm,size_t lindx=0);
 
  /**/
  static bool extractData(const EventEtDensity& data,
              std::vector<double>& lparm,size_t lindx=0);
  /**/
  bool extractData(std::vector<double>& lparm,size_t lindx=0) const;

  static EventEtDensity createStore(const std::vector<double>& lparm,
                    size_t lindx=0);
  static void clear(EventEtDensity& data);
  void clear();
 
 
  static void reset(EventEtDensity& data);
  void reset();
 
  bool inRange(double eta) const;
 
// private:
 
//   /*! @brief Cached @f$ \rho @f$ */
//   double m_rho;
//   /*! @brief Cached @f$ \sigma @f$ */
//   double m_sigma;
//   /*! @brief Cached integrated jet catchment area */
//   double m_area;
//   /*! @brief Cached lower @f$ \eta @f$ range limit */
//   double m_etaLow;
//   /*! @brief Cached upper @f$ \eta @f$ range limit */
//   double m_etaHigh;
};

 
inline EventEtDensity::EventEtDensity()
{
  rhoValue     = 0.;
  rhoSigma     = 0.;
  etaRangeLow  = -5.;
  etaRangeHigh = 5.;
  areaValue    = (etaRangeHigh-etaRangeLow)*2.*M_PI;
}
inline EventEtDensity::EventEtDensity(double rho,double sigma,double area,
                      double etaLow,double etaHigh)
{
  rhoValue     = rho;
  rhoSigma     = sigma;
  areaValue    = area;
  etaRangeLow  = etaLow;
  etaRangeHigh = etaHigh;
}
inline EventEtDensity::EventEtDensity(double rho,double sigma,double area,
                      const fastjet::RangeDefinition& range)
{ 
  rhoValue  = rho;
  rhoSigma  = sigma;
  areaValue = area;
  range.get_rap_limits(etaRangeLow,etaRangeHigh); 
}
 
inline EventEtDensity::~EventEtDensity()
{ }
 
inline fastjet::RangeDefinition EventEtDensity::range() const
{ return fastjet::RangeDefinition(etaRangeLow,etaRangeHigh); }
 
inline fastjet::RangeDefinition EventEtDensity::range(const EventEtDensity&
                              data)
{ return data.range(); }
 
inline bool
EventEtDensity::copyData(const std::vector<double>& lparm,EventEtDensity& data,
             size_t lindx)
{ return data.copyData(lparm,lindx); }                   
 
inline 
EventEtDensity EventEtDensity::createStore(const std::vector<double>& lparm,
                       size_t                     lindx)
{ EventEtDensity data; copyData(lparm,data,lindx); return data; }
 
inline 
bool EventEtDensity::extractData(const EventEtDensity& data,
                 std::vector<double>& lparm,size_t lindx)
{ return data.extractData(lparm,lindx); }
 
inline void EventEtDensity::clear(EventEtDensity& data)
{ data.clear(); }
 
inline size_t EventEtDensity::capacity()
{ return (size_t)MAXINDEX; }
 
inline void EventEtDensity::reset(EventEtDensity& data) 
{ data.reset(); }
 
inline bool EventEtDensity::inRange(double eta) const 
{ return eta >= etaRangeLow && eta <= etaRangeHigh; }
 
// inline EventEtDensity::EventEtDensity()
//   : m_rho(0.), m_sigma(0.), m_area(40.*asin(1.)),
//   , m_etaLow(-5.), m_etaHigh(5.)
// { }
 
// inline EventEtDensity::EventEtDensity(double rho,double sigma,double area,
//                    double etaLow,double etaHigh)
//   : m_rho(rho), m_sigma(sigma), m_area(area),
//   , m_etaLow(etaLow), m_etaHigh(etaHigh)
// { }
// inline EventEtDensity::EventEtDensity(double rho,double sigma,double area,
//                    const fastjet::RangeDefinition range)
//   : m_rho(rho), m_sigma(sigma), m_area(area)
// { range.get_rap_limits(m_etaLow,m_etaHigh); }
 
// inline EventEtDensity::EventEtDensity(const EventEtDensity& data)
//   : m_rho(data.m_rho), m_sigma(data.m_sigma), m_area(data.m_area)
//   , m_etaLow(data.m_etaLow), m_etaHigh(data.m_etaHigh)
// { }
 
// inline double EventEtDensity::rho() const
// { return m_rho; } 
 
// inline double EventEtDensity::sigma() const
// { return m_sigma; }
 
// inline double EventEtDensity::area() const
// { return m_area; }
 
// inline double EventEtDensity::etaLow() const
// { return m_etaLow; }
 
// inline double EventEtDensity::etaHigh() const
// { return m_etaHigh; }
#endif