This wiki is now closed and kept for historical purposes. Please visit the new wiki at


Jump to navigationJump to search

Welcome to the water Cherenkov simulation page!


  • Oh how exciting this box is
  • Surely there will be fantastic news in this box soon.
  • Working on this page
  • Please add info to the page!

Event Display

WCSim is a general Geant4 based water Cherenkov simulation package.
This Wiki page describes the current WC simulation effort.
For more information you can contact Conveners Ion Stancu and Chris Walter
The group Mailing list is here:

Under heavy construction - CWW

How to Get WCSim

The code is stored in SVN and requires ROOT and Geant4.
For instructions on how to check out and compile the code, see WCSim.

Simulation Detector Information

Reference designs include a right rectangular prism ("mailbox") and a cylindrical detector volume. These water volumes are further subdivided by Phototube coverage and efficiency.



Here is info on the ~5 reference designs and their configurations. For the 150 kton Mailbox:

The FV is 68 x 40 x 55 m^3 (149.6 metric tons) , and the total water volume (just add 2*2.5m = 5m to each dimension) =73 x 45 x 60 m^3 = 197.1 metric tons (equating volume entirely to water).

Visualization Software

A event display is available from the svn repository:

svn co event3D

  • The glui-dev librairies are needed for this installation. On Ubuntu these can be found in the Synaptic Package Manager or by just typing:

sudo apt-get libglui-dev

  • In MacOS 10.6, the GLUI librairies involve a little more work. They can be found here The GLUI Framework for Mac OS X installs into /Library/Frameworks by default. You should copy GLUI.framework from there into /System/Library/Frameworks for a successful installation.
  • There is currently a known bug with people with Intel architecture in Linux where the screen flickers. Solutions are currently being sought.

M. Askins Event Display

Simulation Physics Information

Geant requires the user to specify which physics will be included in the simulation.

See WCSim_Physics_Lists for a table of the included physics particles and processes.

Where to find data and vector files

Lists of current samples simulated at Fermilab for the 200 kton reference design, with and without light collectors, as well as older designs with beam neutrinos, atmospheric neutrinos, cosmic rays, and single particles, are available at:

Fermilab WCSim samples

The vector files for Mono-energetic data can be found at:


Monoergetic Isotropic Samples : Electron, Muon, Pion(pizero), Piplus, Piminus

  • WCSim svn release number = 281
  • Geometry:
/WCSim/WCgeom DUSEL_100kton_10inch_HQE_15perCent
  • WCSim Rootfiles available at ssh -Y -A -t <username> ssh lbne0001 in the directory
  • Input vector files are in
  • File naming convention:
[particle type]_vtx_[energy in MeV][section number from 1 to 10].root


File Structure:


  • Particle_name : muon, pions, electron
  • Energy in MeV : 400 to 4000 in steps of 400 MeV
  • Filenumber: 1-10

Each file contains 1000 events. There are 10 such files in each sub directory. In total 300K events were simulated and can be accessed at FNAL and BNL.

Geant4 Simulated Data at BNL :

WCSim Vesrion: 260

WCSim parameters :

Geometry : Cylinder

/WCSim/WCgeom DUSEL_100kton_10inch_HQE_30perCent

Login into BNL machine:

ssh -Y -A -t <username> ssh lbne0001

You would find data in directory :

ParticleName || Energy range || Location

Electron || 400 Mev - 4 GeV || /data1/lbne/pgupta/electron

Pizeros || 400 Mev - 4 GeV || /data1/lbne/pgupta/pion

Muon || 400 Mev - 4 GeV || /data1/lbne/pgupta/muon

File Structure:


Geant4 Simulated Data at FNAL :

Geant4 Cosmic Muon Simulated Data:

A sample of cosmic muons has been simulated using a 3-dimensional Geant4 model of the Homestake overburden. The detector is taken to be at an (easting,northing) of (3618.4,-2499.7), where (0,0) corresponds to 44.359086 degrees latitude and -103.76592 degrees longitude. The documentation for this simulation can be found on DocDB #3144-v3.

The muon flux at a point 3 meters above the cavern ceiling, an elevation of 491.2 feet, is estimated to be 6.85x10^-5 Hz/m^2. Kinematic information for the surviving muons has been written to a root file which can be downloaded from :

In the ntuple x,y,z (in km),E (in TeV), costheta, phi are muon surface parameters. xdet,ydet,zdet (in m), px,py,pz (momentums in TeV) are parameters at the detector depth, and charge is the muon charge (+1/-1). This sample represents an effective lifetime of approximately 55,108 seconds.

Vector and root file descriptions

Input vector file description

Output root file description

Documents, information and talks

Important WCSim related DOCDB references

  • Johannes geometry document
  • describing root output format
  • spreadsheet by John
  • Link to all documents tagged WCSim

Links to video meetings on the DOCDB [Can I embed?]

Link to other relevant documents

  • Kadel geometry document
  • Hamamatsu tube info

Todo lists, whiteboards, etc.


For communicating we are using EVO.... BNL mailing lists

  • How to connect to an EVO meeting (link to elsewhere on Wiki?)
  • Link to mailing list page
  • Link to archives

Useful external links