davidchall

# Editor: https://shinylive.io/py/editor/#gist=697a9f02c3c000de6ae4b10ac1e725eb
# App: https://shinylive.io/py/app/#gist=697a9f02c3c000de6ae4b10ac1e725eb

from pathlib import Path
from shiny import ui, render, App, Inputs, Outputs, Session


calories_target = 200

app_ui = ui.page_fluid(

davidchall

<?xml version="1.0" encoding="utf-8"?>
<style xmlns="http://purl.org/net/xbiblio/csl" class="note" version="1.0" demote-non-dropping-particle="display-and-sort" page-range-format="chicago">
  <info>
    <title>Chicago Manual of Style 17th edition (full note, including subsequent)</title>
    <id>http://www.zotero.org/styles/chicago-fullnote-bibliography-including-subsequent</id>
    <link href="http://www.zotero.org/styles/chicago-fullnote-bibliography" rel="template"/>
    <link href="http://www.chicagomanualofstyle.org/tools_citationguide.html" rel="documentation"/>
    <author>
      <name>Julian Onions</name>
      <email>julian.onions@gmail.com</email>

davidchall

#!/bin/bash

# aliases
alias e='emacs -nw'
alias nb='jupyter notebook'
alias lab='jupyter lab'
alias chrome="open -a Google\ Chrome.app"
alias finder="open -a Finder"
alias preview="open -a Preview"
quicklook () { qlmanage -p "$*" >& /dev/null; }

davidchall

Index: COMIX/Main/Single_Process.C
===================================================================
diff --git a/branches/COMIX/Main/Single_Process.C b/branches/COMIX/Main/Single_Process.C
--- a/branches/COMIX/Main/Single_Process.C	(revision 31007)
+++ b/branches/COMIX/Main/Single_Process.C	(revision 31008)
@@ -265,8 +265,8 @@
 	  for (size_t j(0);j<nfmap;++j) {
 	    long int src, dest;
 	    *map>>src>>dest;
-	    Flavour ft((kf_code)(abs(src)),src<0);

davidchall

diff -Naur geant4.10.02.p01/source/geometry/solids/specific/src/G4TessellatedSolid.cc geant4.10.02.p01-topas3/source/geometry/solids/specific/src/G4TessellatedSolid.cc
--- geant4.10.02.p01/source/geometry/solids/specific/src/G4TessellatedSolid.cc	2016-03-02 03:24:12.000000000 -0500
+++ geant4.10.02.p01-topas3/source/geometry/solids/specific/src/G4TessellatedSolid.cc	2016-05-13 11:42:21.000000000 -0400
@@ -1674,7 +1674,7 @@
                                           const G4ThreeVector& v)const
 {
   G4double dist = DistanceToInCore(p,v,kInfinity);
-#ifdef G4SPECSDEBUG
+//#ifdef G4SPECSDEBUG
   if (dist < kInfinity)

davidchall

#!/bin/sh

# This script will (partially) convert TOPAS parameter files from v2.0 to v3.0.
# A backup of the original file is created with a .bak extension.
# WARNING: using this script twice on the same file will delete the backup file.

if [ "$#" -eq 0 ]; then
    echo "Usage: $0 filename"
fi

davidchall

Given a set of covariance matrices pcov of shape (n,m,m) and a set of parameter estimates popt of shape (n,m), the mean parameter estimates are:

weight = np.linalg.inv(pcov)
pcov_mean = np.linalg.inv(np.sum(weight, axis=0))
popt_wgt = np.einsum('ijk,ik->ij', weight, popt)

popt_mean = np.dot(pcov_mean, np.sum(popt_wgt, axis=0))

davidchall

#!/usr/bin/env python

###############################################################
##  Author:   David Hall
##  Email:    david(dot)hall(dot)physics(at)gmail(dot)com
##  Date:     March 1st, 2016
##  Program:  topas_batch_uncertainty
##
##  This program reads a set of TOPAS result files and uses 
##  the batch method to compute the relative statistical

davidchall

#!/usr/bin/env python

###############################################################
##  Author:   David Hall
##  Email:    David(dot)Hall(at)physics(dot)ox(dot)ac(dot)uk
##  Date:     15th August 2013
##  Program:  top2root
##
##  This program reads in a topdrawer file and converts the
##  plots contained within to the ROOT histogram format.