wirecell.sigproc.minisim.Minisim Class Reference

Inheritance diagram for wirecell.sigproc.minisim.Minisim:

Public Member Functions
def	__init__ (self, pib, kwds)
def	__getattr__ (self, key)
def	wire_space (self, hits)
def	apply_block (self, block, iplane, ich, time)
def	clear_frame (self)
def	add_hits (self, hits)

Public Attributes
	pib
	wire_yz_dir
	wire_yz_pit
	wire_yz_off
	frame

Static Public Attributes
	defaults

Private Attributes
	_cfg

Detailed Description

Simulate the response to a distribution of charge.

todo:

    - for each hit, find all rf's for the hits impact, use each
      region number as offset to hit's wire number, fill frame.

Definition at line 17 of file minisim.py.

Constructor & Destructor Documentation

def wirecell.sigproc.minisim.Minisim.__init__	(		self,
			pib,
			kwds
	)

Create a mini simulation.

- pib :: a response.PlaneImpactBlocks object

Definition at line 35 of file minisim.py.

    def __init__(self, pib, **kwds):
        '''
        Create a mini simulation.

        - pib :: a response.PlaneImpactBlocks object

        '''
        self.pib = pib
        self._cfg = dict(self.defaults, **kwds)

        self.wire_yz_dir = numpy.asarray((dir_yz(30), dir_yz(150), dir_yz(90)))
        self.wire_yz_pit = numpy.asarray((dir_yz(-60), dir_yz(60), dir_yz(0)))
        self.wire_yz_off = numpy.asarray(((0.0, 0.0),
                                          (0.0, 0.0),
                                          (0.0, 0.5*self.pitch)))

        self.frame = None       # becomes 3-tuple
        self.clear_frame()      # prime

Member Function Documentation

def wirecell.sigproc.minisim.Minisim.__getattr__	(		self,
			key
	)

Definition at line 54 of file minisim.py.

    def __getattr__(self, key):
        return self._cfg[key]

def wirecell.sigproc.minisim.Minisim.add_hits	(		self,
			hits
	)

Add response to hits to current frame.

Definition at line 110 of file minisim.py.

    def add_hits(self, hits):
        '''
        Add response to hits to current frame.
        '''
        nearest_chimps = self.wire_space(hits)
        for iplane, chimps in enumerate(nearest_chimps):
            for ihit, (ch,imp) in enumerate(chimps):
                x,y,z,q,t = hit = hits[ihit]
                block = self.pib.region_block("uvw"[iplane], imp)
                self.apply_block(block, iplane, ch, (t-self.pib.tmin) + (x-self.pib.xstart)/self.velocity)
            
            
            

def wirecell.sigproc.minisim.Minisim.apply_block	(		self,
			block,
			iplane,
			ich,
			time
	)

Apply block array to plane array with ich and itick offsets

Definition at line 76 of file minisim.py.

    def apply_block(self, block, iplane, ich, time):
        '''
        Apply block array to plane array with ich and itick offsets
        '''
        dch = block.shape[0]//2
        chm = int(max(ich-dch, 0))
        chp = int(min(ich+dch, self.nwires[iplane]-1))
        #print 'ch:',dch, chm, chp

        tfinebin_min = int(time/self.pib.tbin)
        tfinebin_jmp = int(self.tick/self.pib.tbin)
        #print 'time:',time,tfinebin_min, tfinebin_jmp

        sampled = block[:,tfinebin_min::tfinebin_jmp]
        #print 'sampled:',sampled.shape

        tsampbin_num = len(sampled[0])
        tsampbin_min = int(time/self.tick)
        tsampbin_max = min(tsampbin_min+tsampbin_num, len(self.frame[iplane][0]))
        print 'samp bin:',tsampbin_num, tsampbin_min, tsampbin_max*self.tick/units.us

        self.frame[iplane][chm:chp+1, tsampbin_min:tsampbin_max] = sampled

        

def wirecell.sigproc.minisim.Minisim.clear_frame

(

self

)

Initialize a new frame.

Definition at line 100 of file minisim.py.

    def clear_frame(self):
        '''
        Initialize a new frame.
        '''
        frame = list()
        for nwires in self.nwires:
            frame.append(numpy.zeros((nwires, self.nticks)))
        self.frame = tuple(frame)
        return

def wirecell.sigproc.minisim.Minisim.wire_space	(		self,
			hits
	)

Return the location of a hit in terms of (wire number, impact).

Definition at line 57 of file minisim.py.

    def wire_space(self, hits):
        '''
        Return the location of a hit in terms of (wire number, impact).
        '''
        ret = list()

        hit_yz = hits[:,1:3]
        for iplane in range(3):
            yz = hit_yz - self.wire_yz_off[iplane]
            pit = numpy.dot(yz, self.wire_yz_pit[iplane]) # distance along pitch from center
            pwidth = self.pitch * self.nwires[iplane]/2
            pit += pwidth       # distance along pitch from wire 0
            wiredist = pit/self.pitch
            wire = numpy.asarray(numpy.round(wiredist), dtype=int)
            impact = numpy.round((wiredist-wire)/self.impact) * self.impact
            ret.append(numpy.vstack((wire, impact)).T)
        return ret
        

Member Data Documentation

wirecell.sigproc.minisim.Minisim._cfg

private

Definition at line 43 of file minisim.py.

wirecell.sigproc.minisim.Minisim.defaults

static

Initial value:

= dict(
        nticks=9600,
        nwires=(2400,2400,3456),
        pitch=3.0*units.mm,
        impact=0.3*units.mm,
        velocity=1.6*units.mm/units.us,
        tick=0.5*units.us,      # digitization tick
        )

Definition at line 27 of file minisim.py.

wirecell.sigproc.minisim.Minisim.frame

Definition at line 51 of file minisim.py.

wirecell.sigproc.minisim.Minisim.pib

Definition at line 42 of file minisim.py.

wirecell.sigproc.minisim.Minisim.wire_yz_dir

Definition at line 45 of file minisim.py.

wirecell.sigproc.minisim.Minisim.wire_yz_off

Definition at line 47 of file minisim.py.

wirecell.sigproc.minisim.Minisim.wire_yz_pit

Definition at line 46 of file minisim.py.

---

The documentation for this class was generated from the following file:

• wire-cell-build/python/wirecell/sigproc/minisim.py

---

Generated by   1.8.11