110ReadTree Namespace Reference

Functions
def	printPrimaryParticle (depth, primaryParticle)
def	printPrimaryVertex (depth, primaryVertex)
def	printTrajectoryPoint (depth, trajectoryPoint)
def	printTrajectory (depth, trajectory)
def	printHitSegment (depth, hitSegment)
def	printSegmentContainer (depth, containerName, hitSegments)
def	main (argv=None)

Function Documentation

def 110ReadTree.main

(

argv = None

)

Definition at line 100 of file 110ReadTree.py.

def main(argv=None):
    if argv is None:
        argv = sys.argv

    # The input file is generated in a previous test (100TestTree.sh).
    inputFile=TFile("100TestTree.root")

    # Get the input tree out of the file.
    inputTree=inputFile.Get("EDepSimEvents")
    print("Class:",inputTree.ClassName())

    # Attach a brach to the events.
    event = TG4Event()
    inputTree.SetBranchAddress("Event",event)

    # Read all of the events.
    entries=inputTree.GetEntriesFast()
    for jentry in xrange(entries):
        nb = inputTree.GetEntry(jentry)
        if nb<=0: continue
        print("Class: ", event.ClassName())
        print("Event number:", event.EventId)
        # Dump the primary vertices
        for primaryVertex in event.Primaries:
            printPrimaryVertex("PP", primaryVertex)
        # Dump the trajectories
        for trajectory in event.Trajectories: printTrajectory("TT",trajectory)
        # Dump the segment containers
        print("Number of segment containers:", event.SegmentDetectors.size())
        for containerName, hitSegments in event.SegmentDetectors:
            printSegmentContainer("HH", containerName, hitSegments)
    

def 110ReadTree.printHitSegment	(		depth,
			hitSegment
	)

Definition at line 75 of file 110ReadTree.py.

def printHitSegment(depth, hitSegment):
    print(depth,"Class: ", hitSegment.ClassName())
    print(depth,"Primary Id:", hitSegment.GetPrimaryId());
    print(depth,"Energy Deposit:",hitSegment.GetEnergyDeposit())
    print(depth,"Secondary Deposit:", hitSegment.GetSecondaryDeposit())
    print(depth,"Track Length:",hitSegment.GetTrackLength())
    print(depth,"Start:", hitSegment.GetStart().X(),
          hitSegment.GetStart().Y(),
          hitSegment.GetStart().Z(),
          hitSegment.GetStart().T())
    print(depth,"Stop:", hitSegment.GetStop().X(),
          hitSegment.GetStop().Y(),
          hitSegment.GetStop().Z(),
          hitSegment.GetStop().T())
    print(depth,"Contributor:", [contributor for contributor in hitSegment.Contrib])

# Print the fields in a single element of the SegmentDetectors map.
# The container name is the key, and the hitSegments is the value (a
# vector of TG4HitSegment objects).

def 110ReadTree.printPrimaryParticle	(		depth,
			primaryParticle
	)

Definition at line 13 of file 110ReadTree.py.

def printPrimaryParticle(depth, primaryParticle):
    print(depth,"Class: ", primaryParticle.ClassName())
    print(depth,"Track Id:", primaryParticle.GetTrackId())
    print(depth,"Name:", primaryParticle.GetName())
    print(depth,"PDG Code:",primaryParticle.GetPDGCode())
    print(depth,"Momentum:",primaryParticle.GetMomentum().X(),
          primaryParticle.GetMomentum().Y(),
          primaryParticle.GetMomentum().Z(),
          primaryParticle.GetMomentum().E(),
          primaryParticle.GetMomentum().P(),
          primaryParticle.GetMomentum().M())

# Print the fields in an TG4PrimaryVertex object

def 110ReadTree.printPrimaryVertex	(		depth,
			primaryVertex
	)

Definition at line 26 of file 110ReadTree.py.

def printPrimaryVertex(depth, primaryVertex):
    print(depth,"Class: ", primaryVertex.ClassName())
    print(depth,"Position:", primaryVertex.GetPosition().X(),
          primaryVertex.GetPosition().Y(),
          primaryVertex.GetPosition().Z(),
          primaryVertex.GetPosition().T())
    print(depth,"Generator:",primaryVertex.GetGeneratorName())
    print(depth,"Reaction:",primaryVertex.GetReaction())
    print(depth,"Filename:",primaryVertex.GetFilename())
    print(depth,"InteractionNumber:",primaryVertex.GetInteractionNumber())
    depth = depth + ".."
    for infoVertex in primaryVertex.Informational:
        printPrimaryVertex(depth,infoVertex)
    for primaryParticle in primaryVertex.Particles:
        printPrimaryParticle(depth,primaryParticle)

# Print the fields in a TG4TrajectoryPoint object

def 110ReadTree.printSegmentContainer	(		depth,
			containerName,
			hitSegments
	)

Definition at line 94 of file 110ReadTree.py.

def printSegmentContainer(depth, containerName, hitSegments):
    print(depth,"Detector: ", containerName, hitSegments.size())
    depth = depth + ".."
    for hitSegment in hitSegments: printHitSegment(depth, hitSegment)

# Read a file and dump it.

def 110ReadTree.printTrajectory	(		depth,
			trajectory
	)

Definition at line 57 of file 110ReadTree.py.

def printTrajectory(depth, trajectory):
    print(depth,"Class: ", trajectory.ClassName())
    depth = depth + ".."
    print(depth,"Track Id/Parent Id:",
          trajectory.GetTrackId(),
          trajectory.GetParentId())
    print(depth,"Name:",trajectory.GetName())
    print(depth,"PDG Code",trajectory.GetPDGCode())
    print(depth,"Initial Momentum:",trajectory.GetInitialMomentum().X(),
          trajectory.GetInitialMomentum().Y(),
          trajectory.GetInitialMomentum().Z(),
          trajectory.GetInitialMomentum().E(),
          trajectory.GetInitialMomentum().P(),
          trajectory.GetInitialMomentum().M())
    for trajectoryPoint in trajectory.Points:
        printTrajectoryPoint(depth,trajectoryPoint)

# Print the fields in a TG4HitSegment object

def 110ReadTree.printTrajectoryPoint	(		depth,
			trajectoryPoint
	)

Definition at line 43 of file 110ReadTree.py.

def printTrajectoryPoint(depth, trajectoryPoint):
    print(depth,"Class: ", trajectoryPoint.ClassName())
    print(depth,"Position:", trajectoryPoint.GetPosition().X(),
          trajectoryPoint.GetPosition().Y(),
          trajectoryPoint.GetPosition().Z(),
          trajectoryPoint.GetPosition().T())
    print(depth,"Momentum:", trajectoryPoint.GetMomentum().X(),
          trajectoryPoint.GetMomentum().Y(),
          trajectoryPoint.GetMomentum().Z(),
          trajectoryPoint.GetMomentum().Mag())
    print(depth,"Process",trajectoryPoint.GetProcess())
    print(depth,"Subprocess",trajectoryPoint.GetSubprocess())

# Print the fields in a TG4Trajectory object