Transcript SoLID simulation - Shandong University

SoLID Simulation
Zhiwen Zhao 赵志文
University of Virginia
Third Workshop on Hadron Physics in China
and Opportunities in US
2011/8/9
• Introduction
• Simulation framework
• Simulation study
1
SoLID - Solenoidal Large Intensity Device
• One of three major new equipments of Hall A 12GeV upgrade
besides Super Bigbite (SBS) and Moller
• General purpose device.
• Physics
approved proposals:
1. PVDIS (Xiaochao Zheng’s talk)
2. SIDIS (Xin Qian’s talk)
Submitted proposals:
1. b1, deuteron tensor structure function
2. proton transversity
Proposals in preparation:
1. g3z, parity violating spin structure function
2. PVRES, parity violation in resonance region
2
Comgeant
The Past
• Geant3 based simulation program.
• geometry/sensitivity/digitization/field as input
files and detached from main code, run
different settings without recompilation.
• Successfully used for PVDIS and SIDIS
proposals.
3
GEMC (GEant4 MonteCarlo)
The Present and Future
http://gemc.jlab.org
https://hallaweb.jlab.org
/wiki/index.php/Solid_si
m_geant4
• C++ program that simulates particles through matter
using the Geant4.
• Successfully used for CLAS12.
• Detector information are stored in mysql database.
configuration changes are immediately available to
users without need of recompiling the code.
• Hit process factory: associate detectors with
external digitization routines at run time.
• perl script I/O to database, no need to know C++ or
Geant4 to build detector and run the simulation.
GEOMETRY,
BANKS,
DIGITIZATION
DATABASE
network
gemc
4
GEMC interface
Batch mode
5
GEMC interface
Detector
Interactive mode
Run Control
Camera
6
SIDIS with BaBar Magnet
Magnet/coil/yoke
7
SIDIS with BaBar Magnet
Target/Beam line
8
SIDIS with BaBar Magnet
GEM
9
SIDIS with BaBar Magnet
EC, large angle
10
SIDIS with BaBar Magnet
Collimator
11
SIDIS with BaBar Magnet
Cherenkov, light gas
12
SIDIS with BaBar Magnet
Scintillator
13
SIDIS with BaBar Magnet
Cherenkov, heavy gas
14
SIDIS with BaBar Magnet
MRPC (Multigap Resistive Plate Chambers)
15
SIDIS with BaBar Magnet
EC, forward angle
16
SIDIS with BaBar Magnet
3D Geant4
2D Geant3
17
SIDIS with BaBar Magnet
3D Geant4
2D Geant3
18
PVDIS with BaBar Magnet
2D Geant3
3D Geant4
Baffle
19
PVDIS with BaBar Magnet
2D Geant3
3D Geant4
Baffle
20
SoLID GEMC Framework
• geometry/sensitivity/digitization/field in
mysql database.
• Customized hit processing for various
detectors.
• Unified individual detector simulation and
the whole SoLID simulation.
• GEMC can be used for other projects.
21
Simulation Study
•
•
•
•
•
•
•
•
Magnet Option
SIDIS Kinematics Study
PVDIS baffle design and FOM
Background rate and GEM tracking
Energy flux and EC
Cherenkov
Neutron background
Other progress
22
Magnet Option
Poisson 2D
BaBar
CDF
CLEO
ZEUS
23
Magnet Comparison
Glue-X
CLEO
ZEUS
CDF
Old
SLAC
Cryostat
Inner
Radius
150 cm
150 cm
86 cm
150 cm
90 cm
Length
345 cm
350cm
245cm
500 cm
350 cm
Central
Field
1.49T
1.5T
1.8T
1.47T
2T
Flux
Return
Iron
Yes
Yes
No
No
No
Cool Icon
Yes
Yes
Yes
No
No
Variation
in Current
density
with z
2x more
in end
than
central
4.2%
more in
end than
central
40% more
in end
than
central
No
Yes
Available
02 June 2011
Probably
Not??
Probably
Probably
Probably
Paul E. Reimer
New
New
Whatever we need
BaBar
Yes
One will be available
?
24
SIDIS Kinematic Coverage@11GeV
• Green area,
large angle coverage
• Black area,
forward angle coverage
25
SIDIS Kinematic Coverage@11GeV
ZEUS
BaBar/CLEO
CDF
Glue-X
x
0.05-0.58
0.05-0.65
0.05-0.64
0.05-0.64
z
0.3-0.7
0.3-0.7
0.3-0.7
0.3-0.7
Q2
1-6
1-9
1-7.2
1-8
W
2.3-4.2
2.3-4.4
2.3-4.2
2.3-4.2
W’
1.6-3.4
1.6-3.5
1.6-3.4
1.6-3.4
PT
0-1.45
0-1.7
0-1.45
0-1.45
PVDIS Baffle
Reduce background by 50
27
PVDIS FOM
28
Background Rate on GEM for SIDIS
Condition: 15uA 11GeV e- beam, 40cm 3He 10amg gas target
Todo: more realistic GEM module description in progress,
borrowed from SBS simulation.
29
Tracking Progressive Method (curved tracks)
PVDIS, based on simulated background on GEM
3/4
3/4
No EC
Add EC, with BG Single/Multi : 97.5/0.27% time: 100 s
30
Energy Flux Rate on EC
for SIDIS with BaBar Magnet
60krad/y
Forward angle
Large angle
Condition: 15uA 11GeV e- beam, 40cm 3He 10amg gas target
Todo: more careful study of hadron energy flux in progress 31
IHEP 2010 module
EC (Shashlik)
•
•
•
•
•
•
•
•
Dimensions
Radiation length
Moliere radius
Radiation thickness
Scintillator thickness
Lead thickness
Radiation hardness
Energy resolution
38.2x38.2 mm2
17.5mm
36mm
22.5 X0
1.5mm
0.8mm
500 krad
6.5%/√E 1%
32
EC (Shashlik) transverse size
Rough numbers only
w/ 50ns ADC gate
4
4.0
3.5
3
2.5
2.6
2.4
Resolution. (cm)
2
Back Ground (%)
1.7
1.5
Cost (M$)
1
1.0
0.9
0.6
0.7
0.6
0.5
0.4
0
0
5
10
15
20
25
block Size (cm)
33
SIDIS Cherenkov: Optics
R
One spherical mirror
xi  incident ray on mirror
2
1
1 xr  reflected ray
cos  (

)
xi
xr   angle between incident ray
and normal to the mirror
SIDIS Cherenkov: Detector
 (Some) Requirements: 1) resistant in magnetic field
