General remarks to the procedure

Note: The input file formats and some general issues are explained in chuck3_manual.pdf and angcor_manual.pdf.


-----------------

0.) Create the input files for CHUCK3

There are some examples in the input file folders.

The Harakeh Dieperink dipole form factors which are needed for isoscalar dipole excitations can be calculated using the program formf in the corresponding directory.
There is an example input file called "input" in the same directory. The program can be used with an input file or via typing the values in the command line.
It basically need the mass and optical potential parameters.


-----------------

1.) Run Chuck3 with the desired input-file
./chuck < input.com

This command creates a "fort.2" file in the chuck directory which contains the relevant output for ANGCOR.
ANGCOR needs the "fort.2" file in the directory where it is started.
One can permantly fix this by placing a symbolic link in this directory.


-----------------

2.) Create input files for ANGCOR
You can use and adjust the c program (make_input) provided in the angcor/input/ directory to create input files for various angles.
These input files are needed to run angcor in its directory.


-----------------

3.) Run ANGCOR for each input file

./angcor < input/example_010.com > output/example_010.lis

for each input file.

I usually use it like this:

for ball in input_ca_*.com; do ../angcor < $ball > ../output/${ball%com}lis; done


-----------------

4.) Extract the angular correlation information from the ANGCOR output files

There is a script in the ANGCOR output directory which reads the *.lis output files and extracts the necessary values from the *.lis files.

The script

./make_final.sh

uses awk sort.awk and generates a 4-column "final.dat" file which contains the relevant information of angles theta_particle,theta_gamma,phi_gamma and W(theta_particle,theta_gamma,phi_gamma).


-----------------

5.) Averaging over the acceptance of the spectrometer

The last step is to use the "final.dat" file to average over the spectrometer acceptance that means the angles that can be detected in the experiment.

W(theta_particle,theta_gamma,phi_gamma)  --->  W(theta_gamma)

I have an old c program but it's very specialized, so it is probably not useful for you. I just put it that you get an idea of how it can be done.
./averaging final.dat


-----------------

VD on 03/31/15

