GooFit::Dim5 Struct Reference

#include <EvalVar.h>

Inheritance diagram for GooFit::Dim5:

Public Member Functions
	Dim5 ()
__host__ __device__ mcbooster::GReal_t	cosHELANG (const mcbooster::Vector4R p, const mcbooster::Vector4R q, const mcbooster::Vector4R d)
__host__ __device__ mcbooster::GReal_t	phi (const mcbooster::Vector4R &p4_p, const mcbooster::Vector4R &p4_d1, const mcbooster::Vector4R &p4_d2, const mcbooster::Vector4R &p4_h1, const mcbooster::Vector4R &p4_h2)
__host__ __device__ void	operator() (const mcbooster::GInt_t n, mcbooster::Vector4R **particles, mcbooster::GReal_t *variables) override

Detailed Description

Definition at line 19 of file EvalVar.h.

Constructor & Destructor Documentation

Dim5()

GooFit::Dim5::Dim5 ( )

inline

Definition at line 20 of file EvalVar.h.

References GooFit::dim.

{ dim = 4; }

Member Function Documentation

cosHELANG()

__host__ __device__ mcbooster::GReal_t GooFit::Dim5::cosHELANG ( const mcbooster::Vector4R  p, const mcbooster::Vector4R  q, const mcbooster::Vector4R  d  )

inline

Definition at line 23 of file EvalVar.h.

Referenced by operator()().

                                                                                               {
        mcbooster::GReal_t pd  = p * d;
        mcbooster::GReal_t pq  = p * q;
        mcbooster::GReal_t qd  = q * d;
        mcbooster::GReal_t mp2 = p.mass2();
        mcbooster::GReal_t mq2 = q.mass2();
        mcbooster::GReal_t md2 = d.mass2();

        return (pd * mq2 - pq * qd) / sqrt((pq * pq - mq2 * mp2) * (qd * qd - mq2 * md2));
    }

operator()()

__host__ __device__ void GooFit::Dim5::operator() ( const mcbooster::GInt_t  n, mcbooster::Vector4R **  particles, mcbooster::GReal_t *  variables  )

inlineoverride

Definition at line 94 of file EvalVar.h.

References cosHELANG(), and phi().

                                                                                                             {
        mcbooster::Vector4R ppip  = *particles[0];
        mcbooster::Vector4R ppim  = *particles[1];
        mcbooster::Vector4R pK    = *particles[2];
        mcbooster::Vector4R ppip2 = *particles[3];

        mcbooster::Vector4R pM    = ppip + ppim + pK + ppip2;
        mcbooster::Vector4R ppipi = ppip + ppim;
        mcbooster::Vector4R pKpi  = pK + ppip2;

        variables[0] = ppipi.mass();
        variables[1] = pKpi.mass();
        variables[2] = cosHELANG(pM, ppipi, ppip);
        variables[3] = cosHELANG(pM, pKpi, pK);
        variables[4] = phi(pM, ppip, ppim, pK, ppip2);
    }

phi()

__host__ __device__ mcbooster::GReal_t GooFit::Dim5::phi ( const mcbooster::Vector4R &  p4_p, const mcbooster::Vector4R &  p4_d1, const mcbooster::Vector4R &  p4_d2, const mcbooster::Vector4R &  p4_h1, const mcbooster::Vector4R &  p4_h2  )

inline

Definition at line 34 of file EvalVar.h.

Referenced by operator()().

                                                                               {
        mcbooster::Vector4R p4_d1p, p4_h1p, p4_h2p, p4_d2p;

        mcbooster::Vector4R d1_perp, d1_prime, h1_perp;
        mcbooster::Vector4R D;

        D                      = p4_d1 + p4_d2;
        mcbooster::Vector4R D2 = p4_h1 + p4_h2;

        d1_perp = p4_d1 - (D.dot(p4_d1) / D.dot(D)) * D;
        h1_perp = p4_h1 - (D2.dot(p4_h1) / D2.dot(D2)) * D2;

        // orthogonal to both D and d1_perp

        d1_prime = D.cross(d1_perp);

        d1_perp  = d1_perp / d1_perp.d3mag();
        d1_prime = d1_prime / d1_prime.d3mag();

        mcbooster::GReal_t x, y;

        x                      = d1_perp.dot(h1_perp);
        y                      = d1_prime.dot(h1_perp);
        mcbooster::GReal_t phi = atan2(y, x);
        // printf("x:%.5g, y%.5g phi %.5g\n", x, y, phi );

        if(phi < 0.0)
            phi += 2.0 * M_PI;

        mcbooster::Vector4R d1n  = p4_d1 / p4_d1.d3mag();
        mcbooster::Vector4R d2n  = p4_d2 / p4_d2.d3mag();
        mcbooster::Vector4R h1n  = p4_h1 / p4_h1.d3mag();
        mcbooster::Vector4R h2n  = p4_h2 / p4_h2.d3mag();
        mcbooster::Vector4R h12n = (p4_h1 + p4_h2);
        h12n *= 1.0 / h12n.d3mag();

        mcbooster::Vector4R n1 = d1n.cross(d2n);
        mcbooster::Vector4R n2 = h1n.cross(h2n);
        n1 *= 1.0 / n1.d3mag();
        n2 *= 1.0 / n2.d3mag();
        mcbooster::Vector4R n3 = n1.cross(n2);

        mcbooster::GReal_t cp = (n1.dot(n2));
        mcbooster::GReal_t sp = (n3.dot(h12n));

        mcbooster::GReal_t phi2 = acos(cp);

        if(sp < 0)
            phi2 *= -1;

        // printf("cp %.5g, sp %.5g, phi2 %.5g\n",cp, sp, phi2 );

        return phi2;
    }

---

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

• /home/travis/build/GooFit/GooFit/include/goofit/PDFs/physics/EvalVar.h