This wiki describe track resolution studies using 2 designs with different outer tracker radius. The plots are done assuming flat distribution of single muons in |y|<3.

Here is a short description on how to check track resolution using CEPC-sidcc2 detector which has the SiD envelope for tracking (R=1.25 m), 25 um pitch size and 3 Tesla field. The data are taken from rfull101 dataset list.

Here are the steps to check the resolution using single muons. You will need this script cepc_track_resolution.py and the input dataset pgun_muon_eta3_pt150gev

bash   #  set bash if you haven't done this before, and make sure Java is installed. 
wget http://atlaswww.hep.anl.gov/asc/jas4pp/download/current.php -O jas4pp.tgz
tar -zvxf jas4pp.tgz
source jas4pp/setup.sh # takes 5 sec for first-time optimization
hs-get pgun_muon_eta3_pt150gev%rfull101 data 2 500 # get 500 files in 2 threads
# download the python script shown above and run it as
fpad cepc_track_resolution.py

This script calculates the tracking resolution of this detector as shown below:

Here is the script that does sigma(1/pt). cepc_track_resolution_overpt.py

A new CEPC-Si detector was created using B=3 and R=1.8m tracker, which can fits to the envelope of CEPC-ILD detector, The new detector, called CEPC-sidcc3, with modified tracker was used to create rfull102 dataset list.

The analysis of this data can be done exactly as before. The only difference is that you will need to download the different dataset:

hs-get pgun_muon_eta3_pt150gev%rfull102 data 2 500 # get 500 files in 2 threads

The track resolutions are shown in

Here is the tracking resolution for 1/pt:

The data file is cepc_track_resolution_overpt.jdat. You can read jdat file in PyROOT using this simple macro script to read jdat files

As before, this was done using B=3 and R=1.8m tracker and CEPC-sidcc3.

Data files are cepc_track_resolution_overpt_barrel.jdat, cepc_track_resolution_overpt_endcap.jdat

The results were created with the scripts cepc_track_resolution_overpt_barrel.py, cepc_track_resolution_overpt_endcap.py

The files with data are:

• cepc_track_resolution_z0_barrel.jdat
• cepc_track_resolution_z0_endcap.jdat
• cepc_track_resolution_d0_barrel.jdat
• cepc_track_resolution_d0_endcap.jdat

Here are a few plots illustrating physics performances. The figures have been done using fully identified muons and electrons (PandoraPFO). The following samples (tag rfull102) have been used:

• gev250ee_pythia6_zpole_mumu - Z→mu+mu-
• gev250ee_pythia8_higgs_mumu - Higgs →mu+mu+ (all Higgs processes)
• gev250ee_pythia6_zpole_ee - Z→ e+e-
• gev250ee_pythia8_zhiggs_nunumumu - Z0+Higgs→ nunu+mu+mu-
• gev250ee_pythia8_zhiggs_nunugg - Z0+Higgs→ nunu+gg

You can get the data (for Z0+Higgs→ nunu+mu+mu-, for example), as:

bash   #  set bash if you haven't done this before 
wget http://atlaswww.hep.anl.gov/hepsim/soft/hs-toolkit.tgz -O - | tar -xz;
source hs-toolkit/setup.sh
hs-get gev250ee_pythia8_zhiggs_nunumumu          # download ProMC (truth) files
hs-get gev250ee_pythia8_zhiggs_nunumumu%rfull102 # download SLCIO files for rfull102 tag

See ProMC manual.

Now you can convert ProMC data to STDHEP or LCIO (MCParticle) with full simulations using promc2stdhep or promc2lcio converters as discussed in ProMC to LCIO/STDHEP section.

Here are a few plots that show the performance using gev250ee_pythia6_zpole_mumu:

Here is the plot for gev250ee_pythia8_higgs_mumu - Higgs →mu+mu+ (all Higgs processes)

Here is the plot for gev250ee_pythia6_zpole_ee - Z→ e+e-

— Sergei Chekanov 2017/09/07 17:04