Page 1 of 1

option OP,RAW

PostPosted: Mon May 09, 2011 5:12 pm
by rfgarciar
Dear Gosia Users

I want to use the option OP,RAW, I not sure how to add this options when we have more than one experiment. Let us suppose we have two experiments (two different scattering angles) and 11 gamma detectors.
If I have only one experiment with no clusters, I understand, the input has this form:

OP,YIEL
(...)
4 !NTAP
0,0
0,0
0,0
0,0
OP,RAW
1
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
4.94122,3.20128,-2.92488,0.765712,0,0,0,0
0
0
OP,INTG
(...)

If I have two experiments (only different scattering angle), with no clusters, with the same gamma detectors, How can I do for adding the second experiment?

Best regards,
Ronald

Re: option OP,RAW

PostPosted: Tue May 10, 2011 2:42 pm
by hayes
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

Re: option OP,RAW

PostPosted: Wed May 11, 2011 12:38 pm
by rfgarciar
Hello Adam,

Thanks a lot, It is all I need!

Best regards,
Ronald

Re: option OP,RAW

PostPosted: Thu May 12, 2011 8:10 pm
by rfgarciar
Hi Adam,
And how the input for the experimental yields is? (input.yld)
The file .yld has two yields in this case? each one containing the sum over five Gecrystals for each experiment?. Seeing your example, it is more or less this form (?):
1 1 72 178. 976.45 1 1.
2 1 1634 13 (dummy values)
2 1 72 178. 976.45 1 1.
2 1 486 23
Or must I leave the entry as it was without using OP, RAW?, as in your example; ten experimental yields (for each Ge crystal), and in the last yield of each experiment to put the sum over the five Ge crystals?

Thanks a lot,
Ronald

Re: option OP,RAW

PostPosted: Fri May 20, 2011 10:49 pm
by hayes
rfgarciar wrote:1 1 72 178. 976.45 1 1.
2 1 1634 13 (dummy values)
2 1 72 178. 976.45 1 1.
2 1 486 23


Hi Ronald,

The format you gave above is correct for two experiments, each with one cluster. One purpose of OP,RAW is to sum the yields in a set of Ge detectors, so you give the summed data for each cluster. In my examples, I have used values of the parameters for the efficiency curve to give a constant efficiency as a function of energy. OP,RAW is meant to accept the different efficiency curves for each Ge in the cluster, but I have not used it that way.

Attached is another example with two experiments and one cluster in each. The clusters are different for each experiment, but you will probably re-use clusters. In that case, you should repeat the entire set:

Code: Select all
1 1 72  178.  976.45 1  1.         
2   1   1634   13


for each cluster in each experiment, so the whole file will be ordered as

Code: Select all
1 1 72  178.  976.45 1  1.   Expt 1, cluster 1
2   1   1634   13             
1 1 72  178.  976.45 1  1.   Expt 1, cluster 2   
2   1   1500 100 

etc.

1 1 72  178.  976.45 1  1.   Expt 2, cluster 1
2   1   1650 100
1 1 72  178.  976.45 1  1.   Expt 2, cluster 2
2   1   1480 100

etc.


Gosia matches the data with the experiment and detector numbers in the order they appear in the gosia input file.


Edit: When you run the OP,CORR step with NTAP=3, you can see if the right data are being assigned to each detector by looking at the .out file. They will be listed by experiment and detector number.


Adam