59 "pi- + p",
"n,dilepton", 1.0 );
68 "pi+ + n",
"p,dilepton", 1.0 );
71 pdmutil->
AddSubGroup(
"pion_beam",
"Elementary pion beam reactions",
"root");
75 Pi_minusBeamAmplitude =
77 "Model from Lutz/Soyeur/Friman", -1);
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");
87 "Model from Lutz/Soyeur/Friman", -1);
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");
95 TF1 *
flat =
new TF1(
"pion_flat",
"1.", 0, 0.99);
98 Pi_minusBeamAmplitude_gen =
99 new PInclusiveModel(
"pi-_flat@pi- + p_brems_n_dilepton/generator",
"Dilepton generator", -1);
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);
107 Pi_minusBeamAmplitude_gen->EnableGenerator();
108 pdmutil->
Add(Pi_minusBeamAmplitude_gen);
110 Pi_plusBeamAmplitude_gen =
111 new PInclusiveModel(
"pi+_flat@pi+ + n_brems_p_dilepton/generator",
"Dilepton generator", -1);
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);
119 Pi_plusBeamAmplitude_gen->EnableGenerator();
120 pdmutil->
Add(Pi_plusBeamAmplitude_gen);
123 double Graph_Rho_x[150] = { 1.200010,
274 double Graph_Rho_y_Re_T11[150] = {10.927144,
425 double Graph_Rho_y_Im_T11[150] = {0.910200,
577 double Graph_Rho_y_Re_T13[150] = {79.999754,
729 double Graph_Rho_y_Im_T13[150] = {0.172797,
881 double Graph_Rho_y_Re_T31[150] = {-18.898530,
1034 double Graph_Rho_y_Im_T31[150] = {5.177131,
1186 double Graph_Rho_y_Re_T33[150] = {-44.165564,
1339 double Graph_Rho_y_Im_T33[150] = {-0.009574,
1491 double Graph_Om_x[150] = { 1.200010,
1649 double Graph_Om_y_Im_T11 [150] = { 1.389068,
1802 double Graph_Om_y_Re_T11 [150] = { 16.676058,
1954 double Graph_Om_y_Re_T13 [150] = { 145.316208,
2107 double Graph_Om_y_Im_T13 [150] = { 0.313878,
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 );
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 );
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);
2286 Pi_minusBeamAmplitude->EnableWeighting();
2287 Pi_minusBeamAmplitude->SetExpectedWeightMean(-1);
2288 Pi_minusBeamAmplitude->SetVersionFlag(512);
2291 pdmutil->
Add(Pi_minusBeamAmplitude);
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);
2308 Pi_plusBeamAmplitude->EnableWeighting();
2309 Pi_plusBeamAmplitude->SetExpectedWeightMean(-1);
2310 Pi_plusBeamAmplitude->SetVersionFlag(512);
2312 pdmutil->
Add(Pi_plusBeamAmplitude);
2314 if (Pi_plusBeamAmplitude && Pi_minusBeamAmplitude)
return kTRUE;