#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>
#include "TVector3.h"
#include "TLorentzVector.h"
#include "TH1D.h"
#include "TH2D.h"
#include "TH3D.h"
#include "TGenPhaseSpace.h"
#include "TCanvas.h"
#include "TRandom3.h"
#include "TFile.h"
#include "TTree.h"
#include "TStyle.h"
#include "TROOT.h"
#include "TMath.h"

Functions
const string	header (" ***********************************\n simulating decay into K- pi+ pi- \n ***********************************\n")
TLorentzVector	makeBeam ()
TVector3	makeVertex ()
bool	writePwa2000Ascii (ostream &out, const TLorentzVector &beam, TGenPhaseSpace &event)
bool	writeComGeantAscii (ostream &out, const TLorentzVector &beam, TGenPhaseSpace &event, const TVector3 &vertexpos, const TLorentzVector &recoilproton, bool formated=true)
ostream &	progressIndicator (const long currentPos, const long nmbTotal, const int nmbSteps=10, ostream &out=cout)
float	GetPhi (const TLorentzVector &beam, TGenPhaseSpace &event, const TLorentzVector &recoilproton)
float	Calc_t_prime (const TLorentzVector &particle_In, const TLorentzVector &particle_Out)
void	genPhaseSpaceData (const double xMassMin=2.100, const double xMassMax=2.140, const TString &outFileName="2100.2140.genbod.evt", const TString &thetaHistFileName="./hTheta.root", const int nmbEvent=2000, const bool plot=true, const TString &outFileNameComGeant="")

Variables
const double	gProtonMass = 0.938272013
const double	gPionMass = 0.13957018
const double	gKaonMass = 0.493677
const int	geantIds [3] = { 8, 9, 12}
const int	charges [3] = {+1, -1, -1}
Double_t	daughterMasses [3] = {gPionMass, gPionMass, gKaonMass}
const double	gbeampartMass2 = daughterMasses[2]*daughterMasses[2]
const double	gTargetZPos = -30.0
const double	gTargetlength = -40.0
const double	gBeamMomSigma = 1.2
const double	gBeamMom = 191
const double	gBeamDxDz = 0.00026
const double	gBeamDxDzSigma = 0.00010
const double	gBeamDyDz = 0.00001
const double	gBeamDyDzSigma = 0.00018
const double	gBeamOffsetX = 0.0
const double	gBeamOffsetY = 0.0
const double	gBeamSpotSizeX = 1.5
const double	gBeamSpotSizeY = 1.5
const double	tMin = 0.001

Function Documentation

float Calc_t_prime	(	const TLorentzVector &	particle_In,
		const TLorentzVector &	particle_Out
	)

Definition at line 235 of file genPhaseSpaceData.C.

References M2.

Referenced by genPhaseSpaceData().

void genPhaseSpaceData	(	const double	xMassMin = `2.100`,
		const double	xMassMax = `2.140`,
		const TString &	outFileName = `"2100.2140.genbod.evt"`,
		const TString &	thetaHistFileName = `"./hTheta.root"`,
		const int	nmbEvent = `2000`,
		const bool	plot = `true`,
		const TString &	outFileNameComGeant = `""`
	)

Definition at line 256 of file genPhaseSpaceData.C.

References Calc_t_prime(), daughterMasses, gbeampartMass2, GetPhi(), gProtonMass, header(), i, int, makeBeam(), makeVertex(), progressIndicator(), q(), tMin, writeComGeantAscii(), and writePwa2000Ascii().

float GetPhi	(	const TLorentzVector &	beam,
		TGenPhaseSpace &	event,
		const TLorentzVector &	recoilproton
	)

Definition at line 213 of file genPhaseSpaceData.C.

References i.

Referenced by genPhaseSpaceData().

const string header ( " ************************************\n * simulating decay into K- pi+ pi- *\n ************************************\n" )

Referenced by genPhaseSpaceData().

TLorentzVector makeBeam ( )

Definition at line 98 of file genPhaseSpaceData.C.

References gBeamDxDz, gBeamDxDzSigma, gBeamDyDz, gBeamDyDzSigma, gBeamMom, gBeamMomSigma, and gbeampartMass2.

Referenced by genPhaseSpaceData().

TVector3 makeVertex ( )

Definition at line 115 of file genPhaseSpaceData.C.

References gBeamOffsetX, gBeamOffsetY, gBeamSpotSizeX, gBeamSpotSizeY, gTargetlength, and gTargetZPos.

