# define fluka_closeinp fluka_closeinp_
# define mcihad mcihad_
# define mpdgha mpdgha_
+# define newplo newplo_
#else
# define flukam FLUKAM
# define fluka_openinp FLUKA_OPENINP
// Prepare input file with the current physics settings
InitPhysics();
-
- cout << "\t* InitPhysics() - Prepare input file was called" << endl;
-
- if (fVerbosityLevel >=2)
- cout << "\t* Changing lfdrtr = (" << (GLOBAL.lfdrtr?'T':'F')
- << ") in fluka..." << endl;
- GLOBAL.lfdrtr = true;
-
- if (fVerbosityLevel >=2)
- cout << "\t* Opening file " << fInputFileName << endl;
+// Open fortran files
const char* fname = fInputFileName;
-
fluka_openinp(lunin, PASSCHARA(fname));
fluka_openout(11, PASSCHARA("fluka.out"));
-
- if (fVerbosityLevel >=2)
- cout << "\t* Calling flukam..." << endl;
+// Read input cards
+ GLOBAL.lfdrtr = true;
flukam(1);
-
- if (fVerbosityLevel >=2)
- cout << "\t* Closing file " << fInputFileName << endl;
+// Close input file
fluka_closeinp(lunin);
-
+// Finish geometry
FinishGeometry();
-
- if (fVerbosityLevel >=3)
- cout << "<== TFluka::Init() called." << endl;
-
- if (fVerbosityLevel >=3)
- cout << "<== TFluka::BuildPhysics() called." << endl;
}
//______________________________________________________________________________
}
mopo->WriteFlukaInputCards();
}
-
+
+// Add RANDOMIZ card
+ fprintf(pFlukaVmcInp,"RANDOMIZ %10.1f%10.0f\n", 1., Float_t(gRandom->GetSeed()));
// Add START and STOP card
fprintf(pFlukaVmcInp,"START %10.1f\n",fEventsPerRun);
fprintf(pFlukaVmcInp,"STOP \n");
icode == kEMFSCOcompton || // Compton scattering
icode == kEMFSCOphotoel || // Photoelectric effect
icode == kKASNEUhadronic || // hadronic interaction
- icode == kKASHEAdray // delta-ray
+ icode == kKASHEAdray // delta-ray
) return 1;
else return 0;
}
//______________________________________________________________________________
void TFluka::Gmtod(Double_t* xm, Double_t* xd, Int_t iflag)
{
+//
+// See Gmtod(Float_t*, Float_t*, Int_t)
+//
if (iflag == 1) gGeoManager->MasterToLocal(xm,xd);
else gGeoManager->MasterToLocalVect(xm,xd);
}
//______________________________________________________________________________
void TFluka::Gdtom(Double_t* xd, Double_t* xm, Int_t iflag)
{
+//
+// See Gdtom(Float_t*, Float_t*, Int_t)
+//
if (iflag == 1) gGeoManager->LocalToMaster(xd,xm);
else gGeoManager->LocalToMasterVect(xd,xm);
}
//______________________________________________________________________________
TObjArray *TFluka::GetFlukaMaterials()
{
+//
+// Get array of Fluka materials
return fGeom->GetMatList();
}