Skip to content
Snippets Groups Projects
Commit d6e510ed authored by deserevi's avatar deserevi
Browse files

propagate from branch 'NPTool.pro' (head bfd52c310dc09fb00c97bc81e1d8758bfba81548) 

            to branch 'NPTool.dev' (head ab0e8c73675c98f6ac3dc2d44ade51ae18a36095)
parents 70b828f4 a594fc74
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1280 additions and 20 deletions
INSTALL
+ 6
− 1
View file @ d6e510ed
...@@ -54,7 +54,9 @@ III) INSTALLATION ...@@ -54,7 +54,9 @@ III) INSTALLATION
b) First, you need to define some environment variables. You can b) First, you need to define some environment variables. You can
source the NPEnvironment.bash or NPEnvironment.tchrc file source the NPEnvironment.bash or NPEnvironment.tchrc file
depending on your shell. depending on your shell or include the relevant lines in your
.bashrc or .tcshrc file. Don't forget to customize these file
with the relevant path to your NPTool installation.
This step is mandatory since the Makefiles used to build the package This step is mandatory since the Makefiles used to build the package
depend on some of these environment variables. depend on some of these environment variables.
...@@ -76,3 +78,6 @@ III) INSTALLATION ...@@ -76,3 +78,6 @@ III) INSTALLATION
cd Template cd Template
make make
This step is not mandatory for the G4 simulation to work. This step is not mandatory for the G4 simulation to work.
If you want to use the ROOT macros in the $NPAnalysis/macros
directory you should execute in a ROOT terminal the commands
you can find in rootlogon.C
Inputs/CrossSection/132Sndp_10A_MeV_2f7_ZR_FRC.lis 0 → 100644
+ 180
− 0
View file @ d6e510ed
0.00 1.3462E+00
1.00 1.3514E+00
2.00 1.3667E+00
3.00 1.3912E+00
4.00 1.4235E+00
5.00 1.4617E+00
6.00 1.5036E+00
7.00 1.5465E+00
8.00 1.5880E+00
9.00 1.6253E+00
10.00 1.6560E+00
11.00 1.6779E+00
12.00 1.6892E+00
13.00 1.6883E+00
14.00 1.6745E+00
15.00 1.6475E+00
16.00 1.6074E+00
17.00 1.5550E+00
18.00 1.4916E+00
19.00 1.4189E+00
20.00 1.3388E+00
21.00 1.2537E+00
22.00 1.1659E+00
23.00 1.0778E+00
24.00 9.9188E-01
25.00 9.1010E-01
26.00 8.3436E-01
27.00 7.6610E-01
28.00 7.0638E-01
29.00 6.5580E-01
30.00 6.1451E-01
31.00 5.8225E-01
32.00 5.5836E-01
33.00 5.4186E-01
34.00 5.3155E-01
35.00 5.2602E-01
36.00 5.2380E-01
37.00 5.2344E-01
38.00 5.2353E-01
39.00 5.2283E-01
40.00 5.2029E-01
41.00 5.1508E-01
42.00 5.0664E-01
43.00 4.9465E-01
44.00 4.7904E-01
45.00 4.5999E-01
46.00 4.3784E-01
47.00 4.1310E-01
48.00 3.8641E-01
49.00 3.5846E-01
50.00 3.2997E-01
51.00 3.0165E-01
52.00 2.7416E-01
53.00 2.4807E-01
54.00 2.2388E-01
55.00 2.0195E-01
56.00 1.8253E-01
57.00 1.6576E-01
58.00 1.5166E-01
59.00 1.4015E-01
60.00 1.3106E-01
61.00 1.2416E-01
62.00 1.1916E-01
63.00 1.1574E-01
64.00 1.1356E-01
65.00 1.1229E-01
66.00 1.1162E-01
67.00 1.1124E-01
68.00 1.1090E-01
69.00 1.1037E-01
70.00 1.0950E-01
71.00 1.0815E-01
72.00 1.0624E-01
73.00 1.0374E-01
74.00 1.0063E-01
75.00 9.6948E-02
76.00 9.2758E-02
77.00 8.8137E-02
78.00 8.3179E-02
79.00 7.7987E-02
80.00 7.2667E-02
81.00 6.7326E-02
82.00 6.2064E-02
83.00 5.6972E-02
84.00 5.2130E-02
85.00 4.7608E-02
86.00 4.3458E-02
87.00 3.9720E-02
88.00 3.6417E-02
89.00 3.3560E-02
90.00 3.1146E-02
91.00 2.9158E-02
92.00 2.7571E-02
93.00 2.6351E-02
94.00 2.5456E-02
95.00 2.4840E-02
96.00 2.4455E-02
97.00 2.4249E-02
98.00 2.4174E-02
99.00 2.4181E-02
100.00 2.4225E-02
101.00 2.4267E-02
102.00 2.4271E-02
103.00 2.4208E-02
104.00 2.4056E-02
105.00 2.3799E-02
106.00 2.3426E-02
107.00 2.2933E-02
108.00 2.2323E-02
109.00 2.1602E-02
110.00 2.0782E-02
111.00 1.9876E-02
112.00 1.8903E-02
113.00 1.7882E-02
114.00 1.6834E-02
115.00 1.5778E-02
116.00 1.4735E-02
117.00 1.3723E-02
118.00 1.2758E-02
119.00 1.1856E-02
120.00 1.1029E-02
121.00 1.0284E-02
122.00 9.6287E-03
123.00 9.0660E-03
124.00 8.5967E-03
125.00 8.2187E-03
126.00 7.9277E-03
127.00 7.7175E-03
128.00 7.5802E-03
129.00 7.5067E-03
130.00 7.4870E-03
131.00 7.5108E-03
132.00 7.5676E-03
133.00 7.6472E-03
134.00 7.7401E-03
135.00 7.8375E-03
136.00 7.9313E-03
137.00 8.0149E-03
138.00 8.0825E-03
139.00 8.1299E-03
140.00 8.1539E-03
141.00 8.1527E-03
142.00 8.1255E-03
143.00 8.0728E-03
144.00 7.9958E-03
145.00 7.8968E-03
146.00 7.7786E-03
147.00 7.6444E-03
148.00 7.4980E-03
149.00 7.3433E-03
150.00 7.1843E-03
151.00 7.0249E-03
152.00 6.8691E-03
153.00 6.7204E-03
154.00 6.5822E-03
155.00 6.4574E-03
156.00 6.3484E-03
157.00 6.2573E-03
158.00 6.1857E-03
159.00 6.1345E-03
160.00 6.1042E-03
161.00 6.0948E-03
162.00 6.1056E-03
163.00 6.1357E-03
164.00 6.1836E-03
165.00 6.2475E-03
166.00 6.3252E-03
167.00 6.4144E-03
168.00 6.5123E-03
169.00 6.6163E-03
170.00 6.7236E-03
171.00 6.8313E-03
172.00 6.9367E-03
173.00 7.0369E-03
174.00 7.1296E-03
175.00 7.2123E-03
176.00 7.2829E-03
177.00 7.3397E-03
178.00 7.3814E-03
179.00 7.4068E-03
Inputs/CrossSection/132Sndp_10A_MeV_3p3_ZR_FRC.lis 0 → 100644
+ 180
− 0
View file @ d6e510ed
0.00 1.0371E+00
1.00 1.0303E+00
2.00 1.0103E+00
3.00 9.7770E-01
4.00 9.3402E-01
5.00 8.8111E-01
6.00 8.2140E-01
7.00 7.5780E-01
8.00 6.9354E-01
9.00 6.3199E-01
10.00 5.7645E-01
11.00 5.2991E-01
12.00 4.9476E-01
13.00 4.7268E-01
14.00 4.6444E-01
15.00 4.6984E-01
16.00 4.8774E-01
17.00 5.1611E-01
18.00 5.5220E-01
19.00 5.9279E-01
20.00 6.3440E-01
21.00 6.7364E-01
22.00 7.0747E-01
23.00 7.3337E-01
24.00 7.4958E-01
25.00 7.5507E-01
26.00 7.4961E-01
27.00 7.3360E-01
28.00 7.0807E-01
29.00 6.7443E-01
30.00 6.3445E-01
31.00 5.9000E-01
32.00 5.4304E-01
33.00 4.9542E-01
34.00 4.4879E-01
35.00 4.0454E-01
36.00 3.6378E-01
37.00 3.2725E-01
38.00 2.9542E-01
39.00 2.6851E-01
40.00 2.4648E-01
41.00 2.2916E-01
42.00 2.1623E-01
43.00 2.0724E-01
44.00 2.0168E-01
45.00 1.9899E-01
46.00 1.9856E-01
47.00 1.9981E-01
48.00 2.0215E-01
49.00 2.0506E-01
50.00 2.0808E-01
51.00 2.1079E-01
52.00 2.1283E-01
53.00 2.1391E-01
54.00 2.1378E-01
55.00 2.1227E-01
56.00 2.0927E-01
57.00 2.0476E-01
58.00 1.9876E-01
59.00 1.9139E-01
60.00 1.8280E-01
61.00 1.7318E-01
62.00 1.6276E-01
63.00 1.5176E-01
64.00 1.4042E-01
65.00 1.2899E-01
66.00 1.1771E-01
67.00 1.0680E-01
68.00 9.6462E-02
69.00 8.6886E-02
70.00 7.8213E-02
71.00 7.0542E-02
72.00 6.3934E-02
73.00 5.8409E-02
74.00 5.3950E-02
75.00 5.0511E-02
76.00 4.8019E-02
77.00 4.6380E-02
78.00 4.5481E-02
79.00 4.5196E-02
80.00 4.5387E-02
81.00 4.5913E-02
82.00 4.6637E-02
83.00 4.7429E-02
84.00 4.8173E-02
85.00 4.8774E-02
86.00 4.9153E-02
87.00 4.9254E-02
88.00 4.9035E-02
89.00 4.8474E-02
90.00 4.7562E-02
91.00 4.6306E-02
92.00 4.4724E-02
93.00 4.2850E-02
94.00 4.0726E-02
95.00 3.8402E-02
96.00 3.5934E-02
97.00 3.3382E-02
98.00 3.0803E-02
99.00 2.8253E-02
100.00 2.5780E-02
101.00 2.3430E-02
102.00 2.1240E-02
103.00 1.9241E-02
104.00 1.7454E-02
105.00 1.5893E-02
106.00 1.4565E-02
107.00 1.3469E-02
108.00 1.2595E-02
109.00 1.1931E-02
110.00 1.1458E-02
111.00 1.1157E-02
112.00 1.1003E-02
113.00 1.0973E-02
114.00 1.1040E-02
115.00 1.1179E-02
116.00 1.1365E-02
117.00 1.1574E-02
118.00 1.1785E-02
119.00 1.1977E-02
120.00 1.2137E-02
121.00 1.2253E-02
122.00 1.2316E-02
123.00 1.2321E-02
124.00 1.2265E-02
125.00 1.2147E-02
126.00 1.1968E-02
127.00 1.1729E-02
128.00 1.1436E-02
129.00 1.1094E-02
130.00 1.0711E-02
131.00 1.0295E-02
132.00 9.8568E-03
133.00 9.4050E-03
134.00 8.9484E-03
135.00 8.4946E-03
136.00 8.0503E-03
137.00 7.6215E-03
138.00 7.2132E-03
139.00 6.8301E-03
140.00 6.4764E-03
141.00 6.1557E-03
142.00 5.8708E-03
143.00 5.6235E-03
144.00 5.4145E-03
145.00 5.2436E-03
146.00 5.1096E-03
147.00 5.0109E-03
148.00 4.9455E-03
149.00 4.9113E-03
150.00 4.9062E-03
151.00 4.9279E-03
152.00 4.9741E-03
153.00 5.0421E-03
154.00 5.1293E-03
155.00 5.2327E-03
156.00 5.3496E-03
157.00 5.4772E-03
158.00 5.6132E-03
159.00 5.7552E-03
160.00 5.9014E-03
161.00 6.0497E-03
162.00 6.1984E-03
163.00 6.3455E-03
164.00 6.4894E-03
165.00 6.6287E-03
166.00 6.7622E-03
167.00 6.8891E-03
168.00 7.0088E-03
169.00 7.1210E-03
170.00 7.2252E-03
171.00 7.3209E-03
172.00 7.4075E-03
173.00 7.4842E-03
174.00 7.5504E-03
175.00 7.6056E-03
176.00 7.6498E-03
177.00 7.6833E-03
178.00 7.7065E-03
179.00 7.7201E-03
Inputs/CrossSection/132Sndp_5A_MeV_2f7_ZR_FRC.lis 0 → 100644
+ 180
− 0
View file @ d6e510ed
0.00 7.9396E-02
1.00 7.9415E-02
2.00 7.9475E-02
3.00 7.9581E-02
4.00 7.9747E-02
5.00 7.9991E-02
6.00 8.0337E-02
7.00 8.0821E-02
8.00 8.1489E-02
9.00 8.2399E-02
10.00 8.3625E-02
11.00 8.5256E-02
12.00 8.7401E-02
13.00 9.0184E-02
14.00 9.3746E-02
15.00 9.8242E-02
16.00 1.0384E-01
17.00 1.1070E-01
18.00 1.1899E-01
19.00 1.2887E-01
20.00 1.4046E-01
21.00 1.5387E-01
22.00 1.6916E-01
23.00 1.8635E-01
24.00 2.0541E-01
25.00 2.2626E-01
26.00 2.4875E-01
27.00 2.7269E-01
28.00 2.9782E-01
29.00 3.2387E-01
30.00 3.5050E-01
31.00 3.7737E-01
32.00 4.0410E-01
33.00 4.3032E-01
34.00 4.5568E-01
35.00 4.7983E-01
36.00 5.0246E-01
37.00 5.2329E-01
38.00 5.4208E-01
39.00 5.5866E-01
40.00 5.7287E-01
41.00 5.8463E-01
42.00 5.9391E-01
43.00 6.0072E-01
44.00 6.0510E-01
45.00 6.0716E-01
46.00 6.0703E-01
47.00 6.0486E-01
48.00 6.0083E-01
49.00 5.9514E-01
50.00 5.8800E-01
51.00 5.7962E-01
52.00 5.7022E-01
53.00 5.6000E-01
54.00 5.4917E-01
55.00 5.3792E-01
56.00 5.2642E-01
57.00 5.1484E-01
58.00 5.0333E-01
59.00 4.9201E-01
60.00 4.8101E-01
61.00 4.7041E-01
62.00 4.6029E-01
63.00 4.5073E-01
64.00 4.4176E-01
65.00 4.3343E-01
66.00 4.2574E-01
67.00 4.1872E-01
68.00 4.1235E-01
69.00 4.0663E-01
70.00 4.0153E-01
71.00 3.9701E-01
72.00 3.9304E-01
73.00 3.8958E-01
74.00 3.8658E-01
75.00 3.8398E-01
76.00 3.8173E-01
77.00 3.7977E-01
78.00 3.7805E-01
79.00 3.7650E-01
80.00 3.7507E-01
81.00 3.7370E-01
82.00 3.7233E-01
83.00 3.7092E-01
84.00 3.6941E-01
85.00 3.6776E-01
86.00 3.6592E-01
87.00 3.6387E-01
88.00 3.6156E-01
89.00 3.5898E-01
90.00 3.5610E-01
91.00 3.5290E-01
92.00 3.4938E-01
93.00 3.4553E-01
94.00 3.4135E-01
95.00 3.3685E-01
96.00 3.3204E-01
97.00 3.2694E-01
98.00 3.2157E-01
99.00 3.1594E-01
100.00 3.1010E-01
101.00 3.0406E-01
102.00 2.9786E-01
103.00 2.9154E-01
104.00 2.8513E-01
105.00 2.7867E-01
106.00 2.7219E-01
107.00 2.6572E-01
108.00 2.5932E-01
109.00 2.5299E-01
110.00 2.4679E-01
111.00 2.4073E-01
112.00 2.3485E-01
113.00 2.2916E-01
114.00 2.2370E-01
115.00 2.1848E-01
116.00 2.1351E-01
117.00 2.0880E-01
118.00 2.0437E-01
119.00 2.0022E-01
120.00 1.9636E-01
121.00 1.9278E-01
122.00 1.8948E-01
123.00 1.8645E-01
124.00 1.8369E-01
125.00 1.8119E-01
126.00 1.7894E-01
127.00 1.7691E-01
128.00 1.7511E-01
129.00 1.7350E-01
130.00 1.7209E-01
131.00 1.7084E-01
132.00 1.6974E-01
133.00 1.6877E-01
134.00 1.6793E-01
135.00 1.6718E-01
136.00 1.6651E-01
137.00 1.6592E-01
138.00 1.6538E-01
139.00 1.6487E-01
140.00 1.6439E-01
141.00 1.6393E-01
142.00 1.6347E-01
143.00 1.6300E-01
144.00 1.6253E-01
145.00 1.6203E-01
146.00 1.6150E-01
147.00 1.6095E-01
148.00 1.6036E-01
149.00 1.5973E-01
150.00 1.5907E-01
151.00 1.5837E-01
152.00 1.5764E-01
153.00 1.5687E-01
154.00 1.5607E-01
155.00 1.5524E-01
156.00 1.5438E-01
157.00 1.5350E-01
158.00 1.5260E-01
159.00 1.5170E-01
160.00 1.5078E-01
161.00 1.4987E-01
162.00 1.4896E-01
163.00 1.4806E-01
164.00 1.4717E-01
165.00 1.4631E-01
166.00 1.4547E-01
167.00 1.4467E-01
168.00 1.4390E-01
169.00 1.4318E-01
170.00 1.4250E-01
171.00 1.4187E-01
172.00 1.4130E-01
173.00 1.4078E-01
174.00 1.4033E-01
175.00 1.3994E-01
176.00 1.3962E-01
177.00 1.3937E-01
178.00 1.3919E-01
179.00 1.3908E-01
Inputs/CrossSection/132Sndp_5A_MeV_3p3_ZR_FRC.lis 0 → 100644
+ 180
− 0
View file @ d6e510ed
0.00 2.7754E-01
1.00 2.7835E-01
2.00 2.8077E-01
3.00 2.8484E-01
4.00 2.9057E-01
5.00 2.9801E-01
6.00 3.0720E-01
7.00 3.1816E-01
8.00 3.3091E-01
9.00 3.4542E-01
10.00 3.6166E-01
11.00 3.7953E-01
12.00 3.9889E-01
13.00 4.1958E-01
14.00 4.4134E-01
15.00 4.6393E-01
16.00 4.8701E-01
17.00 5.1025E-01
18.00 5.3330E-01
19.00 5.5580E-01
20.00 5.7741E-01
21.00 5.9782E-01
22.00 6.1680E-01
23.00 6.3419E-01
24.00 6.4991E-01
25.00 6.6398E-01
26.00 6.7651E-01
27.00 6.8769E-01
28.00 6.9777E-01
29.00 7.0702E-01
30.00 7.1572E-01
31.00 7.2416E-01
32.00 7.3257E-01
33.00 7.4116E-01
34.00 7.5012E-01
35.00 7.5956E-01
36.00 7.6959E-01
37.00 7.8028E-01
38.00 7.9167E-01
39.00 8.0374E-01
40.00 8.1648E-01
41.00 8.2982E-01
42.00 8.4365E-01
43.00 8.5785E-01
44.00 8.7226E-01
45.00 8.8670E-01
46.00 9.0100E-01
47.00 9.1499E-01
48.00 9.2850E-01
49.00 9.4139E-01
50.00 9.5353E-01
51.00 9.6480E-01
52.00 9.7510E-01
53.00 9.8435E-01
54.00 9.9247E-01
55.00 9.9939E-01
56.00 1.0050E+00
57.00 1.0094E+00
58.00 1.0123E+00
59.00 1.0139E+00
60.00 1.0140E+00
61.00 1.0127E+00
62.00 1.0100E+00
63.00 1.0059E+00
64.00 1.0004E+00
65.00 9.9367E-01
66.00 9.8565E-01
67.00 9.7645E-01
68.00 9.6613E-01
69.00 9.5476E-01
70.00 9.4242E-01
71.00 9.2918E-01
72.00 9.1512E-01
73.00 9.0030E-01
74.00 8.8481E-01
75.00 8.6873E-01
76.00 8.5213E-01
77.00 8.3510E-01
78.00 8.1773E-01
79.00 8.0011E-01
80.00 7.8231E-01
81.00 7.6444E-01
82.00 7.4658E-01
83.00 7.2880E-01
84.00 7.1120E-01
85.00 6.9384E-01
86.00 6.7680E-01
87.00 6.6014E-01
88.00 6.4392E-01
89.00 6.2818E-01
90.00 6.1297E-01
91.00 5.9833E-01
92.00 5.8428E-01
93.00 5.7085E-01
94.00 5.5804E-01
95.00 5.4586E-01
96.00 5.3433E-01
97.00 5.2342E-01
98.00 5.1314E-01
99.00 5.0347E-01
100.00 4.9439E-01
101.00 4.8588E-01
102.00 4.7792E-01
103.00 4.7046E-01
104.00 4.6350E-01
105.00 4.5698E-01
106.00 4.5088E-01
107.00 4.4516E-01
108.00 4.3979E-01
109.00 4.3473E-01
110.00 4.2994E-01
111.00 4.2539E-01
112.00 4.2105E-01
113.00 4.1688E-01
114.00 4.1285E-01
115.00 4.0894E-01
116.00 4.0511E-01
117.00 4.0134E-01
118.00 3.9762E-01
119.00 3.9391E-01
120.00 3.9020E-01
121.00 3.8648E-01
122.00 3.8273E-01
123.00 3.7895E-01
124.00 3.7512E-01
125.00 3.7124E-01
126.00 3.6732E-01
127.00 3.6334E-01
128.00 3.5931E-01
129.00 3.5523E-01
130.00 3.5111E-01
131.00 3.4695E-01
132.00 3.4276E-01
133.00 3.3854E-01
134.00 3.3431E-01
135.00 3.3007E-01
136.00 3.2584E-01
137.00 3.2161E-01
138.00 3.1741E-01
139.00 3.1324E-01
140.00 3.0911E-01
141.00 3.0503E-01
142.00 3.0100E-01
143.00 2.9705E-01
144.00 2.9317E-01
145.00 2.8937E-01
146.00 2.8566E-01
147.00 2.8205E-01
148.00 2.7854E-01
149.00 2.7513E-01
150.00 2.7184E-01
151.00 2.6865E-01
152.00 2.6559E-01
153.00 2.6264E-01
154.00 2.5982E-01
155.00 2.5711E-01
156.00 2.5453E-01
157.00 2.5207E-01
158.00 2.4974E-01
159.00 2.4753E-01
160.00 2.4544E-01
161.00 2.4347E-01
162.00 2.4161E-01
163.00 2.3988E-01
164.00 2.3826E-01
165.00 2.3675E-01
166.00 2.3536E-01
167.00 2.3407E-01
168.00 2.3289E-01
169.00 2.3181E-01
170.00 2.3083E-01
171.00 2.2996E-01
172.00 2.2918E-01
173.00 2.2850E-01
174.00 2.2791E-01
175.00 2.2742E-01
176.00 2.2702E-01
177.00 2.2671E-01
178.00 2.2648E-01
179.00 2.2635E-01
Inputs/CrossSection/CS_Ep10MeV_sn134pt_gs_1h9demi.dat 0 → 100644
+ 90
− 0
View file @ d6e510ed
0.00 0.16569E-06
2.00 0.16942E-06
4.00 0.18063E-06
6.00 0.19928E-06
8.00 0.22517E-06
10.00 0.25765E-06
12.00 0.29547E-06
14.00 0.33661E-06
16.00 0.37833E-06
18.00 0.41727E-06
20.00 0.44978E-06
22.00 0.47231E-06
24.00 0.48189E-06
26.00 0.47646E-06
28.00 0.45530E-06
30.00 0.41908E-06
32.00 0.36999E-06
34.00 0.31146E-06
36.00 0.24791E-06
38.00 0.18432E-06
40.00 0.12569E-06
42.00 0.76637E-07
44.00 0.40899E-07
46.00 0.21007E-07
48.00 0.18072E-07
50.00 0.31701E-07
52.00 0.60080E-07
54.00 0.10018E-06
56.00 0.14810E-06
58.00 0.19942E-06
60.00 0.24972E-06
62.00 0.29489E-06
64.00 0.33156E-06
66.00 0.35735E-06
68.00 0.37099E-06
70.00 0.37239E-06
72.00 0.36249E-06
74.00 0.34313E-06
76.00 0.31677E-06
78.00 0.28621E-06
80.00 0.25427E-06
82.00 0.22360E-06
84.00 0.19639E-06
86.00 0.17431E-06
88.00 0.15835E-06
90.00 0.14892E-06
92.00 0.14581E-06
94.00 0.14834E-06
96.00 0.15548E-06
98.00 0.16599E-06
100.00 0.17853E-06
102.00 0.19180E-06
104.00 0.20464E-06
106.00 0.21608E-06
108.00 0.22541E-06
110.00 0.23215E-06
112.00 0.23609E-06
114.00 0.23723E-06
116.00 0.23576E-06
118.00 0.23202E-06
120.00 0.22642E-06
122.00 0.21944E-06
124.00 0.21156E-06
126.00 0.20325E-06
128.00 0.19492E-06
130.00 0.18692E-06
132.00 0.17954E-06
134.00 0.17299E-06
136.00 0.16740E-06
138.00 0.16285E-06
140.00 0.15935E-06
142.00 0.15687E-06
144.00 0.15536E-06
146.00 0.15472E-06
148.00 0.15483E-06
150.00 0.15559E-06
152.00 0.15687E-06
154.00 0.15855E-06
156.00 0.16052E-06
158.00 0.16267E-06
160.00 0.16491E-06
162.00 0.16716E-06
164.00 0.16934E-06
166.00 0.17138E-06
168.00 0.17325E-06
170.00 0.17489E-06
172.00 0.17628E-06
174.00 0.17738E-06
176.00 0.17818E-06
178.00 0.17867E-06
Inputs/CrossSection/CS_Ep15MeV_sn134pt_gs_1h9demi.dat 0 → 100644
+ 90
− 0
View file @ d6e510ed
0.00 0.74986E-06
2.00 0.72457E-06
4.00 0.65371E-06
6.00 0.55132E-06
8.00 0.43767E-06
10.00 0.33534E-06
12.00 0.26506E-06
14.00 0.24199E-06
16.00 0.27332E-06
18.00 0.35716E-06
20.00 0.48312E-06
22.00 0.63395E-06
24.00 0.78829E-06
26.00 0.92376E-06
28.00 0.10203E-05
30.00 0.10631E-05
32.00 0.10449E-05
34.00 0.96677E-06
36.00 0.83837E-06
38.00 0.67611E-06
40.00 0.50059E-06
42.00 0.33336E-06
44.00 0.19355E-06
46.00 0.95102E-07
48.00 0.44996E-07
50.00 0.42738E-07
52.00 0.81134E-07
54.00 0.14805E-06
56.00 0.22881E-06
58.00 0.30865E-06
60.00 0.37494E-06
62.00 0.41876E-06
64.00 0.43567E-06
66.00 0.42567E-06
68.00 0.39261E-06
70.00 0.34296E-06
72.00 0.28456E-06
74.00 0.22530E-06
76.00 0.17196E-06
78.00 0.12955E-06
80.00 0.10083E-06
82.00 0.86394E-07
84.00 0.84918E-07
86.00 0.93674E-07
88.00 0.10914E-06
90.00 0.12757E-06
92.00 0.14556E-06
94.00 0.16034E-06
96.00 0.17009E-06
98.00 0.17389E-06
100.00 0.17174E-06
102.00 0.16432E-06
104.00 0.15283E-06
106.00 0.13871E-06
108.00 0.12345E-06
110.00 0.10841E-06
112.00 0.94699E-07
114.00 0.83117E-07
116.00 0.74118E-07
118.00 0.67843E-07
120.00 0.64166E-07
122.00 0.62763E-07
124.00 0.63175E-07
126.00 0.64884E-07
128.00 0.67365E-07
130.00 0.70138E-07
132.00 0.72798E-07
134.00 0.75033E-07
136.00 0.76634E-07
138.00 0.77491E-07
140.00 0.77577E-07
142.00 0.76937E-07
144.00 0.75666E-07
146.00 0.73887E-07
148.00 0.71735E-07
150.00 0.69340E-07
152.00 0.66816E-07
154.00 0.64254E-07
156.00 0.61722E-07
158.00 0.59262E-07
160.00 0.56902E-07
162.00 0.54656E-07
164.00 0.52535E-07
166.00 0.50553E-07
168.00 0.48728E-07
170.00 0.47089E-07
172.00 0.45670E-07
174.00 0.44511E-07
176.00 0.43650E-07
178.00 0.43119E-07
Inputs/CrossSection/CS_Ep20MeV_sn134pt_gs_1h9demi.dat 0 → 100644
+ 90
− 0
View file @ d6e510ed
0.00 0.32271E-05
2.00 0.31495E-05
4.00 0.29259E-05
6.00 0.25828E-05
8.00 0.21597E-05
10.00 0.17028E-05
12.00 0.12592E-05
14.00 0.87079E-06
16.00 0.56980E-06
18.00 0.37542E-06
20.00 0.29185E-06
22.00 0.30811E-06
24.00 0.39985E-06
26.00 0.53346E-06
28.00 0.67184E-06
30.00 0.78100E-06
32.00 0.83609E-06
34.00 0.82530E-06
36.00 0.75101E-06
38.00 0.62785E-06
40.00 0.47846E-06
42.00 0.32808E-06
44.00 0.19927E-06
46.00 0.10786E-06
48.00 0.60824E-07
50.00 0.56181E-07
52.00 0.84677E-07
54.00 0.13263E-06
56.00 0.18522E-06
58.00 0.22945E-06
60.00 0.25634E-06
62.00 0.26191E-06
64.00 0.24710E-06
66.00 0.21667E-06
68.00 0.17755E-06
70.00 0.13712E-06
72.00 0.10164E-06
74.00 0.75298E-07
76.00 0.59816E-07
78.00 0.54627E-07
80.00 0.57467E-07
82.00 0.65166E-07
84.00 0.74426E-07
86.00 0.82450E-07
88.00 0.87327E-07
90.00 0.88163E-07
92.00 0.84998E-07
94.00 0.78573E-07
96.00 0.70036E-07
98.00 0.60642E-07
100.00 0.51524E-07
102.00 0.43529E-07
104.00 0.37160E-07
106.00 0.32574E-07
108.00 0.29655E-07
110.00 0.28097E-07
112.00 0.27507E-07
114.00 0.27486E-07
116.00 0.27687E-07
118.00 0.27854E-07
120.00 0.27826E-07
122.00 0.27528E-07
124.00 0.26951E-07
126.00 0.26122E-07
128.00 0.25090E-07
130.00 0.23902E-07
132.00 0.22601E-07
134.00 0.21225E-07
136.00 0.19809E-07
138.00 0.18389E-07
140.00 0.17013E-07
142.00 0.15735E-07
144.00 0.14618E-07
146.00 0.13719E-07
148.00 0.13087E-07
150.00 0.12743E-07
152.00 0.12683E-07
154.00 0.12864E-07
156.00 0.13213E-07
158.00 0.13633E-07
160.00 0.14018E-07
162.00 0.14267E-07
164.00 0.14306E-07
166.00 0.14095E-07
168.00 0.13638E-07
170.00 0.12986E-07
172.00 0.12225E-07
174.00 0.11465E-07
176.00 0.10819E-07
178.00 0.10386E-07
Inputs/DetectorConfiguration/gaspardTestSpheric.detector
+ 4
− 3
View file @ d6e510ed
%Fichier de configuration manip E225
%%%%%%%%%%%Target%%%%%%%%%%%%%%%%%%%1 %%%%%%%%%%%Target%%%%%%%%%%%%%%%%%%%1
%Thickness in micrometer %Thickness in micrometer
%Radius in mm %Radius in mm
...@@ -9,14 +8,16 @@ ...@@ -9,14 +8,16 @@
GeneralTarget GeneralTarget
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Target Target
THICKNESS= 10.3 THICKNESS= 1.0
ANGLE= 0 ANGLE= 0
RADIUS= 12 RADIUS= 12
MATERIAL= CD2 MATERIAL= CH2
NBLAYERS= 50 NBLAYERS= 50
X= 0 X= 0
Y= 0 Y= 0
Z= 0 Z= 0
%%%%%%%%%%Detector%%%%%%%%%%%%%%%%%%% %%%%%%%%%%Detector%%%%%%%%%%%%%%%%%%%
%%Position and R given in mm %%Position and R given in mm
%%Angle given in degree %%Angle given in degree
......
Inputs/EventGenerator/10He.reaction
+ 1
− 1
View file @ d6e510ed
...@@ -19,7 +19,7 @@ TransfertToResonance ...@@ -19,7 +19,7 @@ TransfertToResonance
ResonanceDecayA= 8 ResonanceDecayA= 8
CrossSectionPath= 11Li(d,3He)10He.txt CrossSectionPath= 11Li(d,3He)10He.txt
ShootLight= 1 ShootLight= 1
ShootHeavy= 0 ShootHeavy= 1
ShootDecayProduct= 0 ShootDecayProduct= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Inputs/EventGenerator/132Sndp.reaction
+ 2
− 1
View file @ d6e510ed
...@@ -14,7 +14,8 @@ Transfert ...@@ -14,7 +14,8 @@ Transfert
SigmaY= 0 SigmaY= 0
SigmaThetaX= 0 SigmaThetaX= 0
SigmaPhiY= 0 SigmaPhiY= 0
CrossSectionPath= sn132dp_gs_10AMeV.txt CrossSectionPath= 132Sndp_10A_MeV_2f7_ZR_FRC.lis
% CrossSectionPath= flat.txt
ShootLight= 1 ShootLight= 1
ShootHeavy= 0 ShootHeavy= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
......
Inputs/EventGenerator/134Snpt.reaction 0 → 100644
+ 21
− 0
View file @ d6e510ed
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%% Reaction file for 60Fe(d,p)61Fe reaction %%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%Beam energy given in MeV ; Excitation in MeV ; emmitance in rad
Transfert
Beam= 134Sn
Target= 1H
Light= 3H
Heavy= 132Sn
ExcitationEnergy= 0.0
BeamEnergy= 1340
BeamEnergySpread= 0
SigmaX= 0
SigmaY= 0
SigmaThetaX= 0
SigmaPhiY= 0
CrossSectionPath= CS_Ep10MeV_sn134pt_gs_1h9demi.dat
ShootLight= 1
ShootHeavy= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NPEnvironment.bash → Misc/NPEnvironment.bash
+ 0
− 0
View file @ d6e510ed
File moved
NPEnvironment.tchrc → Misc/NPEnvironment.tchrc
+ 0
− 0
View file @ d6e510ed
File moved
Misc/rootlogon.C 0 → 100644
+ 5
− 0
View file @ d6e510ed
{
gROOT->ProcessLine(".include /YourPathToNPTool/NPTool/NPLib/include");
gSystem->Load("/YourPathToNPTool/NPTool/NPLib/lib/libInitialConditions.so");
gSystem->Load("/YourPathToNPTool/NPTool/NPLib/lib/libInteractionCoordinates.so");
}
NPAnalysis/Gaspard/RunToTreat.txt
+ 3
− 1
View file @ d6e510ed
...@@ -2,4 +2,6 @@ TTreeName ...@@ -2,4 +2,6 @@ TTreeName
SimulatedTree SimulatedTree
RootFileName RootFileName
../../Outputs/Simulation/mySimul.root ../../Outputs/Simulation/mySimul.root
% ../../Outputs/Simulation/sn132dp_10MeVA_T1_B0_E0.root % ../../Outputs/Simulation/134Snpt_1h9_10MeVA_T0_B0_E0_S2mm.root
% ../../Outputs/Simulation/132Sndp_3p3_10MeVA_T0_B1_E0_S05mm.root
% ../../Outputs/Simulation/134Snpt_1h9_10MeVA_T1_B1_E0_S05mm.root
NPAnalysis/Gaspard/include/ObjectManager.hh
+ 11
− 3
View file @ d6e510ed
...@@ -23,7 +23,7 @@ ...@@ -23,7 +23,7 @@
#include <TFile.h> #include <TFile.h>
#include <TLeaf.h> #include <TLeaf.h>
#include <TVector3.h> #include <TVector3.h>
#include <TRandom.h> #include <TRandom3.h>
// NPL // NPL
#include "TMust2Data.h" #include "TMust2Data.h"
...@@ -32,6 +32,7 @@ ...@@ -32,6 +32,7 @@
#include "RootInput.h" #include "RootInput.h"
#include "RootOutput.h" #include "RootOutput.h"
#include "TInteractionCoordinates.h" #include "TInteractionCoordinates.h"
#include "TInitialConditions.h"
// Use CLHEP System of unit and Physical Constant // Use CLHEP System of unit and Physical Constant
#include "CLHEP/Units/GlobalSystemOfUnits.h" #include "CLHEP/Units/GlobalSystemOfUnits.h"
...@@ -105,8 +106,15 @@ namespace ENERGYLOSS ...@@ -105,8 +106,15 @@ namespace ENERGYLOSS
{ {
// Declare your Energy loss here // Declare your Energy loss here
// EnergyLoss LightTargetCD2 = EnergyLoss("proton_cd2.txt", 100, 1, 1); // LISE++ // EnergyLoss LightTargetCD2 = EnergyLoss("proton_cd2.txt", 100, 1, 1); // LISE++
EnergyLoss LightTarget = EnergyLoss("proton_CD2.G4table", "G4Table", 1000); // G4 // For 132Sn(d,p)
EnergyLoss BeamTarget = EnergyLoss("Sn132[0.0]_CD2.G4table", "G4Table", 1000); // G4 // EnergyLoss LightTarget = EnergyLoss("proton_CD2.G4table", "G4Table", 1000); // G4
// EnergyLoss BeamTarget = EnergyLoss("Sn132[0.0]_CD2.G4table", "G4Table", 1000); // G4
// For 134Sn(p,t)
EnergyLoss LightTarget = EnergyLoss("triton_CH2.G4table", "G4Table", 1000); // G4
EnergyLoss BeamTarget = EnergyLoss("Sn134[0.0]_CH2.G4table", "G4Table", 1000); // G4
// For 132Sn(d,t)
// EnergyLoss LightTarget = EnergyLoss("triton_CD2.G4table", "G4Table", 1000); // G4
// EnergyLoss BeamTarget = EnergyLoss("Sn132[0.0]_CD2.G4table", "G4Table", 1000); // G4
} }
using namespace ENERGYLOSS ; using namespace ENERGYLOSS ;
......
NPAnalysis/Gaspard/macros/DisplayInputCrossSection.C 0 → 100644
+ 148
− 0
View file @ d6e510ed
void DisplayInputCrossSection()
{
// Path to cross-section files
TString path = gSystem->Getenv("NPTOOL");
path += "/Inputs/CrossSection/";
// Read cross-sections 132Sn(d,p)
/* TGraph *gr1 = new TGraph(path + "132Sndp_5A_MeV_3p3_ZR_FRC.lis");
TGraph *gr2 = new TGraph(path + "132Sndp_5A_MeV_2f7_ZR_FRC.lis");
TGraph *gr3 = new TGraph(path + "132Sndp_10A_MeV_3p3_ZR_FRC.lis");
TGraph *gr4 = new TGraph(path + "132Sndp_10A_MeV_2f7_ZR_FRC.lis");*/
// Jacques
TGraph *gr1 = ReadCrossSection("132Sndp_5A_MeV_3p3_ZR_FRC.lis");
TGraph *gr2 = ReadCrossSection("132Sndp_5A_MeV_2f7_ZR_FRC.lis");
TGraph *gr3 = ReadCrossSection("132Sndp_10A_MeV_3p3_ZR_FRC.lis");
TGraph *gr4 = ReadCrossSection("132Sndp_10A_MeV_2f7_ZR_FRC.lis");
// Didier
TGraph *gr5 = new TGraph(path + "sn132dp_gs_10AMeV.txt");
// Read cross-section 134Sn(p,t)
// Didier
TGraph *grpt1 = ReadCrossSectionPT("CS_Ep10MeV_sn134pt_gs_1h9demi.dat");
TGraph *grpt2 = ReadCrossSectionPT("CS_Ep15MeV_sn134pt_gs_1h9demi.dat");
TGraph *grpt3 = ReadCrossSectionPT("CS_Ep20MeV_sn134pt_gs_1h9demi.dat");
// Read cross-section 132Sn(d,d)
// Angel
TGraph *grdd = new TGraph(path + "132Sndd_10A_MeV_ruth.dat");
TGraph *grpp = new TGraph(path + "132Snpp_10A_MeV_ruth.dat");
// Draw cross-sections
TCanvas *can = new TCanvas("can");
can->SetLogy();
can->Draw();
// TH2F *hframe = new TH2F("hframe", "^{2}H(^{132}Sn,p)^{133}Sn", 180, 0, 180, 100, 1e-2, 100);
// TH2F *hframe = new TH2F("hframe", "", 180, 0, 180, 100, 1e-2, 100);
// TH2F *hframe = new TH2F("hframe", "^{1}H(^{134}Sn,t)^{132}Sn_{g.s.}", 180, 0, 180, 100, 1e-8, 1e-5);
// TH2F *hframe = new TH2F("hframe", "", 180, 0, 180, 100, 1e-8, 1e-5);
TH2F *hframe = new TH2F("hframe", "", 180, 0, 180, 100, 1e-3, 10);
hframe->Draw();
hframe->SetXTitle("#Theta_{c.m.} [deg]");
// hframe->SetYTitle("d#sigma/d#Omega [mb/sr]");
hframe->SetYTitle("d#sigma/d#Omega / (d#sigma/d#Omega)_{R}");
// hframe->SetYTitle("#propto d#sigma/d#Omega");
/* grpt1->SetLineColor(kRed); grpt1->Draw("l");
grpt2->SetLineColor(kMagenta); grpt2->Draw("l");
grpt3->SetLineColor(kBlue); grpt3->Draw("l");*/
/* gr1->SetLineColor(kRed); gr1->SetLineStyle(2); gr1->Draw("l");
gr2->SetLineColor(kRed); gr2->Draw("l");
gr3->SetLineColor(kBlue); gr3->SetLineStyle(2); gr3->Draw("l");
gr4->SetLineColor(kBlue); gr4->Draw("l");*/
// gr5->Draw("l");
// grdd->Draw("l");
grpp->Draw("l");
// TLegend
TLegend *leg = new TLegend(0.50, 0.64, 0.82, 0.84);
// leg->AddEntry(grdd, "10 MeV/u", "l");
leg->AddEntry(grpp, "10 MeV/u", "l");
/* TLegend *leg = new TLegend(0.50, 0.64, 0.82, 0.84);
leg->AddEntry(grpt1, "1h9/2 10 MeV/u", "l");
leg->AddEntry(grpt2, "1h9/2 15 MeV/u", "l");
leg->AddEntry(grpt3, "1h9/2 20 MeV/u", "l");*/
/* TLegend *leg = new TLegend(0.16, 0.17, 0.48, 0.37);
leg->AddEntry(gr1, "3p3/2 5 MeV/u", "l");
leg->AddEntry(gr2, "2f7/2 5 MeV/u", "l");
leg->AddEntry(gr3, "3p3/2 10 MeV/u", "l");
leg->AddEntry(gr4, "2f7/2 10 MeV/u", "l");*/
leg->SetBorderSize(1);
leg->Draw();
/* TMultiGraph *mgr = new TMultiGraph();
mgr->Add(gr1, "lp");
mgr->Add(gr2, "lp");
mgr->Add(gr3, "lp");
mgr->Add(gr4, "lp");
mgr->Draw("a*");*/
// gr1->Draw("alp");
}
TGraph* ReadCrossSection(const char* fname)
{
// Path to cross-section files
TString path = gSystem->Getenv("NPTOOL");
path += "/Inputs/CrossSection/";
// Open file
ifstream fich;
fich.open(path + fname);
if (!fich) cout << "Probleme d'ouverture dans le fichier " << fname << endl;
// Read file
Double_t angle, sigma;
Int_t nlines = 0;
TGraph *gr = new TGraph();
while (fich >> angle >> sigma) {
gr->SetPoint(nlines++, angle, sigma * 15); // 15: fm^2 -> mb + D0^2
}
// Close file
fich.close();
// TGraph name
gr->SetTitle(fname);
// test pour savoir si on a bien rempli le TGraph
if (gr->GetN() == 0)
cout << "Mauvaise lecture du fichier --> probablement mauvais format" << endl;
return gr;
}
TGraph* ReadCrossSectionPT(const char* fname)
{
// Path to cross-section files
TString path = gSystem->Getenv("NPTOOL");
path += "/Inputs/CrossSection/";
// Open file
ifstream fich;
fich.open(path + fname);
if (!fich) cout << "Probleme d'ouverture dans le fichier " << fname << endl;
// Read file
Double_t angle, sigma, dum;
Int_t nlines = 0;
TGraph *gr = new TGraph();
// while (fich >> angle >> sigma >> dum >> dum >> dum >> dum >> dum >> dum >> dum >> dum >> dum) {
while (fich >> angle >> sigma) {
gr->SetPoint(nlines++, angle, sigma);
}
// Close file
fich.close();
// TGraph name
gr->SetTitle(fname);
// test pour savoir si on a bien rempli le TGraph
if (gr->GetN() == 0)
cout << "Mauvaise lecture du fichier --> probablement mauvais format" << endl;
return gr;
}
NPAnalysis/Gaspard/src/Analysis.cc
+ 74
− 8
View file @ d6e510ed
...@@ -30,6 +30,17 @@ int main(int argc,char** argv) ...@@ -30,6 +30,17 @@ int main(int argc,char** argv)
NPA::DetectorManager* myDetector = new DetectorManager; NPA::DetectorManager* myDetector = new DetectorManager;
myDetector->ReadConfigurationFile(detectorfileName); myDetector->ReadConfigurationFile(detectorfileName);
// Calculate beam energy at target middle
// Get nominal beam energy
Double_t BeamEnergyNominal = myReaction->GetBeamEnergy() * MeV;
cout << BeamEnergyNominal << endl;
// Slow beam at target middle
Double_t BeamEnergy = BeamEnergyNominal - BeamTarget.Slow(BeamEnergyNominal, myDetector->GetTargetThickness()/2 * micrometer, 0);
// Double_t BeamEnergy = 1293.56 * MeV;
cout << BeamEnergy << endl;
// Set energy beam at target middle
myReaction->SetBeamEnergy(BeamEnergy);
// nominal beam energy // nominal beam energy
Double_t BeamEnergyNominal = myReaction->GetBeamEnergy() * MeV; Double_t BeamEnergyNominal = myReaction->GetBeamEnergy() * MeV;
cout << BeamEnergyNominal << endl; cout << BeamEnergyNominal << endl;
...@@ -51,9 +62,14 @@ int main(int argc,char** argv) ...@@ -51,9 +62,14 @@ int main(int argc,char** argv)
// Get GaspardTracker pointer // Get GaspardTracker pointer
GaspardTracker* GPDTrack = (GaspardTracker*) myDetector->m_Detector["GASPARD"]; GaspardTracker* GPDTrack = (GaspardTracker*) myDetector->m_Detector["GASPARD"];
// Get the TChain and treat it // Get the input TChain and treat it
TChain* chain = RootInput:: getInstance() -> GetChain(); TChain* chain = RootInput:: getInstance() -> GetChain();
// Connect TInitialConditions branch
TInitialConditions *initCond = 0;
chain->SetBranchAddress("InitialConditions", &initCond);
chain->SetBranchStatus("InitialConditions", 1);
// Connect TInteractionCoordinates branch // Connect TInteractionCoordinates branch
TInteractionCoordinates *interCoord = 0; TInteractionCoordinates *interCoord = 0;
chain->SetBranchAddress("InteractionCoordinates", &interCoord); chain->SetBranchAddress("InteractionCoordinates", &interCoord);
...@@ -63,6 +79,16 @@ int main(int argc,char** argv) ...@@ -63,6 +79,16 @@ int main(int argc,char** argv)
int nentries = chain->GetEntries(); int nentries = chain->GetEntries();
cout << "Number of entries to be analysed: " << nentries << endl; cout << "Number of entries to be analysed: " << nentries << endl;
// Default initialization
double XTarget = 0;
double YTarget = 0;
double BeamTheta = 0;
double BeamPhi = 0;
// random generator
TRandom3 *gene = new TRandom3();
// Loop on all events
for (int i = 0; i < nentries; i ++) { for (int i = 0; i < nentries; i ++) {
if (i%10000 == 0 && i!=0) cout << "\r" << i << " analyzed events" << flush; if (i%10000 == 0 && i!=0) cout << "\r" << i << " analyzed events" << flush;
chain -> GetEntry(i); chain -> GetEntry(i);
...@@ -75,22 +101,62 @@ int main(int argc,char** argv) ...@@ -75,22 +101,62 @@ int main(int argc,char** argv)
double E = GPDTrack->GetEnergyDeposit(); double E = GPDTrack->GetEnergyDeposit();
// if there is a hit in the detector array, treat it. // if there is a hit in the detector array, treat it.
double Theta, ThetaStrip; double Theta, ThetaStrip, angle, ThetaCM;
double DetecX, DetecY, DetecZ;
double r;
TVector3 A; TVector3 A;
if (E > -1000) { if (E > -1000) {
// Get exact scattering angle from TInteractionCoordinates object // Get c.m. angle
Theta = interCoord->GetDetectedAngleTheta(0) * deg; ThetaCM = initCond->GetICEmittedAngleThetaCM(0) * deg;
// Get interaction coordinates taking into account the strips // Get exact scattering angle from TInteractionCoordinates object
// Theta = interCoord->GetDetectedAngleTheta(0) * deg;
DetecX = interCoord->GetDetectedPositionX(0);
DetecY = interCoord->GetDetectedPositionY(0);
DetecZ = interCoord->GetDetectedPositionZ(0);
TVector3 Detec(DetecX, DetecY, DetecZ);
// Get interaction position in detector
// This takes into account the strips
A = GPDTrack->GetPositionOfInteraction(); A = GPDTrack->GetPositionOfInteraction();
// Calculate scattering angle
ThetaStrip = ThetaCalculation (A ,TVector3(0,0,1));
// Get beam interaction coordinates on target (from initial condition)
XTarget = initCond->GetICPositionX(0);
YTarget = initCond->GetICPositionY(0);
// cout << XTarget << " " << YTarget << endl;
BeamTheta = initCond->GetICIncidentAngleTheta(0)*deg;
BeamPhi = initCond->GetICIncidentAnglePhi(0)*deg;
TVector3 BeamDirection = TVector3(cos(BeamPhi)*sin(BeamTheta), sin(BeamPhi)*sin(BeamTheta), cos(BeamTheta));
// cout << BeamDirection.X() << " " << BeamDirection.Y() << " " << BeamDirection.Z() << endl;
// Hit direction taking into account beam position on target
TVector3 HitDirection = A - TVector3(XTarget, YTarget, 0);
// cout << "A: " << A.X() << " " << A.Y() << " " << A.Z() << endl;
// cout << "HitDirection: " << HitDirection.X() << " " << HitDirection.Y() << " " << HitDirection.Z() << endl;
// Calculate scattering angle w.r.t. optical beam axis (do not take into account beam position on target)
ThetaStrip = ThetaCalculation(A, TVector3(0,0,1));
Theta = ThetaCalculation(Detec, TVector3(0, 0, 1));
// Calculate scattering angle w.r.t. beam (ideal case)
// ThetaStrip = ThetaCalculation(HitDirection, BeamDirection);
// Theta = ThetaCalculation(Detec - TVector3(XTarget, YTarget, 0), BeamDirection);
// Calculate scattering angle w.r.t. beam (finite spatial resolution)
/* double resol = 800; // in micrometer
angle = gene->Rndm() * 2*3.14;
r = fabs(gene->Gaus(0, resol)) * micrometer;
ThetaStrip = ThetaCalculation(A - TVector3(XTarget + r*cos(angle), YTarget + r*sin(angle), 0), BeamDirection);
Theta = ThetaCalculation(Detec - TVector3(XTarget + r*cos(angle), YTarget + r*sin(angle), 0), BeamDirection);
*/
// Correct for energy loss in the target // Correct for energy loss in the target
E = LightTarget.EvaluateInitialEnergy(E, myDetector->GetTargetThickness()/2 * micrometer, ThetaStrip); E = LightTarget.EvaluateInitialEnergy(E, myDetector->GetTargetThickness()/2 * micrometer, ThetaStrip);
// Calculate excitation energy // Calculate excitation energy
if (Theta/deg > 90) { // if (Theta/deg > 150 && Theta/deg < 180) {
// if (Theta/deg < 60 && ThetaCM/deg < 90) {
// if (Theta/deg > 35 && Theta/deg < 45 && E/MeV < 17) {
// if (Theta/deg < 45) {
if (E/MeV < 38) {
// if (Theta/deg > 90) {
ExNoStrips = myReaction->ReconstructRelativistic(E, Theta / rad); ExNoStrips = myReaction->ReconstructRelativistic(E, Theta / rad);
Ex = myReaction->ReconstructRelativistic(E, ThetaStrip); Ex = myReaction->ReconstructRelativistic(E, ThetaStrip);
} }
......
NPAnalysis/macros/ControlSimu.C
+ 15
− 2
View file @ d6e510ed
...@@ -37,7 +37,6 @@ ...@@ -37,7 +37,6 @@
#include "TH1F.h" #include "TH1F.h"
#include "TH2F.h" #include "TH2F.h"
#include "TGaspardTrackerData.h"
#include "TInitialConditions.h" #include "TInitialConditions.h"
#include "TInteractionCoordinates.h" #include "TInteractionCoordinates.h"
...@@ -91,7 +90,7 @@ void ControlSimu(const char * fname = "mySimul") ...@@ -91,7 +90,7 @@ void ControlSimu(const char * fname = "mySimul")
hEmittanceXY -> Fill(initCond->GetICPositionX(0), initCond->GetICPositionY(0)); hEmittanceXY -> Fill(initCond->GetICPositionX(0), initCond->GetICPositionY(0));
hIncidentZ -> Fill(initCond->GetICPositionZ(0)); hIncidentZ -> Fill(initCond->GetICPositionZ(0));
hEmittanceXTheta -> Fill(initCond->GetICPositionX(0), initCond->GetICIncidentEmittanceTheta(0)); hEmittanceXTheta -> Fill(initCond->GetICPositionX(0), initCond->GetICIncidentEmittanceTheta(0));
hEmittanceYPhi -> Fill(initCond->GetICPositionX(0), initCond->GetICIncidentEmittancePhi(0)); hEmittanceYPhi -> Fill(initCond->GetICPositionY(0), initCond->GetICIncidentEmittancePhi(0));
hIncidentTheta -> Fill(initCond->GetICIncidentAngleTheta(0)); hIncidentTheta -> Fill(initCond->GetICIncidentAngleTheta(0));
hIncidentPhi -> Fill(initCond->GetICIncidentAnglePhi(0)); hIncidentPhi -> Fill(initCond->GetICIncidentAnglePhi(0));
// ejected particle // ejected particle
...@@ -141,4 +140,18 @@ void ControlSimu(const char * fname = "mySimul") ...@@ -141,4 +140,18 @@ void ControlSimu(const char * fname = "mySimul")
hEmittedPhiWF->Draw(); hEmittedPhiWF->Draw();
// hControlPhi->SetLineColor(kRed); // hControlPhi->SetLineColor(kRed);
// hControlPhi->Draw("same"); // hControlPhi->Draw("same");
// Display histograms
TCanvas *canvasxi32 = new TCanvas("canvas3", "Emitted particle properties", 300, 600);
canvas3->Divide(1,3);
canvas3->cd(1);
hEmittedThetaCM->SetXTitle("#Theta_{c.m.}");
hEmittedThetaCM->Draw();
canvas3->cd(2);
hEmittedThetaWF->Draw();
hEmittedThetaWF->SetXTitle("#Theta_{lab}");
canvas3->cd(3);
hEmittedETheta->SetXTitle("#Theta_{lab}");
hEmittedETheta->SetYTitle("E [MeV]");
hEmittedETheta->Draw();
} }
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment