Calculates the reaction position using telescope and tracking data. More...

#include <TTGTIKProcessor.h>

Inheritance diagram for art::crib::TTGTIKProcessor:

Collaboration diagram for art::crib::TTGTIKProcessor:

Public Member Functions
	TTGTIKProcessor ()
	Constructor.

	~TTGTIKProcessor () override
	Default destructor.

void	Init (TEventCollection *col) override
	Initialize the processor with an event collection.

void	Process () override
	Main processing function.

Private Member Functions
	TTGTIKProcessor (const TTGTIKProcessor &rhs)=delete

Double_t	bisection (const TTrack track, const TTelescopeData data)
	Bisection method for calculating reaction position.

	ClassDefOverride (TTGTIKProcessor, 1)
	ROOT class definition macro.

Double_t	GetCMAngle (Double_t ELab, Double_t Ecm, Double_t ALab)
	Recalculate the LAB angle after reconstruction.

Double_t	GetCustomExcitedEnergy (Int_t telID, Double_t Etotal)
	Generate a custom excited state energy.

Double_t	GetEcm_classic_kinematics (Double_t energy, Double_t theta)
	Calculate the center-of-mass energy using classical kinematics.

Double_t	GetEcm_kinematics (Double_t energy, Double_t theta, Double_t low_e, Double_t high_e)
	Calculate the center-of-mass energy using relativistic kinematics.

Double_t	GetEcmFromBeam (Double_t z, const TTrack *track)
	Calculate the center-of-mass energy from beam data.

Double_t	GetEcmFromDetectParticle (Double_t z, const TTrack track, const TTelescopeData data)
	Calculate the center-of-mass energy from detected particle data.

std::pair< Double_t, Double_t >	GetELabALabPair (Double_t z, const TTrack track, const TTelescopeData data)
	Calculate the LAB energy and LAB angle from detected particle data.

Double_t	GetReactionPosition (const TTrack track, const TTelescopeData data)
	Calculate the reaction position along the Z-axis.

Double_t	newton (const TTrack track, const TTelescopeData data)
	(Deprecated) Newton method for calculating reaction position.

TTGTIKProcessor &	operator= (const TTGTIKProcessor &rhs)=delete

Double_t	TargetFunction (Double_t z, const TTrack track, const TTelescopeData data)
	Target function for the bisection (and Newton) method.

Private Attributes
TString	fDetectorParameterName
	Name of the detector parameter (TDetectorParameter)

TClonesArray *	fDetectorPrm
	! Pointer to detector parameter objects (TClonesArray of TDetectorParameter)

Bool_t	fDoCenterPos
	Flag to use the detector center position.

Bool_t	fDoCustom
	Flag to enable custom processing.

Double_t	fExcitedEnergy
	Excited state energy (MeV)

TClonesArray **	fInData
	! Pointer to the input telescope data (TClonesArray of TTelescopeData)

Double_t	fInitialBeamEnergy
	Beam energy immediately after the window (MeV)

TString	fInputColName
	Name of the input telescope data collection (TTelescopeData)

TString	fInputTrackColName
	Name of the input tracking data collection (TTrack)

TClonesArray **	fInTrackData
	! Pointer to the input tracking data (TClonesArray of TTrack)

TClonesArray *	fOutData
	! Pointer to the output reaction information (TClonesArray of TReactionInfo)

TString	fOutputColName
	Name of the output reaction information collection (TReactionInfo)

IntVec_t	fParticleAArray
	Array of mass numbers for reaction particles.

IntVec_t	fParticleZArray
	Array of atomic numbers for reaction particles.

Double_t	fPressure
	Gas pressure in Torr.

TString	fTargetName
	Name of the gas target used in TSrim calculations.

TString	fTargetParameterName
	Name of the target parameter (TTargetParameter)

TClonesArray *	fTargetPrm
	! Pointer to target parameter objects (TClonesArray of TTargetParameter)

Double_t	fTemperature
	Gas temperature in Kelvin.

const Double_t	kEpsilon = 1.0e-3
	Convergence threshold for the bisection method.

const Double_t	kInitialMax = 1000.0
	Initial maximum value for bisection method (mm)

const Double_t	kInitialMin = -250.0
	Initial minimum value for bisection method (mm)

