# Generate some data, fill the static histograms and NTuple.
for i in range(10000):
x = gauss( mean, sigma )
sthist.addValues ( (x, ) )
ntuple.addRow ( (x, ) )
if i < 1000 :
sthists.addValues ( (x, ) ) # only fill with first 1000
# Print some statistics from static histogram
# Could do same for dynamic
datarep = sthist.getDataRep()
print "Histogram :"
print " Title : " + sthist.getTitle()
print " Entries : %i" % sthist.numberOfEntries()
print " Mean = %f" % datarep.getMean ( 'x' )
print " Rms = %f" % datarep.getRMS ( 'x' )
# Print the average X value on the display
canvas.selectDisplay ( sthist)
canvas.addTextRep ( sthist, 'averagex' )
canvas.selectDisplay ( dyhist)
canvas.addTextRep ( dyhist, 'averagex' )
# Get the contents of the bins as a DataArray
high = sthist.createDataArray()
low = sthists.createDataArray ()
# Take difference with high statistics one scaled down, and a column
# to the low one.
#
# If GemConditionsWord == 7., then fill element with `nbrTkrTriggered',
# otherwise with -1
#
daSvac [ label_TowerTkrTrigGemCond ] = \
numarray.choose ( t, ( nbrTkrTriggered, -1 ) )
daSvac [ label_TowerCalLeTrigGemCond ] = \
numarray.choose ( t, ( nbrCalLeTriggered, -1 ) )
tend = time.time()
print "Took %f seconds create the 4 new columns with 500,000 rows each" % \
(tend -tstart)
#
# The rest is standard procedure
#
tkrtrighist = Display ("Histogram", daSvac, (label_TkrTriggered,) )
canvas.addDisplay ( tkrtrighist )
tkrtrighist.setLog ( 'y', True)
calletrighist = Display ("Histogram", daSvac, (label_CalLeTriggered, ) )
canvas.addDisplay ( calletrighist )
calletrighist.setLog ( 'y', True)
tkrtrighist_gemcond = Display ( "Histogram", daSvac,
(label_TowerTkrTrigGemCond, ) )
canvas.addDisplay ( tkrtrighist_gemcond )
tkrtrighist_gemcond.setRange ( 'x', 0, 16)
tkrtrighist_gemcond.setLog ( 'y', True)
print tkrtrighist_gemcond.numberOfEntries()
