Muon system simulation

From Charm-Tau Detector
(Difference between revisions)
Jump to: navigation, search
(Created page with "* This page is dedicated to the sketch standalone simulation of the SCT detector based on pure Geant4")
 
m (→‎Version 0.2)
 
(15 intermediate revisions by one user not shown)
Line 1: Line 1:
−
* This page is dedicated to the sketch standalone simulation of the SCT detector based on pure Geant4
+
== This page is dedicated to the sketch standalone simulation of the SCT detector based on pure Geant4 ==
 +
 
 +
One of the most important and urgent tasks for muon system technology choice and further optimization is to create a fast and reliable simulation tool.
 +
Optimal toolkit in this case is pure Geant4 simulation with simplified geometry and physics lists, able to produce fast results and estimate main parameters and characteristics of the detector.
 +
 
 +
== Version 0.1 ==
 +
=== Tasks for the 0.1 version of the simulation tools ===
 +
 
 +
* Create geant4-based geometry description reflecting all main features influencing muon-related physics. The model are to be extendible to describe muon detector in detail if needed.
 +
* Estimate basic features of the interacting muons (despite of the detector technology choice):
 +
** Muon smearing due to the multiple scattering (depending of the muon energy, direction and/or specific production process; magnetic field ON/OFF)
 +
** Desirable thickness of the muon system to reach maximal muon detection efficiency
 +
** Muon/pion separation: decay point (detector layer) for the muons and pions of the same momentum, feasibility do distinguish processes based on this information
 +
** Feasibility to distinguish pion kink: distinguish muons originated from the primary vertex and produced in the pion decays. 
 +
 +
 
 +
=== The geometry used in version 0.1 simulation ===
 +
