Gabriel N. Perdue Fermilab - Genie
Contents
1. e Computing the total cross section is computationally expensive but it doesn t change event to event so do it once for all the materials in your experiment and store the results e Generate interactions e Run validation apps Gabriel N Perdue Fermilab NuSTEC School Liverpool n Qs y ar gmksp lL e GENIE Make Spline e Example gmkspl p 14 14 t 1000060120 o ccge carbon splines xml event generator list CCQE e Set the projectile PDG code the target PDG Ion Code the Output and the event generator List e See GENIE config EventGeneratorListAssembler xml for a list of valid event generator lists e You will almost always want the default but you may want to produce a narrower set of channels for specific studies e Event generators themselves are defined in GENIE config UserPhysicsOptions xml You will see this file a lot Gabriel N Perdue Fermilab NuSTEC School Liverpool Qs Yeas EventGeneratorListAssembler xml pe slide lt param_set name CCQE gt lt param type int name NGenerators gt 1 lt param gt lt param type alg name Generator 0 gt genie EventGenerator QEL CC lt param gt lt param_set gt A Defined in EventGenerator xml UserPhysicsOptions xml oo NINININI OSD ODDS OD NINI NINININI OD OD OAD OD OSD ODORS OD OD OND OD OD OD NINININI NI Specify which cross section model is to be used by each GENIE event generation thread The p2. e LwlynSmithQELCCPXSec Go to GENIE src S find name LwlynSmithOELCCPXSec LlewellynSmith LwlynSmithQELCCPXSec cxx LlewellynSmith LwlynSmithOELCCPXSec h Gabriel N Perdue Fermilab 22 NuSTEC School Liverpool LoadConf ig S void LwlynSmithQELCCPXSec LoadConfig void AlgConfigPool confp AlgConfigPool Instance const Registry x gc confp gt GlobalParameterList double thc fConfig gt GetDoubleDef CabbiboAngle gc gt GetDouble CabbiboAngle fCos8c2 TMath Power TMath Cos thc 2 load QEL form factors model fFormFactorsModel dynamic_cast lt const QELFormFactorsModell x gt this gt SubAlg FormFactorsAlg assert fFormFactorsModel fFormFactors SetModel fFormFactorsModel lt attach algorithm load XSec Integrator fXSecIntegrator dynamic_cast lt const XSecIntegratorI zs this gt SubAlg XSec Integrator assert fXSecIntegrator 23 Gabriel N Perdue Fermilab NuSTEC School Liverpool Jt af If we grep the config directory for FormFactorsAlg LwlynSmithQELCCPXSec LwlynSmithQELCCPXSec LwlynSmithQELCCPXSec LwlynSmithQELCCPXSec Xml Xml xml Xml lt param type lt param type lt param type lt param type alg alg alg alg name FormFactorsAlg name FormFactorsAlg name FormFactorsAlg name FormFactorsAlg Vv M MM genie genie genie genie we LwlynSmithFFCC Default Lwl
3. 1000080160 016 gt path length 1 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx ComputeProbScales 719 gt Computing the max interaction probability probability scale 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx ComputeProbScales 810 gt Probability scale 6 78377e 15 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx Configure 487 gt Finished configuring GMCJDriver 30 Gabriel N Perdue Fermilab NuSTEC School Liverpool Je n w Flux Neutrinos e Once we know the probability scale we begin drawing flux neutrinos and see if they interact most will not but thanks to the scaling it will not require an overly Large number of neutrinos to get an interaction 1396387998 NOTICE gevgen n lt gEvGen cxx GenerateEventsUsingFluxOrTgtMix 328 gt zk Generating event H 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateEvent 823 gt Generating next event 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateFluxNeutrino 1022 gt Generating a flux neutrino 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateFluxNeutrino 1036 gt Generated flux neutrino o PDG code 14 o 4 momentum E 0 327306 px 0 py 0 pz 0 327306 M 2 P 0 327306 o 4 position t 0 x 0 y z 0 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateEventiTry 875 gt Computing intera
4. lt E gt lt xsec gt 4 103713443e 12 lt xsec gt lt knot gt lt knot gt lt E gt 0 2303 lt E gt lt xsec gt 8 36 6966 79e 12 lt xsec gt lt knot gt lt knot gt lt E gt 0 29257 lt E gt lt xsec gt 1 290439962e 11 lt xsec gt lt knot gt lt knot gt lt E gt 0 37168 lt xsec gt lt knot gt Gabriel N Perdue Fermilab lt E gt 1 Ee 11 Cross Section lt xsec gt NuSTEC School Liverpool Oy Ni 3 Gr e Algorithm and class names follow the convention of naming the file after the class so they are easy to find S cd SGENIE src find name LwlynSmithOELCCPXSec LlewellynSmith LwlynSmithQELCCPXSec cxx LlewellynSmith LwlynSmithOELCCPXSec h When running find and grep you may want to mask certain directories e g svn with find name xfilex grep v svn etc Gabriel N Perdue Fermilab NuSTEC School Liverpool n VK e More on Splines e Generating splines takes a significant amount of CPU time recall we need a computation over dozens of energy positions dozens of models and dozens of targets e You can find some useful pre built splines here e https www hepforge org archive genie data 2 8 0 e And at Fermilab read e https cdcvs fnal gov redmine projects genie wiki UPS products genie phyopt and genie xsec Gabriel N Perdue Fermilab 7 NuSTEC School Liverpool Other Spline Utilities Comments e See the Physi
5. may have built ROOT without gdml enabled It is a library I didn t know was needed Goto ROOT Add enable gdml to the config status file make distclean then configure cat config status and make Then goto GENIE and make clean amp amp make If you installed after May 11 or so the Scripts included this library in the configure Fermilab NuSTEC School Liverpool n LVL y ar Installation e Installation on Linux systems with recent versions of GCC is very straightforward e We do not exercise builds on Windows very often and don t really support it e Installation on a newer Mac or older Mac with newer 5 x Xcode is tricky because Apple has stopped supporting the GNU toolchain and GENIE does not build using the version of clang that comes with Xcode 5 x The fix for clang is possibly understood but not implemented yet e One fix is to install GCC yourself You will need gfortran also See the following link and check for the subsection on building on OSX for ideas on how to do this we plan on a more official fix this summer e http home fnal gov perdue GENTEDeve LopersManual html _building genie Gabriel N Perdue Fermilab NuSTEC School Liverpool Oy Ni 3 Gr Running GENIE e Three basic modes e Generate interaction splines These are files that store the total cross section as a function of channel and neutrino energy on different targets
6. tgt 1000080160 N 2212 proc Weak CC RES res 0 6 1741 vee etc 3 Gabriel N Perdue Fermilab NuSTEC School Liverpool n Ve Gr II II e Discussed programatically later But how DD ne E E E GENIE GHEP Event Record print level 31 a nnn nnn enden ne nnn nei en Idx Name Ist PDG Mother Daughter Px Py Pz E m Keeser numu ol 14 1 1 4 4 0 000 0 000 3 000 3 000 0 000 1 C12 1000060120 1 1 2 3 0 000 0 000 0 000 11 179 11 179 2 neutron 11 2112 1 1 5 5 0 054 0 160 0 080 0 934 940 M 0 915 3 C11 2 1000060110 1 1 12 12 0 054 160 0 080 10 245 10 243 4 mu 1 13 0 1 1 1 0 133 0 656 0 513 0 849 0 106 P 0 158 0 778 0 608 5 HadrSyst 12 2000000001 2 1 6 7 0 188 0 495 2 407 3 084 000 M 1 854 6 neutron 14 2112 5 1 8 19 0 670 0 529 1 086 1 670 0 940 FSI 7 7 pi 14 211 5 1 11 11 0 482 0 034 1 321 1 414 0 140 FSI 1 8 neutron 1 2112 6 1 1 1 0 014 0 186 0 749 1 216 0 940 9 proton 1 2212 6 1 1 1 0 653 0 408 0 296 1 249 0 938 10 pi 1 211 6 1 1 1 0 002 0 057 0 024 0 153 0 140 11 pi 1 211 7 1 1 1 0 482 0 034 1 321 1 414 0 140 12 HadrBlob 15
7. 11 7 1 1 2000000002 3 1 1 nu_mu x 0 00000 m y Err flag bits 15 gt 0 0000000000000000 Ist set n R Validation Apps Compare predictions from GENIE to data e There is an extensive database of experimental results stored as ASCII plain text files distributed with GENIE along with citation information More on this ina Later Lecture Gabriel N Perdue Fermilab 36 NuSTEC School Liverpool e Install GENIE using instructions in the Physics and User s Manual or a script package like lamp e gmkspl is the primary app for producing cross section summary splines These files are used for event selection during event generation e gevgen creates events e Running either app without arguments will produce a usage guide e Both apps feature detailed logs With default verbosity you may study the details of interactions during event generation Gabriel N Perdue Fermilab 37 NuSTEC School Liverpool Gabriel N Perdue Fermilab Thanks NuSTEC School Liverpool Gabriel N Perdue Fermilab Back up NuSTEC School Liverpool gmksp lL e What happens when we run gmkspl e See instructions later in this presentation for a gdb walkthrough during event generation S gmkspl p 14 14 t 1000060120 o test xml event generator list CCCOH e For test runs it is a good idea to use a single element for your target and limit the spec
8. 2000000002 3 1 1 1 0 049 0 038 0 096 9 297 0 000 M 9 297 a _ Fin Init 0 000 0 000 0 000 0 000 ee a Vertex nu_mu x 0 00000 m y 0 00000 m z 0 00000 m t Os ee Err flag bits 15 gt 0 0000000000000000 1st set none Err mask bits 15 gt 0 1111111111111111 Is unphysical NO Accepted YES EE HH sig Ev 8 1385e 38 cm 2 d2sig x y E dxdy 2 4893e 37 cm 2 Weight 1 ee EE EE 32 Gabriel N Perdue Fermilab NuSTEC School Liverpool Qe A Re gt GHepRecords e We also get an interaction summary more on this in a later lecture Init State gt probe PDG code gt nucl target Z 6 A gt hit nucleon PDC Code 14 nu_mu 12 PDG Code 1000060120 C12 2112 neutron gt hit quark PDC Code 1 d valence gt probe 4P E 3 Px 0 Py 0 Pz 3 gt target 4P CG ZE 11 179 Px 0 Py 0 Pz 0 gt nucleon 4P E 0 933909 Px 542517 Py 0 160284 Pz 0 0800058 Process Info gt Interaction Weak CC gt Scattering DIS Kinematics gt Running Hadronic invariant mass W gt Selected Bjorken x 0 435806 gt Selected Inelasticity y 0 757872 gt Selected Momentum transfer Q2 gt 0 gt Selected Hadronic invariant mass W Exclusive Process Info gt charm prod false gt f
9. 9 297 memmen E i E a ee ee ee ee SS Fin Init 0 000 0 000 0 000 gt Vertex nu_mu x 0 00000 m y 0 00000 m z 0 00000 m t et ee Se PP ee ee ee ee Fa Err flag bitg A5 gt 0 0000000000000000 1st set Err mask bit6f15 gt 0 1111111111111111 Is unphysical NO Accepted YES sig Ev 8 1385e 38 cm 2 d2sig x y E dxdy 2 4893e 37 cm 2 Weight 1 s go to GEANT See Not enumerated Src GHEP GHepStatus h 34 Gabriel N Perdue Fermilab NuSTEC School Liverpool Jt af typedef enum EGHepStatus kIStUndef ined kIStInitialState kIStStableFinalState kIStIntermediateState kIStDecayedState KIStCorrelatedNuc leon KIStNucleonTarget kIStDISPreFragmHadronicState kIStPreDecayResonantState KIStHadronInTheNucleus KIStFinalStateNuclearRemnant KIStNucleonClusterTarget GHepStatus_t Gabriel N Perdue Fermilab H generator level initial state 1 4 generator level final state particles to be tracked by detector level MC x 14 N hadrons inside the nucleus marked for hadron transport modules to act on 15 low energy nuclear fragments entering the record as a pseudo particle 16 3 PDG Mother Daug 14 1 1 4 1000060120 1 1 2 2112 1 1 5 1000060110 1 1 12 13 0 1 1 2000000001 2 1 6 2112 5 1 8 211 5 1 11 2112 6 1 1 2212 6 1 1 211 6 1 1 2
10. ExcitationSim Default HadronTransporter Default NucBindEnergyAggregator Default Unstab leParticleDecayer AfterHadronTransport QELInteractionListGenerator CC Default Ks Jl D y as lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt lt param gt When running a QEL CC event these are the modules we step through Note that this is specified in text not in code and may be changed without re compiling and even on the fly during a run Where do we find them NuSTEC School Liverpool Qs y as JE En CJ You can often find a lookup table in master_config xml lt genie_config gt lt xkkkkk CONFIGURATION FOR EVENT GENERATION MODULES kkk gt lt config alg genie EventGenerator gt EventGenerator xml lt config gt lt config alg genie FermiMover gt FermiMover xml lt config gt lt config alg genie HadronTransporter gt HadronTransporter xml lt config gt lt config alg genie QELKinematicsGenerator gt QELKinematicsGenerator xml lt config gt And extra information if you follow the table here QELKinematicsGenerator xml etc Typically though these files also correspond to class names by convention ls 1 src QEL Q h cxx src QEL QELHadronicSystemGenerator cxx src QEL QELHadronicSystemGenerator h src QEL QELIn
11. Introduction to GENIE Gabriel N Perdue Fermilab 2 Fermilab ENERGY Oy Ni 3 Gr Overview Getting help Installation e Run the basic applications e Simple output checks The goal for this lecture install GENIE make a cross section spline file generate some events A greenfield build of GENIE will take too long if we are only starting now but start the build anyway Hopefully you already have it in which case focus on making splines and interactions Gabriel N Perdue Fermilab NuSTEC School Liverpool Prefatory Remarks e These presentations attempt to be consistent in how they refer to files functions and directories e When I remember to I use back ticks to mark functions classes compiled programs and code snippets and single ticks to mark files and directories e Directory paths are generally understood to begin with GENIE the installation root environment var e Directory paths into the source code sometimes start with src and sometimes don t If it isn t clear ask e Configuration xml will usually be assumed to be found in GENIE config e Bash shell expansion will be used in the text so File h cxx is File h File cxx as would resolve under Let Commands will often but not always start with a dollar prompt and use a different Courier font Gabriel N Perdue Fermilab NuSTEC School Liverpool Getting Help e http genie hepforg
12. arameter name is build as XSecModel name of thread gt lt param type alg name XSecModel genie EventGenerator QEL CC gt genie LwlynSmithQELCCPxXSec DefauLt lt param gt Gabriel N Perdue Fermilab 2 NuSTEC School Liverpool Jt af To drill down further SCE conf ig EventGenerator xml lt param set name QEL CC gt type string name VldContext gt lt param gt lt param lt param lt param lt param lt param lt param lt param lt param lt param lt param lt param lt param lt param lt param lt param type int type alg type alg type alg type alg type alg type alg type alg type alg type alg type alg type alg type alg type alg lt param_set gt Gabriel N Perdue Fermilab name NModules gt 12 name Modu le 0 gt name Modu le 1 gt name Module 2 gt genie name Module 3 gt genie name Module 4 gt genie name Modu Le 5 gt name Module 6 gt genie name Module 7 gt genie name Module 8 gt genie name Module 9 gt name Module 10 gt genie name Module 11 gt genie name ILstGen gt genie genie genie genie genie InitialStateAppender Default VertexGenerator Default FermiMover Default QELKinematicsGenerator CC Default QELPrimaryLeptonGenerator Default QELHadronicSystemGenerator Default PauliBlocker Default UnstableParticleDecayer BeforeHadronTransport NucDe
13. cs and User s Manual for details e gspladd merge spline files e gspl2root convert XML cross section files in a ROOT file e Currently the full cross section is loaded at the beginning of a run There is some work on going for making lazy loading possible but it is not available yet Gabriel N Perdue Fermilab 18 NuSTEC School Liverpool n Jt af eS More on Configuration e Recall the List of modules that appeared in EventGenerator xml for QEL CC e How do we configure the associated parameters Example How to turn off intranuclear rescattering Go to GENIE src Ss find name HadronTransporter EVGModules HadronTransporter cxx EVGModules HadronTransporter h Gabriel N Perdue Fermilab NuSTEC School Liverpool n JC af LoadConf ig GE void HadronTransporter LoadConfig void AlgConfigPool confp AlgConfigPool Instance const Registry x gc confp gt GlobalParameterList fHadTranspModel Q fEnabled fConfig gt GetBoolDef enable gc gt GetBool HadronTransp Enable LOG HadTransp pDEBUG lt lt Hadron transport was lt lt fEnabled not lt lt enabled if fEnabled RgAlg hadtransp_model fConfig gt GetAlgDef model gc gt GetAlg HadronTransp Model LOG HadTransp pDEBUG lt lt Loading the hadron transport model lt lt hadtransp_model Note this relies on the fact that if a miss
14. ction probabilities for max path lengths 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx ComputeInteractionProbabilities 1092 gt Computing relative interaction probabilities for each material 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateEventiTry 880 gt The no interaction probability max path lengths is 97 2556 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateEventiTry 890 gt xx Rejecting current flux neutrino 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateEvent 839 gt Flux neutrino didn t interact Trying the next one e Eventually we get an interaction and the actual channel is drawn from a cross section weighted list of all possible reactions unphysical reactions will have cross section weight Q 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GenerateEventKinematics 1210 gt Asking the selected GEVGDriver object to generate an event 1396387998 NOTICE IntSel n lt PhysInteractionSelector cxx SelectInteraction 98 gt 1396387998 NOTICE IntSel n lt PhysInteractionSelector cxx SelectInteraction 100 gt Computing xsecs for all relevant modeled interactions 1396387998 NOTICE IntSel n lt PhysInteractionSelector cxx SelectInteraction 170 gt nu 14 tgt 1000080160 N 2112 proc Weak CC QES 7 41357 nu 14 tgt 1000080160 N 2212 proc Weak NC QES 1 24454 nu 14 tgt 1000080160 N 2112 proc Weak NC QES 1 7118 nu 14
15. e org e For most user and physics questions check the Physics and Users Manual http genie hepforge org manuals GENIE PhysicsAndUserManual 20130615 pdf e When the manual is updated just check the left hand side of the GENIE main page e There is also a mailing list https www 1iscmail ac uk cgi bin webadmin AQ NEUTRINO MC SUPPORT e We are also working on a Developer s Manual aimed at implementing tuning models http home fnal gov perdue GENIEDevelopersManual html e When it is done it will live on the GENIE main page e You can look at the code it is written in AsciiDoc which is a mostly literate markup language that compiles to HTML and PDF and modify it and submit pull requests on GitHub https github com GENIEMC GENIEDeve LopersManual Gabriel N Perdue Fermilab NuSTEC School Liverpool Framework Internals e We won t discuss the framework internals here In Gabriel N Perdue Fermilab general you don t need to know much about them to use GENIE and develop and tune physics models e There are a couple of organizing patterns we need to understand We ll go over those in this school If you are really interested in framework internals see e Section Five of the NIMx e http projects docdb fnal gov cgi bin ShowDocument docid 2922 Andreopoulos C and Bell A and Bhattacharya D and Cavanna F and Dobson J and others The GENIE Neutrino Monte Carlo Generat
16. fig gt GetDoubleDef WeinbergAngle gc gt GetDouble WeinbergAngle fSin28w TMath Power TMath Sin thw 2 Finally going back to the config directory we grep for QEL Ma 25 Gabriel N Perdue Fermilab NuSTEC School Liverpool n VK Sei J y ar In UserPhysicsOptions xml ql Axial and vector masses for quasi elastic scattering gt lt param type double name QEL Ma gt 0 990 lt param gt lt param type double name QEL Mv gt 0 840 lt param gt Whew e It can sometimes take a bit of find ing and grep ing to get what you are after but everything is there and organized this way for a reason e Note though that we could have found this quickly by searching UserPhysicsOptions for MA and Ma case sensitive For example suppose you wanted to implement a QE model that used different form factors Would it make sense for Ma to be located in QELKinematics Similarly would you want two copies of a parameter CC amp NC when there should be only one Gabriel N Perdue Fermilab 26 NuSTEC School Liverpool n R gevgen GENIE Event Generation e Example gevgen n 10 p 14 t 1000080160 e 1 10 r 100 f x exp x seed 398819 cross sections SPATH gxspl vA 2 8 0 xml e Many other options of course use a geometry more on this in a later lecture weight the events be careful specify a generator list etc Gabriel N Perdue Fer
17. ies it can take a very long time to iterate through all the reaction channels and complicated geometries may feature dozens of nuclear targets Gabriel N Perdue Fermilab 40 NuSTEC School Liverpool
18. ing registry entry was found above then it would be filled up by the default choice so that when the registry is looked up again the same key would point to the correct value Probably it would be better to convert the RgAlg to AlgId and query the AlgFactory directly fHadT ranspModel dynamic_cast lt const EventRecordVisitorI x gt this gt SubAlg model assert fHadTranspModel We see this pattern of reading from the AlgConfigPool and Registry quite frequently Gabriel N Perdue Fermilab 20 NuSTEC School Liverpool We search GENIE config ae grep HadronTransp Model HadronTransporter xml model alg Yes GPL HadronTransp Model UserPhysicsOptions xml lt param type alg name HadronTransp Model gt genie HAIntranuke Default lt param gt Drilling down in UserPhysicsOptions xml Si Intranuclear rescattering Use the HadronTransportEnable option to toggle intranuclear rescattering on off Also set the preferred hadron transport model Options include genie HAIntranuke Defau lt genie HNIntranuke Defau lt gt lt param type bool name HadronTransp Enable gt true lt param gt lt param type alg name HadronTransp Model gt genie HAIntranuke Default lt param gt All we need to do then is set the enable flag to false Gabriel N Perdue Fermilab 21 NuSTEC School Liverpool Oy N 3 Gr e Another example Ma for CCQE e Recall the cross section algorithm is
19. ist Gabriel N Perdue Fermilab 29 NuSTEC School Liverpool Jt af 1396387997 NOTICE GMCJDriver 1396387997 NOTICE GMCJDriver 1396387998 NOTICE GMCJDriver O gt LVL y ar Path length Scaling First GENIE will look at all the interactions it has cross section information for in the spline files and build a total cross section model Then it will look at the geometry and find the largest possible path length here it is 1 for a single oxygen nucleus Then it will scale the interaction probability so that a neutrino crossing the largest path has an interaction probability of 1 This scaling keeps the MC efficient and alters the rate we use up neutrinos when we are drawing from flux file representing a certain number of protons on target lt GMCJDriver cxx BootstrapXSecSplineSummation 675 gt Summing up splines to get total cross section for each init state lt GMCJDriver cxx BootstrapXSecSplineSummation 683 gt k Summing xsec splines for init state nu pdg 14 tgt pdg 1000080160 lt GMCJDriver cxx GetMaxPathLengthList 587 gt Querying the geometry driver to compute the max path length list n 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx BootstrapXSecSplineSummation 702 gt Finished summing all interaction xsec splines per initial state n n 1396387998 NOTICE GMCJDriver n lt GMCJDriver cxx GetMaxPathLengthList 592 gt Maximum path length list o code
20. milab 27 NuSTEC School Liverpool n CETA y ar gevgen What happens in detail when we run gevgen It is very worthwhile walking through with a debugger e g gdb See e http home fnal gov perdue GENTEDeve LopersManual html _running_the_code_the_ccqge_walkthrough gdb is also very useful for understanding new models under development without resorting to dozens of debug print statements The Developer s Manual still under construction also has some useful resources like gdbinit files and basic gdb and valgrind usage etc clone the repository to get all the files 28 Gabriel N Perdue Fermilab NuSTEC School Liverpool O gt CNE y ar Output How do we know the code succeeded e GENIE leaves voluminous log files use GENIE config Messenger xmLl to control verbosity the default setting may be too talkative for your needs e Under default settings there are pretty printed GHepRecords we may examine e Failures are usually clearly explained 1396638461 NOTICE Rndm n lt RandomGen cxx SetSeed 92 gt Setting default random number seed at random number generator initialization 65539 1396638461 NOTICE Rndm n lt RandomGen cxx SetSeed 92 gt Setting random number seed 2989819 1396638461 FATAL AppInit n lt AppInit cxx XSecTable 64 gt Input cross section file minerva app users perdue GENIE data carbon_oxygen_splines xml does not ex
21. or Nucl Instrum Meth A614 87 104 2010 gt NuSTEC School Liverpool Gabriel N Perdue n VK e Installation There are instructions in the Physics and Users manual Additionally there are instructions on the web e http genie hepforge org instaLlation html The only tricky parts are third party Libraries and compiler and OS support Generally GCC on xnix systems will work Fermilab NuSTEC School Liverpool n Ns In Installation S e If you download the GENIE source code you will find installation Scripts for the third party libraries in e GENIE src scripts build ext Also for LINUX ONLYx there is a full installation script hosted on our UNOFFICIAL scratch project site on GitHub e https github com GENIEMC lamp git clone https github com GENIEMC lamp git rub the lamp sh no arguments to see help rub the lamp sh p 6 lt this is probably what you want e Note that lamp uses slightly different scripts to install the third party packages than those included with GENIE A to do for us is to unify these e ROOT takes a while to build but installation is otherwise fairly quick The total time is less than one hour For Macs see http home fnal gov perdue GENIEDeve LopersManual html Gabriel N Perdue Fermilab 7 NuSTEC School Liverpool Gabriel N Perdue n R Short Side Note If you installed GENIE via lamp last week you
22. s nucleons N p N n gt f s pions N pi 0 N pi N pi gt resonance not set 1 85428 2 00929 1 85428 Gabriel N Perdue Fermilab 33 NuSTEC School Liverpool Sor A A gt Jt bin bg GHepRecords NEE These specify ranges State GENIE GHEP Event Record print level Idx Name Ist PDG 14 1 1 4 4 0 000 0 000 3 000 3 000 0 000 1 C12 o 1000060120 1 1 2 3 0 000 0 000 0 000 11 179 11 179 2 neutron 11 2112 1 1 5 5 0 054 0 160 0 080 0 934 940 M 0 915 3 C11 2 1000060110 1 1 12 12 0 054 0 160 0 080 10 245 10 243 4 mu 41 13 ol 1 1 1 0 133 0 656 0 513 0 849 0 106 P 0 158 0 778 0 608 5 HadrSyst 12 2000000001 2 1 6 7 0 188 0 495 2 407 3 084 0 000 M 1 854 6 neutron 14 2112 5 1 8 10 0 670 0 529 1 086 1 670 0 940 FSI 7 7 pi f _14 211 5 1 11 11 0 482 0 034 1 321 1 414 0 140 FSI 1 8 neutron 1 2112 6 1 1 1 0 014 0 186 0 749 1 0 9 1 2212 6 1 1 1 0 653 0 408 0 296 1 0 10 1 211 6 1 1 1 0 002 0 057 0 024 0 0 11 1 211 7 1 1 1 0 482 0 034 1 321 1 0 12 T5 2000000002 3 1 1 1 0 049 0 038 0 096 9 0
23. teractionListGenerator cxx src QEL QELInteractionListGenerator h src QEL QELKinematicsGenerator cxx src QEL QELKinematicsGenerator h src QEL QELPrimaryLeptonGenerator cxx src QEL QELPrimaryLeptonGenerator h Gabriel N Perdue Fermilab 14 NuSTEC School Liverpool Jt af Example Cross Section Spline gt 4 In general the spline you need for your experiment will use both gt neutrinos and antineutrinos and have many different models and targets lt xml version 1 0 encoding IS0 8859 1 gt lt generated by genie XSecSplineList SaveSplineList gt lt genie_xsec_spline_ list version 2 00 ES Le Cr Foy lt spline Cross Section Algorithn Danes genie S Tgt 1000060120 gem Target Carbon d proc Weak CC 0ES nknots 44 gt lt knot gt lt E gt 0 01000 lt E gt lt xsec gt Q lt xsec gt lt knot gt lt knot gt lt E gt 0 030466 lt E gt lt xsec gt 0 lt xsec gt lt knot gt lt knot gt lt E gt 0 050932 lt E gt lt xsec gt Q lt xsec gt lt knot gt lt knot gt lt E gt 0 071399 lt E gt lt xsec gt Q lt xsec gt lt knot gt lt knot gt lt E gt 0 091865 lt E gt lt xsec gt Q lt xsec gt lt knot gt lt knot gt lt E gt 0 11233 lt E gt lt xsec gt Q lt xsec gt lt knot gt lt knot gt lt E gt 0 1427 lt E gt lt xsec gt 8 487371629e 13 lt xsec gt lt knot gt lt knot gt lt E gt 0 18129
24. ynSmithFFCC Dipole LwlynSmithFFCC BBA2003 LwlynSmithFFCC BBA2005 Then searching src for LwlynSmithFFCC we find LlewellynSmith LwlynSmithFFCC cxx LlewellynSmith LwlynSmithFFCC h find lt param gt lt param gt lt param gt lt param gt There we find the computations are done in LwlynSmithFF Gabriel N Perdue Fermilab 24 NuSTEC School Liverpool n LoadConfig Ka void LwlynSmithFF LoadConfig void Load configuration data from its configuration Registry or global defaults to private data members fElFFModel AlgConfigPool confp AlgConfigPool Instance const Registry gc confp gt GlobalParameterList load elastic form factors model RgAlg form_factors_model fConfig gt GetAlgDef ElasticFormFactorsModel gc gt GetAlg ElasticFormFactorsModel fELFFModel dynamic_cast lt const ELFormFactorsModell zz this gt SubAlg ElasticFormFactorsModel assert fELFFModel fELFF SetModel fELFFModel axial mass and Fa q2 0 fMa fConfig gt GetDoubleDef Ma gc gt GetDouble QEL Ma Axial mass fFAQ fConfig gt GetDoubleDef FAQ gc gt GetDouble QEL FAQ FA q2 0 fMa2 TMath Power fMa 2 anomalous magnetic moments fMuP fConfig gt GetDoubleDef MuP gc gt GetDouble AnomMagnMoment P fMUN fConfig gt GetDoubleDef MuN gc gt GetDouble AnomMagnMoment N weinberg angle double thw fCon
Download Pdf Manuals
Related Search