Skip to content

Combined SuperSim and G4DMC simulation of 4mm test wafer in SLAC B33 cryostat

Notifications You must be signed in to change notification settings

kelseymh/SLACB33_G115

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

16 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

From Noah Kurinsky:

I think we should aim to simulate G115, a 3”x4mm Ge detector with the SNOLAB
QETs, and compare the simulation to the data taken in the Bluefors at
SLAC. Below is a list of tasks and timeline for getting this done, and I'm
hoping you can take the lead on the G4CMP part of the simulation. We can
either try to put the G4CMP part into the SuperSim simulation Mike already
has, or have mike build the input spectrum for the G4CMP run and modify the
example geometry to run this crystal.

Once we get the list of phonon hits, we will pass off to Jon and Jorge; Jon
will run the TES simulation and Jorge will do post-processing. Xuji and I
are taking/analyzing the SLAC data. The idea is to have this done or mostly
done by the collaboration meeting as a first simulation milestone for
SNOLAB detectors and a first sanity check with detectors of the G4CMP.

This, plus the charge validation I've done at Stanford, constitutes the
first full test of G4CMP charge and phonon simulations, and will allow us to
go forward with larger detectors at voltage when we have COMSOL models and
data from test facilities (we already have a lot of HV data from UMN, just
with high noise).

-----
Timeline for 0V Simulation

Model Backgrounds - identify sources - DONE (Mike, Noah)

SuperSim - Need to generate ~100000 detector hits given background spectrum,
hopefully done by 8/1 – Mike (should be quick)

Data - need to take pulse data to build up moderate statistics sample,
hopefully ~10,000 pulses free of background/muons data taken by 8/4 - Noah

G4CMP - need to implement G115 geometry, just needs to be a cylinder 3" in
diameter 4mm in thickness, located at top of stack - Tyler (with help from
Mike for SuperSim)

Will take care of charge and phonon absorption

Aim to be done by 8/4, run jobs over next ~3 days

TES Simulation + Postprocessing - Goal of being done by collaboration
 meeting plenary (8/13)

Modify TES simulator to read in example crystal outputs, segment into iZIP
channel layout. G115 looks like a CDMSlite detector - Jon

Do CDMSlite postprocessing on data (Xuji + Noah) and DMC (Jorge)
-------

		  ========================================

Building and running for end-to-end simulation

0)  Configure your session with Geant4 environment variables.

	$ source <G4-path>/share/Geant4-10.x.y/geant4make/geant4make.csh

1)  Check out and build G4CMP "develop" branch, source g4cmp_env.* script

	$ git clone https://github.com/kelseymh/G4CMP
	$ make all
	$ source g4cmp_env.csh	(or .sh for Bash)

    You make also build using CMake, following the README instructions.

2)  Check out and build SuperSim "develop" branch, source CDMSbuild/g4setup.*

	https://confluence.slac.stanford.edu/display/CDMS/SuperCDMS+Simulation

	$ git clone ssh://titus.stanford.edu/data/git/Simulations/supersim
	$ CDMSbuild/setup
	$ source CDMSbuild/g4setup.csh
	$ make all

3)  Switch back to this directory.

4)  Generate raw gammas from B33 floor and process for rethrowing.  Edit
    file to select either G115 or G124 (Zip).

	$ raw_gammas.csh

5)  Move rethrowing file into data/ directory

	$ mv atG115_gamma.dat data/
	$ mv atG124_gamma.dat data/

6)  Generate SuperSim hits output using rethrowing spectrum

	$ CDMS_SimDraw G115_rethrow.mac >&! G115_rethrow.log
	$ CDMS_SimDraw G124_rethrow.mac >&! G124_rethrow.log

7)  Edit DMC job (G115_DMC.mac) to point to new SuperSim hits file, then run

	$ CDMS_G4DMC G115_DMC.mac >&! G115_DMC.log
	$ CDMS_G4DMC G124_DMC.mac >&! G124_DMC.log


SuperSim macros for G115 simulation
-----------------------------------

SLACB33_G115.mac:	Defines cryostat geometry with custom ZIP detector
SLACB33_G124.mac:	As above, but makes regular G124 active, G115 not.
G115_zip.mac:		Define custom ZIP detector dimensions (3" x 4 mm)
G115_draw.mac:		Cutaway view of B33 geometry, SLACB33_G115.eps
gamma-floor.mac:	Draws gamma spectrum from floor toward ZIPs

G115_gammas.mac:	Throw gammas from floor for event production (low eff)

raw_gammas.mac:		Generate spectrum of gammas at G115 or G124 from floor
G115_rethrow.mac:	Rethrow gammas at G115 using raw_gammas.mac result
G124_rethrow.mac:	Rethrow gammas at G124 using raw_gammas.mac result
G115_DMC.mac:		Use G115_rethrow output to generate G4DMC results
G124_DMC.mac:		Use G115_rethrow output to generate G4DMC results

raw_gammas.csh:		Runs raw_gammas.mac and CDMSspectrum.C on output

hit_plot.C
event_sum.C
lowE_events.C
gamma_flux.C:		Produce single-page plots from SuperSim mczip6 TTree

About

Combined SuperSim and G4DMC simulation of 4mm test wafer in SLAC B33 cryostat

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published