Esempi in Python per CalibrateMods

Esempio n. 1

File: findPeaks.py Progetto: Ak1508/PyROOT-scripts

import pickle

usage = 'Usage: ./findPeaks DataSetName \n'

if len(sys.argv) < 2:
  print usage
  sys.exit()


#ROOT.TROOT.gApplication.ExecuteFile("~/rootlogon.C")

setName=sys.argv[1]
parentDirectory = '/Users/rob/Development/AnalysisToolkit/WIPPnRuns/CuShield/'

file = ROOT.TFile(setName+'Calibrated.root')
tree = file.Get('T')
tree.Draw('fEnergyCal>>h1(5000,0,5000)')
h1 = ROOT.gDirectory.Get('h1')
#h1.Sumw2()
h1.Draw("hist")

peakInfo = []

peakInfo = CalibrateMods.fitPeaks(h1)

output = open('Pickled'+setName+'Peaks.txt','w')
pickle.dump(peakInfo, output)
output.close()

Esempio n. 2

File: plotSigma.py Progetto: robj137/PyROOT-scripts

#!/usr/bin/python

import ROOT
import CalibrateMods
import sys
import os
import pickle

usage = 'Usage: '+str(sys.argv[0])+ 'DataSetName \n'

if len(sys.argv) < 2:
  print usage
  sys.exit()


#ROOT.TROOT.gApplication.ExecuteFile("~/rootlogon.C")

setName=sys.argv[1]
parentDirectory = '/Users/rob/Development/AnalysisToolkit/WIPPnRuns/CuShield/'

f = open('PickledCalibration'+setName+'.txt')
peakInfo = pickle.load(f)
CalibrateMods.resolvePoints(peakInfo, setName)

Esempio n. 3

File: calibrate.py Progetto: Ak1508/PyROOT-scripts

if len(sys.argv) < 2:
    print usage
    sys.exit()

    ratio = 1.0
if len(sys.argv) > 3:
    ratio = float(sys.argv[3])

#ROOT.TROOT.gApplication.ExecuteFile("~/rootlogon.C")

setName = sys.argv[1]
parentDirectory = '/Users/rob/Development/AnalysisToolkit/WIPPnRuns/CuShield/' + setName + '/'
mcaDirectory = parentDirectory + 'MCA/'
print mcaDirectory
if (sys.argv[2] == 'find'):
    filename = CalibrateMods.createMCAFile(setName, mcaDirectory)

    histo = CalibrateMods.getMCASumHist(filename, setName)
    peaks = [
        129.065, 209.253, 238.632, 270.245, 277.371, 300.087, 328.000, 338.320,
        463.004, 583.187, 609.32, 661.657, 727.330, 794.947, 860.557, 911.204,
        1173.226, 1332.492, 1460.822, 1588.20, 1592.511, 1764.491, 2103.511,
        2614.511
    ]
    peakInfo = []
    for peak in peaks:
        peakFit = CalibrateMods.fitPeak(histo, peak, ratio)
        if len(peakFit) > 0:
            peakFit.append(peak)
            peakInfo.append(peakFit)

Esempio n. 4

File: calibrate.py Progetto: robj137/PyROOT-scripts

if len(sys.argv) < 2:
  print usage
  sys.exit()

  ratio = 1.0
if len(sys.argv) >3:
  ratio = float(sys.argv[3])

#ROOT.TROOT.gApplication.ExecuteFile("~/rootlogon.C")

setName=sys.argv[1]
parentDirectory = '/Users/rob/Development/AnalysisToolkit/WIPPnRuns/CuShield/'+setName+'/'
mcaDirectory = parentDirectory+'MCA/'
print mcaDirectory
if(sys.argv[2] == 'find'):
  filename = CalibrateMods.createMCAFile(setName,mcaDirectory)

  histo = CalibrateMods.getMCASumHist(filename, setName)
  peaks = [ 129.065, 209.253, 238.632, 270.245, 277.371, 300.087, 328.000, 
            338.320, 463.004, 583.187, 609.32, 661.657, 727.330, 
            794.947, 860.557, 911.204, 1173.226, 1332.492, 1460.822, 
            1588.20, 1592.511, 1764.491, 2103.511, 2614.511]
  peakInfo = []
  for peak in peaks:
    peakFit = CalibrateMods.fitPeak(histo, peak, ratio)
    if len(peakFit) > 0:
      peakFit.append(peak)
      peakInfo.append(peakFit)

  output = open('PickledCalibration'+setName+'.txt','w')
  pickle.dump(peakInfo, output)

Esempio n. 5

File: plotSigma.py Progetto: Ak1508/PyROOT-scripts

#!/usr/bin/python

import ROOT
import CalibrateMods
import sys
import os
import pickle

usage = 'Usage: ' + str(sys.argv[0]) + 'DataSetName \n'

if len(sys.argv) < 2:
    print usage
    sys.exit()

#ROOT.TROOT.gApplication.ExecuteFile("~/rootlogon.C")

setName = sys.argv[1]
parentDirectory = '/Users/rob/Development/AnalysisToolkit/WIPPnRuns/CuShield/'

f = open('PickledCalibration' + setName + '.txt')
peakInfo = pickle.load(f)
CalibrateMods.resolvePoints(peakInfo, setName)