Skip to content

Instantly share code, notes, and snippets.

View JeremyMcCormick's full-sized avatar
:octocat:

Jeremy McCormick JeremyMcCormick

:octocat:
View GitHub Profile
from lsst.rsp import get_tap_service
from IPython.display import display, HTML
def print_test(msg, color="black"):
display(HTML(f'<span style="color: {color}">{msg}</span>'))
def print_error(msg, object_name, exception):
display(HTML(f'<span style="color: red; font-weight: bold;">Error {msg}: {object_name}</span>'))
display(HTML(f'<span style="color: red;">{exception}</span>'))
"""Example of using Great Expectations to validate ci_imsim test
data.
"""
import lsst.daf.butler as dafButler
import great_expectations as gx
from great_expectations.data_context import DataContext
from great_expectations.core.batch import BatchRequest
# Get the test data frame from the Butler
@JeremyMcCormick
JeremyMcCormick / gist:602e424159452aa1f67aee125d122e30
Created December 8, 2022 04:19
RimWorld Research Advancing Stacktrace
[Tech Advancing] [Information] Flushing old Research Manager values.
UnityEngine.StackTraceUtility:ExtractStackTrace ()
Verse.Log:Message (string)
TechAdvancing.LogOutput:WriteLogMessage (TechAdvancing.Errorlevel,string)
TechAdvancing.TA_ResearchManager:FlushCfg ()
TechAdvancing.WorldCompSaveHandler:ExposeData ()
Verse.ScribeExtractor:SaveableFromNode<RimWorld.Planet.WorldComponent> (System.Xml.XmlNode,object[])
Verse.Scribe_Collections:Look<RimWorld.Planet.WorldComponent> (System.Collections.Generic.List`1<RimWorld.Planet.WorldComponent>&,bool,string,Verse.LookMode,object[])
Verse.Scribe_Collections:Look<RimWorld.Planet.WorldComponent> (System.Collections.Generic.List`1<RimWorld.Planet.WorldComponent>&,string,Verse.LookMode,object[])
RimWorld.Planet.World:ExposeComponents ()
@JeremyMcCormick
JeremyMcCormick / cv_model_selection.R
Created October 23, 2021 20:27
Model selection using CV
# snippet: not all code included
# revised loop to do cross-validation
for (m in 1:length(allModels)) {
allpredicted = rep(NA,n) # storage for honest predictions
for (ii in 1: nfolds) { # ii is an easier string to search for index
groupii = (cvgroups == ii)
trainset = bodyfat[!groupii,] # all data EXCEPT for group ii
testset = bodyfat[groupii, ] # data in group ii
data = read.csv("mydata.csv")
n = nrow(data)
nfolds = 5 # 5-fold CV
groups = rep(1:nfolds, length=n)
set.seed(2)
cvgroups = sample(groups, n)
allpredicted = rep(NA, n)
"""
Connect to an MS SQL Server and execute a query
"""
import pyodbc
conn = pyodbc.connect('Driver={SQL Server};'
'Server=THINKSGIVING\SQL;'
'Database=AdventureWorks2017;'
'Trusted_Connection=yes;')
@JeremyMcCormick
JeremyMcCormick / print_loc.sh
Created September 9, 2020 21:24
Print non-blank LoC
for f in $(find src/ include/ hps_eve.cxx -type f); do r=$(cat $f | awk 'NF' | wc -l); echo $f:$r; done
auto vec = step->GetPreStepPoint()->GetPosition();
const double pos[] = {vec[0], vec[1], vec[2]};
double bfield[] = {0, 0, 0, 0};
G4StepPoint* pre = step->GetPreStepPoint();
G4FieldManager* globalFieldMgr = G4TransportationManager::GetTransportationManager()-> GetFieldManager();
globalFieldMgr->GetDetectorField()->GetFieldValue(pos, bfield);
std::cout << ">>>> Field val: " << pre->GetPosition() << " -> ("
<< bfield[0] << ", " << bfield[1] << ", " << bfield[2] << ") " << std::endl;
#include <vector>
#include <iostream>
#include "TEveManager.h"
#include "TEveBrowser.h"
#include "TEveGeoNode.h"
#include "TEveScene.h"
#include "TEveViewer.h"
#include "TGLViewer.h"
#include <vector>
void MakeGeometry() {
TGeoManager *geom = new TGeoManager ("simple1", "Simple geometry");
//--- define some materials
TGeoMaterial *matVacuum = new TGeoMaterial ("Vacuum", 0, 0, 0);
TGeoMaterial *matAl = new TGeoMaterial ("Al", 26.98, 13, 2.7);
//--- define some media
TGeoMedium *Vacuum = new TGeoMedium ("Vacuum", 1, matVacuum);