Changeset 1055 for trunk/source/processes/hadronic/models/chiral_inv_phase_space/interface/src/G4QANuENuclearCrossSection.cc
- Timestamp:
- May 28, 2009, 4:26:57 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/chiral_inv_phase_space/interface/src/G4QANuENuclearCrossSection.cc
r1007 r1055 25 25 // 26 26 // 27 // $Id: G4QANuENuclearCrossSection.cc,v 1. 2 2007/11/01 16:09:38mkossov Exp $28 // GEANT4 tag $Name: geant4-09-0 2$27 // $Id: G4QANuENuclearCrossSection.cc,v 1.4 2009/05/08 15:16:26 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03-beta-cand-01 $ 29 29 // 30 30 // 31 // G4 Physics class: G4QANuENuclearCrossSection for (anti_n i_e,e+) cross sections31 // G4 Physics class: G4QANuENuclearCrossSection for (anti_nu_e,e+) cross sections 32 32 // Created: M.V. Kossov, CERN/ITEP(Moscow), 24-SEP-2007 33 33 // The last update: M.V. Kossov, CERN/ITEP (Moscow) 24-SEP-2007 … … 37 37 // ******* DO NOT MAKE ANY CHANGE! With time it'll move back to photolepton...(M.K.) ****** 38 38 // **************************************************************************************** 39 //========================================================================================= ======39 //========================================================================================= 40 40 41 41 //#define debug … … 90 90 <<"("<<lastN<<"),PDG="<<pPDG<<"("<<lastPDG<<"), T="<<pEn<<"("<<lastTH<<")"<<",Sz=" 91 91 <<colN.size()<<G4endl; 92 92 //CalculateCrossSection(fCS,-27,j,lastPDG,lastZ,lastN,pMom); // DUMMY TEST 93 93 #endif 94 94 if(pPDG!=-12) … … 111 111 j = 0; // A#0f records found in DB for this projectile 112 112 if(lastI) for(G4int i=0; i<lastI; i++) if(colPDG[i]==pPDG) // The partType is found 113 113 { // The nucleus with projPDG is found in AMDB 114 114 if(colN[i]==tgN && colZ[i]==tgZ) 115 115 { 116 116 lastI=i; 117 117 lastTH =colTH[i]; // Last THreshold (A-dependent) … … 163 163 #endif 164 164 j++; // Increment a#0f records found in DB for this pPDG 165 166 167 165 } 166 if(!in) // This nucleus has not been calculated previously 167 { 168 168 #ifdef pdebug 169 169 G4cout<<"G4QAENCS::GetCrosSec:CalcNew P="<<pMom<<",f="<<fCS<<",lstI="<<lastI<<G4endl; … … 172 172 lastCS=CalculateCrossSection(fCS,0,j,lastPDG,lastZ,lastN,pMom); //calculate & create 173 173 if(lastCS<=0.) 174 174 { 175 175 lastTH = ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last 176 176 #ifdef pdebug … … 184 184 lastTH=pEn; 185 185 } 186 186 } 187 187 #ifdef pdebug 188 188 G4cout<<"G4QAENCS::GetCrosSec:New CS="<<lastCS<<",lZ="<<lastN<<",lN="<<lastZ<<G4endl; … … 200 200 #endif 201 201 return lastCS*millibarn; 202 202 } // End of creation of the new set of parameters 203 203 else 204 204 { 205 205 #ifdef pdebug 206 206 G4cout<<"G4QAENCS::GetCS: Update lastI="<<lastI<<",j="<<j<<G4endl; … … 245 245 G4double G4QANuENuclearCrossSection::ThresholdEnergy(G4int Z, G4int N, G4int) 246 246 { 247 //static const G4double mNeut = G4NucleiProperties::GetNuclearMass(1 .,0.)/GeV;248 //static const G4double mProt = G4NucleiProperties::GetNuclearMass(1 .,1.)/GeV;249 //static const G4double mDeut = G4NucleiProperties::GetNuclearMass(2 .,1.)/GeV/2.;247 //static const G4double mNeut = G4NucleiProperties::GetNuclearMass(1,0)/GeV; 248 //static const G4double mProt = G4NucleiProperties::GetNuclearMass(1,1)/GeV; 249 //static const G4double mDeut = G4NucleiProperties::GetNuclearMass(2,1)/GeV/2.; 250 250 static const G4double mN=.931494043;// Nucleon mass (inside nucleus, AtomicMassUnit, GeV) 251 251 static const G4double dmN=mN+mN; // Doubled nucleon mass (2*AtomicMassUnit, GeV) … … 294 294 { 295 295 if(F<0) // This isotope was found in DAMDB =========> RETRIEVE 296 296 { 297 297 lastTX =TX[I]; // Pointer to the prepared TX function (same isotope) 298 298 lastQE =QE[I]; // Pointer to the prepared QE function (same isotope) 299 300 301 299 } 300 else // This isotope wasn't calculated previously => CREATE 301 { 302 302 if(first) 303 303 { … … 315 315 TX.push_back(lastTX); 316 316 QE.push_back(lastQE); 317 317 } // End of creation of the new set of parameters 318 318 } // End of parameters udate 319 319 // ============================== NOW Calculate the Cross Section ===================== … … 330 330 G4int sep=ran; // as a result = an index of the left edge of the interval 331 331 while(ran>=2) 332 332 { 333 333 G4int newran=ran/2; 334 334 G4double oldL=lastEN[sep]; … … 386 386 static const G4int nE=65; // !! If change this, change it in CalculateCrossSection() !! 387 387 static const G4double nuEn[nE]={thresh, 388 .00051331,.00053602,.00056078,.00058783,.00061743,.00064990,.00068559,.00072492,389 .00076834,.00081641,.00086975,.00092912,.00099536,.00106950,.00115273,.00124646,390 .00135235,.00147241,.00160901,.00176503,.00194392,.00214986,.00238797,.00266448,391 .00298709,.00336531,.00381094,.00433879,.00496745,.00572047,.00662785,.00772806,392 .00907075,.01072050,.01276190,.01530660,.01850330,.02255110,.02771990,.03437780,393 .04303240,.05438970,.06944210,.08959920,.11688400,.15423600,.20597200,.27851200,394 .38153100,.52979600,.74616300,1.0665200,1.5480900,2.2834800,3.4251100,5.2281000,395 8.1270200,12.875900,20.808500,34.331200,57.877800,99.796200,176.16300,318.68200};388 .00051331,.00053602,.00056078,.00058783,.00061743,.00064990,.00068559,.00072492, 389 .00076834,.00081641,.00086975,.00092912,.00099536,.00106950,.00115273,.00124646, 390 .00135235,.00147241,.00160901,.00176503,.00194392,.00214986,.00238797,.00266448, 391 .00298709,.00336531,.00381094,.00433879,.00496745,.00572047,.00662785,.00772806, 392 .00907075,.01072050,.01276190,.01530660,.01850330,.02255110,.02771990,.03437780, 393 .04303240,.05438970,.06944210,.08959920,.11688400,.15423600,.20597200,.27851200, 394 .38153100,.52979600,.74616300,1.0665200,1.5480900,2.2834800,3.4251100,5.2281000, 395 8.1270200,12.875900,20.808500,34.331200,57.877800,99.796200,176.16300,318.68200}; 396 396 static const G4double TOTX[nE]={0., 397 .00046923,.00160693,.00229560,.00288772,.00344344,.00398831,.00453725,.00510087,398 .00568797,.00630663,.00696489,.00767121,.00843477,.00926586,.01017610,.01117910,399 .01229040,.01352820,.01491430,.01647420,.01823820,.02024280,.02253130,.02515600,400 .02818010,.03167950,.03574640,.04049260,.04605330,.05739330,.07028190,.08396290,401 .09969400,.11756200,.13812500,.16207000,.19099600,.22447700,.26434900,.31051900,402 .36357400,.42242900,.48500000,.54692800,.60140500,.63996800,.65519300,.64339200,403 .60784300,.55734400,.50212600,.45116800,.40962600,.37854800,.35702100,.34330900,404 .33576200,.33300100,.33387100,.33737700,.34235200,.34880100,.35507500,.35961200};397 .00046923,.00160693,.00229560,.00288772,.00344344,.00398831,.00453725,.00510087, 398 .00568797,.00630663,.00696489,.00767121,.00843477,.00926586,.01017610,.01117910, 399 .01229040,.01352820,.01491430,.01647420,.01823820,.02024280,.02253130,.02515600, 400 .02818010,.03167950,.03574640,.04049260,.04605330,.05739330,.07028190,.08396290, 401 .09969400,.11756200,.13812500,.16207000,.19099600,.22447700,.26434900,.31051900, 402 .36357400,.42242900,.48500000,.54692800,.60140500,.63996800,.65519300,.64339200, 403 .60784300,.55734400,.50212600,.45116800,.40962600,.37854800,.35702100,.34330900, 404 .33576200,.33300100,.33387100,.33737700,.34235200,.34880100,.35507500,.35961200}; 405 405 static const G4double QELX[nE]={0., 406 2.40856e-7,8.61338e-7,1.28732e-6,1.69747e-6,2.12608e-6,2.59201e-6,3.11071e-6,3.69770e-6,407 4.37028e-6,5.14877e-6,6.05773e-6,7.12746e-6,8.39567e-6,9.90987e-6,1.17304e-5,1.39343e-5,408 1.66209e-5,1.99191e-5,2.39974e-5,2.90775e-5,3.54536e-5,4.35191e-5,5.38039e-5,6.70278e-5,409 8.41765e-5,.000106611,.000136228,.000175689,.000228768,.000328317,.000465818,.000648870,410 .000904299,.001260330,.001762740,.002480740,.003534040,.005062210,.007327720,.010675000,411 .015645500,.022975800,.033679400,.049004300,.070294800,.098706100,.134951000,.179193000,412 .231297000,.287287000,.344760000,.404397000,.465915000,.527060000,.584085000,.632748000,413 .671346000,.699744000,.719355000,.732541000,.740996000,.746249000,.749265000,.751160000};406 2.40856e-7,8.61338e-7,1.28732e-6,1.69747e-6,2.12608e-6,2.59201e-6,3.11071e-6,3.69770e-6, 407 4.37028e-6,5.14877e-6,6.05773e-6,7.12746e-6,8.39567e-6,9.90987e-6,1.17304e-5,1.39343e-5, 408 1.66209e-5,1.99191e-5,2.39974e-5,2.90775e-5,3.54536e-5,4.35191e-5,5.38039e-5,6.70278e-5, 409 8.41765e-5,.000106611,.000136228,.000175689,.000228768,.000328317,.000465818,.000648870, 410 .000904299,.001260330,.001762740,.002480740,.003534040,.005062210,.007327720,.010675000, 411 .015645500,.022975800,.033679400,.049004300,.070294800,.098706100,.134951000,.179193000, 412 .231297000,.287287000,.344760000,.404397000,.465915000,.527060000,.584085000,.632748000, 413 .671346000,.699744000,.719355000,.732541000,.740996000,.746249000,.749265000,.751160000}; 414 414 // -------------------------------- 415 415 G4int first=0; 416 416 if(z<0.) 417 417 { 418 418 first=1; 419 419 z=-z; … … 434 434 t[k]=TOTX[k]*nuEn[k]*(za+za)/ta+QELX[k]*(dz+dz-da)/ta; // TotalCrossSection 435 435 q[k]=QELX[k]*dz/a; // QuasiElasticCrossSection 436 436 } 437 437 return first; 438 438 } … … 454 454 static const G4double Xl[nQ2]={5.20224e-16, 455 455 .006125,.0137008,.0218166,.0302652,.0389497,.0478144,.0568228,.0659497,.0751768,.0844898, 456 .093878, .103332, .112844, .122410, .132023, .141680, .151376, .161109, .170875,.180672,457 .190499, .200352, .210230, .220131, .230055, .239999, .249963, .259945, .269944,.279960,458 .289992, .300039, .310099, .320173, .330260, .340359, .350470, .360592, .370724,.380867,459 .391019, .401181, .411352, .421531, .431719, .441915, .452118, .462329, .472547,.482771,460 .493003, .503240, .513484, .523734, .533989, .544250, .554517, .564788, .575065,.585346,461 .595632, .605923, .616218, .626517, .636820, .647127, .657438, .667753, .678072,.688394,462 .698719, .709048, .719380, .729715, .740053, .750394, .760738, .771085, .781434,.791786,463 .802140, .812497, .822857, .833219, .843582, .853949, .864317, .874687, .885060,.895434,464 .905810, .916188, .926568, .936950, .947333, .957719, .968105, .978493, .988883,.999275};465 456 .093878, .103332, .112844, .122410, .132023, .141680, .151376, .161109, .170875, .180672, 457 .190499, .200352, .210230, .220131, .230055, .239999, .249963, .259945, .269944, .279960, 458 .289992, .300039, .310099, .320173, .330260, .340359, .350470, .360592, .370724, .380867, 459 .391019, .401181, .411352, .421531, .431719, .441915, .452118, .462329, .472547, .482771, 460 .493003, .503240, .513484, .523734, .533989, .544250, .554517, .564788, .575065, .585346, 461 .595632, .605923, .616218, .626517, .636820, .647127, .657438, .667753, .678072, .688394, 462 .698719, .709048, .719380, .729715, .740053, .750394, .760738, .771085, .781434, .791786, 463 .802140, .812497, .822857, .833219, .843582, .853949, .864317, .874687, .885060, .895434, 464 .905810, .916188, .926568, .936950, .947333, .957719, .968105, .978493, .988883, .999275}; 465 // Direct table 466 466 static const G4double Xmax=Xl[lQ2]; 467 467 static const G4double Xmin=Xl[0]; 468 468 static const G4double dX=(Xmax-Xmin)/lQ2; // step in X(Q2, GeV^2) 469 469 static const G4double inl[nQ2]={0, 470 1.52225, 2.77846, 3.96651, 5.11612, 6.23990, 7.34467, 8.43466, 9.51272, 10.5809,11.6406,471 12.6932, 13.7394, 14.7801, 15.8158, 16.8471, 17.8743, 18.8979, 19.9181, 20.9353,21.9496,472 22.9614, 23.9707, 24.9777, 25.9826, 26.9855, 27.9866, 28.9860, 29.9837, 30.9798,31.9745,473 32.9678, 33.9598, 34.9505, 35.9400, 36.9284, 37.9158, 38.9021, 39.8874, 40.8718,41.8553,474 42.8379, 43.8197, 44.8007, 45.7810, 46.7605, 47.7393, 48.7174, 49.6950, 50.6718,51.6481,475 52.6238, 53.5990, 54.5736, 55.5476, 56.5212, 57.4943, 58.4670, 59.4391, 60.4109,61.3822,476 62.3531, 63.3236, 64.2937, 65.2635, 66.2329, 67.2019, 68.1707, 69.1390, 70.1071,71.0748,477 72.0423, 73.0095, 73.9763, 74.9429, 75.9093, 76.8754, 77.8412, 78.8068, 79.7721,80.7373,478 81.7022, 82.6668, 83.6313, 84.5956, 85.5596, 86.5235, 87.4872, 88.4507, 89.4140,90.3771,479 91.3401, 92.3029, 93.2656, 94.2281, 95.1904, 96.1526, 97.1147, 98.0766, 99.0384,100.000};470 1.52225, 2.77846, 3.96651, 5.11612, 6.23990, 7.34467, 8.43466, 9.51272, 10.5809, 11.6406, 471 12.6932, 13.7394, 14.7801, 15.8158, 16.8471, 17.8743, 18.8979, 19.9181, 20.9353, 21.9496, 472 22.9614, 23.9707, 24.9777, 25.9826, 26.9855, 27.9866, 28.9860, 29.9837, 30.9798, 31.9745, 473 32.9678, 33.9598, 34.9505, 35.9400, 36.9284, 37.9158, 38.9021, 39.8874, 40.8718, 41.8553, 474 42.8379, 43.8197, 44.8007, 45.7810, 46.7605, 47.7393, 48.7174, 49.6950, 50.6718, 51.6481, 475 52.6238, 53.5990, 54.5736, 55.5476, 56.5212, 57.4943, 58.4670, 59.4391, 60.4109, 61.3822, 476 62.3531, 63.3236, 64.2937, 65.2635, 66.2329, 67.2019, 68.1707, 69.1390, 70.1071, 71.0748, 477 72.0423, 73.0095, 73.9763, 74.9429, 75.9093, 76.8754, 77.8412, 78.8068, 79.7721, 80.7373, 478 81.7022, 82.6668, 83.6313, 84.5956, 85.5596, 86.5235, 87.4872, 88.4507, 89.4140, 90.3771, 479 91.3401, 92.3029, 93.2656, 94.2281, 95.1904, 96.1526, 97.1147, 98.0766, 99.0384, 100.000}; 480 480 G4double Enu=lastE; // Get energy of the last calculated cross-section 481 481 G4double dEnu=Enu+Enu; // doubled energy of nu/anu … … 544 544 // Reversed table 545 545 static const G4double X0[nX]={5.21412e-05, 546 .437860, .681908, .891529, 1.08434, 1.26751, 1.44494, 1.61915, 1.79198, 1.96493,2.13937,547 2.31664, 2.49816, 2.68559, 2.88097, 3.08705, 3.30774, 3.54917, 3.82233,4.15131, 4.62182};546 .437860, .681908, .891529, 1.08434, 1.26751, 1.44494, 1.61915, 1.79198, 1.96493, 2.13937, 547 2.31664, 2.49816, 2.68559, 2.88097, 3.08705, 3.30774, 3.54917, 3.82233, 4.15131, 4.62182}; 548 548 static const G4double X1[nX]={.00102591, 549 1.00443, 1.55828, 2.03126, 2.46406, 2.87311, 3.26723, 3.65199, 4.03134, 4.40835,4.78561,550 5.16549, 5.55031, 5.94252, 6.34484, 6.76049, 7.19349, 7.64917, 8.13502,8.66246, 9.25086};549 1.00443, 1.55828, 2.03126, 2.46406, 2.87311, 3.26723, 3.65199, 4.03134, 4.40835, 4.78561, 550 5.16549, 5.55031, 5.94252, 6.34484, 6.76049, 7.19349, 7.64917, 8.13502, 8.66246, 9.25086}; 551 551 static const G4double X2[nX]={.0120304, 552 2.59903, 3.98637, 5.15131, 6.20159, 7.18024, 8.10986, 9.00426, 9.87265, 10.7217,11.5564,553 12.3808, 13.1983, 14.0116, 14.8234, 15.6359, 16.4515, 17.2723, 18.1006,18.9386, 19.7892};552 2.59903, 3.98637, 5.15131, 6.20159, 7.18024, 8.10986, 9.00426, 9.87265, 10.7217, 11.5564, 553 12.3808, 13.1983, 14.0116, 14.8234, 15.6359, 16.4515, 17.2723, 18.1006, 18.9386, 19.7892}; 554 554 static const G4double X3[nX]={.060124, 555 5.73857, 8.62595, 10.9849, 13.0644, 14.9636, 16.7340, 18.4066, 20.0019, 21.5342,23.0142,556 24.4497, 25.8471, 27.2114, 28.5467, 29.8564, 31.1434, 32.4102, 33.6589,34.8912, 36.1095};555 5.73857, 8.62595, 10.9849, 13.0644, 14.9636, 16.7340, 18.4066, 20.0019, 21.5342, 23.0142, 556 24.4497, 25.8471, 27.2114, 28.5467, 29.8564, 31.1434, 32.4102, 33.6589, 34.8912, 36.1095}; 557 557 static const G4double X4[nX]={.0992363, 558 8.23746, 12.1036, 15.1740, 17.8231, 20.1992, 22.3792, 24.4092, 26.3198, 28.1320,29.8615,559 31.5200, 33.1169, 34.6594, 36.1536, 37.6044, 39.0160, 40.3920, 41.7353,43.0485, 44.3354};558 8.23746, 12.1036, 15.1740, 17.8231, 20.1992, 22.3792, 24.4092, 26.3198, 28.1320, 29.8615, 559 31.5200, 33.1169, 34.6594, 36.1536, 37.6044, 39.0160, 40.3920, 41.7353, 43.0485, 44.3354}; 560 560 static const G4double X5[nX]={.0561127, 561 7.33661, 10.5694, 13.0778, 15.2061, 17.0893, 18.7973, 20.3717, 21.8400, 23.2211,24.5291,562 25.7745, 26.9655, 28.1087, 29.2094, 30.2721, 31.3003, 32.2972, 33.2656,34.2076, 35.1265};561 7.33661, 10.5694, 13.0778, 15.2061, 17.0893, 18.7973, 20.3717, 21.8400, 23.2211, 24.5291, 562 25.7745, 26.9655, 28.1087, 29.2094, 30.2721, 31.3003, 32.2972, 33.2656, 34.2076, 35.1265}; 563 563 static const G4double X6[nX]={.0145859, 564 4.81774, 6.83565, 8.37399, 9.66291, 10.7920, 11.8075, 12.7366, 13.5975, 14.4025,15.1608,565 15.8791, 16.5628, 17.2162, 17.8427, 18.4451, 19.0259, 19.5869, 20.1300,20.6566, 21.1706};564 4.81774, 6.83565, 8.37399, 9.66291, 10.7920, 11.8075, 12.7366, 13.5975, 14.4025, 15.1608, 565 15.8791, 16.5628, 17.2162, 17.8427, 18.4451, 19.0259, 19.5869, 20.1300, 20.6566, 21.1706}; 566 566 static const G4double X7[nX]={.00241155, 567 2.87095, 4.02492, 4.89243, 5.61207, 6.23747, 6.79613, 7.30433, 7.77270, 8.20858,8.61732,568 9.00296, 9.36863, 9.71682, 10.0495, 10.3684, 10.6749, 10.9701, 11.2550,11.5306, 11.7982};567 2.87095, 4.02492, 4.89243, 5.61207, 6.23747, 6.79613, 7.30433, 7.77270, 8.20858, 8.61732, 568 9.00296, 9.36863, 9.71682, 10.0495, 10.3684, 10.6749, 10.9701, 11.2550, 11.5306, 11.7982}; 569 569 static const G4double X8[nX]={.000316863, 570 1.76189, 2.44632, 2.95477, 3.37292, 3.73378, 4.05420, 4.34415, 4.61009, 4.85651,5.08666,571 5.30299, 5.50738, 5.70134, 5.88609, 6.06262, 6.23178, 6.39425, 6.55065,6.70149, 6.84742};570 1.76189, 2.44632, 2.95477, 3.37292, 3.73378, 4.05420, 4.34415, 4.61009, 4.85651, 5.08666, 571 5.30299, 5.50738, 5.70134, 5.88609, 6.06262, 6.23178, 6.39425, 6.55065, 6.70149, 6.84742}; 572 572 static const G4double X9[nX]={3.73544e-05, 573 1.17106, 1.61289, 1.93763, 2.20259, 2.42976, 2.63034, 2.81094, 2.97582, 3.12796,3.26949,574 3.40202, 3.52680, 3.64482, 3.75687, 3.86360, 3.96557, 4.06323, 4.15697,4.24713, 4.33413};573 1.17106, 1.61289, 1.93763, 2.20259, 2.42976, 2.63034, 2.81094, 2.97582, 3.12796, 3.26949, 574 3.40202, 3.52680, 3.64482, 3.75687, 3.86360, 3.96557, 4.06323, 4.15697, 4.24713, 4.33413}; 575 575 static const G4double XA[nX]={4.19131e-06, 576 .849573, 1.16208, 1.38955, 1.57379, 1.73079, 1.86867, 1.99221, 2.10451, 2.20770,2.30332,577 2.39252, 2.47622, 2.55511, 2.62977, 2.70066, 2.76818, 2.83265, 2.89437,2.95355, 3.01051};576 .849573, 1.16208, 1.38955, 1.57379, 1.73079, 1.86867, 1.99221, 2.10451, 2.20770, 2.30332, 577 2.39252, 2.47622, 2.55511, 2.62977, 2.70066, 2.76818, 2.83265, 2.89437, 2.95355, 3.01051}; 578 578 static const G4double XB[nX]={4.59981e-07, 579 .666131, .905836, 1.07880, 1.21796, 1.33587, 1.43890, 1.53080, 1.61399, 1.69011,1.76040,580 1.82573, 1.88682, 1.94421, 1.99834, 2.04959, 2.09824, 2.14457, 2.18878,2.23107, 2.27162};579 .666131, .905836, 1.07880, 1.21796, 1.33587, 1.43890, 1.53080, 1.61399, 1.69011, 1.76040, 580 1.82573, 1.88682, 1.94421, 1.99834, 2.04959, 2.09824, 2.14457, 2.18878, 2.23107, 2.27162}; 581 581 static const G4double XC[nX]={4.99861e-08, 582 .556280, .752730, .893387, 1.00587, 1.10070, 1.18317, 1.25643, 1.32247, 1.38269,1.43809,583 1.48941, 1.53724, 1.58203, 1.62416, 1.66391, 1.70155, 1.73728, 1.77128,1.80371, 1.83473};582 .556280, .752730, .893387, 1.00587, 1.10070, 1.18317, 1.25643, 1.32247, 1.38269, 1.43809, 583 1.48941, 1.53724, 1.58203, 1.62416, 1.66391, 1.70155, 1.73728, 1.77128, 1.80371, 1.83473}; 584 584 static const G4double XD[nX]={5.40832e-09, 585 .488069, .657650, .778236, .874148, .954621, 1.02432, 1.08599, 1.14138, 1.19172,1.23787,586 1.28049, 1.32008, 1.35705, 1.39172, 1.42434, 1.45514, 1.48429, 1.51197,1.53829, 1.56339};585 .488069, .657650, .778236, .874148, .954621, 1.02432, 1.08599, 1.14138, 1.19172, 1.23787, 586 1.28049, 1.32008, 1.35705, 1.39172, 1.42434, 1.45514, 1.48429, 1.51197, 1.53829, 1.56339}; 587 587 static const G4double XE[nX]={5.84029e-10, 588 .445057, .597434, .705099, .790298, .861468, .922865, .976982, 1.02542, 1.06930,1.10939,589 1.14630, 1.18050, 1.21233, 1.24208, 1.27001, 1.29630, 1.32113, 1.34462,1.36691, 1.38812};588 .445057, .597434, .705099, .790298, .861468, .922865, .976982, 1.02542, 1.06930, 1.10939, 589 1.14630, 1.18050, 1.21233, 1.24208, 1.27001, 1.29630, 1.32113, 1.34462, 1.36691, 1.38812}; 590 590 static const G4double XF[nX]={6.30137e-11, 591 .418735, .560003, .659168, .737230, .802138, .857898, .906854, .950515, .989915,1.02580,592 1.05873, 1.08913, 1.11734, 1.14364, 1.16824, 1.19133, 1.21306, 1.23358,1.25298, 1.27139};591 .418735, .560003, .659168, .737230, .802138, .857898, .906854, .950515, .989915, 1.02580, 592 1.05873, 1.08913, 1.11734, 1.14364, 1.16824, 1.19133, 1.21306, 1.23358, 1.25298, 1.27139}; 593 593 static const G4double XG[nX]={6.79627e-12, 594 .405286, .539651, .633227, .706417, .766929, .818642, .863824, .903931, .939963,.972639,595 1.00250, 1.02995, 1.05532, 1.07887, 1.10082, 1.12134, 1.14058, 1.15867,1.17572, 1.19183};594 .405286, .539651, .633227, .706417, .766929, .818642, .863824, .903931, .939963, .972639, 595 1.00250, 1.02995, 1.05532, 1.07887, 1.10082, 1.12134, 1.14058, 1.15867, 1.17572, 1.19183}; 596 596 static const G4double XH[nX]={7.32882e-13, 597 .404391, .535199, .625259, .695036, .752243, .800752, .842823, .879906, .912994,.942802,598 .969862, .994583, 1.01729, 1.03823, 1.05763, 1.07566, 1.09246, 1.10816,1.12286, 1.13667};597 .404391, .535199, .625259, .695036, .752243, .800752, .842823, .879906, .912994, .942802, 598 .969862, .994583, 1.01729, 1.03823, 1.05763, 1.07566, 1.09246, 1.10816, 1.12286, 1.13667}; 599 599 static const G4double XI[nX]={7.90251e-14, 600 .418084, .548382, .636489, .703728, .758106, .803630, .842633, .876608, .906576,.933269,601 .957233, .978886, .998556, 1.01651, 1.03295, 1.04807, 1.06201, 1.07489,1.08683, 1.09792};600 .418084, .548382, .636489, .703728, .758106, .803630, .842633, .876608, .906576, .933269, 601 .957233, .978886, .998556, 1.01651, 1.03295, 1.04807, 1.06201, 1.07489, 1.08683, 1.09792}; 602 602 static const G4double XJ[nX]={8.52083e-15, 603 .447299, .579635, .666780, .731788, .783268, .825512, .861013, .891356, .917626,.940597,604 .960842, .978802, .994820, 1.00917, 1.02208, 1.03373, 1.04427, 1.05383,1.06253, 1.07046};603 .447299, .579635, .666780, .731788, .783268, .825512, .861013, .891356, .917626, .940597, 604 .960842, .978802, .994820, 1.00917, 1.02208, 1.03373, 1.04427, 1.05383, 1.06253, 1.07046}; 605 605 // Direct table 606 606 static const G4double Xmin[nE]={X0[0],X1[0],X2[0],X3[0],X4[0],X5[0],X6[0],X7[0],X8[0], … … 615 615 {X0,X1,X2,X3,X4,X5,X6,X7,X8,X9,XA,XB,XC,XD,XE,XF,XG,XH,XI,XJ}; 616 616 static const G4double I0[nX]={0, 617 .354631, 1.08972, 2.05138, 3.16564, 4.38343, 5.66828, 6.99127, 8.32858, 9.65998,10.9680,618 12.2371, 13.4536, 14.6050, 15.6802, 16.6686, 17.5609, 18.3482, 19.0221, 19.5752,20.0000};617 .354631, 1.08972, 2.05138, 3.16564, 4.38343, 5.66828, 6.99127, 8.32858, 9.65998, 10.9680, 618 12.2371, 13.4536, 14.6050, 15.6802, 16.6686, 17.5609, 18.3482, 19.0221, 19.5752, 20.0000}; 619 619 static const G4double I1[nX]={0, 620 .281625, .877354, 1.67084, 2.60566, 3.64420, 4.75838, 5.92589, 7.12829, 8.34989,9.57708,621 10.7978, 12.0014, 13.1781, 14.3190, 15.4162, 16.4620, 17.4496, 18.3724, 19.2245,20.0000};620 .281625, .877354, 1.67084, 2.60566, 3.64420, 4.75838, 5.92589, 7.12829, 8.34989, 9.57708, 621 10.7978, 12.0014, 13.1781, 14.3190, 15.4162, 16.4620, 17.4496, 18.3724, 19.2245, 20.0000}; 622 622 static const G4double I2[nX]={0, 623 .201909, .642991, 1.24946, 1.98463, 2.82370, 3.74802, 4.74263, 5.79509, 6.89474,8.03228,624 9.19947, 10.3889, 11.5938, 12.8082, 14.0262, 15.2427, 16.4527, 17.6518, 18.8356,20.0000};623 .201909, .642991, 1.24946, 1.98463, 2.82370, 3.74802, 4.74263, 5.79509, 6.89474, 8.03228, 624 9.19947, 10.3889, 11.5938, 12.8082, 14.0262, 15.2427, 16.4527, 17.6518, 18.8356, 20.0000}; 625 625 static const G4double I3[nX]={0, 626 .140937, .461189, .920216, 1.49706, 2.17728, 2.94985, 3.80580, 4.73758, 5.73867,6.80331,627 7.92637, 9.10316, 10.3294, 11.6013, 12.9150, 14.2672, 15.6548, 17.0746, 18.5239,20.0000};626 .140937, .461189, .920216, 1.49706, 2.17728, 2.94985, 3.80580, 4.73758, 5.73867, 6.80331, 627 7.92637, 9.10316, 10.3294, 11.6013, 12.9150, 14.2672, 15.6548, 17.0746, 18.5239, 20.0000}; 628 628 static const G4double I4[nX]={0, 629 .099161, .337358, .694560, 1.16037, 1.72761, 2.39078, 3.14540, 3.98768, 4.91433,5.92245,630 7.00942, 8.17287, 9.41060, 10.7206, 12.1010, 13.5500, 15.0659, 16.6472, 18.2924,20.0000};629 .099161, .337358, .694560, 1.16037, 1.72761, 2.39078, 3.14540, 3.98768, 4.91433, 5.92245, 630 7.00942, 8.17287, 9.41060, 10.7206, 12.1010, 13.5500, 15.0659, 16.6472, 18.2924, 20.0000}; 631 631 static const G4double I5[nX]={0, 632 .071131, .255084, .543312, .932025, 1.41892, 2.00243, 2.68144, 3.45512, 4.32283,5.28411,633 6.33859, 7.48602, 8.72621, 10.0590, 11.4844, 13.0023, 14.6128, 16.3158, 18.1115,20.0000};632 .071131, .255084, .543312, .932025, 1.41892, 2.00243, 2.68144, 3.45512, 4.32283, 5.28411, 633 6.33859, 7.48602, 8.72621, 10.0590, 11.4844, 13.0023, 14.6128, 16.3158, 18.1115, 20.0000}; 634 634 static const G4double I6[nX]={0, 635 .053692, .202354, .443946, .778765, 1.20774, 1.73208, 2.35319, 3.07256, 3.89177,4.81249,636 5.83641, 6.96528, 8.20092, 9.54516, 10.9999, 12.5670, 14.2486, 16.0466, 17.9630,20.0000};635 .053692, .202354, .443946, .778765, 1.20774, 1.73208, 2.35319, 3.07256, 3.89177, 4.81249, 636 5.83641, 6.96528, 8.20092, 9.54516, 10.9999, 12.5670, 14.2486, 16.0466, 17.9630, 20.0000}; 637 637 static const G4double I7[nX]={0, 638 .043065, .168099, .376879, .672273, 1.05738, 1.53543, 2.10973, 2.78364, 3.56065,4.44429,639 5.43819, 6.54610, 7.77186, 9.11940, 10.5928, 12.1963, 13.9342, 15.8110, 17.8313,20.0000};638 .043065, .168099, .376879, .672273, 1.05738, 1.53543, 2.10973, 2.78364, 3.56065, 4.44429, 639 5.43819, 6.54610, 7.77186, 9.11940, 10.5928, 12.1963, 13.9342, 15.8110, 17.8313, 20.0000}; 640 640 static const G4double I8[nX]={0, 641 .036051, .143997, .327877, .592202, .941572, 1.38068, 1.91433, 2.54746, 3.28517,4.13277,642 5.09574, 6.17984, 7.39106, 8.73568, 10.2203, 11.8519, 13.6377, 15.5854, 17.7033,20.0000};641 .036051, .143997, .327877, .592202, .941572, 1.38068, 1.91433, 2.54746, 3.28517, 4.13277, 642 5.09574, 6.17984, 7.39106, 8.73568, 10.2203, 11.8519, 13.6377, 15.5854, 17.7033, 20.0000}; 643 643 static const G4double I9[nX]={0, 644 .030977, .125727, .289605, .528146, .846967, 1.25183, 1.74871, 2.34384, 3.04376,3.85535,645 4.78594, 5.84329, 7.03567, 8.37194, 9.86163, 11.5150, 13.3430, 15.3576, 17.5719,20.0000};644 .030977, .125727, .289605, .528146, .846967, 1.25183, 1.74871, 2.34384, 3.04376, 3.85535, 645 4.78594, 5.84329, 7.03567, 8.37194, 9.86163, 11.5150, 13.3430, 15.3576, 17.5719, 20.0000}; 646 646 static const G4double IA[nX]={0, 647 .027129, .111420, .258935, .475812, .768320, 1.14297, 1.60661, 2.16648, 2.83034,3.60650,648 4.50394, 5.53238, 6.70244, 8.02569, 9.51488, 11.1841, 13.0488, 15.1264, 17.4362,20.0000};647 .027129, .111420, .258935, .475812, .768320, 1.14297, 1.60661, 2.16648, 2.83034, 3.60650, 648 4.50394, 5.53238, 6.70244, 8.02569, 9.51488, 11.1841, 13.0488, 15.1264, 17.4362, 20.0000}; 649 649 static const G4double IB[nX]={0, 650 .024170, .100153, .234345, .433198, .703363, 1.05184, 1.48607, 2.01409, 2.64459,3.38708,651 4.25198, 5.25084, 6.39647, 7.70319, 9.18708, 10.8663, 12.7617, 14.8968, 17.2990,20.0000};650 .024170, .100153, .234345, .433198, .703363, 1.05184, 1.48607, 2.01409, 2.64459, 3.38708, 651 4.25198, 5.25084, 6.39647, 7.70319, 9.18708, 10.8663, 12.7617, 14.8968, 17.2990, 20.0000}; 652 652 static const G4double IC[nX]={0, 653 .021877, .091263, .214670, .398677, .650133, .976322, 1.38510, 1.88504, 2.48555,3.19709,654 4.03129, 5.00127, 6.12184, 7.40989, 8.88482, 10.5690, 12.4888, 14.6748, 17.1638,20.0000};653 .021877, .091263, .214670, .398677, .650133, .976322, 1.38510, 1.88504, 2.48555, 3.19709, 654 4.03129, 5.00127, 6.12184, 7.40989, 8.88482, 10.5690, 12.4888, 14.6748, 17.1638, 20.0000}; 655 655 static const G4double ID[nX]={0, 656 .020062, .084127, .198702, .370384, .606100, .913288, 1.30006, 1.77535, 2.34912,3.03253,657 3.83822, 4.78063, 5.87634, 7.14459, 8.60791, 10.2929, 12.2315, 14.4621, 17.0320,20.0000};656 .020062, .084127, .198702, .370384, .606100, .913288, 1.30006, 1.77535, 2.34912, 3.03253, 657 3.83822, 4.78063, 5.87634, 7.14459, 8.60791, 10.2929, 12.2315, 14.4621, 17.0320, 20.0000}; 658 658 static const G4double IE[nX]={0, 659 .018547, .078104, .185102, .346090, .567998, .858331, 1.22535, 1.67824, 2.22735,2.88443,660 3.66294, 4.57845, 5.64911, 6.89637, 8.34578, 10.0282, 11.9812, 14.2519, 16.8993,20.0000};659 .018547, .078104, .185102, .346090, .567998, .858331, 1.22535, 1.67824, 2.22735, 2.88443, 660 3.66294, 4.57845, 5.64911, 6.89637, 8.34578, 10.0282, 11.9812, 14.2519, 16.8993, 20.0000}; 661 661 static const G4double IF[nX]={0, 662 .017143, .072466, .172271, .323007, .531545, .805393, 1.15288, 1.58338, 2.10754,2.73758,663 3.48769, 4.37450, 5.41770, 6.64092, 8.07288, 9.74894, 11.7135, 14.0232, 16.7522,20.0000};662 .017143, .072466, .172271, .323007, .531545, .805393, 1.15288, 1.58338, 2.10754, 2.73758, 663 3.48769, 4.37450, 5.41770, 6.64092, 8.07288, 9.74894, 11.7135, 14.0232, 16.7522, 20.0000}; 664 664 static const G4double IG[nX]={0, 665 .015618, .066285, .158094, .297316, .490692, .745653, 1.07053, 1.47479, 1.96931,2.56677,666 3.28205, 4.13289, 5.14068, 6.33158, 7.73808, 9.40133, 11.3745, 13.7279, 16.5577,20.0000};665 .015618, .066285, .158094, .297316, .490692, .745653, 1.07053, 1.47479, 1.96931, 2.56677, 666 3.28205, 4.13289, 5.14068, 6.33158, 7.73808, 9.40133, 11.3745, 13.7279, 16.5577, 20.0000}; 667 667 static const G4double IH[nX]={0, 668 .013702, .058434, .139923, .264115, .437466, .667179, .961433, 1.32965, 1.78283,2.33399,669 2.99871, 3.79596, 4.74916, 5.88771, 7.24937, 8.88367, 10.8576, 13.2646, 16.2417,20.0000};668 .013702, .058434, .139923, .264115, .437466, .667179, .961433, 1.32965, 1.78283, 2.33399, 669 2.99871, 3.79596, 4.74916, 5.88771, 7.24937, 8.88367, 10.8576, 13.2646, 16.2417, 20.0000}; 670 670 static const G4double II[nX]={0, 671 .011264, .048311, .116235, .220381, .366634, .561656, .813132, 1.13008, 1.52322,2.00554,672 2.59296, 3.30542, 4.16834, 5.21490, 6.48964, 8.05434, 9.99835, 12.4580, 15.6567,20.0000};671 .011264, .048311, .116235, .220381, .366634, .561656, .813132, 1.13008, 1.52322, 2.00554, 672 2.59296, 3.30542, 4.16834, 5.21490, 6.48964, 8.05434, 9.99835, 12.4580, 15.6567, 20.0000}; 673 673 static const G4double IJ[nX]={0, 674 .008628, .037206, .089928, .171242, .286114, .440251, .640343, .894382, 1.21208,1.60544,675 2.08962, 2.68414, 3.41486, 4.31700, 5.44048, 6.85936, 8.69067, 11.1358, 14.5885,20.0000};674 .008628, .037206, .089928, .171242, .286114, .440251, .640343, .894382, 1.21208, 1.60544, 675 2.08962, 2.68414, 3.41486, 4.31700, 5.44048, 6.85936, 8.69067, 11.1358, 14.5885, 20.0000}; 676 676 static const G4double* Il[nE]= 677 677 {I0,I1,I2,I3,I4,I5,I6,I7,I8,I9,IA,IB,IC,ID,IE,IF,IG,IH,II,IJ}; 678 678 static const G4double lE[nE]={ 679 -1.98842,-1.58049,-1.17256,-.764638,-.356711, .051215, .459141, .867068, 1.27499,1.68292,680 2.09085, 2.49877, 2.90670, 3.31463, 3.72255, 4.13048, 4.53840, 4.94633, 5.35426,5.76218};679 -1.98842,-1.58049,-1.17256,-.764638,-.356711, .051215, .459141, .867068, 1.27499, 1.68292, 680 2.09085, 2.49877, 2.90670, 3.31463, 3.72255, 4.13048, 4.53840, 4.94633, 5.35426, 5.76218}; 681 681 static const G4double lEmi=lE[0]; 682 682 static const G4double lEma=lE[nE-1]; 683 683 static const G4double dlE=(lEma-lEmi)/bE; 684 684 //*************************************************************************************** 685 685 G4double Enu=lastE; // Get energy of the last calculated cross-section 686 686 G4double lEn=std::log(Enu); // log(E) for interpolation … … 772 772 } 773 773 774 // This class can provide only virtual exchange pi + (a substitute for W+boson)774 // This class can provide only virtual exchange pi- (a substitute for W- boson) 775 775 G4int G4QANuENuclearCrossSection::GetExchangePDGCode() {return -211;}
Note: See TracChangeset
for help on using the changeset viewer.