Detailed Description

Calculate elastic scattering count rates in the detectors of multidetector arrays ,.

Created by KVClassFactory on Fri Nov 20 2015 Author: John Frankland

Use this class to calculate the count rates & energy losses of elastically scattered nuclei in the detectors of a multidetector array corresponding to a given angular range.

For details on setting up and fine-tuning the simulation, see class KVElasticScatter.

Running the calculation, obtaining the results

To perform a calculation, the minimum you need to do is:

KVMultiDetArray::MakeMultiDetector("name_of_dataset"); // build geometry
 
KVElasticCountRates ecr(run_number); // choose run => define projectile, target etc.
ecr.SetAngularRange(theta_min,theta_max,phi_min,phi_max); // all detectors in range will be sampled
ecr.CaclulateScattering(10000);

Make sure the number of sampled points (default=10000) is large enough for accurate determination of count rates.

At the end of calculation, we print infos for all detectors hit by elastic particles, these are the count rates for a nominal beam intensity of 10**7 particles per second, the mean energy loss in the detector, the total integrated cross-section per detector in barn, the number of particles per second per unit area hitting each detector, and the total energy per second per unit area deposited in each detector:

esa.CalculateScattering(10000);
 
SI_0101  :  N=  204.42/sec.    <E>=  890.3 MeV   Tot.Xsec=6.687E+01 barn    fluence=2.711E+01/sec./cm**2     dissip.=2.413E+04 MeV/sec./cm**2
SI_0201  :  N=   61.37/sec.    <E>=  709.7 MeV   Tot.Xsec=2.008E+01 barn    fluence=1.413E+01/sec./cm**2     dissip.=1.002E+04 MeV/sec./cm**2
SI_0301  :  N=   15.29/sec.    <E>=  706.4 MeV   Tot.Xsec=5.002E+00 barn    fluence=2.390E+00/sec./cm**2     dissip.=1.688E+03 MeV/sec./cm**2
SI_0401  :  N=    6.33/sec.    <E>=  699.2 MeV   Tot.Xsec=2.070E+00 barn    fluence=1.582E+00/sec./cm**2     dissip.=1.106E+03 MeV/sec./cm**2

If you want to see the count rates for a different beam intensity, call

ecr.PrintResults(5.e+06)

with the required number of projectiles per second.

For each detector hit we fill three histograms:

OBJ: TH1F  SI_0223  Eloss in SI_0223 : 0 at: 0x10898980
OBJ: TH1F  SI_0223_dW  Solid angle of SI_0223 : 0 at: 0x10899540
OBJ: TH2F  SI_0223_map Map of SI_0223 : 0 at: 0x10e2e7d0

Each histogram can be obtained using

ecr.GetHistos().FindObject("SI_0223_map")

The different informations for each detector can be obtained using

ecr.GetDetector("SI_0201").count_rate
(double)6.13702995974404857e+01
 
ecr.GetDetector("SI_0201").mean_energy
(double)7.09659910827807721e+02
 
ecr.GetDetector("SI_0201").intXsec
(double)2.00758282694385635e+01

See KVElasticCountRate class for details.

Warning: The values returned correspond to the last beam intensity value given to PrintResults().

Definition at line 115 of file KVElasticCountRates.h.

#include <KVElasticCountRates.h>

Inheritance diagram for KVElasticCountRates:

[legend]

Public Member Functions
	KVElasticCountRates (int run)

void	FillHistograms (const KVNameValueList *, double theta, double phi, double xsec)

KVElasticCountRate	GetDetector (const std::string &name)

const KVHashList &	GetHistos () const

void	PrintResults (Double_t beam_intensity=1.e+07)
	Print mean energy deposit & counting rate for given beam intensity in particles per second. More...

KVNameValueList	PutResultsInList (Double_t beam_intensity=1.e+07)
	Print mean energy deposit & counting rate for given beam intensity in particles per second. More...

void	SetAngularRange (Double_t theta_min=0, Double_t theta_max=180, Double_t phi_min=0, Double_t phi_max=360)

void	SetBeamIntensity (Double_t bi)

