AnalysisTest.py

Go to the documentation of this file.

#--Initialize ------------------------------------------------
execfile('../options/AnalysisInit.py')
#-------------------------------------------------------------

#-------------------------------------------------------------
def electron(c):
#-------------------------------------------------------------
    if c : return 1
    return 0

#-------------------------------------------------------------
def initialize():
#-------------------------------------------------------------
    global h0, h1, h2
    try:
        h0 = his.histo('h0')
        h1 = his.histo('h1')
        h2 = his.histo('h2')
    except:
        h0 = his.book('h0', 'Histogram 1D for random tests', 35, -10., 10.) 
        h1 = his.book('h1', 'Histogram 1D for tests', 20, 0., 40.) 
        h2 = his.book('h2', 'Histogram 2D for test2', 20, 0., 40., 20, 0., 10.) 

#-------------------------------------------------------------
def execute(): 
#-------------------------------------------------------------
  global cands
    cands = evt['Anal/AxParticleCandidates']

    electrons = select(cands,electron)
    print 'MyTracks collection contains ' + `len(tks)` + ' tracks and ' + `len(electrons)` + ' electrons'


#-------End of declarative part-----------------------------------------------------

#g.run(10)
#excel.plot(h1)

g.run(1)
sel = evt['Phys/PhysSel']
sel.decayIsInMCTree('B0PP')