const Int_t	kMaxIteration = 1000
	Maximum number of iterations for the bisection method.

Double_t	M1

Double_t	M2

Double_t	M3

Double_t	M3_default

Double_t	M4

TSrim *	srim
	! TSrim object to calculate energy loss

Detailed Description

Calculates the reaction position using telescope and tracking data.

This class read two input data from a TClonesArray("art::crib::TTelescopeData") and a TClonesArray("art::TTrack"). It applies Thick Gas Target Method and get reaction position and Ecm.

Example Steering File

Processor:
  - name: MyTTGTIKProcessor
    type: art::crib::TTGTIKProcessor
    parameter:
      DetectorParameter: prm_detectors  # [TString] Name of the telescope parameter collection (detector parameters)
      ExcitedEnergy: -1  # [Double_t] Excited state energy (MeV); use a negative value if not applicable
      InitialBeamEnergy: 0  # [Double_t] Initial beam energy (in MeV) immediately after the exit window
      InputCollection: tel  # [TString] Input collection of telescope data objects (derived from TTelescopeData)
      InputTrackCollection: track  # [TString] Input collection of tracking data objects (derived from TTrack)
      OutputCollection: result  # [TString] Output collection containing reaction reconstruction information using the TGTIK method
      OutputTransparency: 0  # [Bool_t] Output is persistent if false (default)
      ParticleAArray: []  # [IntVec_t] Array of mass numbers (A) for reaction particles
      ParticleZArray: []  # [IntVec_t] Array of atomic numbers (Z) for reaction particles
      TargetName: ""  # [TString] Name of the target gas (used in TSrim calculation)
      TargetParameter: prm_targets  # [TString] Name of the target parameter collection
      TargetPressure: 0  # [Double_t] Target gas pressure in Torr
      TargetTemperature: 0  # [Double_t] Target gas temperature in Kelvin
      UseCenterPosition: 0  # [Bool_t] Flag to use the detector's center position (useful when the DSSSD is not operational)
      UseCustomFunction: 0  # [Bool_t] Flag to enable custom processing functions for additional corrections
      Verbose: 1  # [Int_t] verbose level (default 1 : non quiet)

Warning: UseCustomFunction is designed for specific analysis, currently for 26Si(a, p) analysis.

Kinematics Calculation

This processor is using classical (non-relativistic) kinematics. Okawa note the relationship. See my master thesis (Japanese) for the details.

Consider Ion1 + Ion2 -> Ion3 + Ion4 reaction, that is, Ion2(Ion1, Ion3)Ion4 reaction.

\( M \): mass
\( v \): velocity
\( v_2 = 0 \): for the target ion
\( {}_{\mathrm{CM}} \): center-of-mass system
\( \theta\): scatter angle of Ion3
\( Q \): Q-value

Based on the relationship between the center-of-mass frame and the laboratory frame,

\[ \begin{align} v_{\mathrm{CM}} &= \frac{M_2}{M_1 + M_2} v_1 \\ v_{1\mathrm{CM}} &= v_1 - v_{\mathrm{CM}} = \frac{M_2}{M_1}v_{\mathrm{CM}} \\ v_{2\mathrm{CM}} &= 0 - v_{\mathrm{CM}} = -v_{\mathrm{CM}} \end{align} \]

The speed of the center of mass, \( v_{\mathrm{CM}} \), and kinetic energy, \( E_{\mathrm{CM}} \) can be expressed as follows based on the definition of the center of mass.

\[ \begin{align} E_{\mathrm{CM}} &= \frac{M_2(M_1 + M_2)}{2M_1} v_{\mathrm{CM}}^2 \\ &= \alpha v_{\mathrm{CM}}^2 \end{align} \]

Here, \( \alpha \) is defined as:

\[ \begin{align} \alpha = \frac{M_2(M_1 + M_2)}{2M_1} \end{align} \]

Based on the relationship between the reaction Q-value and the conservation laws,

