ATLAS Support Center

This wiki explains how to view electron-proton collision events for CEPC. To analyse events for such studies, install Jas4pp program (assumes “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

Description is available from https://wiki.inside.anl.gov/inside/HEP/CEPC.

Let's calculate Z peak from particle-flow objects after full reconstruction using Pandora. You do not run “hs-get” command if you have done this before.

cd examples/slic/
hs-get gev250ee_pythia6_zpole_ee%rfull002 gev250ee_pythia6_zpole_ee # download all files
fpad mc_pflow.py

you can also do do some simple fit:

fpad mc_pflow_fit.py

Let's calculate Z peak from reconstructed tracks (no Pandora and no particle identification) You do not run “hs-get” command if you have done this before.

fpad mc_tracks.py

Please correct on line 47 the variable Bfield=5 to Bfield=4 Tesla!

Let's calculate the Higgs peak from reconstructed photons after the PFA algorithm. This time you will need H→ gamma+gamma sample gev250ee_pythia6_higgs_gamgam: We set pT(γ)>10 GeV.

hs-get gev250ee_pythia6_higgs_gamgam%rfull002 gev250ee_pythia6_higgs_gamgam 
fpad mc_gammagamma.py

Let's calculate Z peak from reconstructed b-jets. Jets are reconstructed using the Durham kT algorithm with ycut=0.05. This time we will download Z→ bbar sample. Look at gev250ee_pythia6_zpole_bbar. We assume pT(jet)>20 GeV.

hs-get gev250ee_pythia6_zpole_bbar%rfull002 gev250ee_pythia6_zpole_bbar
fpad mc_jets.py

No any optimization for jet reconstruction was used. Consider other approaches to make more realistic dijet mass.

Let's calculate Higgs peak from reconstructed leptons in the channel H(125)→ ZZ* → 4L (electrons and muons). We will use particle flow objects to reconstruct and identify electrons and muons. This time we will download Higgs→ 4 lepton sample. Look at gev250ee_pythia6_higgs_zz_4l.

hs-get gev250ee_pythia6_higgs_zz_4l%rfull002 gev250ee_pythia6_higgs_zz_4l  # download 100 files in 2 threads
fpad mc_higgs_zz.py

This example uses the Jade algorithm to force 2 jets, and uses H→ bbar sample:

hs-get gev250ee_pythia6_higgs_bbar%rfull002 gev250ee_pythia6_higgs_bbar  
fpad mc_higgs_bbar.py

This example shows the jet energy resolution in 2 pT(truth) regions. It is based on inclusive jet sample gev250ee_pythia6_qcd_all. We use Durham jets with pT(truth)>20 GeV. Jet are constructed from PFA objects. Truth jet and reco jets are matched using dR<0.05 in phi and rapidity.

hs-get gev250ee_pythia6_qcd_all%rfull002 gev250ee_pythia6_qcd_all
fpad mc_jets_resol.py

— Sergei Chekanov 2015/10/08 09:31