ScifiCTR.Scifi_TimeOfTracks Class Reference

Inheritance diagram for ScifiCTR.Scifi_TimeOfTracks:

Collaboration diagram for ScifiCTR.Scifi_TimeOfTracks:

Public Member Functions
	Init (self, monitor)
	ExecuteEvent (self, event)
	Plot (self)

Public Attributes
	M
	trackTask
	V
	tdcScifiStationCalib

Detailed Description

Definition at line 426 of file ScifiCTR.py.

Member Function Documentation

ExecuteEvent()

ScifiCTR.Scifi_TimeOfTracks.ExecuteEvent	(		self,
			event
	)

Definition at line 438 of file ScifiCTR.py.

   def ExecuteEvent(self,event):
       h = self.M.h
       for aTrack in self.M.Reco_MuonTracks:
          if not aTrack.GetUniqueID()==1: continue
          fitStatus = aTrack.getFitStatus()
          if not fitStatus.isFitConverged(): continue
          state = aTrack.getFittedState()
          mom = state.getMom()
          slopeX = mom.X()/mom.Z()
          slopeY = mom.Y()/mom.Z()
          sortedClusters={}
          DetID2Key={}
          for nHit in range(event.Digi_ScifiHits.GetEntries()):
            DetID2Key[event.Digi_ScifiHits[nHit].GetDetectorID()] = nHit
#
          pos = {}
          for s in range(1, nStations+1):  sortedClusters[s] = {'H':[],'V':[]}
          for nM in range(aTrack.getNumPointsWithMeasurement()):
              state = aTrack.getFittedState(nM)
              Meas = aTrack.getPointWithMeasurement(nM)
              W      = Meas.getRawMeasurement()
              clkey = W.getHitId()
              aCl   = self.trackTask.clusScifi[clkey]
              aHit = event.Digi_ScifiHits[DetID2Key[aCl.GetFirst()]]
              s = aCl.GetFirst()//1000000
              aCl.GetPosition(A,B)
              mat = (aCl.GetFirst()//10000)%10
 
              if aHit.isVertical(): 
                        L = B[1]-state.getPos()[1]
                        sortedClusters[s]['V'].append( [clkey,L,B[0],state.getPos()[1],mat,(A[2]+B[2])/2.,state.getPos()[0],aCl.GetFirst(),aCl.GetTime()] )
              else:  
                        L = A[0]-state.getPos()[0]
                        sortedClusters[s]['H'].append( [clkey,L,A[1],state.getPos()[0],mat,(A[2]+B[2])/2.,state.getPos()[1],aCl.GetFirst(),aCl.GetTime()] )
          stationTimes = {}
          for s in range(1, nStations+1):
             if not (len(sortedClusters[s]['V']) * len(sortedClusters[s]['H']) ) ==1: continue
             sTime = 0
             for proj in ['V','H']:
                clkey = sortedClusters[s][proj][0][0]
                aCl    = self.trackTask.clusScifi[clkey]
                L = sortedClusters[s][proj][0][1]
                time =  aCl.GetTime()   # Get time in ns, use fastest TDC of cluster
                mat  =  sortedClusters[s][proj][0][4]
                sipm_array = int(sortedClusters[s][proj][0][7]/1000)%10
                sipm_channel = sortedClusters[s][proj][0][7]%1000
                if channelTimeAlignment==1 and s==2 and proj=='H':# apply the corrections per channel
                   time-=  self.tdcScifiStationCalib[s][1][proj][1][sipm_array][sipm_channel]
                time-=  self.tdcScifiStationCalib[s][1][proj][mat]  # correct as function of station / projection / mat
                time-=  self.tdcScifiStationCalib[s][0]                  # internal station calibration
                time-=  abs(L)/self.V
 
# cross check 
               # corTime = self.M.Scifi.GetCorrectedTime(aCl.GetFirst(),aCl.GetTime(),abs(L) )
               # print('test',aCl.GetTime(),time,corTime)
 
                sTime += time
             stationTimes[s] = [sTime/2.,(sortedClusters[s]['H'][0][5] + sortedClusters[s]['V'][0][5])/2.]
          # require station 1 to be present which defines T0
          if not 1 in stationTimes: continue
 
          T0 = stationTimes[1][0]
          Z0 = stationTimes[1][1]
          for s in stationTimes:
             if s == 1: continue
             dZ = stationTimes[s][1]-Z0
             dT = stationTimes[s][0]-T0
             dL = dZ * ROOT.TMath.Sqrt( slopeX**2+slopeY**2+1 )
             # if slopeY>0.1:  dL = -dL     # cosmics from the back
             # dT -= dL / u.speedOfLight
             rc = h['dTvsZ'].Fill(dL,dT)
             if abs(slopeX)<0.1 and abs(slopeY)<0.1:
                  rc = h['dTvsZ_beam'].Fill(dL,dT)
 
 
 

Init()

ScifiCTR.Scifi_TimeOfTracks.Init	(		self,
			monitor
	)

Definition at line 428 of file ScifiCTR.py.

   def Init(self,monitor):
       self.M = monitor
       self.trackTask = self.M.FairTasks['simpleTracking']
       h = self.M.h
       self.V = M.Scifi.GetConfParF("Scifi/signalSpeed")
       ut.bookHist(M.h,'dTvsZ','dT versus dL; dL/0.5cm [cm];dT/100ps [ns]',140,0.,70.,120,-6.,6.)
       ut.bookHist(M.h,'dTvsZ_beam','dT versus dL, beam dL/0.5cm [cm];dT/100ps [ns]',140,0.,70.,120,-6.,6.)
       with open('ScifiTimeAlignment_v2', 'rb') as fh:
           self.tdcScifiStationCalib = pickle.load(fh)
 

Plot()

ScifiCTR.Scifi_TimeOfTracks.Plot

(

self

)

Definition at line 514 of file ScifiCTR.py.

   def Plot(self):
      h = self.M.h
      ut.bookCanvas(h,'Tscifi','',1800,1200,2,1)
      h['Tscifi'].cd(1)
      h['dTvsZ'].Draw('colz')
      h['Tscifi'].cd(2)
      h['dTvsZ_beam'].Draw('colz')
 

Member Data Documentation

M

ScifiCTR.Scifi_TimeOfTracks.M

Definition at line 429 of file ScifiCTR.py.

tdcScifiStationCalib

ScifiCTR.Scifi_TimeOfTracks.tdcScifiStationCalib

Definition at line 436 of file ScifiCTR.py.

trackTask

ScifiCTR.Scifi_TimeOfTracks.trackTask

Definition at line 430 of file ScifiCTR.py.

V

ScifiCTR.Scifi_TimeOfTracks.V

Definition at line 432 of file ScifiCTR.py.

---

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

• shipLHC/scripts/ScifiCTR.py