HepSim

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
{{indexmenu_n>6}}

[[:|<< back to HepSim manual]]

====== Detector visualization ======

Much of this topic has been discussed in the Section [[hepsim:jas4pp|Jas4pp Manual]]. See also topics discussed in Sections [[fcs:sid:intro|SiD introduction]] and [[fcs:fcchh:intro|FCC introduction]].



====== Event 1 ======

Here is one example of event visualization: A hypothetical  particle (Z') with the mass 40 TeV decays to W+ and W-  bosons. 
Such heavy particles can only be produced at the FCC-hh
experiment with the center-of-mass of energy of 100 TeV for colliding protons.
The transverse momenta of these bosons from the Z' decay is about 20 TeV. 
Here are a few images illustrating a single event Z'(40 TeV) → W+W- → hadrons:

{{:hepsim:z40_qqbar1.png?310|}}
{{:hepsim:z40_qqbar2.png?320|}}
{{:hepsim:z40_qqbar4.png?275|}}

{{:hepsim:z40_qqbar6.png?310|}}
{{:hepsim:z40_qqbar8.png?310|}}
{{:hepsim:z40_qqbar9.png?300|}}


This simulation was created using a FCC-like detector (SiFCC) developed at ANL/Fermilab/Duke. 
This detector is described [[http://atlaswww.hep.anl.gov/hepsim/soft/detectors/sifcch4/sifcch4.html|here]]. The Monte Carlo sample used
for this simulation  is available from the [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=206|HepSim repository]].

You can see two high-pT jets (~20 TeV) from each W decay, and many back-splash particles.
This event has many secondary interaction in the detector material. These images show  Geant4 truth-level particles (straight lines) created by secondary interactions
due to "back-splash" events.

{{:hepsim:z40_sifcc4_v3.png?300|}}
{{:hepsim:z40_sifcc4_v7.png?310|}}

A real pp (100 TeV) → Z'(40)  → W+W- is expected to be more complicated due to 
the presence of the protons in the initial state and pileup events (up to 1000 MinBias events).



This simulation is very challenging since
the [[http://atlaswww.hep.anl.gov/hepsim/soft/detectors/sifcch4/sifcch4.html|SiFCC(v4) detector]] has high-granularity calorimeters with more than 150 million cells.  For example, the HCAL (TileCal) of the ATLAS experiment has about 5000 cells (3 layers), while SiFCC has 50 million cells (64 layers).  The simulation of a single event like this in the [[http://atlaswww.hep.anl.gov/hepsim/soft/detectors/sifcch4/sifcch4.html|SiFCC(v4) detector]]
takes 4 hours (8 CPU cores, Xeon(R) CPU E5-2450 2.50GHz) with 16 GB memory on the ANL/HEP cluster. 
The simulation of a typical pp-collision event at the LHC takes 15 min.

The event shown above has:

  * 10,000 hits in ECAL
  * 46,000 hits in HCAL
  * 12,000 hits in the outer tracker
  * 1,000  hits in the inner tracker


====== Event 2 ======

This is another Z'(40 TeV) → W+ W- example, where W+ decays to hadrons, and W- decays to electron and neutrino.  
The neutrino is not detected, but the electron can be well reconstructed as shown here:


{{:hepsim:z40_siffcc_n1.png?280|}}
{{:hepsim:z40_sifcc_n2.png?300|}}
{{:hepsim:z40_sifcc_n3.png?350|}}

This event has 25,000 hits in the HCAL+ECAL, and about 6,000 in the tracker.

====== How to make such images ======


Use [[hepsim:jas4pp|Jas4pp]] as:

<code bash>
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
wget  http://mc.hep.anl.gov/asc/hepsim/events/ee/250gev/pythia6_zpole_bbar/rfull002/gev250_pythia6_zpole_bbar_0001_pandora.slcio
jaspp gev250_pythia6_zpole_bbar_0001_pandora.slcio
</code>
Then press the small "run" (blue arrow) button on the toolbar. Then go to File → New → Wired 4 View. The you should wait for 5-10 before the image will appear. The first event from this file is what you see in the above images.
You will need a computer with at least 8 GB RAM. This example is shown for e+e- events at 250 GeV. Similarly, use other files with the tag "rfull" (full simulations).

You can download all reconstructed files as 

<code>
hs-get tev40mumu_pythia6_zprime_ww%rfull006
</code>


====== Read more ======

Physics studies used this detector were discussed in the talk [[https://indico.cern.ch/event/438866/contributions/1085149/attachments/1256884/1855819/FCCweek_2016_Rome_April.pdf|Simulation of a high-granular hadronic calorimeter for multi-TeV physics]].



 --- //[[chekanov@anl.gov|Sergei Chekanov]] 2016/04/29 16:23//