WCSim Update Erin OSullivan, on behalf on the WCSim group HK Open - - PowerPoint PPT Presentation

WCSim Update

Erin O’Sullivan, on behalf on the WCSim group HK Open Meeting January 30, 2014

An Introduction to WCSim

What is WCSim?

WCSim: An open-sourced, GEANT4-based simulation code for water Cherenkov detectors.

• User can define a detector

configuration by choosing a PMT type, cylinder size, option to add Gadolinium loading. Also, the HK detector is defined here.

• Run particles in your detector

using either a user-defined vector file or using GEANT4 particle generation.

How do I get WCSim?

• 1. git clone the development branch from
• GitHub. This will give you the most up-to-

date version of the code.

• 2. Download a zip file of a release (currently

v1.2.0). This will give you a stable version of the code.

• 3. Download it as part of the Hyper-K

software package. 3 Different Methods

git clone commands for https, ssh or svn Get the zip or tar file of the newest release (v1.2.0)

https://github.com/WCSim

Fork a copy of the code

https://github.com/WCSim

Browse the code Get emails when there are new pull requests or issues See ongoing issues or report bugs you find See any outstanding pull requests

How do I implement a change in WCSim?

• Make a new branch on your local forked copy
• Push your new branch to GitHub
• Make a pull request and discuss!
• Modify your pull request based on feedback
• Your pull request is merged

How do I implement a change in WCSim?

Where can I get help with WCSim?

• If you find a bug or have a suggestion

for how to enhance the code, you can

• pen an issue on GitHub.
• If you have a general code question,

you can email the HK software mailing list: hkwg-software@suketto.icrr.u- tokyo.ac.jp

v1.2.0: Changes to the code structure

Motivation

• Separating the code allows for different

developers to work on the individual components

• Is a move toward making the electronics in

WCSim more flexible by making it easier in the code to implement a new type of PMT or digitizer.

Separating the PMT and Digitizer in WCSim

True Info: Contains info about which PMTs were hit, hit time, etc

Separating the PMT and Digitizer in WCSim

Old setup in WCSim: PMT and digitizer together. “Digitizer”: Contains info about the PMT (ex. charge conversion), the digitizer (ex. gating, threshold), and information about the PMT and digitizer together (ex. electronics response)

Separating the PMT and Digitizer in WCSim

New setup in WCSim: PMT and digitizer are separate. PMT: Converts a hit to the corresponding charge. True Info: Contains info about which PMTs were hit, hit time, etc Digitizer: Contains info specific to the digitizer

Making the PMTs in WCSim more modular

Motivation

• Implementing this change will allow users to

“plug in” any of the predefined PMT types when they create a detector geometry.

• This will also allow all of the PMT

information to be in one place, making it more straight forward for users who want to implement a different type of PMT.

Making the PMTs in WCSim more modular

WCSimPMTObject Base Class:

20” PMT 8” PMT

GetPMTName GetPMTRadius GetPMTExposeHeight

10” PMT 10” HQEPMT 12” HQEPMT 20” HQEHPD

GetPMTGlassThickness GettimingResolution GetQEWavelength GetQE GetmaxQE Getqpe0 Derived Classes: Member Functions:

Making the PMTs in WCSim more modular

WCSimDetectorConfigs.cc: Possibilities for CreatePMTObject are PMT20inch, PMT8inch, PMT10inch, PMT10inchHQE, PMT12inchHQE, HPD20inchHQE

v1.2.0: New features in WCSim

Adding HPDs to WCSim

• 1 p.e. distribution: measured value of 8 inch HPD

in EGADS tank (champion data)

• Quantum Efficiency: QE of High QE SK PMTs

which were the candidates of PMTs in EGADS.

• Timing resolution: measured value of 8 inch HPDs
• Photodetector shape: based on the design drawing

Information from Okajima-san

New Detector Configuration for Hyper-K w/HPD

Some small add-ons

• Alex has added a python script which

generates vectors for monoenergetic single particles in the SK geometry

• There is now a crude verification script

which looks at the basic outputs of WCSim (charge, time, number of tubes hit) and compares them to the output of a clean check out of the code.

Ongoing work: More changes in the electronics scheme

FADC in WCSim

Information from Thomas Lindner

Making the digitizer more flexible in the code and implementing the FADC to see if there are any benefits to using an FADC over the TDC/ADC electronics

How can I get involved?

Simulate the SK-IV Detector in WCSim

Other Open Issues on GitHub

• Document and internally comment the HK geometry

section of the code

• Check/validate the Gd model
• Check/validate the physics list
• Check if code is compatible with newest Geant4

release

• Add a validation suite

Conclusion

• WCSim is now on GitHub
• Many changes to the electronics in the code to

make it easier to swap out different PMT types

• HPDs are now available in the code and the Hyper-

K with HPD geometry is available for use.

• Please get involved!

Wave Length [nm] 300 350 400 450 500 550 600 650 700 QE [%] 5 10 15 20 25 30 35 40

ZP0007 ZP0012 ZP0014 ZP0015 ZP0021 ZP0022 ZP0024 ZP0025 Normal SK PMT

High-QE R3600

Normal-QE R3600

×1.5

Charge [p.e.] 5 10 15 20 25 30 Timing resolution [ns] 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

HPD timing resolution (measured value) HPD timing resolution (WCSim parameter) PMT timing resolution (WCSim parameter)