Classes
class	DrawTracks

class	DrawVetoDigi

class	EventLoop

class	IO

class	Rulers

Functions
bool	_file_accessible (str path)

None	evExit ()

None	printMCTrack (int n, MCTrack)

None	dump (float pcut=0)

None	printFittedTracks ()

None	printParticles ()

None	projection ()

None	projection_prescale ()

None	storeCameraSetting (str fname="camSetting.root")

None	readCameraSetting (str fname="camSetting.root")

None	speedUp ()

None	DisplayNuDetector ()

None	switchOff (tag)

None	switchOn (tag)

None	hidePlasticScintillator ()

None	switchOffRock ()

None	switchOffAll (exc)

None	switchOnAll (exc)

def	select (pattern)

None	search (lvdict, tag)

None	rename (str name="ship.TGeant4.root")

None	mydebug ()

None	debugStraw (n)

None	DrawCharmTracks ()

None	DrawSimpleMCTracks ()

None	positionText (r, float x, float y, z, angle, str txt, int size=200, color=ROOT.kBlue, mode=ROOT.TGLFont.kExtrude, light=ROOT.kTRUE)

None	PRVersion ()

Variables
None	fMan = None

None	fRun = None

ROOT	pdg = ROOT.TDatabasePDG.Instance()

ROOT	g = ROOT.gROOT

ROOT	gEnv = ROOT.gEnv

str	mcEngine = "TGeant4"

str	simEngine = "Pythia8"

bool	withGeo = False

bool	withMCTracks = True

dict	transparentMaterials

ArgumentParser	parser = ArgumentParser()

	dest

	help

	required

	default

	action

ArgumentParser	options = parser.parse_args()

	recoFile

	InputFile

	else :

ArgumentParser	candidate = options.InputFile.replace(".root", "_rec.root")

ROOT	inFile = ROOT.FairFileSource(options.InputFile)

	OutputFile

None	rtdb = fRun.GetRuntimeDb()

ROOT	parInput1 = ROOT.FairParRootFileIo()

int	verbose = 0

ROOT	Track = ROOT.FairMCTracksDraw("Monte-Carlo Tracks")

ROOT	GTrack = ROOT.FairGeoTracksDraw("GeoTracks")

load_from_root_file	ShipGeo = load_from_root_file(fRun.GetGeoFile(), "ShipGeo")

ROOT	_tmpFile = ROOT.TFile.Open(options.InputFile)

ROOT	_tmpTree = _tmpFile.Get("cbmsim") if _tmpFile else None

dict	mcHits = {}

dict	_candidates

ROOT	source = ROOT.VectorMCPointSource[branch](branch)

ROOT	fRman = ROOT.FairRootManager.Instance()

ROOT	sTree = fRman.GetInChain()

ROOT	lsOfGlobals = ROOT.gROOT.GetListOfGlobals()

ROOT	sGeo = ROOT.gGeoManager

ROOT	top = sGeo.GetTopVolume()

ROOT	gEve = ROOT.gEve

geomGeant4	fieldMaker = geomGeant4.addVMCFields(ShipGeo, "", True, withVirtualMC=False)

ROOT	bfield = ROOT.genfit.FairShipFields()

ROOT	bellField

ROOT	compField = ROOT.ShipCompField("compfield", bellField)

ROOT	geoMat = ROOT.genfit.TGeoMaterialInterface()

ROOT	fM = ROOT.genfit.FieldManager.getInstance()

ROOT	br = gEve.GetBrowser()

Rulers	rulers = Rulers()

EventLoop	SHiPDisplay = EventLoop()

Function Documentation

◆ _file_accessible()

bool eventDisplay._file_accessible ( str path )

protected

Check whether a file exists (local or remote via XRootD).

Definition at line 26 of file eventDisplay.py.

def _file_accessible(path: str) -> bool:
    """Check whether a file exists (local or remote via XRootD)."""
    if path.startswith("root://"):
        from urllib.parse import urlparse
 
        from XRootD import client as xrd_client
 
        parsed = urlparse(path)
        fs = xrd_client.FileSystem(f"{parsed.scheme}://{parsed.netloc}")
        status, _ = fs.stat(parsed.path)
        return status.ok
    return os.path.exists(path)
 
 

◆ debugStraw()

None eventDisplay.debugStraw ( n )

Definition at line 1152 of file eventDisplay.py.

def debugStraw(n) -> None:
    sGeo = ROOT.gGeoManager
    vols = sGeo.GetListOfVolumes()
    sTree = g.FindObjectAny("cbmsim")
    sTree.GetEntry(n)
    for s in sTree.strawtubesPoint:
        print(vols[s.GetDetectorID() - 1].GetName())
 
 
