diff --git a/Projects/Vendeta/Analysis.cxx b/Projects/Vendeta/Analysis.cxx index 5eb10ffa62772783514e9dd337e6af58f045e643..c1a183235c1b5164d57a5f6d0eeb88099889c13e 100644 --- a/Projects/Vendeta/Analysis.cxx +++ b/Projects/Vendeta/Analysis.cxx @@ -88,7 +88,7 @@ void Analysis::TreatEvent(){ double DT = Time_Vendeta - Time_FC;// + ToF_Shift_Vendlg[DetNbr-1]; - if(DT>0){ + if(DT>-500){ double DeltaTheta = atan(63.5/Rdet); double Theta_Vendeta = DetPos.Theta(); @@ -125,7 +125,7 @@ void Analysis::TreatEvent(){ double DT = Time_Vendeta - Time_FC;// + ToF_Shift_Vendhg[DetNbr-1]; - if(DT>0){ + if(DT>-500){ double DeltaTheta = atan(63.5/Rdet); double Theta_Vendeta = DetPos.Theta(); double Theta_random = ra.Uniform(Theta_Vendeta-DeltaTheta,Theta_Vendeta+DeltaTheta); diff --git a/Projects/Vendeta/calibration/Vendeta_Time.cal b/Projects/Vendeta/calibration/Vendeta_Time.cal index dbd6a3af84d0c29e489d3daa1062590b40267189..5bd6814598b3a4bd533c791bdf37a39e289f6ec4 100755 --- a/Projects/Vendeta/calibration/Vendeta_Time.cal +++ b/Projects/Vendeta/calibration/Vendeta_Time.cal @@ -1,144 +1,144 @@ -Vendeta_DET1_LG_ANODE6_TIMEOFFSET -59.6746 -Vendeta_DET1_HG_ANODE6_TIMEOFFSET -44.8191 -Vendeta_DET2_LG_ANODE6_TIMEOFFSET -57.8197 -Vendeta_DET2_HG_ANODE6_TIMEOFFSET -49.0571 -Vendeta_DET3_LG_ANODE6_TIMEOFFSET -62.3393 -Vendeta_DET3_HG_ANODE6_TIMEOFFSET -55.258 -Vendeta_DET4_LG_ANODE6_TIMEOFFSET -62.742 -Vendeta_DET4_HG_ANODE6_TIMEOFFSET -59.7595 -Vendeta_DET5_LG_ANODE6_TIMEOFFSET -67.6211 -Vendeta_DET5_HG_ANODE6_TIMEOFFSET -50.2892 -Vendeta_DET6_LG_ANODE6_TIMEOFFSET -60.7652 -Vendeta_DET6_HG_ANODE6_TIMEOFFSET -45.226 -Vendeta_DET7_LG_ANODE6_TIMEOFFSET -62.8352 -Vendeta_DET7_HG_ANODE6_TIMEOFFSET -50.6045 -Vendeta_DET8_LG_ANODE6_TIMEOFFSET -61.0707 -Vendeta_DET8_HG_ANODE6_TIMEOFFSET -48.8541 -Vendeta_DET9_LG_ANODE6_TIMEOFFSET -55.6891 -Vendeta_DET9_HG_ANODE6_TIMEOFFSET -52.6181 -Vendeta_DET10_LG_ANODE6_TIMEOFFSET -56.8744 -Vendeta_DET10_HG_ANODE6_TIMEOFFSET -52.0033 -Vendeta_DET11_LG_ANODE6_TIMEOFFSET -56.35 -Vendeta_DET11_HG_ANODE6_TIMEOFFSET -59.0317 -Vendeta_DET12_LG_ANODE6_TIMEOFFSET -61.8728 -Vendeta_DET12_HG_ANODE6_TIMEOFFSET -60.5685 -Vendeta_DET13_LG_ANODE6_TIMEOFFSET -58.6669 -Vendeta_DET13_HG_ANODE6_TIMEOFFSET -54.4535 -Vendeta_DET14_LG_ANODE6_TIMEOFFSET -55.848 -Vendeta_DET14_HG_ANODE6_TIMEOFFSET -51.7572 -Vendeta_DET15_LG_ANODE6_TIMEOFFSET -59.6486 -Vendeta_DET15_HG_ANODE6_TIMEOFFSET -51.003 -Vendeta_DET16_LG_ANODE6_TIMEOFFSET -61.0214 -Vendeta_DET16_HG_ANODE6_TIMEOFFSET -52.5003 -Vendeta_DET17_LG_ANODE6_TIMEOFFSET -63.2736 -Vendeta_DET17_HG_ANODE6_TIMEOFFSET -56.756 -Vendeta_DET18_LG_ANODE6_TIMEOFFSET -65.1051 -Vendeta_DET18_HG_ANODE6_TIMEOFFSET -56.6478 -Vendeta_DET19_LG_ANODE6_TIMEOFFSET -60.6122 -Vendeta_DET19_HG_ANODE6_TIMEOFFSET -53.7386 -Vendeta_DET20_LG_ANODE6_TIMEOFFSET -60.5737 -Vendeta_DET20_HG_ANODE6_TIMEOFFSET -53.6525 -Vendeta_DET21_LG_ANODE6_TIMEOFFSET -57.2636 +Vendeta_DET1_LG_ANODE6_TIMEOFFSET -51.7104 +Vendeta_DET1_HG_ANODE6_TIMEOFFSET -40.9152 +Vendeta_DET2_LG_ANODE6_TIMEOFFSET -53.8816 +Vendeta_DET2_HG_ANODE6_TIMEOFFSET -45.077 +Vendeta_DET3_LG_ANODE6_TIMEOFFSET -50.2365 +Vendeta_DET3_HG_ANODE6_TIMEOFFSET -51.2556 +Vendeta_DET4_LG_ANODE6_TIMEOFFSET -58.7939 +Vendeta_DET4_HG_ANODE6_TIMEOFFSET -55.7588 +Vendeta_DET5_LG_ANODE6_TIMEOFFSET -3.28871 +Vendeta_DET5_HG_ANODE6_TIMEOFFSET -50.2452 +Vendeta_DET6_LG_ANODE6_TIMEOFFSET 3.54608 +Vendeta_DET6_HG_ANODE6_TIMEOFFSET -45.2883 +Vendeta_DET7_LG_ANODE6_TIMEOFFSET 5.4746 +Vendeta_DET7_HG_ANODE6_TIMEOFFSET -46.5201 +Vendeta_DET8_LG_ANODE6_TIMEOFFSET 7.30702 +Vendeta_DET8_HG_ANODE6_TIMEOFFSET -44.8232 +Vendeta_DET9_LG_ANODE6_TIMEOFFSET 8.65419 +Vendeta_DET9_HG_ANODE6_TIMEOFFSET -44.6394 +Vendeta_DET10_LG_ANODE6_TIMEOFFSET 7.42644 +Vendeta_DET10_HG_ANODE6_TIMEOFFSET -44.008 +Vendeta_DET11_LG_ANODE6_TIMEOFFSET 8.0436 +Vendeta_DET11_HG_ANODE6_TIMEOFFSET -55.0314 +Vendeta_DET12_LG_ANODE6_TIMEOFFSET 10.4369 +Vendeta_DET12_HG_ANODE6_TIMEOFFSET -56.5682 +Vendeta_DET13_LG_ANODE6_TIMEOFFSET 5.65558 +Vendeta_DET13_HG_ANODE6_TIMEOFFSET -50.5655 +Vendeta_DET14_LG_ANODE6_TIMEOFFSET 8.43334 +Vendeta_DET14_HG_ANODE6_TIMEOFFSET -47.7949 +Vendeta_DET15_LG_ANODE6_TIMEOFFSET 0.650771 +Vendeta_DET15_HG_ANODE6_TIMEOFFSET -51.1356 +Vendeta_DET16_LG_ANODE6_TIMEOFFSET 3.33674 +Vendeta_DET16_HG_ANODE6_TIMEOFFSET -52.4834 +Vendeta_DET17_LG_ANODE6_TIMEOFFSET 4.8467 +Vendeta_DET17_HG_ANODE6_TIMEOFFSET -48.8195 +Vendeta_DET18_LG_ANODE6_TIMEOFFSET 3.04881 +Vendeta_DET18_HG_ANODE6_TIMEOFFSET -48.6242 +Vendeta_DET19_LG_ANODE6_TIMEOFFSET 3.72647 +Vendeta_DET19_HG_ANODE6_TIMEOFFSET -53.6117 +Vendeta_DET20_LG_ANODE6_TIMEOFFSET 7.71313 +Vendeta_DET20_HG_ANODE6_TIMEOFFSET -53.6624 +Vendeta_DET21_LG_ANODE6_TIMEOFFSET -53.2306 Vendeta_DET21_HG_ANODE6_TIMEOFFSET 0 -Vendeta_DET22_LG_ANODE6_TIMEOFFSET -53.1483 -Vendeta_DET22_HG_ANODE6_TIMEOFFSET -47.8539 -Vendeta_DET23_LG_ANODE6_TIMEOFFSET -54.9371 -Vendeta_DET23_HG_ANODE6_TIMEOFFSET -49.2088 -Vendeta_DET24_LG_ANODE6_TIMEOFFSET -57.0301 -Vendeta_DET24_HG_ANODE6_TIMEOFFSET -49.1904 -Vendeta_DET25_LG_ANODE6_TIMEOFFSET -61.0623 -Vendeta_DET25_HG_ANODE6_TIMEOFFSET -46.9562 -Vendeta_DET26_LG_ANODE6_TIMEOFFSET -63.0825 -Vendeta_DET26_HG_ANODE6_TIMEOFFSET -49.0353 -Vendeta_DET27_LG_ANODE6_TIMEOFFSET -58.8964 -Vendeta_DET27_HG_ANODE6_TIMEOFFSET -48.3858 -Vendeta_DET28_LG_ANODE6_TIMEOFFSET -62.5929 -Vendeta_DET28_HG_ANODE6_TIMEOFFSET -56.2581 -Vendeta_DET29_LG_ANODE6_TIMEOFFSET -55.7432 -Vendeta_DET29_HG_ANODE6_TIMEOFFSET -53.4057 -Vendeta_DET30_LG_ANODE6_TIMEOFFSET -57.9021 -Vendeta_DET30_HG_ANODE6_TIMEOFFSET -51.8945 -Vendeta_DET31_LG_ANODE6_TIMEOFFSET -54.7518 -Vendeta_DET31_HG_ANODE6_TIMEOFFSET -54.2489 -Vendeta_DET32_LG_ANODE6_TIMEOFFSET -51.8877 -Vendeta_DET32_HG_ANODE6_TIMEOFFSET -140.63 -Vendeta_DET33_LG_ANODE6_TIMEOFFSET -63.1306 -Vendeta_DET33_HG_ANODE6_TIMEOFFSET -49.753 -Vendeta_DET34_LG_ANODE6_TIMEOFFSET -64.3207 -Vendeta_DET34_HG_ANODE6_TIMEOFFSET -55.0256 -Vendeta_DET35_LG_ANODE6_TIMEOFFSET -61.6654 -Vendeta_DET35_HG_ANODE6_TIMEOFFSET -56.0665 -Vendeta_DET36_LG_ANODE6_TIMEOFFSET -64.8874 -Vendeta_DET36_HG_ANODE6_TIMEOFFSET -55.3767 -Vendeta_DET37_LG_ANODE6_TIMEOFFSET -56.513 -Vendeta_DET37_HG_ANODE6_TIMEOFFSET -43.5398 -Vendeta_DET38_LG_ANODE6_TIMEOFFSET -69.4906 -Vendeta_DET38_HG_ANODE6_TIMEOFFSET -50.343 -Vendeta_DET39_LG_ANODE6_TIMEOFFSET -63.5212 -Vendeta_DET39_HG_ANODE6_TIMEOFFSET -46.0627 -Vendeta_DET40_LG_ANODE6_TIMEOFFSET -59.6885 -Vendeta_DET40_HG_ANODE6_TIMEOFFSET -46.179 -Vendeta_DET41_LG_ANODE6_TIMEOFFSET -62.0332 -Vendeta_DET41_HG_ANODE6_TIMEOFFSET -51.9157 -Vendeta_DET42_LG_ANODE6_TIMEOFFSET -65.0402 -Vendeta_DET42_HG_ANODE6_TIMEOFFSET -51.0274 -Vendeta_DET43_LG_ANODE6_TIMEOFFSET -49.6916 -Vendeta_DET43_HG_ANODE6_TIMEOFFSET -49.4824 -Vendeta_DET44_LG_ANODE6_TIMEOFFSET -58.8582 -Vendeta_DET44_HG_ANODE6_TIMEOFFSET -56.4819 -Vendeta_DET45_LG_ANODE6_TIMEOFFSET -59.3703 -Vendeta_DET45_HG_ANODE6_TIMEOFFSET -50.2117 -Vendeta_DET46_LG_ANODE6_TIMEOFFSET -59.4935 -Vendeta_DET46_HG_ANODE6_TIMEOFFSET -50.5453 -Vendeta_DET47_LG_ANODE6_TIMEOFFSET -53.8864 -Vendeta_DET47_HG_ANODE6_TIMEOFFSET -54.4564 -Vendeta_DET48_LG_ANODE6_TIMEOFFSET -56.2152 -Vendeta_DET48_HG_ANODE6_TIMEOFFSET -52.7713 -Vendeta_DET49_LG_ANODE6_TIMEOFFSET -61.4114 -Vendeta_DET49_HG_ANODE6_TIMEOFFSET -50.4108 -Vendeta_DET50_LG_ANODE6_TIMEOFFSET -52.6468 -Vendeta_DET50_HG_ANODE6_TIMEOFFSET -45.3396 -Vendeta_DET51_LG_ANODE6_TIMEOFFSET -67.2692 -Vendeta_DET51_HG_ANODE6_TIMEOFFSET -51.5164 -Vendeta_DET52_LG_ANODE6_TIMEOFFSET -67.5185 -Vendeta_DET52_HG_ANODE6_TIMEOFFSET -51.6829 -Vendeta_DET53_LG_ANODE6_TIMEOFFSET -61.1497 -Vendeta_DET53_HG_ANODE6_TIMEOFFSET -48.5667 -Vendeta_DET54_LG_ANODE6_TIMEOFFSET -63.3634 -Vendeta_DET54_HG_ANODE6_TIMEOFFSET -50.8508 -Vendeta_DET55_LG_ANODE6_TIMEOFFSET -58.1545 -Vendeta_DET55_HG_ANODE6_TIMEOFFSET -57.2594 -Vendeta_DET56_LG_ANODE6_TIMEOFFSET -58.9024 -Vendeta_DET56_HG_ANODE6_TIMEOFFSET -55.8663 -Vendeta_DET57_LG_ANODE6_TIMEOFFSET -58.0078 -Vendeta_DET57_HG_ANODE6_TIMEOFFSET -67.2925 -Vendeta_DET58_LG_ANODE6_TIMEOFFSET -54.1625 -Vendeta_DET58_HG_ANODE6_TIMEOFFSET -43.688 -Vendeta_DET59_LG_ANODE6_TIMEOFFSET -52.3144 -Vendeta_DET59_HG_ANODE6_TIMEOFFSET -50.1517 -Vendeta_DET60_LG_ANODE6_TIMEOFFSET -71.1129 -Vendeta_DET60_HG_ANODE6_TIMEOFFSET -89.1563 -Vendeta_DET61_LG_ANODE6_TIMEOFFSET -58.3139 -Vendeta_DET61_HG_ANODE6_TIMEOFFSET -63.612 -Vendeta_DET62_LG_ANODE6_TIMEOFFSET -67.4705 -Vendeta_DET62_HG_ANODE6_TIMEOFFSET -296.824 -Vendeta_DET63_LG_ANODE6_TIMEOFFSET -60.612 -Vendeta_DET63_HG_ANODE6_TIMEOFFSET -48.091 -Vendeta_DET64_LG_ANODE6_TIMEOFFSET -54.2712 -Vendeta_DET64_HG_ANODE6_TIMEOFFSET -131.131 -Vendeta_DET65_LG_ANODE6_TIMEOFFSET -54.1577 -Vendeta_DET65_HG_ANODE6_TIMEOFFSET -53.6022 -Vendeta_DET66_LG_ANODE6_TIMEOFFSET -50.9856 -Vendeta_DET66_HG_ANODE6_TIMEOFFSET -119.727 -Vendeta_DET67_LG_ANODE6_TIMEOFFSET -62.2156 -Vendeta_DET67_HG_ANODE6_TIMEOFFSET -55.7379 -Vendeta_DET68_LG_ANODE6_TIMEOFFSET -50.3268 -Vendeta_DET68_HG_ANODE6_TIMEOFFSET -55.35 -Vendeta_DET69_LG_ANODE6_TIMEOFFSET -54.2274 -Vendeta_DET69_HG_ANODE6_TIMEOFFSET -49.0134 -Vendeta_DET70_LG_ANODE6_TIMEOFFSET -56.6098 -Vendeta_DET70_HG_ANODE6_TIMEOFFSET -51.9793 -Vendeta_DET71_LG_ANODE6_TIMEOFFSET -55.9297 -Vendeta_DET71_HG_ANODE6_TIMEOFFSET -93.1255 -Vendeta_DET72_LG_ANODE6_TIMEOFFSET -60.1854 -Vendeta_DET72_HG_ANODE6_TIMEOFFSET -51.1992 +Vendeta_DET22_LG_ANODE6_TIMEOFFSET -53.1936 +Vendeta_DET22_HG_ANODE6_TIMEOFFSET -43.976 +Vendeta_DET23_LG_ANODE6_TIMEOFFSET -62.8682 +Vendeta_DET23_HG_ANODE6_TIMEOFFSET -49.2518 +Vendeta_DET24_LG_ANODE6_TIMEOFFSET -57.1077 +Vendeta_DET24_HG_ANODE6_TIMEOFFSET -49.368 +Vendeta_DET25_LG_ANODE6_TIMEOFFSET -57.11 +Vendeta_DET25_HG_ANODE6_TIMEOFFSET -51.0045 +Vendeta_DET26_LG_ANODE6_TIMEOFFSET -55.1781 +Vendeta_DET26_HG_ANODE6_TIMEOFFSET -53.0488 +Vendeta_DET27_LG_ANODE6_TIMEOFFSET -54.9169 +Vendeta_DET27_HG_ANODE6_TIMEOFFSET -40.474 +Vendeta_DET28_LG_ANODE6_TIMEOFFSET -62.7126 +Vendeta_DET28_HG_ANODE6_TIMEOFFSET -48.2743 +Vendeta_DET29_LG_ANODE6_TIMEOFFSET -56.6157 +Vendeta_DET29_HG_ANODE6_TIMEOFFSET -45.8951 +Vendeta_DET30_LG_ANODE6_TIMEOFFSET -49.9337 +Vendeta_DET30_HG_ANODE6_TIMEOFFSET -41.2286 +Vendeta_DET31_LG_ANODE6_TIMEOFFSET -59.4607 +Vendeta_DET31_HG_ANODE6_TIMEOFFSET -46.3677 +Vendeta_DET32_LG_ANODE6_TIMEOFFSET -56.5875 +Vendeta_DET32_HG_ANODE6_TIMEOFFSET -47.4646 +Vendeta_DET33_LG_ANODE6_TIMEOFFSET -55.0935 +Vendeta_DET33_HG_ANODE6_TIMEOFFSET -53.8766 +Vendeta_DET34_LG_ANODE6_TIMEOFFSET -56.4013 +Vendeta_DET34_HG_ANODE6_TIMEOFFSET -51.0475 +Vendeta_DET35_LG_ANODE6_TIMEOFFSET -61.7078 +Vendeta_DET35_HG_ANODE6_TIMEOFFSET -48.2275 +Vendeta_DET36_LG_ANODE6_TIMEOFFSET -60.9581 +Vendeta_DET36_HG_ANODE6_TIMEOFFSET -47.5593 +Vendeta_DET37_LG_ANODE6_TIMEOFFSET -52.6198 +Vendeta_DET37_HG_ANODE6_TIMEOFFSET -43.5903 +Vendeta_DET38_LG_ANODE6_TIMEOFFSET -57.6652 +Vendeta_DET38_HG_ANODE6_TIMEOFFSET -50.6965 +Vendeta_DET39_LG_ANODE6_TIMEOFFSET -55.6432 +Vendeta_DET39_HG_ANODE6_TIMEOFFSET -46.2478 +Vendeta_DET40_LG_ANODE6_TIMEOFFSET -51.7419 +Vendeta_DET40_HG_ANODE6_TIMEOFFSET -46.3196 +Vendeta_DET41_LG_ANODE6_TIMEOFFSET -62.0913 +Vendeta_DET41_HG_ANODE6_TIMEOFFSET -48.0791 +Vendeta_DET42_LG_ANODE6_TIMEOFFSET -57.0765 +Vendeta_DET42_HG_ANODE6_TIMEOFFSET -47.1489 +Vendeta_DET43_LG_ANODE6_TIMEOFFSET -54.0486 +Vendeta_DET43_HG_ANODE6_TIMEOFFSET -45.6618 +Vendeta_DET44_LG_ANODE6_TIMEOFFSET -55.0131 +Vendeta_DET44_HG_ANODE6_TIMEOFFSET -52.7097 +Vendeta_DET45_LG_ANODE6_TIMEOFFSET -59.4656 +Vendeta_DET45_HG_ANODE6_TIMEOFFSET -54.4277 +Vendeta_DET46_LG_ANODE6_TIMEOFFSET -59.7746 +Vendeta_DET46_HG_ANODE6_TIMEOFFSET -54.7329 +Vendeta_DET47_LG_ANODE6_TIMEOFFSET -58.192 +Vendeta_DET47_HG_ANODE6_TIMEOFFSET -48.8151 +Vendeta_DET48_LG_ANODE6_TIMEOFFSET -52.4046 +Vendeta_DET48_HG_ANODE6_TIMEOFFSET -49.042 +Vendeta_DET49_LG_ANODE6_TIMEOFFSET -1.19671 +Vendeta_DET49_HG_ANODE6_TIMEOFFSET -46.533 +Vendeta_DET50_LG_ANODE6_TIMEOFFSET 3.66615 +Vendeta_DET50_HG_ANODE6_TIMEOFFSET -41.6431 +Vendeta_DET51_LG_ANODE6_TIMEOFFSET 0 +Vendeta_DET51_HG_ANODE6_TIMEOFFSET -47.5382 +Vendeta_DET52_LG_ANODE6_TIMEOFFSET 0.553368 +Vendeta_DET52_HG_ANODE6_TIMEOFFSET -47.7952 +Vendeta_DET53_LG_ANODE6_TIMEOFFSET -49.2352 +Vendeta_DET53_HG_ANODE6_TIMEOFFSET -48.9065 +Vendeta_DET54_LG_ANODE6_TIMEOFFSET -55.485 +Vendeta_DET54_HG_ANODE6_TIMEOFFSET -51.0403 +Vendeta_DET55_LG_ANODE6_TIMEOFFSET -54.3748 +Vendeta_DET55_HG_ANODE6_TIMEOFFSET -45.3554 +Vendeta_DET56_LG_ANODE6_TIMEOFFSET -51.0772 +Vendeta_DET56_HG_ANODE6_TIMEOFFSET -44.1187 +Vendeta_DET57_LG_ANODE6_TIMEOFFSET -59.008 +Vendeta_DET57_HG_ANODE6_TIMEOFFSET -52.1784 +Vendeta_DET58_LG_ANODE6_TIMEOFFSET -54.9421 +Vendeta_DET58_HG_ANODE6_TIMEOFFSET -48.062 +Vendeta_DET59_LG_ANODE6_TIMEOFFSET -53.0501 +Vendeta_DET59_HG_ANODE6_TIMEOFFSET -49.7265 +Vendeta_DET60_LG_ANODE6_TIMEOFFSET -59.0688 +Vendeta_DET60_HG_ANODE6_TIMEOFFSET -55.7496 +Vendeta_DET61_LG_ANODE6_TIMEOFFSET -62.9952 +Vendeta_DET61_HG_ANODE6_TIMEOFFSET -53.1383 +Vendeta_DET62_LG_ANODE6_TIMEOFFSET -63.9224 +Vendeta_DET62_HG_ANODE6_TIMEOFFSET -58.0663 +Vendeta_DET63_LG_ANODE6_TIMEOFFSET -56.604 +Vendeta_DET63_HG_ANODE6_TIMEOFFSET -44.1934 +Vendeta_DET64_LG_ANODE6_TIMEOFFSET -59.142 +Vendeta_DET64_HG_ANODE6_TIMEOFFSET -42.7817 +Vendeta_DET65_LG_ANODE6_TIMEOFFSET -50.8804 +Vendeta_DET65_HG_ANODE6_TIMEOFFSET -41.9318 +Vendeta_DET66_LG_ANODE6_TIMEOFFSET -55.8415 +Vendeta_DET66_HG_ANODE6_TIMEOFFSET -44.9224 +Vendeta_DET67_LG_ANODE6_TIMEOFFSET -59.3876 +Vendeta_DET67_HG_ANODE6_TIMEOFFSET -52.1028 +Vendeta_DET68_LG_ANODE6_TIMEOFFSET -58.9977 +Vendeta_DET68_HG_ANODE6_TIMEOFFSET -51.7337 +Vendeta_DET69_LG_ANODE6_TIMEOFFSET -55.1883 +Vendeta_DET69_HG_ANODE6_TIMEOFFSET -44.9601 +Vendeta_DET70_LG_ANODE6_TIMEOFFSET -62.2589 +Vendeta_DET70_HG_ANODE6_TIMEOFFSET -48.1123 +Vendeta_DET71_LG_ANODE6_TIMEOFFSET -58.7874 +Vendeta_DET71_HG_ANODE6_TIMEOFFSET -51.9722 +Vendeta_DET72_LG_ANODE6_TIMEOFFSET -52.3003 +Vendeta_DET72_HG_ANODE6_TIMEOFFSET -51.6784 diff --git a/Projects/Vendeta/macro/FillTOFHisto.C b/Projects/Vendeta/macro/FillTOFHisto.C new file mode 100644 index 0000000000000000000000000000000000000000..171b775bcc23f406d27e46f4f912d88a6a35a284 --- /dev/null +++ b/Projects/Vendeta/macro/FillTOFHisto.C @@ -0,0 +1,104 @@ +TChain* chain; + +int NumberOfDetectors= 72; +int NumberOfAnodes= 1; +int nentries=1e6; + +///////////////////////////////////// +void LoadRootFile(){ + chain = new TChain("PhysicsTree"); + chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_1.root"); + //chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_2.root"); + //chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_3.root"); + //chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_4.root"); +} + +///////////////////////////////////// +void FillTOFHisto(){ + + LoadRootFile(); + nentries = chain->GetEntries(); + cout << "Number of entries: " << nentries << endl; + + TFile* ofile = new TFile("histo_tof_file.root","recreate"); + TH1F* hLG[792]; + TH1F* hHG[792]; + + vector<double>* FC_Q1 = new vector<double>(); + + vector<double>* LG_Tof = new vector<double>(); + vector<int>* LG_ID = new vector<int>(); + vector<int>* LG_Anode_ID = new vector<int>(); + + vector<double>* HG_Tof = new vector<double>(); + vector<int>* HG_ID = new vector<int>(); + vector<int>* HG_Anode_ID = new vector<int>(); + + TFissionChamberPhysics* FC = new TFissionChamberPhysics(); + chain->SetBranchAddress("FissionChamber",&FC); + + chain->SetBranchAddress("FC_Q1",&FC_Q1); + chain->SetBranchAddress("LG_Tof",&LG_Tof); + chain->SetBranchAddress("LG_ID",&LG_ID); + chain->SetBranchAddress("LG_Anode_ID",&LG_Anode_ID); + chain->SetBranchAddress("HG_Tof",&HG_Tof); + chain->SetBranchAddress("HG_ID",&HG_ID); + chain->SetBranchAddress("HG_Anode_ID",&HG_Anode_ID); + +for(int i=0; i<NumberOfDetectors; i++){ + for(int j=0; j<NumberOfAnodes; j++){ + int anode = 6; + //int index = (i+1) * (j+1); + int index = (i+1) * anode; + //TString histo_name = Form("hLG_Det%i_Anode%i",i+1,j+1); + TString histo_name = Form("hLG_Det%i_Anode%i",i+1,anode); + hLG[index-1] = new TH1F(histo_name,histo_name,4000,-100,300); + //histo_name = Form("hHG_Det%i_Anode%i",i+1,j+1); + histo_name = Form("hHG_Det%i_Anode%i",i+1,anode); + hHG[index-1] = new TH1F(histo_name,histo_name,3000,0,300); + } + } + + for(int i=0; i<nentries; i++){ + chain->GetEntry(i); + + if(i%100000==0){ + cout << "\033[34m\r Processing tree..." << (double)i/nentries*100 << "\% done" << flush; + } + //int FC_mult = FC->AnodeNumber.size(); + //cout << FC_mult << endl; + int mysize = LG_Tof->size(); + for(int j=0; j<mysize; j++){ + // LG // + int index_LG = LG_ID->at(j) * LG_Anode_ID->at(j); + if(LG_ID->at(j)>0 && LG_Anode_ID->at(j)>0){ + hLG[index_LG-1]->Fill(LG_Tof->at(j)); + } + } + + mysize = HG_Tof->size(); + for(int j=0; j<mysize; j++){ + // HG // + int index_HG = HG_ID->at(j) * HG_Anode_ID->at(j); + if(HG_ID->at(j)>0 && HG_Anode_ID->at(j)>0){ + hHG[index_HG-1]->Fill(HG_Tof->at(j)); + } + } + } + + /* for(int i=0; i<NumberOfDetectors; i++){ */ + /* for(int j=0; j<NumberOfAnodes; j++){ */ + /* int anode = 6; */ + /* //int index = (i+1) * (j+1); */ + /* int index = (i+1) * anode; */ + + /* hLG[index]->Write(); */ + /* hHG[index]->Write(); */ + /* } */ + /* } */ + + //hLG[6]->Draw(); + ofile->Write(); + ofile->Close(); + +} diff --git a/Projects/Vendeta/macro/FitTofGammaPeak.C b/Projects/Vendeta/macro/FitTofGammaPeak.C index 061b71adc19dce179dca47abe11583b357f6538e..610acbe975c878435a272ad0682541438d557793 100644 --- a/Projects/Vendeta/macro/FitTofGammaPeak.C +++ b/Projects/Vendeta/macro/FitTofGammaPeak.C @@ -10,7 +10,8 @@ bool Finder(TH1F* h, Double_t *mean, Double_t *sigma); ///////////////////////////////////////////////////// void OpenRootFile(){ - ifile = new TFile("histo_tof_file_q1_80ns.root"); + //ifile = new TFile("histo_tof_file_q1_80ns.root"); + ifile = new TFile("histo_tof_file.root"); } ///////////////////////////////////////////////////// diff --git a/Projects/Vendeta/macro/GenerateTOFHisto.C b/Projects/Vendeta/macro/GenerateTOFHisto.C index 6a66a83d3b06b818a0a55fa9720d8cb3019e4110..34eecb304802851cd92fb0bb6043dc01edd9b743 100644 --- a/Projects/Vendeta/macro/GenerateTOFHisto.C +++ b/Projects/Vendeta/macro/GenerateTOFHisto.C @@ -2,24 +2,30 @@ TChain* chain; int NumberOfDetectors= 72; int NumberOfAnodes= 1; - +int nentries=1e6; ///////////////////////////////////// void LoadRootFile(){ chain = new TChain("PhysicsTree"); - chain->Add("/home/faster/nptool/Outputs/Analysis/test_cf4_q1_160ns.root"); + //chain->Add("/home/faster/nptool/Outputs/Analysis/test_cf4_q1_160ns.root"); + chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_1.root"); + chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_2.root"); + chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_3.root"); + chain->Add("/home/faster/nptool/Outputs/Analysis/test_sampler_qdc_cf_4.root"); } ///////////////////////////////////// void GenerateTOFHisto(){ LoadRootFile(); - TFile* ofile = new TFile("histo_tof_file_q1_160ns.root","recreate"); + nentries = chain->GetEntries(); + TFile* ofile = new TFile("histo_tof_file.root","recreate"); + for(int i=0; i<NumberOfDetectors; i++){ for(int j=0; j<NumberOfAnodes; j++){ j=5; // LG // - TString to_draw = Form("LG_Tof>>hLG_Det%i_Anode%i(3000,0,300)",i+1,j+1); + TString to_draw = Form("LG_Tof>>hLG_Det%i_Anode%i(5000,-200,300)",i+1,j+1); TString condition = Form("LG_ID==%i && LG_Anode_ID==%i && FC_Q1>5500",i+1,j+1); TString histo_name_LG = Form("hLG_Det%i_Anode%i",i+1,j+1); chain->Draw(to_draw,condition,"",1e6); diff --git a/Projects/Vendeta/macro/RawSpectra/FillRawFCSpectra.C b/Projects/Vendeta/macro/RawSpectra/FillRawFCSpectra.C new file mode 100644 index 0000000000000000000000000000000000000000..cb4fd2e44a5bb22b9762db2b3ac53aa6ab5ab950 --- /dev/null +++ b/Projects/Vendeta/macro/RawSpectra/FillRawFCSpectra.C @@ -0,0 +1,66 @@ +TChain* chain; + +int NumberOfAnodes= 1; +int entries= 1e6; + +////////////////////////////////////////////////// +void OpenRootFile(){ + chain = new TChain("RawTree"); + chain->Add("/home/faster/fastertonptool/data/rootfiles/V4B_SAMPLING_6_0001.root"); +} + +////////////////////////////////////////////////// +void ShowRawFCSpectra(string Nucleus){ + OpenRootFile(); + nentries = chain->GetEntries(); + + TFile* ofile = new TFile("FC_Raw_spectra.root","recreate"); + TH1F* Q1[11]; + TH1F* Q2[11]; + TH1F* Qmax[11]; + TH1F* Q2vsQ1[11]; + TH1F* QmaxvsQ1[11]; + + TFissionChamberData* FC = new TFissionChamberDaat(); + chain->SetBranchAddress("FissionChamber",&FC); + + TString histo_name; + for(int i=0; i< NumberOfAnodes; i++){ + histo_name = Form("Q1_Anode%i",i+1); + Q1[i] = new TH1F(histo_name, histo_name, 500,0,100e3); + + histo_name = Form("Q2_Anode%i",i+1); + Q2[i] = new TH1F(histo_name, histo_name, 500,0,20e3); + + histo_name = Form("Qmax_Anode%i",i+1); + Qmax[i] = new TH1F(histo_name, histo_name, 500,0,10e3); + + histo_name = Form("Q2vsQ1_Anode%i",i+1); + Q2vsQ1[i] = new TH2F(histo_name, histo_name, 500,0,100e3,500,0,20e3); + + histo_name = Form("QmaxvsQ1_Anode%i",i+1); + QmaxvsQ1[i] = new TH2F(histo_name, histo_name, 500,0,100e3,500,0,10e3); + } + + for(int i=0; i<nentries; i++){ + chain->GetEntry(i); + + int mult = FC->GetMultiplicity(); + for(int j=0; i<mult; j++){ + int anode = FC->GetAnodeNbr(j); + double Q1val = FC->GetQ1(j); + double Q2val = FC->GetQ2(j); + double Qmaxval = FC->GetQmax(j); + + Q1[anode-1]->Fill(Q1val); + Q2[anode-1]->Fill(Q2val); + Qmax[anode-1]->Fill(Qmaxval); + Q2vsQ1[anode-1]->Fill(Q1,Q2); + QmaxvsQ1[anode-1]->Fill(Q1,Q2); + + } + } + + ofile->Write(); + ofile->Close(); +} diff --git a/Projects/Vendeta/macro/RawSpectra/FillRawPSD.C b/Projects/Vendeta/macro/RawSpectra/FillRawPSD.C new file mode 100644 index 0000000000000000000000000000000000000000..f79733f9ddceb470c91be673cd14fb2a4e73c35f --- /dev/null +++ b/Projects/Vendeta/macro/RawSpectra/FillRawPSD.C @@ -0,0 +1,72 @@ +TChain* chain; + +int NumberOfDetectors= 72; +int nentries=1e6; + +////////////////////////////////////////////////// +void OpenRootFile(){ + chain = new TChain("RawTree"); + chain->Add("/home/faster/fastertonptool/data/rootfiles/test_sampler_qdc_cf_1.root"); + chain->Add("/home/faster/fastertonptool/data/rootfiles/test_sampler_qdc_cf_2.root"); + chain->Add("/home/faster/fastertonptool/data/rootfiles/test_sampler_qdc_cf_3.root"); + chain->Add("/home/faster/fastertonptool/data/rootfiles/test_sampler_qdc_cf_4.root"); +} + +////////////////////////////////////////////////// +void FillRawPSD(){ + OpenRootFile(); + nentries = chain->GetEntries(); + + TFile* ofile = new TFile("PSD_histo.root","recreate"); + TH2F* hLG[72]; + TH2F* hHG[72]; + + TString histo_name; + for(int i=0; i<NumberOfDetectors;i++){ + histo_name = Form("hLG_det%i",i+1); + hLG[i] = new TH2F(histo_name,histo_name,500,0,500e3,500,0,1.1); + histo_name = Form("hHG_det%i",i+1); + hHG[i] = new TH2F(histo_name,histo_name,500,0,900e3,500,0,1.1); + } + + TVendetaData* Vendeta = new TVendetaData(); + chain->SetBranchAddress("Vendeta",&Vendeta); + + for(int i=0; i<nentries; i++){ + chain->GetEntry(i); + if(i%100000==0){ + cout << "\033[34m\r Processing tree..." << (double)i/nentries*100 << "\% done" << flush; + } + + // LG // + int mult_LG = Vendeta->GetLGMultEnergy(); + if(mult_LG>0){ + for(int j=0; j<mult_LG; j++){ + double Q1 = Vendeta->GetLGQ1(j); + double Q2 = Vendeta->GetLGQ2(j); + double PSD_LG = Q2/Q1; + int det = Vendeta->GetLGDetectorNbr(j); + + hLG[det-1]->Fill(Q1,PSD_LG); + } + } + + // HG // + int mult_HG = Vendeta->GetHGMultEnergy(); + if(mult_HG>0){ + for(int j=0; j<mult_HG; j++){ + double Q1 = Vendeta->GetHGQ1(j); + double Q2 = Vendeta->GetHGQ2(j); + double PSD_HG = Q2/Q1; + int det = Vendeta->GetHGDetectorNbr(j); + + hHG[det-1]->Fill(Q1,PSD_HG); + } + } + + + } + + ofile->Write(); + ofile->Close(); +} diff --git a/Projects/Vendeta/macro/RawSpectra/ShowRawPSD.C b/Projects/Vendeta/macro/RawSpectra/ShowRawPSD.C index 8dd2b98161fdcfa61e8f23672570583ff5672eba..505d31af419e395abbb3e266be28f931473ae813 100644 --- a/Projects/Vendeta/macro/RawSpectra/ShowRawPSD.C +++ b/Projects/Vendeta/macro/RawSpectra/ShowRawPSD.C @@ -8,14 +8,15 @@ int NumberOfEvents= 1e7; void OpenRootFile(){ chain = new TChain("RawTree"); //chain->Add("/home/faster/fastertonptool/data/rootfiles/V4B_SAMPLING_6_0001.root"); - chain->Add("/home/faster/fastertonptool/data/rootfiles/test_cf4_q1_160ns.root"); + //chain->Add("/home/faster/fastertonptool/data/rootfiles/test_cf4_q1_160ns.root"); + chain->Add("/home/faster/fastertonptool/data/rootfiles/test_sampler_qdc_cf_1.root"); } ////////////////////////////////////////////////// void ShowRawPSD(){ OpenRootFile(); - TFile* ofile = new TFile("PSD_Q1_160ns.root","recreate"); + TFile* ofile = new TFile("PSD_histo.root","recreate"); // Canvas Definition for Low Gain // TCanvas* cLG_RI = new TCanvas("Det-LG RI","Det-LG RI",1800,1800);