PYTHIA8/pythia8140/include/LHAFortran.h

   1 // LHAFortran.h is a part of the PYTHIA event generator.
   2 // Copyright (C) 2010 Torbjorn Sjostrand.
   3 // PYTHIA is licenced under the GNU GPL version 2, see COPYING for details.
   4 // Please respect the MCnet Guidelines, see GUIDELINES for details.
   5
   6 // Header file for Fortran Les Houches Accord user process information.
   7 // LHAupFortran: derived class with the HEPRUP and HEPEUP Fortran info.
   8 // You are expected to supply the fillHepRup and fillHepEup methods.
   9
  10 #ifndef Pythia8_LHAFortran_H
  11 #define Pythia8_LHAFortran_H
  12
  13 #include "PythiaStdlib.h"
  14
  15 namespace Pythia8 {
  16
  17 //==========================================================================
  18
  19 // Give access to the HEPRUP and HEPEUP Fortran commonblocks.
  20
  21 #ifdef _WIN32
  22   #define heprup_ HEPRUP
  23   #define hepeup_ HEPEUP
  24 #endif
  25
  26 extern "C" {
  27
  28   extern struct {
  29     int idbmup[2];
  30     double ebmup[2];
  31     int pdfgup[2], pdfsup[2], idwtup, nprup;
  32     double xsecup[100], xerrup[100], xmaxup[100];
  33     int lprup[100];
  34   } heprup_;
  35
  36   extern struct {
  37     int nup, idprup;
  38     double xwgtup, scalup, aqedup, aqcdup;
  39     int idup[500], istup[500], mothup[500][2], icolup[500][2];
  40     double pup[500][5], vtimup[500],spinup[500];
  41   } hepeup_;
  42
  43 }
  44
  45 //==========================================================================
  46
  47 // A derived class with initialization information from the HEPRUP
  48 // Fortran commonblock and event information from the HEPEUP one.
  49
  50 class LHAupFortran : public LHAup {
  51
  52 public:
  53
  54   // Constructor.
  55   LHAupFortran() {}
  56
  57   // Routine for doing the job of setting initialization info.
  58   bool setInit() {
  59     // Call the routine that does the job.
  60     if (!fillHepRup()) return false;
  61     // Store beam and strategy info.
  62     setBeamA(heprup_.idbmup[0], heprup_.ebmup[0], heprup_.pdfgup[0],
  63       heprup_.pdfsup[0]);
  64     setBeamB(heprup_.idbmup[1], heprup_.ebmup[1], heprup_.pdfgup[1],
  65       heprup_.pdfsup[1]);
  66     setStrategy(heprup_.idwtup);
  67     // Store process info. Protect against vanishing cross section.
  68     for (int ip = 0; ip < heprup_.nprup; ++ip) {
  69       double xsec = max( 1e-10, heprup_.xsecup[ip]);
  70       addProcess( heprup_.lprup[ip], xsec, heprup_.xerrup[ip],
  71         heprup_.xmaxup[ip] );
  72     }
  73     // Store the beam energies to calculate x values later.
  74     eBeamA = heprup_.ebmup[0];
  75     eBeamB = heprup_.ebmup[1];
  76     // Done.
  77     return true;
  78   }
  79
  80   // Routine for doing the job of setting info on next event.
  81   bool setEvent(int idProcIn = 0) {
  82     // In some strategies the type of the next event has been set.
  83     hepeup_.idprup = idProcIn;
  84     // Call the routine that does the job.
  85     if (!fillHepEup()) return false;
  86     // Store process info.
  87     setProcess(hepeup_.idprup, hepeup_.xwgtup, hepeup_.scalup,
  88       hepeup_.aqedup, hepeup_.aqcdup);
  89     // Store particle info.
  90     for (int ip = 0; ip < hepeup_.nup; ++ip) addParticle(hepeup_.idup[ip],
  91       hepeup_.istup[ip], hepeup_.mothup[ip][0], hepeup_.mothup[ip][1],
  92       hepeup_.icolup[ip][0], hepeup_.icolup[ip][1], hepeup_.pup[ip][0],
  93       hepeup_.pup[ip][1], hepeup_.pup[ip][2], hepeup_.pup[ip][3],
  94       hepeup_.pup[ip][4], hepeup_.vtimup[ip], hepeup_.spinup[ip]) ;
  95     // Store x values (here E = pup[ip][3]), but note incomplete info.
  96     setPdf( hepeup_.idup[0], hepeup_.idup[1], hepeup_.pup[0][3]/eBeamA,
  97             hepeup_.pup[1][3]/eBeamB, 0., 0., 0., false);
  98     // Done.
  99     return true;
 100   }
 101
 102 private:
 103
 104   // Save beam energies to calculate x values.
 105   double eBeamA, eBeamB;
 106
 107   // User-written routine that does the intialization and fills heprup.
 108   bool fillHepRup();
 109
 110   // User-written routine that does the event generation and fills hepeup.
 111   bool fillHepEup();
 112
 113 };
 114
 115 //==========================================================================
 116
 117 } // end namespace Pythia8
 118
 119 #endif // Pythia8_LHAFortran_H