\[ \begin{align} \frac{1}{2}M_1 v_{1\mathrm{CM}}^2 + \frac{1}{2}M_2 v_{2\mathrm{CM}}^2 + Q &= \frac{1}{2}M_3 v_{3\mathrm{CM}}^2 + \frac{1}{2}M_4 v_{4\mathrm{CM}}^2 \\ M_1 v_{1\mathrm{CM}} &= M_2 v_{2\mathrm{CM}} \\ M_3 v_{3\mathrm{CM}} &= M_4 v_{4\mathrm{CM}} \\ \end{align} \]

Cancelling \( v_{4\mathrm{CM}} \), we obtain:

\[ \begin{align} \alpha v_{\mathrm{CM}}^2 + Q &= \frac{M_3(M_3 + M_4)}{2M_4}v_{3\mathrm{CM}}^2 \\ \alpha v_{\mathrm{CM}}^2 + Q &= \beta v_{3\mathrm{CM}}^2 \end{align} \]

Here, \( \beta \) is defined as:

\[ \begin{align} \beta = \frac{M_3(M_3+M_4)}{2M_4} \end{align} \]

For Ion 3, the velocity component perpendicular to the z-axis remains unchanged between the LAB and CM frames, while the parallel component can be transformed using vcm. Therefore,

\[ \begin{align} v_{3\mathrm{CM}}\sin{\theta_{\mathrm{CM}}} &= v_3\sin{\theta} \\ v_{3\mathrm{CM}}\cos{\theta_{\mathrm{CM}}} &= v_3\cos{\theta} - v_{\mathrm{CM}} \end{align} \]

Cancelling \( \theta_{\mathrm{CM}} \) from these two equations, we obtain:

\[ \begin{align} v_{3\mathrm{cm}}^2 &= v_3^2 -2v_3 v_{\mathrm{CM}}\cos\theta + v_{\mathrm{cm}}^2 \end{align} \]

\( v_3 \) and \( \theta \) can be considered observable quantities. Cancelling \( v_{3\mathrm{CM}} \), we obtain:

\[ \begin{align} \alpha v_{\mathrm{CM}}^2 + Q = \beta\left(v_3^2 -2v_3v_{\mathrm{cm}}\cos\theta + v_{\mathrm{CM}}^2\right) \\ (\alpha - \beta)v_{\mathrm{CM}}^2 + 2(\beta v_3 \cos\theta)v_{\mathrm{CM}} + (Q - \beta v_3^2) = 0 \end{align} \]

By solving this quadratic equation, various quantities can be obtained. These relationships are utilized in the methods of this class.

In inverse kinematics, ion 4 is detected. Therefore, the indices 3 and 4 should be swapped, and the transformation \( \theta_{\mathrm{CM}} \rightarrow \pi - \theta_{\mathrm{CM}} \) should be applied.

Definition at line 162 of file TTGTIKProcessor.h.

Constructor & Destructor Documentation

◆ TTGTIKProcessor() [1/2]

art::crib::TTGTIKProcessor::TTGTIKProcessor ( )

Constructor.

Definition at line 27 of file TTGTIKProcessor.cc.

◆ ~TTGTIKProcessor()

art::crib::TTGTIKProcessor::~TTGTIKProcessor ( )

override

Default destructor.

Definition at line 68 of file TTGTIKProcessor.cc.

◆ TTGTIKProcessor() [2/2]

art::crib::TTGTIKProcessor::TTGTIKProcessor ( const TTGTIKProcessor & rhs )

privatedelete

Member Function Documentation

◆ bisection()

Double_t art::crib::TTGTIKProcessor::bisection	(	const TTrack *	track,
		const TTelescopeData *	data )

private

Bisection method for calculating reaction position.

Parameters

track	Pointer to the tracking data (TTrack).
data	Pointer to the telescope data (TTelescopeData).

Returns: Calculated reaction Z position (mm).

This function first performs a rough search over the interval [kInitialMin, kInitialMax] to identify a valid subinterval where the target function is defined and exhibits a sign change. It then applies the standard bisection method within that subinterval to converge on a zero.

Definition at line 256 of file TTGTIKProcessor.cc.

◆ ClassDefOverride()

art::crib::TTGTIKProcessor::ClassDefOverride	(	TTGTIKProcessor	,
		1	)

private

ROOT class definition macro.

◆ GetCMAngle()

Double_t art::crib::TTGTIKProcessor::GetCMAngle	(	Double_t	ELab,
		Double_t	Ecm,
		Double_t	ALab )

