vetoTimeCalibration.vetoTDCchannelCalibration Class Reference

Inheritance diagram for vetoTimeCalibration.vetoTDCchannelCalibration:

Collaboration diagram for vetoTimeCalibration.vetoTDCchannelCalibration:

Public Member Functions
	Init (self, options, monitor)
	ExecuteEvent (self, event)
	Finalize (self)

Public Attributes
	M
	options
	xCheck
	tag
	tdcVetoCalib

Detailed Description

Definition at line 216 of file vetoTimeCalibration.py.

Member Function Documentation

ExecuteEvent()

vetoTimeCalibration.vetoTDCchannelCalibration.ExecuteEvent	(		self,
			event
	)

Definition at line 238 of file vetoTimeCalibration.py.

   def ExecuteEvent(self,event):
       h = self.M.h
       s = 1
       vetoHits = {10:[],11:[]}
       for aHit in event.Digi_MuFilterHits:
          if not aHit.isValid(): continue
          detID = aHit.GetDetectorID()//1000
          if not detID<20: continue
          vetoHits[detID].append(aHit)
       if not(len(vetoHits[10])==1) and not(len(vetoHits[11])==1): return
 
       for aTrack in event.Reco_MuonTracks:
          if not aTrack.GetUniqueID()==1: continue
          tdc = {10:{},11:{}}
          S = aTrack.getFitStatus()
          if S.isFitConverged(): continue
          mom    = aTrack.getFittedState().getMom()
          pos      = aTrack.getFittedState().getPos()
 
          for l in range(2):
             zEx = self.M.zPos['MuFilter'][1*10+l]
             lam      = (zEx-pos.z())/mom.z()
             yEx        = pos.y()+lam*mom.y()
             xEx        = pos.x()+lam*mom.x()  # needed for correction of signal propagation
             for aHit in vetoHits[10+l]:
                detID = aHit.GetDetectorID()
                self.M.MuFilter.GetPosition(detID,A,B)
                D = (A[1]+B[1])/2. - yEx
                if abs(D)<5:
                   bar = (detID%1000)
                   tdc[10+l][bar] = {'L':{},'R':{}}
                   for k in range(16):
                      qdc = aHit.GetSignal(k)
                      if qdc <0: continue
                      kx = k
                      side = 'L'
                      if not k < 8 : 
                           kx = k - 8
                           side = 'R'
                      tdc[10+l][bar][side][kx] = aHit.GetTime(k)
                      if self.xCheck:
                          cor = 0
                          if kx in self.tdcVetoCalib[s*10+l][side][bar]:
                              cor = self.tdcVetoCalib[s*10+l][side][bar][kx][0]
                              tdc[10+l][bar][side][kx] -= cor
      # make internal alignment
          for l in range(2):
            for bar in tdc[10+l]:
             for side in tdc[10+l][bar]:
                if not refChannel in tdc[10+l][bar][side]: continue
                t0 = tdc[10+l][bar][side][refChannel]
                for k in tdc[10+l][bar][side]:
                   dt = tdc[10+l][bar][side][k]-t0
                   key = self.M.sdict[s]+str(10+l)+'_'+str(bar)+side+str(k)+self.tag
                   rc = h['dtChan_'+key].Fill(dt)
 

Finalize()

vetoTimeCalibration.vetoTDCchannelCalibration.Finalize

(

self

)

