phsp_unweighting.inl File Reference

#include <iostream>
#include <assert.h>
#include <time.h>
#include <vector>
#include <array>
#include <chrono>
#include <tclap/CmdLine.h>
#include <hydra/Types.h>
#include <hydra/Vector4R.h>
#include <hydra/PhaseSpace.h>
#include <hydra/Function.h>
#include <hydra/Lambda.h>
#include <hydra/Algorithm.h>
#include <hydra/Tuple.h>
#include <hydra/host/System.h>
#include <hydra/device/System.h>
#include <hydra/Decays.h>
#include <hydra/DenseHistogram.h>
#include <hydra/Range.h>

Include dependency graph for phsp_unweighting.inl:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Macros
#define	PHSP_UNWEIGHTING_INL_

Functions
cmd	add (NArg)
	catch (TCLAP::ArgException &e)
	declarg (A, hydra::Vector4R) declarg(B
hydra::Vector4R	declarg (C, hydra::Vector4R) int main(int argv
phsp	Generate (Parent, Events)
TCLAP::ValueArg< size_t >	NArg ("n", "nevents", "Number of events to generate. Default is [ 10e6 ].", true, 10e6, "unsigned long")
hydra::Vector4R	Parent (P_mass, 0.0, 0.0, 0.0)
cmd	parse (argv, argc)

Variables
double	A_mass = 0.13957061
hydra::Vector4R char **	argc
double	B_mass = 0.13957061
double	C_mass = 0.13957061
auto	dalitz_calculator
std::chrono::duration< double, std::milli >	elapsed = end - start
auto	end = std::chrono::high_resolution_clock::now()
double	masses [3] {A_mass, B_mass, C_mass }
	nentries = NArg.getValue()
double	P_mass = 0.493677
hydra::PhaseSpace< 3 >	phsp {P_mass, masses}
auto	start = std::chrono::high_resolution_clock::now()
	try

Macro Definition Documentation

PHSP_UNWEIGHTING_INL_

#define PHSP_UNWEIGHTING_INL_

Function Documentation

add()

cmd add

(

NArg

)

catch()

catch

(

TCLAP::ArgException &

e

)

declarg() [1/2]

declarg	(	A	,
		hydra::Vector4R
	)

Examples:

phsp_unweighting.inl.

declarg() [2/2]

hydra::Vector4R declarg	(	C	,
		hydra::Vector4R
	)

Generate()

phsp Generate	(	Parent	,
		Events
	)

NArg()

TCLAP::ValueArg<size_t> NArg	(	"n"	,
		"nevents"	,
		"Number of events to generate. Default is ."	[10e6],
		true	,
		10e6	,
		"unsigned long"
	)

Examples:

phsp_unweighting.inl.

Parent()

hydra::Vector4R Parent	(	P_mass	,
		0.	0,
		0.	0,
		0.	0
	)

Examples:

phsp_averaging_functor.inl, phsp_basic.inl, phsp_reweighting.inl, phsp_unweighting.inl, and phsp_unweighting_functor.inl.

parse()

cmd parse	(	argv	,
		argc
	)

Variable Documentation

A_mass

double A_mass = 0.13957061

Examples:

phsp_averaging_functor.inl, phsp_basic.inl, phsp_reweighting.inl, phsp_unweighting.inl, and phsp_unweighting_functor.inl.

argc

hydra::Vector4R char** argc

Initial value:

{


    size_t  nentries   = 0

Examples:

phsp_unweighting.inl.

B_mass

double B_mass = 0.13957061

Examples:

phsp_averaging_functor.inl, phsp_basic.inl, phsp_reweighting.inl, phsp_unweighting.inl, and phsp_unweighting_functor.inl.

C_mass

double C_mass = 0.13957061

Examples:

phsp_averaging_functor.inl, phsp_basic.inl, phsp_reweighting.inl, phsp_unweighting.inl, and phsp_unweighting_functor.inl.

dalitz_calculator

auto dalitz_calculator

Initial value:

= hydra::wrap_lambda(
            [] __hydra_dual__ (A a, B b, C c) {

        return hydra::make_tuple( (a + b).mass2(), (b + c).mass2());
    })

Examples:

phsp_unweighting.inl.

elapsed

std::chrono::duration<double, std::milli> elapsed = end - start

Examples:

phsp_unweighting.inl.

end

auto end = std::chrono::high_resolution_clock::now()

masses

double masses[3] {A_mass, B_mass, C_mass }

Examples:

dalitz_plot.C, dalitz_plot.inl, phsp_averaging_functor.inl, phsp_basic.inl, phsp_reweighting.inl, phsp_unweighting.inl, phsp_unweighting_functor.inl, and range_semantics.inl.

nentries

nentries = NArg.getValue()

Examples:

phsp_unweighting.inl.

P_mass

double P_mass = 0.493677

Examples:

phsp_averaging_functor.inl, phsp_basic.inl, phsp_reweighting.inl, phsp_unweighting.inl, and phsp_unweighting_functor.inl.

phsp

hydra::PhaseSpace<3> phsp {P_mass, masses}

Examples:

phsp_unweighting.inl.

start

auto start = std::chrono::high_resolution_clock::now()

Examples:

phsp_unweighting.inl.

try

try

Initial value:

{

        TCLAP::CmdLine cmd("Command line arguments for PHSP B0 -> J/psi K pi", '=')