mgmarino

TGaxis* rescale_frame(TPad *canvas, RooPlot *frame, Double_t scale, const char *title)
{
 
 //   HACK
 //   Takes a frame and rescales it to arbitrary coordinates, writing over the graphical axis.
 //   This is helpful when dealing with RooPlot to get the axes
 //   correct.  Returns axis in case anything else needs to be done.

    TAxis* yaxis = frame.GetYaxis();
    yaxis->SetTickLength(0);

mgmarino

import ROOT
ROOT.gSystem.Load("libEXOROOT")
t1 = ROOT.TChain("tree")
t2 = ROOT.TChain("tree")


t1.Add("recon00001620-000.root")
t1.Add("recon00001683-000.root")

t2.Add("run00001620-000.root")

mgmarino

import ROOT
ROOT.gSystem.Load("libEXOROOT")
t1 = ROOT.TChain("tree")
t2 = ROOT.TChain("tree")

ED1 = ROOT.EXOEventData();
ED2 = ROOT.EXOEventData();

t1.Add("recon00001620-000.root")
t1.Add("recon00001683-000.root")

mgmarino

from load_tree import load_tree, load_limited_tree

# use the following if you want a lot of events (~800000)
atree = load_tree()

# use the following if you want fewer events (~30000)
atree = load_limited_tree()

# I would use the second to do tests (comment out the first)
# and then the first to do "production" data runs.  

mgmarino

import ROOT
import random
import math
import os
import glob

ROOT.gSystem.Load("libEXOROOT")

def make_wfs():

mgmarino

#include "EXOUtilities/EXOEventData.hh"
#include "EXOUtilities/EXOChargeCluster.hh"
#include "EXOUtilities/EXOScintillationCluster.hh"
#include "EXOUtilities/EXOWaveform.hh"
#include "TFile.h"
#include "TTree.h"
#include "TChain.h"
#include <iostream>
#include <TMath.h>
#include <TH1D.h>

mgmarino

import ROOT

# These are the data that are read in, there should be two versions for SS and
# MS
adc_value = ROOT.RooRealVar("adc_value", "ADC", 500, 5000) 
sigma = ROOT.RooRealVar("sigma", "sigma", 0, 10) 
true_energy = ROOT.RooRealVar("true_energy", "true_energy", 0, 5000) 

# These guys need to be close to the correct values otherwise Minuit has a hard
# time

mgmarino

    SS_nll = ROOT.RooNLLVar("SS_nll", "SS_nll", modelSS, dataSS, 
                            ROOT.RooFit.Extended(), ROOT.RooFit.NumCPU(cpus))
    MS_nll = ROOT.RooNLLVar("MS_nll", "MS_nll", modelMS, dataMS, 
                            ROOT.RooFit.Extended(), ROOT.RooFit.NumCPU(cpus))
    nll_add = []
    for floater_one, var_one, err_one, float_two, var_two, err_two, rho in [ 
        (OffsetSS, offsetSSMean, offsetSSMeanSigma,                                     
         SlopeSS, slopeSSMean, slopeSSMeanSigma, SSCovarianceTerm),
        (OffsetMS, offsetMSMean, offsetMSMeanSigma,
         SlopeMS, slopeMSMean, slopeMSMeanSigma, MSCovarianceTerm),

mgmarino

#!/usr/bin/python
import ROOT

ROOT.gSystem.Load("libEXOUtilities")
f = ROOT.TFile("/media/750GB/FiberLight/run00002225-000.root")
t = f.Get("tree")
A = t.GetUserInfo().First()

while A:
    #A.Print()

mgmarino

import pyWIMP.utilities.utilities as util 
import glob
import os
import sys

def do_child_work(afile, force):
    output_filename = "out_" + os.path.basename(afile)
    if os.path.exists(output_filename) and not force: 
        print "Refusing to overwrite: ", output_filename
        return