Detailed Description

Generate momenta for an event using microcanonical phase space sampling.

Created by KVClassFactory on Fri Apr 17 10:19:02 2015 Author: John Frankland,,,

We use ROOT class TGenPhaseSpace to sample the phase space for each event. This corresponds to sampling the microcanonical phase space for a gas composed of the nuclei of the break-up channel with a given total kinetic energy determined by energy balance.

Example of use

 
KVNucleus CN(6,12);
CN.SetExcitEnergy(20); // minimum = 7.27 MeV
 
KVEvent decay;
KVNucleus* n = decay.AddParticle();
n->SetZandA(2,4);
n = decay.AddParticle();
n->SetZandA(2,4);
n = decay.AddParticle();
n->SetZandA(2,4);
 
KVGenPhaseSpace gps;
if(!gps.SetBreakUpChannel(CN, &decay))
{
  printf("Break-up channel is not allowed\n");
  return;
}
 
Double_t etot = gps.GetAvailableEnergy();
TH1F* h1 = new TH1F("h1","Kinetic energy of alpha particle 3",200,0,etot*2./3.);
h1->Sumw2();
 
Int_t nevents = 10000;
while(nevents--){
  wgt = gps.Generate();
  h1->Fill(n->GetKE(),wgt);
}
 
h1->Draw();

Examples: KVGenPhaseSpace_example.C, and KVGenPhaseSpace_example2.C.

Definition at line 61 of file KVGenPhaseSpace.h.

#include <KVGenPhaseSpace.h>

Inheritance diagram for KVGenPhaseSpace:

[legend]

Public Member Functions
	KVGenPhaseSpace ()
	Default constructor. More...

	KVGenPhaseSpace (const Char_t name, const Char_t title="")
	Constructor with name and title. More...

virtual	~KVGenPhaseSpace ()
	Destructor. More...

void	Copy (TObject &obj) const

virtual Double_t	Generate ()

Double_t	GetAvailableEnergy () const

Bool_t	SetBreakUpChannel (const KVNucleus &CN, KVEvent *e)

Public Member Functions inherited from KVBase
	KVBase ()
	Default constructor. More...

	KVBase (const Char_t name, const Char_t title="")
	Ctor for object with given name and type. More...

	KVBase (const KVBase &)
	copy ctor More...

virtual	~ KVBase ()

virtual void	Clear (Option_t *opt="")
	Clear object properties : name, type/title, number, label. More...

const Char_t *	GetLabel () const

UInt_t	GetNumber () const

UInt_t	GetNumberOfObjects () const

virtual TObject *	GetObject () const

virtual const Char_t *	GetType () const

Bool_t	HasLabel () const

virtual Bool_t	IsCalled (const Char_t *name) const

Bool_t	IsLabelled (const Char_t *l) const

virtual Bool_t	IsType (const Char_t *typ) const

virtual void	List ()

KVBase &	operator= (const KVBase &)
	copy assignment operator More...

virtual void	Print (Option_t *option="") const

Double_t	ProtectedGetX (const TF1 *func, Double_t val, int &status, Double_t xmin=0.0, Double_t xmax=0.0) const

void	SetLabel (const Char_t *lab)

virtual void	SetNumber (UInt_t num)

virtual void	SetType (const Char_t *str)

Private Member Functions
Bool_t	CheckBreakUpChannel ()

void	init ()
	default initialisations More...

virtual void	InitialiseMCSampler ()

Private Attributes
KVNucleus	fCompound
	initial nucleus which undergoes break-up More...

Double_t	fEtot
	available kinetic energy for decay More...

KVEvent *	fEvent
	break-up channel More...

std::vector< Double_t >	fMasses
	masses of nuclei in break-up channel More...

TObject *	fMCSampler
	Monte-Carlo phase space sampler. More...

Int_t	fMult
	multiplicity of channel More...

Bool_t	fOK
	ready to calculate or not More...

Additional Inherited Members
Public Types inherited from KVBase
enum	EKaliVedaBits { kIsKaliVedaObject = BIT(23) }

Static Public Member Functions inherited from KVBase
static Bool_t	AreEqual (Double_t x, Double_t y, Long64_t maxdif=1)
	Comparison between two 64-bit floating-point values. More...

static void	BackupFileWithDate (const Char_t *path)

static void	CombineFiles (const Char_t file1, const Char_t file2, const Char_t *newfilename, Bool_t keep=kTRUE)

