截面六边形的楔形罩子（test）

chenqy commented 11 years ago

将pmt外部圆柱改成截面六边形的楔形块

chenqy commented 11 years ago

现有的摆法楔形会overlap，需要转动phi? 5pmts

目前的转动是对于一个固定的全局坐标系，所以必须先绕Z轴转动，然后再绕Y，再绕Z转，第一步绕Z转30度如下：

即使转动了，pmt仍旧不够密排。。。需要机械组提供具体的x,y,z坐标和转动角度，必需的是六边形的转动角度！

chenqy commented 11 years ago

刘子然同学的密排法目前未完成：示意图如下 22222222

121212

与五边形相接的六边形形状会略微偏离正六边形有12个五边形，12*5个不是正六边形？

chenqy commented 11 years ago

两段六边楔形，后面可加液闪？ new new2 new3

chenqy commented 11 years ago

修改pmt位置和转动角度的相关代码，使其能读取文件信息：id,x,y,z,r,theta,phi（以前输出过一个文件）。用ID：0-19做测试奇怪，输出的total竟然是21？读取完了可能多读了一行。。。加上旋转角度phi0尝试： 001

chenqy commented 11 years ago

需要机械组提供： id, x, y, z, r, theta, phi, phi0 id: pmt号 x, y, z： pmt坐标 r：摆放球面半径，可无 theta, phi: 相对球心的旋转 phi0: 六边形绕自身对称轴的旋转

chenqy commented 11 years ago

一个面的n个六边形好排，如何弄个正三角形把他们装进去呢？Polyhedra为什么总是偏那么一点呢？偏了之后就不容易得到另外19个面的信息了。。。 tosuccess why

chenqy commented 11 years ago

找到偏的原因了。。。计算错误。。。现在可以学习正二十面体的性质，通过一个基本单元的旋转操作，得到正二十面体了。将来用于刘子然提供的一个球面三角形旋转。 1of20

chenqy commented 11 years ago

200px-icosahedron-golden-rectangles svg

chenqy commented 11 years ago

1,12,5 1,5,6 1,6,9 1,9,2 1,2,12 2,12,8 2,8,7 2,7,9 3,8,11 3,11,4 3,4,10 3,10,7 3,7,8 4,10,6 4,6,5 4,5,11 5,11,12 6,10,9 7,10,9 11,12,8

chenqy commented 11 years ago

如何得到旋转矩阵？ tmp 若已知两个矢量，求解变换矩阵的三个参数不知道行不行。。 http://proj-clhep.web.cern.ch/proj-clhep/manual/UserGuide/VectorDefs/node64.html

chenqy commented 11 years ago

什么是四元数？什么叉乘求旋转轴，点乘、反余弦求夹角。。。

chenqy commented 11 years ago

完成了正二十面体摆放。需要计算参数。。。

chenqy commented 11 years ago