# ----Load the default libraries------

◆ DisplayNuDetector()

None eventDisplay.DisplayNuDetector ( )

Definition at line 876 of file eventDisplay.py.

def DisplayNuDetector() -> None:
    for x in ["Wall"]:
        xvol = sGeo.GetVolume(x)
        if not xvol:
            continue
        xvol.SetVisDaughters(0)
        xvol.SetVisibility(1)
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    gEve.ElementChanged(geoscene, True, True)
 
 

◆ DrawCharmTracks()

None eventDisplay.DrawCharmTracks ( )

Definition at line 1299 of file eventDisplay.py.

def DrawCharmTracks() -> None:
    i = -1
    for aTrack in sTree.MCTrack:
        i += 1
        if i < 2:
            continue
        if aTrack.GetMotherId() == 1:
            pa = pdg.GetParticle(sTree.MCTrack[i].GetPdgCode())
            if pa.Lifetime() > 1.0e-12:
                print(sTree.MCTrack[i])
                SHiPDisplay.tracks.DrawMCTrack(i)
 
 

◆ DrawSimpleMCTracks()

None eventDisplay.DrawSimpleMCTracks ( )

Definition at line 1312 of file eventDisplay.py.

def DrawSimpleMCTracks() -> None:
    comp = SHiPDisplay.tracks.comp
    comp.OpenCompound()
    n = -1
    ntot = 0
    fPos = ROOT.TVector3()
    fMom = ROOT.TVector3()
    delZ = 10 * u.m
    for fT in sTree.MCTrack:
        n += 1
        DTrack = ROOT.TEveLine()
        DTrack.SetPickable(ROOT.kTRUE)
        DTrack.SetTitle(fT.__repr__())
        fT.GetStartVertex(fPos)
        fT.GetMomentum(fMom)
        hitlist = {}
        hitlist[fPos.Z()] = [fPos.X(), fPos.Y()]
        z = fPos.Z() + delZ
        slx, sly = fMom.X() / fMom.Z(), fMom.Y() / fMom.Z()
        hitlist[z] = [fPos.X() + slx * delZ, fPos.Y() + sly * delZ]
        for z in hitlist:
            DTrack.SetNextPoint(hitlist[z][0], hitlist[z][1], z)
        p = pdg.GetParticle(fT.GetPdgCode())
        if p:
            pName = p.GetName()
        else:
            pName = str(fT.GetPdgCode())
        DTrack.SetName("MCTrack_" + str(n) + "_" + pName)
        c = ROOT.kYellow
        DTrack.SetMainColor(c)
        DTrack.SetLineWidth(3)
        comp.AddElement(DTrack)
        ntot += 1
    comp.CloseCompound()
    gEve.ElementChanged(SHiPDisplay.tracks.evscene, True, True)
 
 

◆ dump()

None eventDisplay.dump ( float pcut = 0 )

Definition at line 153 of file eventDisplay.py.

def dump(pcut: float = 0) -> None:
    print("   #         pid   px    py      pz     vx      vy       vz      mid")
    n = -1
    for mcp in sTree.MCTrack:
        n += 1
        if mcp.GetP() / u.GeV < pcut:
            continue
        printMCTrack(n, sTree.MCTrack)
 
 

◆ evExit()

None eventDisplay.evExit ( )

Prevent double delete due to a FairRoot bug.

Definition at line 40 of file eventDisplay.py.

def evExit() -> None:
    """Prevent double delete due to a FairRoot bug."""
    # Check whether the Eve window was closed/destructed
    if ROOT.addressof(ROOT.gEve) == 0:
        # Prevent the FairEventManager destructor from being called
        ROOT.SetOwnership(fMan, False)
 
 

◆ hidePlasticScintillator()

None eventDisplay.hidePlasticScintillator ( )

Definition at line 911 of file eventDisplay.py.

def hidePlasticScintillator() -> None:
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    v = sGeo.FindVolumeFast("vleft")
    v.SetVisibility(0)
    v.SetVisDaughters(0)
    for v in sGeo.GetListOfVolumes():
        if v.GetName().find("wallVeto") > 0:
            v.SetVisibility(0)
            v.SetVisDaughters(0)
    gEve.ElementChanged(geoscene, True, True)
 
 
# switch off drawing of rock

◆ mydebug()

None eventDisplay.mydebug ( )

Definition at line 1103 of file eventDisplay.py.