Definition at line 294 of file vetoTimeCalibration.py.

   def Finalize(self):
       h = self.M.h
       s = 1
       h['tdcCalib'] = {10:{'L':{},'R':{}},11:{'L':{},'R':{}}}
       ut.bookCanvas(h,'dummy','',900,600,1,1)
       h['dummy'].cd()
       for l in range(self.M.systemAndPlanes[s]):
          for side in ['L','R']:
             tname = 'v'+str(l)+side+self.tag
             ut.bookCanvas(h,tname,tname,2048,2048,8,8)
             for bar in range(self.M.systemAndBars[s]):
                 h['tdcCalib'][s*10+l][side][bar] = {}
                 for k in range(0,8):
                    h[tname].cd(bar*8+k+1)
                    key = self.M.sdict[s]+str(10+l)+'_'+str(bar)+side+str(k)+self.tag
                    dt = h['dtChan_'+key]
                    if not k == refChannel:  
                        rc = dt.Fit('gaus','SQ')
                        result = rc.Get()
                        if result:    h['tdcCalib'][s*10+l][side][bar][k]=[result.Parameter(1),result.Parameter(2)]
                    h['dtChan_'+key].Draw()
 
       ut.bookHist(h,'TDCshiftsMean'+self.tag,'TDCshiftsMean '+self.tag,100,-0.2,0.2)
       ut.bookHist(h,'TDCshiftsSigma'+self.tag,'TDCshiftsSigma '+self.tag,100,0.,0.25)
       for l in h['tdcCalib']:
           for side in h['tdcCalib'][l]:
              for bar in h['tdcCalib'][l][side]:
                 for k in h['tdcCalib'][l][side][bar]:
                     if k==refChannel: continue
                     rc = h['TDCshiftsMean'+self.tag].Fill(h['tdcCalib'][l][side][bar][k][0])
                     rc = h['TDCshiftsSigma'+self.tag].Fill(h['tdcCalib'][l][side][bar][k][1])
       ut.bookCanvas(h,'ChannelSummary'+self.tag,'Channel Summary',1200,600,2,1)
       tc = h['ChannelSummary'+self.tag].cd(1)
       h['TDCshiftsMean'+self.tag].Draw()
       tc = h['ChannelSummary'+self.tag].cd(2)
       h['TDCshiftsSigma'+self.tag].Draw()
 
       if not self.xCheck:
           with open('tdcVetoInternalCalibration', 'wb') as fh:
               pickle.dump(h['tdcCalib'], fh)
           ut.writeHists(h,'tdcCalib.root')
       else:
           ut.writeHists(h,'tdcCalibCor.root')
 

Init()

vetoTimeCalibration.vetoTDCchannelCalibration.Init	(		self,
			options,
			monitor
	)

Definition at line 218 of file vetoTimeCalibration.py.

   def Init(self,options,monitor):
       self.M = monitor
       self.options = options
       self.xCheck = False
       self.tag = ""
       if options.xCheck : 
           self.xCheck = True
           self.tag = "X"
       h = self.M.h
       s = 1
       for l in range(self.M.systemAndPlanes[s]):
          for bar in range(self.M.systemAndBars[s]):
             for side in ['L','R']:
                for k in range(8):
                   key = self.M.sdict[s]+str(s*10+l)+'_'+str(bar)+side+str(k)+self.tag
                   ut.bookHist(h,'dtChan_'+key,'dt '+key+";dt [timestamp]",200,-1.,1.)
       if self.xCheck:
            with open('tdcVetoInternalCalibration', 'rb') as fh:
               self.tdcVetoCalib = pickle.load(fh)
 

Member Data Documentation

M

vetoTimeCalibration.vetoTDCchannelCalibration.M

Definition at line 219 of file vetoTimeCalibration.py.

options

vetoTimeCalibration.vetoTDCchannelCalibration.options

Definition at line 220 of file vetoTimeCalibration.py.

tag

vetoTimeCalibration.vetoTDCchannelCalibration.tag

Definition at line 222 of file vetoTimeCalibration.py.

tdcVetoCalib

vetoTimeCalibration.vetoTDCchannelCalibration.tdcVetoCalib

Definition at line 236 of file vetoTimeCalibration.py.

xCheck

vetoTimeCalibration.vetoTDCchannelCalibration.xCheck

Definition at line 221 of file vetoTimeCalibration.py.

---

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

• shipLHC/scripts/vetoTimeCalibration.py