User Tools

Site Tools


community:hepsim:usage_full
This manual is outdated. Please use http://atlaswww.hep.anl.gov/hepsim/doc/doku.php

<< back to HepSim manual

Full detector simulations

Currently, the Geant4-based detector simulation can be done using the SLIC software. You will need to convert ProMC files to LCIO files and use these files for the “slic” Geant4-based program. To process SLCIO files, you will need Jas4pp program.

The files with full simulations are located under the “rfullNNN” tags (“NNN” is a number) Here how you can find all samples that have “rfull” tag for download. Click here: http://atlaswww.hep.anl.gov/hepsim/list.php?find=rfull. The description of the detector tags is given in Detector description. If you need to find all files for a give tag, say rfull005, use this command:

hs-find tev10mumu%rfull005

which searches for all samples with the dataset name “tev10mumu” (mu+mu-, 40 TeV samples). Then download a selected sample as:

hs-get tev10mumu_pythia6_zprime10tev_qqbar%rfull005 

The detector repository contains the description of available detectors. Each “rfullNNN” tag corresponds to the detector given in the detector repository.

The physics performance studies based on full simulation files from HepSim are listed here.

Processing LCIO files

To perform data analysis using LCIO files, use the Jas4pp data-analysis package. The installation does not have external dependencies besides Java. Make sure Java7 and above is installed. Then run these commands to install the package using Linux/Mac with the “bash” shell:

wget http://atlaswww.hep.anl.gov/asc/jas4pp/download/current.php -O jas4pp.tgz
tar -zvxf jas4pp.tgz
cd jas4pp
source ./setup.sh # takes 5 sec for first-time optimization

The last command optimizes Java packages for the Python language. Now you are ready to run over any file with truth-level and datector-simulation files, such as LCIO and ProMC. This program can be used for:

  • Downloading and searching HepSim data
  • Processing ProMC files from HepSim
  • Running over SLCIO files with Geant4 simulated / reconstructed events.
  • Data analysis (jets, physics vectors, histogram packages)
  • Visualisation of reconstructed events using Wired4 display

You can find more details in Jas4pp manual. The examples showing how to run over SLCIO files using Jas4pp are located in the “examples/slic” directory of the Jas4pp package.

The standard set of containers is available from the LCIO files posted on HepSim. Read Manual. The API is described in LCIO Java API description.

Look the Java API definition of the needed objects in org.lcsim.lcio. The data containers are shown below:

Click to display ⇲

Click to hide ⇱

BeamCalHits
CalorimeterHitRelations
EM_BARREL
EM_ENDCAP
EcalBarrelHits
EcalEndcapHits
HAD_BARREL
HAD_ENDCAP
HcalBarrelHits
HcalEndcapHits
HelicalTrackHitRelations
HelicalTrackHits
HelicalTrackMCRelations
LumiCalHits
MCParticle
MUON_BARREL
MUON_ENDCAP
MuonBarrelHits
MuonEndcapHits
PandoraPFOCollection
ReconClusters
SiTrackerBarrelHits
SiTrackerEndcapHits
SiTrackerForwardHits
SiVertexBarrelHits
SiVertexEndcapHits
StateAtECal
StateAtEnd
StateAtStart
TKR_RawTrackerHits
TKR_TrackerHits
Tracks
VXD_RawTrackerHits
VXD_TrackerHits

Using C++/ROOT with LCIO

If you need to read LCIO files in C++ code with ROOT/FastJet, use this example package:

wget http://atlaswww.hep.anl.gov/hepsim/soft/lcio-cpp.tgz -O - | tar -xz;
cd lcio-cpp; make 

It will compile “example.cpp” file linking FastJet package. Look at README file to see how to install the “LCIO” package required by this example. Your LCIO“ files from HepSim should be in the directory “data”.

Visualizing events with Jas+Wired

You can run Jas+Wired to visualize the simulated events. The Wired program is included in Jas4pp, so you simply run it as:

./jaspp

This will start a Jas3-like environment with all needed plugins. Then copy the detector geometry file from the detector repository.

wget http://atlaswww.hep.anl.gov/hepsim/soft/detectors/sidloi3.tgz -O - | tar -xz;

This will create a directory “sidloi3”. This detector corresponds to “rfull001” tag used for the reconstruction of pythia6_zpole_ee (Z→e+e-).

Now we can visualize the detector as [File]-[Open data source]-[HepRep] XML and select the file “sidloi3.heprep” from sidloi3.tgz: This is how to do this using the command line:

./jaspp  sidloi3.heprep

You will see the detector layout:

Whired4 event display

Now, we will read the event: Open any *.slcio file you copied from HepSim as [File]-[Open data source]-[LCIO] file. Then click a small button [Go] (top menu bar). It will process events. Then select again [File]-[New]-[Wired 4 view]. You will get an image in the Wired4 display. Now press [Go] again to look at next event.

If you want to see how data records are organized inside the slcio file, do this [File]-[New]-[LCSim Event browser]

Generally, when you open a SLCIO file from the HepSim, you do not need to load the geometry file. The geometry file will be downloaded automatically from HepSim and will be put to ”.lcsim/cache/“ inside your home directory.

Visualizing geometry using ROOT

(contribution from N.Nikiforou):

You can also visualize the detector using ROOT. Here are a few steps: Locate the file detector.gdml inside the file detector.zip (where “detector” shows the name of the detector).

show more details here

show more details here

This file is created using this command during the detector design stage

slic -g detector.lcdd -G detector.gdml

Once you have the GDML file, you can use ROOT to visualize it. You just need to make sure the ROOT installation has openGL/GDML support (CERN AFS installations have it for sure). At the ROOT prompt do:

TGeoManager::Import("detector.gdml");
gGeoManager->GetTopVolume()->Draw("ogl");

This should popup an OpenGL display with the detector which you can clip, pan, rotate etc. You will the image as shown here:

People

  • S.Chekanov (ANL, main developer)
  • E.May (ANL, debugging, Wired4)

Additional help

community/hepsim/usage_full.txt · Last modified: 2017/09/30 14:35 by asc