def mydebug() -> None:
    t = g.FindObjectAny("cbmsim")
    nev = t.GetEntriesFast()
    t.GetEntry(0)
    # Loop over Geo tracks
    for i in range(min(5, nev)):
        t.GetEntry(i)
        for gTr in t.GeoTracks:
            gTr.Print()
            gTr.GetParticle()
            lorv = ROOT.TLorentzVector()
            print(
                "xyz E pxpypz",
                gTr.GetPoint(0)[0],
                gTr.GetPoint(0)[1],
                gTr.GetPoint(0)[2],
                lorv.E(),
                lorv.Px(),
                lorv.Py(),
                lorv.Pz(),
            )
    # Loop over MC tracks
    for i in range(min(5, nev)):
        t.GetEntry(i)
        for gMCTr in t.MCTrack:
            gMCTr.Print()
            print(gMCTr.GetPdgCode(), gMCTr.GetMass(), gMCTr.GetP())
    # MC event header
    for i in range(nev):
        t.GetEntry(i)
        print(
            t.MCEventHeader.GetEventID(),
            t.MCEventHeader.GetRunID(),
            t.MCEventHeader.GetZ(),
        )
    # Load geometry
    sGeo = ROOT.gGeoManager
    cave = sGeo.GetTopVolume()
    cave.Draw("ogl")
    # eve
    gEve = ROOT.gEve
    #
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    topnode = geoscene.FindChild("cave_1")
    topnode.SetVisLevel(4)
    gEve.ElementChanged(geoscene, True, True)
 
 

◆ positionText()

None eventDisplay.positionText	(		r,
		float	x,
		float	y,
			z,
			angle,
		str	txt,
		int	size = `200`,
			color = `ROOT.kBlue`,
			mode = `ROOT.TGLFont.kExtrude`,
			light = `ROOT.kTRUE`
	)

Definition at line 1349 of file eventDisplay.py.

) -> None:
    tt = ROOT.TEveText(txt)
    tt.SetFontSize(size)
    tt.RefMainTrans().SetPos(x, y, z)
    tt.PtrMainTrans().RotateLF(1, 3, angle)
    tt.SetMainColor(color)
    tt.SetFontMode(mode)
    tt.SetLighting(light)
    r.AddElement(tt)
 
 

◆ printFittedTracks()

None eventDisplay.printFittedTracks ( )

Definition at line 163 of file eventDisplay.py.

def printFittedTracks() -> None:
    print("  # converged Ndf chi2/Ndf    P      Pt      MCid")
    n = -1
    for ft in sTree.FitTracks:
        n += 1
        fitStatus = ft.getFitStatus()
        fitState = ft.getFittedState()
        mom = fitState.getMom()
        print(
            "%3i %6i   %4i %6.3F   %6.3F %6.3F %6i "
            % (
                n,
                fitStatus.isFitConverged(),
                fitStatus.getNdf(),
                fitStatus.getChi2() / fitStatus.getNdf(),
                mom.Mag() / u.GeV,
                mom.Pt() / u.GeV,
                sTree.fitTrack2MC[n],
            )
        )
 
 

◆ printMCTrack()

None eventDisplay.printMCTrack	(	int	n,
			MCTrack
	)

Definition at line 135 of file eventDisplay.py.

def printMCTrack(n: int, MCTrack) -> None:
    mcp = MCTrack[n]
    print(
        " %6i %7i %6.3F %6.3F %7.3F %7.3F %7.3F %7.3F %6i "
        % (
            n,
            mcp.GetPdgCode(),
            mcp.GetPx() / u.GeV,
            mcp.GetPy() / u.GeV,
            mcp.GetPz() / u.GeV,
            mcp.GetStartX() / u.m,
            mcp.GetStartY() / u.m,
            mcp.GetStartZ() / u.m,
            mcp.GetMotherId(),
        )
    )
 
 

◆ printParticles()

None eventDisplay.printParticles ( )

Definition at line 185 of file eventDisplay.py.

def printParticles() -> None:
    print("  #    P    Pt[GeV/c]   DOCA[mm]    Rsq    Vz[m]    d1    d2")
    n = -1
    for aP in sTree.Particles:
        n += 1
        doca = -1.0
        if aP.GetMother(1) == 99:  # DOCA is set
            doca = aP.T()
        Rsq = (aP.Vx() / (2.45 * u.m)) ** 2 + (aP.Vy() / ((10.0 / 2.0 - 0.05) * u.m)) ** 2
        print(
            "%3i %6.3F  %6.3F  %9.3F    %6.3F   %6.3F %4i  %4i "
            % (
                n,
                aP.P() / u.GeV,
                aP.Pt() / u.GeV,
                doca / u.mm,
                Rsq,
                aP.Vz() / u.m,
                aP.GetDaughter(0),
                aP.GetDaughter(1),
            )
        )
 
 