Referenced by genPhaseSpaceData().

ostream& progressIndicator	(	const long	currentPos,
		const long	nmbTotal,
		const int	nmbSteps = `10`,
		ostream &	out = `cout`
	)

Definition at line 197 of file genPhaseSpaceData.C.

References int.

Referenced by calcNewAmps(), createDiffTree(), rpwa::fillTreeFromAmp(), rpwa::fillTreeFromEvt(), fillUdstDataIntoMassBins_example(), genPhaseSpaceData(), rpwa::ampIntegralMatrix::integrate(), main(), plotAmpDiffs(), rpwa::processTree(), rpwa::writeAmpFromTree(), and rpwa::writeEvtFromTree().

bool writeComGeantAscii	(	ostream &	out,
		const TLorentzVector &	beam,
		TGenPhaseSpace &	event,
		const TVector3 &	vertexpos,
		const TLorentzVector &	recoilproton,
		bool	formated = `true`
	)

Definition at line 156 of file genPhaseSpaceData.C.

References geantIds, and i.

Referenced by genPhaseSpaceData().

bool writePwa2000Ascii	(	ostream &	out,
		const TLorentzVector &	beam,
		TGenPhaseSpace &	event
	)

Definition at line 127 of file genPhaseSpaceData.C.

References charges, geantIds, and i.

Referenced by genPhaseSpaceData().

Variable Documentation

const int charges[3] = {+1, -1, -1}

Definition at line 55 of file genPhaseSpaceData.C.

Referenced by writePwa2000Ascii().

Double_t daughterMasses[3] = {gPionMass, gPionMass, gKaonMass}

Definition at line 56 of file genPhaseSpaceData.C.

Referenced by rpwa::diffractivePhaseSpace::BuildDaughterList(), compareAmplitudes(), genPhaseSpaceData(), and rpwa::nBodyPhaseSpaceGen::setDecay().

const double gBeamDxDz = 0.00026

Definition at line 73 of file genPhaseSpaceData.C.

Referenced by makeBeam().

const double gBeamDxDzSigma = 0.00010

Definition at line 74 of file genPhaseSpaceData.C.

Referenced by makeBeam().

const double gBeamDyDz = 0.00001

Definition at line 75 of file genPhaseSpaceData.C.

Referenced by makeBeam().

const double gBeamDyDzSigma = 0.00018

Definition at line 76 of file genPhaseSpaceData.C.

Referenced by makeBeam().

const double gBeamMom = 191

Definition at line 71 of file genPhaseSpaceData.C.

Referenced by makeBeam().

const double gBeamMomSigma = 1.2

Definition at line 67 of file genPhaseSpaceData.C.

Referenced by makeBeam().

const double gBeamOffsetX = 0.0

Definition at line 84 of file genPhaseSpaceData.C.

Referenced by makeVertex().

const double gBeamOffsetY = 0.0

Definition at line 85 of file genPhaseSpaceData.C.

Referenced by makeVertex().

const double gbeampartMass2 = daughterMasses[2]*daughterMasses[2]

Definition at line 59 of file genPhaseSpaceData.C.

Referenced by genPhaseSpaceData(), and makeBeam().

const double gBeamSpotSizeX = 1.5

Definition at line 86 of file genPhaseSpaceData.C.

Referenced by makeVertex().

const double gBeamSpotSizeY = 1.5

Definition at line 87 of file genPhaseSpaceData.C.

Referenced by makeVertex().

const int geantIds[3] = { 8, 9, 12}

Definition at line 54 of file genPhaseSpaceData.C.

Referenced by writeComGeantAscii(), and writePwa2000Ascii().

const double gKaonMass = 0.493677

Definition at line 37 of file genPhaseSpaceData.C.

const double gPionMass = 0.13957018

Definition at line 35 of file genPhaseSpaceData.C.

const double gProtonMass = 0.938272013

Definition at line 34 of file genPhaseSpaceData.C.

Referenced by genPhaseSpaceData().

const double gTargetlength = -40.0

Definition at line 64 of file genPhaseSpaceData.C.

Referenced by makeVertex().

const double gTargetZPos = -30.0

Definition at line 63 of file genPhaseSpaceData.C.

Referenced by makeVertex().

const double tMin = 0.001

Definition at line 90 of file genPhaseSpaceData.C.

Referenced by genPhaseSpaceData(), and rpwa::diffractivePhaseSpace::SettMin().

Functions

Variables

Function Documentation

Variable Documentation