Public Member Functions inherited from KVElasticScatter
	KVElasticScatter (Int_t run)
	Default constructor. More...

virtual	~ KVElasticScatter ()

void	CalculateScattering (Int_t N)
	Perform scattering 'N' times. More...

TH1F *	GetDepth ()
	Return pointer to histogram of 'depth' of scattering point in target (in mg/cm2) More...

Int_t	GetEbinning (void)
	Returns the number of bins of the GetEnergy histogram. More...

TH1F *	GetEnergy ()
	Return pointer to energy loss histogram for chosen detector (in MeV) More...

TH1F *	GetEnergy (const Char_t *label)
	Energy loss in detector with given 'label' through which scattered particle passes. More...

TH1F *	GetEnergy (Int_t index)

auto	GetKinematics ()

Int_t	GetNDets () const
	Returns the number of detectors crossed by the scattered particle. More...

void	GetTargetNucleusFromTargetLayer (Bool_t yes=kTRUE)

TH1F *	GetTheta ()
	Return pointer to polar angle distribution of scattered particle (in degrees) More...

void	Print () const
	Print details of calculation to be performed. More...

void	SetBeamDirection (const TVector3 &beam_dir)

void	SetDetector (const Char_t *det)

void	SetDetectorResolution (double fwhm, double eref)

void	SetEbinning (Int_t nbins)

void	SetExcitation (Double_t ex)

void	SetKinematicSolution (KV2Body::kinematic_solution ik)

void	SetNucleusOfInterest (KV2Body::nucleus_of_interest noi)

void	SetOutGoingProjectileLike (const KVNucleus &opl)

void	SetProjectileNucleus (const KVNucleus &pn)

void	SetTargetNucleusForScattering (const KVNucleus &tn)

void	SetTargetScatteringLayer (const Char_t *name)

Constructor & Destructor Documentation

◆ KVElasticCountRates()

KVElasticCountRates::KVElasticCountRates ( int run )

inline

Definition at line 132 of file KVElasticCountRates.h.

Member Function Documentation

◆ FillHistograms()

void KVElasticCountRates::FillHistograms	(	const KVNameValueList *	dets,
		double	theta,
		double	phi,
		double	xsec
	)

parse the list dets fill histograms with energy loss for all detectors clear the detector energy losses delete the list

Definition at line 104 of file KVElasticCountRates.cpp.

◆ GetDetector()

KVElasticCountRate KVElasticCountRates::GetDetector ( const std::string & name )

inline

Definition at line 150 of file KVElasticCountRates.h.

◆ GetHistos()

const KVHashList& KVElasticCountRates::GetHistos ( ) const

inline

Definition at line 142 of file KVElasticCountRates.h.

◆ PrintResults()

void KVElasticCountRates::PrintResults ( Double_t beam_intensity = 1.e+07 )

Print mean energy deposit & counting rate for given beam intensity in particles per second.

Definition at line 202 of file KVElasticCountRates.cpp.

◆ PutResultsInList()

KVNameValueList KVElasticCountRates::PutResultsInList ( Double_t beam_intensity = 1.e+07 )

Print mean energy deposit & counting rate for given beam intensity in particles per second.

Definition at line 243 of file KVElasticCountRates.cpp.

◆ SetAngularRange()

void KVElasticCountRates::SetAngularRange	(	Double_t	theta_min = `0`,
		Double_t	theta_max = `180`,
		Double_t	phi_min = `0`,
		Double_t	phi_max = `360`
	)

Definition at line 33 of file KVElasticCountRates.cpp.

◆ SetBeamIntensity()

void KVElasticCountRates::SetBeamIntensity ( Double_t bi )

inline

Definition at line 137 of file KVElasticCountRates.h.

Detailed Description

Running the calculation, obtaining the results

Public Member Functions

Constructor & Destructor Documentation

◆ KVElasticCountRates()

Member Function Documentation

◆ FillHistograms()

◆ GetDetector()

◆ GetHistos()

◆ PrintResults()

◆ PutResultsInList()

◆ SetAngularRange()

◆ SetBeamIntensity()