Search/Navigation:
Related:
SLAC/EPP/HPS Public
Jefferson Lab/Hall B/HPS Run Wiki
S30XL-LESA/LDMX
This guide will explain how to generate physics events in an HPS setup including the ECal and the tracking detectors.
The Maven 2 tool must be installed locally (not covered here), as well as a Java JDK.
You must first set CVSROOT to SLAC CVS (example in bash shell).
No Format |
---|
export CVSROOT=:pserver:anonymous@cvs.freehep.org:/cvs/lcd
|
Other shells like csh have not been tested.
Create a work area where packages will be installed.
No Format |
---|
mkdir hps
cd hps
|
In order to get a working version of slic, follow the
SimDist Build Instructions.
SLAC users with NFS access may use this script.
No Format |
---|
/nfs/slac/g/lcd/mc/prj/sw/dist/SimDist/pro/scripts/slic.sh
|
In your work area, checkout the hps-detectors module from cvs and install it.
No Format |
---|
cvs co hps-detectors
cd hps-detectors
mvn
cd ..
|
You will need:
You need a copy of SLIC.
See the SLIC installation instructions.
The lcdd detector description files are included in the hps-java distribution.
No Format |
---|
hps-java/detector-data/detectors/
|
This will locally install the HPS detector data.
The following example shows how to generate some LCIO events using slic .with an stdhep file as input
No Format |
---|
./SimDist/scripts/slic.sh -g ./hps-detectorsjava/detector-data/detectors/HPS-EcalTest-v1pt0/HPS-EcalTest-v1pt0.lcdd \ -i /nfs/slac/g/hps/mgraham/DarkPhoton/SignalEvents/ap6.6gev100mevsel_10_20u_beamspot_gammactau_0cm.stdhep \ -x -O -o hpsEcalTest -r 10000 |
If you do not have SLAC NFS access, sample signal events can be obtained from this ftp space.
...
Note that a full listing of the command line options for slic can be printed out with the following command:
No Format |
---|
> slic -h
Wed Oct 16 11:53:20 2019 :: OKAY :: SlicApplication :: SLIC is starting.
SLIC
**************
* SLIC Usage *
**************
Command line usage:
slic [single_macro_path]
slic [options]
Interactive usage:
slic -n
************************
* Command Line Options *
************************
Option Full Name Macro Command Description
------------------------------------------------------------
-h --help /slic/usage Print SLIC usage.
-? --help /slic/usage Print SLIC usage.
-n --interactive /control/interactive Start a Geant4 interactive session.
-v --version /slic/version Print SLIC version.
-m --macro /control/execute Execute Geant4 commands from a file.
-g --lcdd-url /lcdd/url Set LCDD geometry file URL.
-i --event-file /generator/filename Set event input file full path.
-o --lcio-file /lcio/filename Set name of LCIO output file.
-p --lcio-path /lcio/path Set directory for LCIO output.
-O --autoname /lcio/autoname Automatically name the LCIO output file.
-x --lcio-delete /lcio/fileExists delete Delete an existing LCIO file.
-r --run-events /run/beamOn Run # of events.
-s --skip-events /generator/skipEvents Set number of events to skip.
-l --physics-list /physics/select Set Geant4 physics list.
-L --log-file /log/filename Set logfile name.
-d --seed /random/seed Set the random seed. (No argument seeds with time.)
-G --dump-gdml /lcdd/dumpGDML Dump geometry to GDML file.
-c --optical /physics/enableOptical Enable optical physics processes.
-P --pdg-file /physics/setPDGFile Set location of particle.tbl file.
-V --lcdd-version /lcdd/version Set the version of the GDML setup tag.
-S --lcdd-setup /lcdd/setupName Set the name of the GDML setup tag.
|
The use of the -O argument is highly recommended as it includes some run metadata into the output file name, which can be very useful in providing some provenance for the content.
One can also generate events internally using the GEANT4 General Particle Source. One can either issue commands on the interactive command line, or invoke them via a Geant4 macro.
No Format |
---|
slic -g detector.lcdd -m G4commands.mac -O
|
Here is a simple example to generate 100k single 4.56 GeV electrons isotropically over a fixed range in theta and uniformly in phi:
No Format |
---|
/generator/select gps
/gps/particle e-
/gps/position 0 0 0 cm
/gps/energy 4.56 GeV
/gps/ang/type iso
/gps/ang/mintheta 170 deg
/gps/ang/maxtheta 180 deg
/gps/ang/minphi 0 deg
/gps/ang/maxphi 360 deg
/run/beamOn 100000 |