def test_scatterercopy(self): '''scatterercopy''' # create a weird shape from mccomposite.geometry import primitives shape = primitives.block((1, 1, 1)) #create pure python representation of scatterer composite composite1 = mccomposite.composite(shape) nprinter = NeutronPrinter(shape) composite1.addElement(nprinter) #create a copy copy = mccomposite.scatterercopy(composite1) #create a larget composite shape = primitives.block((1, 1, 2)) composite = mccomposite.composite(shape) composite.addElement(composite1, (0, 0, -0.5)) composite.addElement(copy, (0, 0, +0.5)) #render the c++ representation ccomposite = mccomposite.scattererEngine(composite) ev = mcni.neutron(r=(0, 0, -5), v=(0, 0, 1)) ccomposite.scatter(ev) return
def test_copy(self): '''copy''' print "This test creates two identical blocks, each of which "\ "does not interact with neutrons. They print the info "\ "about the neutrons passing thru them, however. "\ "This test then send one neutron through these two "\ "blocks, so we should see two printings of neutron info, "\ "differing only on time-of-flight." # create a shape from mccomposite.geometry import primitives smallblock = primitives.block((1, 1, 1)) #create pure python representation of scatterer composite composite1 = mccomposite.composite(smallblock) import UseNeutronPrinter2 nprinter = UseNeutronPrinter2.NeutronPrinter(smallblock) composite1.addElement(nprinter) #create a copy copy = Copy.Copy(composite1) #create a larget composite largeblock = primitives.block((1, 1, 2)) composite = mccomposite.composite(largeblock) composite.addElement(composite1, (0, 0, -0.5)) #composite.addElement( nprinter, (0,0,-0.5) ) composite.addElement(copy, (0, 0, +0.5)) #render the c++ representation ccomposite = mccomposite.scattererEngine(composite) ev = mcni.neutron(r=(0, 0, -5), v=(0, 0, 1)) ccomposite.scatter(ev) return
def test_scatterercopy(self): '''scatterercopy''' # create a weird shape from mccomposite.geometry import primitives shape = primitives.block( (1,1,1) ) #create pure python representation of scatterer composite composite1 = mccomposite.composite( shape ) nprinter = NeutronPrinter( shape ) composite1.addElement( nprinter ) #create a copy copy = mccomposite.scatterercopy( composite1 ) #create a larget composite shape = primitives.block( (1,1,2) ) composite = mccomposite.composite( shape ) composite.addElement( composite1, (0,0,-0.5) ) composite.addElement( copy, (0,0,+0.5) ) #render the c++ representation ccomposite = mccomposite.scattererEngine( composite ) ev = mcni.neutron( r = (0,0,-5), v = (0,0,1) ) ccomposite.scatter(ev) return
def test_copy(self): '''copy''' print "This test creates two identical blocks, each of which "\ "does not interact with neutrons. They print the info "\ "about the neutrons passing thru them, however. "\ "This test then send one neutron through these two "\ "blocks, so we should see two printings of neutron info, "\ "differing only on time-of-flight." # create a shape from mccomposite.geometry import primitives smallblock = primitives.block( (1,1,1) ) #create pure python representation of scatterer composite composite1 = mccomposite.composite( smallblock ) import UseNeutronPrinter2 nprinter = UseNeutronPrinter2.NeutronPrinter( smallblock ) composite1.addElement( nprinter ) #create a copy copy = Copy.Copy( composite1 ) #create a larget composite largeblock = primitives.block( (1,1,2) ) composite = mccomposite.composite( largeblock ) composite.addElement( composite1, (0,0,-0.5) ) #composite.addElement( nprinter, (0,0,-0.5) ) composite.addElement( copy, (0,0,+0.5) ) #render the c++ representation ccomposite = mccomposite.scattererEngine( composite ) ev = mcni.neutron( r = (0,0,-5), v = (0,0,1) ) ccomposite.scatter(ev) return
def test(self): # create a weird shape from mccomposite.geometry import primitives block = primitives.block((1, 1, 1)) sphere = primitives.sphere(1) cylinder = primitives.cylinder(2, 2.001) from mccomposite.geometry import operations dilated = operations.dilate(sphere, 2) translated = operations.translate(block, (0, 0, 0.5)) united = operations.unite(dilated, translated) rotated = operations.rotate(united, (90, 0, 0)) intersect = operations.intersect(rotated, cylinder) difference = operations.subtract(intersect, sphere) print mccomposite.scattererEngine(difference) shape = difference #shape = block #shape = dilated #shape = united #shape = intersect #shape = operations.rotate(block, (90,0,0) ) #shape = rotated #shape = sphere #shape = operations.subtract(sphere, block) #shape = operations.subtract( primitives.cylinder(1, 2.1), sphere ) #create pure python representation of scatterer composite composite = mccomposite.composite(shape) nprinter = NeutronPrinter(shape) composite.addElement(nprinter) #render the c++ representation ccomposite = mccomposite.scattererEngine(composite) ev = mcni.neutron(r=(0, 0, -5), v=(0, 0, 1)) ccomposite.scatter(ev) return
def test(self): # create a weird shape from mccomposite.geometry import primitives block = primitives.block( (1,1,1) ) sphere = primitives.sphere( 1 ) cylinder = primitives.cylinder( 2,2.001 ) from mccomposite.geometry import operations dilated = operations.dilate( sphere, 2 ) translated = operations.translate( block, (0,0,0.5) ) united = operations.unite( dilated, translated ) rotated = operations.rotate( united, (90,0,0) ) intersect = operations.intersect( rotated, cylinder ) difference = operations.subtract( intersect, sphere ) print mccomposite.scattererEngine( difference ) shape = difference #shape = block #shape = dilated #shape = united #shape = intersect #shape = operations.rotate(block, (90,0,0) ) #shape = rotated #shape = sphere #shape = operations.subtract(sphere, block) #shape = operations.subtract( primitives.cylinder(1, 2.1), sphere ) #create pure python representation of scatterer composite composite = mccomposite.composite( shape ) nprinter = NeutronPrinter( shape ) composite.addElement( nprinter ) #render the c++ representation ccomposite = mccomposite.scattererEngine( composite ) ev = mcni.neutron( r = (0,0,-5), v = (0,0,1) ) ccomposite.scatter(ev) return
def test(self): '''HollowCylinder''' print "This test creates a hollow cylinder."\ "It does not interact with neutrons. It prints the info "\ "about the neutrons passing thru them, however. "\ "This test then send one neutron through this hollow cylinder, "\ "and we should see two printings of neutron info, "\ "differing only on time-of-flight." # create a shape shape = HollowCylinder.HollowCylinder( 1, 1.2, 1 ) import UseNeutronPrinter2 nprinter = UseNeutronPrinter2.NeutronPrinter( shape ) #render the c++ representation cinstance = mccomposite.scattererEngine( nprinter ) ev = mcni.neutron( r = (0,0,-5), v = (0,0,1) ) cinstance.scatter(ev) return