source: fact/tools/pyscripts/sandbox/vogler/DRS_End_Noise_1.py

Last change on this file was 14173, checked in by vogler, 12 years ago
inital filling of my sandbox
  • Property svn:executable set to *
File size: 4.0 KB
Line 
1#!/usr/bin/python -tt
2# ********************************
3# Test script for the CalFits class
4#
5# written by Thomas Kraehenbuehl, ETH Zurich
6# tpk@phys.ethz.ch, +41 44 633 3973
7# April 2012
8# ********************************
9#
10# modified and adapted py Patrick Vogler
11#
12# ################################
13from ROOT import gSystem
14gSystem.Load("calfactfits_h.so") # according to new naming scheme
15from ROOT import *
16
17
18
19
20
21
22
23#define filenames
24
25# 2011 12 05
26calibfilename = '/fact/raw/2011/12/05/20111205_011.drs.fits.gz' # NROI 300, pedestal
27
28datafilename = '/fact/raw/2011/12/05/20111205_022.fits.gz' # NROI 300, 5 minutes physics data with
29
30
31import numpy as np
32from scipy import weave
33from scipy.weave import converters
34
35from plotters import Plotter # ADC display
36# from plotters import CamPlotter # event display
37#from drs_spikes import DRSSpikes
38
39
40print "Testing object creation: "
41caltest = CalFactFits(datafilename, calibfilename)
42npcalevent = np.empty( caltest.npix * caltest.nroi, np.float64) #.reshape(caltest.npix ,caltest.nroi)
43caltest.SetNpcaldataPtr(npcalevent)
44
45
46numroi = np.int64(caltest.nroi)
47numpix = np.int64(caltest.npix)
48numevents = np.int64(caltest.nevents)
49
50num_LP_ev = 0 # number of ext Lightpulser events
51num_ped = 0 # number of pedestal events
52
53
54print "Common variables run information: "
55print "ROI: ", numroi
56print "#Pix: ", numpix
57print "Number of events: ", numevents
58print
59
60
61event = np.zeros((numpix, numroi)) # create an array to store an event in the format numpix * numroi (2-dim array)
62print "Array event created"
63run_average = np.zeros((numpix, numroi)) # create an array to store the "average event" of a run
64data_average_run = np.zeros(numroi)
65print "Arrays created "
66
67
68print "... looping..."
69
70while caltest.GetCalEvent(): # Loop ueber alle events
71 # print npcalevent ## Daten, Array 1 dim Laenge caltest.npix * caltest.nroi 1 Event
72
73
74 if (caltest.event_triggertype == 4): #ext LP event
75 print 'event id:', caltest.event_id, ' Trigger Type:', caltest.event_triggertype
76 num_LP_ev = num_LP_ev + 1
77 event = np.reshape(npcalevent, (numpix, -1)) # bring the event the shape numpix * numroi (2-dim array)
78 run_average += event
79
80print "Looped ... "
81
82
83
84
85#for i in range (0, (numpix - 1) ):
86# data_average_run += run_average[i]
87#print "Looped second loop "
88
89
90del caltest
91print "caltest deleted "
92
93
94
95print "Common variables run information: "
96print "ROI: ", numroi
97print "#Pix: ", numpix
98print "Number of events: ", numevents
99print "Number of external Lightpulser events: ", num_LP_ev
100print "Number of pedestal: ", num_ped
101print
102
103
104
105#####################################################################################################################
106print "creating plotters ..."
107
108#make a Plotter class ... this is an easy way for plotting ... but there are
109# many was to plot data ... without this class ... it was written for convenience, but there is no strong reason to use it...
110#myplotter = Plotter('titel of the plot', xlabel='time in slices', ylabel='amplitude calibrated data ... in mV')
111
112myplotter = Plotter('titel of the plot', xlabel='time in slices', ylabel='amplitude calibrated data ... in mV')
113myplotter2 = Plotter('titel of the plot', xlabel='time in slices', ylabel='amplitude calibrated data ... in mV')
114myplotter3 = Plotter('titel of the plot', xlabel='time in slices', ylabel='amplitude calibrated data ... in mV')
115
116myplotter4 = Plotter('titel of the plot', xlabel='time in slices', ylabel='amplitude calibrated data ... in mV')
117
118print "Plotters created "
119###################################################################################################################
120
121myplotter((run_average[20]) / num_LP_ev, 'pix 20, average' )
122
123myplotter2((run_average[21]) / num_LP_ev, 'pix 21, average' )
124
125myplotter3((run_average[22]) / num_LP_ev, 'pix 22, average' )
126
127myplotter4( data_average_run / (numpix * num_LP_ev) , 'average signal (all pixel) over the whole run' )
128
129answer = raw_input('type "quit" to quit ')
130while not 'quit' in answer:
131 answer = raw_input('type "quit" to quit ')
132
Note: See TracBrowser for help on using the repository browser.