========>chenqy: edge length of icosahedron is 21.7468 m; 22.2304 m; 22.8404 m. ========>chenqy: radius of circumscribed sphere is 20.6825 m; 21.1423 m; 21.7225 m. ========>chenqy: radius of inscribed sphere is 16.4354 m; 16.8009 ; 17261.9 m. 上面的不对。。。 ========>chenqy: edge length of icosahedron (inner) is 21.7468 m. ========>chenqy: radius of circumscribed sphere (outer) is 21.7899 m. ========>chenqy: radius of inscribed sphere is 16.4354 m; 16.8245 ; 17.3154 m. 。。。 ========>chenqy: edge length of icosahedron (inner, middle) is 21.7468 m; 22.2617 m. ========>chenqy: radius of circumscribed sphere (outer) is 21.7899 m. ========>chenqy: radius of inscribed sphere (inner, middle, outer) is 16.4354 m; 16.8245 ; 17.3154 m. ==> x1 = 9713.65; y1 = 9713.65; z1 = 9713.65; r1 = 16824.5; theta1 = 0.955317; phi1 = 0.785398 ==> x2 = 0; y2 = 15717; z2 = 6003.36; r2 = 16824.5; theta2 = 1.20593; phi2 = 1.5708 ==> x3 = -9713.65; y3 = 9713.65; z3 = 9713.65; r3 = 16824.5; theta3 = 0.955317; phi3 = -0.785398 ==> x4 = -6003.36; y4 = 0; z4 = 15717; r4 = 16824.5; theta4 = 0.364864; phi4 = -0 ==> x5 = 6003.36; y5 = 0; z5 = 15717; r5 = 16824.5; theta5 = 0.364864; phi5 = 0 ==> x6 = 9713.65; y6 = -9713.65; z6 = 9713.65; r6 = 16824.5; theta6 = 0.955317; phi6 = -0.785398 ==> x7 = 0; y7 = -15717; z7 = 6003.36; r7 = 16824.5; theta7 = 1.20593; phi7 = -1.5708 ==> x8 = -9713.65; y8 = -9713.65; z8 = 9713.65; r8 = 16824.5; theta8 = 0.955317; phi8 = 0.785398 ==> x9 = 9713.65; y9 = -9713.65; z9 = -9713.65; r9 = 16824.5; theta9 = 2.18628; phi9 = -0.785398 ==> x10 = 6003.36; y10 = 0; z10 = -15717; r10 = 16824.5; theta10 = 2.77673; phi10 = 0 ==> x11 = -6003.36; y11 = 0; z11 = -15717; r11 = 16824.5; theta11 = 2.77673; phi11 = -0 ==> x12 = -9713.65; y12 = -9713.65; z12 = -9713.65; r12 = 16824.5; theta12 = 2.18628; phi12 = 0.785398 ==> x13 = 0; y13 = -15717; z13 = -6003.36; r13 = 16824.5; theta13 = 1.93566; phi13 = -1.5708 ==> x14 = -9713.65; y14 = 9713.65; z14 = -9713.65; r14 = 16824.5; theta14 = 2.18628; phi14 = -0.785398 ==> x15 = 0; y15 = 15717; z15 = -6003.36; r15 = 16824.5; theta15 = 1.93566; phi15 = 1.5708 ==> x16 = 9713.65; y16 = 9713.65; z16 = -9713.65; r16 = 16824.5; theta16 = 2.18628; phi16 = 0.785398 ==> x17 = 15717; y17 = 6003.36; z17 = 0; r17 = 16824.5; theta17 = 1.5708; phi17 = 0.364864 ==> x18 = -15717; y18 = 6003.36; z18 = 0; r18 = 16824.5; theta18 = 1.5708; phi18 = -0.364864 ==> x19 = -15717; y19 = -6003.36; z19 = 0; r19 = 16824.5; theta19 = 1.5708; phi19 = 0.364864 ==> x20 = 15717; y20 = -6003.36; z20 = 0; r20 = 16824.5; theta20 = 1.5708; phi20 = -0.364864

chenqy commented 11 years ago

六边截面柱之间有缝隙，如下：更改三角形尺寸避免三角形overlap：

41 861*20 17220pmts

chenqy commented 11 years ago

有个尖尖出来了。（1）一样大的正二十面体，减少pmt，由41到40,17220到16xxx，太多了。（2）把正二十面体放大，应该比较好。

chenqy commented 11 years ago

新的球面划分？

chenqy commented 11 years ago

[x] 目前20个三角形在平移旋转后有Overlap，可以通过改变三角形的大小等方法解决，如下图，但是需要更详细的参数--------------直接摆放PMT，不用三角形了，免去了考虑三角形overlap的麻烦
[ ] 新的PMT检查无误后，应该用新的PMT
[x] pmtID不对，比如一个面有903个PMT，第二个面如何从903往后计数----------已解决，并且可以不用考虑了
[ ] 刘子然对一个球面三角形的逆施奈德投影还没完成，完成后可以尝试用类似的方法得到所有的球面三角形
[x] 本底产生子还需要考虑----------同上
[x] 若daughter不能超出mother，考虑先用overlap的运行得到数据输入到txt，然后直接用读取的办法摆放，省掉构建不必要的物理体----------同上

chenqy commented 11 years ago

机械组工程上不容易实现是什么意思？做出来大的框架，小的部分，只要能在相应的坐标焊上放PMT的架子就可以了，（球面）三角形怎么分都行的，不是么

新奥尔良城中的名为NOAH (New Orleans Arcology Habitat)或者Arco(short for arcology) noah http://jz.zhulong.com/jz_proj_detail.asp?id=64124