◆ projection()

None eventDisplay.projection ( )

Definition at line 795 of file eventDisplay.py.

def projection() -> None:
    # if 1>0:
    # camera
    s = ROOT.gEve.SpawnNewScene("Projected Event")
    ROOT.gEve.GetDefaultViewer().AddScene(s)
    v = ROOT.gEve.GetDefaultGLViewer()
    v.SetCurrentCamera(ROOT.TGLViewer.kCameraOrthoXOY)
    cam = v.CurrentCamera()
    cam.SetZoomMinMax(0.2, 20)
    # projections
    mng = ROOT.TEveProjectionManager(ROOT.TEveProjection.kPT_RPhi)
    s.AddElement(mng)
    axes = ROOT.TEveProjectionAxes(mng)
    axes.SetTitle("TEveProjections demo")
    s.AddElement(axes)
    ROOT.gEve.AddToListTree(axes, ROOT.kTRUE)
    ROOT.gEve.AddToListTree(mng, ROOT.kTRUE)
 
 

◆ projection_prescale()

None eventDisplay.projection_prescale ( )

Definition at line 814 of file eventDisplay.py.

def projection_prescale() -> None:
    # if 1>0:
    v = gEve.GetViewers()
    vw = v.FindChild("Viewer 1")
    if vw:
        vw.SetName("3d")
    sev = ROOT.gEve.SpawnNewViewer("Scaled 2D")
    smng = ROOT.TEveProjectionManager(ROOT.TEveProjection.kPP_Plane)
    sp = smng.GetProjection()
    sp.SetUsePreScale(ROOT.kTRUE)
    sp.AddPreScaleEntry(2, 100000000.0, 0.1)
    ss = ROOT.gEve.SpawnNewScene("Scaled Geom")
    sev.AddScene(ss)
    ss.AddElement(smng)
    N = sGeo.GetTopNode()
    TNod = ROOT.TEveGeoTopNode(sGeo, N, 1, 3, 10)
    ss.AddElement(TNod)
    eventscene = ROOT.gEve.SpawnNewScene("Scaled event")
    eventscene.AddElement(ROOT.FairEventManager.Instance())
    sev.AddScene(eventscene)
    eventscene.AddElement(smng)
    ROOT.gEve.GetBrowser().GetTabRight().SetTab(1)
    ROOT.gEve.FullRedraw3D(ROOT.kTRUE)
 
 

◆ PRVersion()

None eventDisplay.PRVersion ( )

Definition at line 1371 of file eventDisplay.py.