private

Recalculate the LAB angle after reconstruction.

Parameters

ELab	LAB energy (MeV).
Ecm	Center-of-mass energy (MeV).
ALab	LAB angle (radian).

Returns: Reconstructed LAB angle (radian).

This function calculates the CM scattering angle (θ_cm) based on the detected particle's laboratory energy (ELab) and laboratory scattering angle (ALab), as well as the center-of-mass energy (Ecm). The calculation employs classical kinematics, using coefficients derived from the masses of the reaction participants. In particular:

α is computed from the beam (particle 1) and target (particle 2) masses.
β is computed from the detected particle (particle 4) and the complementary fragment (particle 3) masses.
q-value represents the mass difference between the entrance and exit channels. The function then computes the classical velocities of the detected particle in the LAB frame (v4), the beam velocity in the CM frame (v_cm), and the detected particle's velocity in the CM frame (v4_cm). Finally, the CM angle is obtained from the cosine relationship: θ_cm = acos((v4*cos(ALab) - v_cm) / v4_cm).

Definition at line 629 of file TTGTIKProcessor.cc.

◆ GetCustomExcitedEnergy()

Double_t art::crib::TTGTIKProcessor::GetCustomExcitedEnergy	(	Int_t	telID,
		Double_t	Etotal )

private

Generate a custom excited state energy.

This function is used for custom processing (e.g., handling excited state effects).

Parameters

telID	Identifier for the telescope.
Etotal	Total measured energy (MeV).

Returns: Generated excited state energy (MeV).

This is used for 26Si(a, p)29P analysis. This function uses a custom random generator ROOT file (specific to this analysis) to assign an excited state energy for 29P based on TALYS simulation data. It reads graphs from a directory corresponding to the given telescope ID, computes the ratio of excited state contributions as a function of energy, and then, using a uniform random number, selects one excited state.

NOTE: This function is not designed for generic use and performs file I/O on an event-by-event basis, which is inefficient.

Definition at line 659 of file TTGTIKProcessor.cc.

◆ GetEcm_classic_kinematics()

Double_t art::crib::TTGTIKProcessor::GetEcm_classic_kinematics	(	Double_t	energy,
		Double_t	theta )

private

Calculate the center-of-mass energy using classical kinematics.

Parameters

energy	LAB energy (MeV).
theta	LAB angle (radian).

Returns: Calculated center-of-mass energy (MeV).

This function calculates the center-of-mass energy (Ecm) from the detected particle's laboratory energy and angle using classical (non-relativistic) kinematics. It is used in the GetEcmFromDetectParticle method. The formulas employed here are based on those detailed in Okawa's master thesis.

Definition at line 557 of file TTGTIKProcessor.cc.

◆ GetEcm_kinematics()

Double_t art::crib::TTGTIKProcessor::GetEcm_kinematics	(	Double_t	energy,
		Double_t	theta,
		Double_t	low_e,
		Double_t	high_e )

private

Calculate the center-of-mass energy using relativistic kinematics.

Parameters

energy	LAB energy (MeV).
theta	LAB angle (radian).
low_e	Lower bound for energy (MeV).
high_e	Upper bound for energy (MeV).

Returns: Calculated center-of-mass energy (MeV).

This function is intended to compute the center-of-mass energy (Ecm) from the detected laboratory energy and angle using relativistic formulas. It is designed to be used in the GetEcmFromDetectParticle method. However, the relativistic kinematics implementation is currently not available, and this function always returns kInvalidD.

Definition at line 546 of file TTGTIKProcessor.cc.

◆ GetEcmFromBeam()

Double_t art::crib::TTGTIKProcessor::GetEcmFromBeam	(	Double_t	z,
		const TTrack *	track )

private

Calculate the center-of-mass energy from beam data.

Parameters

z	Reaction position (mm).
track	Pointer to the tracking data (TTrack).

Returns: Calculated center-of-mass energy (MeV).

This function calculates the center-of-mass energy (Ecm) based on the beam's kinematics and its energy loss when traversing the target material. The beam is considered to be the first particle (Z1) of the reaction system [Z1, Z2, Z3, Z4].

The procedure is as follows:

