MufiCTR Namespace Reference

Functions
	start ()
	start0 (plist=['tmp4705p9'], s=3)
	xCheck ()
	execute (s=3)
	timingOfEvent (makeCluster=False, debug=False)

Variables
	parser = ArgumentParser()
	dest
	help
	default
	options = parser.parse_args()
	f = ROOT.TFile(options.file)
dict	h = {}
dict	S = {1:[1800,800,2,1],2:[1800,1500,2,3],3:[1800,1800,2,4]}
dict	channelsPerPlane = {1:7*16,3:60*2}
dict	Planes = {1:2,3:7}
	runN = options.file[options.file.find(".root")-6:options.file.find(".root")]
dict	par0 = {}
dict	alignTpar = {}
dict	alignTparB2 = {}

Function Documentation

execute()

MufiCTR.execute

(

s = 3

)

Definition at line 137 of file MufiCTR.py.

def execute(s=3):
  detector = 'mufi-'
  system = 'DS'
  if s==1: system = 'Veto'
  ut.bookCanvas(h,'c1','',640,480,1,1)
  h['c1'].cd()
  fitLimits ={1:[-15,-5],11:[-15,-5],3:[-10,-2],13:[-25,-15]}
  for l in range(S[s][2]*S[s][3]):
    if s==3 and l>6: continue
    tag = str(s*10+l)
    hist = h[detector+'dT_'+tag]
    alignTpar[s*10+l] = {}
    alignTparB2[s*10+l] = {}
    for i in range(hist.GetNbinsY()):
        if s==1 and i > 8*16: continue
        tagi = str(i*1000+s*10+l)
        alignTpar[s*10+l][i] = [100,-100]
        h['tmp'+tagi] = hist.ProjectionX('tmp'+tagi,i+1,i+1)
        tmp = h['tmp'+tagi]
        rc = tmp.Fit('gaus','SQ','',fitLimits[s][0],fitLimits[s][1])
        fitres = rc.Get()
        if fitres:
             alignTpar[s*10+l][i]=[fitres.Parameter(1),fitres.Parameter(2)]
        alignTparB2[s*10+l][i] = [100,-100]
        try:
          tmp = h[detector+'dT_'+tag+'B2noB1'].ProjectionX('tmp',i+1,i+1)
          rc = tmp.Fit('gaus','SQ','',fitLimits[s+10][0],fitLimits[s+10][1])
          fitres = rc.Get()
          if fitres:
               alignTparB2[s*10+l][i]=[fitres.Parameter(1),fitres.Parameter(2)]
        except: pass
  nbins = channelsPerPlane[s] * Planes[s]
  ut.bookHist(h,'talign','time alignment B1;channel;offset [ns]',nbins,-0.5,nbins-0.5)
  ut.bookHist(h,'talign2','time alignment B2;channel;offset [ns]',nbins,-0.5,nbins-0.5)
  ut.bookHist(h,'tres','time diff res;channel;sigma [ns]',nbins,-0.5,nbins-0.5)
  ut.bookHist(h,'t12','time diff B1 B2;channel;sigma [ns]',nbins,-0.5,nbins-0.5)
  for l in range(S[s][2]*S[s][3]):
    if s==3 and l>6: continue
    for i in range(hist.GetNbinsY()):
       k = l*channelsPerPlane[s]
       rc = h['talign'].Fill(k+i,alignTpar[s*10+l][i][0])
       rc = h['talign2'].Fill(k+i,alignTparB2[s*10+l][i][0])
       rc = h['tres'].Fill(k+i,abs(alignTpar[s*10+l][i][1]))
       if alignTpar[s*10+l][i][0]<99 and alignTparB2[s*10+l][i][0]<99:
           rc = h['t12'].Fill(k+i,alignTparB2[s*10+l][i][0]-alignTpar[s*10+l][i][0])
       else: rc = h['t12'].Fill(k+i,100)
  if s==3:
     h['talign'].SetMinimum(-10)
     h['talign'].SetMaximum(1)
  else:
     h['talign'].SetMinimum(-15)
     h['talign'].SetMaximum(-5)
  h['talign2'].SetMinimum(-25)
  h['talign2'].SetMaximum(1)
  h['tres'].SetMinimum(0)
  h['tres'].SetMaximum(2.5)
  h['talign'].SetMarkerStyle(30)
  h['talign'].SetMarkerColor(ROOT.kRed)
  h['talign2'].SetMarkerStyle(32)
  h['talign2'].SetMarkerColor(ROOT.kBlue)
  h['tres'].SetMarkerStyle(31)
  h['tres'].SetMarkerColor(ROOT.kGreen)
  h['tres'].SetStats(0)
  h['talign'].SetStats(0)
  h['talign2'].SetStats(0)
  ut.bookCanvas(h,'results',system+' time alignment',2400,1200,1,2)
  tc = h['results'].cd(1)
  h['talign'].Draw('phist')
  tc = h['results'].cd(2)
  h['tres'].Draw('phist')
  h['results'].Print(system+'timeAlign-run'+runN+'.png')
  h['t12'].SetMinimum(-30)
  h['t12'].SetMaximum(0)
  h['t12'].SetStats(0)
  h['t12'].SetMarkerColor(ROOT.kMagenta)
  h['t12'].SetMarkerStyle(53)
 
  ut.bookCanvas(h,'results2','time diff between FW and BW tracks',1200,600,1,1)
  h ['results2'].cd()
  h['t12'].Draw('phist')
  h['results2'].Print(system+'timeAlignFWBW-run'+runN+'.png')
  
  for k in range(0,4):
    h['talignReorder'+str(k)] = h['talign'].Clone('talignReorder'+str(k))
    h['talignReorder'+str(k)].Reset()
  if s==3: 
   slope = 0.0833
   for i in range(h['talign'].GetNbinsX()):
     if i<120 and i%2==0: p=0
     elif i<120 and i%2==1: p=1
     elif i<240 and i%2==0: p=2
     elif i<360 and i%2==0: p=3
     elif i<360 and i%2==1: p=4
     elif i<480 and i%2==0: p=5
     elif i<600 and i%2==0: p=6
     elif i<600 and i%2==1: p=7
     elif i<720 and i%2==0: p=8
     elif i<840 and i%2==0: p=9
     else: continue
     if i%2==0: j=(i//2)%60
     if i%2==1: j=((i-1)//2)%60
     X = h['talign'].GetBinContent(i+1)
     h['talignReorder0'].SetBinContent(j+1+p*60,X)
     h['talignReorder1'].SetBinContent(j+1+p*60,X)
     channel = (j+1+p*60)%60
     if channel<30: X+= channel*slope - 15*slope
     else: X-=(channel-30)*slope - 15*slope
     h['talignReorder2'].SetBinContent(j+1+p*60,X)
 
   for l in range(0,3):
    c = 'talignReorder'+str(l)
    ut.bookCanvas(h,'results'+c,system+' time alignment '+str(l),1200,600,1,1)
    h['results'+c].SetGridy(1)
    h['results'+c].cd()
    par0[l]={}
    for x in h[c].GetListOfFunctions(): x.Delete()
    if l==2:
       par1 = 0
       planes = 20
       iv = 30
       h[c].SetMinimum(-7)
       h[c].SetMaximum(-4)
    else: 
       par1 = slope
       planes = 20
       iv = 30
       h[c].SetMinimum(-20)
       h[c].SetMaximum(3)
    for p in range(planes):
      if p==10:
        imin = 0.+p*iv+7
        imax = 0.+imin+iv-3
      elif p==3:
        imin = 0.+p*iv
        imax = 0.+imin+iv-10
      elif p==6:
        imin = 0.+p*iv+5
        imax = 0.+imin+iv-5
      else:
        imin = 0.+p*iv
        imax = 0.+imin+iv
      theFun =  'fun'+str(l*100+p)
      h[theFun]=ROOT.TF1(theFun,'[0]+x*[1]')
      sign = 1
      if p%2==0: sign = -1
      h[theFun].SetParameter(0, -6)
      if l>0: h[theFun].FixParameter(1, sign*par1)
      rc = h[c].Fit(h[theFun],'SQw+','',imin+2,imax-2)
      res = rc.Get()
      ix = 2
      while res.Chi2() > 1.5:
          ix+=2
          if imax-ix < imin+2: 1/0
          if h[c].GetFunction(theFun): h[c].GetFunction(theFun).Delete()
          if p==3: rc = h[c].Fit(h[theFun],'SQw+','',imin+ix,imax-ix)
          else   : rc = h[c].Fit(h[theFun],'SQw+','',imin+2,imax-ix)  # vertical no statistics
          res = rc.Get()
      par0[l][p]=[res.Parameter(0),res.Parameter(1)]
      print("%i slope %5.3F b = %5.2F"%(p,res.Parameter(1),res.Parameter(0)))
    h[c].GetXaxis().SetRange(1,600)
    h[c].Draw()
    h['results'+c].Print(system+'timeAlignFit'+str(l)+'-run'+runN+'.png')
    
   c = 'talignReorder3'
   ut.bookHist(h,'resT','t residuals',100,-3.,3.)
   for i in range(h[c].GetNbinsX()):
     p = i//30
     if p<20:
         res = h[c.replace('3','2')].GetBinContent(i+1)-par0[2][p][0]
         h[c].SetBinContent(i+1,res)
         rc = h['resT'].Fill(res)
   ut.bookCanvas(h,'results'+c,system+' time alignment 3',1200,600,1,1)
   h['results'+c].SetGridy(1)
   h['results'+c].cd()
   h[c].GetXaxis().SetRange(1,600)
   h[c].SetMinimum(-1)
   h[c].SetMaximum(1)
   h[c].Draw()
   h['results'+c].Print(system+'timeAlignFit3'+'-run'+runN+'.png')
   h['c1'].cd()
   rc = h['resT'].Fit('gaus','S','')
   h['c1'].Update()
   stats = h['resT'].FindObject('stats')
   stats.SetOptFit(1)
   stats.SetOptStat(1000000000)
   stats.SetX1NDC(0.595611)
   stats.SetY1NDC(0.614035)
   stats.SetX2NDC(0.979624)
   stats.SetY2NDC(0.936404)
   h['c1'].Print(system+'timeResiudals'+'-run'+runN+'.png')
      
   h['c1'].cd()
   for I in [ [10.5,89.5],[249.5,328.5],[492.5,571.5] ]:
      p+=1
      rc = h['t12'].Fit('pol1','SQ','',I[0],I[1])
      res = rc.Get()
      print("%i delta t = %5.2F"%(p,res.Parameter(0)))
   m=0
   for i in par0[0]:
    m+=abs(par0[0][i][1])
   print('average slope = ',m/len(par0[0]))
  
# golden nu event
#  $EOSSHIP/eos/experiment/sndlhc/convertedData/physics/2022/run_005056/sndsw_raw-0101.root -g geofile_sndlhc_TI18_V7_22November2022.root
#  849600

start()

MufiCTR.start

(

)

Definition at line 18 of file MufiCTR.py.

def start():
 for x in ['dTScifiDS','dTcorScifiDS']:
  for d in ['','B2noB1','B1only']:
   dd=''
   if d!='': dd = d+'/'
   # safety net if B2noB1 and B1Only plots don't exist
   try:
     h['mufi-'+x+d] = f.mufilter.Get(dd+'mufi-'+x+d).Clone('mufi-'+x+d)
   except:
       try:
         #2024: changed structure of offline monitoring file with shifter and expert subdirs
         h['mufi-'+x+d] = f.mufilter.Get('expert/'+dd+'mufi-'+x+d).Clone('mufi-'+x+d)
       except:
         continue
   h['mufi-'+x+d].Draw()
   h['mufi-'+x+d].Print('mufi-'+x+d+'.png')
   for pad in h['mufi-'+x+d].GetListOfPrimitives():
      for y in pad.GetListOfPrimitives():
        if not y.ClassName().find('TH')<0:
           hname = y.GetName()
           h[hname] = y.Clone(hname)
 

start0()

MufiCTR.start0	(		plist = ['tmp4705p9'],
			s = 3
	)

Definition at line 40 of file MufiCTR.py.

def start0(plist=['tmp4705p9'],s=3):
   detector = 'mufi-'
   wanted = []
   system = 'DS'
   if s==1: system = 'Veto'
   for x in ['','cor']:
    for d in ['','B2noB1','B1only']:
     ut.bookCanvas(h,detector+'dT'+x+'Scifi'+system+d,'dt rel to scifi '+system,S[s][0],S[s][1],S[s][2],S[s][3])
     for l in range(7):
      tag = str(10*s+l)+d
      wanted.append(detector+'dT'+x+'_'+tag)
   for p in plist:
      ut.readHists(h,p,wanted=wanted)
   for x in ['','cor']:
    for d in ['','B2noB1','B1only']:
     for l in range(S[s][2]*S[s][3]):
              if s==3 and l>6: continue
              n = l+1
              tag = str(10*s+l)+d
              if s==3 and n==7: n=8
              # safety net if B2noB1 and B1Only plots don't exist
              try:
                tc = h[detector+'dT'+x+'Scifi'+system+d].cd(n)
                h[detector+'dT'+x+'_'+tag].Draw('colz')
              except:
                 continue
     h[detector+'dT'+x+'Scifi'+system+d].Print('mufi-dT'+x+'Scifi'+system+d+'.png')
 

timingOfEvent()

MufiCTR.timingOfEvent	(		makeCluster = False,
			debug = False
	)

Definition at line 342 of file MufiCTR.py.

def timingOfEvent(makeCluster=False,debug=False):
   firstScifi_z = 300*u.cm
   TDC2ns = 1E9/160.316E6
   ut.bookHist(h,'evTimeDS','cor time of hits;[ns]',100,-20.,20)
   ut.bookHist(h,'evTimeScifi','cor time of hits blue DS red Scifi;[ns]',100,-20.,20)
   ut.bookCanvas(h,'tevTime','cor time of hits',1024,768,1,1)
   h['evTimeScifi'].SetLineColor(ROOT.kRed)
   h['evTimeDS'].SetLineColor(ROOT.kBlue)
   h['evTimeScifi'].SetStats(0)
   h['evTimeDS'].SetStats(0)
   if makeCluster: trackTask.scifiCluster()
   meanXY = {}
   for siCl in trackTask.clusScifi:
       detID = siCl.GetFirst()
       s = detID//1000000
       isVertical = detID%1000000//100000
       siCl.GetPosition(A,B)
       z=(A[2]+B[2])/2.
       pos = (A[1]+B[1])/2.
       L = abs(A[0]-B[0])/2.
       if isVertical:
          pos = (A[0]+B[0])/2.
          L = abs(A[1]-B[1])/2.
       corTime = geo.modules['Scifi'].GetCorrectedTime(detID, siCl.GetTime(), L) - (z-firstScifi_z)/u.speedOfLight
       h['evTimeScifi'].Fill(corTime)
       if debug: print(detID,corTime,pos)
   for aHit in eventTree.Digi_MuFilterHits:
       detID = aHit.GetDetectorID()
       if not detID//10000==3: continue
       if aHit.isVertical(): nmax = 1
       else: nmax=2
       geo.modules['MuFilter'].GetPosition(detID,A,B)
       z=(A[2]+B[2])/2.
       pos = (A[1]+B[1])/2.
       L = abs(A[0]-B[0])/2.
       if isVertical: 
          pos = (A[0]+B[0])/2.
          L = abs(A[1]-B[1])/2.
       for i in range(nmax):
            corTime = geo.modules['MuFilter'].GetCorrectedTime(detID, i, aHit.GetTime(i)*TDC2ns, L)- (z-firstScifi_z)/u.speedOfLight
            h['evTimeDS'].Fill(corTime)
            if debug: print(detID,i,corTime,pos)
   tc=h['tevTime'].cd()
   h['evTimeScifi'].Draw()
   h['evTimeDS'].Draw('same')
 
 
 

xCheck()

MufiCTR.xCheck

(

)

Definition at line 68 of file MufiCTR.py.

def xCheck():
   ut.bookCanvas(h,'xCheckH',' ',1200,600,1,1)
   ut.bookCanvas(h,'xCheckV',' ',1200,600,1,1)
   tcolors = [0,ROOT.kRed,ROOT.kRed,ROOT.kBlue,ROOT.kBlue,ROOT.kGreen,ROOT.kGreen,ROOT.kCyan,ROOT.kCyan]
   h['meanAndSig'] = {}
   for b in ['','B2noB1']:
    j=0
    try:
      tc = h['mufi-dTcorScifiDS'+b]
    except:
       continue
    h['meanAndSig'][b] = {}
    xmin,xmax = -5.,5.
    if b=='B2noB1': xmin,xmax = -20.,-8.
    for pad in tc.GetListOfPrimitives():
      j+=1
      if j==7: j==8
      tpad = tc.cd(j)
      for x in pad.GetListOfPrimitives():
         if not x.ClassName().find('TH')<0:
            if j%2==1: tcX = h['xCheckH'].cd()
            else: tcX = h['xCheckV'].cd()
            hname = x.GetName()
            h[hname] = x.Clone(hname)
            h[hname+'B1'] = h[hname].ProjectionX(hname+'B1')
            h[hname+'B1'].SetStats(0)
            h[hname+'B1'].SetLineColor(tcolors[j])
            h[hname+'B1'].SetMarkerColor(tcolors[j])
            h[hname+'B1'].SetMarkerStyle(30+j)
            rc = h[hname+'B1'].Fit('gaus','SQ','',xmin,xmax)
            fitRes = rc.Get()
            h[hname+'B1'].GetFunction('gaus').SetLineColor(tcolors[j])
            h['meanAndSig'][b][j] = [fitRes.Parameter(1),fitRes.Parameter(2)]
   plots = {'H':[0,2,4],'V':[1,3,5,6]}
   b=''
   txt = ROOT.TLatex()
   txt.SetTextFont(42)
   txt.SetTextSize(0.04)
   for c in plots:
      tcX = h['xCheck'+c].cd()
      tcX.SetLogy(1)
      j=0
      for x in plots[c]:
         if x<2: h['mufi-dTcor_3'+str(x)+'B1'].Draw()
         else: h['mufi-dTcor_3'+str(x)+'B1'].Draw('same')
         if x==6: m,s = h['meanAndSig'][b][x+2]
         else: m,s = h['meanAndSig'][b][x+1]
         h['txt'+str(x)] = txt.DrawLatexNDC(0.15,0.85-j*0.05,"DS%i%s  m=%5.2Fns  #sigma=%5.2Fns"%(x//2,c,m,s))
         h['txt'+str(x)].SetTextColor(h['mufi-dTcor_3'+str(x)+'B1'].GetLineColor())
         j+=1
      j+=1
      try:
        if c=='H':
         h['txtB2H'] = txt.DrawLatexNDC(0.15,0.85-j*0.05,"BW tracks DS0H: %5.2Fns DS1H: %5.2Fns DS2H: %5.2Fns  "%(
           h['meanAndSig']['B2noB1'][1][0],h['meanAndSig']['B2noB1'][3][0],h['meanAndSig']['B2noB1'][5][0]))
         h['txtB2H'].SetTextColor(ROOT.kMagenta)
        else:
         h['txtB2V'] = txt.DrawLatexNDC(0.15,0.85-j*0.05,"BW tracks DS0V: %5.2Fns DS1V: %5.2Fns DS2V: %5.2Fns  DS3V: %5.2Fns  "%(
           h['meanAndSig']['B2noB1'][2][0],h['meanAndSig']['B2noB1'][4][0],h['meanAndSig']['B2noB1'][6][0],h['meanAndSig']['B2noB1'][8][0]))
         h['txtB2V'].SetTextColor(ROOT.kMagenta)
      except: pass
      h['xCheck'+c].Update()
      h['xCheck'+c].Print('mufi-dTScifiDS-xCheck'+c+'.png')
 
 

Variable Documentation

alignTpar

dict MufiCTR.alignTpar = {}

Definition at line 134 of file MufiCTR.py.

alignTparB2

dict MufiCTR.alignTparB2 = {}

Definition at line 135 of file MufiCTR.py.

channelsPerPlane

dict MufiCTR.channelsPerPlane = {1:7*16,3:60*2}

Definition at line 13 of file MufiCTR.py.

default

MufiCTR.default

Definition at line 6 of file MufiCTR.py.

dest

MufiCTR.dest

Definition at line 6 of file MufiCTR.py.

f

MufiCTR.f = ROOT.TFile(options.file)

Definition at line 9 of file MufiCTR.py.

h

dict MufiCTR.h = {}

Definition at line 11 of file MufiCTR.py.

help

MufiCTR.help

Definition at line 6 of file MufiCTR.py.

options

MufiCTR.options = parser.parse_args()

Definition at line 7 of file MufiCTR.py.

par0

dict MufiCTR.par0 = {}

Definition at line 133 of file MufiCTR.py.

parser

MufiCTR.parser = ArgumentParser()

Definition at line 5 of file MufiCTR.py.

Planes

dict MufiCTR.Planes = {1:2,3:7}

Definition at line 14 of file MufiCTR.py.

runN

MufiCTR.runN = options.file[options.file.find(".root")-6:options.file.find(".root")]

Definition at line 16 of file MufiCTR.py.

S

dict MufiCTR.S = {1:[1800,800,2,1],2:[1800,1500,2,3],3:[1800,1800,2,4]}

Definition at line 12 of file MufiCTR.py.