def PRVersion() -> None:
    readCameraSetting()
    for x in [
        "moreShieldingSide",
        "moreShieldingTopBot",
        "CoatWall",
        "CoatVol",
        "AbsorberVol",
    ]:
        vol = ROOT.gGeoManager.FindVolumeFast(x)
        vol.SetVisibility(0)
    ROOT.gGeoManager.GetMaterial("Concrete").SetTransparency(0)
    r = rulers.ruler
    ticks = 5
    r.DestroyElements()
    r.OpenCompound()
    xpos, ypos = -500.0, -1500.0
    zstart = ShipGeo.target.z0
    zlength = ShipGeo.MuonStation3.z - zstart + 10 * u.m
    z = zstart
    for i in range(int(zlength / 100 / ticks)):
        m = ROOT.TEveLine()
        m.SetNextPoint(xpos, ypos, z)
        m.SetNextPoint(xpos - 1 * u.m, ypos, z)
        m.SetMainColor(ROOT.kRed)
        m.SetLineWidth(5)
        r.AddElement(m)
        t1 = ROOT.TEveText(str(i * ticks) + "m")
        t1.SetMainColor(ROOT.kGray + 3)
        t1.SetFontSize(5)
        t1.RefMainTrans().SetPos(xpos - 0.1 * u.m, ypos + 0.2 * u.m, z)
        r.AddElement(t1)
        z += ticks * u.m
    xpos, ypos = 0.0, 0.0
    z = ShipGeo.MuonStation3.z + 6 * u.m
    ylength = 7 * u.m
    ypos = -ylength
    for i in range(-int(ylength / 100), int(ylength / 100), 1):
        m = ROOT.TEveLine()
        m.SetNextPoint(xpos, ypos, z)
        m.SetNextPoint(xpos + 0.05 * u.m, ypos, z)
        m.SetMainColor(ROOT.kRed)
        m.SetLineWidth(3)
        r.AddElement(m)
        t1 = ROOT.TEveText(str(i) + "m")
        t1.SetMainColor(ROOT.kGray + 3)
        t1.SetFontSize(5)
        t1.RefMainTrans().SetPos(xpos - 0.5 * u.m, ypos, z)
        r.AddElement(t1)
        ypos += 1 * u.m
    ty = ROOT.TEveText("y-axis")
    ty.SetFontSize(10)
    ty.RefMainTrans().SetPos(0.0, ypos + 1 * u.m, z)
    ty.SetMainColor(ROOT.kRed - 2)
    r.AddElement(ty)
    xpos, ypos = 0.0, 0.0
    z = ShipGeo.MuonStation3.z + 10 * u.m
    xlength = 3 * u.m
    xpos = -xlength
    for i in range(-int(xlength / 100), int(xlength / 100), 1):
        m = ROOT.TEveLine()
        m.SetNextPoint(xpos, ypos, z)
        m.SetNextPoint(xpos, ypos - 0.05 * u.m, z)
        m.SetMainColor(ROOT.kRed)
        m.SetLineWidth(3)
        r.AddElement(m)
        t1 = ROOT.TEveText(str(i) + "m")
        t1.SetMainColor(ROOT.kGray + 3)
        t1.SetFontSize(5)
        t1.RefMainTrans().SetPos(xpos, ypos - 0.1 * u.m, z)
        r.AddElement(t1)
        xpos += 1 * u.m
    tx = ROOT.TEveText("x-axis")
    tx.SetFontSize(10)
    tx.RefMainTrans().SetPos(xpos + 1 * u.m, 0.0, z)
    tx.SetMainColor(ROOT.kRed - 2)
    r.AddElement(tx)
    rotAngle = ROOT.TMath.Pi() + ROOT.TMath.PiOver2() * 5.0 / 2.0
    positionText(
        r,
        0.0,
        900.0,
        ShipGeo.TrackStation1.z - 20 * u.m,
        rotAngle,
        "SHiP",
        200,
        ROOT.kOrange - 2,
    )
    positionText(
        r,
        0.0,
        750.0,
        ShipGeo.TrackStation1.z - 40 * u.m,
        rotAngle,
        "Vacuum decay vessel",
        200,
        ROOT.kGray + 1,
    )
    positionText(r, 0.0, 100.0, ShipGeo.target.z - 6 * u.m, rotAngle, "Target", 200, ROOT.kBlue)
    positionText(
        r,
        0.0,
        600.0,
        ShipGeo.muShield.z - 10 * u.m,
        rotAngle,
        "Active muon shield",
        200,
        ROOT.kGreen - 2,
    )
    positionText(
        r,
        0.0,
        600.0,
        ShipGeo.tauMudet.zMudetC - 10 * u.m,
        rotAngle,
        "Tau neutrino detector",
        200,
        ROOT.kRed - 2,
    )
    positionText(
        r,
        0.0,
        900.0,
        ShipGeo.Bfield.z - 5 * u.m,
        rotAngle,
        "Dipole Magnet",
        200,
        ROOT.kBlue + 2,
    )
    positionText(
        r,
        -1500.0,
        -800.0,
        ShipGeo.TrackStation3.z - 2 * u.m,
        rotAngle,
        "Strawtracker",
        200,
        ROOT.kRed + 2,
    )
    positionText(r, 0.0, 700.0, ShipGeo.MuonFilter2.z, rotAngle, "Muon", 200, ROOT.kGreen + 2)
    r.CloseCompound()
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    gEve.ElementChanged(geoscene, True, True)

◆ readCameraSetting()

None eventDisplay.readCameraSetting ( str fname = "camSetting.root" )

Definition at line 846 of file eventDisplay.py.

def readCameraSetting(fname: str = "camSetting.root") -> None:
    f = ROOT.TFile.Open(fname)
    cam = ROOT.gEve.GetDefaultGLViewer().CurrentCamera()
    f.GetKey(cam.ClassName()).Read(cam)
    cam.IncTimeStamp()
    gEve.GetDefaultGLViewer().RequestDraw()
    f.Close()
 
 

◆ rename()

None eventDisplay.rename ( str name = "ship.TGeant4.root" )

Definition at line 983 of file eventDisplay.py.

def rename(name: str = "ship.TGeant4.root") -> None:
    f = ROOT.TFile(name, "UPDATE")
    t = f.Get("cbmsim")
    for x in t.GetListOfBranches():
        nm = x.GetName().replace("_1", "")
        x.SetName(nm)
    t.Write()
    f.Close()
 
 

◆ search()

None eventDisplay.search	(	lvdict,
		tag
	)

Definition at line 977 of file eventDisplay.py.

def search(lvdict, tag) -> None:
    for x in lvdict:
        if not x.find(tag) < 0:
            print(x)
 
 

◆ select()

def eventDisplay.select ( pattern )

Definition at line 968 of file eventDisplay.py.

def select(pattern):
    exc = []
    for v in sGeo.GetListOfVolumes():
        vname = v.GetName()
        if not vname.find(pattern) < 0:
            exc.append(vname)
    return exc
 
 

◆ speedUp()

None eventDisplay.speedUp ( )

Definition at line 855 of file eventDisplay.py.

def speedUp() -> None:
    for x in ["wire", "gas", "rockD", "rockS", "rockSFe"]:
        xvol = sGeo.GetVolume(x)
        if xvol:
            xvol.SetVisibility(0)
    for k in range(1, 7):
        va = sGeo.GetVolume("T" + str(k))
        if not va:
            continue
        for x in va.GetNodes():
            nm = x.GetName()
            if not nm.find("Inner") < 0 and k < 3:
                x.SetVisDaughters(False)
                x.SetVisibility(False)
            if not nm.find("LiSc") < 0:
                x.SetVisDaughters(False)
            if not nm.find("RibPhi") < 0:
                x.SetVisDaughters(False)
 
 
# set display properties for tau nu target

◆ storeCameraSetting()

None eventDisplay.storeCameraSetting ( str fname = "camSetting.root" )

Definition at line 839 of file eventDisplay.py.

def storeCameraSetting(fname: str = "camSetting.root") -> None:
    f = ROOT.TFile.Open(fname, "RECREATE")
    cam = ROOT.gEve.GetDefaultGLViewer().CurrentCamera()
    cam.Write()
    f.Close()
 
 

◆ switchOff()

None eventDisplay.switchOff ( tag )

Definition at line 888 of file eventDisplay.py.

def switchOff(tag) -> None:
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    for v in top.GetNodes():
        vname = v.GetName()
        if not vname.find(tag) < 0:
            v.SetVisibility(0)
            v.SetVisDaughters(0)
    gEve.ElementChanged(geoscene, True, True)
 
 

◆ switchOffAll()

None eventDisplay.switchOffAll ( exc )

Definition at line 934 of file eventDisplay.py.

def switchOffAll(exc) -> None:
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    for v in top.GetNodes():
        vname = v.GetName()
        if not vname.find("cave") < 0:
            continue
        todo = True
        for tag in exc:
            if not tag.find(vname) < 0:
                todo = False
        if todo:
            v.SetVisibility(0)
            v.SetVisDaughters(0)
    gEve.ElementChanged(geoscene, True, True)
 
 

◆ switchOffRock()

None eventDisplay.switchOffRock ( )

Definition at line 925 of file eventDisplay.py.

def switchOffRock() -> None:
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    for x in ["rockD", "rockS"]:
        v = sGeo.FindVolumeFast(x)
        v.SetVisibility(0)
    gEve.ElementChanged(geoscene, True, True)
 
 

◆ switchOn()

None eventDisplay.switchOn ( tag )

Definition at line 899 of file eventDisplay.py.

def switchOn(tag) -> None:
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    for v in top.GetNodes():
        vname = v.GetName()
        if not vname.find(tag) < 0:
            print("switch on ", vname)
            v.SetVisibility(1)
            v.SetVisDaughters(1)
    gEve.ElementChanged(geoscene, True, True)
 
 

◆ switchOnAll()

None eventDisplay.switchOnAll ( exc )

Definition at line 951 of file eventDisplay.py.

def switchOnAll(exc) -> None:
    sc = gEve.GetScenes()
    geoscene = sc.FindChild("Geometry scene")
    for v in top.GetNodes():
        vname = v.GetName()
        if not vname.find("cave") < 0:
            continue
        todo = True
        for tag in exc:
            if not tag.find(vname) < 0:
                todo = False
        if todo:
            v.SetVisibility(1)
            v.SetVisDaughters(1)
    gEve.ElementChanged(geoscene, True, True)
 
 

Variable Documentation

◆ _candidates

dict eventDisplay._candidates

protected

Initial value:

=  {
    "VetoPoints": ("vetoPoint", ROOT.kBlue, ROOT.kFullDiamond),
    "TimeDetPoints": ("TimeDetPoint", ROOT.kBlue, ROOT.kFullDiamond),
    "StrawPoints": ("strawtubesPoint", ROOT.kGreen, ROOT.kFullCircle),
    "TargetPoints": ("TargetPoint", ROOT.kRed, ROOT.kFullSquare),
    "MTCDetPoint": ("MTCDetPoint", ROOT.kGreen, ROOT.kFullSquare),
    "SiliconTargetPoint": ("SiliconTargetPoint", ROOT.kCyan, ROOT.kFullSquare),
}

