Run

From solidwiki
Revision as of 20:55, 3 December 2020 by Saw (Talk | contribs) (Created 01/20/11(Zwzhao), Last modified 09/05/13(Zwzhao))

Jump to: navigation, search

There are files under "script" at SoLID_GEMC SVN to use.

edit "set_solidGEMC_farm" or "set_solidGEMC_local" according to your system, see instruction in the file

run "source set_solidGEMC_farm" or "source set_solidGEMC_local" to setup environment.

then do ssh tunnel if you are not on jlab network, refer to Solid mysql database

If you are on local machine, put a copy of SoLID CLEO magnetic field map file solenoid_CLEO.dat at where defined by the FIELD_DIR option.

If you are on ifarm, the field map are there already

then run the commands below in a terminal

run_solid_CLEO_SIDIS_3he
run_solid_CLEO_PVDIS

These are examples of gemc running interactive mode. you can change options within to your need.

keep -DBHOST=$GEMC_HOST unchanged to accommodate mysql server access.

change "-USE_PHYSICSL=gemc" to use other physics list, but the current default is "QGSP_BERT". Geant4 group recommend "FTFP_BERT".

To record passby particles which deposit zero energy in the sensitive area, turn it on by "-RECORD_PASSBY=1"

turn on optical physics by "-OPT_PH=1"

use -USE_QT=0 to run in batch mode

use option -geometry=1000x1000 to get larger screen instead of default 600x600

use -EXEC_MACRO=$SoLID_GEMC/script/solid_slice.vis to see the cut-half plot.

combine gcard and command line options, you can control manyfeature of your simulation.

