Exemple #1
0
def create_geometry(file_name=None, pixel_width=None, pixel_height=None):
    inst_name = "REF_L"

    if pixel_width is None:
        pixel_width = PIXEL_WIDTH
    if pixel_height is None:
        pixel_height = PIXEL_HEIGHT

    if file_name is None:
        xml_outfile = inst_name + "_Definition.xml"
    else:
        xml_outfile = file_name

    det = MantidGeom(inst_name)
    det.addSnsDefaults()
    det.addComment("SOURCE AND SAMPLE POSITION")
    det.addModerator(-13.63)  #was 13.601
    det.addSamplePosition()
    det.addComment("MONITORS")
    det.addMonitors(names=["monitor1"], distance=["-0.23368"])
    det.addComment("DETECTOR")
    id_str = "detector1"
    det.addComponent("detector1", id_str)
    doc_handle = det.makeTypeElement(id_str)

    det.addPixelatedTube("tube",
                         NUM_PIXELS_PER_TUBE,
                         PIXEL_HEIGHT * NUM_PIXELS_PER_TUBE,
                         type_name="pixel")

    det.addCylinderPixel("pixel", (0.0, 0.0, 0.0), (0.0, 1.0, 0.0),
                         pixel_width / 2.0, pixel_height)

    for i in range(0, NUM_TUBES):
        det.addComponent("tube%d" % i, root=doc_handle)
        det.addDetector(str((i - NUM_TUBES / 2 + 0.5) * pixel_width), "0",
                        "+1.28", "0", "0", "0", "tube%d" % i, "tube")

    det.addComment("FIXME: Do something real here.")
    det.addDummyMonitor(0.01, 0.03)

    id_list = [0, NUM_TUBES * NUM_PIXELS_PER_TUBE - 1, None]
    det.addDetectorIds(id_str, id_list)
    det.addMonitorIds(["-1"])

    det.writeGeom(xml_outfile)
Exemple #2
0
def create_geometry(file_name=None, pixel_width=None, pixel_height=None):
    inst_name = "REF_L"
    short_name = "REF_L"

    if pixel_width is None: pixel_width = PIXEL_WIDTH
    if pixel_height is None: pixel_height = PIXEL_HEIGHT

    if file_name is None:
        xml_outfile = inst_name+"_Definition.xml"
    else:
        xml_outfile = file_name

    det = MantidGeom(inst_name)
    det.addSnsDefaults()
    det.addComment("SOURCE AND SAMPLE POSITION")
    det.addModerator(-13.601)
    det.addSamplePosition()
    det.addComment("MONITORS")
    det.addMonitors(names=["monitor1"], distance=["-0.23368"])

    id_str = "detector1"
    det.addComponent("detector1", id_str)
    doc_handle = det.makeTypeElement(id_str)

    det.addPixelatedTube("tube", NUM_PIXELS_PER_TUBE, pixel_height*NUM_PIXELS_PER_TUBE, type_name="pixel")
    
    det.addCylinderPixel("pixel", (0.0, 0.0, 0.0), (0.0, 1.0, 0.0),
                         pixel_width/2.0, pixel_height)
    
    for i in range(0, NUM_TUBES):
        det.addComponent("tube%d" % i, root=doc_handle)
        det.addDetector(str((i-NUM_TUBES/2+0.5)*pixel_width), "0", "0", "0", "0", "0", "tube%d"%i, "tube")

    det.addComment("FIXME: Do something real here.")
    det.addDummyMonitor(0.01, 0.03)

    id_list = [0, NUM_TUBES*NUM_PIXELS_PER_TUBE-1, None]
    det.addDetectorIds(id_str, id_list)
    det.addMonitorIds(["-1"])


    det.writeGeom(xml_outfile)
