def get_b_field(canvas, z = None): global bfield_hist, bfield_graph if not maus_globals.has_instance(): configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) maus_globals.birth(configuration) z_list = [] bz_list = [] for z_pos in range(-8000, 8001, 50): z_pos = float(z_pos) (bx_field, by_field, bz_field, ex_field, ey_field, ez_field) = \ field.get_field_value(0., 0., z_pos, 0.) z_list.append(z_pos) # z in mm bz_list.append(bz_field*1e3) # bz in T bfield_hist, bfield_graph = common.make_root_graph("bz vs z", z_list, "z [m]", bz_list, "B_{z} [T]") canvas.cd() bfield_hist.Draw() bfield_graph.Draw('l') if z != None: print "Plotting z", z marker_graph = ROOT.TGraph(2) marker_graph.SetPoint(0, z, -1000) marker_graph.SetPoint(1, z, +1000) marker_graph.Draw('l') root_list.append(marker_graph)
def get_b_field(canvas, z=None): global bfield_hist, bfield_graph if not maus_globals.has_instance(): configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) maus_globals.birth(configuration) z_list = [] bz_list = [] for z_pos in range(-8000, 8001, 50): z_pos = float(z_pos) (bx_field, by_field, bz_field, ex_field, ey_field, ez_field) = \ field.get_field_value(0., 0., z_pos, 0.) z_list.append(z_pos) # z in mm bz_list.append(bz_field * 1e3) # bz in T bfield_hist, bfield_graph = common.make_root_graph( "bz vs z", z_list, "z [m]", bz_list, "B_{z} [T]") canvas.cd() bfield_hist.Draw() bfield_graph.Draw('l') if z != None: print "Plotting z", z marker_graph = ROOT.TGraph(2) marker_graph.SetPoint(0, z, -1000) marker_graph.SetPoint(1, z, +1000) marker_graph.Draw('l') root_list.append(marker_graph)
def main(): """ Make a plot of z, bz """ # set up datacards print "Welcome to MAUS field map maker example" configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) # initialise field maps and geometry print "Building field maps and other initialisation (this can take a while)" maus_globals.birth(configuration) # make lists of z, bz points print "Getting field values" canvas = None graph_list = [] for r_pos, line_color in [(0., 1), (150., 2)]: z_list = [float(z_pos) for z_pos in range(12000, 18001, 50)] bz_list = [] for z_pos in z_list: (bx_field, by_field, bz_field, ex_field, ey_field, ez_field) = \ field.get_field_value(r_pos, 0., z_pos, 0.) bz_list.append(bz_field * 1e3) # bz in T print 'z:', z_pos, ' ** b:', bx_field, by_field, bz_field, \ 'e:', ex_field, ey_field, ez_field # now make a ROOT graph of bz against z print "Graphing field values" hist, graph = xboa.common.make_root_graph("x=" + str(r_pos) + " mm", z_list, "z [m]", bz_list, "B_{z} [T]") graph.SetLineColor(line_color) if canvas == None: canvas = xboa.common.make_root_canvas("bz vs z") hist.Draw() graph_list.append(graph) graph.Draw('l') canvas.Update() legend = xboa.common.make_root_legend(canvas, graph_list) print legend.GetX1NDC(), legend.GetX2NDC() legend.SetX1NDC(0.7) legend.SetX2NDC(0.9) legend.SetBorderSize(1) #legend.Draw() canvas.Print('plots/bfield_vs_z.png') # Clean up maus_globals.death()
def main(): """ Spit out material information for the example geometry. """ # set up datacards print "Welcome to MAUS geometry navigator example" configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) conf_json = json.loads(configuration) conf_json["simulation_geometry_filename"] = "GeometryNavigatorTest.dat" configuration = json.dumps(conf_json) # initialise field maps and geometry print "Building field maps and other initialisation (this can take a while)" maus_globals.birth(configuration) print "Testing Functionality" print maus_cpp.material.set_position( 0.0, 0.0, -50.0 ) print "Current Position:", maus_cpp.material.get_position() print for i in range( 100 ) : maus_cpp.material.step( 0.0, 0.0, 1.0 ) print (-50.0 + float(i)), maus_cpp.material.get_material_data() print print print "Final Position:", maus_cpp.material.get_position() print # Clean up maus_globals.death() # Finished print "Finished"
def main(): """ Make a plot of z, bz """ # set up datacards print "Welcome to MAUS field map maker example" configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) # initialise field maps and geometry print "Building field maps and other initialisation (this can take a while)" maus_globals.birth(configuration) # make lists of z, bz points print "Getting field values" z_list = [] bz_list = [] for z_pos in range(-8000, 8001, 50): z_pos = float(z_pos) (bx_field, by_field, bz_field, ex_field, ey_field, ez_field) = \ field.get_field_value(0., 0., z_pos, 0.) z_list.append(z_pos*1e-3) # z in metres bz_list.append(bz_field*1e3) # bz in T print 'z:', z_pos, ' ** b:', bx_field, by_field, bz_field, \ 'e:', ex_field, ey_field, ez_field # now make a ROOT graph of bz against z print "Graphing field values" canvas = xboa.common.make_root_canvas("bz vs z") hist, graph = xboa.common.make_root_graph("bz vs z", z_list, "z [m]", bz_list, "B_{z} [T]") hist.Draw() graph.Draw('l') canvas.Update() canvas.Print('bfield_vs_z.png') # Clean up maus_globals.death() # Finished print "Finished - press <Carriage Return> to continue"
def main(): """ Make a plot of z, bz """ # set up datacards print "Welcome to MAUS field map maker example" configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) # initialise field maps and geometry print "Building field maps and other initialisation (this can take a while)" maus_globals.birth(configuration) # make lists of z, bz points print "Getting field values" z_list = [] bz_list = [] for z_pos in range(-8000, 8001, 50): z_pos = float(z_pos) (bx_field, by_field, bz_field, ex_field, ey_field, ez_field) = \ field.get_field_value(0., 0., z_pos, 0.) z_list.append(z_pos * 1e-3) # z in metres bz_list.append(bz_field * 1e3) # bz in T print 'z:', z_pos, ' ** b:', bx_field, by_field, bz_field, \ 'e:', ex_field, ey_field, ez_field # now make a ROOT graph of bz against z print "Graphing field values" canvas = xboa.Common.make_root_canvas("bz vs z") hist, graph = xboa.Common.make_root_graph("bz vs z", z_list, "z [m]", bz_list, "B_{z} [T]") hist.Draw() graph.Draw('l') canvas.Update() canvas.Print('bfield_vs_z.png') # Clean up maus_globals.death() # Finished print "Finished"
def main(): """ Spit out material information for the example geometry. """ # set up datacards print "Welcome to MAUS geometry navigator example" configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) conf_json = json.loads(configuration) conf_json["simulation_geometry_filename"] = "GeometryNavigatorTest.dat" configuration = json.dumps(conf_json) # initialise field maps and geometry print "Building field maps and other initialisation (this can take a while)" maus_globals.birth(configuration) print "Testing Functionality" print maus_cpp.material.set_position(0.0, 0.0, -50.0) print "Current Position:", maus_cpp.material.get_position() print for i in range(100): maus_cpp.material.step(0.0, 0.0, 1.0) print(-50.0 + float(i)), maus_cpp.material.get_material_data() print print print "Final Position:", maus_cpp.material.get_position() print # Clean up maus_globals.death() # Finished print "Finished"
def initialise_maus(): configuration = Configuration.Configuration().\ getConfigJSON(command_line_args=True) maus_globals.birth(configuration)