55 using std::ostringstream;
60 using namespace genie;
76 LOG(
"testINukeHadroData",
pINFO) <<
"Reading Event Energy";
79 LOG(
"testINukeHadroData",
pINFO) <<
"Unspecified energy, write to file";
98 string data_dir = (gSystem->Getenv(
"GINUKEHADRONDATA")) ?
99 string(gSystem->Getenv(
"GINUKEHADRONDATA")) :
100 string(gSystem->Getenv(
"GENIE")) +
string(
"/data/intranuke/");
103 <<
"Saving INTRANUKE hadron x-section data to file: " <<
filename;
105 TFile
f(filename.c_str(),
"UPDATE");
109 double cstep = 2.0 / (numPoints);
113 const int hN_kpNelas_nfiles = 18;
114 const int hN_kpNelas_points_per_file = 37;
115 const int hN_kpNelas_npoints = hN_kpNelas_points_per_file * hN_kpNelas_nfiles;
117 int hN_kpNelas_energies[hN_kpNelas_nfiles] = {
118 100, 200, 300, 400, 500, 600, 700, 800, 900, 1000,
119 1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800
122 TGraph2D * hN_kpNelas_graph =
new TGraph2D(hN_kpNelas_npoints);
123 TGraph2D * hN_kpNelas_int =
new TGraph2D(hN_kpNelas_nfiles*numPoints);
128 for(
int ifile = 0; ifile < hN_kpNelas_nfiles; ifile++) {
130 ostringstream hN_datafile;
131 int ke = hN_kpNelas_energies[ifile];
132 hN_datafile << data_dir <<
"/diff_ang/kpn/kpn" << ke <<
".txt";
134 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
135 for(
int i=0; i<buff->GetN(); i++) {
136 buff -> GetPoint(i,x,y);
137 hN_kpNelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
153 for (
int ifile=0;ifile<hN_kpNelas_nfiles;ifile++) {
154 int subdiv = numPoints*ifile;
155 int ke = hN_kpNelas_energies[ifile];
156 for(
int j=0;j<numPoints; j++)
158 hN_kpNelas_int -> SetPoint(subdiv+j,
double(ke),-1+j*cstep,
163 hN_kpNelas_graph ->
Write(
"d_kpN_elas_data");
164 if (test_intbounce) {
165 hN_kpNelas_int ->
Write(
"d_kpN_elas_int");
167 delete hN_kpNelas_graph;
168 delete hN_kpNelas_int;
173 const int hN_kpPelas_nfiles = 18;
174 const int hN_kpPelas_points_per_file = 37;
175 const int hN_kpPelas_npoints = hN_kpPelas_points_per_file * hN_kpPelas_nfiles;
177 int hN_kpPelas_energies[hN_kpPelas_nfiles] = {
178 100, 200, 300, 400, 500, 600, 700, 800, 900, 1000,
179 1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800
182 TGraph2D * hN_kpPelas_graph =
new TGraph2D(hN_kpPelas_npoints);
184 TGraph2D * hN_kpPelas_int =
new TGraph2D(hN_kpPelas_nfiles*numPoints);
189 for(
int ifile = 0; ifile < hN_kpPelas_nfiles; ifile++) {
191 ostringstream hN_datafile;
192 int ke = hN_kpPelas_energies[ifile];
193 hN_datafile << data_dir <<
"/diff_ang/kpp/kpp" << ke <<
".txt";
195 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
196 for(
int i=0; i<buff->GetN(); i++) {
197 buff -> GetPoint(i,x,y);
198 hN_kpPelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
205 int subdiv = numPoints*ifile;
207 for(
int j=0;j<numPoints; j++)
212 hN_kpPelas_int -> SetPoint(subdiv+j,
double(ke),-1.+
double(j)*cstep,
228 hN_kpPelas_graph ->
Write(
"d_kpP_elas_data");
231 hN_kpPelas_int ->
Write(
"d_kpP_elas_int");
234 delete hN_kpPelas_graph;
235 delete hN_kpPelas_int;
241 const int hN_ppelas_nfiles = 20;
242 const int hN_ppelas_points_per_file = 21;
243 const int hN_ppelas_npoints = hN_ppelas_points_per_file * hN_ppelas_nfiles;
245 int hN_ppelas_energies[hN_ppelas_nfiles] = {
246 50, 100, 150, 200, 250, 300, 350, 400, 450, 500,
247 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000
250 TGraph2D * hN_ppelas_graph =
new TGraph2D(hN_ppelas_npoints);
251 TGraph2D * hN_ppelas_int =
new TGraph2D(hN_ppelas_nfiles*numPoints);
257 for(
int ifile = 0; ifile < hN_ppelas_nfiles; ifile++) {
259 ostringstream hN_datafile;
260 int ke = hN_ppelas_energies[ifile];
261 hN_datafile << data_dir <<
"/diff_ang/pp/pp" << ke <<
".txt";
263 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
264 for(
int i=0; i<buff->GetN(); i++) {
265 buff -> GetPoint(i,x,y);
266 hN_ppelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
273 int subdiv = numPoints*ifile;
274 for(
int j=0;j<numPoints; j++)
276 hN_ppelas_int -> SetPoint(subdiv+j,
double(ke),-1+j*cstep,
282 hN_ppelas_graph ->
Write(
"d_pp_elas_data");
285 hN_ppelas_int ->
Write(
"d_pp_elas_int");
288 delete hN_ppelas_graph;
289 delete hN_ppelas_int;
295 const int hN_npelas_nfiles = 20;
296 const int hN_npelas_points_per_file = 21;
297 const int hN_npelas_npoints = hN_npelas_points_per_file * hN_npelas_nfiles;
299 int hN_npelas_energies[hN_npelas_nfiles] = {
300 50, 100, 150, 200, 250, 300, 350, 400, 450, 500,
301 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000
304 TGraph2D * hN_npelas_graph =
new TGraph2D(hN_npelas_npoints);
305 TGraph2D * hN_npelas_int =
new TGraph2D(hN_npelas_nfiles*numPoints);
311 for(
int ifile = 0; ifile < hN_npelas_nfiles; ifile++) {
313 ostringstream hN_datafile;
314 int ke = hN_npelas_energies[ifile];
315 hN_datafile << data_dir <<
"/diff_ang/pn/pn" << ke <<
".txt";
317 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
318 for(
int i=0; i<buff->GetN(); i++) {
319 buff -> GetPoint(i,x,y);
320 hN_npelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
327 int subdiv = numPoints*ifile;
328 for(
int j=0;j<numPoints; j++)
330 hN_npelas_int -> SetPoint(subdiv+j,
double(ke),-1+j*cstep,
336 hN_npelas_graph ->
Write(
"d_np_elas_data");
339 hN_npelas_int ->
Write(
"d_np_elas_int");
342 delete hN_npelas_graph;
343 delete hN_npelas_int;
349 const int hN_pipNelas_nfiles = 60;
350 const int hN_pipNelas_points_per_file = 21;
351 const int hN_pipNelas_npoints = hN_pipNelas_points_per_file * hN_pipNelas_nfiles;
353 int hN_pipNelas_energies[hN_pipNelas_nfiles] = {
354 10, 20, 30, 40, 50, 60, 70, 80, 90,
355 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
356 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
357 300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
358 700, 740, 780, 820, 860, 900, 940, 980,
359 1020, 1060, 1100, 1140, 1180, 1220, 1260,
360 1300, 1340, 1380, 1420, 1460, 1500
363 TGraph2D * hN_pipNelas_graph =
new TGraph2D(hN_pipNelas_npoints);
364 TGraph2D * hN_pipNelas_int =
new TGraph2D(hN_pipNelas_nfiles*numPoints);
370 for(
int ifile = 0; ifile < hN_pipNelas_nfiles; ifile++) {
372 ostringstream hN_datafile;
373 int ke = hN_pipNelas_energies[ifile];
374 hN_datafile << data_dir <<
"/diff_ang/pip/pip" << ke <<
".txt";
376 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
377 for(
int i=0; i<buff->GetN(); i++) {
378 buff -> GetPoint(i,x,y);
379 hN_pipNelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
386 int subdiv = numPoints*ifile;
387 for(
int j=0;j<numPoints; j++)
389 hN_pipNelas_int -> SetPoint(subdiv+j,
double(ke),-1+j*cstep,
397 hN_pipNelas_graph ->
Write(
"d_pipN_elas_data");
400 hN_pipNelas_int ->
Write(
"d_pipN_elas_int");
403 delete hN_pipNelas_graph;
404 delete hN_pipNelas_int;
410 const int hN_pi0Nelas_nfiles = 60;
411 const int hN_pi0Nelas_points_per_file = 21;
412 const int hN_pi0Nelas_npoints = hN_pi0Nelas_points_per_file * hN_pi0Nelas_nfiles;
414 int hN_pi0Nelas_energies[hN_pi0Nelas_nfiles] = {
415 10, 20, 30, 40, 50, 60, 70, 80, 90,
416 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
417 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
418 300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
419 700, 740, 780, 820, 860, 900, 940, 980,
420 1020, 1060, 1100, 1140, 1180, 1220, 1260,
421 1300, 1340, 1380, 1420, 1460, 1500
424 TGraph2D * hN_pi0Nelas_graph =
new TGraph2D(hN_pi0Nelas_npoints);
425 TGraph2D * hN_pi0Nelas_int =
new TGraph2D(hN_pi0Nelas_nfiles*numPoints);
431 for(
int ifile = 0; ifile < hN_pi0Nelas_nfiles; ifile++) {
433 ostringstream hN_datafile;
434 int ke = hN_pi0Nelas_energies[ifile];
435 hN_datafile << data_dir <<
"/diff_ang/pip/pip" << ke <<
".txt";
437 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
438 for(
int i=0; i<buff->GetN(); i++) {
439 buff -> GetPoint(i,x,y);
440 hN_pi0Nelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
447 int subdiv = numPoints*ifile;
448 for(
int j=0;j<numPoints; j++)
450 hN_pi0Nelas_int -> SetPoint(subdiv+j,
double(ke),-1+j*cstep,
458 hN_pi0Nelas_graph ->
Write(
"d_pi0N_elas_data");
461 hN_pi0Nelas_int ->
Write(
"d_pi0N_elas_int");
464 delete hN_pi0Nelas_graph;
465 delete hN_pi0Nelas_int;
471 const int hN_pimNelas_nfiles = 60;
472 const int hN_pimNelas_points_per_file = 21;
473 const int hN_pimNelas_npoints = hN_pimNelas_points_per_file * hN_pimNelas_nfiles;
475 int hN_pimNelas_energies[hN_pimNelas_nfiles] = {
476 10, 20, 30, 40, 50, 60, 70, 80, 90,
477 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
478 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
479 300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
480 700, 740, 780, 820, 860, 900, 940, 980,
481 1020, 1060, 1100, 1140, 1180, 1220, 1260,
482 1300, 1340, 1380, 1420, 1460, 1500
485 TGraph2D * hN_pimNelas_graph =
new TGraph2D(hN_pimNelas_npoints);
490 for(
int ifile = 0; ifile < hN_pimNelas_nfiles; ifile++) {
492 ostringstream hN_datafile;
493 int ke = hN_pimNelas_energies[ifile];
494 hN_datafile << data_dir <<
"/diff_ang/pim/pim" << ke <<
".txt";
496 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
497 for(
int i=0; i<buff->GetN(); i++) {
498 buff -> GetPoint(i,x,y);
499 hN_pimNelas_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
505 hN_pimNelas_graph ->
Write(
"d_pimN_elas_data");
506 delete hN_pimNelas_graph;
511 const int hN_piNcex_nfiles = 60;
512 const int hN_piNcex_points_per_file = 21;
513 const int hN_piNcex_npoints = hN_piNcex_points_per_file * hN_piNcex_nfiles;
515 int hN_piNcex_energies[hN_piNcex_nfiles] = {
516 10, 20, 30, 40, 50, 60, 70, 80, 90,
517 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
518 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
519 300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
520 700, 740, 780, 820, 860, 900, 940, 980,
521 1020, 1060, 1100, 1140, 1180, 1220, 1260,
522 1300, 1340, 1380, 1420, 1460, 1500
525 TGraph2D * hN_piNcex_graph =
new TGraph2D(hN_piNcex_npoints);
530 for(
int ifile = 0; ifile < hN_piNcex_nfiles; ifile++) {
532 ostringstream hN_datafile;
533 int ke = hN_piNcex_energies[ifile];
534 hN_datafile << data_dir <<
"/diff_ang/pie/pie" << ke <<
".txt";
536 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
537 for(
int i=0; i<buff->GetN(); i++) {
538 buff -> GetPoint(i,x,y);
539 hN_piNcex_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
545 hN_piNcex_graph ->
Write(
"d_piN_cex_data");
546 delete hN_piNcex_graph;
551 const int hN_piNabs_nfiles = 19;
552 const int hN_piNabs_points_per_file = 21;
553 const int hN_piNabs_npoints = hN_piNabs_points_per_file * hN_piNabs_nfiles;
555 int hN_piNabs_energies[hN_piNabs_nfiles] = {
556 50, 75, 100, 125, 150, 175, 200, 225, 250, 275,
557 300, 325, 350, 375, 400, 425, 450, 475, 500
560 TGraph2D * hN_piNabs_graph =
new TGraph2D(hN_piNabs_npoints);
565 for(
int ifile = 0; ifile < hN_piNabs_nfiles; ifile++) {
567 ostringstream hN_datafile;
568 int ke = hN_piNabs_energies[ifile];
569 hN_datafile << data_dir <<
"/diff_ang/pid2p/pid2p" << ke <<
".txt";
571 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
572 for(
int i=0; i<buff->GetN(); i++) {
573 buff -> GetPoint(i,x,y);
574 hN_piNabs_graph -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
580 hN_piNabs_graph ->
Write(
"d_piN_abs_data");
581 delete hN_piNabs_graph;
586 const int hN_gampi0pInelas_nfiles = 29;
587 const int hN_gampi0pInelas_points_per_file = 37;
588 const int hN_gampi0pInelas_npoints = hN_gampi0pInelas_points_per_file * hN_gampi0pInelas_nfiles;
590 int hN_gampi0pInelas_energies[hN_gampi0pInelas_nfiles] = {
591 160, 180, 200, 220, 240, 260, 280, 300, 320, 340,
592 360, 380, 400, 450, 500, 550, 600, 650, 700, 750,
593 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200
596 TGraph2D * fhN2dXSecGamPi0P_Inelas =
new TGraph2D(hN_gampi0pInelas_npoints);
601 for(
int ifile = 0; ifile < hN_gampi0pInelas_nfiles; ifile++) {
603 ostringstream hN_datafile;
604 int ke = hN_gampi0pInelas_energies[ifile];
605 hN_datafile << data_dir <<
"/diff_ang/gampi0p/" << ke <<
"-pi0p.txt";
607 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
608 for(
int i=0; i<buff->GetN(); i++) {
609 buff -> GetPoint(i,x,y);
610 fhN2dXSecGamPi0P_Inelas -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
616 fhN2dXSecGamPi0P_Inelas ->
Write(
"d_gam_pi0p_inel_data");
617 delete fhN2dXSecGamPi0P_Inelas;
622 const int hN_gampi0nInelas_nfiles = 29;
623 const int hN_gampi0nInelas_points_per_file = 37;
624 const int hN_gampi0nInelas_npoints = hN_gampi0nInelas_points_per_file * hN_gampi0nInelas_nfiles;
626 int hN_gampi0nInelas_energies[hN_gampi0nInelas_nfiles] = {
627 160, 180, 200, 220, 240, 260, 280, 300, 320, 340,
628 360, 380, 400, 450, 500, 550, 600, 650, 700, 750,
629 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200
632 TGraph2D * fhN2dXSecGamPi0N_Inelas =
new TGraph2D(hN_gampi0nInelas_npoints);
637 for(
int ifile = 0; ifile < hN_gampi0nInelas_nfiles; ifile++) {
639 ostringstream hN_datafile;
640 int ke = hN_gampi0nInelas_energies[ifile];
641 hN_datafile << data_dir <<
"/diff_ang/gampi0n/" << ke <<
"-pi0n.txt";
643 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
644 for(
int i=0; i<buff->GetN(); i++) {
645 buff -> GetPoint(i,x,y);
646 fhN2dXSecGamPi0N_Inelas -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
652 fhN2dXSecGamPi0N_Inelas ->
Write(
"d_gam_pi0n_inel_data");
653 delete fhN2dXSecGamPi0N_Inelas;
658 const int hN_gampipnInelas_nfiles = 29;
659 const int hN_gampipnInelas_points_per_file = 37;
660 const int hN_gampipnInelas_npoints = hN_gampipnInelas_points_per_file * hN_gampipnInelas_nfiles;
662 int hN_gampipnInelas_energies[hN_gampipnInelas_nfiles] = {
663 160, 180, 200, 220, 240, 260, 280, 300, 320, 340,
664 360, 380, 400, 450, 500, 550, 600, 650, 700, 750,
665 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200
668 TGraph2D * fhN2dXSecGamPipN_Inelas =
new TGraph2D(hN_gampipnInelas_npoints);
673 for(
int ifile = 0; ifile < hN_gampipnInelas_nfiles; ifile++) {
675 ostringstream hN_datafile;
676 int ke = hN_gampipnInelas_energies[ifile];
677 hN_datafile << data_dir <<
"/diff_ang/gampi+n/" << ke <<
"-pi+n.txt";
679 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
680 for(
int i=0; i<buff->GetN(); i++) {
681 buff -> GetPoint(i,x,y);
682 fhN2dXSecGamPipN_Inelas -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
688 fhN2dXSecGamPipN_Inelas ->
Write(
"d_gam_pipn_inel_data");
689 delete fhN2dXSecGamPipN_Inelas;
694 const int hN_gampimpInelas_nfiles = 29;
695 const int hN_gampimpInelas_points_per_file = 37;
696 const int hN_gampimpInelas_npoints = hN_gampimpInelas_points_per_file * hN_gampimpInelas_nfiles;
698 int hN_gampimpInelas_energies[hN_gampimpInelas_nfiles] = {
699 160, 180, 200, 220, 240, 260, 280, 300, 320, 340,
700 360, 380, 400, 450, 500, 550, 600, 650, 700, 750,
701 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200
704 TGraph2D * fhN2dXSecGamPimP_Inelas =
new TGraph2D(hN_gampimpInelas_npoints);
709 for(
int ifile = 0; ifile < hN_gampimpInelas_nfiles; ifile++) {
711 ostringstream hN_datafile;
712 int ke = hN_gampimpInelas_energies[ifile];
713 hN_datafile << data_dir <<
"/diff_ang/gampi-p/" << ke <<
"-pi-p.txt";
715 TGraph * buff =
new TGraph(hN_datafile.str().c_str());
716 for(
int i=0; i<buff->GetN(); i++) {
717 buff -> GetPoint(i,x,y);
718 fhN2dXSecGamPimP_Inelas -> SetPoint(ipoint,(
double)ke,TMath::Cos(x*
kPi/180.),y);
724 fhN2dXSecGamPimP_Inelas ->
Write(
"d_gam_pimp_inel_data");
725 delete fhN2dXSecGamPimP_Inelas;
736 <<
"Saving INTRANUKE hadron x-section splines to file: " <<
filename;
753 inukd -> XSecPn_Tot () -> SaveAsROOT(filename,
"pn_tot",
true);
754 inukd -> XSecPn_Elas () -> SaveAsROOT(filename,
"pn_elas",
false);
755 inukd -> XSecPn_Reac () -> SaveAsROOT(filename,
"pn_reac",
false);
756 inukd -> XSecNn_Tot () -> SaveAsROOT(filename,
"nn_tot",
false);
757 inukd -> XSecNn_Elas () -> SaveAsROOT(filename,
"nn_elas",
false);
758 inukd -> XSecNn_Reac () -> SaveAsROOT(filename,
"nn_reac",
false);
759 inukd -> XSecPp_Tot () -> SaveAsROOT(filename,
"pp_tot",
true);
760 inukd -> XSecPp_Elas () -> SaveAsROOT(filename,
"pp_elas",
false);
761 inukd -> XSecPp_Reac () -> SaveAsROOT(filename,
"pp_reac",
false);
763 inukd -> XSecPipn_Tot () -> SaveAsROOT(filename,
"pipn_tot",
false);
764 inukd -> XSecPipn_CEx () -> SaveAsROOT(filename,
"pipn_cex",
false);
765 inukd -> XSecPipn_Elas () -> SaveAsROOT(filename,
"pipn_elas",
false);
766 inukd -> XSecPipn_Reac () -> SaveAsROOT(filename,
"pipn_reac",
false);
767 inukd -> XSecPi0n_Tot () -> SaveAsROOT(filename,
"pi0n_tot",
false);
768 inukd -> XSecPi0n_CEx () -> SaveAsROOT(filename,
"pi0n_cex",
false);
769 inukd -> XSecPi0n_Elas () -> SaveAsROOT(filename,
"pi0n_elas",
false);
770 inukd -> XSecPi0n_Reac () -> SaveAsROOT(filename,
"pi0n_reac",
false);
771 inukd -> XSecPipp_Tot () -> SaveAsROOT(filename,
"pipp_tot",
false);
772 inukd -> XSecPipp_CEx () -> SaveAsROOT(filename,
"pipp_cex",
false);
773 inukd -> XSecPipp_Elas () -> SaveAsROOT(filename,
"pipp_elas",
false);
774 inukd -> XSecPipp_Reac () -> SaveAsROOT(filename,
"pipp_reac",
false);
775 inukd -> XSecPipp_Abs () -> SaveAsROOT(filename,
"pipp_abs",
false);
776 inukd -> XSecPi0p_Tot () -> SaveAsROOT(filename,
"pi0p_tot",
false);
777 inukd -> XSecPi0p_CEx () -> SaveAsROOT(filename,
"pi0p_cex",
false);
778 inukd -> XSecPi0p_Elas () -> SaveAsROOT(filename,
"pi0p_elas",
false);
779 inukd -> XSecPi0p_Reac () -> SaveAsROOT(filename,
"pi0p_reac",
false);
780 inukd -> XSecPi0p_Abs () -> SaveAsROOT(filename,
"pi0p_abs",
false);
782 inukd -> XSecKpn_Elas () -> SaveAsROOT(filename,
"kpn_elas",
false);
783 inukd -> XSecKpp_Elas () -> SaveAsROOT(filename,
"kpp_elas",
false);
784 inukd -> XSecGamp_fs () -> SaveAsROOT(filename,
"gamp_fs",
false);
785 inukd -> XSecGamn_fs () -> SaveAsROOT(filename,
"gamn_fs",
false);
786 inukd -> XSecGamN_Tot () -> SaveAsROOT(filename,
"gamN_tot",
false);
788 inukd -> FracPA_Tot () -> SaveAsROOT(filename,
"pA_tot",
false);
789 inukd -> FracPA_Elas () -> SaveAsROOT(filename,
"pA_elas",
false);
790 inukd -> FracPA_Inel () -> SaveAsROOT(filename,
"pA_inel",
false);
791 inukd -> FracPA_CEx () -> SaveAsROOT(filename,
"pA_cex",
false);
792 inukd -> FracPA_Abs () -> SaveAsROOT(filename,
"pA_abs",
false);
793 inukd -> FracPA_Pipro () -> SaveAsROOT(filename,
"pA_pipro",
false);
794 inukd -> FracNA_Tot () -> SaveAsROOT(filename,
"nA_tot",
false);
795 inukd -> FracNA_Elas () -> SaveAsROOT(filename,
"nA_elas",
false);
796 inukd -> FracNA_Inel () -> SaveAsROOT(filename,
"nA_inel",
false);
797 inukd -> FracNA_CEx () -> SaveAsROOT(filename,
"nA_cex",
false);
798 inukd -> FracNA_Abs () -> SaveAsROOT(filename,
"nA_abs",
false);
799 inukd -> FracNA_Pipro () -> SaveAsROOT(filename,
"nA_pipro",
false);
800 inukd -> FracPipA_Tot () -> SaveAsROOT(filename,
"pipA_tot",
false);
801 inukd -> FracPipA_Elas () -> SaveAsROOT(filename,
"pipA_elas",
false);
802 inukd -> FracPipA_Inel () -> SaveAsROOT(filename,
"pipA_inel",
false);
803 inukd -> FracPipA_CEx () -> SaveAsROOT(filename,
"pipA_cex",
false);
804 inukd -> FracPipA_Abs () -> SaveAsROOT(filename,
"pipA_abs",
false);
805 inukd -> FracPimA_Tot () -> SaveAsROOT(filename,
"pimA_tot",
false);
806 inukd -> FracPimA_Elas () -> SaveAsROOT(filename,
"pimA_elas",
false);
807 inukd -> FracPimA_Inel () -> SaveAsROOT(filename,
"pimA_inel",
false);
808 inukd -> FracPimA_CEx () -> SaveAsROOT(filename,
"pimA_cex",
false);
809 inukd -> FracPimA_Abs () -> SaveAsROOT(filename,
"pimA_abs",
false);
810 inukd -> FracPi0A_Tot () -> SaveAsROOT(filename,
"pi0A_tot",
false);
811 inukd -> FracPi0A_Elas () -> SaveAsROOT(filename,
"pi0A_elas",
false);
812 inukd -> FracPi0A_Inel () -> SaveAsROOT(filename,
"pi0A_inel",
false);
813 inukd -> FracPi0A_CEx () -> SaveAsROOT(filename,
"pi0A_cex",
false);
814 inukd -> FracPi0A_Abs () -> SaveAsROOT(filename,
"pi0A_abs",
false);
820 <<
"Printing out INTRANUKE hadron x-sections";
827 <<
"\n hN mode x-sections:" 828 <<
"\n XSec[pn/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecPn_Tot () -> Evaluate(ke) <<
" mbarn" 829 <<
"\n XSec[pn/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecPn_Elas () -> Evaluate(ke) <<
" mbarn" 830 <<
"\n XSec[pn/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecPn_Reac () -> Evaluate(ke) <<
" mbarn" 831 <<
"\n XSec[nn/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecNn_Tot () -> Evaluate(ke) <<
" mbarn" 832 <<
"\n XSec[nn/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecNn_Elas () -> Evaluate(ke) <<
" mbarn" 833 <<
"\n XSec[nn/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecNn_Reac () -> Evaluate(ke) <<
" mbarn" 834 <<
"\n XSec[pp/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecPp_Tot () -> Evaluate(ke) <<
" mbarn" 835 <<
"\n XSec[pp/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecPp_Elas () -> Evaluate(ke) <<
" mbarn" 836 <<
"\n XSec[pp/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecPp_Reac () -> Evaluate(ke) <<
" mbarn" 837 <<
"\n XSec[pipn/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecPipn_Tot () -> Evaluate(ke) <<
" mbarn" 838 <<
"\n XSec[pipn/cex] (K=" << ke <<
" MeV) = " << inukd -> XSecPipn_CEx () -> Evaluate(ke) <<
" mbarn" 839 <<
"\n XSec[pipn/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecPipn_Elas () -> Evaluate(ke) <<
" mbarn" 840 <<
"\n XSec[pipn/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecPipn_Reac () -> Evaluate(ke) <<
" mbarn" 841 <<
"\n XSec[pi0n/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0n_Tot () -> Evaluate(ke) <<
" mbarn" 842 <<
"\n XSec[pi0n/cex] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0n_CEx () -> Evaluate(ke) <<
" mbarn" 843 <<
"\n XSec[pi0n/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0n_Elas () -> Evaluate(ke) <<
" mbarn" 844 <<
"\n XSec[pi0n/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0n_Reac () -> Evaluate(ke) <<
" mbarn" 845 <<
"\n XSec[pipp/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecPipp_Tot () -> Evaluate(ke) <<
" mbarn" 846 <<
"\n XSec[pipp/cex] (K=" << ke <<
" MeV) = " << inukd -> XSecPipp_CEx () -> Evaluate(ke) <<
" mbarn" 847 <<
"\n XSec[pipp/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecPipp_Elas () -> Evaluate(ke) <<
" mbarn" 848 <<
"\n XSec[pipp/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecPipp_Reac () -> Evaluate(ke) <<
" mbarn" 849 <<
"\n XSec[pipp/abs] (K=" << ke <<
" MeV) = " << inukd -> XSecPipp_Abs () -> Evaluate(ke) <<
" mbarn" 850 <<
"\n XSec[pi0p/tot] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0p_Tot () -> Evaluate(ke) <<
" mbarn" 851 <<
"\n XSec[pi0p/cex] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0p_CEx () -> Evaluate(ke) <<
" mbarn" 852 <<
"\n XSec[pi0p/elas] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0p_Elas () -> Evaluate(ke) <<
" mbarn" 853 <<
"\n XSec[pi0p/reac] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0p_Reac () -> Evaluate(ke) <<
" mbarn" 854 <<
"\n XSec[pi0p/abs] (K=" << ke <<
" MeV) = " << inukd -> XSecPi0p_Abs () -> Evaluate(ke) <<
" mbarn" 856 <<
"\n hA mode x-sections:" 857 <<
"\n Frac[pA/tot] (K=" << ke <<
" MeV) = " << inukd -> FracPA_Tot () -> Evaluate(ke) <<
" mbarn" 858 <<
"\n Frac[pA/elas] (K=" << ke <<
" MeV) = " << inukd -> FracPA_Elas () -> Evaluate(ke) <<
" mbarn" 859 <<
"\n Frac[pA/inel] (K=" << ke <<
" MeV) = " << inukd -> FracPA_Inel () -> Evaluate(ke) <<
" mbarn" 860 <<
"\n Frac[pA/cex] (K=" << ke <<
" MeV) = " << inukd -> FracPA_CEx () -> Evaluate(ke) <<
" mbarn" 861 <<
"\n Frac[pA/abs] (K=" << ke <<
" MeV) = " << inukd -> FracPA_Abs () -> Evaluate(ke) <<
" mbarn" 862 <<
"\n Frac[pA/pipro] (K=" << ke <<
" MeV) = " << inukd -> FracPA_Pipro () -> Evaluate(ke) <<
" mbarn" 863 <<
"\n Frac[nA/tot] (K=" << ke <<
" MeV) = " << inukd -> FracNA_Tot () -> Evaluate(ke) <<
" mbarn" 864 <<
"\n Frac[nA/elas] (K=" << ke <<
" MeV) = " << inukd -> FracNA_Elas () -> Evaluate(ke) <<
" mbarn" 865 <<
"\n Frac[nA/inel] (K=" << ke <<
" MeV) = " << inukd -> FracNA_Inel () -> Evaluate(ke) <<
" mbarn" 866 <<
"\n Frac[nA/reac] (K=" << ke <<
" MeV) = " << inukd -> FracNA_CEx () -> Evaluate(ke) <<
" mbarn" 867 <<
"\n Frac[nA/abs] (K=" << ke <<
" MeV) = " << inukd -> FracNA_Abs () -> Evaluate(ke) <<
" mbarn" 868 <<
"\n Frac[nA/pipro] (K=" << ke <<
" MeV) = " << inukd -> FracNA_Pipro () -> Evaluate(ke) <<
" mbarn" 869 <<
"\n Frac[pipA/tot] (K=" << ke <<
" MeV) = " << inukd -> FracPipA_Tot () -> Evaluate(ke) <<
" mbarn" 870 <<
"\n Frac[pipA/elas] (K=" << ke <<
" MeV) = " << inukd -> FracPipA_Elas () -> Evaluate(ke) <<
" mbarn" 871 <<
"\n Frac[pipA/inel] (K=" << ke <<
" MeV) = " << inukd -> FracPipA_Inel () -> Evaluate(ke) <<
" mbarn" 872 <<
"\n Frac[pipA/cex] (K=" << ke <<
" MeV) = " << inukd -> FracPipA_CEx () -> Evaluate(ke) <<
" mbarn" 873 <<
"\n Frac[pipA/abs] (K=" << ke <<
" MeV) = " << inukd -> FracPipA_Abs () -> Evaluate(ke) <<
" mbarn" 874 <<
"\n Frac[pimA/tot] (K=" << ke <<
" MeV) = " << inukd -> FracPimA_Tot () -> Evaluate(ke) <<
" mbarn" 875 <<
"\n Frac[pimA/elas] (K=" << ke <<
" MeV) = " << inukd -> FracPimA_Elas () -> Evaluate(ke) <<
" mbarn" 876 <<
"\n Frac[pimA/inel] (K=" << ke <<
" MeV) = " << inukd -> FracPimA_Inel () -> Evaluate(ke) <<
" mbarn" 877 <<
"\n Frac[pimA/cex] (K=" << ke <<
" MeV) = " << inukd -> FracPimA_CEx () -> Evaluate(ke) <<
" mbarn" 878 <<
"\n Frac[pimA/abs] (K=" << ke <<
" MeV) = " << inukd -> FracPimA_Abs () -> Evaluate(ke) <<
" mbarn" 879 <<
"\n Frac[pi0A/tot] (K=" << ke <<
" MeV) = " << inukd -> FracPi0A_Tot () -> Evaluate(ke) <<
" mbarn" 880 <<
"\n Frac[pi0A/elas] (K=" << ke <<
" MeV) = " << inukd -> FracPi0A_Elas () -> Evaluate(ke) <<
" mbarn" 881 <<
"\n Frac[pi0A/inel] (K=" << ke <<
" MeV) = " << inukd -> FracPi0A_Inel () -> Evaluate(ke) <<
" mbarn" 882 <<
"\n Frac[pi0A/cex] (K=" << ke <<
" MeV) = " << inukd -> FracPi0A_CEx () -> Evaluate(ke) <<
" mbarn" 883 <<
"\n Frac[pi0A/abs] (K=" << ke <<
" MeV) = " << inukd -> FracPi0A_Abs () -> Evaluate(ke) <<
" mbarn";
THE MAIN GENIE PROJECT NAMESPACE
void PrintOutForInputKE(double ke)
int main(int argc, char **argv)
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
static INukeHadroData * Instance(void)
void SaveDataToRootFile(bool test_intbounce)
void SaveSplinesToRootFile(void)
Singleton class to load & serve hadron x-section splines used by GENIE's version of the INTRANUKE cascade...
Command line argument parser.
Most commonly used PDG codes. A set of utility functions to handle PDG codes is provided in PDGUtils...
bool OptionExists(char opt)
was option set?
QTextStream & endl(QTextStream &s)