Warning, /eic.github.io/_software/djangoh.md is written in an unsupported language. File is not indexed.
0001 ---
0002 title: DJANGOH
0003 name: djangoh
0004 category: djangoh
0005 layout: default
0006 ---
0007 {% include layouts/title.md %}
0008
0009 * TOC
0010 {:toc}
0011
0012 The event generator [DJANGOH](http://www.desy.de/~hspiesb/mcp.html) simulates deep inelastic lepton-proton scattering for both NC and CC events including both QED and QCD radiative effects.
0013
0014 DJANGOH contains the Monte Carlo program HERACLES and an interface of HERACLES to LEPTO. The use of HERACLES allows to take into account the complete one-loop electroweak radiative corrections and radiative scattering.
0015 The LUND string fragmentation as implemented in the event simulation program [JETSET](http://home.thep.lu.se/~torbjorn/Pythia.html) is used to obtain the complete hadronic final state. At low hadronic mass, [SOPHIA](http://ebl.stanford.edu/) is used instead of LEPTO.
0016 DJANGOH comprises the programs (formerly kept separately) DJANGO6 and HERACLES. The interface is to version 6.5.1 of [LEPTO](http://www.isv.uu.se/thep/lepto/).
0017
0018 For the EIC, DJANGOH was upgraded to use nuclear PDFs as available in [LHAPDF](https://lhapdf.hepforge.org/).
0019 From version 4.6.10 on, DJANGOH simulates also longitudinal polarised deep inelastic lepton-proton scattering for both NC and CC events including both QED and QCD radiative effects.
0020
0021
0022 ##### Manuals
0023 Detailed documentation on DJANGOH is available from the following sources:
0024 * {% include documents/doc.md name='Djangoh_m.pdf' tag='Manual for DJANGOH 4.6.8' -%}
0025 * {% include documents/doc.md name='Djangoh_Updateh-4.6.10.pdf' tag='How are nuclear PDFs and polarised effects included in Djangoh' -%}
0026 * {% include documents/doc.md name='Lepto-6.5.1.pdf' tag='LEPTO Manual 6.5.1' -%}
0027
0028
0029 ##### Output file structure
0030
0031 The output file is in a text format which has the following structure:
0032 * 1st line: <tt>DJANGOH EVENT FILE</tt>
0033 * 2nd line: <tt>============================================</tt>
0034 * 3rd line: Information on event wise variables stored in the file:
0035
0036 | I: | 0 \(line index\) |
0037 |---------------------------------------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
0038 | ievent: | eventnumber running from 1 to XXX |
0039 | IChannel: | contains information on the origin of the event from the channels of HERACLES \(radiative corrections\) |
0040 | process \(LST 23\): | specifies process simulated \(=2: weak charged current; =4: neutral current\) |
0041 | subprocess \(LST 24\): | information about first\-order QCD process in current event \(=1: q\-event, =2: qg\-event, =3: q¯q\-event\) |
0042 | nucleon \(LST 22\): | information about chosen target nucleon in current event |
0043 | struckparton \(LST 25\): | information about flavor of the struck quark in current event |
0044 | partontrack \(LST 26\): | entry line in event record of outgoing struck quark\. In parton shower case, quark at boson vertex before final state shower; |
0045 | y, Q2, x, W2, Nu: | are the kinematic variables of the event as calculated from the scattered lepton for NC events\. For CC events it is calculated from the electrons and neutrino\. Includes effects from radiative corrections |
0046 | truey, trueQ2, truex, trueW2, trueNu: | are the kinematic variables of the event at the hard scattering vertex\. Uses the hadronic information \(CC: W\+neutrino\)\. Radiative effects are fully excluded\. |
0047 | SIGtot, errSIGtot | total cross section and its uncertainty |
0048 | D | depolarisation factor |
0049 | F1NC, F3NC | the unpolarised structure functions F1NC~\(q\+antiq\) and F3NC~\(q\-antiq\) for neutral current events |
0050 | G1NC, G3NC, A1NC | the polarised structure functions G1NC~\(dq\+dantiq\) and G3NC~\(dq\-dantiq\) and the simulated asymmetry A1NC=F1NC/G1NC for neutral current events |
0051 | F1CC, F3CC | the unpolarised structure functions for charge current events with W\-: F1CC=\(uptype\+antidowntype\) and F3CC=2\*\(uptype\-antidowntype\); with up <\-> down for W\+ |
0052 | G1CC, G5CC | the polarised structure functions for charge current events with W\-: G1CC=\(duptype\+dantidowntype\) and G5CC=\(\-duptype\+dantidowntype\); with up <\-> down for W\+ |
0053 | nrTracks: | number of tracks in this event, includes also virtual particles |
0054 {:.table-bordered}
0055 {:.table-striped}
0056
0057 * 4th line: <tt>============================================</tt>
0058 * 5th line: Information on track-wise variables stored in the file:
0059
0060 | I: | line index, runs from 1 to nrTracks |
0061 |-----------|-----------------------------------------------------------------------------------------------------------|
0062 | K\(I,1\): | status code KS \(1: stable particles 11: particles which decay 55; radiative photon\) |
0063 | K\(I,2\): | particle KF code \(211: pion, 2112:n, \.\.\.\.\) |
0064 | K\(I,3\): | line number of parent particle |
0065 | K\(I,4\): | normally the line number of the first daughter; it is 0 for an undecayed particle or unfragmented parton |
0066 | K\(I,5\): | normally the line number of the last daughter; it is 0 for an undecayed particle or unfragmented parton\. |
0067 | P\(I,1\): | px of particle |
0068 | P\(I,2\): | py of particle |
0069 | P\(I,3\): | pz of particle |
0070 | P\(I,4\): | Energy of particle |
0071 | P\(I,5\): | mass of particle |
0072 | V\(I,1\): | x vertex information |
0073 | V\(I,2\): | y vertex information |
0074 | V\(I,3\): | z vertex information |
0075 {:.table-bordered}
0076 {:.table-striped}
0077 <br />
0078
0079 * 6th line: <tt>============================================</tt>
0080 * 7th line: event information for first event
0081 * 8th line: <tt>============================================</tt>
0082 * 9th to X-1 line: track-wise info of 1st event
0083 * Xth line <tt>============================================</tt>
0084
0085 **For each subsequent event, lines 7 through X repeat analogously.**
0086
0087 ##### Unpolarised Example Input
0088
0089 Here is an example input file with settings for unpolarised inelastic neutral-current interactions with radiative corrections turned on.The manual explains the meaning of each input parameter and the associated values.
0090
0091 The differences between DJANGOH 4.6.9 and 4.6.10 are described in detail * {% include documents/doc.md name='Djangoh_Updateh-4.6.10.pdf' tag='here' %}
0092
0093
0094 ```
0095 OUTFILENAM
0096 djangoh_20x250_ep.Rad=1.NLO.NC
0097 TITLE
0098 DJANGOH 4.6.10 for eRHIC for PR, NLO at 20x250, Wmin=1.4
0099 EL-BEAM
0100 20D0 0.0D0 -1
0101 PR-BEAM
0102 250D0 0.0D0
0103 GSW-PARAM
0104 2 1 3 1 1 1 2 1 1 1 1
0105 KINEM-CUTS
0106 3 0D0 1.00D0 0.01D0 0.95D0 1.0D2 1.0D5 1.4D0
0107 EGAM-MIN
0108 0D0
0109 INT-OPT-NC
0110 1 18 18 18 18 0 0 0 0
0111 INT-OPT-CC
0112 0 0 0 0
0113 INT-ONLY
0114 1
0115 INT-POINTS
0116 30000
0117 SAM-OPT-NC
0118 1 1 1 1 1 0 0 0 0
0119 SAM-OPT-CC
0120 0 0 0 0
0121 NUCLEUS
0122 250D0 1D0 1D0
0123 NUCL-MOD
0124 0
0125 STRUCTFUNC
0126 0 2 10150
0127 POLPDF
0128 0
0129 LHAPATH
0130 /afs/rhic.bnl.gov/eic/bin/LHAPDF-5.8.6/share/lhapdf/PDFsets
0131 FLONG
0132 111 0.01 0.03
0133 ALFAS
0134 1 0 0.20 0.235
0135 FLAVORS
0136 0 5
0137 IOUNITS
0138 6 10 11
0139 RNDM-SEEDS
0140 -1 -1
0141 START
0142 100000
0143 SOPHIA
0144 1.5
0145 OUT-LEP
0146 1
0147 FRAG
0148 1
0149 CASCADES
0150 12
0151 MAX-VIRT
0152 5
0153 CONTINUE
0154 ```
0155
0156 To turn radiative corrections off, the following change needs to be made:
0157 ```
0158 SAM-OPT-NC
0159 1 0 0 0 0 0 0 0 0
0160 INT-OPT-NC
0161 1 0 0 0 0 0 0 0 0
0162 GSW-PARAM
0163 2 0 3 1 0 0 2 1 1 1 1
0164 ```
0165
0166 To run for nuclear targets the following changes need to be made. Any nucleus available in a specific nuclear pdf can be simulated, i.e. D, He, ..., Au,
0167 ```
0168 NUCLEUS
0169 100.0D0 197 79
0170 NUCL-MOD
0171 4201
0172 ```
0173
0174 For elastic (e + p → e + p) and quasi-elastic (e + p → e + p + γ) neutral current interactions the ```-OPT-NC``` settings should be set as follows:
0175 ```
0176 INT-OPT-NC
0177 0 0 0 0 0 1 20 20 20
0178 SAM-OPT-NC
0179 0 0 0 0 0 1 1 1 1
0180 ```
0181
0182 For charged current running, the ```-OPT-NC``` and ```-OPT-CC``` settings should be set as follows:
0183 ```
0184 INT-OPT-NC
0185 0 0 0 0 0 0 0 0 0
0186 INT-OPT-CC
0187 1 20 0 20
0188 SAM-OPT-NC
0189 0 0 0 0 0 0 0 0 0
0190 SAM-OPT-CC
0191 1 1 0 1
0192 ```
0193
0194 ##### Polarised Example Input
0195
0196 Here is an example input file with settings for polarized inelastic neutral-current interactions with radiative corrections turned on.
0197 ATTENTION: <br>
0198 to generate W+ exchange: ```EL-BEAM``` needs to be Energy +1.0D0 +1<br>
0199 to generate W- exchange: ```EL-BEAM``` needs to be Energy -1.0D0 -1<br>
0200 the hadron beam can have both helicity states + and -
0201
0202 ```
0203 OUTFILENAM
0204 djangoh_20x250_proton
0205 TITLE
0206 DJANGOH 4.6.10 for eRHIC for PR, LO at 20x250, Wmin=1.4
0207 EL-BEAM
0208 20D0 -1.0D0 -1
0209 PR-BEAM
0210 250D0 -1.0D0
0211 GSW-PARAM
0212 2 1 3 1 1 1 2 1 1 1 1
0213 KINEM-CUTS
0214 3 0D0 1.00D0 0.01D0 0.95D0 1.0D2 1.0D5 1.4D0
0215 EGAM-MIN
0216 0D0
0217 INT-OPT-NC
0218 1 18 18 18 18 0 0 0 0
0219 INT-OPT-CC
0220 0 0 0 0
0221 INT-ONLY
0222 1
0223 INT-POINTS
0224 30000
0225 SAM-OPT-NC
0226 1 1 1 1 1 0 0 0 0
0227 SAM-OPT-CC
0228 0 0 0 0
0229 NUCLEUS
0230 250D0 1 1
0231 NUCL-MOD
0232 0
0233 STRUCTFUNC
0234 0 2 21000
0235 POLPDF
0236 100
0237 LHAPATH
0238 /afs/rhic.bnl.gov/eic/bin/LHAPDF-5.8.6/share/lhapdf/PDFsets
0239 FLONG
0240 0 0.01 0.03
0241 ALFAS
0242 1 0 0.20 0.235
0243 NFLAVORS
0244 0 5
0245 IOUNITS
0246 6 10 11
0247 RNDM-SEEDS
0248 -1 -1
0249 START
0250 100000
0251 SOPHIA
0252 1.5
0253 OUT-LEP
0254 1
0255 FRAG
0256 1
0257 CASCADES
0258 1
0259 MAX-VIRT
0260 5
0261 CONTINUE
0262 ```
0263
0264 There is also the option to simulate a polarized neutron target by using
0265 ```
0266 NUCLEUS
0267 250D0 1D0 0D0
0268 ```
