gdml2G4viz.cc

Go to the documentation of this file.

/*
* 
* Copyright: Mathieu Trocme, mathieu.trocme@gmail.com (01 october 2009)
* 
* This software is a computer program whose purpose is to: convert any GDML file to any G4 visualization driver file. 
* 
* This software is governed by the CeCILL license under French law and abiding by the rules of distribution of free software. 
* You can use, modify and/or redistribute the software under the terms of the CeCILL license as circulated by CEA, CNRS and INRIA at the following URL "http://www.cecill.info".
* As a counterpart to the access to the source code and rights to copy, modify and redistribute granted by the license, users are provided only with a limited warranty and the software's author, the holder of the economic rights, and the successive licensors have only limited liability. 
* In this respect, the user's attention is drawn to the risks associated with loading, using, modifying and/or developing or reproducing the software by the user in light of its specific status of free software, that may mean that it is complicated to manipulate, and that also therefore means that it is reserved for developers and experienced professionals having in-depth computer knowledge. 
* Users are therefore encouraged to load and test the software's suitability as regards their requirements in conditions enabling the security of their systems and/or data to be ensured and, more generally, to use and operate it in the same conditions as regards security. 
* 
* The fact that you are presently reading this means that you have had knowledge of the CeCILL license and that you accept its terms.
* 
*/

// 
// Program       : gdml2G4viz.cc
// Task          : Convert any GDML file to any G4 visualization driver file
// Location      : $G4MYSIMU (See https://sbgtrac.in2p3.fr/projects/ramses/wiki/G4_MyTuto for the list of environment variables used)
// Compilation   : 'make' ou 'make -j $NPROC'  # export NPROC=$(cat /proc/cpuinfo | grep "processor.*:" | wc -l)
// Execution     : 'cd $G4BINDIR ; [ ! -a gdml2G4viz.mac ] && ln -s $G4MYSIMU/gdml2G4viz.mac ; ./gdml2G4viz myGDMLfile.gdml myMacroFile.mac'
// Output        : It's up to you. Default will be: $G4BINDIR/myGDMLfile.wrl
// Documentation : See $G4MYSIMU/doc.html
// Convention    : Root coding convention: TMyClass for classes, fMyVariable for class data members, kMyVariable for constants, gMyVariable for global variables
// Help          : http://www.lcsim.org/software/geant4/doxygen/html/classes.html
// Author        : Mathieu Trocme (mathieu.trocme@ires.in2p3.fr)
// Date          : 29 September 2009
// Framework     : PhD thesis, CNRS/IPHC/DRS/RAMSES, Strasbourg, France
// 


#include "G4RunManager.hh"
#include "G4VisExecutive.hh"
#include "G4UImanager.hh"

#include "G4GDMLParser.hh"
#include "TGDML_DetectorConstruction.hh"

#include "TGPS_PrimaryGeneratorAction.hh"  /// cannot be skipped unless someday... http://hypernews.slac.stanford.edu/HyperNews/geant4/get/eventtrackmanage/895.html
#include "QGSP.hh"                         /// or any other physics list (See $G4SOURCE/physics_lists/lists/src/G4PhysListFactory.cc). This is just a fake for G4.



G4String RenameG4VizFile ( G4String gdmlFileName, G4String g4VizDefaultFileName )
{
   G4String extension = g4VizDefaultFileName ; extension.replace(0,extension.first('.')+1,"") ;
   G4String g4VizFileName = gdmlFileName ; g4VizFileName.replace(g4VizFileName.last('.')+1,4,extension.data()) ;
   system (G4String("mv "+g4VizDefaultFileName+" "+g4VizFileName).data()) ; 
   return (g4VizFileName) ;
}


/////