3) decent size
2) “quiet”
Used by PHENIX successfully
Gaseous Electron Multiplier + CsI
• GEMs + CsI: resistant in magnetic field, size is not a problem
• Consists of 3 layers of GEMs,
first coated with CsI which acts as
a photocathode
• First GEM metallic surface
overlayed with Ni and Au to ensure
stability of CsI (CsI not stable on
Cu)
The simulation shows good collection
efficiency.
SIDIS L.-G. Cherenkov: Photon Detector
 (Some) Requirements: 1) resistant in magnetic field
3) decent size
2) “quiet”
Photomultiplier Tubes
• Multi-anode 2” PMT: fairly resistant in magnetic field; it can be
tiled (data from Hamamatsu)
1.93” effective area (94%) Square shaped and 94%
effective area: ideal for
tiling
2.05”
Initial test shows we can
safely run at less than 70G
The simulation gives us the
guidance of local magenetic
field where the PMT is located.
Neutron Background
Damage function
FLUKA
Shielding:
Polyethylene
37
Neutron Background
Shielding reduces neutron flux
in half at two test locations
38
Other Progress
• GEM
A small prototype was tested at Jlab. Combined Efforts
from UVa/INFN/Jlab/China are in coincidence with
GEM R&D for the SuperBigbite & EIC. Several large
prototypes are being built in US and China.
• MRPC
Chinese collaborators will come onsite for beam test
later this year.
• DAQ
Collaborating with Hall D.
39
Summary
• SoLID collaboration has successfully adopted
GEMC as its Geant4 simulation framework
and joined in GEMC development. The
simulation is ready to be used for various
studies to help detector design.
• A lot of subsystem design and simulation
progresses have been made. More studies are
under way.
• In preparation for the director review.
40
Thanks
• Maurizio Ungaro (GEMC)
• Paul Reimer (Magnet, Calorimeter)
• Seamus Riordan (Baffle, PVDIS FOM)
• Lorenzo Zana (Neutron BG)
• Simona Malace, Eric Fuchey, Yi Qiang (Cherenkov)
• Jin Huang, Mehdi Meziane (Calorimeter)
• Yang Zhang (SIDIS kinematics)
• Eugene Chudakov (Comgeant PVDIS, Baffle)
• Xin Qian (Comgeant SIDIS, tracking)
• SoLID Collaboration
41
Backup
42
How To: new detector, hits
$detector{"pos"}
$detector{"rotation"}
$detector{"color"}
$detector{"type"}
$detector{"dimensions"}
$detector{"material"}
$detector{"mfield"}
$detector{"ncopy"}
$detector{"pMany"}
$detector{"exist"}
$detector{"visible"}
$detector{"style"}
$detector{"sensitivity"}
$detector{"hit_type"}
$detector{"identifiers"}
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
”10*cm 20*cm 305*mm";
"90*deg 25*deg 0*deg";
"66bbff";
"Trd";
”1*cm 2*cm 3*cm 4*cm 5*cm";
"Scintillator";
"no";
12;
16th: Bank
1;
1;
1;
1;
"CTOF";
17th: Digitization Routine
"CTOF";
"paddle manual 2";
In general, 1 bank  1 digitization routine… but not necessary
43
Factory Method for Hit Processes
Hit Process,
Digitizations
External
Routines
SVT
CTOF
Automatic Process
Routines Still External
gemc DC
gemc
FTOF
Easy to:
• add new routine
• debug
• modify
44
Digitization
Available For every G4 step
• Hit Position
• Volume Local Hit Position
• Deposited energy
• Time of the hit
• Momentum of the Track
• Energy of the track
• Primary Vertex of track
• Particle ID
• Identifier
• Mother Particle ID
• Mother Vertex
Hit Process Example
Average (x,y,z)
Average (lx, ly, lz)
Total E
Average t
Average p (final p)
Energy
Primary Vertex of track
Particle ID
Strip, Layer, Sector
Event Generation
1)Particle gun built in, two luminosity beams
can be added
2)LUND Format (txt) for physics events
Data Output
1)evio, bank alike binary format by Jlab DAQ
group
2)Root tree, convert from evio
3)text
46
SoLID Event Generator
• DIS e- and pion generators are ready in C++
• e- and pion coincidence generator is ready in
C++
SoLID Hit Processing
• FLUX, raw, EC …
47
SoLID Simulation Databasre
soliddb.jlab.org
• Mysql 5 cluster server. It is highly efficient and
has minimum downtime.
• Flexible development structure.
48
Documentation
• gemc.jlab.org
• https://hallaweb.jlab.org/
wiki/index.php/Solid_sim_
geant4
49
Netpbm
50
CLAS12 SVT
51
GEMC update
Progress
• Mirrors, done in the “identifiers” entry of the geometry, control
optical property on fly.
• Right click to output geometry in GDML format.
• Mother particle tracking becoming optional to optimize speed.
Todo list
•
•
•
•
Move material definition into database also.
Move svn repository out of clas12svn and restructure.
Improve database I/O.
Adapt to Geant4.9.4.
52
SoLID GEMC update
Progress
• Add “solid” HIT_PROCESS_LIST
• More database added in soliddb.jlab.org to allow for the full SoLID, its
subsystems simulation. Also database for individual developers.
• PVDIS and SIDIS yoke designs and field maps are unified
• More materials added for our setup.
• More instructions on wiki
• Rewrote many geometry to avoid overlap and added more
• EC simulation in GEMC is under work.
• Baffle redesign for various magnets
• Event generators updated for PVDIS and SIDIS
• Study configuration with ZEUS magnet.
Todo list
• Move subsystem simulation to GEMC
• Customize hit routine
• Direct root output
53
Compare geant4 to geant3 results
Progress
• SIDIS kinematics and angle distribution
• SIDIS and PVDIS low energy background rate.
Todo list
• Acceptance
• Detector resolution
54