# 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()