int main (int argc, char **argv)
{
   /// Parsing the input line

   G4String macroFileName="", gdmlFileName="" ; FILE* fileTest = NULL ;
   if ( argc != 3 ) 
      { G4Exception("\n!!!\n!!! Error in main(): USAGE is './gdml2G4viz myMacroFile.mac myGDMLfile.gdml'.\n!!!") ; }
   else
   {
      gdmlFileName = G4String(argv[2]) ; macroFileName = G4String(argv[1]) ; 
      /// checking the parameter order
      if ( !gdmlFileName.contains(".gdml") ) { G4Exception("\n!!!\n!!! Error in main(): USAGE is './gdml2G4viz myMacroFile.mac myGDMLfile.gdml'. Mind the order!\n!!!") ; }
      /// checking if the gdml file exists
      fileTest = fopen (gdmlFileName.data(),"r") ;  
         if ( fileTest == NULL ) { G4Exception( G4String("\n!!!\n!!! Error in main(): The gdml file '" + gdmlFileName + "' cannot be opened or read. It may not exist. \n!!!").data() ) ; }  
      fclose (fileTest) ;
      /// checking if the macro file exists
      fileTest = fopen (macroFileName.data(),"r") ;  
         if ( fileTest == NULL ) { G4Exception( G4String("\n!!!\n!!! Error in main(): The macro file '" + macroFileName + "' cannot be opened or read. It may not exist. \n!!!").data() ) ; }  
      fclose (fileTest) ;
   }

   /// Reading the GDML file

   G4GDMLParser parser ;
   parser.Read (gdmlFileName.data()) ; 

   /// Run manager

   G4RunManager* runManager = new G4RunManager() ;
   runManager->SetUserInitialization ( new TGDML_DetectorConstruction(parser.GetWorldVolume()) ) ;
   runManager->SetUserInitialization ( new QGSP() ) ; 
   runManager->SetUserAction ( new TGPS_PrimaryGeneratorAction() ) ;
   runManager->Initialize() ;

   /// Visualization manager

   G4VisManager *visManager = new G4VisExecutive () ;
   visManager->Initialize () ;

   /// UI manager 

   G4UImanager* UI = G4UImanager::GetUIpointer() ;  
   UI->ApplyCommand("/control/execute gdml2G4viz.mac") ; 

   /// Cleaning

   delete visManager ; 
   delete runManager ;  

   /// Renaming the visualization files (if default name is used) with the gdml file root name (e.g. foo.gdml -> foo.heprep.zip, foo.wrl, ...)

   G4String g4VizDefaultFileName="", g4VizFileName="" ;
   g4VizDefaultFileName="gdml2G4viz.heprep.zip" ; if ( system (G4String("[ -a "+g4VizDefaultFileName+" ]").data()) == 0 ) { g4VizFileName = RenameG4VizFile(gdmlFileName,g4VizDefaultFileName) ; }  /// WIRED4, See my post http://forum.freehep.org/index.php?t=msg&th=936&start=0&rid=0&S=7a7fcf2ac3b842d246213ec1f395643b
   g4VizDefaultFileName="G4Data0.heprep"        ; if ( system (G4String("[ -a "+g4VizDefaultFileName+" ]").data()) == 0 ) { system (G4String("[ -a HepRApp.jar ] && java -jar HepRApp.jar -file "+RenameG4VizFile(gdmlFileName,g4VizDefaultFileName)+" &").data()) ; }  /// WIRED3, COPY HepRApp.jar in $G4BINDIR
   g4VizDefaultFileName="g4_00.wrl"             ; if ( system (G4String("[ -a "+g4VizDefaultFileName+" ]").data()) == 0 ) { system (G4String("[ -a firefox ]     && ./firefox "                  +RenameG4VizFile(gdmlFileName,g4VizDefaultFileName)+" &").data()) ; }  /// VRML  , link the firefox executable in $G4BINDIR
   g4VizDefaultFileName="g4_00.prim"            ; if ( system (G4String("[ -a "+g4VizDefaultFileName+" ]").data()) == 0 ) { system (G4String("[ -a dawn ]        && ./dawn "                     +RenameG4VizFile(gdmlFileName,g4VizDefaultFileName)+" &").data()) ; }  /// DAWN  , link the dawn executable in $G4BINDIR

   return (0) ;  
}

Copyright: Mathieu Trocmé, Sun 4 Oct 2009 (CeCILL license)