JSIRENE::JFoxWolframMoments Struct Reference

Class for computing Fox-Wolfram moments. More...

#include <JEventShapeVariables.hh>

Inheritance diagram for JSIRENE::JFoxWolframMoments:

Public Types
typedef std::vector< double > ::reference	reference

Public Member Functions
	JFoxWolframMoments ()
	Default constructor. More...
	JFoxWolframMoments (const std::vector< Trk >::const_iterator __begin, const std::vector< Trk >::const_iterator __end, const JCylinder3D &can=getMaximumContainmentVolume(), const int N=4)
	Constructor. More...
	JFoxWolframMoments (const Evt &event, const JCylinder3D &can=getMaximumContainmentVolume(), const int N=4)
	Constructor. More...

Private Member Functions
double	getContribution (const Vec &p0, const Vec &p1, const size_t n)
	Get contribution to Fox-Wolfram moment from two momenta. More...

Detailed Description

Class for computing Fox-Wolfram moments.

The Fox-Wolfram moments are calculated as defined in section 15.2.3 of arXiv:hep-ph/0603175v2.

Definition at line 138 of file JEventShapeVariables.hh.

Member Typedef Documentation

typedef std::vector<double>::reference JSIRENE::JFoxWolframMoments::reference

Definition at line 141 of file JEventShapeVariables.hh.

Constructor & Destructor Documentation

JSIRENE::JFoxWolframMoments::JFoxWolframMoments ( )

inline

Default constructor.

Definition at line 146 of file JEventShapeVariables.hh.

                         :
      std::vector<double>(0)      
    {}

JSIRENE::JFoxWolframMoments::JFoxWolframMoments ( const std::vector< Trk >::const_iterator  __begin, const std::vector< Trk >::const_iterator  __end, const JCylinder3D &  can = getMaximumContainmentVolume(), const int  N = 4  )

inline

Constructor.

Parameters

__begin	beginning of track data
__end	end of track data
can	detector can
N	maximum order

Definition at line 159 of file JEventShapeVariables.hh.

                                                                       :
      std::vector<double>(N+1)
    {
      using namespace std;
      using namespace JPP;

      const double Evis = getVisibleEnergy(__begin, __end, can);

      if (Evis > 0.0) {

        for (vector<Trk>::const_iterator i = __begin; i != __end; ++i) {

          if (is_finalstate(*i)) {

            const Vec p0 = getKineticEnergy(*i) * i->dir;
        
            for (vector<Trk>::const_iterator j = __begin; j != __end; ++j) {
        
              if (is_finalstate(*j)) {

                const Vec p1 = getKineticEnergy(*j) * j->dir;

                for (int k = 0; k < (int) this->size(); ++k) {
                  (*this)[k] += getContribution(p0, p1, k) / Evis / Evis;
                }
              }
            }
          }
        }
      }
    }

JSIRENE::JFoxWolframMoments::JFoxWolframMoments ( const Evt &  event, const JCylinder3D &  can = getMaximumContainmentVolume(), const int  N = 4  )

inline

Constructor.

Parameters

event	event
can	can
N	maximum order

Definition at line 202 of file JEventShapeVariables.hh.

                                                   :
      JFoxWolframMoments(event.mc_trks.begin(),
                         event.mc_trks.end(),
                         can,
                         N)
    {}

Member Function Documentation

double JSIRENE::JFoxWolframMoments::getContribution ( const Vec &  p0, const Vec &  p1, const size_t  n  )

inlineprivate

Get contribution to Fox-Wolfram moment from two momenta.

Parameters

p0	first momentum vector
p1	second momentum vector
n	Legendre polynomial index

Returns

Fox-Wolfram moment contribution

Definition at line 221 of file JEventShapeVariables.hh.

    {
      using namespace JPP;
      
      return p0.len() * p1.len() * legendre(n, p0.dot(p1) / p0.len() / p1.len());
    }

---

The documentation for this struct was generated from the following file:

• JEventShapeVariables.hh