Run
From solidwiki
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 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'
