EIC code displayed by LXR master/​include/​Rivet/​Projections/​JetShape.hh	Source navigation Diff markup Identifier search General search
	Version: master test Revert   Change

File indexing completed on 2025-05-12 09:05:02

 // -*- C++ -*-
 #ifndef RIVET_JetShape_HH
 #define RIVET_JetShape_HH
 
 #include "Rivet/Config/RivetCommon.hh"
 #include "Rivet/Projection.hh"
 #include "Rivet/Projections/JetFinder.hh"
 #include "Rivet/Particle.hh"
 #include "Rivet/Event.hh"
 #include "Rivet/Tools/Utils.hh"
 
 namespace Rivet {
 
 
   /// @brief Calculate transverse jet profiles
   ///
   /// Calculate the differential and integral jet shapes in \f$P_{\perp}\f$ for a given
   /// set of jets. This particular jet shape projection calculates jet shapes relative
   /// to jet centroids, using only the particles associated to each jet, for the hardest
   /// \f$ n \f$ jets.
   ///
   /// The rapidity scheme (\f$ \eta \f$ or \f$ y \f$) has to be specified when
   /// invoking the constructor.
   ///
   /// The differential jet shape around a given jet axis at distance interval
   /// \f$ r \pm \delta{r}/2 \f$ is defined as
   /// \f[
   /// \rho(r) =
   ///   \frac{1}{\delta r} \frac{1}{N_\mathrm{jets}}
   ///   \sum_\mathrm{jets} \frac{P_\perp(r - \delta r/2, r+\delta r/2)}{p_\perp(0, R)}
   /// \f]
   /// with \f$ 0 \le r \le R \f$ and \f$ P_\perp(r_1, r_2) = \sum_{\in [r_1, r_2)} p_\perp \f$.
   ///
   /// The integral jet shape around a given jet axes until distance \f$ r \f$ is defined as
   /// \f[
   /// \Psi(r) =
   ///   \frac{1}{N_\mathrm{jets}}
   ///   \sum_\mathrm{jets} \frac{P_\perp(0, r)}{p_\perp(0, R)}
   /// \f]
   /// with \f$ 0 \le r \le R \f$ and \f$ P_\perp(r_1, r_2) = \sum_{\in [r_1, r_2)} p_\perp \f$.
   ///
   /// The constructor expects also the binning in radius \f$ r \f$ to be supplied.
   ///
   class JetShape : public Projection {
   public:
 
     /// @name Constructors etc.
     /// @{
 
     /// Constructor from histo range and number of bins.
     JetShape(const JetFinder& jetalg,
              double rmin, double rmax, size_t nbins,
              double ptmin=0, double ptmax=DBL_MAX,
              double absrapmin=-DBL_MAX, double absrapmax=-DBL_MAX,
              RapScheme rapscheme=RAPIDITY);
 
     /// Constructor from vector of bin edges.
     JetShape(const JetFinder& jetalg, vector<double> binedges,
              double ptmin=0, double ptmax=DBL_MAX,
              double absrapmin=-DBL_MAX, double absrapmax=-DBL_MAX,
              RapScheme rapscheme=RAPIDITY);
 
     /// Clone on the heap.
     RIVET_DEFAULT_PROJ_CLONE(JetShape);
 
     /// @}
 
     /// Import to avoid warnings about overload-hiding
     using Projection::operator =;
 
 
     /// Reset projection between events.
     void clear();
 
 
     /// Do the calculation directly on a supplied collection of Jet objects.
     void calc(const Jets& jets);
 
 
   public:
 
 
     /// Number of equidistant radius bins.
     size_t numBins() const {
       return _binedges.size() - 1;
     }
 
     /// Number of jets which passed cuts.
     size_t numJets() const {
       return _diffjetshapes.size();
     }
 
     /// \f$ r_\text{min} \f$ value.
     double rMin() const {
       return _binedges.front();
     }
 
     /// \f$ r_\text{max} \f$ value.
     double rMax() const {
       return _binedges.back();
     }
 
     /// \f$ p_\perp^\text{min} \f$ value.
     double ptMin() const {
       return _ptcuts.first;
     }
 
     /// \f$ p_\perp^\text{max} \f$ value.
     double ptMax() const {
       return _ptcuts.second;
     }
 
     /// Central \f$ r \f$ value for bin @a rbin.
     double rBinMin(size_t rbin) const {
       assert(inRange(rbin, 0u, numBins()));
       return _binedges[rbin];
     }
 
     /// Central \f$ r \f$ value for bin @a rbin.
     double rBinMax(size_t rbin) const {
       assert(inRange(rbin, 0u, numBins()));
       return _binedges[rbin+1];
     }
 
     /// Central \f$ r \f$ value for bin @a rbin.
     double rBinMid(size_t rbin) const {
       assert(inRange(rbin, 0u, numBins()));
       //cout << _binedges << '\n';
       return (_binedges[rbin] + _binedges[rbin+1])/2.0;
     }
 
     /// Return value of differential jet shape profile histo bin.
     double diffJetShape(size_t ijet, size_t rbin) const {
       assert(inRange(ijet, 0u, numJets()));
       assert(inRange(rbin, 0u, numBins()));
       return _diffjetshapes[ijet][rbin];
     }
 
     /// Return value of integrated jet shape profile histo bin.
     double intJetShape(size_t ijet, size_t rbin) const {
       assert(inRange(ijet, 0u, numJets()));
       assert(inRange(rbin, 0u, numBins()));
       double rtn  = 0;
       for (size_t i = 0; i <= rbin; ++i) {
         rtn += _diffjetshapes[ijet][i];
       }
       return rtn;
     }
 
     /// @todo Provide int and diff jet shapes with some sort of area normalisation?
 
     // /// Return value of \f$ \Psi \f$ (integrated jet shape) at given radius for a \f$ p_T \f$ bin.
     // /// @todo Remove this external indexing thing
     // double psi(size_t pTbin) const {
     //   return _PsiSlot[pTbin];
     // }
 
 
   protected:
 
     /// Apply the projection to the event.
     void project(const Event& e);
 
     /// Compare projections.
     CmpState compare(const Projection& p) const;
 
 
   protected:
 
     /// @name Jet-shape parameters
     /// @{
 
     /// Vector of radius bin edges
     vector<double> _binedges;
 
     /// Lower and upper cuts on contributing jet \f$ p_\perp \f$.
     pair<double, double> _ptcuts;
 
     /// Lower and upper cuts on contributing jet (pseudo)rapidity.
     pair<double, double> _rapcuts;
 
     /// Rapidity scheme
     RapScheme _rapscheme;
 
     /// @}
 
 
     /// @name The projected jet shapes
     /// @{
 
     /// Jet shape histo -- first index is jet number, second is r bin
     vector< vector<double> > _diffjetshapes;
 
     /// @}
 
   };
 
 
 }
 
 #endif

This page was automatically generated by the 2.3.7 LXR engine. The LXR team