forked from CAGRA-GrandRaiden/GRAnalyzer
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathhist.def
1921 lines (1920 loc) · 145 KB
/
hist.def
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
# Created by /home/g01/daq/bin/hm on aino-1, 21:13:25 17-Oct-2015
1 com Definition of Parameters
2. com Constants
3 def FALSE = 0
3 def TRUE = 1
3 def pi = 3.14159265358979323846
3 def d2r = pi/180
3 def r2d = 180/pi
3 def cos45 = 0.70710678118654752440
3 def sin45 = 0.70710678118654752440
3 com Velocity of light in vaccuum = 299792458 [m/s] (exact)
4 def c = 2.99792458e+8
3 com Fine structure constant = 1/137.0359895(61)
4 def alpha = 1/137.0359895
3 com Atomic mass unit = 931.49432(28) [MeV/c2]
4 def AMU = 931.49432
3 com Electron charge = 1.60217733(49)e-19 [C]
4 def e = 1.60217733e-19
3 com hbarc = 197.327053(59) [MeVfm]
4 def hbc = 197.327053
3 com Electron mass = 0.51099907(15) [MeV/c2]
4 def M_e = 0.51099907
3 com Proton mass = 938.27231(28) [MeV/c2]
4 def M_p = 938.27231
3 com Neutron mass = 939.56563(28) [MeV/c2]
4 def M_n = 939.56563
3 com Deuteron mass = 1875.61339(57) [MeV/c2]
4 def M_d = 1875.61339
3. com Nuclear masses [MeV/c2]
4 def M_3He = 2808.392
4 def M_4He = 3728.4
4 def M_10B = 10*AMU+12.0508-5*M_e
4 def M_10B = 10*AMU+12.0508-5*M_e
4 def M_11B = 11*AMU+8.6680-5*M_e
4 def M_12C = 12*AMU+0.000-6*M_e
4 def M_12N = 12*AMU+17.3381-7*M_e
4 def M_24Mg = 24*AMU-13.933-12*M_e
4 def M_26Mg = 26*AMU-16.214-14*M_e
4 def M_28Si = 28*AMU-21.49283-14*M_e
4 def M_40Ca = 40*AMU-34.846-20*M_e
4 def M_48Ca = 48*AMU-44.215-20*M_e
4 def M_58Ni = 58*AMU-60.2230-28*M_e
4 def M_90Zr = 90*AMU-86.4457-40*M_e
4 def M_90Y = 90*AMU-86.4807-39*M_e
4 def M_90Y = 197*AMU-31.157-79*M_e
4 def M_197Au = 197*AMU-31.157-79*M_e
4 def M_208Pb = 208*AMU-21.764-82*M_e
3. com Debug Variables
-4 def DEBUG_DUMP_VARIABLES = 1
-4 def DEBUG_DUMP_DEBUG = 1
4 def DEBUG_DUMP_SHOW_NONE = 1
-4 def DEBUG = GR_ADC
-4 def DEBUG_GR_N_OUTPUTRAY = 200
-4 def DEBUG_LAS_N_OUTPUTRAY = 1000
2 com Analysis Parameters
3. com Analyzing System
-4. com Blocks (ANA_BLK_START <= BLOCK_N <= ANA_BLK_END)
5 def ANA_BLK_START = 23000
5 def ANA_BLK_END = 24000
-4. com Blocks (ANA_BLK_START <= BLOCK_N <= ANA_BLK_END)
5 def ANA_BLK_START = 24000
5 def ANA_BLK_END = 25000
-4. com Blocks (ANA_BLK_START <= BLOCK_N <= ANA_BLK_END)
5 def ANA_BLK_START = 34000
5 def ANA_BLK_END = 35000
-4. com Blocks (ANA_BLK_START <= BLOCK_N <= ANA_BLK_END)
5 def ANA_BLK_START = 10000
5 def ANA_BLK_END = 20000
3. com GR
4 def ANALYZE_GR = !BLP
-4 def ANALYZE_GR = TRUE
4 def ANALYZE_GR_VDC = TRUE
-4 def ANALYZE_GR_VDC = GR_PI && GR_RF
4 def ANALYZE_GR_MATRIX = TRUE
4 com Cluster Analysis
5 def GR_VDC_ANA_MULTI_CLUST = TRUE # Analyze Multi-Cluster Event
-5 def GR_VDC_ANA_2HIT_CLUST = TRUE
5 def GR_VDC_MAX_CLUST_SIZE = 10 # Default Value
5 def GR_VDC_MIN_CLUST_SIZE = 2 # Default Value
5 def GR_VDC_MAX_CHI2 = 50 # Maximum allowed chi2 when counting number of combinations
3. com LAS
4 def ANALYZE_LAS = TRUE
-4 def ANALYZE_LAS_VDC = TRUE
4 def ANALYZE_LAS_VDC = LAS_RF_G
-4 def ANALYZE_LAS_VDC = LAS_PI
4 com Cluster Analysis
5 def LAS_VDC_ANA_MULTI_CLUST = TRUE # Analyze Multi-Cluster Event
-5 def LAS_VDC_ANA_2HIT_CLUST = TRUE
5 def LAS_VDC_MAX_CLUST_SIZE = 10 # Default Value
5 def LAS_VDC_MIN_CLUST_SIZE = 2 # Default Value
5 def LAS_VDC_MAX_CHI2 = 50. # Maximum allowed chi2 when counting number of combinations
-3. com GR-FPP
4 def ANALYZE_GR_FPP = TRUE
-4 def ANALYZE_GR_FPP = GF_ANA_G
3. com Multi-Cluster Analysis (Method#1)
-4. com Set#1 (see Analysis book pp.33)
5 def GF_ANA_MULTI_CLUST = TRUE # Analyze Multi-Cluster Event
5 def GF_V_POS = {1759.4, 0.0, 0.0} # Assumed vertex position plane in mm {z, theta, phi} (the center of the carbon scatterer)
5 def GF_MAX_CHI2 = 9 # Maximum allowed chi2 when counting number of combinations
5 def GF_MAX_NCOMB = 4 # Maximum allowed number of combinations in wire positions
5 def GF_MAX_NCLUST = 2 # Maximum allowed number of clusters in a plane
-5 def GF_VD_RES = 10 # Vertex distance resolution (in mm)
5 def GF_VX_RES = 7 # Vertex X position resolution (in mm)
5 def GF_VY_RES = 7 # Vertex Y position resolution (in mm)
4 com Set#2 (see Analysis book pp.33)
5 def GF_ANA_MULTI_CLUST = TRUE # Analyze Multi-Cluster Event
5 def GF_V_POS = {1759.4, 0.0, 0.0} # Assumed vertex position plane in mm {z, theta, phi} (the center of the carbon scatterer)
5 def GF_MAX_CHI2 = 9 # Maximum allowed chi2 when counting number of combinations
5 def GF_MAX_NCOMB = 9 # Maximum allowed number of combinations in wire positions
5 def GF_MAX_NCLUST = 4 # Maximum allowed number of clusters in a plane
-5 def GF_VD_RES = 10 # Vertex distance resolution (in mm)
5 def GF_VX_RES = 7 # Vertex X position resolution (in mm)
5 def GF_VY_RES = 7 # Vertex Y position resolution (in mm)
-3. com Multi-Cluster Analysis (Method#2)
4 def GF_ANA_MULTI_CLUST = TRUE # Analyze Multi-Cluster Event
4 def GF_V_POS = {1759.4, 0.0, 0.0} # Assumed vertex position plane in mm {z, theta, phi} (the center of the carbon scatterer)
-3. com Scaler Output
-4 def PRINT_SCALER = TRUE
-4 def PRINT_SCALER = FALSE
4 def PRINT_SCALER = EVENT_BEND && BLOCK_N[0]%100==0
3. com DST Output
4 def DST_CONDITION = 1 # No Condition
4. com DST Step #2 (VDC-Ray-Analyzed, GR data)
5 def DST_VAR = "QTC_TRAILING_CH:QTC_TRAILING_TDC:QTC_LEADING_CH:QTC_LEADING_TDC:GR_RAYID:GR_RF:GR_X:GR_Y:GR_TH:GR_PH:GR_ADC1:GR_ADC2:GR_ADC3:GR_ADC4:GR_ADC5:GR_ADC6:GR_TDC1:GR_TDC2:GR_TDC3:GR_TDC4:GR_TDC5:GR_TDC6:GR_TPOS1:GR_TPOS2:GR_MADC1:GR_MADC2:GR_MADC3:GR_WIRE_X1:GR_WIRE_U1:GR_WIRE_X2:GR_WIRE_U2:GR_WTDC_X1:GR_WTDC_U1:GR_WTDC_X2:GR_WTDC_U2:GR_TDCR_X1:GR_TDCR_U1:GR_TDCR_X2:GR_TDCR_U2:LAS_RAYID:LAS_RF1:LAS_RF2:LAS_RF3:LAS_X:LAS_Y:LAS_TH:LAS_PH:LAS_ADC1:LAS_ADC2:LAS_ADC3:LAS_ADC4:LAS_ADC5:LAS_ADC6:LAS_ADC7:LAS_ADC8:LAS_ADC9:LAS_ADC10:LAS_ADC11:LAS_ADC12:LAS_TDC1:LAS_TDC2:LAS_TDC3:LAS_TDC4:LAS_TDC5:LAS_TDC6:LAS_TDC7:LAS_TDC8:LAS_TDC9:LAS_TDC10:LAS_TDC11:LAS_TDC12:LAS_TPOS1:LAS_TPOS2:LAS_TPOS3:LAS_TPOS4:LAS_TPOS5:LAS_TPOS6:LAS_MADC1:LAS_MADC2:LAS_MADC3:LAS_MADC4:LAS_MADC5:LAS_MADC6:LAS_WIRE_X1:LAS_WIRE_U1:LAS_WIRE_V1:LAS_WIRE_X2:LAS_WIRE_U2:LAS_WIRE_V2:LAS_TDCR_X1:LAS_TDCR_U1:LAS_TDCR_V1:LAS_TDCR_X2:LAS_TDCR_U2:LAS_TDCR_V2"
#"QTC_TRAILING_CH:QTC_TRAILING_TDC:QTC_LEADING_CH:QTC_LEADING_TDC:GR_RAYID:GR_RF:GR_X:GR_Y:GR_TH:GR_PH:GR_ADC1:GR_ADC2:GR_ADC3:GR_ADC4:GR_ADC5:GR_ADC6:GR_TDC1:GR_TDC2:GR_TDC3:GR_TDC4:GR_TDC5:GR_TDC6:GR_TPOS1:GR_TPOS2:GR_MADC1:GR_MADC2:GR_MADC3:GR_WIRE_X1:GR_WIRE_U1:GR_WIRE_X2:GR_WIRE_U2:GR_WTDC_X1:GR_WTDC_U1:GR_WTDC_X2:GR_WTDC_U2:GR_TDCR_X1:GR_TDCR_U1:GR_TDCR_X2:GR_TDCR_U2:LAS_RAYID:LAS_RF1:LAS_RF2:LAS_RF3:LAS_X:LAS_Y:LAS_TH:LAS_PH:LAS_ADC1:LAS_ADC2:LAS_ADC3:LAS_ADC4:LAS_ADC5:LAS_ADC6:LAS_ADC7:LAS_ADC8:LAS_ADC9:LAS_ADC10:LAS_ADC11:LAS_ADC12:LAS_TDC1:LAS_TDC2:LAS_TDC3:LAS_TDC4:LAS_TDC5:LAS_TDC6:LAS_TDC7:LAS_TDC8:LAS_TDC9:LAS_TDC10:LAS_TDC11:LAS_TDC12:LAS_TPOS1:LAS_TPOS2:LAS_TPOS3:LAS_TPOS4:LAS_TPOS5:LAS_TPOS6:LAS_MADC1:LAS_MADC2:LAS_MADC3:LAS_MADC4:LAS_MADC5:LAS_MADC6:LAS_WIRE_X1:LAS_WIRE_U1:LAS_WIRE_V1:LAS_WIRE_X2:LAS_WIRE_U2:LAS_WIRE_V2:LAS_TDCR_X1:LAS_TDCR_U1:LAS_TDCR_V1:LAS_TDCR_X2:LAS_TDCR_U2:LAS_TDCR_V2:BLP_ADC2:BLP_ADC3:BLP_ADC4:BLP_ADC5:BLP_ADC6:BLP_ADC7:BLP_ADC8:BLP_TDC1:BLP_TDC2:BLP_TDC3:BLP_TDC4:BLP_TDC5:BLP_TDC6:BLP_TDC7:BLP_TDC8:MADC32_0_CH:MADC32_0_ADC:MADC32_1_CH:MADC32_1_ADC:MADC32_2_CH:MADC32_2_ADC:SSD_TRAI_CH:SSD_TRAI_TDC:SSD_LEAD_CH:SSD_LEAD_TDC:GR_TIMESTAMP"
5. com Variables for ROOT conversion (arrays reduced to single variables)
5 def GR_ADC1 = GR_FERA_ADC[0x10]
5 def GR_ADC2 = GR_FERA_ADC[0x11]
5 def GR_ADC3 = GR_FERA_ADC[0x12]
5 def GR_ADC4 = GR_FERA_ADC[0x13]
5 def GR_ADC5 = GR_FERA_ADC[0x20]
5 def GR_ADC6 = GR_FERA_ADC[0x28]
5 def GR_TDC1 = GR_FERA_TDC[0x16]
5 def GR_TDC2 = GR_FERA_TDC[0x11]
5 def GR_TDC3 = GR_FERA_TDC[0x12]
5 def GR_TDC4 = GR_FERA_TDC[0x13]
5 def GR_TDC5 = GR_FERA_TDC[0x20]
5 def GR_TDC6 = GR_FERA_TDC[0x28]
5 def LAS_ADC1 = LAS_FERA_ADC[0x10]
5 def LAS_ADC2 = LAS_FERA_ADC[0x11]
5 def LAS_ADC3 = LAS_FERA_ADC[0x12]
5 def LAS_ADC4 = LAS_FERA_ADC[0x13]
5 def LAS_ADC5 = LAS_FERA_ADC[0x14]
5 def LAS_ADC6 = LAS_FERA_ADC[0x15]
5 def LAS_ADC7 = LAS_FERA_ADC[0x16]
5 def LAS_ADC8 = LAS_FERA_ADC[0x17]
5 def LAS_ADC9 = LAS_FERA_ADC[0x18]
5 def LAS_ADC10= LAS_FERA_ADC[0x19]
5 def LAS_ADC11= LAS_FERA_ADC[0x1a]
5 def LAS_ADC12= LAS_FERA_ADC[0x1b]
5 def LAS_TDC1 = LAS_FERA_TDC[0x10]
5 def LAS_TDC2 = LAS_FERA_TDC[0x12]
5 def LAS_TDC3 = LAS_FERA_TDC[0x14]
5 def LAS_TDC4 = LAS_FERA_TDC[0x16]
5 def LAS_TDC5 = LAS_FERA_TDC[0x18]
5 def LAS_TDC6 = LAS_FERA_TDC[0x1a]
5 def LAS_TDC7 = LAS_FERA_TDC[0x1c]
5 def LAS_TDC8 = LAS_FERA_TDC[0x1e]
5 def LAS_TDC9 = LAS_FERA_TDC[0x20]
5 def LAS_TDC10= LAS_FERA_TDC[0x22]
5 def LAS_TDC11= LAS_FERA_TDC[0x24]
5 def LAS_TDC12= LAS_FERA_TDC[0x26]
5 def LAS_RF1 = LAS_FERA_TDC[0x2a]
5 def LAS_RF2 = LAS_FERA_TDC[0x2c]
5 def LAS_RF3 = LAS_FERA_TDC[0x2e]
5 def GR_TPOS1 = GR_TDC[2]-GR_TDC[1]
5 def GR_TPOS2 = GR_TDC[4]-GR_TDC[3]
5 def GR_MADC1 = sqrt(GR_ADC[1]*GR_ADC[2])
5 def GR_MADC2 = sqrt(GR_ADC[3]*GR_ADC[4])
5 def GR_MADC3 = sqrt(GR_ADC[5]*GR_ADC[6])
5 def BLP_ADC1 = GR_ADC_OLD[0]
5 def BLP_ADC2 = GR_ADC_OLD[1]
5 def BLP_ADC3 = GR_ADC_OLD[2]
5 def BLP_ADC4 = GR_ADC_OLD[3]
5 def BLP_ADC5 = GR_ADC_OLD[4]
5 def BLP_ADC6 = GR_ADC_OLD[5]
5 def BLP_ADC7 = GR_ADC_OLD[6]
5 def BLP_ADC8 = GR_ADC_OLD[7]
5 def BLP_TDC1 = GR_TDC_OLD[0]
5 def BLP_TDC2 = GR_TDC_OLD[1]
5 def BLP_TDC3 = GR_TDC_OLD[2]
5 def BLP_TDC4 = GR_TDC_OLD[3]
5 def BLP_TDC5 = GR_TDC_OLD[4]
5 def BLP_TDC6 = GR_TDC_OLD[5]
5 def BLP_TDC7 = GR_TDC_OLD[6]
5 def BLP_TDC8 = GR_TDC_OLD[7]
5 def LAS_TPOS1 = LAS_TDC[2]-LAS_TDC[1]
5 def LAS_TPOS2 = LAS_TDC[4]-LAS_TDC[3]
5 def LAS_TPOS3 = LAS_TDC[6]-LAS_TDC[5]
5 def LAS_TPOS4 = LAS_TDC[8]-LAS_TDC[7]
5 def LAS_TPOS5 = LAS_TDC[10]-LAS_TDC[9]
5 def LAS_TPOS6 = LAS_TDC[12]-LAS_TDC[11]
5 def LAS_MADC1 = sqrt(LAS_ADC[1]*LAS_ADC[2])
5 def LAS_MADC2 = sqrt(LAS_ADC[3]*LAS_ADC[4])
5 def LAS_MADC3 = sqrt(LAS_ADC[5]*LAS_ADC[6])
5 def LAS_MADC4 = sqrt(LAS_ADC[7]*LAS_ADC[8])
5 def LAS_MADC5 = sqrt(LAS_ADC[9]*LAS_ADC[10])
5 def LAS_MADC6 = sqrt(LAS_ADC[11]*LAS_ADC[12])
3. com Detector Configurations
4 com GR Chamber Configuration
5 def VDC_U_DX = 0
5 def GR_PL_X1 = {573.0+VDC_DX,0,0, -6, 0.00, 10.}
5 def GR_PL_U1 = {602.9+VDC_DX+VDC_U_DX,0,20, -4, -48.19, 10.}
5 def GR_PL_X2 = {323.0+VDC_DX,0,250, -6, 0.00, 10}
5 def GR_PL_U2 = {352.7+VDC_DX+VDC_U_DX,0,270, -4, -48.19, 10.}
-5 def GR_VDC_TILT = {0,45.25*d2r,0} # E316 09nov
5 com Angular Resolution of each Plane [deg]
6 def GR_ARES = {1.05, 1.97, 1.97, 1.05}
4 com LAS Chamber Configuration (UV)
5 def LAS_VDC_DY = +38.671 # The LAS detector is lifted by ~40mm (by Noro group). E282 08Nov
-5 def LAS_PL_X1 = {-757.76,LAS_VDC_DY,0, 6., 0.0, 10.}
5 def LAS_PL_U1 = {-825.76,LAS_VDC_DY,20, 6.,-31.0, 10.}
5 def LAS_PL_V1 = {-821.26,LAS_VDC_DY,40, 6.,+31.0, 10.} # E249 05Jul
-5 def LAS_PL_X2 = {-508.36,LAS_VDC_DY,164, 6., 0.0, 10.}
5 def LAS_PL_U2 = {-576.36,LAS_VDC_DY,184, 6.,-31.0, 10.}
5 def LAS_PL_V2 = {-571.86,LAS_VDC_DY,204, 6.,+31.0, 10.}
5 def LAS_VDC_TILT = {0,-54.04*d2r,0}
5 com Angular Resolution of each Plane [deg]
6 def LAS_ARES = {1.05, 1.97, 1.97, 1.05}
3. com GR-FPP Chamber Configuration
-4 com E316 09nov (by offline calibration of Run#9049)
5 def GF3_DX = 357.69 + 0.08
5 def GF3_DY = 4.05 - 0.73
5 def GF4_DX = GF3_DX + 0.25 - 2.97
5 def GF4_DY = GF3_DY - 0.52 - 0.58
-4 com E282&E316 Common
5 def GF_PL_X1 = { 470.00,0,1210.5, -2, 0, 1.}
5 def GF_PL_X2 = { 469.97,0,1466.1, -2, 0, 1.}
5 def GF_PL_U3 = {-903.68+GF3_DX,GF3_DY,1873.0, 2, -45, 1.}
5 def GF_PL_V3 = { 903.68+GF3_DX,GF3_DY,1885.0, -2, 45, 1.}
5 def GF_PL_U4 = {-994.19+GF4_DX,GF4_DY,2233.0, 2, -45, 1.}
5 def GF_PL_V4 = { 994.19+GF4_DX,GF4_DY,2245.0, -2, 45, 1.}
3. com Projection Plane {z,theta,phi}
4 def GR_N_VPLANE = 8
5 com at Z=0
6 def GR_PLANE_0 = {0,GR_VDC_TILT[1],0}
5 com at X1
6 def GR_PLANE_1 = {GR_PL_X1[2],GR_VDC_TILT[1],0}
5 com at X2
6 def GR_PLANE_2 = {GR_PL_X2[2],GR_VDC_TILT[1],0}
5 com at PS1
6 def GR_PLANE_3 = {452,GR_VDC_TILT[1],0}
5 com at PS3
6 def GR_PLANE_4 = {1730,0,0}
5 com at Carbon Block
6 def GR_PLANE_5 = {1790,0,0}
5 com at FPP-U3
6 def GR_PLANE_6 = {1873,0,0}
5 com at FPP-U4
6 def GR_PLANE_7 = {2233,0,0}
4 def LAS_N_VPLANE = 2
5 def LAS_PLANE_0 = {0,LAS_VDC_TILT[1],0}
5 def LAS_PLANE_1 = {LAS_PL_X1[2],LAS_VDC_TILT[1],0}
4 def GF_OUT_N_VPLANE = 2
5 com at FPP-U3
6 def GF_OUT_PLANE_0 = {1873,0,0}
5 com at FPP-U4
6 def GF_OUT_PLANE_1 = {2233,0,0}
2. com GR Trace-Back Matrix
-3 def GR_EN = GR_XC*0.0265+9.391 # 8deg
3 def GR_EN = GR_XC*0.0247+18.147-0.4
-3 def GR_YC = GR_Y+700.*GR_PH-6.4 # for natC E421 E422
3 def GR_YC = GR_Y+700.*GR_PH-(6.4-4.7) # for main target E421 E422
3 def GR_PHC = GR_PH-0.000000001*GR_Y-0.002 # for main target E421 E422
-3 def GR_PHC = GR_PH-0.0001*GR_Y-0.002 # for main target E421 E422, 94Zr
2. com XC for e421online
3 def GR_XC = GR_X - polynom(GR_TH*r2d, GR_XC_C)
-3 def GR_XC = GR_X - polynom(GR_THC*r2d, GR_XC_C)
-3 def GR_XC_C = { 0,14.488,-1.2052} # 0deg
-3 def GR_XC_C = { 0,14.373,-0.70152} # Au elstic 8deg
3 def GR_XC_C = { 0,14.373-0.886295,-0.70152-0.00248843} # Au elstic 8deg on 7-June-2015
-3 def GR_XC_C = { 0,11.45,-0.5225} # C elstic 8deg
-3 def GR_XC_C = { -0.375, -7.54} # Au(d,d0) at 10 deg
-3 def GR_XC_C = { 0.00, +6} # Au(d,d0) at 10 deg
-2. com Matrix for E282-mag2 070507 for NIM
-4. com Matrix for A (for E282-mag2 070507 for NIM)
5 def GR_AMAT_01 = { 0000, 1.61034e-02}
5 def GR_AMAT_02 = { 0100, -4.17069e-01}
5 def GR_AMAT_03 = { 1000, 2.36614e-05}
-4. com Matrix for B (for E282-mag2 070507 for NIM)
5 def GR_BMAT_01 = { 0000, -6.80045e-03}
5 def GR_BMAT_02 = { 0000, -1.75957e-03*RUN_LAS_YG}
5 def GR_BMAT_03 = { 0001, 2.46027e+00}
5 def GR_BMAT_04 = { 0010, -3.26673e-03}
5 def GR_BMAT_05 = { 0011, -4.29038e-03}
5 def GR_BMAT_06 = { 0100, -2.76651e-02}
5 def GR_BMAT_07 = { 0101, 2.37656e+01}
5 def GR_BMAT_08 = { 0110, 1.20404e-02}
5 def GR_BMAT_09 = { 0111, 8.86801e-02}
5 def GR_BMAT_10 = { 1000, 7.38007e-06}
5 def GR_BMAT_11 = { 1000, -1.16108e-07*RUN_LAS_YG}
5 def GR_BMAT_12 = { 1001, -1.75519e-03}
5 def GR_BMAT_13 = { 1010, -1.86523e-06}
5 def GR_BMAT_14 = { 1011, 2.90648e-06}
5 def GR_BMAT_15 = { 1100, 4.98760e-05}
5 def GR_BMAT_16 = { 1101, 1.02579e-02}
5 def GR_BMAT_17 = { 1110, 6.02465e-06}
5 def GR_BMAT_18 = { 1111, 4.02911e-04}
-3 com Matrix for E316 09Nov
4. com Matrix for A (not calibratied yet)
5 def GR_AMAT_01 = { 0000, 1.61034e-02}
5 def GR_AMAT_02 = { 0100, -4.17069e-01}
5 def GR_AMAT_03 = { 1000, 2.36614e-05}
5 def GR_AMAT_04 = { 0000, 1.09e-03} # Offset due to mis-alignment of the sieve-slit, see log note vol1-p.36, analysis note p.58
4. com Matrix for B (not calibrated yet)
-5 def GR_BMAT_01 = { 0000, -9.65234e-03+0.008}
5 def GR_BMAT_01 = { 0000, -9.65234e-03}
5 def GR_BMAT_02 = { 0000, -1.07086e-03*RUN_LAS_YG}
5 def GR_BMAT_03 = { 0001, 4.17063e+00}
5 def GR_BMAT_04 = { 0010, -1.99037e-03}
5 def GR_BMAT_05 = { 1000, 6.16186e-06}
5 def GR_BMAT_06 = { 1000, 6.17463e-07*RUN_LAS_YG}
5 def GR_BMAT_07 = { 0000, +(2.15e-04*RUN_LAS_YG)} # calculated ideal slope, see analysis note p.43
4. com Matrix for YC (not calibrated yet)
5 def GR_YMAT_01 = { 0010, 1.00000e+00}
5 def GR_YMAT_02 = { 0000, -3.74081e-01+2.720}
5 def GR_YMAT_03 = { 0000, 1.05732e+00*RUN_LAS_YG}
5 def GR_YMAT_04 = { 0001, 7.17335e+02+700}
5 def GR_YMAT_05 = { 0100, -5.20327e+01}
5 def GR_YMAT_06 = { 0101, 9.00266e+03}
5 def GR_YMAT_07 = { 1000, 1.24845e-02}
5 def GR_YMAT_08 = { 1001, -3.22387e+00}
5 def GR_YMAT_09 = { 1101, -1.37672e+01}
4 com Matrix for Xc
-5. com no correction
6 def GR_XMAT_01 = { 1000, 1.000e-00 }
5. com xc-th (not calibrated yet)
6 def GR_XMAT_01 = { 0000, 0.00000e+00}
6 def GR_XMAT_02 = { 0100, -4.95293e+02}
6 def GR_XMAT_03 = { 0200, 5.35956e+03}
6 def GR_XMAT_04 = { 0300, -3.27460e+04}
6 def GR_XMAT_05 = { 0400, 2.02434e+05}
6 def GR_XMAT_06 = { 1000, 1.00000e+00}
6 def GR_XMAT_07 = { 1100, -7.37824e-01}
6 def GR_XMAT_08 = { 1200, 2.79054e+00}
6 def GR_XMAT_09 = { 1300, -2.76952e+01}
6 def GR_XMAT_10 = { 1400, 4.44175e+02}
4. com Matrix for PHC (not calibrated yet)
5 def GR_PHMAT_01 = { 0000, 0.000 }
5 def GR_PHMAT_02 = { 0001, 1.000 }
5 def GR_PHMAT_03 = { 0010, -2.27E-4}
2. com LAS Trace-Back Matrix
3 com converted on 27-JUL-2000 from jul93las_xmat.par by Noro-group
4 com akv1=0.2e-4 k-value of this matrix
2. com Matrix for X
3 def LAS_XMAT_01 = { 0000, 0.22806E-01 }
3 def LAS_XMAT_02 = { 1000, 0.23392E-03 }
3 def LAS_XMAT_03 = { 0100, -0.29214E-01 }
3 def LAS_XMAT_04 = { 0010, 0.90449E-05 }
3 def LAS_XMAT_05 = { 0001, 0.61528E-02 }
3 def LAS_XMAT_06 = { 2000, 0.65014E-07 }
3 def LAS_XMAT_07 = { 1100, -0.33089E-05 }
3 def LAS_XMAT_08 = { 1010, -0.10777E-07 }
3 def LAS_XMAT_09 = { 1001, 0.11167E-04 }
3 def LAS_XMAT_10 = { 0200, 0.16362E-01 }
3 def LAS_XMAT_11 = { 0110, -0.76469E-05 }
3 def LAS_XMAT_12 = { 0101, -0.14984E-01 }
3 def LAS_XMAT_13 = { 0020, 0.89652E-07 }
3 def LAS_XMAT_14 = { 0011, -0.22450E-04 }
3 def LAS_XMAT_15 = { 0002, -0.10874E-02 }
3 def LAS_XMAT_16 = { 3000, 0.49676E-11 }
3 def LAS_XMAT_17 = { 2100, -0.48154E-07 }
3 def LAS_XMAT_18 = { 2010, -0.44922E-11 }
3 def LAS_XMAT_19 = { 2001, -0.21891E-08 }
3 def LAS_XMAT_20 = { 1200, -0.77687E-05 }
3 def LAS_XMAT_21 = { 1110, 0.61315E-08 }
3 def LAS_XMAT_22 = { 1101, -0.83385E-05 }
3 def LAS_XMAT_23 = { 1020, 0.40743E-10 }
3 def LAS_XMAT_24 = { 1011, 0.20944E-07 }
3 def LAS_XMAT_25 = { 1002, 0.17309E-05 }
3 def LAS_XMAT_26 = { 0300, -0.48981E-02 }
3 def LAS_XMAT_27 = { 0210, 0.20841E-05 }
3 def LAS_XMAT_28 = { 0201, 0.53219E-02 }
3 def LAS_XMAT_29 = { 0120, 0.32648E-08 }
3 def LAS_XMAT_30 = { 0111, 0.11717E-04 }
3 def LAS_XMAT_31 = { 0102, -0.73120E-03 }
3 def LAS_XMAT_32 = { 0021, 0.12566E-06 }
3 def LAS_XMAT_33 = { 0012, -0.52859E-04 }
3 def LAS_XMAT_34 = { 0003, 0.35165E-01 }
3 def LAS_XMAT_35 = { 4000, -0.35735E-13 }
3 def LAS_XMAT_36 = { 0400, -0.13701E-03 }
3 def LAS_XMAT_37 = { 0040, -0.66081E-11 }
3 def LAS_XMAT_38 = { 0004, 0.10286E-01 }
3 def LAS_XMAT_39 = { 3100, 0.30037E-10 }
3 def LAS_XMAT_40 = { 2200, -0.32596E-08 }
3 def LAS_XMAT_41 = { 1300, 0.75084E-05 }
2. com Matrix for A
3 def LAS_AMAT_01 = { 0000, -0.90882E+03}
3 def LAS_AMAT_02 = { 1000, -0.98402E+00}
3 def LAS_AMAT_03 = { 0100, 0.21078E+04}
3 def LAS_AMAT_04 = { 0010, -0.81389E-02}
3 def LAS_AMAT_05 = { 0001, -0.15176E+03}
3 def LAS_AMAT_06 = { 2000, -0.47535E-03}
3 def LAS_AMAT_07 = { 1100, 0.20539E+01}
3 def LAS_AMAT_08 = { 1010, 0.82019E-05}
3 def LAS_AMAT_09 = { 1001, -0.10719E+00}
3 def LAS_AMAT_10 = { 0200, -0.20867E+04}
3 def LAS_AMAT_11 = { 0110, -0.39076E-02}
3 def LAS_AMAT_12 = { 0101, 0.23063E+03}
3 def LAS_AMAT_13 = { 0020, 0.17441E-02}
3 def LAS_AMAT_14 = { 0011, 0.17144E+00}
3 def LAS_AMAT_15 = { 0002, 0.46881E+02}
3 def LAS_AMAT_16 = { 3000, -0.12898E-06}
3 def LAS_AMAT_17 = { 2100, 0.72953E-03}
3 def LAS_AMAT_18 = { 2010, -0.50616E-09}
3 def LAS_AMAT_19 = { 2001, 0.72687E-05}
3 def LAS_AMAT_20 = { 1200, -0.15317E+01}
3 def LAS_AMAT_21 = { 1110, -0.72434E-05}
3 def LAS_AMAT_22 = { 1101, 0.76453E-01}
3 def LAS_AMAT_23 = { 1020, 0.38152E-06}
3 def LAS_AMAT_24 = { 1011, 0.44978E-04}
3 def LAS_AMAT_25 = { 1002, 0.41038E-01}
3 def LAS_AMAT_26 = { 0300, 0.10045E+04}
3 def LAS_AMAT_27 = { 0210, 0.33079E-02}
3 def LAS_AMAT_28 = { 0201, -0.76794E+02}
3 def LAS_AMAT_29 = { 0120, -0.94281E-03}
3 def LAS_AMAT_30 = { 0111, -0.38992E-01}
3 def LAS_AMAT_31 = { 0102, -0.10189E+01}
3 def LAS_AMAT_32 = { 0021, -0.62901E-03}
3 def LAS_AMAT_33 = { 0012, -0.92175E-01}
3 def LAS_AMAT_34 = { 0003, -0.21562E+03}
3 def LAS_AMAT_35 = { 4000, -0.13113E-10}
3 def LAS_AMAT_36 = { 0400, -0.18417E+03}
3 def LAS_AMAT_37 = { 0040, -0.23016E-07}
3 def LAS_AMAT_38 = { 0004, -0.30098E+03}
3 def LAS_AMAT_39 = { 3100, 0.97679E-07}
3 def LAS_AMAT_40 = { 2200, -0.27604E-03}
3 def LAS_AMAT_41 = { 1300, 0.38034E+00}
2. com Reaction
3 com Run Dependent Default Value
4 def RUN_dEX = 0.000 # Excitation Energy Shift
4 def RUN_dYC = 0.000 # YC Shift (mm)
4 def RUN_dPHC = 0.000 # PHC Shift (rad)
4 def RUN_LAS_YG = 0.000 # LAS Y Central Position for (-50<=LAS_X<=50mm)
4 def M2 = M_12C # Target Mass [MeV/c2]
4 def GR_ANG = 0.0
3 com Common
4 com Projectile Kinetic Energy [MeV]
5 def T1 = 65
4 com Projectile Mass [MeV/c2]
5 def M1 = M_p
4 com Scattered Particle Mass [MeV/c2]
5 def GR_M3 = M_p
4 com Recoil Particle Mass (g.s.) [MeV/c2]
5 def GR_M4 = M2
4 com Charge of the particle [eC]
5 def GR_q = 1
5 def LAS_q = 1
4 com Spectrometer Bending Radius [m]
5 def GR_Rho = 3.00
5 def LAS_Rho = 1.75
4 com Residual Particle Mass for GR+LAS [MeV/c2]
5 def M5 = 0.0
4 com Energy loss of the ejectile in the target [MeV]
5 def GR_DE = 0
3 com faint
4 com Target Mass [MeV/c2]
5 def M2 = M_p
4 com Magnetic field of the Dipoles [mT]
5 def GR_B = 868.3 # 0deg for oct06
-3. com p+12C inelastic
4 com Target Mass [MeV/c2]
5 def M2 = M_12C
4 com Magnetic field of the Dipoles [mT]
5 def GR_B = 868.3 # 0deg for oct06
-3. com p+27Al inelastic
4 com Target Mass [MeV/c2]
5 def M2 = 25133.150 # 27Al mass MeV/c^2
-3. com p+58Ni inelastic
4 com Target Mass [MeV/c2]
5 def M2 = M_58Ni
4 com Magnetic field of the Dipoles [mT]
5 def GR_B = 871.5 # 0deg for oct06
-3 com p+120Sn inelastic
-3. com a+90Zr inelastic
5 com Magnetic field of the Dipoles [mT]
6 def GR_B = 950.621 # 2.8deg 90Zr
4. com Target Mass [MeV/c2]
5 def M2 = 83745.7 # 90Zr mass MeV/c^2
-3. com p+154Sm inelastic
4 com Target Mass [MeV/c2]
5 def M2 = 143377.660 # 154Sm mass MeV/c^2
4 com Magnetic field of the Dipoles [mT]
-3. com p+197Au elastic
4 com Projectile Kinetic Energy [MeV]
5 def T1 = 301.84
4 com Energy loss of the ejectile in the target [MeV]
5 def GR_DE = 0
4 com Target Mass [MeV/c2]
5 def M2 = 197*AMU
4 com Spectrometer angle
5 def GR_ANG = 8.0
5 def LAS_ANG = 33.7
4 com Magnetic field of the Dipoles [mT]
5 def GR_B = 889.187
5 def LAS_B = 1511.325
-3. com p+208Pb inelastic
4 com Target Mass [MeV/c2]
5 def M2 = M_208Pb
4 com Magnetic field of the Dipoles [mT]
1 com Definition of Variables
2. com Spin
3 def SPIN_UP = !(IPR&0x0002)
3 def SPIN_DN = !(IPR&0x0004)
3 def SPIN_DOWN = SPIN_DN
3 def SPIN = SPIN_UP+SPIN_DN*2
2. com BLP Mode
3 def BLP1_IN = !(IPR&0x0008)
3 def BLP2_IN = !(IPR&0x0010)
3 def BLP = BLP1_IN+BLP2_IN*2
2. com Event
3. com Bit
4 def EVENT_GR = bit(IPR,6-1)
4 def EVENT_LAS = bit(IPR,8-1)
4 def EVENT_COIN = bit(IPR,9-1)
4 def EVENT_2FPP = bit(IPR,10-1)
4 def EVENT_SLEVEL = bit(IPR,11-1)
4 def EVENT_FPP = bit(IPR,12-1)
4 def SLEVEL_ACC = bit(IPR,13-1)
4 def SLEVEL_REJ = bit(IPR,14-1)
4 def EVENT_BLP1 = bit(IPR,13-1)
4 def EVENT_BLP2 = bit(IPR,14-1)
4 def EVENT_BEND = bit(IPR,16-1)
3 def EVENT_GR_COIN = EVENT_GR || EVENT_COIN
3 def EVENT_LAS_COIN = EVENT_LAS || EVENT_COIN
-3 def EVENT_GR_FPP = EVENT_GR_COIN && EVENT_FPP
-3 def EVENT_GR_FPP = EVENT_GR_COIN
3 def EVENT_1ST = EVENT_LAS_COIN && EVENT_2FPP
-3 def EVENT_2ND = EVENT_GR && EVENT_2FPP && EVENT_SLEVEL
-3 def EVENT_2ND_SAMP = EVENT_GR && EVENT_2FPP && EVENT_SLEVEL
3 def EVENT_2ND = EVENT_SLEVEL
3 def EVENT_2ND_SAMP = !EVENT_SLEVEL
3 def EVENT = EVENT_SLEVEL | (EVENT_FPP<<1) | (EVENT_GR<<2) | (EVENT_LAS<<3) | (EVENT_COIN<<4)
3 def SLEVEL = SLEVEL_ACC | (SLEVEL_REJ<<1)
2 com GR ADC
3 def GR_ADC[1] = GR_FERA_ADC[0x10]
3 def GR_ADC[2] = GR_FERA_ADC[0x11]
3 def GR_ADC[3] = GR_FERA_ADC[0x12]
3 def GR_ADC[4] = GR_FERA_ADC[0x13]
3 def GR_ADC[5] = GR_FERA_ADC[0x20]
3 def GR_ADC[6] = GR_FERA_ADC[0x28]
3 def GR_MADC[1] = sqrt(GR_ADC[1]*GR_ADC[2])
3 def GR_MADC[2] = sqrt(GR_ADC[3]*GR_ADC[4])
3 def GR_MADC[3] = sqrt(GR_ADC[5]*GR_ADC[6])
2. com GR TDC
-3 def GR_TDC[1] = GR_FERA_TDC[0x10]
3 def GR_TDC[1] = GR_FERA_TDC[0x16] # due to broke channel
3 def GR_TDC[2] = GR_FERA_TDC[0x11]
3 def GR_TDC[3] = GR_FERA_TDC[0x12]
3 def GR_TDC[4] = GR_FERA_TDC[0x13]
3 def GR_TDC[5] = GR_FERA_TDC[0x20]
3 def GR_TDC[6] = GR_FERA_TDC[0x28]
3 def GR_RF = GR_FERA_TDC[0x14]
3 def GR_TLAS = GR_FERA_TDC[0x15]
3 def GR_TPOS[1] = GR_TDC[2]-GR_TDC[1]
3 def GR_TPOS[2] = GR_TDC[4]-GR_TDC[3]
2. com GR Optics
3. com RFC
4 def GR_RFC = GR_RF[0]-0.067*GR_X[0]-31.0*GR_THC[0]*r2d-2.23e-3*GR_X[0]-1.67e-5*GR_X[0]*GR_X[0]
-4 def GR_RFC = GR_RF[0]-0.067*GR_X[0]-31.0*GR_THC[0]*r2d
-4 def GR_RFC = GR_RF[0]-0.067*GR_X[0]
3 com Ray
4 com No B.G. Correction
5 def GR_X = GR_XP[0]
5 def GR_Y = GR_YP[0]
5 def GR_TH = GR_RAY_TH[0]
5 def GR_PH = GR_RAY_PH[0]
3 com Particle momentum [MeV/c]
3 com Horizontal scattering angle [rad]
-4 def GR_THC = -0.4191*(GR_TH-8.770e-3-3.591e-5*GR_X)
4 def GR_THC = -0.4191*(GR_TH-6.770e-2-9.591e-5*GR_X)
-4 def GR_THC = GR_TH-3.591e-5*GR_X
4 def GR_EX_BI_CORR = {0.0} # Normal
4 def GR_EX_C = polynom(GR_EX[0], GR_EX_C_EX)+RUN_dEX # 06Nov
2. com LAS ADC
3 def LAS_ADC[1] = LAS_FERA_ADC[0x10]
3 def LAS_ADC[2] = LAS_FERA_ADC[0x11]
3 def LAS_ADC[3] = LAS_FERA_ADC[0x12]
3 def LAS_ADC[4] = LAS_FERA_ADC[0x13]
3 def LAS_ADC[5] = LAS_FERA_ADC[0x14]
3 def LAS_ADC[6] = LAS_FERA_ADC[0x15]
3 def LAS_ADC[7] = LAS_FERA_ADC[0x16]
3 def LAS_ADC[8] = LAS_FERA_ADC[0x17]
3 def LAS_ADC[9] = LAS_FERA_ADC[0x18]
3 def LAS_ADC[10] = LAS_FERA_ADC[0x19]
3 def LAS_ADC[11] = LAS_FERA_ADC[0x1a]
3 def LAS_ADC[12] = LAS_FERA_ADC[0x1b]
3 def LAS_MADC[1] = sqrt(LAS_ADC[1]*LAS_ADC[2])
3 def LAS_MADC[2] = sqrt(LAS_ADC[3]*LAS_ADC[4])
3 def LAS_MADC[3] = sqrt(LAS_ADC[5]*LAS_ADC[6])
3 def LAS_MADC[4] = sqrt(LAS_ADC[7]*LAS_ADC[8])
3 def LAS_MADC[5] = sqrt(LAS_ADC[9]*LAS_ADC[10])
3 def LAS_MADC[6] = sqrt(LAS_ADC[11]*LAS_ADC[12])
2. com LAS TDC
3 def LAS_TDC[1] = LAS_FERA_TDC[0x10]
3 def LAS_TDC[2] = LAS_FERA_TDC[0x12]
3 def LAS_TDC[3] = LAS_FERA_TDC[0x14]
3 def LAS_TDC[4] = LAS_FERA_TDC[0x16]
3 def LAS_TDC[5] = LAS_FERA_TDC[0x18]
3 def LAS_TDC[6] = LAS_FERA_TDC[0x1a]
3 def LAS_TDC[7] = LAS_FERA_TDC[0x1c]
3 def LAS_TDC[8] = LAS_FERA_TDC[0x1e]
3 def LAS_TDC[9] = LAS_FERA_TDC[0x20]
3 def LAS_TDC[10] = LAS_FERA_TDC[0x22]
3 def LAS_TDC[11] = LAS_FERA_TDC[0x24]
3 def LAS_TDC[12] = LAS_FERA_TDC[0x26]
3 def LAS_RF[1] = LAS_FERA_TDC[0x2a]
3 def LAS_RF[2] = LAS_FERA_TDC[0x2c]
3 def LAS_RF[3] = LAS_FERA_TDC[0x2e]
3 def LAS_TPOS[1] = LAS_TDC[2]-LAS_TDC[1]
3 def LAS_TPOS[2] = LAS_TDC[4]-LAS_TDC[3]
3 def LAS_TPOS[3] = LAS_TDC[6]-LAS_TDC[5]
3 def LAS_TPOS[4] = LAS_TDC[8]-LAS_TDC[7]
3 def LAS_TPOS[5] = LAS_TDC[10]-LAS_TDC[9]
3 def LAS_TPOS[6] = LAS_TDC[12]-LAS_TDC[11]
2. com EVENT LAS
3 def EVENT_LAS_TRUE = LAS_PI_P && LAS_RAYID[0]==0 && LAS_CHI2 < 10
-3 def EVENT_LAS_TRUE = LAS_PI_P && LAS_RAYID[0]==0
2. com LAS Optics
3 com Particle momentum [MeV/c]
4 def LAS_P = (LAS_DP+1)*0.3*LAS_B*LAS_Rho*LAS_q
3 com Scattering angle [red]
4 def LAS_TH3 = LAS_ANG*d2r+LAS_AI
4 def LAS_EXC = LAS_EX
2. com GR-FPP Hodoscope ADC
3 def GF_ADC_XU[1] = GR_FERA_ADC[0x20]
3 def GF_ADC_XU[2] = GR_FERA_ADC[0x21]
3 def GF_ADC_XU[3] = GR_FERA_ADC[0x22]
3 def GF_ADC_XU[4] = GR_FERA_ADC[0x23]
3 def GF_ADC_XU[5] = GR_FERA_ADC[0x24]
3 def GF_ADC_XU[6] = GR_FERA_ADC[0x25]
3 def GF_ADC_XU[7] = GR_FERA_ADC[0x26]
3 def GF_ADC_XU[8] = GR_FERA_ADC[0x27]
3 def GF_ADC_XD[1] = GR_FERA_ADC[0x28]
3 def GF_ADC_XD[2] = GR_FERA_ADC[0x29]
3 def GF_ADC_XD[3] = GR_FERA_ADC[0x2a]
3 def GF_ADC_XD[4] = GR_FERA_ADC[0x2b]
3 def GF_ADC_XD[5] = GR_FERA_ADC[0x2c]
3 def GF_ADC_XD[6] = GR_FERA_ADC[0x2d]
3 def GF_ADC_XD[7] = GR_FERA_ADC[0x2e]
3 def GF_ADC_XD[8] = GR_FERA_ADC[0x2f]
3 def GF_ADC_YL[1] = GR_FERA_ADC[0x16]
3 def GF_ADC_YL[2] = GR_FERA_ADC[0x17]
3 def GF_ADC_YL[3] = GR_FERA_ADC[0x18]
3 def GF_ADC_YL[4] = GR_FERA_ADC[0x19]
3 def GF_ADC_YL[5] = GR_FERA_ADC[0x1a]
3 def GF_ADC_YR[1] = GR_FERA_ADC[0x1b]
3 def GF_ADC_YR[2] = GR_FERA_ADC[0x1c]
3 def GF_ADC_YR[3] = GR_FERA_ADC[0x1d]
3 def GF_ADC_YR[4] = GR_FERA_ADC[0x1e]
3 def GF_ADC_YR[5] = GR_FERA_ADC[0x1f]
3 def GF_MADC_X[1] = sqrt(GF_ADC_XU[1]*GF_ADC_XD[1])
3 def GF_MADC_X[2] = sqrt(GF_ADC_XU[2]*GF_ADC_XD[2])
3 def GF_MADC_X[3] = sqrt(GF_ADC_XU[3]*GF_ADC_XD[3])
3 def GF_MADC_X[4] = sqrt(GF_ADC_XU[4]*GF_ADC_XD[4])
3 def GF_MADC_X[5] = sqrt(GF_ADC_XU[5]*GF_ADC_XD[5])
3 def GF_MADC_X[6] = sqrt(GF_ADC_XU[6]*GF_ADC_XD[6])
3 def GF_MADC_X[7] = sqrt(GF_ADC_XU[7]*GF_ADC_XD[7])
3 def GF_MADC_X[8] = sqrt(GF_ADC_XU[8]*GF_ADC_XD[8])
3 def GF_MADC_Y[1] = sqrt(GF_ADC_YL[1]*GF_ADC_YR[1])
3 def GF_MADC_Y[2] = sqrt(GF_ADC_YL[2]*GF_ADC_YR[2])
3 def GF_MADC_Y[3] = sqrt(GF_ADC_YL[3]*GF_ADC_YR[3])
3 def GF_MADC_Y[4] = sqrt(GF_ADC_YL[4]*GF_ADC_YR[4])
3 def GF_MADC_Y[5] = sqrt(GF_ADC_YL[5]*GF_ADC_YR[5])
2. com GR-FPP Hodoscope TDC
3 def GF_TDC_XU[1] = GR_FERA_TDC[0x20]
3 def GF_TDC_XU[2] = GR_FERA_TDC[0x21]
3 def GF_TDC_XU[3] = GR_FERA_TDC[0x22]
3 def GF_TDC_XU[4] = GR_FERA_TDC[0x23]
3 def GF_TDC_XU[5] = GR_FERA_TDC[0x24]
3 def GF_TDC_XU[6] = GR_FERA_TDC[0x25]
3 def GF_TDC_XU[7] = GR_FERA_TDC[0x26]
3 def GF_TDC_XU[8] = GR_FERA_TDC[0x27]
3 def GF_TDC_XD[1] = GR_FERA_TDC[0x28]
3 def GF_TDC_XD[2] = GR_FERA_TDC[0x29]
3 def GF_TDC_XD[3] = GR_FERA_TDC[0x2a]
3 def GF_TDC_XD[4] = GR_FERA_TDC[0x2b]
3 def GF_TDC_XD[5] = GR_FERA_TDC[0x2c]
3 def GF_TDC_XD[6] = GR_FERA_TDC[0x2d]
3 def GF_TDC_XD[7] = GR_FERA_TDC[0x2e]
3 def GF_TDC_XD[8] = GR_FERA_TDC[0x2f]
3 def GF_TDC_YL[1] = GR_FERA_TDC[0x16]
3 def GF_TDC_YL[2] = GR_FERA_TDC[0x17]
3 def GF_TDC_YL[3] = GR_FERA_TDC[0x18]
3 def GF_TDC_YL[4] = GR_FERA_TDC[0x19]
3 def GF_TDC_YL[5] = GR_FERA_TDC[0x1a]
3 def GF_TDC_YR[1] = GR_FERA_TDC[0x1b]
3 def GF_TDC_YR[2] = GR_FERA_TDC[0x1c]
3 def GF_TDC_YR[3] = GR_FERA_TDC[0x1d]
3 def GF_TDC_YR[4] = GR_FERA_TDC[0x1e]
3 def GF_TDC_YR[5] = GR_FERA_TDC[0x1f]
3 def GF_TPOS_X[1] = GF_TDC_XD[1]-GF_TDC_XU[1]
3 def GF_TPOS_X[2] = GF_TDC_XD[2]-GF_TDC_XU[2]
3 def GF_TPOS_X[3] = GF_TDC_XD[3]-GF_TDC_XU[3]
3 def GF_TPOS_X[4] = GF_TDC_XD[4]-GF_TDC_XU[4]
3 def GF_TPOS_X[5] = GF_TDC_XD[5]-GF_TDC_XU[5]
3 def GF_TPOS_X[6] = GF_TDC_XD[6]-GF_TDC_XU[6]
3 def GF_TPOS_X[7] = GF_TDC_XD[7]-GF_TDC_XU[7]
3 def GF_TPOS_X[8] = GF_TDC_XD[8]-GF_TDC_XU[8]
3 def GF_TPOS_Y[1] = GF_TDC_YR[1]-GF_TDC_YL[1]
3 def GF_TPOS_Y[2] = GF_TDC_YR[2]-GF_TDC_YL[2]
3 def GF_TPOS_Y[3] = GF_TDC_YR[3]-GF_TDC_YL[3]
3 def GF_TPOS_Y[4] = GF_TDC_YR[4]-GF_TDC_YL[4]
3 def GF_TPOS_Y[5] = GF_TDC_YR[5]-GF_TDC_YL[5]
2. com GR-FPP Optics
3 def GF_IN_X = GF_IN_XP[0]
3 def GF_IN_Y = GF_IN_YP[0]
3 def GF_IN_TH = GF_IN_RAY_TH[0]
3 def GF_IN_PH = GF_IN_RAY_PH[0]
3 def GF_OUT_X = GF_OUT_XP[0]
3 def GF_OUT_Y = GF_OUT_YP[0]
3 def GF_OUT_TH = GF_OUT_RAY_TH[0]
3 def GF_OUT_PH = GF_OUT_RAY_PH[0]
3 def GF_X = GF_OUT_X[0]
3 def GF_Y = GF_OUT_P[0]
2. com BLP
3. com BLP ADC
4 def BLP_ADC[1] = GR_ADC_OLD[0]
4 def BLP_ADC[2] = GR_ADC_OLD[1]
4 def BLP_ADC[3] = GR_ADC_OLD[2]
4 def BLP_ADC[4] = GR_ADC_OLD[3]
4 def BLP_ADC[5] = GR_ADC_OLD[4]
4 def BLP_ADC[6] = GR_ADC_OLD[5]
4 def BLP_ADC[7] = GR_ADC_OLD[6]
4 def BLP_ADC[8] = GR_ADC_OLD[7]
3. com BLP TDC
4 def BLP_TDC[1] = GR_TDC_OLD[0]
4 def BLP_TDC[2] = GR_TDC_OLD[1]
4 def BLP_TDC[3] = GR_TDC_OLD[2]
4 def BLP_TDC[4] = GR_TDC_OLD[3]
4 def BLP_TDC[5] = GR_TDC_OLD[4]
4 def BLP_TDC[6] = GR_TDC_OLD[5]
4 def BLP_TDC[7] = GR_TDC_OLD[6]
4 def BLP_TDC[8] = GR_TDC_OLD[7]
3. com BLP TDC True Range
-4 def BLP1_T = {2300,2700}
-4 def BLP2_T = {1600,1900}
4 def BLP1_T = {0,4096}
4 def BLP2_T = {0,4096}
1 com Definition of Gates
2 com GR Particle ID
3 def WIREX1_notgood = GR_PI_ADC && GR_RFGate
3 def GR_PI = GR_PI_ADC && GR_RFGate
-3 def GR_PI_ADC = gate(GR_MADC[1],102,227) && gate(GR_MADC[2],142,251)
3 def GR_PI_ADC = gate(GR_MADC[1],184,276) && gate(GR_MADC[2],194,302)
-3 def GR_RFGate = gate(GR_RF,526,658) || gate(GR_RF,1045,1205)
-3 def GR_RFGate = gate(GR_RF,454,512) # E421 Online test (1st Peak)
-3 def GR_RFGate = gate(GR_RF,512,555) # E421 Online test (2nd Peak)
-3 def GR_RFGate = gate(GR_RF,555,613) # E421 Online test (3rd Peak)
3 def GR_RFGate = gate(GR_RF,454,628) || gate(GR_RF,997,1184) # E421 Online until June5
-3 def GR_RFGate = gate(GR_RF,460,620) || gate(GR_RF,991,1159)
3 def GR_BLK = gate(BLOCK_N[0],0,30000)
-3 def GR_PI = TRUE
2. com GR Optics
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,-0.5, 0.5)
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,-0.6, 1.0)
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,0.80, 1.2)
3 def GR_TH_NAR = gate(GR_THC[0]*r2d,-0.5,0.5) # E421 and E422 for 0degree
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,-1.0,-0.5) # E421 and E422 for 0degree
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,0.0,1.0) # E421 and E422 for 2.5 and 4.5 angle
-3 def GR_TH_NAR = gate(GR_TH[0]*r2d,3.0,5.0) # E421 and E422 Au(p,p0) at 8deg
-3 def GR_TH_NAR = gate(GR_TH[0]*r2d,1.0, 3.0) # E421 and E422 for faint beam
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,0.1,0.4)
3 def GR_PH_NAR = gate(GR_PHC[0]*r2d,-0.25,0.25)
-3 def GR_TH_NAR1 = gate(GR_THC[0]*r2d,-0.2,0.2)
-3 def GR_TH_NAR2 = gate(GR_THC[0]*r2d,-0.6,-0.2)
-3 def GR_TH_NAR3 = gate(GR_THC[0]*r2d,-1.0,-0.6)
-3 def GR_TH_NAR4 = gate(GR_THC[0]*r2d,-1.4,-1.0)
-3 def GR_TH_NAR = gate(GR_THC[0]*r2d,-0.8,0.8)
-3. com GR_X Gate 08Nov
4 def G_X = gate(GR_X[0], 17, 46) # Run#8050 12C 15.1 MeV
-4 def G_X = gate(GR_X[0], -295, -273) # Run#8050 12C 7.6 MeV
-4 def G_X = gate(GR_X[0], -72, -64) # Run#8050 12C 12.7 MeV
-4 def G_X = gate(GR_X[0], 174, 256) # Run#8050 12C SDR
3 def G_ELASTIC = gate(GR_X[0],-500,-300)
3 def GR_EX_G = gate(GR_EX,-2,2)
3 def GR_CHI2_G = gate(GR_CHI2, 0, 50)
3 def GR_THC_G1 = gate(GR_THC[0]*r2d,-0.5,0.5)
3 def GR_Y_G = gate(GR_Y[0],-10.0,10.0)
-3 def GR_YC_G = TRUE # full open
3 def GR_YC_BGG = gate(GR_YC[0], -30.0, -10.0) # for calibrated YC, 2015 May E421 E422
-3 def GR_YC_BGGM = gate(GR_YC[0], -15.0, -5.0) # for calibrated YC, 2015 May E421 E422
3 def GR_YC_BGGM = gate(GR_YC[0], -14.0, -4.0) # for calibrated YC, 2015 May E421 E422, 96Zr
-3 def GR_YC_BGGM = gate(GR_YC[0], -12.0, -2.0) # for calibrated YC, 2015 May E421 E422, 94Zr
-3 def GR_YC_BGGP = gate(GR_YC[0], 5.0, 15.0) # for calibrated YC, 2015 May E421 E422
3 def GR_YC_BGGP = gate(GR_YC[0], 6.0, 16.0) # for calibrated YC, 2015 May E421 E422, 96Zr
-3 def GR_YC_BGGP = gate(GR_YC[0], 8.0, 18.0) # for calibrated YC, 2015 May E421 E422, 94Zr
-3 def GR_YC_G = gate(GR_YC[0], -5.0, 5.0) # for calibrated YC, 2015 May E421 E422
3 def GR_YC_G = gate(GR_YC[0], -4.0, 6.0) # for calibrated YC, 2015 May E421 E422, 96Zr
-3 def GR_YC_G = gate(GR_YC[0], -2.0, 8.0) # for calibrated YC, 2015 May E421 E422, 94Zr
3 def GR_YC_B = gate(GR_YC[0], -20, -10) || gate(GR_YC[0], 10, 20) # for calibrated YC
-3 def GR_TRUE = GR_YC_G && GR_PHC_G
3 def GR_TPOSP = gate(GR_TPOS[1],-50,-4)
3 def GR_TPOSC = GR_TPOS[2]-GR_XP[3]*0.045-14
3 def GR_TPOSC2 = GR_TPOS[2]+GR_POS_X1*0.3-42
3. com AI Gates
4 def GR_AI_G = GR_AI_G17
4 def GR_AI_G0 = gate(GR_AI[0]*r2d, -0.2, 0.2)
4 def GR_AI_G1 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G2 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G3 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G4 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G5 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G7 = gate(GR_AI[0]*r2d, -0.5, 0.0)
4 def GR_AI_G8 = gate(GR_AI[0]*r2d, 0.0, 0.5)
4 def GR_AI_G9 = gate(GR_AI[0]*r2d, -1.0,-0.5)
4 def GR_AI_G10 = gate(GR_AI[0]*r2d, 0.5, 1.0)
4 def GR_AI_G11 = gate(GR_AI[0]*r2d, -0.5, 0.0)
4 def GR_AI_G12 = gate(GR_AI[0]*r2d, 0.0, 0.5)
4 def GR_AI_G13 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G14 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G15 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G16 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G17 = gate(GR_AI[0]*r2d, -1.0, 1.0)
4 def GR_AI_G20 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G21 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G22 = gate(GR_AI[0]*r2d, -0.5, 0.5)
4 def GR_AI_G23 = gate(GR_AI[0]*r2d, -0.5, 0.5)
3. com BI Gates
4 def GR_BI_G = GR_BI_G17
4 def GR_BI_G0 = gate(GR_BI[0]*r2d, -0.5, 0.5)
4 def GR_BI_G1 = gate(GR_BI[0]*r2d, -2.5, -1.5)
4 def GR_BI_G2 = gate(GR_BI[0]*r2d, -1.5, -0.5)
4 def GR_BI_G3 = gate(GR_BI[0]*r2d, -0.5, 0.5)
4 def GR_BI_G4 = gate(GR_BI[0]*r2d, 0.5, 1.5)
4 def GR_BI_G5 = gate(GR_BI[0]*r2d, 1.5, 2.5)
4 def GR_BI_G7 = gate(GR_BI[0]*r2d, -2.0, 2.0)
4 def GR_BI_G8 = gate(GR_BI[0]*r2d, -2.0, 2.0)
4 def GR_BI_G9 = gate(GR_BI[0]*r2d, -1.0, 1.0)
4 def GR_BI_G10 = gate(GR_BI[0]*r2d, -1.0, 1.0)
4 def GR_BI_G12 = TRUE
4 def GR_BI_G13 = gate(GR_BI[0]*r2d, -0.5, 0.5)
4 def GR_BI_G14 = gate(GR_BI[0]*r2d, -2.0, 2.0)
4 def GR_BI_G15 = gate(GR_BI[0]*r2d, -1.0, 1.0)
4 def GR_BI_G16 = gate(GR_BI[0]*r2d, -2.5, 2.5)
4 def GR_BI_G17 = gate(GR_BI[0]*r2d, -2.5, 2.5)
4 def GR_BI_G20 = gate(abs(GR_BI[0]*r2d), 0.0, 0.8)
4 def GR_BI_G21 = gate(abs(GR_BI[0]*r2d), 0.8, 1.6)
4 def GR_BI_G22 = gate(abs(GR_BI[0]*r2d), 1.6, 2.4)
4 def GR_BI_G23 = gate(abs(GR_BI[0]*r2d), 0.0, 2.5)
2. com LAS Particle ID
3 def LAS_RF_G = gate(LAS_RF[2], 293, 488) || gate(LAS_RF[2], 604, 807)
3 def LAS_DE_G = gate(LAS_MADC[2], 120, 700) && gate(LAS_MADC[5], 100, 500)
-3 def LAS_PI = LAS_RF_G && LAS_DE_G
3 def LAS_PI = LAS_DE_G
-3 def LAS_PI = TRUE
2. com LAS Optics
3 def LAS_X = LAS_XP[0]
3 def LAS_Y = LAS_YP[0]
3 def LAS_TH = LAS_RAY_TH[0]
3 def LAS_PH = LAS_RAY_PH[0]
2. com GR-FPP
-3 def GF_ANA_G = GR_RAYID==0
-3 def GF_ANA_G = GR_YC_G
3 def GF_ANA_G = GR_YC_G && GR_PHC_G && GR_AI_G && GR_BI_G
3 com 2nd Scattering Angle
-4 def GF_SCAT_TH_MIN = 4.5
4 def GF_SCAT_TH_MIN = 5
4 def GF_SCAT_TH_MAX = 20
4 def GF_SCAT_PH_MAX = 66.8
4 def GF_SCAT_TH_G = gate(GF_SCAT_TH[0]*r2d, GF_SCAT_TH_MIN, GF_SCAT_TH_MAX)
4 def GF_SCAT_PH_L_G = gate(GF_SCAT_PH[0]*r2d, -GF_SCAT_PH_MAX, GF_SCAT_PH_MAX)
4 def GF_SCAT_PH_R_G = abs(GF_SCAT_PH[0]*r2d)>180-GF_SCAT_PH_MAX
4 def GF_SCAT_PH_U_G = gate(GF_SCAT_PH[0]*r2d, +90-GF_SCAT_PH_MAX, +90+GF_SCAT_PH_MAX)
4 def GF_SCAT_PH_D_G = gate(GF_SCAT_PH[0]*r2d, -90-GF_SCAT_PH_MAX, -90+GF_SCAT_PH_MAX)
4 def GF_SCAT_L_G = GF_SCAT_TH_G && GF_SCAT_PH_L_G
4 def GF_SCAT_R_G = GF_SCAT_TH_G && GF_SCAT_PH_R_G
4 def GF_SCAT_U_G = GF_SCAT_TH_G && GF_SCAT_PH_U_G
4 def GF_SCAT_D_G = GF_SCAT_TH_G && GF_SCAT_PH_D_G
-4 def GF_SCAT_G = GF_SCAT_L_G || GF_SCAT_R_G
4 def GF_SCAT_G = GF_SCAT_U_G || GF_SCAT_D_G
4 def GF_VX_G = gate(GF_VERTEX_X[0], -500, 600)
4 def GF_VY_G = gate(GF_VERTEX_Y[0], -100, 100)
4 def GF_VZ_G = gate(GF_VERTEX_Z[0], 1647, 1857)
4 def GF_VD_G = gate(GF_DIST[0], 0, 20)
4 def GF_VERTEX_G = GF_VX_G && GF_VY_G && GF_VZ_G && GF_VD_G
2 def EVENT_TRUE = TRUE
1 com E427 2014Nov Online
2 com Analysis Parameters
3 com Detector Configurations
4 com Drift time to distance
5 com GR
6 com G01 Jan-2016
7 def GR_WTDC_FILE = "gr_wtdc.hst" # GRAF
-7 def GR_WTDC_FILE = "gr_wtdcPAW.hst" # testing stuffs for the sake of the lol
-7 def GR_WTDC_FILE = "gr_wtdcROOT.hst" # testing stuffs for the sake of the lol
7 def GR_WTDC_X1 = {110,350}
7 def GR_WTDC_U1 = {108,352}
7 def GR_WTDC_X2 = {108,350}
7 def GR_WTDC_U2 = {110,347}
-4. com GR Chamber Configuration
5 def VDC_DX = +105.0 # E282 08nov (DSR+) offline calibration for Run#8129 (see analysis note#1 pp.92-)
5 def GR_VDC_TILT = {0,42.2245*d2r,0} # E282 08nov(DSR+) offline calibration for Run#8129
4. com GR Chamber Configuration (E387)
5 def VDC_DX = 0
5 def GR_VDC_TILT = {0,45*d2r,0}
4. com GR-FPP Chamber Configuration
5 com E282 08nov (DSR+, by offline calibration of Run#8129)
6 def GF3_DX = 357.69
6 def GF3_DY = 4.05 + 0.07
6 def GF4_DX = GF3_DX + 0.25 - 17.7
6 def GF4_DY = GF3_DY - 0.52 - 0.19
3. com GR Trace-Back Matrix
4 com Matrix for E316 09Nov (only the center, offsets, etc were adjsted based on the parameters of E282 06Nov)
5 com Matrix for A,B,Yc (for E282-06nov-mag2 070109, see analysis note p.44,58)
6. com Matrix for A (for E282-mag2 070109#1 = Standard)
7 def GR_AMAT_01 = { 0000, 1.61034e-02}
7 def GR_AMAT_02 = { 0100, -4.17069e-01}
7 def GR_AMAT_03 = { 1000, 2.36614e-05}
7 def GR_AMAT_04 = { 0000, -0.0235} # offset from E282 06nov
6. com Matrix for B (for E282-mag2 070109#2 = Standard)
7 def b_factor = 0.850
7 def GR_BMAT_01 = { 0000, -9.65234e-03}
7 def GR_BMAT_02 = { 0000, -1.07086e-03*RUN_LAS_YG}
7 def GR_BMAT_03 = { 0001, 4.17063e+00}
7 def GR_BMAT_04 = { 0010, -1.99037e-03*b_factor}
7 def GR_BMAT_05 = { 1000, 6.16186e-06}
7 def GR_BMAT_06 = { 1000, 6.17463e-07*RUN_LAS_YG}
7 def GR_BMAT_07 = { 0000, +(2.15e-04*RUN_LAS_YG)} # calculated ideal slope, see analysis note p.43
7 def GR_BMAT_08 = { 0000, -0.005}
-6. com Matrix for YC (for E282-mag2 070109#3 = Standard)
7 def GR_YMAT_01 = { 0010, 1.00000e+00}
7 def GR_YMAT_02 = { 0000, -3.74081e-01}
7 def GR_YMAT_03 = { 0000, 1.05732e+00*RUN_LAS_YG}
7 def GR_YMAT_04 = { 0001, 7.17335e+02}
7 def GR_YMAT_05 = { 0100, -5.20327e+01}
7 def GR_YMAT_06 = { 0101, 9.00266e+03}
7 def GR_YMAT_07 = { 1000, 1.24845e-02}
7 def GR_YMAT_08 = { 1001, -3.22387e+00}
7 def GR_YMAT_09 = { 1101, -1.37672e+01}
-6. com Matrix for YC (for E316-09nov Roughly Adjusted based on E282-08nov, 20121111)
7 def GR_YMAT_01 = { 0010, 1.00000e+00}
7 def GR_YMAT_02 = { 0000, 1.54185e+00}
7 def GR_YMAT_03 = { 0001, 1.36192e+03}
7 def GR_YMAT_04 = { 0100, -1.51211e+01}
7 def GR_YMAT_05 = { 0101, 4.43422e+03}
7 def GR_YMAT_06 = { 1000, 6.19236e-03}
7 def GR_YMAT_07 = { 1001, -1.34048e+00}
7 def GR_YMAT_08 = { 1101, -6.95930e+00}
7 def GR_YMAT_09 = { 0000, 5.61352e-01*RUN_LAS_YG}
7 com Matrix for YC (displacement, adjusted for Run#9063)
8 def GR_YMAT_10 = { 0000, -7.78+RUN_dYC }
8 def GR_YMAT_11 = { 0001, 0.600e+03}
8 def GR_YMAT_12 = { 1001, -1.2}
8 def GR_YMAT_12 = { 1000, 0.005}
-6. com Matrix for YC (for E316-09nov from Run9049 2nd order for X,TH 20121113)
7 def GR_YMAT_01 = { 0010, 1.00000e+00}
7 def GR_YMAT_02 = { 0000, -5.85904e+00}
7 def GR_YMAT_03 = { 0001, 2.18630e+03}
7 def GR_YMAT_04 = { 0100, -3.95577e+01}
7 def GR_YMAT_05 = { 0101, 1.13246e+04}
7 def GR_YMAT_06 = { 0200, 8.46496e+01}
7 def GR_YMAT_07 = { 0201, -4.04091e+04}
7 def GR_YMAT_08 = { 1000, 1.69274e-02}
7 def GR_YMAT_09 = { 1001, -4.79066e+00}
7 def GR_YMAT_10 = { 1101, 2.45517e+01}
7 def GR_YMAT_11 = { 1200, -2.15860e+00}
7 def GR_YMAT_12 = { 1201, 1.20621e+03}
7 def GR_YMAT_13 = { 2000, -5.58063e-06}
7 def GR_YMAT_14 = { 2001, 4.67959e-03}
7 def GR_YMAT_15 = { 2100, 3.07743e-04}
7 def GR_YMAT_16 = { 2101, -1.08785e-01}
7 def GR_YMAT_17 = { 2200, 6.13039e-03}
7 def GR_YMAT_18 = { 2201, -1.22234e+00}
7 def GR_YMAT_19 = { 0000, 7.931}
-7 def GR_YMAT_09 = { 0000, 5.61352e-01*RUN_LAS_YG}
7 com Matrix for YC (displacement, adjusted for Run#9063)
8 def GR_YMAT_10 = { 0000, -7.78+RUN_dYC }
-6. com Matrix for YC (for E316-09nov from Run9049 2nd order for X 20121114)
7 def GR_YMAT_01 = { 0010, 1.00000e+00}
7 def GR_YMAT_02 = { 0000, -5.92181e+00}
7 def GR_YMAT_03 = { 0001, 2.20656e+03}
7 def GR_YMAT_04 = { 0100, -4.56032e+01}
7 def GR_YMAT_05 = { 0101, 1.38724e+04}
7 def GR_YMAT_06 = { 1000, 1.80252e-02}
7 def GR_YMAT_07 = { 1001, -5.22284e+00}
7 def GR_YMAT_08 = { 1100, 1.36728e-01}
7 def GR_YMAT_09 = { 1101, -4.25944e+01}
7 def GR_YMAT_10 = { 2000, -8.24725e-06}
7 def GR_YMAT_11 = { 2001, 5.80501e-03}
7 def GR_YMAT_12 = { 2100, -6.56356e-05}
7 def GR_YMAT_13 = { 2101, 2.41542e-04}
-7 def GR_YMAT_14 = { 0000, 5.61352e-01*RUN_LAS_YG}
6. com Matrix for YC (for E316-09nov from Run9049 1st order for all 20121114)
7 def GR_YMAT_01 = { 0010, 1.00000e+00}
7 def GR_YMAT_02 = { 0000, -5.95045e+00}
7 def GR_YMAT_03 = { 0001, 2.36363e+03}
7 def GR_YMAT_04 = { 0100, -4.08823e+01}
7 def GR_YMAT_05 = { 0101, 1.25884e+04}
7 def GR_YMAT_06 = { 1000, 1.46217e-02}
7 def GR_YMAT_07 = { 1001, -3.75711e+00}
7 def GR_YMAT_08 = { 1100, 1.01238e-01}
7 def GR_YMAT_09 = { 1101, -2.26667e+01}
-7 def GR_YMAT_10 = { 0000, 5.61352e-01*RUN_LAS_YG}
5 com Matrix for Xc by Tamii 01-JUNE-2006 with kinematical correction
-6. com no correction
7 def GR_XMAT_01 = { 1000, 1.000e-00 }
-6. com xc-th by Tamii for oct04 on 05-JUNE-2006 with kin. corr. (X:3,TH:4) up to HighEx
7 com Fitted from xc_ai_calib_2.dat
7 def GR_XMAT_01 = { 0000, 1.078e+01 }