static void	Deprecated (const char method, const char advice)

static Bool_t	FindClassSourceFiles (const Char_t class_name, KVString &imp_file, KVString &dec_file, const Char_t dir_name=".")

static Bool_t	FindExecutable (TString &exec, const Char_t *path="$(PATH)")

static const Char_t *	FindFile (const Char_t *search, TString &wfil)

static const Char_t *	GetBINDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetDATABASEFilePath ()

static const Char_t *	GetDATADIRFilePath (const Char_t *namefile="")

static Bool_t	GetDataSetEnv (const Char_t dataset, const Char_t type, Bool_t defval)

static const Char_t *	GetDataSetEnv (const Char_t dataset, const Char_t type, const Char_t *defval)

static Double_t	GetDataSetEnv (const Char_t dataset, const Char_t type, Double_t defval)

static const Char_t *	GetETCDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetExampleFilePath (const Char_t library, const Char_t namefile)
	Return full path to example file for given library (="KVMultiDet", "BackTrack", etc.) More...

static const Char_t *	GetINCDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetKVBuildDate ()
	Returns KaliVeda build date. More...

static const Char_t *	GetKVBuildDir ()
	Returns top-level directory used for build. More...

static const Char_t *	GetKVBuildTime ()
	Returns KaliVeda build time. More...

static const Char_t *	GetKVBuildType ()
	Returns KaliVeda build type (cmake build: Release, Debug, RelWithDebInfo, ...) More...

static const Char_t *	GetKVBuildUser ()
	Returns username of person who performed build. More...

static const Char_t *	GetKVSourceDir ()
	Returns top-level directory of source tree used for build. More...

static const Char_t *	GetKVVersion ()
	Returns KaliVeda version string. More...

static const Char_t *	GetLIBDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetListOfPlugins (const Char_t *base)

static const Char_t *	GetListOfPluginURIs (const Char_t *base)

static const Char_t *	GetPluginURI (const Char_t base, const Char_t plugin)

static void	GetTempFileName (TString &base)

static const Char_t *	GetTEMPLATEDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetWORKDIRFilePath (const Char_t *namefile="")

static const Char_t *	gitBranch ()
	Returns git branch of sources. More...

static const Char_t *	gitCommit ()
	Returns last git commit of sources. More...

static void	InitEnvironment ()

static bool	is_gnuinstall ()

static Bool_t	IsThisAPlugin (const TString &uri, TString &base)

static TPluginHandler *	LoadPlugin (const Char_t base, const Char_t uri="0")

static Bool_t	OpenContextMenu (const char method, TObject obj, const char *alt_method_name="")

static void	OpenTempFile (TString &base, std::ofstream &fp)

static void	PrintSplashScreen ()
	Prints welcome message and infos on version etc. More...

static Bool_t	SearchAndOpenKVFile (const Char_t name, KVSQLite::database &dbfile, const Char_t kvsubdir="")

static Bool_t	SearchAndOpenKVFile (const Char_t name, std::ifstream &file, const Char_t kvsubdir="", KVLockfile *locks=0)

static Bool_t	SearchAndOpenKVFile (const Char_t name, std::ofstream &file, const Char_t kvsubdir="", KVLockfile *locks=0)

static Bool_t	SearchKVFile (const Char_t name, TString &fullpath, const Char_t kvsubdir="")

static const Char_t *	WorkingDirectory ()

Constructor & Destructor Documentation

◆ KVGenPhaseSpace() [1/2]

KVGenPhaseSpace::KVGenPhaseSpace ( )

Default constructor.

Definition at line 49 of file KVGenPhaseSpace.cpp.

◆ KVGenPhaseSpace() [2/2]

KVGenPhaseSpace::KVGenPhaseSpace	(	const Char_t *	name,
		const Char_t *	title = `""`
	)

Constructor with name and title.

Definition at line 61 of file KVGenPhaseSpace.cpp.

◆ ~KVGenPhaseSpace()

KVGenPhaseSpace::~KVGenPhaseSpace ( )

virtual

Destructor.

Definition at line 72 of file KVGenPhaseSpace.cpp.

Member Function Documentation

◆ CheckBreakUpChannel()

Bool_t KVGenPhaseSpace::CheckBreakUpChannel ( )

private

We check that:

mass and charge are conserved
the excitation energy of the compound is greater than the channel Q-value (taking into account any excitation energy of nuclei in exit channel)

Definition at line 180 of file KVGenPhaseSpace.cpp.

◆ Copy()

void KVGenPhaseSpace::Copy ( TObject & obj ) const

virtual

This method copies the current state of 'this' object into 'obj' You should add here any member variables, for example: (supposing a member variable KVGenPhaseSpace::fToto) CastedObj.fToto = fToto; or CastedObj.SetToto( GetToto() );

Reimplemented from KVBase.

Definition at line 89 of file KVGenPhaseSpace.cpp.

◆ Generate()

Double_t KVGenPhaseSpace::Generate ( )

virtual

Generate 1 event for this break-up channel. The weight associated to the event is returned. We also store the value in the list of parameters of the KVEvent object pointed to by fEvent: fEvent->GetParameters()->SetValue("PHASESPACE_WEIGHT", weight);

NB: in order to get e.g. correct kinetic energy distributions for particles, this weight MUST be used as e.g. the filling weight for any histograms (see example in Class Description)

Examples: KVGenPhaseSpace_example.C, and KVGenPhaseSpace_example2.C.

Definition at line 148 of file KVGenPhaseSpace.cpp.

◆ GetAvailableEnergy()

Double_t KVGenPhaseSpace::GetAvailableEnergy ( ) const

inline

Return total kinetic energy available for break-up particles

Examples: KVGenPhaseSpace_example2.C.

Definition at line 87 of file KVGenPhaseSpace.h.

◆ init()

void KVGenPhaseSpace::init ( void )

private

default initialisations

Definition at line 15 of file KVGenPhaseSpace.cpp.

◆ InitialiseMCSampler()

void KVGenPhaseSpace::InitialiseMCSampler ( )

privatevirtual

initialise the TGenPhaseSpace for the decay channel We store the name of the generator in the list of parameters of the KVEvent object pointed to by fEvent: fEvent->GetParameters()->SetValue("PHASESPACE_GENERATOR","TGenPhaseSpace");

Definition at line 32 of file KVGenPhaseSpace.cpp.

◆ SetBreakUpChannel()

Bool_t KVGenPhaseSpace::SetBreakUpChannel	(	const KVNucleus &	CN,
		KVEvent *	e
	)

Define the break-up channel, i.e. the initial compound nucleus with its kinematics & excitation energy, and a list of nuclei to be produced by the decay. Nuclei produced by decay may be excited.

Returns kTRUE if decay is energetically allowed, kFALSE if not.

Examples: KVGenPhaseSpace_example.C, and KVGenPhaseSpace_example2.C.

Definition at line 111 of file KVGenPhaseSpace.cpp.

Member Data Documentation

◆ fCompound

KVNucleus KVGenPhaseSpace::fCompound

private

initial nucleus which undergoes break-up

Definition at line 67 of file KVGenPhaseSpace.h.

◆ fEtot

Double_t KVGenPhaseSpace::fEtot

private

available kinetic energy for decay

Definition at line 72 of file KVGenPhaseSpace.h.

◆ fEvent

KVEvent* KVGenPhaseSpace::fEvent

private

break-up channel

Definition at line 68 of file KVGenPhaseSpace.h.

◆ fMasses

std::vector<Double_t> KVGenPhaseSpace::fMasses

private

masses of nuclei in break-up channel

Definition at line 70 of file KVGenPhaseSpace.h.

◆ fMCSampler

TObject* KVGenPhaseSpace::fMCSampler

private

Monte-Carlo phase space sampler.

Definition at line 66 of file KVGenPhaseSpace.h.

◆ fMult

Int_t KVGenPhaseSpace::fMult

private

multiplicity of channel

Definition at line 69 of file KVGenPhaseSpace.h.

◆ fOK

Bool_t KVGenPhaseSpace::fOK

private

ready to calculate or not

Definition at line 71 of file KVGenPhaseSpace.h.

Detailed Description

Example of use

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ KVGenPhaseSpace() [1/2]

◆ KVGenPhaseSpace() [2/2]

◆ ~KVGenPhaseSpace()

Member Function Documentation

◆ CheckBreakUpChannel()

◆ Copy()

◆ Generate()

◆ GetAvailableEnergy()

◆ init()

◆ InitialiseMCSampler()

◆ SetBreakUpChannel()

Member Data Documentation

◆ fCompound

◆ fEtot

◆ fEvent

◆ fMasses

◆ fMCSampler

◆ fMult

◆ fOK