Exemple #3
0
def create_geometry(file_name=None,
                    tube_width=TUBE_WIDTH,
                    tube_length=TUBE_SIZE):
    """
        Create a geometry file
        @param file_name: name of the output file [optional]
        @param tube_width: width of the tubes (pixel width)
        @param tube_length: length of the tubes (number of pixels times pixel height)
    """
    # Get geometry information file
    inst_name = "EQ-SANS"
    if file_name is None:
        xml_outfile = inst_name + "_Definition.xml"
    else:
        xml_outfile = file_name

    det = MantidGeom(inst_name)
    det.addSnsDefaults()
    det.addComment("SOURCE AND SAMPLE POSITION")
    det.addModerator(-14.122)
    det.addSamplePosition()
    det.addComment("MONITORS")
    det.addMonitors(names=["monitor1"], distance=["-0.23368"])

    id_str = "detector1"
    det.addComponent(id_str, id_str)
    doc_handle = det.makeTypeElement(id_str)

    det.addCylinderPixel("pixel", (0.0, 0.0, 0.0), (0.0, 1.0, 0.0),
                         (tube_width / 2.0),
                         (tube_length / NUM_PIXELS_PER_TUBE))

    for i in range(0, NUM_BANKS / 2):
        i_low_bank = i
        i_high_bank = i + NUM_BANKS / 2
        low_bank = "bank" + str(i_low_bank + 1)
        high_bank = "bank" + str(i_high_bank + 1)

        # FRONT plane
        for j in range(NUM_TUBES_PER_BANK):
            tube_id = "bank" + str(i_low_bank + 1) + "_tube" + str(j + 1)
            x, y, z = get_position(i_low_bank, j, tube_width)
            det.addComponent(tube_id, root=doc_handle)
            det.addPixelatedTube('_' + tube_id,
                                 NUM_PIXELS_PER_TUBE,
                                 tube_length,
                                 type_name="pixel")
            det.addDetector(x, y, z, "0", "0", "0", tube_id, '_' + tube_id)

        # BACK plane
        for j in range(NUM_TUBES_PER_BANK):
            tube_id = "bank" + str(i_high_bank + 1) + "_tube" + str(j + 1)
            x, y, z = get_position(i_high_bank, j, tube_width)
            det.addComponent(tube_id, root=doc_handle)
            det.addPixelatedTube('_' + tube_id,
                                 NUM_PIXELS_PER_TUBE,
                                 tube_length,
                                 type_name="pixel")
            det.addDetector(x, y, z, "0", "0", "0", tube_id, '_' + tube_id)

    det.addComment("MONITOR SHAPE")
    det.addComment("FIXME: Do something real here.")
    det.addDummyMonitor(0.01, 0.03)

    det.addComment("DETECTOR IDs")
    id_list = [
        0, NUM_BANKS * NUM_TUBES_PER_BANK * NUM_PIXELS_PER_TUBE - 1, None
    ]
    det.addDetectorIds(id_str, id_list)

    det.addComment("MONITOR IDs")
    det.addMonitorIds(["-1"])

    det.writeGeom(xml_outfile)
Exemple #4
0
def create_geometry(file_name=None, tube_width=TUBE_WIDTH, tube_length=TUBE_SIZE):
    """
        Create a geometry file
        @param file_name: name of the output file [optional]
        @param tube_width: width of the tubes (pixel width)
        @param tube_length: length of the tubes (number of pixels times pixel height)
    """
    # Get geometry information file
    inst_name = "EQ-SANS"
    if file_name is None:
        xml_outfile = inst_name+"_Definition.xml"
    else:
        xml_outfile = file_name

    det = MantidGeom(inst_name)
    det.addSnsDefaults()
    det.addComment("SOURCE AND SAMPLE POSITION")
    det.addModerator(-14.122)
    det.addSamplePosition()
    det.addComment("MONITORS")
    det.addMonitors(names=["monitor1"], distance=["-0.23368"])

    id_str = "detector1"
    det.addComponent(id_str, id_str)
    doc_handle = det.makeTypeElement(id_str)

    det.addCylinderPixel("pixel", (0.0, 0.0, 0.0), (0.0, 1.0, 0.0),
                        (tube_width/2.0),
                        (tube_length/NUM_PIXELS_PER_TUBE))

    for i in range(0, NUM_BANKS/2):
        i_low_bank = i
        i_high_bank = i+NUM_BANKS/2
        low_bank = "bank"+str(i_low_bank+1)
        high_bank = "bank"+str(i_high_bank+1)

        # FRONT plane
        for j in range(NUM_TUBES_PER_BANK):
            tube_id = "bank"+str(i_low_bank+1)+"_tube"+str(j+1)
            x, y, z = get_position(i_low_bank, j, tube_width)
            det.addComponent(tube_id, root=doc_handle)
            det.addPixelatedTube('_'+tube_id, NUM_PIXELS_PER_TUBE, tube_length, type_name="pixel")
            det.addDetector(x, y, z, "0", "0", "0", tube_id, '_'+tube_id)

        # BACK plane
        for j in range(NUM_TUBES_PER_BANK):
            tube_id = "bank"+str(i_high_bank+1)+"_tube"+str(j+1)
            x, y, z = get_position(i_high_bank, j, tube_width)
            det.addComponent(tube_id, root=doc_handle)
            det.addPixelatedTube('_'+tube_id, NUM_PIXELS_PER_TUBE, tube_length, type_name="pixel")
            det.addDetector(x, y, z, "0", "0", "0", tube_id, '_'+tube_id)

    det.addComment("MONITOR SHAPE")
    det.addComment("FIXME: Do something real here.")
    det.addDummyMonitor(0.01, 0.03)

    det.addComment("DETECTOR IDs")
    id_list = [0, NUM_BANKS*NUM_TUBES_PER_BANK*NUM_PIXELS_PER_TUBE-1, None]
    det.addDetectorIds(id_str, id_list)

    det.addComment("MONITOR IDs")
    det.addMonitorIds(["-1"])

    det.writeGeom(xml_outfile)