Skip to content

Instantly share code, notes, and snippets.

View pdeperio's full-sized avatar

Patrick de Perio pdeperio

7 followers · 1 following
  • Kavli IPMU
  • Japan
View GitHub Profile
@pdeperio
pdeperio / naive_pvalue.py
Last active August 6, 2018 17:29
Naive p-value calculation
import numpy as np
poisson_mean = 7.56
poisson_sigma = 0.61
observation = 14
trials = int(1e7)
expectation = np.random.normal(poisson_mean, poisson_sigma, trials)
counts = np.random.poisson(expectation)
@pdeperio
pdeperio / conda_clone_6.8.0_err.log
Created August 2, 2018 21:47
conda_clone_6.8.0_err
$ conda create --name pax_v6.8.0.1 --clone pax_v6.8.0
An unexpected error has occurred.
Please consider posting the following information to the
conda GitHub issue tracker at:
https://github.com/conda/conda/issues
@pdeperio
pdeperio / po212_tracking.log
Created July 25, 2018 17:48
Po212 tracking log
*********************************************************************************************************
* G4Track Information: Particle = Po212, Track ID = 1, Parent ID = 0
*********************************************************************************************************
Step# X(mm) Y(mm) Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume ProcName
0 108 204 -149 0 0 0 0 LXe initStep
G4RadioactiveDecay::GetMeanLifeTime()
KineticEnergy: 0 GeV, Mass: 197.467 GeV, Life time: 431.366 ns
mean life time: 4.31366e-07 s
G4DecayTable: Po212
@pdeperio
pdeperio / kr83m_41kev_tracking.log
Created July 25, 2018 16:36
Kr83m 41 keV tracking log
*********************************************************************************************************
* G4Track Information: Particle = Kr83[41.557], Track ID = 1, Parent ID = 0
*********************************************************************************************************
Step# X(mm) Y(mm) Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume ProcName
0 281 -59.8 -252 0 0 0 0 LXe initStep
G4RadioactiveDecay::GetMeanLifeTime()
KineticEnergy: 0 GeV, Mass: 77.2341 GeV, Life time: 9.50448e+12 ns
mean life time: 9504.48 s
G4DecayTable: Kr83[41.557]
@pdeperio
pdeperio / kr83m_9kev_tracking.log
Last active July 25, 2018 16:35
Kr83m 9 keV tracking log
*********************************************************************************************************
* G4Track Information: Particle = Kr83[9.405], Track ID = 2, Parent ID = 1
*********************************************************************************************************
Step# X(mm) Y(mm) Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume ProcName
0 281 -59.8 -252 1.2e-07 0 0 0 LXe initStep
G4RadioactiveDecay::GetMeanFreePath()
KineticEnergy: 1.200241968035698e-10 GeV, Mass: 77.21570376305901 GeV, tau: 0 ns
G4Decay::GetMeanFreePath: Kr83[9.405] stops, kinetic energy = 0.0001200241968035698 keV
mean free path: 2.225073858506956e-311 m
@pdeperio
pdeperio / rn220.log
Created July 25, 2018 15:44
214BiPo Decay
******************************************************************************************************************************************
* Row * Instance * trackid * type * creaproc * edproc * time *
******************************************************************************************************************************************
* 9 * 41 * 9 * Bi212 * RadioactiveDecay * ionIoni * 13910.2592316261344 *
* 9 * 42 * 9 * Bi212 * RadioactiveDecay * Scintillation * 19550.3251246495711 *
* 9 * 43 * 16 * e- * RadioactiveDecay * eIoni * 19550.3251246495711 *
* 9 * 44 * 16 * e- * RadioactiveDecay * eIoni * 19550.3251246495711 *
* 9 * 45 * 16 *
@pdeperio
pdeperio / kr83m.log
Created July 25, 2018 15:41
Kr83m Decay
******************************************************************************************************************************************
* Row * Instance * trackid * type * creaproc * edproc * time *
******************************************************************************************************************************************
* 0 * 0 * 1 * Kr83[41.557] * Null * Scintillation * 5303.11377129179255 *
* 0 * 1 * 9 * e- * RadioactiveDecay * eIoni * 5303.11377129179255 *
* 0 * 2 * 9 * e- * RadioactiveDecay * Scintillation * 5303.11377129179255 *
* 0 * 3 * 8 * e- * RadioactiveDecay * eIoni * 5303.11377129179255 *
* 0 * 4 * 8 *
@pdeperio
pdeperio / g4.10.4.2_crash.log
Created July 24, 2018 21:51
G4.10.4.2 Crash
*** glibc detected *** xenon1t_G4p10: corrupted double-linked list: 0x00000000035be670 ***
======= Backtrace: =========
/lib64/libc.so.6[0x3722675f3e]
/lib64/libc.so.6[0x37226763c3]
/lib64/libc.so.6[0x3722679d08]
/lib64/libc.so.6(__libc_malloc+0x5c)[0x372267abfc]
/cvmfs/xenon.opensciencegrid.org/releases/anaconda/2.4/envs/mc/lib/libstdc++.so.6(_Znwm+0x18)[0x7f3d4c6aebe8]
/cvmfs/xenon.opensciencegrid.org/releases/anaconda/2.4/envs/mc/lib/libstdc++.so.6(_Znam+0x9)[0x7f3d4c6aec89]
/cvmfs/xenon.opensciencegrid.org/software/root/6.06.08/lib/libCore.so(_ZN7TBufferC2ENS_5EModeEi+0x67)[0x7f3d55ff30c7]
/cvmfs/xenon.opensciencegrid.org/software/root/6.06.08/lib/libRIO.so(_ZN11TBufferFileC2EN7TBuffer5EModeEi+0x9)[0x7f3d55a37f89]
@pdeperio
pdeperio / waveform_stitch.sh
Last active February 10, 2018 15:54
Waveform stitching
# Make sure to use same pax version as processing
source activate pax_v6.8.0
# Output directories
mkdir -p s1 s2 joined
RAWDIR=/project/lgrandi/xenon1t/raw_for_waveforms
# Specify SR to use same maps as processing
SCIENCE_RUN=1
@pdeperio
pdeperio / plot_bbf_overlay.C
Last active January 19, 2018 14:49
plot_bbf_overlay.C
{
TString lax_version = "1.2.3";
const int nHists = 4;
TString hname = "hbkg_py0_1.00_rf0_1.00";
TString hname_append[nHists] = {"_data", "", "_debug1", "_debug2"};
TString htitle[nHists] = {"Data", "MC (new, ef8ae8a)", "MC (debug1)", "MC (debug2)"};
enum henum {data, mc, mc_old};
int icolor[nHists] = {kBlue, kBlack, kRed, kGreen+2};