phsp_unweighting_functor_and_fit.inl File Reference

#include <iostream>
#include <assert.h>
#include <time.h>
#include <vector>
#include <array>
#include <chrono>
#include <tclap/CmdLine.h>
#include <hydra/device/System.h>
#include <hydra/host/System.h>
#include <hydra/Types.h>
#include <hydra/Vector4R.h>
#include <hydra/Decays.h>
#include <hydra/PhaseSpace.h>
#include <hydra/PhaseSpaceIntegrator.h>
#include <hydra/Function.h>
#include <hydra/FunctorArithmetic.h>
#include <hydra/Lambda.h>
#include <hydra/Tuple.h>
#include <hydra/Range.h>
#include <hydra/Random.h>
#include <hydra/LogLikelihoodFCN.h>
#include <hydra/Parameter.h>
#include <hydra/UserParameters.h>
#include <hydra/Pdf.h>
#include <hydra/Distance.h>
#include <hydra/functions/BreitWignerNR.h>
#include <hydra/DenseHistogram.h>
#include <hydra/SparseHistogram.h>
#include "Minuit2/FunctionMinimum.h"
#include "Minuit2/MnUserParameterState.h"
#include "Minuit2/MnPrint.h"
#include "Minuit2/MnMigrad.h"
#include "Minuit2/MnMinimize.h"
#include "Minuit2/MnMinos.h"
#include "Minuit2/MnContours.h"
#include "Minuit2/CombinedMinimizer.h"
#include "Minuit2/MnPlot.h"
#include "Minuit2/MinosError.h"
#include "Minuit2/ContoursError.h"
#include "Minuit2/VariableMetricMinimizer.h"

Include dependency graph for phsp_unweighting_functor_and_fit.inl:

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

Go to the source code of this file.

Macros
#define	PHSP_UNWEIGHTING_FUNCTOR_AND_FIT_INL_

Functions
cmd	add (NArg)
hydra::Vector4R	B0 (B0_mass, 0.0, 0.0, 0.0)
	catch (TCLAP::ArgException &e)
	declarg (P1, hydra::Vector4R) declarg(P2
hydra::Vector4R	declarg (P3, hydra::Vector4R) int main(int argv
phsp	Generate (B0, Events.begin(), Events.end())
TCLAP::ValueArg< size_t >	NArg ("n", "nevents", "Number of events to generate. Default is [ 10e6 ].", true, 10e6, "unsigned long")
cmd	parse (argv, argc)

Variables
hydra::Vector4R char **	argc
double	B0_mass = 5.27955
auto	breit_wigner
auto	dalitz_calculator
std::chrono::duration< double, std::milli >	elapsed = end - start
auto	end = std::chrono::high_resolution_clock::now()
double	Jpsi_mass = 3.0969
double	K_mass = 0.493677
hydra::Parameter	M0
auto	model
	nentries = NArg.getValue()
hydra::PhaseSpace< 3 >	phsp (B0_mass,{Jpsi_mass, K_mass, pi_mass })
double	pi_mass = 0.13957061
auto	start = std::chrono::high_resolution_clock::now()
	try
hydra::Parameter	W0

Macro Definition Documentation

PHSP_UNWEIGHTING_FUNCTOR_AND_FIT_INL_

#define PHSP_UNWEIGHTING_FUNCTOR_AND_FIT_INL_

Function Documentation

add()

cmd add

(

NArg

)

B0()

hydra::Vector4R B0	(	B0_mass	,
		0.	0,
		0.	0,
		0.	0
	)

Examples:

dalitz_plot.C, dalitz_plot.inl, and phsp_chain.inl.

catch()

catch

(

TCLAP::ArgException &

e

)

declarg() [1/2]

declarg	(	P1	,
		hydra::Vector4R
	)

declarg() [2/2]

hydra::Vector4R declarg	(	P3	,
		hydra::Vector4R
	)

Generate()

phsp Generate	(	B0	,
		Events.	begin(),
		Events.	end()
	)

NArg()

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

Examples:

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

parse()

cmd parse	(	argv	,
		argc
	)

Variable Documentation

argc

hydra::Vector4R char** argc

Initial value:

{


    size_t  nentries   = 0

B0_mass

double B0_mass = 5.27955

Examples:

phsp_chain.inl.

breit_wigner

auto breit_wigner

Initial value:

= hydra::wrap_lambda(
            [] __hydra_dual__ (unsigned int npar, const hydra::Parameter* params, P1 p1, P2 p2, P3 p3 ){

        double mass  = params[0];
        double width = params[1];

        auto   m = (p2+p3).mass();

        double dmass2 = (m-mass)*(m-mass);
        double width2   = width*width;

        double denominator = dmass2 + width2/4.0;
        return ((width2)/4.0)/denominator;

    }, M0, W0 )

Examples:

phsp_reweighting.inl, and phsp_unweighting_functor.inl.

dalitz_calculator

auto dalitz_calculator

Initial value:

= hydra::wrap_lambda( [] __hydra_dual__ ( P1 p1, P2 p2, P3 p3){

            double   MSq_12 = (p1+p2).mass2();
            double   MSq_13 = (p1+p3).mass2();
            double   MSq_23 = (p2+p3).mass2();

            return hydra::make_tuple(MSq_12, MSq_13, MSq_23);
    })

Examples:

dalitz_plot.C, dalitz_plot.inl, phsp_basic.inl, phsp_reweighting.inl, and phsp_unweighting_functor.inl.

elapsed

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

end

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

Examples:

adaptive_gauss_kronrod.inl, dalitz_plot.C, dalitz_plot.inl, gauss_kronrod.inl, and sample_distribution.inl.

Jpsi_mass

double Jpsi_mass = 3.0969

Examples:

phsp_chain.inl.

K_mass

double K_mass = 0.493677

Examples:

phsp_chain.inl.

M0

hydra::Parameter M0

Initial value:

= hydra::Parameter::Create()
        .Name("mass" ).Value(0.8).Error(0.001).Limits(0.7, 0.9)

model

auto model

Initial value:

= hydra::make_pdf( breit_wigner,
            hydra::PhaseSpaceIntegrator<3, hydra::device::sys_t>(B0_mass, {Jpsi_mass, K_mass, pi_mass }, 100000))

nentries

nentries = NArg.getValue()

phsp

hydra::PhaseSpace<3> phsp(B0_mass,{Jpsi_mass, K_mass, pi_mass })

Examples:

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

pi_mass

double pi_mass = 0.13957061

Examples:

phsp_chain.inl.

start

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

Examples:

adaptive_gauss_kronrod.inl, basic_fit.inl, dalitz_plot.C, dalitz_plot.inl, fit_gaussian.C, gauss_kronrod.inl, phsp_averaging_functor.inl, phsp_basic.inl, phsp_chain.inl, phsp_reweighting.inl, phsp_unweighting_functor.inl, plain_mc.inl, pseudo_experiment.inl, sample_distribution.inl, and simultaneous_fit.inl.

try

try

Initial value:

{

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

W0

hydra::Parameter W0

Initial value:

= hydra::Parameter::Create()
        .Name("width").Value(0.05).Error(0.001).Limits(0.04, 0.06)