DWBA and Monte Carlo Simulation

Reaction:

Beam Energy	MeV/u
Beam	J^π	Ex:	MeV
Target
Light recoil
Number of events

Detector:

HELIOS
SOLARIS
ISS

B-field		T (minus sign = field point to upstream)
Array Pos.		mm (negative for upstream)
Recoil Pos.		mm (negative for upstream)

DWBA and E_x List:

For 2-nucleon transfer, Orbital take the form NL=X, where N is number of node, X is momentum number. n and L are related by Σ_i (2n_i+l_i) = 2N + X + 2n + l, where n_i and l_i are the principle number and orbital angular momentum of the each transfered nucleon, and n and l are the internal quanta of the 2-nucleon. e.g. (t,p) reaction to 0f orbtial, the left-hand side would be n_i = 0 and l_i = 3 and the sum is 3+3 = 6 = 2N + X + 2n+l. Assume n = l = 0, we have 6 = 2N+L. Thus, 3L=0, 2L=2,1L=4, 0L=6.

TODO: guess the orbital
Positive parity Negative parity Unknown parity
Isotope: Max Ex: MeV

E_x [MeV]	J^π	Orbital

Cal. DWBA

Incoming Channel
Outgoing Channel

Only DWBA and Don't Sim. Angle range (for only DWBA)

Min [deg]:

Max [deg]:

Plot config:

E vs Z	Ex (cal.)	ThetaCM
ThetaCM vs Z	Recoil X vs Y	Recoil-R vs ThetaCM
Recoil R vs Z	Time diff vs Z	Array X vs Y

Array Hit
Loop = 1
Loop = 2
Loop <= 2
ThetaCM > 10 deg
rhoRecoil > 10 mm
rhoRecoil < 40 mm

Advanced control

Download Sample files:	Reaction File
	DetectorGeo (SOLARIS) File
	DetectorGeo (HELIOS) File
	Ex File
	DWBA File
	Plot Config File

File name can be customized.
For kinematic simulation, only the reactionConfig.txt, detectorGeo.txt, and Ex.txt are needed.
For DWBA calculation, only the DWBA file is needed.
If reactionConfig.txt, detectorGeo.txt, and DWBA file are presented, will do DWBA and use the result for simulation.
When the DWBA file is presented, the kinematic simulation will use the DWBA result and also the excitation energy. i.e. the user provide Ex.txt will not be used.
User can use a customs in File for DWBA calculation. Once the in File exist, it ignores the DWBA file.
To change DWBA angular range, download the in file, edit it. But becareful, DWBA for kinematic simulation must be 0 -180 deg.

The source code for calculation can be found in Here

Reaction File
DetectorGeo File
Ex File
DWBA File ^
*.in File ^
Plot Config File #
^ can be alone # can be omitted