Changeset 14175 for fact/tools/pyscripts


Ignore:
Timestamp:
06/13/12 16:46:02 (12 years ago)
Author:
neise
Message:
can now create trees
File:
1 edited

Legend:

Unmodified
Added
Removed
  • fact/tools/pyscripts/sandbox/dneise/spikes/ana.py

    r14171 r14175  
    110110
    111111event = run.next()
    112 bt_old = event['board_times'].copy
     112bt_old = event['board_times'].copy()
    113113
    114114for event in run:
    115115    data = event['data']
    116     s, d = mars_spikes(data)
     116    singles, doubles = mars_spikes(data)
    117117    sc = event['start_cells']
    118     bt = event['board_times'].copy()
     118    board_time = event['board_times'].copy()
    119119   
    120     if len(s) >0 :
    121         chid.value  = s[0][0]
    122         startcell.value = sc[ chid.value ]
    123         number_of_singles.value = len(s)
    124         time.value = bt[s[0][0]/9]-bt_old[s[0][0]/9]
    125         for i in s:
    126             log = i[1]
    127             phys  = (startcell.value+log)%1024
    128             position_of_spikes_in_logical_pipeline[i] = log
    129             position_of_spikes_in_physical_pipeline[i] = phys
    130         baum.Fill()
    131    
     120    ss = []
     121    for i in range(1440):
     122        ss.append([])
     123    for s in singles:
     124        ss[s[0]].append(s[1])
     125       
     126    for pixel_id,pixel in enumerate(ss):
     127        if len(pixel) >0 :
     128            chid.value  = pixel_id
     129            startcell.value = sc[ chid.value ]
     130            number_of_singles.value = len(pixel)
     131            time.value = board_time[pixel_id/36]-bt_old[pixel_id/36]
     132            for num,s in enumerate(pixel):
     133                log = s
     134                phys  = (startcell.value+log)%1024
     135                position_of_spikes_in_logical_pipeline[num] = log
     136                position_of_spikes_in_physical_pipeline[num] = phys
     137            baum.Fill()
     138
    132139    bt_old = event['board_times'].copy()
    133140
Note: See TracChangeset for help on using the changeset viewer.