CrossSections

Revision as of 04:49, 17 Dec 2006
213.42.21.77 (Talk | contribs)
http://morgan-john.realitsen.info
← Go to previous diff
Revision as of 12:39, 18 Dec 2006
217.20.214.8 (Talk | contribs)

Go to next diff →
Line 1: Line 1:
-[http://licata-riccardo.topo20.org licata, riccardo] - licata, riccardo licata, riccardo licata, riccardo [http://licata-riccardo.topo20.org licata, riccardo], [http://hell.fifa2007it.org hell] - hell hell hell [http://hell.fifa2007it.org hell], [http://ackermann-max.trevisos.org ackermann, max] - ackermann, max ackermann, max ackermann, max [http://ackermann-max.trevisos.org ackermann, max], [http://qixia-si.realitsen.info qixia si] - qixia si qixia si qixia si [http://qixia-si.realitsen.info qixia si], [http://alex-and.design4italy.org alex and emma] - alex and emma alex and emma alex and emma [http://alex-and.design4italy.org alex and emma], [http://capelle-jan.win-for-italy.info capelle, jan van de-] - capelle, jan van de- capelle, jan van de- capelle, jan van de- [http://capelle-jan.win-for-italy.info capelle, jan van de-], [http://bikom.trevisos.org bikom] - bikom bikom bikom [http://bikom.trevisos.org bikom], [http://paraylemonial.realitsen.info paray-le-monial] - paray-le-monial paray-le-monial paray-le-monial [http://paraylemonial.realitsen.info paray-le-monial], [http://mitsuyoshi-tosa.realitsen.info mitsuyoshi tosa] - mitsuyoshi tosa mitsuyoshi tosa mitsuyoshi tosa [http://mitsuyoshi-tosa.realitsen.info mitsuyoshi tosa], [http://four-brothers.fifa2007it.org four brothers - quattro fratelli] - four brothers - quattro fratelli four brothers - quattro fratelli four brothers - quattro fratelli [http://four-brothers.fifa2007it.org four brothers - quattro fratelli], [http://baumhauer-joseph.trevisos.org baumhauer, joseph] - baumhauer, joseph baumhauer, joseph baumhauer, joseph [http://baumhauer-joseph.trevisos.org baumhauer, joseph], [http://bentley-john.trevisos.org bentley, john francis] - bentley, john francis bentley, john francis bentley, john francis [http://bentley-john.trevisos.org bentley, john francis], [http://konarak.topo20.org konarak] - konarak konarak konarak [http://konarak.topo20.org konarak], [http://beauvais.trevisos.org beauvais] - beauvais beauvais beauvais [http://beauvais.trevisos.org beauvais], [http://il-viaggio.motnolado.org il viaggio di felicia] - il viaggio di felicia il viaggio di felicia il viaggio di felicia [http://il-viaggio.motnolado.org il viaggio di felicia], [http://fanum-voltumnae.win-for-italy.info fanum voltumnae] - fanum voltumnae fanum voltumnae fanum voltumnae [http://fanum-voltumnae.win-for-italy.info fanum voltumnae], [http://halle.topo20.org halle] - halle halle halle [http://halle.topo20.org halle], [http://divisionismo.win-for-italy.info divisionismo] - divisionismo divisionismo divisionismo [http://divisionismo.win-for-italy.info divisionismo], [http://vivin-louis.design4italy.org vivin, louis] - vivin, louis vivin, louis vivin, louis [http://vivin-louis.design4italy.org vivin, louis], [http://thule-springs.design4italy.org thule springs] - thule springs thule springs thule springs [http://thule-springs.design4italy.org thule springs], [http://abercrombie-sir.trevisos.org abercrombie, sir patrick] - abercrombie, sir patrick abercrombie, sir patrick abercrombie, sir patrick [http://abercrombie-sir.trevisos.org abercrombie, sir patrick], [http://tessenow-heinrich.design4italy.org tessenow heinrich] - tessenow heinrich tessenow heinrich tessenow heinrich [http://tessenow-heinrich.design4italy.org tessenow heinrich], [http://ebrach.win-for-italy.info ebrach] - ebrach ebrach ebrach [http://ebrach.win-for-italy.info ebrach], [http://piattabanda.realitsen.info piattabanda] - piattabanda piattabanda piattabanda [http://piattabanda.realitsen.info piattabanda], [http://piles-roger.realitsen.info piles, roger de-] - piles, roger de- piles, roger de- piles, roger de- [http://piles-roger.realitsen.info piles, roger de-], [http://charlotte-gray.fifa2007it.org charlotte gray] - charlotte gray charlotte gray charlotte gray [http://charlotte-gray.fifa2007it.org charlotte gray], [http://mills-robert.realitsen.info mills, robert] - mills, robert mills, robert mills, robert [http://mills-robert.realitsen.info mills, robert], [http://belfagor-.fifa2007it.org belfagor - il fantasma del louvre] - belfagor - il fantasma del louvre belfagor - il fantasma del louvre belfagor - il fantasma del louvre [http://belfagor-.fifa2007it.org belfagor - il fantasma del louvre], [http://cenerentola-a.fifa2007it.org cenerentola a parigi] - cenerentola a parigi cenerentola a parigi cenerentola a parigi [http://cenerentola-a.fifa2007it.org cenerentola a parigi], [http://la-spirale.motnolado.org la spirale della vendetta] - la spirale della vendetta la spirale della vendetta la spirale della vendetta [http://la-spirale.motnolado.org la spirale della vendetta], [http://praetexta.realitsen.info praetexta] - praetexta praetexta praetexta [http://praetexta.realitsen.info praetexta], [http://fashion.win-for-italy.info fashion] - fashion fashion fashion [http://fashion.win-for-italy.info fashion], [http://hilton-roger.topo20.org hilton, roger] - hilton, roger hilton, roger hilton, roger [http://hilton-roger.topo20.org hilton, roger], [http://lavedan-pierre.topo20.org lavedan, pierre] - lavedan, pierre lavedan, pierre lavedan, pierre [http://lavedan-pierre.topo20.org lavedan, pierre], [http://pontecagnano-faiano.realitsen.info pontecagnano faiano] - pontecagnano faiano pontecagnano faiano pontecagnano faiano [http://pontecagnano-faiano.realitsen.info pontecagnano faiano], [http://susanna.design4italy.org susanna] - susanna susanna susanna [http://susanna.design4italy.org susanna], [http://bambole-russe.fifa2007it.org bambole russe] - bambole russe bambole russe bambole russe [http://bambole-russe.fifa2007it.org bambole russe], [http://commedia-sexy.fifa2007it.org commedia sexy] - commedia sexy commedia sexy commedia sexy [http://commedia-sexy.fifa2007it.org commedia sexy], [http://bukezukuri.trevisos.org buke-zukuri] - buke-zukuri buke-zukuri buke-zukuri [http://bukezukuri.trevisos.org buke-zukuri], [http://costantino-i.win-for-italy.info costantino i (imperatore romano)] - costantino i (imperatore romano) costantino i (imperatore romano) costantino i (imperatore romano) [http://costantino-i.win-for-italy.info costantino i (imperatore romano)], [http://pawdawmu.realitsen.info pawdawmu] - pawdawmu pawdawmu pawdawmu [http://pawdawmu.realitsen.info pawdawmu], [http://gioco-a.fifa2007it.org gioco a due] - gioco a due gioco a due gioco a due [http://gioco-a.fifa2007it.org gioco a due], [http://leptis-magna.topo20.org leptis magna] - leptis magna leptis magna leptis magna [http://leptis-magna.topo20.org leptis magna], [http://killing-zoe.motnolado.org killing zoe] - killing zoe killing zoe killing zoe [http://killing-zoe.motnolado.org killing zoe], [http://puy-jean.realitsen.info puy, jean] - puy, jean puy, jean puy, jean [http://puy-jean.realitsen.info puy, jean], [http://jerry-maguire.motnolado.org jerry maguire] - jerry maguire jerry maguire jerry maguire [http://jerry-maguire.motnolado.org jerry maguire], [http://brouet-auguste.trevisos.org brouet, auguste] - brouet, auguste brouet, auguste brouet, auguste [http://brouet-auguste.trevisos.org brouet, auguste], [http://pileta-la.realitsen.info pileta, la-] - pileta, la- pileta, la- pileta, la- [http://pileta-la.realitsen.info pileta, la-], [http://campanile.win-for-italy.info campanile] - campanile campanile campanile [http://campanile.win-for-italy.info campanile], [http://arrieta-pedro.trevisos.org arrieta, pedro de-] - arrieta, pedro de- arrieta, pedro de- arrieta, pedro de- [http://arrieta-pedro.trevisos.org arrieta, pedro de-], [http://cambridge-usa.win-for-italy.info cambridge (u.s.a.)] - cambridge (u.s.a.) cambridge (u.s.a.) cambridge (u.s.a.) [http://cambridge-usa.win-for-italy.info cambridge (u.s.a.)], [http://arellano-juan.trevisos.org arellano, juan de-] - arellano, juan de- arellano, juan de- arellano, juan de- [http://arellano-juan.trevisos.org arellano, juan de-], [http://la-cameriera.motnolado.org la cameriera] - la cameriera la cameriera la cameriera [http://la-cameriera.motnolado.org la cameriera], [http://honegger-david.topo20.org honegger, david] - honegger, david honegger, david honegger, david [http://honegger-david.topo20.org honegger, david], [http://hals-frans.topo20.org hals, frans] - hals, frans hals, frans hals, frans [http://hals-frans.topo20.org hals, frans], [http://black-simphony.fifa2007it.org black simphony] - black simphony black simphony black simphony [http://black-simphony.fifa2007it.org black simphony], [http://last-night.motnolado.org last night] - last night last night last night [http://last-night.motnolado.org last night], [http://spencer-stanley.design4italy.org spencer, stanley] - spencer, stanley spencer, stanley spencer, stanley [http://spencer-stanley.design4italy.org spencer, stanley], [http://gertrud.fifa2007it.org gertrud] - gertrud gertrud gertrud [http://gertrud.fifa2007it.org gertrud], [http://tachi.design4italy.org tachi] - tachi tachi tachi [http://tachi.design4italy.org tachi], [http://tisi-benvenuto.design4italy.org tisi, benvenuto] - tisi, benvenuto tisi, benvenuto tisi, benvenuto [http://tisi-benvenuto.design4italy.org tisi, benvenuto], [http://novilara.realitsen.info novilara] - novilara novilara novilara [http://novilara.realitsen.info novilara], [http://dragonfly.fifa2007it.org dragonfly] - dragonfly dragonfly dragonfly [http://dragonfly.fifa2007it.org dragonfly], [http://ricostruire.design4italy.org ricostruire] - ricostruire ricostruire ricostruire [http://ricostruire.design4italy.org ricostruire], [http://base-architettura.trevisos.org base (architettura)] - base (architettura) base (architettura) base (architettura) [http://base-architettura.trevisos.org base (architettura)], [http://muybridge-eadweard.realitsen.info muybridge, eadweard] - muybridge, eadweard muybridge, eadweard muybridge, eadweard [http://muybridge-eadweard.realitsen.info muybridge, eadweard], [http://kahn-louis.topo20.org kahn, louis isadore] - kahn, louis isadore kahn, louis isadore kahn, louis isadore [http://kahn-louis.topo20.org kahn, louis isadore], [http://looy-jacobus.topo20.org looy, jacobus van-] - looy, jacobus van- looy, jacobus van- looy, jacobus van- [http://looy-jacobus.topo20.org looy, jacobus van-], [http://confessions-a.fifa2007it.org confessions of a teenage drama queen] - confessions of a teenage drama queen confessions of a teenage drama queen confessions of a teenage drama queen [http://confessions-a.fifa2007it.org confessions of a teenage drama queen], [http://liberty.topo20.org liberty] - liberty liberty liberty [http://liberty.topo20.org liberty], [http://brinckmann-albert.trevisos.org brinckmann, albert erich] - brinckmann, albert erich brinckmann, albert erich brinckmann, albert erich [http://brinckmann-albert.trevisos.org brinckmann, albert erich], [http://carne-tremula.fifa2007it.org carne tremula] - carne tremula carne tremula carne tremula [http://carne-tremula.fifa2007it.org carne tremula], [http://a-good.fifa2007it.org a good lawyer's wife] - a good lawyer's wife a good lawyer's wife a good lawyer's wife [http://a-good.fifa2007it.org a good lawyer's wife], [http://bourgeois-victor.trevisos.org bourgeois, victor] - bourgeois, victor bourgeois, victor bourgeois, victor [http://bourgeois-victor.trevisos.org bourgeois, victor], [http://gods-and.fifa2007it.org gods and generals] - gods and generals gods and generals gods and generals [http://gods-and.fifa2007it.org gods and generals], [http://fragile.fifa2007it.org fragile] - fragile fragile fragile [http://fragile.fifa2007it.org fragile], [http://kitaj-ronald.topo20.org kitaj, ronald b.] - kitaj, ronald b. kitaj, ronald b. kitaj, ronald b. [http://kitaj-ronald.topo20.org kitaj, ronald b.], [http://foster-sir.win-for-italy.info foster, sir norman] - foster, sir norman foster, sir norman foster, sir norman [http://foster-sir.win-for-italy.info foster, sir norman], [http://vos-cornelis.design4italy.org vos, cornelis de-] - vos, cornelis de- vos, cornelis de- vos, cornelis de- [http://vos-cornelis.design4italy.org vos, cornelis de-], [http://swatch.design4italy.org swatch] - swatch swatch swatch [http://swatch.design4italy.org swatch], [http://venet-bernard.design4italy.org venet, bernard] - venet, bernard venet, bernard venet, bernard [http://venet-bernard.design4italy.org venet, bernard], [http://top.design4italy.org top] - top top top [http://top.design4italy.org top], [http://mora.realitsen.info mora] - mora mora mora [http://mora.realitsen.info mora], [http://big-fish.fifa2007it.org big fish] - big fish big fish big fish [http://big-fish.fifa2007it.org big fish], [http://bosch-y.trevisos.org bosch y gimpera, pedro] - bosch y gimpera, pedro bosch y gimpera, pedro bosch y gimpera, pedro [http://bosch-y.trevisos.org bosch y gimpera, pedro], [http://pailou.realitsen.info pailou] - pailou pailou pailou [http://pailou.realitsen.info pailou], [http://nylon.realitsen.info nylon] - nylon nylon nylon [http://nylon.realitsen.info nylon], [http://carlitos-way.fifa2007it.org carlito's way] - carlito's way carlito's way carlito's way [http://carlitos-way.fifa2007it.org carlito's way], [http://loggiato.topo20.org loggiato] - loggiato loggiato loggiato [http://loggiato.topo20.org loggiato], [http://avedon-richard.trevisos.org avedon, richard] - avedon, richard avedon, richard avedon, richard [http://avedon-richard.trevisos.org avedon, richard], [http://orchard-william.realitsen.info orchard, william] - orchard, william orchard, william orchard, william [http://orchard-william.realitsen.info orchard, william], [http://gremlins.fifa2007it.org gremlins] - gremlins gremlins gremlins [http://gremlins.fifa2007it.org gremlins], [http://accademismo.trevisos.org accademismo] - accademismo accademismo accademismo [http://accademismo.trevisos.org accademismo], [http://candela-felix.win-for-italy.info candela, felix] - candela, felix candela, felix candela, felix [http://candela-felix.win-for-italy.info candela, felix], [http://azzolini-tito.trevisos.org azzolini, tito] - azzolini, tito azzolini, tito azzolini, tito [http://azzolini-tito.trevisos.org azzolini, tito], [http://kiva.topo20.org kiva] - kiva kiva kiva [http://kiva.topo20.org kiva], [http://robert-luzarches.design4italy.org robert de luzarches] - robert de luzarches robert de luzarches robert de luzarches [http://robert-luzarches.design4italy.org robert de luzarches], [http://del-perduto.fifa2007it.org del perduto amore] - del perduto amore del perduto amore del perduto amore [http://del-perduto.fifa2007it.org del perduto amore], [http://shumman-kubo.design4italy.org shumman, kubo] - shumman, kubo shumman, kubo shumman, kubo [http://shumman-kubo.design4italy.org shumman, kubo], [http://morgan-john.realitsen.info morgan, john pierpont] - morgan, john pierpont morgan, john pierpont morgan, john pierpont [http://morgan-john.realitsen.info morgan, john pierpont], morgan, john pierpont+==Как построить графики сечений, используемых Geant'ом==
 + 
 +Этот вопрос, судя по форумам g4 на сайте ЦЕРН, будоражит многие умы. К сожалению, универсального метода извлечения сечений для конкретных процессов по всей видимости не существует (можно, конечно, смоделировать маленький опыт Резерфорда, но это не очень идеальный вариант). Однако, если вас интересуют сечения стандартных электромагнитных взаимодействий, то есть обнадеживающее обстоятельство. Классы этих процессов содержат статический метод <code>ComputeCrossSectionPerAtom</code>, который и даст необходимые сечения.
 + 
 +Вот код программы, которая создает 100 файлов, в которых содержатся сечения в барнах для энергий 0&mdash;10&nbsp;MeV.
 + #include "G4Material.hh"
 + #include "G4PhotoElectricEffect.hh"
 + #include "G4ComptonScattering.hh"
 + #include "G4GammaConversion.hh"
 +
 + #include "globals.hh"
 + #include "G4UnitsTable.hh"
 + #include <string.h>
 + #include <fstream>
 + #include "main.hh"
 + using namespace std;
 +
 + struct material {
 + int z;
 + G4Material *g4mat;
 + };
 + const double dm3=1000*cm3;
 + int doall() {
 +
 + G4UnitDefinition::BuildUnitsTable();
 +
 + material mats[100];
 + mats[0].z=1;mats[0].g4mat=new G4Material
 + ("Hydrogen", 1, 1.00794*g/mole, 0.08988*g/dm3);
 + mats[1].z=2; mats[1].g4mat=new G4Material
 + ("Helium", 2, 4.002602*g/mole, 0.1785*g/dm3);
 + mats[2].z=3; mats[2].g4mat=new G4Material
 + ("Lithium", 3, 6.941*g/mole, 0.534*g/cm3);
 + mats[3].z=4; mats[3].g4mat=new G4Material
 + ("Beryllium", 4, 9.012182*g/mole, 1.848*g/cm3);
 + mats[4].z=5; mats[4].g4mat=new G4Material
 + ("Boron", 5, 10.811*g/mole, 2.37*g/cm3);
 + mats[5].z=6; mats[5].g4mat=new G4Material
 + ("Carbon", 6, 12.01115*g/mole, 1.8*g/cm3);
 + mats[6].z=7; mats[6].g4mat=new G4Material
 + ("Nitrogen", 7, 14.00674*g/mole, 0.8081*g/cm3);
 + mats[7].z=8; mats[7].g4mat=new G4Material
 + ("Oxygen", 8, 15.9994*g/mole, 1.429*g/dm3);
 + mats[8].z=9; mats[8].g4mat=new G4Material
 + ("Flourine", 9, 18.9984032*g/mole, 1.69*g/dm3);
 + mats[9].z=10; mats[9].g4mat=new G4Material
 + ("Neon", 10, 20.1797*g/mole, 0.899*g/dm3);
 + mats[10].z=11; mats[10].g4mat=new G4Material
 + ("Sodium", 11, 22.989768*g/mole, 0.971*g/cm3);
 + mats[11].z=12; mats[11].g4mat=new G4Material
 + ("Magnesium", 12, 24.3050*g/mole, 1.783*g/cm3);
 + mats[12].z=13; mats[12].g4mat=new G4Material
 + ("Aluminium", 13, 26.981539*g/mole, 2.6989*g/cm3);
 + mats[13].z=14; mats[13].g4mat=new G4Material
 + ("Silicon", 14, 28.0855*g/mole, 2.33*g/cm3);
 + mats[14].z=15; mats[14].g4mat=new G4Material
 + ("Phosphorus", 15, 30.973762*g/mole, 1.82*g/cm3);
 + mats[15].z=16; mats[15].g4mat=new G4Material
 + ("Sulphur", 16, 32.066*g/mole, 1.957*g/cm3);
 + mats[16].z=17; mats[16].g4mat=new G4Material
 + ("Chlorine", 17, 35.4527*g/mole, 3.214*g/dm3);
 + mats[17].z=18; mats[17].g4mat=new G4Material
 + ("Argon", 18, 39.948*g/mole, 1.783*g/dm3);
 + mats[18].z=19; mats[18].g4mat=new G4Material
 + ("Potassium", 19, 39.0983*g/mole, 0.862*g/cm3);
 + mats[19].z=20; mats[19].g4mat=new G4Material
 + ("Calcium", 20, 40.078*g/mole, 1.55*g/cm3);
 + mats[20].z=21; mats[20].g4mat=new G4Material
 + ("Scandium", 21, 44.955910*g/mole, 2.989*g/cm3);
 + mats[21].z=22; mats[21].g4mat=new G4Material
 + ("Titanium", 22, 47.88*g/mole, 4.54*g/cm3);
 + mats[22].z=23; mats[22].g4mat=new G4Material
 + ("Vanadium", 23, 50.9415*g/mole, 6.11*g/cm3);
 + mats[23].z=24; mats[23].g4mat=new G4Material
 + ("Chromium", 24, 51.9961*g/mole, 7.14*g/cm3);
 + mats[24].z=25; mats[24].g4mat=new G4Material
 + ("Manganese", 25, 54.93805*g/mole, 7.44*g/cm3);
 + mats[25].z=26; mats[25].g4mat=new G4Material
 + ("Iron", 26, 55.847*g/mole, 7.874*g/cm3);
 + mats[26].z=27; mats[26].g4mat=new G4Material
 + ("Cobalt", 27, 58.93320*g/mole, 8.9*g/cm3);
 + mats[27].z=28; mats[27].g4mat=new G4Material
 + ("Nickel", 28, 58.6934*g/mole, 8.902*g/cm3);
 + mats[28].z=29; mats[28].g4mat=new G4Material
 + ("Copper", 29, 63.546*g/mole, 8.96*g/cm3);
 + mats[29].z=30; mats[29].g4mat=new G4Material
 + ("Zinc", 30, 65.39*g/mole, 7.133*g/cm3);
 + mats[30].z=31; mats[30].g4mat=new G4Material
 + ("Gallium", 31, 69.723*g/mole, 6.095*g/cm3);
 + mats[31].z=32; mats[31].g4mat=new G4Material
 + ("Germanium", 32, 72.61*g/mole, 5.323*g/cm3);
 + mats[32].z=33; mats[32].g4mat=new G4Material
 + ("Arsenic", 33, 74.92159*g/mole, 5.727*g/cm3);
 + mats[33].z=34; mats[33].g4mat=new G4Material
 + ("Selenium", 34, 78.96*g/mole, 4.28*g/cm3);
 + mats[34].z=35; mats[34].g4mat=new G4Material
 + ("Bromine", 35, 79.904*g/mole, 3.119*g/cm3);
 + mats[35].z=36; mats[35].g4mat=new G4Material
 + ("Krypton", 36, 83.80*g/mole, 3.73*g/dm3);
 + mats[36].z=37; mats[36].g4mat=new G4Material
 + ("Rubidium", 37, 85.4678*g/mole, 1.532*g/cm3);
 + mats[37].z=38; mats[37].g4mat=new G4Material
 + ("Strontium", 38, 87.62*g/mole, 2.54*g/cm3);
 + mats[38].z=39; mats[38].g4mat=new G4Material
 + ("Yttrium", 39, 88.90585*g/mole, 4.469*g/cm3);
 + mats[39].z=40; mats[39].g4mat=new G4Material
 + ("Zirconium", 40, 91.224*g/mole, 6.506*g/cm3);
 + mats[40].z=41; mats[40].g4mat=new G4Material
 + ("Niobium", 41, 92.90638*g/mole, 8.57*g/cm3);
 + mats[41].z=42; mats[41].g4mat=new G4Material
 + ("Molybdenum", 42, 95.94*g/mole, 10.22*g/cm3);
 + mats[42].z=43; mats[42].g4mat=new G4Material
 + ("Technetium", 43,
 + (98.9063)*g/mole,
 + (11.50)*g/cm3);
 + mats[43].z=44; mats[43].g4mat=new G4Material
 + ("Ruthenium", 44, 101.07*g/mole, 12.41*g/cm3);
 + mats[44].z=45; mats[44].g4mat=new G4Material
 + ("Rhodium", 45, 102.90550*g/mole, 12.41*g/cm3);
 + mats[45].z=46; mats[45].g4mat=new G4Material
 + ("Palladium", 46, 106.42*g/mole, 12.02*g/cm3);
 + mats[46].z=47; mats[46].g4mat=new G4Material
 + ("Silver", 47, 107.8682*g/mole, 10.49*g/cm3);
 + mats[47].z=48; mats[47].g4mat=new G4Material
 + ("Cadmium", 48, 112.411*g/mole, 8.642*g/cm3);
 + mats[48].z=49; mats[48].g4mat=new G4Material
 + ("Indium", 49, 114.818*g/mole, 7.31*g/cm3);
 + mats[49].z=50; mats[49].g4mat=new G4Material
 + ("Tin", 50, 118.710*g/mole, 7.31*g/cm3);
 + mats[50].z=51; mats[50].g4mat=new G4Material
 + ("Antimony", 51, 121.757*g/mole, 6.697*g/cm3);
 + mats[51].z=52; mats[51].g4mat=new G4Material
 + ("Tellurium", 52, 127.60*g/mole, 6.24*g/cm3);
 + mats[52].z=53; mats[52].g4mat=new G4Material
 + ("Iodine", 53, 126.90447*g/mole, 4.93*g/cm3);
 + mats[53].z=54; mats[53].g4mat=new G4Material
 + ("Xenon", 54, 131.29*g/mole, 5.887*g/dm3);
 + mats[54].z=55; mats[54].g4mat=new G4Material
 + ("Caesium", 55, 132.90543*g/mole, 1.873*g/cm3);
 + mats[55].z=56; mats[55].g4mat=new G4Material
 + ("Barium", 56, 137.327*g/mole, 3.51*g/cm3);
 + mats[56].z=57; mats[56].g4mat=new G4Material
 + ("Lanthanum", 57, 138.9055*g/mole, 6.145*g/cm3);
 + mats[57].z=72; mats[57].g4mat=new G4Material
 + ("Hafnium", 72, 178.49*g/mole, 13.28*g/cm3);
 + mats[58].z=73; mats[58].g4mat=new G4Material
 + ("Tantalum", 73, 180.9479*g/mole, 16.654*g/cm3);
 + mats[59].z=74; mats[59].g4mat=new G4Material
 + ("Tungsten", 74, 183.84*g/mole, 19.35*g/cm3);
 + mats[60].z=75; mats[60].g4mat=new G4Material
 + ("Rhenium", 75, 186.207*g/mole, 21.02*g/cm3);
 + mats[61].z=76; mats[61].g4mat=new G4Material
 + ("Osmium", 76, 190.23*g/mole, 22.57*g/cm3);
 + mats[62].z=77; mats[62].g4mat=new G4Material
 + ("Iridium", 77, 192.22*g/mole, 22.421*g/cm3);
 + mats[63].z=78; mats[63].g4mat=new G4Material
 + ("Platinum", 78, 195.08*g/mole, 21.45*g/cm3);
 + mats[64].z=79; mats[64].g4mat=new G4Material
 + ("Gold", 79, 196.96654*g/mole, 19.3*g/cm3);
 + mats[65].z=80; mats[65].g4mat=new G4Material
 + ("Mercury", 80, 200.59*g/mole, 13.546*g/cm3);
 + mats[66].z=81; mats[66].g4mat=new G4Material
 + ("Thallium", 81, 204.3833*g/mole, 11.85*g/cm3);
 + mats[67].z=82; mats[67].g4mat=new G4Material
 + ("Lead", 82, 207.2*g/mole, 11.3437*g/cm3);
 + mats[68].z=83; mats[68].g4mat=new G4Material
 + ("Bismuth", 83, 208.98037*g/mole, 9.747*g/cm3);
 + mats[69].z=84; mats[69].g4mat=new G4Material
 + ("Polonium", 84,
 + (209.9828)*g/mole, 9.32*g/cm3);
 + mats[70].z=86; mats[70].g4mat=new G4Material
 + ("Radon", 86,
 + (222.0176)*g/mole, 9.73*g/dm3);
 + mats[71].z=87; mats[71].g4mat=new G4Material
 + ("Francium", 87,
 + (223.0197)*g/mole, 1.87*g/cm3);
 + mats[72].z=88; mats[72].g4mat=new G4Material
 + ("Radium", 88,
 + (226.0254)*g/mole, 5.5*g/cm3);
 + mats[73].z=89; mats[73].g4mat=new G4Material
 + ("Actinium", 89,
 + (227.0278)*g/mole, 10.07*g/cm3);
 + mats[74].z=58; mats[74].g4mat=new G4Material
 + ("Cerium", 58, 140.115*g/mole, 6.773*g/cm3);
 + mats[75].z=59; mats[75].g4mat=new G4Material
 + ("Praseodymium", 59, 140.90765*g/mole, 6.773*g/cm3);
 + mats[76].z=60; mats[76].g4mat=new G4Material
 + ("Neodymium", 60, 144.24*g/mole, 7.008*g/cm3);
 + mats[77].z=61; mats[77].g4mat=new G4Material
 + ("Promethium", 61,
 + (144.9127)*g/mole, 7.26*g/cm3);
 + mats[78].z=62; mats[78].g4mat=new G4Material
 + ("Samarium", 62, 150.36*g/mole, 7.520*g/cm3);
 + mats[79].z=63; mats[79].g4mat=new G4Material
 + ("Europium", 63, 151.965*g/mole, 5.244*g/cm3);
 + mats[80].z=64; mats[80].g4mat=new G4Material
 + ("Gadolinium", 64, 157.25*g/mole, 7.901*g/cm3);
 + mats[81].z=65; mats[81].g4mat=new G4Material
 + ("Terbium", 65, 158.92534*g/mole, 8.230*g/cm3);
 + mats[82].z=66; mats[82].g4mat=new G4Material
 + ("Dysprosium", 66, 162.50*g/mole, 8.551*g/cm3);
 + mats[83].z=67; mats[83].g4mat=new G4Material
 + ("Holmium", 67, 164.93032*g/mole, 8.795*g/cm3);
 + mats[84].z=68; mats[84].g4mat=new G4Material
 + ("Erbium", 68, 167.26*g/mole, 9.066*g/cm3);
 + mats[85].z=69; mats[85].g4mat=new G4Material
 + ("Thulium", 69, 168.93421*g/mole, 9.321*g/cm3);
 + mats[86].z=70; mats[86].g4mat=new G4Material
 + ("Ytterbium", 70, 173.04*g/mole, 6.966*g/cm3);
 + mats[87].z=71; mats[87].g4mat=new G4Material
 + ("Lutetium", 71, 174.967*g/mole, 9.841*g/cm3);
 + mats[88].z=90; mats[88].g4mat=new G4Material
 + ("Thorium", 90, 232.0381*g/mole, 11.72*g/cm3);
 + mats[89].z=91; mats[89].g4mat=new G4Material
 + ("Protactinium", 91,
 + (231.03588)*g/mole, 15.37*g/cm3);
 + mats[90].z=92; mats[90].g4mat=new G4Material
 + ("Uranium", 92, 238.0508*g/mole, 18.05*g/cm3);
 + mats[91].z=93; mats[91].g4mat=new G4Material
 + ("Neptunium", 93,
 + (237.0482)*g/mole, 20.25*g/cm3);
 + mats[92].z=94; mats[92].g4mat=new G4Material
 + ("Plutonium", 94,
 + (244.0642)*g/mole, 19.84*g/cm3);
 + mats[93].z=95; mats[93].g4mat=new G4Material
 + ("Americium", 95,
 + (243.0614)*g/mole, 13.67*g/cm3);
 + mats[94].z=96; mats[94].g4mat=new G4Material
 + ("Curium", 96,
 + (247.0703)*g/mole, 13.51*g/cm3);
 + mats[95].z=97; mats[95].g4mat=new G4Material
 + ("Berkelium", 97, (247.0703)*g/mole, 14*g/cm3);
 +
 + // initialise processes
 + //
 + G4PhotoElectricEffect* phot = new G4PhotoElectricEffect();
 + G4ComptonScattering* comp = new G4ComptonScattering();
 + G4GammaConversion* conv = new G4GammaConversion();
 +
 + // print cross section per atom
 + //
 + G4double Emin = 5*keV, Emax = 10*MeV, dE = 10*keV;
 + material *iter = mats;
 + double Z;
 + for(int i=0; i<96; i++) {
 + char buf[100];
 + sprintf(buf, "%02d.txt", mats[i].z);
 + Z = mats[i].z;
 + ofstream file(buf);
 + for (G4double Energy = Emin; Energy <= Emax; Energy += dE) {
 + file << Energy/MeV << " " <<
 + comp->ComputeCrossSectionPerAtom(Energy,Z)/barn +
 + phot->ComputeCrossSectionPerAtom(Energy,Z)/barn +
 + conv->ComputeCrossSectionPerAtom(Energy,Z)/barn
 + << endl;
 + }
 + }
 +
 + return EXIT_SUCCESS;
 + }
 + 
 +Для работы потребовалось определить все элементы, задействованные в вычислениях, т. к. это необходимо для работы <code>G4PhotoElectricProcess</code>.

Revision as of 12:39, 18 Dec 2006

Как построить графики сечений, используемых Geant'ом

Этот вопрос, судя по форумам g4 на сайте ЦЕРН, будоражит многие умы. К сожалению, универсального метода извлечения сечений для конкретных процессов по всей видимости не существует (можно, конечно, смоделировать маленький опыт Резерфорда, но это не очень идеальный вариант). Однако, если вас интересуют сечения стандартных электромагнитных взаимодействий, то есть обнадеживающее обстоятельство. Классы этих процессов содержат статический метод ComputeCrossSectionPerAtom, который и даст необходимые сечения.

Вот код программы, которая создает 100 файлов, в которых содержатся сечения в барнах для энергий 0—10 MeV.

#include "G4Material.hh"
#include "G4PhotoElectricEffect.hh"
#include "G4ComptonScattering.hh"
#include "G4GammaConversion.hh"

#include "globals.hh"
#include "G4UnitsTable.hh"
#include <string.h>
#include <fstream>
#include "main.hh"
using namespace std;

struct material {
	int z;
	G4Material *g4mat;
};
const double dm3=1000*cm3;
int doall() {
 
	G4UnitDefinition::BuildUnitsTable(); 

	material mats[100];
	mats[0].z=1;mats[0].g4mat=new G4Material
			    ("Hydrogen", 1, 1.00794*g/mole, 0.08988*g/dm3);
	mats[1].z=2; mats[1].g4mat=new G4Material
			     ("Helium", 2, 4.002602*g/mole, 0.1785*g/dm3);
	mats[2].z=3; mats[2].g4mat=new G4Material
			     ("Lithium", 3, 6.941*g/mole, 0.534*g/cm3);
	mats[3].z=4; mats[3].g4mat=new G4Material
			     ("Beryllium", 4, 9.012182*g/mole, 1.848*g/cm3);
	mats[4].z=5; mats[4].g4mat=new G4Material
			     ("Boron", 5, 10.811*g/mole, 2.37*g/cm3);
	mats[5].z=6; mats[5].g4mat=new G4Material
			     ("Carbon", 6, 12.01115*g/mole, 1.8*g/cm3);
	mats[6].z=7; mats[6].g4mat=new G4Material
			     ("Nitrogen", 7, 14.00674*g/mole, 0.8081*g/cm3);
	mats[7].z=8; mats[7].g4mat=new G4Material
			     ("Oxygen", 8, 15.9994*g/mole, 1.429*g/dm3);
	mats[8].z=9; mats[8].g4mat=new G4Material
			     ("Flourine", 9, 18.9984032*g/mole, 1.69*g/dm3);
	mats[9].z=10; mats[9].g4mat=new G4Material
			      ("Neon", 10, 20.1797*g/mole, 0.899*g/dm3);
	mats[10].z=11; mats[10].g4mat=new G4Material
			       ("Sodium", 11, 22.989768*g/mole, 0.971*g/cm3);
	mats[11].z=12; mats[11].g4mat=new G4Material
			       ("Magnesium", 12, 24.3050*g/mole, 1.783*g/cm3);
	mats[12].z=13; mats[12].g4mat=new G4Material
			       ("Aluminium", 13, 26.981539*g/mole, 2.6989*g/cm3);
	mats[13].z=14; mats[13].g4mat=new G4Material
			       ("Silicon", 14, 28.0855*g/mole, 2.33*g/cm3);
	mats[14].z=15; mats[14].g4mat=new G4Material
			       ("Phosphorus", 15, 30.973762*g/mole, 1.82*g/cm3);
	mats[15].z=16; mats[15].g4mat=new G4Material
			       ("Sulphur", 16, 32.066*g/mole, 1.957*g/cm3);
	mats[16].z=17; mats[16].g4mat=new G4Material
			       ("Chlorine", 17, 35.4527*g/mole, 3.214*g/dm3);
	mats[17].z=18; mats[17].g4mat=new G4Material
			       ("Argon", 18, 39.948*g/mole, 1.783*g/dm3);
	mats[18].z=19; mats[18].g4mat=new G4Material
			       ("Potassium", 19, 39.0983*g/mole, 0.862*g/cm3);
	mats[19].z=20; mats[19].g4mat=new G4Material
			       ("Calcium", 20, 40.078*g/mole, 1.55*g/cm3);
	mats[20].z=21; mats[20].g4mat=new G4Material
			       ("Scandium", 21, 44.955910*g/mole, 2.989*g/cm3);
	mats[21].z=22; mats[21].g4mat=new G4Material
			       ("Titanium", 22, 47.88*g/mole, 4.54*g/cm3);
	mats[22].z=23; mats[22].g4mat=new G4Material
			       ("Vanadium", 23, 50.9415*g/mole, 6.11*g/cm3);
	mats[23].z=24; mats[23].g4mat=new G4Material
			       ("Chromium", 24, 51.9961*g/mole, 7.14*g/cm3);
	mats[24].z=25; mats[24].g4mat=new G4Material
			       ("Manganese", 25, 54.93805*g/mole, 7.44*g/cm3);
	mats[25].z=26; mats[25].g4mat=new G4Material
			       ("Iron", 26, 55.847*g/mole, 7.874*g/cm3);
	mats[26].z=27; mats[26].g4mat=new G4Material
			       ("Cobalt", 27, 58.93320*g/mole, 8.9*g/cm3);
	mats[27].z=28; mats[27].g4mat=new G4Material
			       ("Nickel", 28, 58.6934*g/mole, 8.902*g/cm3);
	mats[28].z=29; mats[28].g4mat=new G4Material
			       ("Copper", 29, 63.546*g/mole, 8.96*g/cm3);
	mats[29].z=30; mats[29].g4mat=new G4Material
			       ("Zinc", 30, 65.39*g/mole, 7.133*g/cm3);
	mats[30].z=31; mats[30].g4mat=new G4Material
			       ("Gallium", 31, 69.723*g/mole, 6.095*g/cm3);
	mats[31].z=32; mats[31].g4mat=new G4Material
			       ("Germanium", 32, 72.61*g/mole, 5.323*g/cm3);
	mats[32].z=33; mats[32].g4mat=new G4Material
			       ("Arsenic", 33, 74.92159*g/mole, 5.727*g/cm3);
	mats[33].z=34; mats[33].g4mat=new G4Material
			       ("Selenium", 34, 78.96*g/mole, 4.28*g/cm3);
	mats[34].z=35; mats[34].g4mat=new G4Material
			       ("Bromine", 35, 79.904*g/mole, 3.119*g/cm3);
	mats[35].z=36; mats[35].g4mat=new G4Material
			       ("Krypton", 36, 83.80*g/mole, 3.73*g/dm3);
	mats[36].z=37; mats[36].g4mat=new G4Material
			       ("Rubidium", 37, 85.4678*g/mole, 1.532*g/cm3);
	mats[37].z=38; mats[37].g4mat=new G4Material
			       ("Strontium", 38, 87.62*g/mole, 2.54*g/cm3);
	mats[38].z=39; mats[38].g4mat=new G4Material
			       ("Yttrium", 39, 88.90585*g/mole, 4.469*g/cm3);
	mats[39].z=40; mats[39].g4mat=new G4Material
			       ("Zirconium", 40, 91.224*g/mole, 6.506*g/cm3);
	mats[40].z=41; mats[40].g4mat=new G4Material
			       ("Niobium", 41, 92.90638*g/mole, 8.57*g/cm3);
	mats[41].z=42; mats[41].g4mat=new G4Material
			       ("Molybdenum", 42, 95.94*g/mole, 10.22*g/cm3);
	mats[42].z=43; mats[42].g4mat=new G4Material
			       ("Technetium", 43,
				(98.9063)*g/mole,
				(11.50)*g/cm3);
	mats[43].z=44; mats[43].g4mat=new G4Material
			       ("Ruthenium", 44, 101.07*g/mole, 12.41*g/cm3);
	mats[44].z=45; mats[44].g4mat=new G4Material
			       ("Rhodium", 45, 102.90550*g/mole, 12.41*g/cm3);
	mats[45].z=46; mats[45].g4mat=new G4Material
			       ("Palladium", 46, 106.42*g/mole, 12.02*g/cm3);
	mats[46].z=47; mats[46].g4mat=new G4Material
			       ("Silver", 47, 107.8682*g/mole, 10.49*g/cm3);
	mats[47].z=48; mats[47].g4mat=new G4Material
			       ("Cadmium", 48, 112.411*g/mole, 8.642*g/cm3);
	mats[48].z=49; mats[48].g4mat=new G4Material
			       ("Indium", 49, 114.818*g/mole, 7.31*g/cm3);
	mats[49].z=50; mats[49].g4mat=new G4Material
			       ("Tin", 50, 118.710*g/mole, 7.31*g/cm3);
	mats[50].z=51; mats[50].g4mat=new G4Material
			       ("Antimony", 51, 121.757*g/mole, 6.697*g/cm3);
	mats[51].z=52; mats[51].g4mat=new G4Material
			       ("Tellurium", 52, 127.60*g/mole, 6.24*g/cm3);
	mats[52].z=53; mats[52].g4mat=new G4Material
			       ("Iodine", 53, 126.90447*g/mole, 4.93*g/cm3);
	mats[53].z=54; mats[53].g4mat=new G4Material
			       ("Xenon", 54, 131.29*g/mole, 5.887*g/dm3);
	mats[54].z=55; mats[54].g4mat=new G4Material
			       ("Caesium", 55, 132.90543*g/mole, 1.873*g/cm3);
	mats[55].z=56; mats[55].g4mat=new G4Material
			       ("Barium", 56, 137.327*g/mole, 3.51*g/cm3);
	mats[56].z=57; mats[56].g4mat=new G4Material
			       ("Lanthanum", 57, 138.9055*g/mole, 6.145*g/cm3);
	mats[57].z=72; mats[57].g4mat=new G4Material
			       ("Hafnium", 72, 178.49*g/mole, 13.28*g/cm3);
	mats[58].z=73; mats[58].g4mat=new G4Material
			       ("Tantalum", 73, 180.9479*g/mole, 16.654*g/cm3);
	mats[59].z=74; mats[59].g4mat=new G4Material
			       ("Tungsten", 74, 183.84*g/mole, 19.35*g/cm3);
	mats[60].z=75; mats[60].g4mat=new G4Material
			       ("Rhenium", 75, 186.207*g/mole, 21.02*g/cm3);
	mats[61].z=76; mats[61].g4mat=new G4Material
			       ("Osmium", 76, 190.23*g/mole, 22.57*g/cm3);
	mats[62].z=77; mats[62].g4mat=new G4Material
			       ("Iridium", 77, 192.22*g/mole, 22.421*g/cm3);
	mats[63].z=78; mats[63].g4mat=new G4Material
			       ("Platinum", 78, 195.08*g/mole, 21.45*g/cm3);
	mats[64].z=79; mats[64].g4mat=new G4Material
			       ("Gold", 79, 196.96654*g/mole, 19.3*g/cm3);
	mats[65].z=80; mats[65].g4mat=new G4Material
			       ("Mercury", 80, 200.59*g/mole, 13.546*g/cm3);
	mats[66].z=81; mats[66].g4mat=new G4Material
			       ("Thallium", 81, 204.3833*g/mole, 11.85*g/cm3);
	mats[67].z=82; mats[67].g4mat=new G4Material
			       ("Lead", 82, 207.2*g/mole, 11.3437*g/cm3);
	mats[68].z=83; mats[68].g4mat=new G4Material
			       ("Bismuth", 83, 208.98037*g/mole, 9.747*g/cm3);
	mats[69].z=84; mats[69].g4mat=new G4Material
			       ("Polonium", 84,
				(209.9828)*g/mole, 9.32*g/cm3);
	mats[70].z=86; mats[70].g4mat=new G4Material
			       ("Radon", 86,
				(222.0176)*g/mole, 9.73*g/dm3);
	mats[71].z=87; mats[71].g4mat=new G4Material
			       ("Francium", 87,
				(223.0197)*g/mole, 1.87*g/cm3);
	mats[72].z=88; mats[72].g4mat=new G4Material
			       ("Radium", 88,
				(226.0254)*g/mole, 5.5*g/cm3);
	mats[73].z=89; mats[73].g4mat=new G4Material
			       ("Actinium", 89,
				(227.0278)*g/mole, 10.07*g/cm3);
	mats[74].z=58; mats[74].g4mat=new G4Material
			       ("Cerium", 58, 140.115*g/mole, 6.773*g/cm3);
	mats[75].z=59; mats[75].g4mat=new G4Material
			       ("Praseodymium", 59, 140.90765*g/mole, 6.773*g/cm3);
	mats[76].z=60; mats[76].g4mat=new G4Material
			       ("Neodymium", 60, 144.24*g/mole, 7.008*g/cm3);
	mats[77].z=61; mats[77].g4mat=new G4Material
			       ("Promethium", 61,
				(144.9127)*g/mole, 7.26*g/cm3);
	mats[78].z=62; mats[78].g4mat=new G4Material
			       ("Samarium", 62, 150.36*g/mole, 7.520*g/cm3);
	mats[79].z=63; mats[79].g4mat=new G4Material
			       ("Europium", 63, 151.965*g/mole, 5.244*g/cm3);
	mats[80].z=64; mats[80].g4mat=new G4Material
			       ("Gadolinium", 64, 157.25*g/mole, 7.901*g/cm3);
	mats[81].z=65; mats[81].g4mat=new G4Material
			       ("Terbium", 65, 158.92534*g/mole, 8.230*g/cm3);
	mats[82].z=66; mats[82].g4mat=new G4Material
			       ("Dysprosium", 66, 162.50*g/mole, 8.551*g/cm3);
	mats[83].z=67; mats[83].g4mat=new G4Material
			       ("Holmium", 67, 164.93032*g/mole, 8.795*g/cm3);
	mats[84].z=68; mats[84].g4mat=new G4Material
			       ("Erbium", 68, 167.26*g/mole, 9.066*g/cm3);
	mats[85].z=69; mats[85].g4mat=new G4Material
			       ("Thulium", 69, 168.93421*g/mole, 9.321*g/cm3);
	mats[86].z=70; mats[86].g4mat=new G4Material
			       ("Ytterbium", 70, 173.04*g/mole, 6.966*g/cm3);
	mats[87].z=71; mats[87].g4mat=new G4Material
			       ("Lutetium", 71, 174.967*g/mole, 9.841*g/cm3);
	mats[88].z=90; mats[88].g4mat=new G4Material
			       ("Thorium", 90, 232.0381*g/mole, 11.72*g/cm3);
	mats[89].z=91; mats[89].g4mat=new G4Material
			       ("Protactinium", 91,
				(231.03588)*g/mole, 15.37*g/cm3);
	mats[90].z=92; mats[90].g4mat=new G4Material
			       ("Uranium", 92, 238.0508*g/mole, 18.05*g/cm3);
	mats[91].z=93; mats[91].g4mat=new G4Material
			       ("Neptunium", 93,
				(237.0482)*g/mole, 20.25*g/cm3);
	mats[92].z=94; mats[92].g4mat=new G4Material
			       ("Plutonium", 94,
				(244.0642)*g/mole, 19.84*g/cm3);
	mats[93].z=95; mats[93].g4mat=new G4Material
			       ("Americium", 95,
				(243.0614)*g/mole, 13.67*g/cm3);
	mats[94].z=96; mats[94].g4mat=new G4Material
			       ("Curium", 96,
				(247.0703)*g/mole, 13.51*g/cm3);
	mats[95].z=97; mats[95].g4mat=new G4Material
			       ("Berkelium", 97, (247.0703)*g/mole, 14*g/cm3);
		
	// initialise processes
	//
	G4PhotoElectricEffect* phot = new G4PhotoElectricEffect();
	G4ComptonScattering*   comp = new G4ComptonScattering();
	G4GammaConversion*     conv = new G4GammaConversion(); 

	// print cross section per atom
	//
	G4double Emin = 5*keV, Emax = 10*MeV, dE = 10*keV;
	material *iter = mats;
	double Z;
	for(int i=0; i<96; i++) {
		char buf[100];
		sprintf(buf, "%02d.txt", mats[i].z);
		Z = mats[i].z;
		ofstream file(buf);
		for (G4double Energy = Emin; Energy <= Emax; Energy += dE) {
			file << Energy/MeV << " " <<
				comp->ComputeCrossSectionPerAtom(Energy,Z)/barn +
	 			phot->ComputeCrossSectionPerAtom(Energy,Z)/barn + 
				conv->ComputeCrossSectionPerAtom(Energy,Z)/barn
			     << endl;
		}
	} 

	return EXIT_SUCCESS;
}

Для работы потребовалось определить все элементы, задействованные в вычислениях, т. к. это необходимо для работы G4PhotoElectricProcess.

Edit page