Functions
matrix	ReadDensityMatrix ("fort.925", 3, kTRUE, 0.49, 0.51)
matrix	SetMatrix (12)
matrix	Do ("mt = sqrt(_z_z + _x_x)")
matrix	Do ("[dilepton]->SetPx(p1)")
matrix	Do ("[dilepton]->SetPy(p2)")
matrix	Do ("[dilepton]->SetPz(p3)")
matrix	Do ("[dilepton]->SetM(_x)")
*	makeDynamicData () -> GetBatchParticle("boost"))=boost
matrix	Do ("foreach(); []->Boost(boost)")
dilepton_decay	AddPrologueBulk (matrix)
my_reaction	Do (histo1, "_x = ([e+]+[e-])->M()")
my_reaction	Do (histo2, "_x = ([e+]+[e-])->Rapidity(); _y = ([e+]+[e-])->Pt()")
my_reaction	Loop (1000000)
histo1	Draw ()

Variables
PParticle *	seed = new PParticle("dilepton", 0., 0., 0.)
PParticle	boost (0, 0, 0, -4.0, 1.0)
PParticle *	ep = new PParticle("e+")
PParticle *	em = new PParticle("e-")
PParticle *	p [] = {seed, em, ep}
PChannel	dilepton_decay (p, 2)
PChannel *	c [] ={&dilepton_decay}
PReaction	my_reaction (c, "output", 1)
TH1F *	histo1 = new TH1F("histo1", "ee mass", 100, 0.01, 1.5)
TH2F *	histo2 = new TH2F("histo2","Y vs. Pt", 100, -4.0, 6.0, 300, 0.,1.)

Function Documentation

◆ AddPrologueBulk()

dilepton_decay AddPrologueBulk ( matrix )

◆ Do() [1/8]

matrix Do ( "->SetM(_x)" [dilepton] )

◆ Do() [2/8]

matrix Do ( "->SetPx(p1)" [dilepton] )

◆ Do() [3/8]

matrix Do ( "->SetPy(p2)" [dilepton] )

◆ Do() [4/8]

matrix Do ( "->SetPz(p3)" [dilepton] )

◆ Do() [5/8]

matrix Do ( "foreach(*); ->Boost(boost)" [ *] )

◆ Do() [6/8]

matrix Do ( )

◆ Do() [7/8]

my_reaction Do ( histo1 )

◆ Do() [8/8]

my_reaction Do ( histo2 )

◆ Draw()

histo1 Draw ( )

◆ Loop()

my_reaction Loop ( 1000000 )

◆ makeDynamicData()

* makeDynamicData ( ) -> GetBatchParticle("boost"))=boost

◆ ReadDensityMatrix()

matrix ReadDensityMatrix	(	"fort.925"	,
		3	,
		kTRUE	,
		0.	49,
		0.	51 )

◆ SetMatrix()

matrix SetMatrix ( 12 )

Variable Documentation

◆ boost

PParticle boost(0, 0, 0, -4.0, 1.0)	(	0	,
		0	,
		0	,
		-4.	0,
		1.	0 )

◆ c

PChannel* c[] ={&dilepton_decay}

Definition at line 52 of file use_DensityMatrix.C.

◆ dilepton_decay

PChannel dilepton_decay(p, 2)	(	p	,
		2	)

◆ em

PParticle* em = new PParticle("e-")

Definition at line 47 of file use_DensityMatrix.C.

◆ ep

PParticle* ep = new PParticle("e+")

Definition at line 46 of file use_DensityMatrix.C.

◆ histo1

TH1F* histo1 = new TH1F("histo1", "ee mass", 100, 0.01, 1.5)

Definition at line 56 of file use_DensityMatrix.C.

◆ histo2

TH2F* histo2 = new TH2F("histo2","Y vs. Pt", 100, -4.0, 6.0, 300, 0.,1.)

Definition at line 58 of file use_DensityMatrix.C.

◆ my_reaction

PReaction my_reaction(c, "output", 1)	(	c	,
		"output"	,
		1	)

◆ p

PParticle* p[] = {seed, em, ep}

Definition at line 48 of file use_DensityMatrix.C.

◆ seed

PParticle* seed = new PParticle("dilepton", 0., 0., 0.)

Definition at line 18 of file use_DensityMatrix.C.

Functions

Variables