The following geometry is based mostly on the  [https://ctd.inp.nsk.su/docs/CTauFactory/Drawings/Detector/Detector_Model.pdf '''this drawings'''].
 +
 
 +
For simplicity, only the barrel part of the detector was simulated. The full simulation is possible, but not urgent for the first estimation.
 +
 
 +
* The following parts of the detector were simulated:
 +
** CsI calorimeter (inner radius 1090mm, thickness 297.6 mm, which corresponds to 16 X0)
 +
** Magnet coil (inner radius 1610mm, thickness 14.4 mm of copper, which corresponds to 1 X0)
 +
** 9 iron absorber layers in octant geometry (see the drawings). The distance to the innermost layer is 1900mm from the beamline, the thickness of the absorber layers 30 mm,30 mm,30 mm,40 mm,40 mm,80 mm,80 mm,80 mm, respectively,  which roughly corresponds to 1.7 X0, 1.7 X0, 1.7 X0, 2.3 X0, 2.3 X0, 4.5 X0, 4.5 X0, 4.5 X0.
 +
** The 30 mm gaps between the absorber layers could be filled with organic scintillator for the energy measurement (not needed for now)
 +
** Internal elements of the detector are estimated to give from 0.35 X0 to 0.6 X0 and are neglected in this study
 +
** Magnetic field is NOT simulated (could be turned on if needed, though)
 +
[[File:Det. 0000.png]] [[File:Det.gif]]
 +
*Green: CsI calorimeter
 +
*Violet: Cu magnet coil
 +
*White: Fe absorber layers
 +
*Yellow: Organic scintillator gaps
 +
 
 +
=== The source ===
 +
 
 +
Muons and pions were simulated as particles originated from the primary vertex according to the various physical distribution, taken from [https://ctd.inp.nsk.su/wiki/index.php/MC_Data_Sets '''these samples'''].
 +
To compare detector response in case of pions and muons, pion and muon of the same sign were generated in the same event with identical momenta. During the event reconstruction all hits were categorized according to the original particle.
 +
 
 +
=== Results ===
 +
 
 +
====Muon track smearing due to the multiple scattering====
 +
 
 +
One of the base characteristics of the muon system is spacial resolution. Since all charged tracks are subjected to the multiple scattering on the material of the inner detectors and absorber, and soft tracks scatters stronger, the use of the muon system as tracking detector is questionable.
 +
For the muon identification  purposes high spacial resolution is not needed, the optimal value should be of order of the scattering spot.
 +
 
 +
With the setup, described above, Monte Carlo simulation of the muon events for the 1st and 7th gap of absorber is shown. The smearing along Z axis is due to the multiple scattering on the calorimeter, magnet and absorber. Upper pictures corresponds to the energy-weighted position of the the hit, lower to the true position of muon/pion, respectively. The simulated physical process
 +
 
 +
<math>e^+e^- \to J/\psi \to \mu^+ \mu^-</math> at E=3.096 GeV.
 +
left: Magnetic field off ; right: magnetic field of 1.5 Tesla on
 +
 
 +
[[File:Mu 3.096 field off smearing.png]][[File:Mu 3.096 field on smearing.png]]
 +
 
 +
 
 +
The same for process
 +
 
 +
<math>e^+e^- \to  \tau^+ \tau^-</math> at E=3.3.55GeV.
 +
left: Magnetic field off ; right: magnetic field of 1.5 Tesla on
 +
 
 +
[[File:Tau 3.55 field off smearing.png]][[File:Tau 3.55 field on smearing.png]]
 +
 
 +
====Muon id feasibility study====
 +
 
 +
With the setup, described above, Monte Carlo simulation of the muon events for the 1st and 7th gap of absorber is shown. The smearing along Z axis is due to the multiple scattering on the calorimeter, magnet and absorber. Upper pictures corresponds to the energy-weighted position of the the hit, lower to the true position of muon/pion, respectively. The simulated physical process
 +
 
 +
<math>e^+e^- \to J/\psi \to \mu^+ \mu^-</math> at E=3.096 GeV.
 +
 
 +
Magnetic field off:
 +
 
 +
[[File:Mu 3.096 field on layers.png]]
 +
 
 +
Magnetic field 1.5 Tesla on:
 +
 
 +
[[File:Mu 3.096 field off layers.png]]
 +
 
 +
 
 +
The same for process <math>e^+e^- \to  \tau^+ \tau^-</math> at E=3.3.55GeV.
 +
 
 +
 
 +
Magnetic field off:
 +
 
 +
[[File:Tau 3.55 field on layers.png]]
 +
 
 +
Magnetic field 1.5 Tesla on:
 +
 
 +
[[File:Tau 3.55 field off layers.png]]
 +
 
 +
 
 +
== Version 0.2 ==
 +
 
 +
=== Tasks for the 0.2 version of the simulation tools ===
 +
 
 +
The main task for the second round of MC simulation is to supply the parametric simulation software with the data suitable for the fast simulation.
 +
The initial particle momenta range is [0.5,1.5] GeV/c, particle type: pi+,  pi-, mu+, and mu-.
 +
 
 +
 
 +
=== Geometry ===
 +
We used the same detector geometry as in the version 0.1, except for the fact, that  endcaps of the same layer structure have been added.
 +
The uniform magnetic field of 1.5T directed along beam line were applied to the volume inside the magnet coil.

Latest revision as of 21:39, 21 May 2019

Contents

[edit] This page is dedicated to the sketch standalone simulation of the SCT detector based on pure Geant4

One of the most important and urgent tasks for muon system technology choice and further optimization is to create a fast and reliable simulation tool. Optimal toolkit in this case is pure Geant4 simulation with simplified geometry and physics lists, able to produce fast results and estimate main parameters and characteristics of the detector.

[edit] Version 0.1

[edit] Tasks for the 0.1 version of the simulation tools

  • Create geant4-based geometry description reflecting all main features influencing muon-related physics. The model are to be extendible to describe muon detector in detail if needed.
  • Estimate basic features of the interacting muons (despite of the detector technology choice):
    • Muon smearing due to the multiple scattering (depending of the muon energy, direction and/or specific production process; magnetic field ON/OFF)
    • Desirable thickness of the muon system to reach maximal muon detection efficiency
    • Muon/pion separation: decay point (detector layer) for the muons and pions of the same momentum, feasibility do distinguish processes based on this information
    • Feasibility to distinguish pion kink: distinguish muons originated from the primary vertex and produced in the pion decays.


[edit] The geometry used in version 0.1 simulation

The following geometry is based mostly on the this drawings.

For simplicity, only the barrel part of the detector was simulated. The full simulation is possible, but not urgent for the first estimation.

  • The following parts of the detector were simulated:
    • CsI calorimeter (inner radius 1090mm, thickness 297.6 mm, which corresponds to 16 X0)
    • Magnet coil (inner radius 1610mm, thickness 14.4 mm of copper, which corresponds to 1 X0)
    • 9 iron absorber layers in octant geometry (see the drawings). The distance to the innermost layer is 1900mm from the beamline, the thickness of the absorber layers 30 mm,30 mm,30 mm,40 mm,40 mm,80 mm,80 mm,80 mm, respectively, which roughly corresponds to 1.7 X0, 1.7 X0, 1.7 X0, 2.3 X0, 2.3 X0, 4.5 X0, 4.5 X0, 4.5 X0.
    • The 30 mm gaps between the absorber layers could be filled with organic scintillator for the energy measurement (not needed for now)
    • Internal elements of the detector are estimated to give from 0.35 X0 to 0.6 X0 and are neglected in this study
    • Magnetic field is NOT simulated (could be turned on if needed, though)

Det. 0000.png Det.gif

  • Green: CsI calorimeter
  • Violet: Cu magnet coil
  • White: Fe absorber layers
  • Yellow: Organic scintillator gaps

[edit] The source

Muons and pions were simulated as particles originated from the primary vertex according to the various physical distribution, taken from these samples. To compare detector response in case of pions and muons, pion and muon of the same sign were generated in the same event with identical momenta. During the event reconstruction all hits were categorized according to the original particle.

[edit] Results

[edit] Muon track smearing due to the multiple scattering

One of the base characteristics of the muon system is spacial resolution. Since all charged tracks are subjected to the multiple scattering on the material of the inner detectors and absorber, and soft tracks scatters stronger, the use of the muon system as tracking detector is questionable. For the muon identification purposes high spacial resolution is not needed, the optimal value should be of order of the scattering spot.

With the setup, described above, Monte Carlo simulation of the muon events for the 1st and 7th gap of absorber is shown. The smearing along Z axis is due to the multiple scattering on the calorimeter, magnet and absorber. Upper pictures corresponds to the energy-weighted position of the the hit, lower to the true position of muon/pion, respectively. The simulated physical process 

e^+e^- \to J/\psi \to \mu^+ \mu^- at E=3.096 GeV.

left: Magnetic field off ; right: magnetic field of 1.5 Tesla on

Mu 3.096 field off smearing.pngMu 3.096 field on smearing.png


The same for process

e^+e^- \to \tau^+ \tau^- at E=3.3.55GeV. left: Magnetic field off ; right: magnetic field of 1.5 Tesla on

Tau 3.55 field off smearing.pngTau 3.55 field on smearing.png

[edit] Muon id feasibility study

With the setup, described above, Monte Carlo simulation of the muon events for the 1st and 7th gap of absorber is shown. The smearing along Z axis is due to the multiple scattering on the calorimeter, magnet and absorber. Upper pictures corresponds to the energy-weighted position of the the hit, lower to the true position of muon/pion, respectively. The simulated physical process

e^+e^- \to J/\psi \to \mu^+ \mu^- at E=3.096 GeV.

Magnetic field off:

Mu 3.096 field on layers.png

Magnetic field 1.5 Tesla on:

Mu 3.096 field off layers.png


The same for process e^+e^- \to \tau^+ \tau^- at E=3.3.55GeV.


Magnetic field off:

Tau 3.55 field on layers.png

Magnetic field 1.5 Tesla on:

Tau 3.55 field off layers.png


[edit] Version 0.2

[edit] Tasks for the 0.2 version of the simulation tools

The main task for the second round of MC simulation is to supply the parametric simulation software with the data suitable for the fast simulation. The initial particle momenta range is [0.5,1.5] GeV/c, particle type: pi+, pi-, mu+, and mu-.


[edit] Geometry

We used the same detector geometry as in the version 0.1, except for the fact, that endcaps of the same layer structure have been added. The uniform magnetic field of 1.5T directed along beam line were applied to the volume inside the magnet coil.

Retrieved from "https://ctd.inp.nsk.su/wiki/index.php?title=Muon_system_simulation&oldid=713"
Personal tools