by hayes » Tue May 10, 2011 2:42 pm
Hi Ronald,
First, I think there is an extra line containing "0" here:
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
0
0
OP,INTG
This might not matter for one cluster, but you should make sure to put the terminating "0" line in only after defining clusters for all experiments. (Maybe I am wrong about the extra "0". Check my thinking.)
Attached is a tar file containing an example of an integration using the same cluster for two experiments. This can be expanded to include more than one cluster for each experiment by "pattern-matching" with the example input.
Things to remember:
Change the paths in the OP,BRIC lines in the attached example input.
Run OP,GDET with the first entry flagged with a negative sign, so that the file number 8 for OP,RAW is created. I call this *.abs.
Here is an annotated input example:
OP,YIEL
0
-1,0
1.0
5,5 Five Ge crystals for each experiment
1, 1, 1, 1, 1 These are all detector type "1" in the "gdt" file (file 9)
142.62, 121.72, 129.93, 129.93, 148.28 Polar angles of the 5 Ge crystals
108.0, 108.0, 85.16, 130.84, 72.0 Azimuthal angles "
1, 1, 1, 1, 1 Similar for experiment 2
142.62, 121.72, 129.93, 129.93, 148.28
108.0, 108.0, 85.16, 130.84, 72.0
5,4 The normalizing transition
1 Only 1 detector defined here for experiment 1 (one cluster)
1000.0 The upper limit on measurements. (Change this to an appropriate value for OP,MINI.)
1.0 The normalization constant for cluster 1, experiment 1.
1 Similar for experiment 2 with one cluster.
1000.0
1.0
0 !NTAP
0,0
0,0
0,0
0,0
OP,RAW
1 Experiment number 1
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0 Efficiency coefficients (dummy values here to produce a flat efficiency curve
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0 for efficiency-corrected data).
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0 One line for each Ge crystal in the cluster.
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
1 One cluster defined.
5 The cluster contains all 5 crystals
1,2,3,4,5 The logical Ge detector numbers in the crystal.
2 Similar definitions for EXPERIMENT 2
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0
1
5
1,2,3,4,5
0 One line containing zero to terminate OP,RAW
OP,INTG
The full input with its associated files produces an output like:
INTEGRATED YIELDS
EXPERIMENT 1 DETECTOR 1
ENERGY RANGE 967.700--- 985.200 MEV SCATTERING ANGLE RANGE 20.000--- 41.670 DEG
NI NF II IF YIELD NORMALIZED YIELD
10 9 18.0 16.0 0.16811E-01 0.82252E-03
24 8 15.0 14.0 0.21758E-03 0.10646E-04
...
INTEGRATED YIELDS
EXPERIMENT 2 DETECTOR 1
ENERGY RANGE 967.700--- 985.200 MEV SCATTERING ANGLE RANGE 20.000--- 41.670 DEG
NI NF II IF YIELD NORMALIZED YIELD
10 9 18.0 16.0 0.20913E-01 0.87559E-03
24 8 15.0 14.0 0.19681E-03 0.82401E-05
...
I hope this helps.
Adam
- Attachments
-
op_raw_example_3.tar
- Extract the files in a new directory using
tar -xvf op_raw_example_3.tar
- (100 KiB) Downloaded 378 times