Compute the beam's flight vector from its initial position (at z = 0) to the assumed reaction position z.
Determine the effective path length through the target by applying a sign based on z.
Use the TSrim library to compute the residual energy after energy loss.
Calculate the beam's velocity (β) from its kinetic energy.
Obtain the beam's direction from the track angles (assumed small so that tan(theta) approximations hold), and compute the normalized beta vector.
Boost a beam TLorentzVector (initially at rest with mass M1) using this beta vector so that its energy becomes M1 + kinetic energy.
Construct a target TLorentzVector at rest (mass M2), sum with the beam, and boost to the center-of-mass frame.
Finally, return the total kinetic energy in the center-of-mass frame, which is the sum of (E - M) for both particles.

Definition at line 388 of file TTGTIKProcessor.cc.

◆ GetEcmFromDetectParticle()

Double_t art::crib::TTGTIKProcessor::GetEcmFromDetectParticle	(	Double_t	z,
		const TTrack *	track,
		const TTelescopeData *	data )

private

Calculate the center-of-mass energy from detected particle data.

Parameters

z	Reaction position (mm).
track	Pointer to the tracking data (TTrack).
data	Pointer to the telescope data (TTelescopeData).

Returns: Calculated center-of-mass energy (MeV).

This function computes the center-of-mass energy (Ecm) using the laboratory energy and angle of the detected particle obtained from the assumed reaction position (z). Currently, it employs classical kinematics for the calculation.

Definition at line 453 of file TTGTIKProcessor.cc.

◆ GetELabALabPair()

std::pair< Double_t, Double_t > art::crib::TTGTIKProcessor::GetELabALabPair	(	Double_t	z,
		const TTrack *	track,
		const TTelescopeData *	data )

private

Calculate the LAB energy and LAB angle from detected particle data.

Parameters

z	Reaction position (mm).
track	Pointer to the tracking data (TTrack).
data	Pointer to the telescope data (TTelescopeData).

Returns: A pair containing the LAB energy (MeV) and LAB angle (radian).

This function calculates the LAB energy (ELab) and LAB angle (ALab) of the detected particle, based on an assumed reaction position (z). It retrieves the detector parameters corresponding to the telescope ID, computes the effective detection position using either the strip information or the detector center position (depending on the flag), and then applies the TSrim library to determine the energy loss. The LAB angle is computed as the angle between the track direction and the vector from the reaction position to the detection position.

Definition at line 475 of file TTGTIKProcessor.cc.

◆ GetReactionPosition()

Double_t art::crib::TTGTIKProcessor::GetReactionPosition	(	const TTrack *	track,
		const TTelescopeData *	data )

private

Calculate the reaction position along the Z-axis.

Parameters

track	Pointer to the tracking data (TTrack).
data	Pointer to the telescope data (TTelescopeData).

Returns: Calculated reaction Z position (mm).

This function computes the reaction position (z-coordinate) by employing the bisection method. Although an alternative Newton method is available, the bisection method is preferred for its robustness. Newton method is not implemented yet.

Definition at line 237 of file TTGTIKProcessor.cc.

◆ Init()

void art::crib::TTGTIKProcessor::Init ( TEventCollection * col )

override

Initialize the processor with an event collection.

Parameters

col	Pointer to the event collection.

This function prepares necessary input and output objects, validates parameters, computes mass values for the reaction, initializes the TSrim object for energy loss calculations, and sets up the output collection.

Definition at line 81 of file TTGTIKProcessor.cc.

◆ newton()

Double_t art::crib::TTGTIKProcessor::newton	(	const TTrack *	track,
		const TTelescopeData *	data )

private

(Deprecated) Newton method for calculating reaction position.

Parameters

track	Pointer to the tracking data (TTrack).
data	Pointer to the telescope data (TTelescopeData).

Returns: Calculated reaction Z position (mm).

Currently it is unavailable, return kInvalidD.

Definition at line 246 of file TTGTIKProcessor.cc.

◆ operator=()

TTGTIKProcessor & art::crib::TTGTIKProcessor::operator= ( const TTGTIKProcessor & rhs )

privatedelete

◆ Process()

void art::crib::TTGTIKProcessor::Process ( )

override

Main processing function.

This function processes the input telescope and tracking data to reconstruct the reaction information.

