Pluto v6.3.0
Loading...
Searching...
No Matches
PPionBeamPlugin.cc
Go to the documentation of this file.
1
2//Plugin for the PPionBeamAmplitude
3//
4//This plugin enables the amplitudes from Ref. 23 for the
5//reaction pi + N -> N e+ e-
6//
7//The energy-dependent matrix elements have been supported by
8//the authors of Ref. 23
9//
10//To use this plugin, type:
11//makeDistributionManager()->Exec("pion_beam");
12//A flat di-electron generator is automatially employed
13//Therefore, the weights of the particles have to be used.
14//
15//
16// Author: Schuldes/Froehlich
17// Written: Jul 2009
19
20#include "PPionBeamPlugin.h"
21
22
23PPionBeamPlugin::PPionBeamPlugin(const Char_t *id, const Char_t *de):
25
26 Rho0Propagator = NULL;
27 OmegaPropagator = NULL;
28
29 Pi_minusBeamAmplitude = NULL;
30 Pi_minusBeamAmplitude_gen = NULL;
31
32 Pi_plusBeamAmplitude = NULL;
33 Pi_plusBeamAmplitude_gen = NULL;
34
35 // RequiresPlugin("elementary");
36}
37
39
41
42 if (!makeStaticData()->IsParticleValid("pi- + p")) {
43 makeStaticData()->AddParticle(14009, "pi- + p", 1.1);
44 }
45 if (!makeStaticData()->IsParticleValid("pi+ + n")) {
46 makeStaticData()->AddParticle(13008, "pi+ + n", 1.1);
47 }
48
49
50 Int_t ipid[3];
51 ipid[0] = makeStaticData()->GetParticleID("pi- + p");
52 ipid[1] = makeStaticData()->GetParticleID("n");
53 ipid[2] = makeStaticData()->GetParticleID("dilepton");
54
55
56 //Add decay
57 if (makeStaticData()->GetDecayKey(ipid, 2) < 0)
58 makeStaticData()->AddDecay(-1, "pi- + p -> n + dilepton",
59 "pi- + p","n,dilepton", 1.0 );
60
61 ipid[0] = makeStaticData()->GetParticleID("pi+ + n");
62 ipid[1] = makeStaticData()->GetParticleID("p");
63 ipid[2] = makeStaticData()->GetParticleID("dilepton");
64
65
66 if (makeStaticData()->GetDecayKey(ipid, 2) < 0)
67 makeStaticData()->AddDecay(-1, "pi+ + n -> p + dilepton",
68 "pi+ + n","p,dilepton", 1.0 );
69
70
71 pdmutil->AddSubGroup("pion_beam", "Elementary pion beam reactions", "root");
72 pdmutil->SetGroup("pion_beam");
73
74
75 Pi_minusBeamAmplitude =
76 new PPionBeamAmplitude("Pi_minusBeamAmplitude@pi- + p_to_n_dilepton",
77 "Model from Lutz/Soyeur/Friman", -1);
78
79 Pi_minusBeamAmplitude->Add("pi-,grandparent,pion");
80 Pi_minusBeamAmplitude->Add("p,grandparent,p_in");
81 Pi_minusBeamAmplitude->Add("dilepton,daughter");
82 Pi_minusBeamAmplitude->Add("n,daughter,p_out");
83 Pi_minusBeamAmplitude->Add("q,parent");
84
85 PPionBeamAmplitude *Pi_plusBeamAmplitude =
86 new PPionBeamAmplitude("Pi_plusBeamAmplitude@pi+ + n_to_p_dilepton",
87 "Model from Lutz/Soyeur/Friman", -1);
88
89 Pi_plusBeamAmplitude->Add("pi+,grandparent,pion");
90 Pi_plusBeamAmplitude->Add("n,grandparent,p_in");
91 Pi_plusBeamAmplitude->Add("dilepton,daughter");
92 Pi_plusBeamAmplitude->Add("p,daughter,p_out");
93 Pi_plusBeamAmplitude->Add("q,parent");
94
95 TF1 *flat = new TF1("pion_flat", "1.", 0, 0.99);
96 //The "PInclusiveModel" can be used as a generator:
97
98 Pi_minusBeamAmplitude_gen =
99 new PInclusiveModel("pi-_flat@pi- + p_brems_n_dilepton/generator", "Dilepton generator", -1);
100
101 //The distribution template:
102 Pi_minusBeamAmplitude_gen->Add("q,parent");
103 Pi_minusBeamAmplitude_gen->Add("n,daughter");
104 Pi_minusBeamAmplitude_gen->Add("dilepton,daughter,primary");
105 Pi_minusBeamAmplitude_gen->SetSampleFunction(flat);
106 //Enable distribution as a generator
107 Pi_minusBeamAmplitude_gen->EnableGenerator();
108 pdmutil->Add(Pi_minusBeamAmplitude_gen);
109
110 Pi_plusBeamAmplitude_gen =
111 new PInclusiveModel("pi+_flat@pi+ + n_brems_p_dilepton/generator", "Dilepton generator", -1);
112
113 //The distribution template:
114 Pi_plusBeamAmplitude_gen->Add("q,parent");
115 Pi_plusBeamAmplitude_gen->Add("p,daughter");
116 Pi_plusBeamAmplitude_gen->Add("dilepton,daughter,primary");
117 Pi_plusBeamAmplitude_gen->SetSampleFunction(flat);
118 //Enable distribution as a generator
119 Pi_plusBeamAmplitude_gen->EnableGenerator();
120 pdmutil->Add(Pi_plusBeamAmplitude_gen);
121
122
123 double Graph_Rho_x[150] = { 1.200010,
124 1.205343,
125 1.210677,
126 1.216010,
127 1.221343,
128 1.226677,
129 1.232010,
130 1.237343,
131 1.242677,
132 1.248010,
133 1.253343,
134 1.258677,
135 1.264010,
136 1.269343,
137 1.274677,
138 1.280010,
139 1.285343,
140 1.290677,
141 1.296010,
142 1.301343,
143 1.306677,
144 1.312010,
145 1.317343,
146 1.322677,
147 1.328010,
148 1.333343,
149 1.338677,
150 1.344010,
151 1.349343,
152 1.354677,
153 1.360010,
154 1.365343,
155 1.370677,
156 1.376010,
157 1.381343,
158 1.386677,
159 1.392010,
160 1.397343,
161 1.402677,
162 1.408010,
163 1.413343,
164 1.418677,
165 1.424010,
166 1.429343,
167 1.434677,
168 1.440010,
169 1.445343,
170 1.450677,
171 1.456010,
172 1.461343,
173 1.466677,
174 1.472010,
175 1.477343,
176 1.482677,
177 1.488010,
178 1.493343,
179 1.498677,
180 1.504010,
181 1.509343,
182 1.514677,
183 1.520010,
184 1.525343,
185 1.530677,
186 1.536010,
187 1.541343,
188 1.546677,
189 1.552010,
190 1.557343,
191 1.562677,
192 1.568010,
193 1.573343,
194 1.578677,
195 1.584010,
196 1.589343,
197 1.594677,
198 1.600010,
199 1.605343,
200 1.610677,
201 1.616010,
202 1.621343,
203 1.626677,
204 1.632010,
205 1.637343,
206 1.642677,
207 1.648010,
208 1.653343,
209 1.658677,
210 1.664010,
211 1.669343,
212 1.674677,
213 1.680010,
214 1.685343,
215 1.690677,
216 1.696010,
217 1.701343,
218 1.706677,
219 1.712010,
220 1.717343,
221 1.722677,
222 1.728010,
223 1.733343,
224 1.738677,
225 1.744010,
226 1.749343,
227 1.754677,
228 1.760010,
229 1.765343,
230 1.770677,
231 1.776010,
232 1.781343,
233 1.786677,
234 1.792010,
235 1.797343,
236 1.802677,
237 1.808010,
238 1.813343,
239 1.818677,
240 1.824010,
241 1.829343,
242 1.834677,
243 1.840010,
244 1.845343,
245 1.850677,
246 1.856010,
247 1.861343,
248 1.866677,
249 1.872010,
250 1.877343,
251 1.882677,
252 1.888010,
253 1.893343,
254 1.898677,
255 1.904010,
256 1.909343,
257 1.914677,
258 1.920010,
259 1.925343,
260 1.930677,
261 1.936010,
262 1.941343,
263 1.946677,
264 1.952010,
265 1.957343,
266 1.962677,
267 1.968010,
268 1.973343,
269 1.978677,
270 1.984010,
271 1.989343,
272 1.994677
273 };
274 double Graph_Rho_y_Re_T11[150] = {10.927144,
275 11.054252,
276 11.183319,
277 11.314396,
278 11.447535,
279 11.582795,
280 11.720238,
281 11.859932,
282 12.001948,
283 12.146364,
284 12.293266,
285 12.442745,
286 12.594900,
287 12.749843,
288 12.907689,
289 13.068568,
290 13.232624,
291 13.400007,
292 13.570890,
293 13.745453,
294 13.923903,
295 14.106458,
296 14.293363,
297 14.484882,
298 14.681305,
299 14.882954,
300 15.090174,
301 15.303348,
302 15.522894,
303 15.749271,
304 15.982983,
305 16.224578,
306 16.474674,
307 16.733924,
308 17.003239,
309 17.283150,
310 17.574746,
311 17.879095,
312 18.197160,
313 18.530794,
314 18.881362,
315 19.250767,
316 19.641255,
317 20.055519,
318 20.496921,
319 20.969657,
320 21.479111,
321 22.032411,
322 22.639301,
323 23.313768,
324 24.077115,
325 24.964850,
326 26.046298,
327 27.514324,
328 23.429189,
329 17.922962,
330 14.284720,
331 11.418496,
332 9.060166,
333 7.095792,
334 5.459748,
335 4.106839,
336 3.001836,
337 2.115744,
338 1.424028,
339 0.905675,
340 0.542946,
341 0.320928,
342 0.227112,
343 0.251003,
344 0.383341,
345 0.615095,
346 0.936011,
347 1.331964,
348 1.780537,
349 2.241361,
350 2.625383,
351 2.419523,
352 2.773354,
353 2.867063,
354 2.729286,
355 2.428496,
356 2.038813,
357 1.616750,
358 1.198547,
359 0.805547,
360 0.449584,
361 0.137842,
362 -0.124826,
363 -0.334509,
364 -0.488020,
365 -0.569560,
366 -0.491498,
367 -0.384744,
368 -0.250865,
369 -0.088530,
370 0.109346,
371 0.357700,
372 0.607211,
373 0.710514,
374 0.734399,
375 0.729459,
376 0.711516,
377 0.686936,
378 0.659246,
379 0.629241,
380 0.597348,
381 0.565423,
382 0.534413,
383 0.504613,
384 0.476116,
385 0.449040,
386 0.423345,
387 0.399003,
388 0.376021,
389 0.354327,
390 0.333904,
391 0.314705,
392 0.296628,
393 0.279665,
394 0.263739,
395 0.248807,
396 0.234862,
397 0.221789,
398 0.209578,
399 0.198175,
400 0.187539,
401 0.177665,
402 0.168470,
403 0.159933,
404 0.152025,
405 0.144707,
406 0.137973,
407 0.131771,
408 0.126073,
409 0.120865,
410 0.116114,
411 0.111810,
412 0.107930,
413 0.104436,
414 0.101328,
415 0.098576,
416 0.096168,
417 0.094107,
418 0.092346,
419 0.090888,
420 0.089713,
421 0.088812,
422 0.088187,
423 0.087804 };
424
425 double Graph_Rho_y_Im_T11[150] = {0.910200,
426 0.954286,
427 0.999784,
428 1.046767,
429 1.095314,
430 1.145507,
431 1.197431,
432 1.251175,
433 1.306835,
434 1.364512,
435 1.424314,
436 1.486355,
437 1.550759,
438 1.617661,
439 1.687204,
440 1.759545,
441 1.834858,
442 1.913330,
443 1.995166,
444 2.080593,
445 2.169870,
446 2.263268,
447 2.361107,
448 2.463729,
449 2.571525,
450 2.684937,
451 2.804456,
452 2.930638,
453 3.064111,
454 3.205591,
455 3.355889,
456 3.515921,
457 3.686747,
458 3.869556,
459 4.065889,
460 4.277049,
461 4.505037,
462 4.752027,
463 5.020286,
464 5.313351,
465 5.634549,
466 5.988255,
467 6.379784,
468 6.815732,
469 7.304491,
470 7.856877,
471 8.487220,
472 9.215156,
473 10.068662,
474 11.089791,
475 12.346527,
476 13.960898,
477 16.194830,
478 19.872674,
479 28.166761,
480 28.334894,
481 27.899177,
482 27.148490,
483 26.189925,
484 25.090924,
485 23.901309,
486 22.659430,
487 21.395293,
488 20.132271,
489 18.888536,
490 17.678337,
491 16.513072,
492 15.402340,
493 14.355056,
494 13.380355,
495 12.488697,
496 11.692907,
497 11.009356,
498 10.459301,
499 10.070405,
500 9.879183,
501 9.937317,
502 10.098127,
503 9.923171,
504 10.131774,
505 10.360310,
506 10.505713,
507 10.533933,
508 10.447338,
509 10.262323,
510 9.997361,
511 9.669314,
512 9.290949,
513 8.873061,
514 8.423727,
515 7.948418,
516 7.407230,
517 6.941536,
518 6.545564,
519 6.197589,
520 5.892170,
521 5.632294,
522 5.444328,
523 5.414461,
524 5.460242,
525 5.471005,
526 5.452900,
527 5.417032,
528 5.370151,
529 5.315913,
530 5.258264,
531 5.196802,
532 5.132822,
533 5.067713,
534 5.002415,
535 4.937646,
536 4.873792,
537 4.811268,
538 4.750327,
539 4.691095,
540 4.633722,
541 4.578238,
542 4.524679,
543 4.473112,
544 4.423478,
545 4.375775,
546 4.329958,
547 4.285930,
548 4.243713,
549 4.203203,
550 4.164347,
551 4.127078,
552 4.091288,
553 4.056955,
554 4.023991,
555 3.992322,
556 3.961887,
557 3.932592,
558 3.904393,
559 3.877226,
560 3.851016,
561 3.825713,
562 3.801244,
563 3.777555,
564 3.754608,
565 3.732333,
566 3.710694,
567 3.689639,
568 3.669109,
569 3.649084,
570 3.629514,
571 3.610366,
572 3.591602,
573 3.573178,
574 3.555080
575 };
576
577 double Graph_Rho_y_Re_T13[150] = {79.999754,
578 81.202991,
579 82.440358,
580 83.713773,
581 85.025272,
582 86.377041,
583 87.771493,
584 89.211189,
585 90.698927,
586 92.237830,
587 93.831169,
588 95.482617,
589 97.196182,
590 98.976049,
591 100.827450,
592 102.755150,
593 104.765348,
594 106.864575,
595 109.059305,
596 111.358697,
597 113.770248,
598 116.305084,
599 118.973334,
600 121.787054,
601 124.760419,
602 127.907485,
603 131.245170,
604 134.791346,
605 138.565223,
606 142.588566,
607 146.884262,
608 151.478873,
609 156.398289,
610 161.676925,
611 167.341747,
612 173.434967,
613 179.996182,
614 187.063240,
615 194.688249,
616 202.906814,
617 211.768987,
618 221.314068,
619 231.557608,
620 242.501176,
621 254.099687,
622 266.228576,
623 278.660484,
624 290.947852,
625 302.405098,
626 311.914614,
627 317.814662,
628 317.820817,
629 309.052088,
630 288.515445,
631 254.048605,
632 205.724301,
633 146.838922,
634 83.613864,
635 23.084481,
636 -29.164873,
637 -70.385655,
638 -100.443169,
639 -120.732715,
640 -133.264735,
641 -140.019930,
642 -142.638508,
643 -142.414483,
644 -140.284675,
645 -136.924826,
646 -132.806143,
647 -128.252667,
648 -123.487439,
649 -118.658135,
650 -113.863949,
651 -109.167613,
652 -104.608551,
653 -100.210174,
654 -95.983230,
655 -91.931735,
656 -88.053837,
657 -84.344307,
658 -80.796371,
659 -77.400747,
660 -74.148347,
661 -71.027939,
662 -68.029386,
663 -65.142647,
664 -62.353245,
665 -59.652588,
666 -57.027252,
667 -54.463349,
668 -51.945514,
669 -49.455479,
670 -46.968581,
671 -44.449700,
672 -41.840615,
673 -39.030936,
674 -35.786804,
675 -32.178549,
676 -29.528051,
677 -27.653129,
678 -26.154901,
679 -24.882372,
680 -23.767237,
681 -22.768257,
682 -21.868339,
683 -21.064489,
684 -20.333940,
685 -19.661180,
686 -19.037930,
687 -18.458681,
688 -17.918412,
689 -17.413810,
690 -16.941920,
691 -16.499816,
692 -16.085387,
693 -15.696490,
694 -15.331194,
695 -14.988424,
696 -14.666358,
697 -14.363838,
698 -14.079714,
699 -13.812463,
700 -13.561723,
701 -13.326213,
702 -13.105173,
703 -12.897808,
704 -12.703081,
705 -12.520644,
706 -12.349726,
707 -12.189612,
708 -12.039848,
709 -11.899627,
710 -11.768584,
711 -11.646277,
712 -11.532149,
713 -11.426017,
714 -11.327059,
715 -11.234826,
716 -11.148986,
717 -11.069003,
718 -10.994557,
719 -10.925266,
720 -10.860516,
721 -10.800336,
722 -10.744141,
723 -10.691750,
724 -10.642770,
725 -10.596803,
726 -10.553743
727 };
728
729 double Graph_Rho_y_Im_T13[150] = {0.172797,
730 0.202131,
731 0.235426,
732 0.273128,
733 0.315725,
734 0.363733,
735 0.417710,
736 0.478276,
737 0.546117,
738 0.622003,
739 0.706793,
740 0.801454,
741 0.907079,
742 1.024912,
743 1.156358,
744 1.303057,
745 1.466877,
746 1.649993,
747 1.855002,
748 2.084838,
749 2.343180,
750 2.634180,
751 2.963042,
752 3.335988,
753 3.760336,
754 4.245326,
755 4.801695,
756 5.442363,
757 6.183616,
758 7.043319,
759 8.043708,
760 9.212110,
761 10.577915,
762 12.181168,
763 14.062523,
764 16.279532,
765 18.897777,
766 21.995371,
767 25.675367,
768 30.055224,
769 35.290852,
770 41.573748,
771 49.136717,
772 58.278429,
773 69.369187,
774 82.861222,
775 99.319087,
776 119.371105,
777 143.764923,
778 173.220093,
779 208.276747,
780 249.049908,
781 294.575917,
782 342.478895,
783 388.404311,
784 426.509539,
785 451.047793,
786 458.502875,
787 449.102019,
788 426.338072,
789 395.230875,
790 360.487558,
791 325.592027,
792 292.650111,
793 262.673758,
794 236.014050,
795 212.576840,
796 192.110301,
797 174.272132,
798 158.713074,
799 145.122917,
800 133.200461,
801 122.713454,
802 113.442718,
803 105.210160,
804 97.871589,
805 91.297458,
806 85.382517,
807 80.038571,
808 75.188341,
809 70.772272,
810 66.731762,
811 63.021465,
812 59.603185,
813 56.439346,
814 53.505943,
815 50.773808,
816 48.217359,
817 45.822452,
818 43.568704,
819 41.442980,
820 39.429563,
821 37.515746,
822 35.691573,
823 33.946246,
824 32.276516,
825 30.702591,
826 29.384301,
827 29.135239,
828 29.803776,
829 30.387191,
830 30.800782,
831 31.084135,
832 31.273437,
833 31.394989,
834 31.471681,
835 31.500655,
836 31.483502,
837 31.429924,
838 31.347404,
839 31.241453,
840 31.116455,
841 30.975392,
842 30.820882,
843 30.655076,
844 30.479455,
845 30.295663,
846 30.104566,
847 29.907259,
848 29.704468,
849 29.496734,
850 29.284797,
851 29.068989,
852 28.849764,
853 28.627389,
854 28.402126,
855 28.174274,
856 27.944018,
857 27.711490,
858 27.476881,
859 27.240225,
860 27.001759,
861 26.761511,
862 26.519548,
863 26.276010,
864 26.031002,
865 25.784291,
866 25.535902,
867 25.285895,
868 25.034338,
869 24.781301,
870 24.526823,
871 24.270938,
872 24.013808,
873 23.755328,
874 23.495667,
875 23.234796,
876 22.972822,
877 22.709809,
878 22.445771
879 };
880
881 double Graph_Rho_y_Re_T31[150] = {-18.898530,
882 -18.954238,
883 -19.008932,
884 -19.062622,
885 -19.115336,
886 -19.167093,
887 -19.217916,
888 -19.267833,
889 -19.316876,
890 -19.365079,
891 -19.412470,
892 -19.459096,
893 -19.504998,
894 -19.550217,
895 -19.594806,
896 -19.638821,
897 -19.682309,
898 -19.725339,
899 -19.767980,
900 -19.810299,
901 -19.852371,
902 -19.894290,
903 -19.936132,
904 -19.978009,
905 -20.020017,
906 -20.062275,
907 -20.104905,
908 -20.148048,
909 -20.191848,
910 -20.236468,
911 -20.282075,
912 -20.328867,
913 -20.377060,
914 -20.426869,
915 -20.478591,
916 -20.532417,
917 -20.588689,
918 -20.647744,
919 -20.709872,
920 -20.775615,
921 -20.845352,
922 -20.919569,
923 -20.998812,
924 -21.083699,
925 -21.174885,
926 -21.273119,
927 -21.379240,
928 -21.494175,
929 -21.618967,
930 -21.754770,
931 -21.902869,
932 -22.064711,
933 -22.241895,
934 -22.436216,
935 -22.649648,
936 -22.884333,
937 -23.142670,
938 -23.427171,
939 -23.740492,
940 -24.085340,
941 -24.464253,
942 -24.879401,
943 -25.332139,
944 -25.822282,
945 -26.347212,
946 -26.900283,
947 -27.468808,
948 -28.031350,
949 -28.554390,
950 -28.989253,
951 -29.270353,
952 -29.317394,
953 -29.044877,
954 -28.380249,
955 -27.288896,
956 -25.795851,
957 -23.990377,
958 -22.009025,
959 -19.999765,
960 -18.088432,
961 -16.359488,
962 -14.852245,
963 -13.572492,
964 -12.503751,
965 -11.618787,
966 -10.887591,
967 -10.280414,
968 -9.771085,
969 -9.336666,
970 -8.955685,
971 -8.602002,
972 -8.373614,
973 -8.381858,
974 -8.305447,
975 -8.198825,
976 -8.075954,
977 -7.943039,
978 -7.803117,
979 -7.658536,
980 -7.510620,
981 -7.360422,
982 -7.208926,
983 -7.056351,
984 -6.903115,
985 -6.750076,
986 -6.596933,
987 -6.444168,
988 -6.291804,
989 -6.140002,
990 -5.988871,
991 -5.838466,
992 -5.688884,
993 -5.540139,
994 -5.392265,
995 -5.245312,
996 -5.099299,
997 -4.954260,
998 -4.810212,
999 -4.667167,
1000 -4.525154,
1001 -4.384183,
1002 -4.244278,
1003 -4.105449,
1004 -3.967714,
1005 -3.831086,
1006 -3.695582,
1007 -3.561217,
1008 -3.427992,
1009 -3.295935,
1010 -3.165056,
1011 -3.035363,
1012 -2.906875,
1013 -2.779587,
1014 -2.653525,
1015 -2.528700,
1016 -2.405113,
1017 -2.282784,
1018 -2.161710,
1019 -2.041901,
1020 -1.923378,
1021 -1.806136,
1022 -1.690193,
1023 -1.575549,
1024 -1.462195,
1025 -1.350166,
1026 -1.239448,
1027 -1.130057,
1028 -1.021995,
1029 -0.915245,
1030 -0.809840
1031 };
1032
1033
1034 double Graph_Rho_y_Im_T31[150] = {5.177131,
1035 5.337784,
1036 5.497719,
1037 5.656976,
1038 5.815590,
1039 5.973591,
1040 6.131000,
1041 6.287828,
1042 6.444084,
1043 6.599785,
1044 6.754926,
1045 6.909510,
1046 7.063533,
1047 7.216990,
1048 7.369875,
1049 7.522178,
1050 7.673890,
1051 7.824999,
1052 7.975490,
1053 8.125350,
1054 8.274567,
1055 8.423120,
1056 8.570994,
1057 8.718172,
1058 8.864631,
1059 9.010357,
1060 9.155329,
1061 9.299518,
1062 9.442903,
1063 9.585461,
1064 9.727158,
1065 9.867959,
1066 10.007833,
1067 10.146736,
1068 10.284625,
1069 10.421424,
1070 10.557077,
1071 10.691512,
1072 10.824623,
1073 10.956336,
1074 11.086503,
1075 11.214980,
1076 11.341581,
1077 11.466103,
1078 11.588274,
1079 11.707783,
1080 11.824250,
1081 11.937219,
1082 12.046145,
1083 12.150328,
1084 12.248938,
1085 12.340956,
1086 12.425123,
1087 12.499879,
1088 12.563336,
1089 12.613054,
1090 12.646085,
1091 12.658707,
1092 12.646276,
1093 12.603044,
1094 12.521635,
1095 12.392996,
1096 12.205773,
1097 11.945939,
1098 11.596455,
1099 11.136726,
1100 10.543048,
1101 9.789461,
1102 8.850032,
1103 7.703896,
1104 6.342526,
1105 4.780726,
1106 3.068518,
1107 1.297886,
1108 -0.401210,
1109 -1.882716,
1110 -3.019820,
1111 -3.736728,
1112 -4.024632,
1113 -3.932383,
1114 -3.542528,
1115 -2.946559,
1116 -2.226248,
1117 -1.446060,
1118 -0.651834,
1119 0.126823,
1120 0.872474,
1121 1.577026,
1122 2.239358,
1123 2.864945,
1124 3.474202,
1125 4.270884,
1126 4.721048,
1127 5.079860,
1128 5.392000,
1129 5.671536,
1130 5.925894,
1131 6.159182,
1132 6.374945,
1133 6.575360,
1134 6.762250,
1135 6.937481,
1136 7.101472,
1137 7.255024,
1138 7.400154,
1139 7.536142,
1140 7.664262,
1141 7.784613,
1142 7.897715,
1143 8.004005,
1144 8.103720,
1145 8.197358,
1146 8.285037,
1147 8.366975,
1148 8.443519,
1149 8.514807,
1150 8.581137,
1151 8.642709,
1152 8.699571,
1153 8.752011,
1154 8.800131,
1155 8.844110,
1156 8.884206,
1157 8.920372,
1158 8.952839,
1159 8.981702,
1160 9.007089,
1161 9.029199,
1162 9.048028,
1163 9.063717,
1164 9.076381,
1165 9.086097,
1166 9.093019,
1167 9.097181,
1168 9.098662,
1169 9.097581,
1170 9.093991,
1171 9.088002,
1172 9.079685,
1173 9.069073,
1174 9.056272,
1175 9.041333,
1176 9.024324,
1177 9.005355,
1178 8.984416,
1179 8.961603,
1180 8.936963,
1181 8.910553,
1182 8.882464,
1183 8.852699
1184 };
1185
1186 double Graph_Rho_y_Re_T33[150] = {-44.165564,
1187 -44.782299,
1188 -45.413780,
1189 -46.060859,
1190 -46.724507,
1191 -47.405805,
1192 -48.105886,
1193 -48.826045,
1194 -49.567696,
1195 -50.332251,
1196 -51.121502,
1197 -51.937242,
1198 -52.781461,
1199 -53.656910,
1200 -54.564928,
1201 -55.509720,
1202 -56.493338,
1203 -57.518695,
1204 -58.592089,
1205 -59.712470,
1206 -60.889741,
1207 -62.121776,
1208 -63.416569,
1209 -64.779063,
1210 -66.206418,
1211 -67.708315,
1212 -69.279943,
1213 -70.916366,
1214 -72.631021,
1215 -74.399175,
1216 -76.218163,
1217 -78.099669,
1218 -79.999924,
1219 -81.952793,
1220 -83.886394,
1221 -85.837848,
1222 -87.787357,
1223 -89.698700,
1224 -91.590760,
1225 -93.410871,
1226 -95.168676,
1227 -96.850432,
1228 -98.419853,
1229 -99.869159,
1230 -101.181770,
1231 -102.333399,
1232 -103.319808,
1233 -104.099973,
1234 -104.679432,
1235 -105.037907,
1236 -105.159607,
1237 -105.044851,
1238 -104.676379,
1239 -104.060251,
1240 -103.194744,
1241 -102.084690,
1242 -100.737332,
1243 -99.165891,
1244 -97.380714,
1245 -95.400311,
1246 -93.244871,
1247 -90.927100,
1248 -88.471985,
1249 -85.899878,
1250 -83.220662,
1251 -80.468165,
1252 -77.645458,
1253 -74.776600,
1254 -71.872527,
1255 -68.940500,
1256 -66.002479,
1257 -63.049322,
1258 -60.105192,
1259 -57.162796,
1260 -54.223449,
1261 -51.296106,
1262 -48.375228,
1263 -45.459236,
1264 -42.548048,
1265 -39.633488,
1266 -36.720837,
1267 -33.796985,
1268 -30.858233,
1269 -27.903517,
1270 -24.922299,
1271 -21.926087,
1272 -18.902438,
1273 -15.856397,
1274 -12.793796,
1275 -9.724366,
1276 -6.661738,
1277 -3.613763,
1278 -0.608983,
1279 2.355058,
1280 5.210696,
1281 7.963255,
1282 10.582684,
1283 13.018834,
1284 15.262771,
1285 17.240957,
1286 18.961181,
1287 20.385818,
1288 21.449639,
1289 22.202664,
1290 22.600870,
1291 22.667623,
1292 22.485054,
1293 21.999008,
1294 21.328073,
1295 20.501632,
1296 19.547779,
1297 18.532754,
1298 17.446044,
1299 16.354439,
1300 15.252220,
1301 14.159107,
1302 13.093216,
1303 12.046872,
1304 11.035352,
1305 10.057886,
1306 9.112631,
1307 8.205265,
1308 7.331877,
1309 6.491906,
1310 5.687883,
1311 4.912460,
1312 4.171532,
1313 3.456577,
1314 2.770329,
1315 2.112600,
1316 1.478096,
1317 0.869541,
1318 0.282205,
1319 -0.282582,
1320 -0.826144,
1321 -1.350564,
1322 -1.855644,
1323 -2.342791,
1324 -2.813130,
1325 -3.265676,
1326 -3.704098,
1327 -4.126980,
1328 -4.534898,
1329 -4.929698,
1330 -5.310221,
1331 -5.679083,
1332 -6.034973,
1333 -6.378910,
1334 -6.711708,
1335 -7.032907
1336 };
1337
1338
1339 double Graph_Rho_y_Im_T33[150] = {-0.009574,
1340 -0.011236,
1341 -0.013128,
1342 -0.015280,
1343 -0.017727,
1344 -0.020566,
1345 -0.023942,
1346 -0.028041,
1347 -0.033098,
1348 -0.039414,
1349 -0.047368,
1350 -0.057439,
1351 -0.070236,
1352 -0.086530,
1353 -0.107291,
1354 -0.133759,
1355 -0.167492,
1356 -0.210456,
1357 -0.265179,
1358 -0.334718,
1359 -0.423091,
1360 -0.535009,
1361 -0.676495,
1362 -0.854774,
1363 -1.078020,
1364 -1.356318,
1365 -1.700164,
1366 -2.120676,
1367 -2.631449,
1368 -3.242145,
1369 -3.964043,
1370 -4.811083,
1371 -5.785905,
1372 -6.906905,
1373 -8.165078,
1374 -9.579957,
1375 -11.155927,
1376 -12.887888,
1377 -14.794431,
1378 -16.858743,
1379 -19.095246,
1380 -21.507339,
1381 -24.079181,
1382 -26.813884,
1383 -29.707263,
1384 -32.743940,
1385 -35.930967,
1386 -39.219763,
1387 -42.625779,
1388 -46.121508,
1389 -49.674060,
1390 -53.290255,
1391 -56.901223,
1392 -60.518525,
1393 -64.107220,
1394 -67.635337,
1395 -71.090745,
1396 -74.436419,
1397 -77.669200,
1398 -80.760091,
1399 -83.692600,
1400 -86.467419,
1401 -89.065923,
1402 -91.483225,
1403 -93.726513,
1404 -95.780370,
1405 -97.660409,
1406 -99.363551,
1407 -100.898416,
1408 -102.273414,
1409 -103.491548,
1410 -104.567807,
1411 -105.504611,
1412 -106.314090,
1413 -107.003055,
1414 -107.577203,
1415 -108.045087,
1416 -108.410391,
1417 -108.678139,
1418 -108.849840,
1419 -108.926569,
1420 -108.909725,
1421 -108.794662,
1422 -108.579441,
1423 -108.252061,
1424 -107.808841,
1425 -107.235865,
1426 -106.519939,
1427 -105.651462,
1428 -104.611737,
1429 -103.391656,
1430 -101.976705,
1431 -100.352962,
1432 -98.518143,
1433 -96.457731,
1434 -94.174990,
1435 -91.670379,
1436 -88.946940,
1437 -86.017365,
1438 -82.906955,
1439 -79.632547,
1440 -76.233013,
1441 -72.772859,
1442 -69.273632,
1443 -65.821844,
1444 -62.468224,
1445 -59.215597,
1446 -56.189297,
1447 -53.345014,
1448 -50.709890,
1449 -48.304668,
1450 -46.084162,
1451 -44.094868,
1452 -42.259156,
1453 -40.597342,
1454 -39.083314,
1455 -37.686660,
1456 -36.417188,
1457 -35.242315,
1458 -34.155783,
1459 -33.152847,
1460 -32.214421,
1461 -31.338977,
1462 -30.519467,
1463 -29.742801,
1464 -29.015464,
1465 -28.319097,
1466 -27.663181,
1467 -27.036550,
1468 -26.433956,
1469 -25.860077,
1470 -25.305275,
1471 -24.773761,
1472 -24.259292,
1473 -23.760774,
1474 -23.278570,
1475 -22.809698,
1476 -22.353615,
1477 -21.909736,
1478 -21.474249,
1479 -21.051088,
1480 -20.636091,
1481 -20.228549,
1482 -19.829581,
1483 -19.436313,
1484 -19.051205,
1485 -18.671257,
1486 -18.296820,
1487 -17.928025,
1488 -17.563446
1489 };
1490
1491 double Graph_Om_x[150] = { 1.200010,
1492 1.205343,
1493 1.210677,
1494 1.216010,
1495 1.221343,
1496 1.226677,
1497 1.232010,
1498 1.237343,
1499 1.242677,
1500 1.248010,
1501 1.253343,
1502 1.258677,
1503 1.264010,
1504 1.269343,
1505 1.274677,
1506 1.280010,
1507 1.285343,
1508 1.290677,
1509 1.296010,
1510 1.301343,
1511 1.306677,
1512 1.312010,
1513 1.317343,
1514 1.322677,
1515 1.328010,
1516 1.333343,
1517 1.338677,
1518 1.344010,
1519 1.349343,
1520 1.354677,
1521 1.360010,
1522 1.365343,
1523 1.370677,
1524 1.376010,
1525 1.381343,
1526 1.386677,
1527 1.392010,
1528 1.397343,
1529 1.402677,
1530 1.408010,
1531 1.413343,
1532 1.418677,
1533 1.424010,
1534 1.429343,
1535 1.434677,
1536 1.440010,
1537 1.445343,
1538 1.450677,
1539 1.456010,
1540 1.461343,
1541 1.466677,
1542 1.472010,
1543 1.477343,
1544 1.482677,
1545 1.488010,
1546 1.493343,
1547 1.498677,
1548 1.504010,
1549 1.509343,
1550 1.514677,
1551 1.520010,
1552 1.525343,
1553 1.530677,
1554 1.536010,
1555 1.541343,
1556 1.546677,
1557 1.552010,
1558 1.557343,
1559 1.562677,
1560 1.568010,
1561 1.573343,
1562 1.578677,
1563 1.584010,
1564 1.589343,
1565 1.594677,
1566 1.600010,
1567 1.605343,
1568 1.610677,
1569 1.616010,
1570 1.621343,
1571 1.626677,
1572 1.632010,
1573 1.637343,
1574 1.642677,
1575 1.648010,
1576 1.653343,
1577 1.658677,
1578 1.664010,
1579 1.669343,
1580 1.674677,
1581 1.680010,
1582 1.685343,
1583 1.690677,
1584 1.696010,
1585 1.701343,
1586 1.706677,
1587 1.712010,
1588 1.717343,
1589 1.722677,
1590 1.728010,
1591 1.733343,
1592 1.738677,
1593 1.744010,
1594 1.749343,
1595 1.754677,
1596 1.760010,
1597 1.765343,
1598 1.770677,
1599 1.776010,
1600 1.781343,
1601 1.786677,
1602 1.792010,
1603 1.797343,
1604 1.802677,
1605 1.808010,
1606 1.813343,
1607 1.818677,
1608 1.824010,
1609 1.829343,
1610 1.834677,
1611 1.840010,
1612 1.845343,
1613 1.850677,
1614 1.856010,
1615 1.861343,
1616 1.866677,
1617 1.872010,
1618 1.877343,
1619 1.882677,
1620 1.888010,
1621 1.893343,
1622 1.898677,
1623 1.904010,
1624 1.909343,
1625 1.914677,
1626 1.920010,
1627 1.925343,
1628 1.930677,
1629 1.936010,
1630 1.941343,
1631 1.946677,
1632 1.952010,
1633 1.957343,
1634 1.962677,
1635 1.968010,
1636 1.973343,
1637 1.978677,
1638 1.984010,
1639 1.989343,
1640 1.994677
1641 };
1642
1643
1644
1645
1646
1647
1648
1649 double Graph_Om_y_Im_T11 [150] = { 1.389068,
1650 1.456658,
1651 1.526448,
1652 1.598560,
1653 1.673117,
1654 1.750252,
1655 1.830102,
1656 1.912813,
1657 1.998538,
1658 2.087443,
1659 2.179702,
1660 2.275503,
1661 2.375050,
1662 2.478561,
1663 2.586272,
1664 2.698442,
1665 2.815355,
1666 2.937317,
1667 3.064668,
1668 3.197779,
1669 3.337072,
1670 3.482994,
1671 3.636067,
1672 3.796850,
1673 3.965981,
1674 4.144179,
1675 4.332243,
1676 4.531073,
1677 4.741692,
1678 4.965255,
1679 5.203077,
1680 5.456638,
1681 5.727655,
1682 6.018059,
1683 6.330312,
1684 6.666592,
1685 7.030124,
1686 7.424451,
1687 7.853343,
1688 8.322490,
1689 8.837434,
1690 9.405405,
1691 10.035236,
1692 10.737938,
1693 11.527571,
1694 12.422345,
1695 13.446496,
1696 14.633362,
1697 16.030622,
1698 17.710142,
1699 19.788316,
1700 22.474312,
1701 26.216790,
1702 32.424445,
1703 46.772593,
1704 47.439029,
1705 47.041099,
1706 46.096837,
1707 44.789221,
1708 43.229954,
1709 41.499436,
1710 39.658386,
1711 37.752824,
1712 35.816881,
1713 33.874883,
1714 31.942785,
1715 30.029712,
1716 28.138978,
1717 26.268922,
1718 24.413804,
1719 22.563957,
1720 20.706500,
1721 18.825883,
1722 16.904263,
1723 14.923372,
1724 12.868877,
1725 10.750605,
1726 10.238991,
1727 10.752178,
1728 10.486870,
1729 10.780760,
1730 11.675000,
1731 13.002159,
1732 14.575257,
1733 16.250461,
1734 17.933620,
1735 19.568770,
1736 21.127969,
1737 22.606847,
1738 24.033553,
1739 25.541496,
1740 27.358317,
1741 26.236157,
1742 25.026124,
1743 23.668561,
1744 22.153756,
1745 20.451159,
1746 18.491888,
1747 16.540051,
1748 15.356865,
1749 14.610553,
1750 14.046350,
1751 13.581249,
1752 13.180783,
1753 12.827425,
1754 12.512975,
1755 12.234395,
1756 11.979383,
1757 11.741607,
1758 11.518226,
1759 11.307183,
1760 11.106938,
1761 10.916210,
1762 10.733968,
1763 10.559365,
1764 10.391660,
1765 10.230250,
1766 10.074594,
1767 9.924207,
1768 9.778691,
1769 9.637670,
1770 9.500835,
1771 9.367901,
1772 9.238590,
1773 9.112678,
1774 8.989949,
1775 8.870220,
1776 8.753314,
1777 8.639065,
1778 8.527328,
1779 8.417973,
1780 8.310875,
1781 8.205922,
1782 8.103007,
1783 8.002030,
1784 7.902903,
1785 7.805546,
1786 7.709879,
1787 7.615831,
1788 7.523333,
1789 7.432321,
1790 7.342739,
1791 7.254531,
1792 7.167650,
1793 7.082044,
1794 6.997669,
1795 6.914482,
1796 6.832444,
1797 6.751522,
1798 6.671677
1799 };
1800
1801
1802 double Graph_Om_y_Re_T11 [150] = { 16.676058,
1803 16.873620,
1804 17.074454,
1805 17.278667,
1806 17.486368,
1807 17.697678,
1808 17.912727,
1809 18.131652,
1810 18.354605,
1811 18.581743,
1812 18.813243,
1813 19.049291,
1814 19.290092,
1815 19.535868,
1816 19.786854,
1817 20.043311,
1818 20.305525,
1819 20.573799,
1820 20.848469,
1821 21.129896,
1822 21.418482,
1823 21.714654,
1824 22.018892,
1825 22.331705,
1826 22.653656,
1827 22.985366,
1828 23.327505,
1829 23.680812,
1830 24.046094,
1831 24.424238,
1832 24.816224,
1833 25.223121,
1834 25.646128,
1835 26.086550,
1836 26.546026,
1837 27.025899,
1838 27.528193,
1839 28.055032,
1840 28.608615,
1841 29.192193,
1842 29.808795,
1843 30.462277,
1844 31.157229,
1845 31.899187,
1846 32.695000,
1847 33.553233,
1848 34.484872,
1849 35.504376,
1850 36.631442,
1851 37.894109,
1852 39.334727,
1853 41.023316,
1854 43.095027,
1855 45.921357,
1856 39.521110,
1857 30.351765,
1858 24.284479,
1859 19.471935,
1860 15.470932,
1861 12.091298,
1862 9.223901,
1863 6.793776,
1864 4.742290,
1865 3.020930,
1866 1.588384,
1867 0.409093,
1868 -0.546716,
1869 -1.303340,
1870 -1.879366,
1871 -2.286958,
1872 -2.531089,
1873 -2.607636,
1874 -2.500004,
1875 -2.173386,
1876 -1.561956,
1877 -0.536450,
1878 1.219817,
1879 5.427883,
1880 5.767103,
1881 7.351354,
1882 9.212795,
1883 10.915352,
1884 12.236154,
1885 13.106596,
1886 13.540076,
1887 13.580716,
1888 13.277968,
1889 12.671440,
1890 11.791697,
1891 10.647630,
1892 9.199456,
1893 5.813249,
1894 2.983655,
1895 0.971801,
1896 -0.634399,
1897 -1.927270,
1898 -2.916273,
1899 -3.479861,
1900 -3.186752,
1901 -2.468163,
1902 -1.900464,
1903 -1.464063,
1904 -1.112943,
1905 -0.819212,
1906 -0.565878,
1907 -0.339191,
1908 -0.142748,
1909 0.027742,
1910 0.178540,
1911 0.313819,
1912 0.436449,
1913 0.548524,
1914 0.651670,
1915 0.747157,
1916 0.835989,
1917 0.918994,
1918 0.996849,
1919 1.070123,
1920 1.139306,
1921 1.204797,
1922 1.266949,
1923 1.326061,
1924 1.382391,
1925 1.436179,
1926 1.487623,
1927 1.536891,
1928 1.584151,
1929 1.629527,
1930 1.673151,
1931 1.715135,
1932 1.755573,
1933 1.794553,
1934 1.832154,
1935 1.868451,
1936 1.903506,
1937 1.937379,
1938 1.970120,
1939 2.001780,
1940 2.032400,
1941 2.062027,
1942 2.090689,
1943 2.118424,
1944 2.145266,
1945 2.171235,
1946 2.196363,
1947 2.220677,
1948 2.244195,
1949 2.266939,
1950 2.288926,
1951 2.310179
1952 };
1953
1954 double Graph_Om_y_Re_T13 [150] = { 145.316208,
1955 148.083055,
1956 150.944508,
1957 153.906334,
1958 156.974935,
1959 160.157338,
1960 163.460977,
1961 166.894114,
1962 170.465700,
1963 174.184960,
1964 178.062835,
1965 182.110576,
1966 186.340514,
1967 190.768185,
1968 195.404279,
1969 200.270501,
1970 205.381185,
1971 210.754440,
1972 216.420894,
1973 222.386461,
1974 228.695727,
1975 235.353877,
1976 242.401765,
1977 249.871843,
1978 257.775024,
1979 266.165408,
1980 275.056849,
1981 284.469814,
1982 294.485659,
1983 305.082245,
1984 316.311351,
1985 328.273510,
1986 340.923025,
1987 354.441140,
1988 368.722597,
1989 383.974977,
1990 400.255158,
1991 417.575487,
1992 436.113490,
1993 455.826814,
1994 476.852325,
1995 499.234042,
1996 522.883651,
1997 547.731147,
1998 573.543413,
1999 599.834811,
2000 625.891083,
2001 650.281582,
2002 671.111683,
2003 685.303913,
2004 688.469236,
2005 674.899026,
2006 637.580290,
2007 569.892821,
2008 467.874905,
2009 333.701202,
2010 177.649521,
2011 16.540020,
2012 -132.003962,
2013 -255.160247,
2014 -347.740621,
2015 -411.019480,
2016 -449.700058,
2017 -469.555335,
2018 -475.915437,
2019 -473.013093,
2020 -464.075398,
2021 -451.365215,
2022 -436.463959,
2023 -420.438963,
2024 -403.990574,
2025 -387.584322,
2026 -371.502351,
2027 -355.923812,
2028 -340.946128,
2029 -326.615863,
2030 -312.947150,
2031 -299.932513,
2032 -287.549580,
2033 -275.769611,
2034 -264.558115,
2035 -253.877010,
2036 -243.689490,
2037 -233.956948,
2038 -224.642118,
2039 -215.708356,
2040 -207.122175,
2041 -198.836439,
2042 -190.825146,
2043 -183.045428,
2044 -175.456039,
2045 -168.007818,
2046 -160.641735,
2047 -153.280419,
2048 -145.810654,
2049 -138.048945,
2050 -129.655317,
2051 -119.951993,
2052 -109.456520,
2053 -102.169916,
2054 -97.141669,
2055 -93.133006,
2056 -89.704576,
2057 -86.667909,
2058 -83.924058,
2059 -81.423497,
2060 -79.172024,
2061 -77.099570,
2062 -75.163161,
2063 -73.342901,
2064 -71.624602,
2065 -69.997257,
2066 -68.452046,
2067 -66.981557,
2068 -65.579433,
2069 -64.240275,
2070 -62.959259,
2071 -61.732225,
2072 -60.555390,
2073 -59.425388,
2074 -58.339212,
2075 -57.294075,
2076 -56.287474,
2077 -55.317106,
2078 -54.380864,
2079 -53.476834,
2080 -52.603212,
2081 -51.758295,
2082 -50.940582,
2083 -50.148609,
2084 -49.381060,
2085 -48.636662,
2086 -47.914216,
2087 -47.212642,
2088 -46.530897,
2089 -45.867933,
2090 -45.222887,
2091 -44.594937,
2092 -43.983259,
2093 -43.387081,
2094 -42.805668,
2095 -42.238315,
2096 -41.684341,
2097 -41.143133,
2098 -40.614095,
2099 -40.096635,
2100 -39.590201,
2101 -39.094293,
2102 -38.608399,
2103 -38.132059
2104 };
2105
2106
2107 double Graph_Om_y_Im_T13 [150] = { 0.313878,
2108 0.368610,
2109 0.431055,
2110 0.502141,
2111 0.582931,
2112 0.674837,
2113 0.779583,
2114 0.899176,
2115 1.035959,
2116 1.192666,
2117 1.372513,
2118 1.579273,
2119 1.817399,
2120 2.092184,
2121 2.409804,
2122 2.777770,
2123 3.204856,
2124 3.701555,
2125 4.280654,
2126 4.956692,
2127 5.747987,
2128 6.674972,
2129 7.762763,
2130 9.040322,
2131 10.540098,
2132 12.300707,
2133 14.363690,
2134 16.775194,
2135 19.589250,
2136 22.859508,
2137 26.649545,
2138 31.034828,
2139 36.090930,
2140 41.925515,
2141 48.639537,
2142 56.386859,
2143 65.338516,
2144 75.695877,
2145 87.732655,
2146 101.748067,
2147 118.149329,
2148 137.428874,
2149 160.171127,
2150 187.126538,
2151 219.209043,
2152 257.508320,
2153 303.372286,
2154 358.207936,
2155 423.674585,
2156 501.192578,
2157 591.537257,
2158 694.244016,
2159 805.850223,
2160 919.402082,
2161 1023.20456,
2162 1102.584176,
2163 1144.255389,
2164 1141.482241,
2165 1097.291368,
2166 1022.350076,
2167 930.208991,
2168 832.775961,
2169 738.301212,
2170 651.388388,
2171 573.925390,
2172 506.193188,
2173 447.542576,
2174 397.002511,
2175 353.486102,
2176 315.960002,
2177 283.518367,
2178 255.347314,
2179 230.793268,
2180 209.283499,
2181 190.349776,
2182 173.609728,
2183 158.733608,
2184 145.453946,
2185 133.545571,
2186 122.818811,
2187 113.120402,
2188 104.309579,
2189 96.275395,
2190 88.921457,
2191 82.162287,
2192 75.931660,
2193 70.165414,
2194 64.801224,
2195 59.802191,
2196 55.121033,
2197 50.722625,
2198 46.570594,
2199 42.633929,
2200 38.887847,
2201 35.313340,
2202 31.917492,
2203 28.800893,
2204 26.552234,
2205 27.660677,
2206 31.165483,
2207 34.081518,
2208 36.317887,
2209 38.063159,
2210 39.455256,
2211 40.590090,
2212 41.550560,
2213 42.324927,
2214 42.924916,
2215 43.384620,
2216 43.730864,
2217 43.982760,
2218 44.155337,
2219 44.259936,
2220 44.305865,
2221 44.300907,
2222 44.251282,
2223 44.162455,
2224 44.038929,
2225 43.884433,
2226 43.702406,
2227 43.495663,
2228 43.266757,
2229 43.018004,
2230 42.751156,
2231 42.468056,
2232 42.170154,
2233 41.858889,
2234 41.535493,
2235 41.200979,
2236 40.856391,
2237 40.502626,
2238 40.140464,
2239 39.770688,
2240 39.393967,
2241 39.010874,
2242 38.622006,
2243 38.227975,
2244 37.829254,
2245 37.426253,
2246 37.019384,
2247 36.609005,
2248 36.195493,
2249 35.779189,
2250 35.360378,
2251 34.939372,
2252 34.516464,
2253 34.091907,
2254 33.665937,
2255 33.238826,
2256 32.810787
2257 };
2258
2259 TGraph *Graph_Rho_Re_T11 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Re_T11 );
2260 TGraph *Graph_Rho_Im_T11 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Im_T11 );
2261 TGraph *Graph_Rho_Re_T13 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Re_T13 );
2262 TGraph *Graph_Rho_Im_T13 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Im_T13 );
2263 TGraph *Graph_Rho_Re_T31 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Re_T31 );
2264 TGraph *Graph_Rho_Im_T31 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Im_T31 );
2265 TGraph *Graph_Rho_Re_T33 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Re_T33 );
2266 TGraph *Graph_Rho_Im_T33 = new TGraph(150, Graph_Rho_x, Graph_Rho_y_Im_T33 );
2267
2268 TGraph *Graph_Om_Re_T11 = new TGraph(150, Graph_Om_x, Graph_Om_y_Re_T11 );
2269 TGraph *Graph_Om_Im_T11 = new TGraph(150, Graph_Om_x, Graph_Om_y_Im_T11 );
2270 TGraph *Graph_Om_Re_T13 = new TGraph(150, Graph_Om_x, Graph_Om_y_Re_T13 );
2271 TGraph *Graph_Om_Im_T13 = new TGraph(150, Graph_Om_x, Graph_Om_y_Im_T13 );
2272
2273 Pi_minusBeamAmplitude->SetGraph_Rho_Re_T11(Graph_Rho_Re_T11);
2274 Pi_minusBeamAmplitude->SetGraph_Rho_Im_T11(Graph_Rho_Im_T11);
2275 Pi_minusBeamAmplitude->SetGraph_Rho_Re_T13(Graph_Rho_Re_T13);
2276 Pi_minusBeamAmplitude->SetGraph_Rho_Im_T13(Graph_Rho_Im_T13);
2277 Pi_minusBeamAmplitude->SetGraph_Rho_Re_T31(Graph_Rho_Re_T31);
2278 Pi_minusBeamAmplitude->SetGraph_Rho_Im_T31(Graph_Rho_Im_T31);
2279 Pi_minusBeamAmplitude->SetGraph_Rho_Re_T33(Graph_Rho_Re_T33);
2280 Pi_minusBeamAmplitude->SetGraph_Rho_Im_T33(Graph_Rho_Im_T33);
2281 Pi_minusBeamAmplitude->SetGraph_Om_Re_T11(Graph_Om_Re_T11);
2282 Pi_minusBeamAmplitude->SetGraph_Om_Im_T11(Graph_Om_Im_T11);
2283 Pi_minusBeamAmplitude->SetGraph_Om_Re_T13(Graph_Om_Re_T13);
2284 Pi_minusBeamAmplitude->SetGraph_Om_Im_T13(Graph_Om_Im_T13);
2285
2286 Pi_minusBeamAmplitude->EnableWeighting();
2287 Pi_minusBeamAmplitude->SetExpectedWeightMean(-1);
2288 Pi_minusBeamAmplitude->SetVersionFlag(512);
2289
2290
2291 pdmutil->Add(Pi_minusBeamAmplitude);
2292
2293
2294
2295 Pi_plusBeamAmplitude->SetGraph_Rho_Re_T11(Graph_Rho_Re_T11);
2296 Pi_plusBeamAmplitude->SetGraph_Rho_Im_T11(Graph_Rho_Im_T11);
2297 Pi_plusBeamAmplitude->SetGraph_Rho_Re_T13(Graph_Rho_Re_T13);
2298 Pi_plusBeamAmplitude->SetGraph_Rho_Im_T13(Graph_Rho_Im_T13);
2299 Pi_plusBeamAmplitude->SetGraph_Rho_Re_T31(Graph_Rho_Re_T31);
2300 Pi_plusBeamAmplitude->SetGraph_Rho_Im_T31(Graph_Rho_Im_T31);
2301 Pi_plusBeamAmplitude->SetGraph_Rho_Re_T33(Graph_Rho_Re_T33);
2302 Pi_plusBeamAmplitude->SetGraph_Rho_Im_T33(Graph_Rho_Im_T33);
2303 Pi_plusBeamAmplitude->SetGraph_Om_Re_T11(Graph_Om_Re_T11);
2304 Pi_plusBeamAmplitude->SetGraph_Om_Im_T11(Graph_Om_Im_T11);
2305 Pi_plusBeamAmplitude->SetGraph_Om_Re_T13(Graph_Om_Re_T13);
2306 Pi_plusBeamAmplitude->SetGraph_Om_Im_T13(Graph_Om_Im_T13);
2307
2308 Pi_plusBeamAmplitude->EnableWeighting();
2309 Pi_plusBeamAmplitude->SetExpectedWeightMean(-1);
2310 Pi_plusBeamAmplitude->SetVersionFlag(512);
2311
2312 pdmutil->Add(Pi_plusBeamAmplitude);
2313
2314 if (Pi_plusBeamAmplitude && Pi_minusBeamAmplitude) return kTRUE;
2315 return kFALSE;
2316}
2317
2320
2321
2322Bool_t PPionBeamPlugin::ExecCommand(const char *, Double_t) {
2323 return kTRUE;
2324}
2325
2326
2327
2329
2330
2331
ClassImp(PBeamLineSimulation) PBeamLineSimulation
PDistributionManagerUtil * makeDistributionManagerUtil()
PStaticData * makeStaticData()
int Add(PDistribution *dist)
Int_t AddSubGroup(const Char_t *id, const Char_t *de, const Char_t *sub)
void SetGroup(const Char_t *id)
Bool_t ExecCommand(const char *command, Double_t value)
Bool_t Activate(void)
PPionBeamPlugin(const Char_t *id, const Char_t *de)
int AddDecay(int didx, const char *name, const char *parent, const char *daughters, double br)
int GetParticleID(const char *id, int warn=1)
int AddParticle(int pid, const char *name, double mass)
TF1 * flat