1mm液闪的影响

[mirguest]

探测器结构：

• 1 mm LS
• 真空
• PMT
• 钢罐

[mirguest]

标签：issue147-1mm

[mirguest]

非常奇怪，设置PMT外面为Vacuum后，运行run.mac，在Event2处停留很长时间。

[mirguest]

难道是发生了全反射？

[mirguest]

外面是否应该再加有机玻璃罐呢？

[mirguest]

邓老师建议，将外部的材料改成与LS的折射率相同，并且吸收长度设为10000m，避免吸收。

[mirguest]

建立tag：issue147-dummymat-1mm

[mirguest]

现在看来，可以将吸收长度变小一点，避免光子不断再传输。

[mirguest]

或者，我们直接用LS，但是不会发光就可以。

[mirguest]

新建issue147-dummymat-10KM-1mm，吸收长度为10000m。

[mirguest]

又建立了1KM的。不过似乎还是不太乐观。

[mirguest]

可以考虑对长度进行限制，大于多少就直接kill。

[mirguest]

如果对原来的探测器进行考察，使用有1mm薄膜的方案。我们额外记录下在薄膜中的沉积能量。

[mirguest]

这样，就可以额外的得到一个粒子在LS薄膜中和总的LS中的沉积比率。

[mirguest]

dummymat-1mm的作业没有跑完就死掉了。

terminate called after throwing an instance of 'std::bad_alloc'
  what():  St9bad_alloc

[mirguest]

如果设为不发光的LS，应该怎么设置光学属性呢？

[mirguest]

issue147-dummyLS-1mm

[DingXuefeng]

DsPhysConsOptical.cc里面有一些开关，不知道是不是设成false就不发光了。

[DingXuefeng]

如果在DsPhysConsOptical.cc中Line 36设置m_useScintillation=false;，那么在104行就会scint->SetNoOp()， 在DsG4Scintillation.h中这个inline函数会设置m_noop为false，进而在DsG4Scintillation的PostStepDoIt()中line 181会直接return，因而不会有新的光子产生。

这里m_useScintillation为private变量，并且在整个DsPhysConsOptical类中只有构造函数中对它进行了赋值，所以只要在line36处修改它为false，就一定能够关闭闪烁过程

[mirguest]

@DingXuefeng 我现在是希望通过设置Material Property Table来管理发光的属性。你这个应该是对物理过程的修改了吧？

[mirguest]

使用dummy LS后，运行速度果然很快了。

[mirguest]

下面是考虑如何考虑1mm LS的影响。

[mirguest]

我们可以计算LS中沉积能量与该粒子总能量间的关系。

[mirguest]

可以考虑在ROOT中直接增加能量这个Branch。？？？

[mirguest]

直接绘制edep_total vs energy

[mirguest]

我们为了简单，先不考虑事例中不同粒子的混合。直接考虑每个粒子的沉积能量，位置，及能量间的关系。

[mirguest]

K40：

[mirguest]

Uranium：

[mirguest]

但此处能量到8MeV，是不是有问题呢？

[mirguest]

Thorium：

[mirguest]

难道我们的产生子有问题？发现能量有些高。

[mirguest]

因为我们使用的能谱计算程序，包含了所有的粒子。没有拆分。

[mirguest]

需要拆分后看看。

[mirguest]

我前面的做法可能不正确，因为alpha粒子和gamma，e-不同。alpha粒子有淬灭效应的影响。因此，如果它沉积的能量为5MeV，但从闪烁光的效果来看，可能只有0.5MeV了。

[mirguest]

需要重新考虑分析的方法。

[mirguest]

我们应该主要考虑totalPE的影响。

[chenqy]

@mirguest 按照我目前的理解，不如简单一点先试试： （1）用1MeV的gamma作为刻度，得到相应的p.e.数，假设线性，根据得到的root文件里的totalPE，推出来沉积能量，和edep_total作一下比较。我们现在的edep_total是分开记录的，比较一下看看差别多大。 （2）用第一步推出来的能量和文献里Fig7比较一下。 （3）让hitTime<300ns，看看这个条件下的totalPE'推出来的能量和Fig6比较有何异同。

[chenqy]

@mirguest 哦，我好像弄错了，hitTime是光子的时间信息？

[chenqy]

而且还与InitTime有关。。。

[mirguest]

@chenqy 我觉得以后的作业还是把本底事例先进行拆分吧。这样相对来说简单一点。只不过会丢失时间信息。

[mirguest]

后面贴出的是K，U，Th的totalPE。我们将事例按照时间的不同进行拆分。

K

U

Th

[mirguest]

此处，应该可以给出total PE cut的限制。我们可以考虑混合所有的本底，这样就可以设置某阈值后，本底的剩余率。

[mirguest]

不过，如果把三个都放在原点混合，似乎就比较极端了。

[mirguest]

初步分析了Th中的alpha, beta和gamma。alpha粒子会基本沉积。而gamma则飞行的比较远，沉积的比例不是太多。

[mirguest]

Th中沉积能量与动能间的关系

alpha

beta

gamma

[mirguest]

Th光子总数和沉积能量间的关系

alpha

beta

gamma