예제 #1
0
 def NextEvent(self, i=-1):
     if i < 0: self.n += 1
     else: self.n = i
     fRun.Run(self.n, self.n + 1)  # go for first event
     # check if tracks are made from real pattern recognition
     if sTree.GetBranch("FitTracks_PR"):
         sTree.FitTracks = sTree.FitTracks_PR
     if sTree.GetBranch("fitTrack2MC_PR"):
         sTree.fitTrack2MC = sTree.fitTrack2MC_PR
     if sTree.GetBranch("Particles_PR"):
         sTree.Particles = sTree.Particles_PR
     if hasattr(self, "tracks"): self.tracks.ExecuteTask()
     if sTree.FindBranch("EcalClusters"):
         if sTree.EcalClusters.GetEntries() > 0:
             self.ecalFiller.Exec('start', sTree.EcalPointLite)
         self.calos.ExecuteTask()
     if sTree.FindBranch("Digi_SBTHits"): self.veto.ExecuteTask()
     if ROOT.gROOT.FindObject('Root Canvas EnergyLoss'):
         evd_fillEnergy.execute()
     print 'Event %i ready' % (self.n)
     # make pointsets pickable
     for x in mcHits:
         p = ROOT.gEve.GetCurrentEvent().FindChild(mcHits[x].GetName())
         if p:
             p.SetPickable(ROOT.kTRUE)
             p.SetTitle(p.__repr__())
예제 #2
0
 def NextEvent(self,i=-1):
   if i<0: self.n+=1
   else  : self.n=i
   fRun.Run(self.n,self.n+1) # go for first event
# check if tracks are made from real pattern recognition
   if sTree.GetBranch("FitTracks_PR"):    sTree.FitTracks = sTree.FitTracks_PR
   if sTree.GetBranch("fitTrack2MC_PR"):  sTree.fitTrack2MC = sTree.fitTrack2MC_PR
   if sTree.GetBranch("Particles_PR"):    sTree.Particles   = sTree.Particles_PR
   if hasattr(self,"tracks"): self.tracks.ExecuteTask()
   if sTree.FindBranch("EcalClusters"):
     if sTree.EcalClusters.GetEntries()>0:
      self.ecalFiller.Exec('start',sTree.EcalPointLite)
     self.calos.ExecuteTask()
   if sTree.FindBranch("Digi_SBTHits"): self.veto.ExecuteTask()
   if ROOT.gROOT.FindObject('Root Canvas EnergyLoss'): evd_fillEnergy.execute()
   print 'Event %i ready'%(self.n)
# make pointsets pickable
   for x in mcHits: 
     p = ROOT.gEve.GetCurrentEvent().FindChild(mcHits[x].GetName())
     if p: 
      p.SetPickable(ROOT.kTRUE)
      p.SetTitle(p.__repr__())
예제 #3
0
            ls.SetSpecularPower(1.1)
            vw.DoDraw()


#
        slot = pack1.NewSlot()
        slot.StartEmbedding()
        can = ROOT.TCanvas(
            "Root Canvas EnergyLoss")  # ROOT.gROOT.FindObject('Root Canvas')
        can.SetTitle('Energy as function of z')
        can.ToggleEditor()
        slot.StopEmbedding()
        ls = ROOT.gROOT.GetListOfGlobals()
        ls.Add(can)
        SHiPDisplay = ls.FindObject('SHiP Displayer')
        SHiPDisplay.transparentMode('on')
if __name__ == "__main__":
    execute()
    v = ROOT.gEve.GetViewers().FindChild('Bar Embedded Viewer side')
    vw = v.GetGLViewer()
    cam = vw.CurrentCamera()
    fr = vw.GetFrame()
    test = ROOT.TGLVertex3(0., 0., 0.)
    vtest = cam.ViewportToWorld(test)
    zmin = vtest.Z()
    test = ROOT.TGLVertex3(fr.GetWidth(), 0., 0.)
    vtest = cam.ViewportToWorld(test)
    zmax = vtest.Z()
    print "?", zmin, zmax
    evd_fillEnergy.execute()
예제 #4
0
   ls.SetFrontPower(0.2)
   ls.SetSidePower(0.4)
   ls.SetSpecularPower(1.1)
   vw.DoDraw()
#
  slot = pack1.NewSlot()
  slot.StartEmbedding()
  can = ROOT.TCanvas("Root Canvas EnergyLoss") # ROOT.gROOT.FindObject('Root Canvas')
  can.SetTitle('Energy as function of z')
  can.ToggleEditor()
  slot.StopEmbedding()
  ls = ROOT.gROOT.GetListOfGlobals()
  ls.Add(can)
  SHiPDisplay = ls.FindObject('SHiP Displayer')
  SHiPDisplay.transparentMode('on')
if __name__=="__main__":
  execute()
  v = ROOT.gEve.GetViewers().FindChild('Bar Embedded Viewer side')
  vw = v.GetGLViewer()
  cam = vw.CurrentCamera()
  fr = vw.GetFrame()
  test = ROOT.TGLVertex3(0.,0.,0.)
  vtest = cam.ViewportToWorld(test)
  zmin = vtest.Z()
  test = ROOT.TGLVertex3(fr.GetWidth(),0.,0.)
  vtest = cam.ViewportToWorld(test)
  zmax = vtest.Z()
  print "?",zmin,zmax
  evd_fillEnergy.execute()