Definition at line 1207 of file eventDisplay.py.

◆ _tmpFile

ROOT eventDisplay._tmpFile = ROOT.TFile.Open(options.InputFile)

protected

Definition at line 1203 of file eventDisplay.py.

◆ _tmpTree

ROOT eventDisplay._tmpTree = _tmpFile.Get("cbmsim") if _tmpFile else None

protected

Definition at line 1204 of file eventDisplay.py.

◆ action

eventDisplay.action

Definition at line 101 of file eventDisplay.py.

◆ bellField

ROOT eventDisplay.bellField

Initial value:

=  ROOT.ShipBellField(
            "bellfield", ShipGeo.Bfield.max, ShipGeo.Bfield.z, 2, ShipGeo.Bfield.y / 2.0 * u.m
        )

Definition at line 1255 of file eventDisplay.py.

◆ bfield

ROOT eventDisplay.bfield = ROOT.genfit.FairShipFields()

Definition at line 1252 of file eventDisplay.py.

◆ br

ROOT eventDisplay.br = gEve.GetBrowser()

Definition at line 1268 of file eventDisplay.py.

◆ candidate

ArgumentParser eventDisplay.candidate = options.InputFile.replace(".root", "_rec.root")

Definition at line 130 of file eventDisplay.py.

◆ compField

ROOT eventDisplay.compField = ROOT.ShipCompField("compfield", bellField)

Definition at line 1258 of file eventDisplay.py.

◆ default

eventDisplay.default

Definition at line 94 of file eventDisplay.py.

◆ dest

eventDisplay.dest

Definition at line 83 of file eventDisplay.py.

◆ else

eventDisplay.else :

Definition at line 129 of file eventDisplay.py.

◆ fieldMaker

geomGeant4 eventDisplay.fieldMaker = geomGeant4.addVMCFields(ShipGeo, "", True, withVirtualMC=False)

Definition at line 1251 of file eventDisplay.py.

◆ fM

ROOT eventDisplay.fM = ROOT.genfit.FieldManager.getInstance()

Definition at line 1263 of file eventDisplay.py.

◆ fMan

ROOT eventDisplay.fMan = None

Definition at line 50 of file eventDisplay.py.

◆ fRman

ROOT eventDisplay.fRman = ROOT.FairRootManager.Instance()

Definition at line 1233 of file eventDisplay.py.

◆ fRun

ROOT eventDisplay.fRun = None

Definition at line 51 of file eventDisplay.py.

◆ g

ROOT eventDisplay.g = ROOT.gROOT

Definition at line 53 of file eventDisplay.py.

◆ gEnv

ROOT eventDisplay.gEnv = ROOT.gEnv

Definition at line 54 of file eventDisplay.py.

◆ geoMat

ROOT eventDisplay.geoMat = ROOT.genfit.TGeoMaterialInterface()

Definition at line 1261 of file eventDisplay.py.

◆ gEve

ROOT eventDisplay.gEve = ROOT.gEve

Definition at line 1243 of file eventDisplay.py.

◆ GTrack

ROOT eventDisplay.GTrack = ROOT.FairGeoTracksDraw("GeoTracks")

Definition at line 1195 of file eventDisplay.py.

◆ help

eventDisplay.help

Definition at line 83 of file eventDisplay.py.

◆ inFile

ROOT eventDisplay.inFile = ROOT.FairFileSource(options.InputFile)

Definition at line 1173 of file eventDisplay.py.

◆ InputFile

eventDisplay.InputFile

Definition at line 128 of file eventDisplay.py.

◆ lsOfGlobals

ROOT eventDisplay.lsOfGlobals = ROOT.gROOT.GetListOfGlobals()

Definition at line 1237 of file eventDisplay.py.

◆ mcEngine

str eventDisplay.mcEngine = "TGeant4"

Definition at line 57 of file eventDisplay.py.

◆ mcHits

dict eventDisplay.mcHits = {}

Definition at line 1206 of file eventDisplay.py.

◆ options

ArgumentParser eventDisplay.options = parser.parse_args()

Definition at line 119 of file eventDisplay.py.

◆ OutputFile

eventDisplay.OutputFile

Definition at line 1176 of file eventDisplay.py.

◆ parInput1

ROOT eventDisplay.parInput1 = ROOT.FairParRootFileIo()