Definition at line 166 of file TTGTIKProcessor.cc.

◆ TargetFunction()

Double_t art::crib::TTGTIKProcessor::TargetFunction	(	Double_t	z,
		const TTrack *	track,
		const TTelescopeData *	data )

private

Target function for the bisection (and Newton) method.

Computes the difference between the beam and detected particle center-of-mass energies.

Parameters

z	Reaction position (mm).
track	Pointer to the tracking data (TTrack).
data	Pointer to the telescope data (TTelescopeData).

Returns: Difference in center-of-mass energy (MeV).

This function calculates two center-of-mass energies based on an assumed reaction position (z):

Ecm(beam): Calculated from beam information.
Ecm(detected): Calculated from detected particle information.

The target function is defined as: f(z) = Ecm(beam) - Ecm(detected)

A zero crossing of this function indicates that the assumed z position corresponds to the true reaction position. The (x, y, z) coordinates are then determined from the tracking data.

Definition at line 361 of file TTGTIKProcessor.cc.

Member Data Documentation

◆ fDetectorParameterName

TString art::crib::TTGTIKProcessor::fDetectorParameterName

private

Name of the detector parameter (TDetectorParameter)

Definition at line 190 of file TTGTIKProcessor.h.

◆ fDetectorPrm

TClonesArray* art::crib::TTGTIKProcessor::fDetectorPrm

private

! Pointer to detector parameter objects (TClonesArray of TDetectorParameter)

Definition at line 199 of file TTGTIKProcessor.h.

◆ fDoCenterPos

Bool_t art::crib::TTGTIKProcessor::fDoCenterPos

private

Flag to use the detector center position.

Definition at line 211 of file TTGTIKProcessor.h.

◆ fDoCustom

Bool_t art::crib::TTGTIKProcessor::fDoCustom

private

Flag to enable custom processing.

Definition at line 210 of file TTGTIKProcessor.h.

◆ fExcitedEnergy

Double_t art::crib::TTGTIKProcessor::fExcitedEnergy

private

Excited state energy (MeV)

Definition at line 209 of file TTGTIKProcessor.h.

◆ fInData

TClonesArray** art::crib::TTGTIKProcessor::fInData

private

! Pointer to the input telescope data (TClonesArray of TTelescopeData)

Definition at line 194 of file TTGTIKProcessor.h.

◆ fInitialBeamEnergy

Double_t art::crib::TTGTIKProcessor::fInitialBeamEnergy

private

Beam energy immediately after the window (MeV)

Definition at line 203 of file TTGTIKProcessor.h.

◆ fInputColName

TString art::crib::TTGTIKProcessor::fInputColName

private

Name of the input telescope data collection (TTelescopeData)

Definition at line 187 of file TTGTIKProcessor.h.

◆ fInputTrackColName

TString art::crib::TTGTIKProcessor::fInputTrackColName

private

Name of the input tracking data collection (TTrack)

Definition at line 188 of file TTGTIKProcessor.h.

◆ fInTrackData

TClonesArray** art::crib::TTGTIKProcessor::fInTrackData

private

! Pointer to the input tracking data (TClonesArray of TTrack)

Definition at line 195 of file TTGTIKProcessor.h.

◆ fOutData

TClonesArray* art::crib::TTGTIKProcessor::fOutData

private

! Pointer to the output reaction information (TClonesArray of TReactionInfo)

Definition at line 196 of file TTGTIKProcessor.h.

◆ fOutputColName

TString art::crib::TTGTIKProcessor::fOutputColName

private

Name of the output reaction information collection (TReactionInfo)

Definition at line 189 of file TTGTIKProcessor.h.

◆ fParticleAArray

IntVec_t art::crib::TTGTIKProcessor::fParticleAArray

private

Array of mass numbers for reaction particles.

Definition at line 208 of file TTGTIKProcessor.h.

◆ fParticleZArray

IntVec_t art::crib::TTGTIKProcessor::fParticleZArray

private

Array of atomic numbers for reaction particles.

Definition at line 207 of file TTGTIKProcessor.h.

◆ fPressure

Double_t art::crib::TTGTIKProcessor::fPressure

private

Gas pressure in Torr.

