]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TFluka/TFluka.cxx
Remove warning message.
[u/mrichter/AliRoot.git] / TFluka / TFluka.cxx
index e3e4be18924a1aaeace909eb1413d905380af7b3..6acbbfc5032317a007621f7befc0e9036dc626c7 100644 (file)
@@ -71,6 +71,7 @@
 # define fluka_closeinp fluka_closeinp_
 # define mcihad mcihad_
 # define mpdgha mpdgha_
+# define newplo newplo_
 #else 
 # define flukam  FLUKAM
 # define fluka_openinp FLUKA_OPENINP
@@ -265,36 +266,17 @@ void TFluka::BuildPhysics() {
     // 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;
 }  
 
 //______________________________________________________________________________ 
@@ -1182,7 +1164,9 @@ void TFluka::InitPhysics()
        }
        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");
@@ -1634,7 +1618,7 @@ Bool_t   TFluka::IsTrackDisappeared() const
       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;
 }
@@ -1958,6 +1942,9 @@ void TFluka::Gmtod(Float_t* xm, Float_t* xd, Int_t iflag)
 //______________________________________________________________________________ 
 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);
 }
@@ -1992,6 +1979,9 @@ void TFluka::Gdtom(Float_t* xd, Float_t* xm, Int_t iflag)
 //______________________________________________________________________________ 
 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);
 }
@@ -1999,6 +1989,8 @@ void TFluka::Gdtom(Double_t* xd, Double_t* xm, Int_t iflag)
 //______________________________________________________________________________
 TObjArray *TFluka::GetFlukaMaterials()
 {
+//
+// Get array of Fluka materials
    return fGeom->GetMatList();
 }