Definition at line 1181 of file eventDisplay.py.

◆ parser

ArgumentParser eventDisplay.parser = ArgumentParser()

Definition at line 81 of file eventDisplay.py.

◆ pdg

ROOT eventDisplay.pdg = ROOT.TDatabasePDG.Instance()

Definition at line 52 of file eventDisplay.py.

◆ recoFile

eventDisplay.recoFile

Definition at line 127 of file eventDisplay.py.

◆ required

eventDisplay.required

Definition at line 83 of file eventDisplay.py.

◆ rtdb

None eventDisplay.rtdb = fRun.GetRuntimeDb()

Definition at line 1180 of file eventDisplay.py.

◆ rulers

Rulers eventDisplay.rulers = Rulers()

Definition at line 1275 of file eventDisplay.py.

◆ sGeo

ROOT eventDisplay.sGeo = ROOT.gGeoManager

Definition at line 1239 of file eventDisplay.py.

◆ SHiPDisplay

eventDisplay.SHiPDisplay = EventLoop()

Definition at line 1276 of file eventDisplay.py.

◆ ShipGeo

load_from_root_file eventDisplay.ShipGeo = load_from_root_file(fRun.GetGeoFile(), "ShipGeo")

Definition at line 1200 of file eventDisplay.py.

◆ simEngine

str eventDisplay.simEngine = "Pythia8"

Definition at line 58 of file eventDisplay.py.

◆ source

ROOT eventDisplay.source = ROOT.VectorMCPointSource[branch](branch)

Definition at line 1217 of file eventDisplay.py.

◆ sTree

ROOT eventDisplay.sTree = fRman.GetInChain()

Definition at line 1234 of file eventDisplay.py.

◆ top

ROOT eventDisplay.top = sGeo.GetTopVolume()

Definition at line 1240 of file eventDisplay.py.

◆ Track

ROOT eventDisplay.Track = ROOT.FairMCTracksDraw("Monte-Carlo Tracks")

Definition at line 1194 of file eventDisplay.py.

◆ transparentMaterials

dict eventDisplay.transparentMaterials

Initial value:

=  {
    "iron": 80,
    "aluminium": 80,
    "mylar": 60,
    "STTmix9010_2bar": 95,
    "steel": 80,
    "Aluminum": 80,
    "Scintillator": 80,
    #                        tau nu detector
    "CoilCopper": 70,
    "copper": 90,
    "HPTgas": 70,
    "Bakelite": 70,
    "RPCgas": 70,
    "TTmedium": 70,
}

Definition at line 63 of file eventDisplay.py.

◆ verbose

int eventDisplay.verbose = 0

Definition at line 1192 of file eventDisplay.py.

◆ withGeo

bool eventDisplay.withGeo = False

Definition at line 59 of file eventDisplay.py.

◆ withMCTracks

bool eventDisplay.withMCTracks = True

Definition at line 60 of file eventDisplay.py.

Classes

Functions

Variables

Function Documentation

◆ _file_accessible()

◆ debugStraw()

◆ DisplayNuDetector()

◆ DrawCharmTracks()

◆ DrawSimpleMCTracks()

◆ dump()

◆ evExit()

◆ hidePlasticScintillator()

◆ mydebug()

◆ positionText()

◆ printFittedTracks()

◆ printMCTrack()

◆ printParticles()

◆ projection()

◆ projection_prescale()

◆ PRVersion()

◆ readCameraSetting()

◆ rename()

◆ search()

◆ select()

◆ speedUp()

◆ storeCameraSetting()

◆ switchOff()

◆ switchOffAll()

◆ switchOffRock()

◆ switchOn()

◆ switchOnAll()

Variable Documentation

◆ _candidates

◆ _tmpFile

◆ _tmpTree

◆ action

◆ bellField

◆ bfield

◆ br

◆ candidate

◆ compField

◆ default

◆ dest

◆ else

◆ fieldMaker

◆ fM

◆ fMan

◆ fRman

◆ fRun

◆ g

◆ gEnv

◆ geoMat

◆ gEve

◆ GTrack

◆ help

◆ inFile

◆ InputFile

◆ lsOfGlobals

◆ mcEngine

◆ mcHits

◆ options

◆ OutputFile

◆ parInput1

◆ parser

◆ pdg

◆ recoFile

◆ required

◆ rtdb

◆ rulers

◆ sGeo

◆ SHiPDisplay

◆ ShipGeo

◆ simEngine

◆ source

◆ sTree

◆ top

◆ Track

◆ transparentMaterials

◆ verbose

◆ withGeo