Definition at line 205 of file TTGTIKProcessor.h.

◆ fTargetName

TString art::crib::TTGTIKProcessor::fTargetName

private

Name of the gas target used in TSrim calculations.

Definition at line 204 of file TTGTIKProcessor.h.

◆ fTargetParameterName

TString art::crib::TTGTIKProcessor::fTargetParameterName

private

Name of the target parameter (TTargetParameter)

Definition at line 191 of file TTGTIKProcessor.h.

◆ fTargetPrm

TClonesArray* art::crib::TTGTIKProcessor::fTargetPrm

private

! Pointer to target parameter objects (TClonesArray of TTargetParameter)

Definition at line 200 of file TTGTIKProcessor.h.

◆ fTemperature

Double_t art::crib::TTGTIKProcessor::fTemperature

private

Gas temperature in Kelvin.

Definition at line 206 of file TTGTIKProcessor.h.

◆ kEpsilon

const Double_t art::crib::TTGTIKProcessor::kEpsilon = 1.0e-3

private

Convergence threshold for the bisection method.

Definition at line 219 of file TTGTIKProcessor.h.

◆ kInitialMax

const Double_t art::crib::TTGTIKProcessor::kInitialMax = 1000.0

private

Initial maximum value for bisection method (mm)

Definition at line 218 of file TTGTIKProcessor.h.

◆ kInitialMin

const Double_t art::crib::TTGTIKProcessor::kInitialMin = -250.0

private

Initial minimum value for bisection method (mm)

Definition at line 217 of file TTGTIKProcessor.h.

◆ kMaxIteration

const Int_t art::crib::TTGTIKProcessor::kMaxIteration = 1000

private

Maximum number of iterations for the bisection method.

Definition at line 220 of file TTGTIKProcessor.h.

◆ M1

Double_t art::crib::TTGTIKProcessor::M1

private

Definition at line 223 of file TTGTIKProcessor.h.

◆ M2

Double_t art::crib::TTGTIKProcessor::M2

private

Definition at line 224 of file TTGTIKProcessor.h.

◆ M3

Double_t art::crib::TTGTIKProcessor::M3

private

Definition at line 226 of file TTGTIKProcessor.h.

◆ M3_default

Double_t art::crib::TTGTIKProcessor::M3_default

private

Definition at line 225 of file TTGTIKProcessor.h.

◆ M4

Double_t art::crib::TTGTIKProcessor::M4

private

Definition at line 227 of file TTGTIKProcessor.h.

◆ srim

TSrim* art::crib::TTGTIKProcessor::srim

private

! TSrim object to calculate energy loss

Definition at line 214 of file TTGTIKProcessor.h.

The documentation for this class was generated from the following files:

src-crib/reconst/TTGTIKProcessor.h
src-crib/reconst/TTGTIKProcessor.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Example Steering File

Kinematics Calculation

Constructor & Destructor Documentation

◆ TTGTIKProcessor() [1/2]

◆ ~TTGTIKProcessor()

◆ TTGTIKProcessor() [2/2]

Member Function Documentation

◆ bisection()

◆ ClassDefOverride()

◆ GetCMAngle()

◆ GetCustomExcitedEnergy()

◆ GetEcm_classic_kinematics()

◆ GetEcm_kinematics()

◆ GetEcmFromBeam()

◆ GetEcmFromDetectParticle()

◆ GetELabALabPair()

◆ GetReactionPosition()

◆ Init()

◆ newton()

◆ operator=()

◆ Process()

◆ TargetFunction()

Member Data Documentation

◆ fDetectorParameterName

◆ fDetectorPrm

◆ fDoCenterPos

◆ fDoCustom

◆ fExcitedEnergy

◆ fInData

◆ fInitialBeamEnergy

◆ fInputColName

◆ fInputTrackColName

◆ fInTrackData

◆ fOutData

◆ fOutputColName

◆ fParticleAArray

◆ fParticleZArray

◆ fPressure

◆ fTargetName

◆ fTargetParameterName

◆ fTargetPrm

◆ fTemperature

◆ kEpsilon

◆ kInitialMax

◆ kInitialMin

◆ kMaxIteration

◆ M1

◆ M2

◆ M3

◆ M3_default

◆ M4

◆ srim