It is currently Thu Mar 28, 2024 8:11 pm


option OP,RAW

For questions about writing Gosia inputs

option OP,RAW

Postby rfgarciar » Mon May 09, 2011 5:12 pm

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
rfgarciar
 
Posts: 6
Joined: Mon May 09, 2011 4:48 pm

Re: option OP,RAW

Postby 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 345 times
hayes
 
Posts: 45
Joined: Tue Feb 08, 2011 2:13 pm
Location: Rochester, NY, USA

Re: option OP,RAW

Postby rfgarciar » Wed May 11, 2011 12:38 pm

Hello Adam,

Thanks a lot, It is all I need!

Best regards,
Ronald
rfgarciar
 
Posts: 6
Joined: Mon May 09, 2011 4:48 pm

Re: option OP,RAW

Postby rfgarciar » Thu May 12, 2011 8:10 pm

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
rfgarciar
 
Posts: 6
Joined: Mon May 09, 2011 4:48 pm

Re: option OP,RAW

Postby hayes » Fri May 20, 2011 10:49 pm

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
Attachments
op_raw_example_2.tar
(720 KiB) Downloaded 284 times
hayes
 
Posts: 45
Joined: Tue Feb 08, 2011 2:13 pm
Location: Rochester, NY, USA


Return to General Usage Questions

Who is online

Users browsing this forum: No registered users and 0 guests

cron