for full command line options, see "Documentation" section at gemc website or use "gemc -help" to get an up2date list, here is a list for your convenience retrieved for GEMC 1.7



    Usage: -Option=<option>

    Options:

   > Option ALIGN_ZAXIS: Align z axis to a custom direction. Options:
      - "beamp"  aligns z axis to the beam directions specified by BEAM_P.
      - "custom, theta*unit, phi*unit" aligns z axis to a custom direction, changes BEAM_P reference frame.
   > Option BANK_DATABASE: Select mysql Bank Database.
   > Option BEAM_P: Beam particle, momentum, angles (in respect of z-axis). 
      Example: -BEAM_P="e-, 6*GeV, 15*deg, 20*deg" sets 6 GeV electrons 15 degrees in theta, 20 degrees in phi. 
      Use -BEAM_P="show_all" to print the list of G4 supported particles.

   > Option BEAM_V: Beam Vertex. Example: -BEAM_V="(0, 0, -20)cm". 
   > Option CATCH: Catch volumes matching the given string.
   > Option CHECK_OVERLAPS: Checks Overlapping Volumes:
      1.  Check Overlaps at Construction Time
      2.  Check Overlaps based on standard lines grid setup
      3.  Check Overlaps by shooting lines according to a cylindrical pattern

   > Option DATABASE: Select mysql Database.
   > Option DAWN_N: Number of events to be displayed with the DAWN driver (also activate the DAWN driver).
   > Option DBHOST: Select mysql server host name.
   > Option DBPSWD: Select mysql password
   > Option DBUSER: Select mysql user name
   > Option DC_MSTAG_R2: Mini Stagger for Region 2. Each layer will alternate +- |this value|
   > Option DC_MSTAG_R3: Mini Stagger for Region 3. Each layer will alternate +- |this value|
   > Option DEFAULT_MATERIAL: Default material for missing material field.

   > Option ENERGY_CUT: Set an energy cut in MeV below which no particle will be tracked further. -1. turns this off.
   > Option EVN: Initial Event Number.
   > Option EXEC_MACRO: Executes commands in macro file.
   > Option FIELD_DIR: Magnetic Field Maps Location. 

   > Option G4P_VERBOSITY: Controls Physical Volumes Construction Log Output.
   > Option G4TRACK_VERBOSITY: Controls Geant4 Track Verbosity.
   > Option GEN_VERBOSITY: Controls Geant4 Generator Verbosity.
   > Option GEO_VERBOSITY: Controls Geometry Construction Log Output.
   > Option GT: Selects Geometry table. This option is overwritten with the gemc read card.
   > Option GUI_VERBOSITY: Controls GUI Construction Log Output.
   > Option HALL_FIELD: Magnetic Field of the Hall. 

   > Option HALL_MATERIAL: Composition of the Experimental Hall. 
            Air normal simulation
            Air_Opt Simulation with Optical Physics (default)
            Vacuum

   > Option HIGH_RES: Use High Resolution Graphics

   > Option HIT_PROCESS_LIST: Registers Hit Process Routines. Can register multiple experiments, separated by space, e.v. "clas12 aprime"
      clas12.  CLAS12 hit process routines (default)
      aprime.  aprime hit process routines
      gluex.   GlueX  hit process routines

   > Option HIT_VERBOSITY: Controls Hits Log Output. 
   > Option INPUT_GEN_FILE: Generator Input. Current availables file formats:
      LUND. 
      xample: -INPUT_GEN_FILE="LUND, input.dat"

   > Option LOG_MSG: Log Messages Header.
   > Option LOW_EM_PHYS: Turn on the low energy Electro-Magnetic physics, down to the X-rat scale.
     Currently only implemented for gemc physics list. 
         0  Turn off [default]  
         1  Standard E&M down to X-rays ( ~1 kEV cutoff) 

   > Option LUMI2_EVENT: Luminosity Beam 2 Parameters: number of Particles/Event, Time Between Bunches. The Time Window is specified with the LUMI_EVENT flag
            Example: -LUMI2_EVENT="10000, 2*ns" simulate 10K particles per event at 2ns intervals. 

   > Option LUMI2_P: Luminosity Beam Particle 2, momentum, angles (in respect of z-axis). 
            Example: -LUMI2_P="proton, 1*GeV, 25*deg, 2*deg" sets 1 GeV protons, 25 degrees in theta, 2 degrees in phi. 
            Use -LUMI2_P="show_all" to print the list of G4 supported particles.

   > Option LUMI2_SPREAD_V: Spread Luminosity Beam 2 Radius, Z position. Example: -SPREAD_V="(0.1, 10)cm". 
   > Option LUMI2_V: Luminosity Beam Particle 2 Vertex. Example: -LUMI2_V="(0, 0, -20)cm". 
   > Option LUMI_EVENT: Luminosity Beam Parameters: number of Particles/Event, Time Window, Time Between Bunches
            Example: -LUMI_EVENT="10000, 120*ns, 2*ns" simulate 10K particles per event distributed over 120 ns, at 2ns intervals. 

   > Option LUMI_P: Luminosity Beam particle, momentum, angles (in respect of z-axis). 
            Example: -LUMI_P="proton, 1*GeV, 25*deg, 2*deg" sets 1 GeV protons, 25 degrees in theta, 2 degrees in phi. 
            Use -LUMI_P="show_all" to print the list of G4 supported particles.

   > Option LUMI_SPREAD_V: Spread Luminosity Beam Radius, Z position. Example: -SPREAD_V="(0.1, 10)cm". 
   > Option LUMI_V: Luminosity Beam Vertex. Example: -LUMI_V="(0, 0, -20)cm". 
   > Option MATERIALSDB: Select Materials DB. 
  Available Databases: 
  * CPP: use the normal geant4 c++ constructor.
  * MYSQL: use the mysql DB.

   > Option MATERIAL_VERBOSITY: Controls Geant4 Material Verbosity.
   > Option MAX_FIELD_STEP: Sets Maximum Acceptable Step in Magnetic Field (in mm).

   > Option MAX_X_POS: Max X Position in millimeters. Beyond this the track will be killed
   > Option MAX_Y_POS: Max Y Position in millimeters. Beyond this the track will be killed
   > Option MAX_Z_POS: Max Z Position in millimeters. Beyond this the track will be killed
   > Option MERGE_FILE: Merge banks from filename. Format supported: 
      - EVIO
   > Option MGN_VERBOSITY: Controls Magnetic Fields Log Output.
   > Option N: Number of events to be simulated.
   > Option NGENP: Max Number of Generated Particles to save in the Output.
   > Option NO_FIELD: Sets Magnetic Field of a volume to zero. "all" means no magnetic field at all. 

   > Option OPT_PH: Activate Optical Physics Processes in gemc Physics List.
   > Option OUTPUT: Type of output, output filename. Supported output: evio, txt. Example: -OUTPUT="evio, out.ev"
   > Option OUT_VERBOSITY: Controls Bank Log Output.
   > Option PARAMETERSDB: Select Parameters DB. 
  Available Databases: 
  * MYSQL: use the mysql DB.

   > Option PARAMETER_VERBOSITY: Controls Parameters Verbosity.
   > Option PHY_VERBOSITY: Controls Physics List Log Output.
   > Option POLAR: Beam particle, polarization percentage and angles  (in respect of z-axis). 
      Example: -POLAR="90, 90*deg, 270*deg" sets 90% polarization 90 degrees in theta, 270 degrees in phi. 
      Use -POLAR="show_all" to print the list of G4 supported particles.

   > Option PRINT_EVENT: -PRINT_EVENT=N: Print Event Number every N events.
   > Option QTSTYLE: Sets the GUI Style. Available options: 
      - QCleanlooksStyle 
      - QMacStyle 
      - QPlastiqueStyle 
      - QWindowsStyle 
      - QMotifStyle
   > Option RANDOM: Random Engine Initialization. The argument (seed) can be an integer or the string TIME.
   > Option RECORD_MIRRORS: Set to one if you want to save mirror hits in the output. Default is 0.

   > Option RECORD_PASSBY: Set to one if you want to save zero energy hits in the output. Default is 0.

   > Option RUNNO: Run Number. Controls the geometry and calibration parameters.

   > Option SAVE_ALL_MOTHERS: Set to 1 to save mother vertex and pid infos in output. High Memory Usage. Default is 0.

   > Option SCALE_FIELD: Scales Magnetic Field by a factor.
      Usage:
      -SCALE_FIELD="fieldname, scalefactor"
      Example: -SCALE_FIELD="srr-solenoid, 0.5"

   > Option SECONDARIES: Threshold (in MeV) to save secondary particles. Negative: don't save any.
   > Option SPREAD_P: Spread beam energy and angles (in respect of z-axis). 
      Example: -SPREAD_P="0*GeV, 10*deg, 20*deg" spreads 10 degrees in theta, 20 degrees in phi. 

   > Option SPREAD_V: Spread Beam Radius, Z position. Example: -SPREAD_V="(0.1, 10)cm". 
   > Option USE_PHYSICSL: Physics List. Avaliable choices: 

            * gemc: comprehensive physics list. Optical Physics may be activated with OPT_PH=1 

            The following is a list of other physics lists. More infos can be found here:

            http://geant4.cern.ch/support/proc_mod_catalog/physics_lists/referencePL.shtml

            * LHEP: This is the main LHEP based physics list, using exclusively parameterised modeling. 
            * LHEP_BERT: Like LHEP, but using Geant4 Bertini cascade for primary protons, neutrons, 
               pions and Kaons below ~10GeV.  
            * LHEP_BERT_HP: Like LHEP_BERT with the addition to use the data driven high precision.
               neutron package (NeutronHP) to transport neutrons below 20 MeV down to thermal energies.
            * QGSP:  Quark-Gluon String model based physics list. 
            * QGSP_BERT: Like QGSP, but using Geant4 Bertini cascade for primary protons, neutrons, 
               pions and Kaons below ~10GeV.  
            * QGSP_BERT_HP: Like QGSP_BERT with the addition to use the data driven high precision.
               neutron package (NeutronHP) to transport neutrons below 20 MeV down to thermal energies.
            * QGSP_BIC: Like QGSP, but using Geant4 Binary cascade for primary protons and neutrons 
               with energies below ~10GeV.  
            * QGSP_BIC_HP: Like QGSP_BIC with the addition to use the data driven high precision.
               neutron package (NeutronHP) to transport neutrons below 20 MeV down to thermal energies.
            * QGSC_BERT: The quark-gluon string (QGS) part handles the formation of strings in the initial
              collision of a hadron with a nucleon in the nucleus. String fragmentation into hadrons is
              handled by the Quark-Gluon String fragmentation model. The Chiral Invariant Phase Space (CHIPS)
              part handles the de-excitation of the remnant nucleus. Uses Geant4 Bertini cascade for primary
              protons, neutrons, pions and Kaons below ~10GeV. 

   > Option USE_QT: QT GUI switch
      0.  Don't use the graphical interface
      1.  QT OpenGL Immediate mode (can interact with picture; sliders works well; slower than Stored mode)
      2.  OpenGL Stored mode (can't interact with picture; sliders works well)
      3.  OpenGL Immediate mode (can't interact with picture; sliders works well; slower than Stored mode)
      4.  QT OpenGL Stored mode (can interact with picture; sliders works but picture needs to be updated by clicking on it)

   > Option gcard: gemc card file.
   > Option geometry: Specify the size of the QT display window. Default '600x600'