close
Warning:
Can't synchronize with repository "(default)" ("(default)" is not readable or not a Git repository.). Look in the Trac log for more information.
- Timestamp:
-
Dec 21, 2016, 11:33:11 AM (8 years ago)
- Author:
-
Clarence Wret
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v5
|
v6
|
|
4 | 4 | |
5 | 5 | |
6 | | == Compiling NEUT and running neutsmpl == |
| 6 | == Compiling NEUT and running neutroot2 == |
7 | 7 | |
8 | 8 | The main neutrino-nucleus interaction generator lives in $NEUT_ROOT/src/neutsmpl/neutroot. NEUT has two dependencies: ROOT and the lesser known CERNLIB. You need to have $ROOTSYS, $CERN and $CERN_LEVEL environment variables set-up to compile NEUT. These are all specified in $NEUT_ROOT/src/neutsmpl/Envmakeneutsmpl.csh. |
… |
… |
|
18 | 18 | == Structure of card files: == |
19 | 19 | |
| 20 | NEUT has a load of defaults set on run-time if the user does not supply them, which means you're running with "reference" NEUT. |
20 | 21 | |
| 22 | Some things are good to specify though, and here's a list for generating all modes on T2K H2O. We can specify how many events we want, the incoming neutrino type, the flux file, number of bound protons and neutrons and free protons, what interaction modes, and so on. |
| 23 | ``` |
| 24 | C Number of events ; EVCT-NEVT |
| 25 | EVCT-NEVT 100000 |
| 26 | C |
| 27 | C Particle Code ; EVCT-IDPT |
| 28 | EVCT-IDPT 14 |
| 29 | EVCT-MPOS 2 |
| 30 | EVCT-RAD 100. |
| 31 | EVCT-MDIR 1 |
| 32 | EVCT-DIR 0. 0. 1. |
| 33 | C |
| 34 | EVCT-MPV 3 |
| 35 | EVCT-FILENM 'nuisance_fluxes/t2kflux_2016_plus250kA.root' |
| 36 | EVCT-HISTNM 'enu_nd280_numu' |
| 37 | EVCT-INMEV 0 |
| 38 | C |
| 39 | C **** TARGET INFORMATION **** |
| 40 | NEUT-NUMBNDN 8 |
| 41 | NEUT-NUMBNDP 8 |
| 42 | NEUT-NUMFREP 2 |
| 43 | NEUT-NUMATOM 16 |
| 44 | C |
| 45 | C **** WHAT INTERACTION MODES DO WE WANT? (0 = ALL) **** |
| 46 | NEUT-MODE 0 |
| 47 | C **** WHAT CCQE MODEL ARE WE USING (2 = RFG WITH BBBA05) *** |
| 48 | NEUT-MDLQE 2 |
| 49 | C |
| 50 | C **** RANDOM NUMBER AND IF WE WANT NEUT TO BE QUIET **** |
| 51 | NEUT-RAND 1 |
| 52 | NEUT-QUIET 2 |
| 53 | ``` |
21 | 54 | |
22 | 55 | = NEUT inputs for NUISANCE = |
… |
… |
|
29 | 62 | }}} |
30 | 63 | |
| 64 | |
31 | 65 | = Using NEUT \ReWeight = |
32 | 66 | |