From 2c7807429fc592f6ce8548a5709b581f99ce7468 Mon Sep 17 00:00:00 2001
From: Morfouace <pierre.morfouace@cea.fr>
Date: Fri, 2 Sep 2022 12:05:05 -0600
Subject: [PATCH] Updating Vendeta macro
---
Projects/Vendeta/Analysis.cxx | 4 +-
Projects/Vendeta/calibration/Vendeta_Time.cal | 286 +++++++++---------
Projects/Vendeta/macro/FillTOFHisto.C | 104 +++++++
Projects/Vendeta/macro/FitTofGammaPeak.C | 3 +-
Projects/Vendeta/macro/GenerateTOFHisto.C | 14 +-
.../macro/RawSpectra/FillRawFCSpectra.C | 66 ++++
.../Vendeta/macro/RawSpectra/FillRawPSD.C | 72 +++++
.../Vendeta/macro/RawSpectra/ShowRawPSD.C | 5 +-
8 files changed, 402 insertions(+), 152 deletions(-)
create mode 100644 Projects/Vendeta/macro/FillTOFHisto.C
create mode 100644 Projects/Vendeta/macro/RawSpectra/FillRawFCSpectra.C
create mode 100644 Projects/Vendeta/macro/RawSpectra/FillRawPSD.C
diff --git a/Projects/Vendeta/Analysis.cxx b/Projects/Vendeta/Analysis.cxx
index 5eb10ffa6..c1a183235 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 dbd6a3af8..5bd681459 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 000000000..171b775bc
--- /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 061b71adc..610acbe97 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 6a66a83d3..34eecb304 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 000000000..cb4fd2e44
--- /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 000000000..f79733f9d
--- /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 8dd2b9816..505d31af4 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);
--
GitLab