TFluka/mgdraw.cxx

   1 #include <Riostream.h>
   2 #include "TVirtualMCApplication.h"
   3 #include "TVirtualMCStack.h"
   4 #include "TFluka.h"
   5 // Fluka include
   6 #include "Fdimpar.h"  //(DIMPAR) fluka include
   7 #include "Fdblprc.h"  //(DBLPRC) fluka common
   8 #include "Ftrackr.h"  //(TRACKR) fluka common
   9
  10 #ifndef WIN32
  11 # define mgdraw mgdraw_
  12 #else
  13 # define mgdraw MGDRAW
  14 #endif
  15
  16 extern "C" {
  17 void mgdraw(Int_t& icode, Int_t& mreg)
  18 {
  19     TFluka* fluka =  (TFluka*) gMC;
  20     Int_t verbosityLevel = fluka->GetVerbosityLevel();
  21 //
  22 //  Make sure that stack has currrent track Id
  23     Int_t trackId = TRACKR.ispusr[mkbmx2-1];
  24     TVirtualMCStack* cppstack = fluka->GetStack();
  25     cppstack->SetCurrentTrack(trackId);
  26 //
  27 //
  28     Int_t oldreg = ((TFluka*) gMC)->GetMreg();
  29     if (oldreg != mreg) {
  30 //
  31 //  Boundary Crossing
  32 //
  33         fluka->SetNewreg(mreg);
  34         if (oldreg == -1) fluka->SetMreg(mreg);
  35         if (verbosityLevel >= 3)
  36             printf("Boundary Crossing %d %d \n", oldreg, mreg);
  37     } else {
  38         fluka->SetMreg(mreg);
  39         fluka->SetNewreg(mreg);
  40         if (verbosityLevel >= 3)
  41             printf("Normal step %d %d \n", oldreg, mreg);
  42     }
  43     fluka->SetIcode(icode);
  44     fluka->SetCaller(4);
  45
  46     if (verbosityLevel >= 3) {
  47         cout << endl << " !!! I am in mgdraw - calling Stepping()" << endl;
  48         cout << endl << " Track Id =" << trackId << endl;
  49     }
  50
  51     fluka->FutoTest();
  52
  53     if (oldreg != mreg) {
  54 //
  55 //  Double step for boundary crossing
  56 //
  57         fluka->SetTrackIsExiting();
  58         (TVirtualMCApplication::Instance())->Stepping();
  59         fluka->SetMreg(mreg);
  60         fluka->SetTrackIsEntering();
  61         (TVirtualMCApplication::Instance())->Stepping();
  62         fluka->SetTrackIsInside();
  63     } else {
  64         (TVirtualMCApplication::Instance())->Stepping();
  65     }
  66 } // end of mgdraw
  67 } // end of extern "C"
  68