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.
	JFoxWolframMoments (const std::vector< Trk >::const_iterator __begin, const std::vector< Trk >::const_iterator __end, const JCylinder3D &can=getMaximumContainmentVolume(), const int N=4)
	Constructor.
	JFoxWolframMoments (const Evt &event, const JCylinder3D &can=getMaximumContainmentVolume(), const int N=4)
	Constructor.

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

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

reference

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

Definition at line 141 of file JEventShapeVariables.hh.

Constructor & Destructor Documentation

JFoxWolframMoments() [1/3]

JSIRENE::JFoxWolframMoments::JFoxWolframMoments ( )

inline

Default constructor.

Definition at line 146 of file JEventShapeVariables.hh.

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

JFoxWolframMoments() [2/3]

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;
                }
              }
            }
          }
        }
      }
    }

JFoxWolframMoments() [3/3]

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

getContribution()

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