diff --git a/macro/beamtime/mcbm2022/mMuchPar.par b/macro/beamtime/mcbm2022/mMuchPar.par
index 81b34c590fa1405f3e68f323da272ce030d84e28..9059ec01e608d42704b5b86bf630bf40208021a0 100644
--- a/macro/beamtime/mcbm2022/mMuchPar.par
+++ b/macro/beamtime/mcbm2022/mMuchPar.par
@@ -1,5 +1,5 @@
##############################################################################
-# Class: CbmMcbm2018MuchPar
+# Class: CbmMcbmUnpackPar
# Context: For mMUCH testing
##############################################################################
[CbmMuchUnpackPar]
@@ -14,14 +14,14 @@ CrobActiveFlag: Int_t \
1 1 1 1 1 1
NrOfFebsInGemA: Int_t 27
nFebsIdsArrayA: Int_t \
- 7 1 -1 10 17 9 -1 14 -1 \ #Id for the connected FEB Position in DPB 1001 GEM Module A
- 8 5 -1 16 13 12 -1 15 -1 \ #Id for the connected FEB Position in DPB 1002 GEM Module A
- 11 8 -1 4 3 -1 -1 0 2 #Id for the connected FEB Position in DPB 1003 GEM Module A
+ 7 8 -1 10 17 9 14 11 -1 \ #Id for the connected FEB Position in DPB 1001 GEM Module A
+ -1 5 -1 16 13 12 -1 15 -1 \ #Id for the connected FEB Position in DPB 1002 GEM Module A
+ -1 -1 -1 4 3 -1 -1 -7 -1 #Id for the connected FEB Position in DPB 1003 GEM Module A (-7 is Pulser)
// GEM2 Number of Febs and location of Febs
NrOfFebsInGemB: Int_t 18
nFebsIdsArrayB: Int_t \
6 16 -1 14 5 8 9 4 15 \ #Id for the connected FEB Position in DPB 2002 GEM Module B
- 10 3 -1 17 -1 11 12 7 13 #Id for the connected FEB Position in DPB 2003 GEM Module B
+ -1 3 -1 17 10 11 12 7 13 #Id for the connected FEB Position in DPB 2003 GEM Module B
// RPC Number of Febs and location of Febs
NrOfFebsInRpc: Int_t 9
nFebsIdsArrayRpc: Int_t \
diff --git a/macro/beamtime/mcbm2022/mMuchParUpto26032022.par b/macro/beamtime/mcbm2022/mMuchParUpto26032022.par
new file mode 100644
index 0000000000000000000000000000000000000000..ff00da31d6990a008a82422b448111ea54a57d7a
--- /dev/null
+++ b/macro/beamtime/mcbm2022/mMuchParUpto26032022.par
@@ -0,0 +1,1049 @@
+##############################################################################
+# Class: CbmMcbm2018MuchPar
+# Context: For mMUCH testing
+##############################################################################
+[CbmMuchUnpackPar]
+//----------------------------------------------------------------------------
+// Total number of MUCH DPBs in system
+NrOfDpbs: Int_t 6
+// Array to hold the unique IDs (equipment ID) for all MUCH DPBs, Array size has to match NrOfDpbs
+DbpIdArray: Int_t \
+ 0x1001 0x1002 0x1003 0x1006 0x2002 0x2003
+// Array to hold the active flag for all CROBs, [ NbDpb * kuNbCrobPerDpb ] = [ NbDpb * 2 ]
+CrobActiveFlag: Int_t \
+ 1 1 1 1 1 1
+NrOfFebsInGemA: Int_t 27
+nFebsIdsArrayA: Int_t \
+ 7 1 -1 10 17 9 -1 14 -1 \ #Id for the connected FEB Position in DPB 1001 GEM Module A
+ 8 5 -1 16 13 12 -1 15 -1 \ #Id for the connected FEB Position in DPB 1002 GEM Module A
+ 11 8 -1 4 3 -1 -1 0 2 #Id for the connected FEB Position in DPB 1003 GEM Module A
+// GEM2 Number of Febs and location of Febs
+NrOfFebsInGemB: Int_t 18
+nFebsIdsArrayB: Int_t \
+ 6 16 -1 14 5 8 9 4 15 \ #Id for the connected FEB Position in DPB 2002 GEM Module B
+ 10 3 -1 17 -1 11 12 7 13 #Id for the connected FEB Position in DPB 2003 GEM Module B
+// RPC Number of Febs and location of Febs
+NrOfFebsInRpc: Int_t 9
+nFebsIdsArrayRpc: Int_t \
+ 3 -1 -1 2 0 -1 1 -1 -1 #Id for the connected FEB Position in GEM Module B
+ChannelsToPadXA: Int_t \
+ 8 6 8 6 8 6 8 5 7 5 7 5 7 5 7 5 \
+ 7 4 6 4 6 4 8 4 9 4 9 3 9 3 9 3 \
+ 9 3 10 3 10 2 10 2 10 2 10 2 11 2 11 1 \
+ 11 1 11 1 11 1 12 1 12 0 12 0 12 0 12 0 \
+ 13 0 13 6 13 5 13 5 13 5 14 5 14 5 14 4 \
+ 14 4 14 4 15 4 15 4 15 3 15 3 15 3 16 3 \
+ 16 3 16 2 16 2 16 2 17 2 17 2 17 1 17 1 \
+ 17 1 18 1 18 1 18 0 18 0 18 0 19 0 19 0 \
+ \
+ -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 -2 -2 -2 17 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 19 -2 20 -2 20 -2 20 21 20 21 20 21 20 21 20 \
+ 21 20 21 21 22 21 22 21 22 -1 22 -1 22 -1 22 -1 \
+ 22 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 \
+ -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 \
+ \
+ 8 6 8 6 8 6 9 7 9 6 9 7 9 7 9 7 \
+ 10 7 10 7 10 7 10 7 10 8 11 8 11 0 11 0 \
+ 12 0 12 11 6 11 6 4 6 4 6 5 7 5 7 5 \
+ 7 4 8 4 8 4 8 4 8 4 8 4 8 3 8 3 \
+ 8 3 9 3 9 3 9 3 9 3 9 2 9 2 9 2 \
+ 9 2 10 2 10 2 10 2 10 2 10 1 10 1 10 1 \
+ 10 1 11 1 11 1 11 1 11 1 11 0 11 0 11 0 \
+ 6 0 6 0 6 0 6 0 5 0 5 0 5 3 5 5 \
+ \
+ 5 4 5 4 5 4 5 4 7 4 6 4 6 4 6 3 \
+ 6 3 6 3 6 3 6 3 6 3 6 3 6 3 7 3 \
+ 7 3 7 2 7 2 7 2 7 2 7 2 7 2 7 2 \
+ 8 2 8 2 8 2 8 1 8 1 8 1 8 1 8 1 \
+ 8 1 8 1 9 1 9 1 9 1 9 0 9 0 9 0 \
+ 4 0 4 0 4 0 4 4 4 4 4 4 4 5 4 5 \
+ 4 5 3 5 3 5 3 5 3 3 4 3 5 3 5 3 \
+ 5 3 5 3 5 2 5 2 5 2 5 2 5 2 5 2 \
+ \
+ 11 9 11 9 11 9 12 10 12 10 12 10 12 10 12 10 \
+ 12 10 12 10 12 10 13 10 13 10 13 11 13 11 13 11 \
+ 14 11 14 11 14 11 14 12 14 12 15 13 15 13 15 13 \
+ 15 13 15 13 16 14 16 14 16 14 16 14 17 14 17 15 \
+ 9 15 9 9 9 8 9 8 9 8 9 8 9 8 9 8 \
+ 9 9 -1 8 -1 8 10 8 10 8 10 7 10 7 10 7 \
+ 10 7 10 7 10 7 10 7 11 7 11 7 11 7 11 6 \
+ 11 6 11 6 11 6 11 6 11 6 11 6 12 6 12 6 \
+ \
+ 17 15 17 15 17 16 17 16 17 16 18 16 -1 16 18 17 \
+ 18 17 18 17 18 18 19 15 19 18 19 18 19 18 19 18 \
+ 19 19 20 19 20 19 20 19 20 20 20 20 21 20 21 20 \
+ 21 20 21 21 21 21 22 21 22 21 22 21 22 22 22 22 \
+ 22 22 22 22 22 22 14 22 14 22 14 13 14 13 14 13 \
+ 14 13 14 13 14 12 14 12 14 12 16 13 15 13 15 13 \
+ 15 13 15 13 15 12 15 12 15 12 15 12 16 12 16 15 \
+ 16 16 16 16 16 16 16 16 17 17 17 17 17 17 17 17 \
+ \
+ 19 17 19 17 19 18 19 18 19 18 19 18 19 18 20 18 \
+ 20 18 20 18 20 18 20 18 21 19 21 19 21 19 21 20 \
+ 21 20 22 20 22 20 22 20 22 21 22 21 22 21 16 21 \
+ 16 21 16 22 16 16 16 16 16 16 17 16 17 16 17 16 \
+ 18 17 18 17 18 17 18 17 18 17 18 17 19 17 19 17 \
+ 19 17 19 18 19 18 19 18 19 18 20 18 20 18 20 19 \
+ 20 19 20 19 20 19 21 19 21 20 21 20 21 20 21 20 \
+ 21 20 21 20 22 21 22 21 22 21 22 21 22 21 22 21 \
+ \
+ 12 11 12 11 12 11 12 11 13 11 13 11 13 11 13 12 \
+ 13 12 13 12 13 12 13 12 14 12 14 12 14 12 14 13 \
+ 14 13 14 13 14 13 15 14 15 14 15 14 15 14 15 14 \
+ 15 15 12 15 12 15 12 15 12 15 12 15 12 11 12 11 \
+ 13 12 13 12 13 12 14 13 14 13 14 13 14 13 14 13 \
+ 14 13 15 13 15 13 15 13 15 14 15 14 15 14 15 14 \
+ 16 14 16 14 16 15 16 15 16 15 16 15 17 15 17 16 \
+ 17 16 11 16 -1 16 -1 16 -1 16 -1 17 -1 17 -1 17 \
+ \
+ 8 7 8 7 8 7 9 7 -2 7 -2 7 -2 7 -2 8 \
+ -2 8 -2 8 9 8 -2 8 -2 8 -2 8 10 8 10 8 \
+ 10 9 10 10 10 10 10 10 11 10 11 10 11 10 8 11 \
+ 8 11 8 7 8 7 8 7 8 7 8 7 8 7 8 8 \
+ 8 7 8 7 9 7 9 7 9 7 9 7 9 6 9 6 \
+ 9 6 10 6 10 6 10 6 10 6 10 9 10 9 11 9 \
+ 11 9 11 9 11 10 11 10 11 10 12 10 -1 10 -1 10 \
+ -1 11 -1 11 -1 11 -1 11 8 12 -1 6 -1 7 -1 6 \
+ \
+ 4 2 4 2 -1 2 4 2 4 1 4 1 3 1 3 1 \
+ 3 1 3 1 3 1 3 1 4 1 4 1 4 0 4 0 \
+ 4 0 4 0 5 0 5 0 5 0 5 0 5 0 5 0 \
+ 5 3 5 3 5 3 5 3 5 3 5 3 6 2 6 2 \
+ 6 2 6 2 6 2 6 2 6 2 6 2 6 2 6 2 \
+ 6 2 6 2 6 1 7 1 7 1 7 1 7 1 7 1 \
+ -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 0 -1 0 \
+ -1 0 -1 0 -1 0 -1 0 -1 0 -1 0 -1 0 -1 0 \
+ \
+ 3 0 3 0 3 2 4 2 4 2 4 2 4 2 4 2 \
+ 4 3 4 3 4 3 4 3 4 3 4 3 4 3 5 3 \
+ 5 3 5 2 5 2 5 2 5 2 5 2 5 2 5 1 \
+ 5 1 5 1 5 1 6 1 6 -1 6 1 6 1 6 1 \
+ 3 1 3 1 3 1 3 0 3 0 3 0 3 0 3 0 \
+ 2 0 2 0 2 0 3 0 3 0 3 0 3 0 4 2 \
+ 4 2 4 2 4 2 4 2 4 2 4 2 4 2 4 2 \
+ 4 1 4 1 4 1 5 1 5 1 5 1 5 1 5 1 \
+ \
+ 8 6 8 6 8 6 8 6 8 6 9 6 9 6 9 6 \
+ 9 6 9 7 9 7 9 7 10 7 10 7 10 8 10 8 \
+ 10 8 10 9 11 9 11 9 11 9 11 9 7 10 7 10 \
+ 6 10 7 10 7 10 6 10 6 11 6 11 6 11 6 6 \
+ 7 6 7 6 7 5 7 5 7 5 7 5 7 5 8 5 \
+ 8 5 8 5 8 5 8 5 8 5 8 5 8 -1 8 -1 \
+ 8 4 8 4 8 4 9 4 9 4 9 6 7 6 7 6 \
+ 7 7 7 6 7 6 7 -2 7 5 6 5 6 5 6 5 \
+ \
+ 12 11 12 11 12 11 12 11 12 12 13 12 13 12 13 12 \
+ 13 12 10 12 10 12 10 13 11 13 11 13 11 13 11 11 \
+ 10 11 10 9 10 9 10 10 11 10 11 10 11 9 11 9 \
+ 11 9 12 10 12 9 12 9 -1 9 12 9 10 12 12 12 \
+ 13 12 13 12 13 12 13 11 13 12 14 13 14 13 14 13 \
+ 14 13 14 13 14 13 15 13 15 14 15 14 10 14 10 14 \
+ 10 14 11 14 11 15 11 9 11 9 10 9 11 9 -1 8 \
+ 11 8 11 8 10 8 10 8 -1 8 9 9 -1 8 10 9 \
+ \
+ 18 17 18 17 18 17 19 17 19 17 19 17 14 18 14 18 \
+ 15 18 15 18 15 18 15 18 15 18 15 18 15 18 15 19 \
+ 14 19 15 19 15 19 15 14 15 14 16 14 16 14 16 14 \
+ 16 14 16 14 16 14 16 14 16 13 16 13 16 13 16 13 \
+ 15 16 15 17 15 17 16 17 16 17 16 17 16 15 16 15 \
+ 16 15 16 15 16 15 13 16 16 16 16 12 17 12 17 12 \
+ 13 12 13 -1 13 -1 14 -1 14 -1 -1 -1 14 -1 14 -1 \
+ 13 -1 14 -1 14 -1 13 -1 13 -1 12 -1 12 -1 13 -1 \
+ \
+ -2 19 22 -2 -2 19 22 19 20 19 20 20 20 20 21 20 \
+ 21 20 21 20 21 21 21 21 22 21 22 21 22 21 22 21 \
+ 22 22 22 22 22 22 18 22 18 22 18 -2 18 17 18 17 \
+ 18 17 18 17 18 17 19 17 19 18 19 18 19 18 19 18 \
+ 20 19 20 19 20 19 21 19 21 19 21 19 21 19 21 20 \
+ 21 20 21 20 22 20 22 20 22 20 22 20 22 20 22 20 \
+ 22 21 22 21 17 21 17 21 17 21 17 22 17 22 17 22 \
+ 17 22 18 -1 18 17 -1 17 18 18 18 18 18 18 18 18 \
+ \
+ 22 22 -1 22 21 22 22 22 22 22 22 22 22 22 22 21 \
+ 21 21 21 21 21 21 21 21 21 21 21 20 20 20 20 20 \
+ 20 20 20 20 20 20 20 19 18 19 22 19 22 19 21 19 \
+ 22 19 17 19 21 19 22 19 21 19 22 19 22 19 22 18 \
+ 16 18 20 18 22 16 21 18 20 17 22 18 17 18 15 18 \
+ 21 18 20 15 19 18 19 17 18 17 18 17 17 16 17 17 \
+ 16 17 15 17 16 16 14 16 14 15 19 14 20 15 21 16 \
+ 22 16 21 16 20 16 19 15 19 15 20 15 20 15 13 15 \
+ \
+ 12 12 -1 19 -1 20 -1 19 -1 21 -1 20 -1 21 -1 19 \
+ -1 11 -1 21 -1 20 -1 20 -1 19 -1 19 -1 -1 -1 -1 \
+ -1 -1 -1 -1 11 -1 13 -1 15 -1 14 -1 10 10 12 -1 \
+ 13 -1 11 -1 12 -1 11 -1 10 -1 10 -1 9 -1 9 9 \
+ -1 -1 -1 -1 -1 -1 -1 8 -1 -1 8 -1 -1 -1 -1 -1 \
+ -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 \
+ -1 7 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 8 -1 8 7 \
+ 7 9 6 7 4 6 5 5 6 6 5 6 5 5 4 5 \
+ \
+ 5 1 -2 1 -2 -2 -2 0 -2 0 -2 0 -2 -1 -2 0 \
+ 4 0 3 -2 3 0 2 0 2 0 2 0 2 0 2 -1 \
+ 2 1 2 1 2 2 1 1 2 1 2 1 -1 0 -2 0 \
+ -2 0 -2 0 -2 0 -2 0 -2 0 -2 0 -2 -1 -2 -1 \
+ -2 -1 3 -1 3 -1 3 -1 3 -1 3 -1 3 -1 3 -1 \
+ 3 -1 4 -1 4 -1 4 -1 4 -1 4 2 4 2 4 2 \
+ 3 2 4 2 2 2 1 1 1 1 1 1 1 1 0 1 \
+ 0 0 0 0 2 0 0 0 3 1 4 0 3 0 2 1
+ChannelsToPadYA: Int_t \
+ 1 2 2 3 3 4 4 0 0 1 1 2 2 3 3 4 \
+ 4 0 0 1 1 2 0 3 4 4 3 0 2 1 1 2 \
+ 0 3 4 4 3 0 2 1 1 2 0 3 4 4 3 0 \
+ 2 1 1 2 0 3 4 4 3 0 2 1 1 2 0 3 \
+ 4 4 3 7 2 5 1 6 0 9 4 8 3 7 2 5 \
+ 1 6 0 9 4 8 3 7 2 5 1 6 0 9 4 8 \
+ 3 7 2 5 1 6 0 9 4 8 3 7 2 5 1 6 \
+ 0 9 4 8 3 7 2 5 1 6 0 9 4 8 3 7 \
+ \
+ -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 -2 -2 -2 11 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 -2 \
+ -2 17 -2 10 -2 11 -2 12 17 13 9 14 8 15 7 16 \
+ 6 17 5 10 15 11 16 12 17 -1 14 -1 13 -1 12 -1 \
+ 11 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 \
+ -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 \
+ \
+ 7 8 6 9 5 5 9 7 8 6 7 8 6 9 5 14 \
+ 9 13 8 12 7 11 6 10 5 9 9 8 8 12 6 11 \
+ 9 10 8 7 13 5 12 11 11 10 10 15 15 16 16 17 \
+ 17 12 10 13 11 14 12 17 13 16 14 15 15 10 16 11 \
+ 17 13 10 14 11 17 12 16 13 15 14 10 15 11 16 12 \
+ 17 13 10 14 11 17 12 16 13 15 14 10 15 11 16 12 \
+ 17 13 10 14 11 17 12 16 13 15 14 17 15 18 16 19 \
+ 14 20 17 21 16 16 15 15 10 14 11 13 12 12 13 14 \
+ \
+ 27 24 26 23 25 22 24 21 21 20 18 19 19 18 20 27 \
+ 21 26 22 25 23 24 24 23 25 22 26 21 27 20 18 19 \
+ 19 18 20 27 22 26 23 25 24 24 25 23 26 22 27 21 \
+ 18 20 19 19 20 18 21 27 22 26 23 25 24 24 25 23 \
+ 26 22 27 21 18 20 19 19 20 18 21 27 22 26 23 25 \
+ 34 24 35 23 36 22 37 25 32 26 31 27 30 18 29 19 \
+ 28 20 28 21 29 22 30 23 31 32 33 37 33 36 34 35 \
+ 35 34 36 33 37 33 32 34 31 35 30 36 29 37 28 32 \
+ \
+ 25 27 26 26 27 25 18 18 19 19 20 20 21 21 22 28 \
+ 24 23 26 24 27 25 19 26 21 27 23 18 25 19 27 20 \
+ 19 21 21 22 22 23 24 23 26 25 18 18 20 20 22 22 \
+ 24 24 26 26 18 18 20 20 22 23 24 25 18 27 24 19 \
+ 30 21 31 24 32 36 37 35 36 34 35 33 34 32 29 37 \
+ 28 33 -1 28 -1 29 29 30 30 31 31 28 37 29 36 30 \
+ 35 31 34 32 33 36 32 35 34 34 35 37 36 33 37 30 \
+ 33 31 32 32 31 35 30 36 29 33 28 34 28 29 29 28 \
+ \
+ 19 25 21 27 23 19 25 21 27 23 18 25 -1 26 21 20 \
+ 23 22 24 26 26 19 18 23 19 20 21 22 23 25 25 27 \
+ 27 20 19 22 21 24 22 26 24 18 26 20 18 23 20 25 \
+ 22 27 24 19 26 21 18 23 20 25 21 27 23 19 24 22 \
+ 32 25 30 26 28 27 28 31 29 29 30 29 37 30 36 31 \
+ 35 28 34 37 33 32 32 33 31 34 27 36 28 35 37 34 \
+ 36 33 35 32 34 35 33 36 32 37 31 31 28 30 37 30 \
+ 35 29 33 36 31 34 30 32 29 28 37 36 35 34 33 32 \
+ \
+ 36 31 35 30 34 28 33 29 32 37 31 36 30 35 29 34 \
+ 36 33 34 32 32 31 30 30 29 28 36 29 35 37 33 28 \
+ 31 37 33 35 39 33 38 31 37 28 36 37 35 34 49 32 \
+ 48 30 47 34 46 38 45 39 44 40 38 41 39 42 40 43 \
+ 39 41 40 42 42 49 48 48 46 47 44 46 38 45 40 44 \
+ 42 43 48 38 46 41 45 49 43 47 39 45 41 43 49 39 \
+ 47 41 45 49 43 47 39 44 41 38 42 40 49 42 47 48 \
+ 45 46 43 44 40 38 41 40 42 50 43 48 49 46 48 44 \
+ \
+ 44 44 45 45 46 46 43 47 45 48 46 49 47 43 48 38 \
+ 49 39 43 40 41 41 39 42 45 49 47 48 49 47 43 44 \
+ 41 42 39 40 38 38 39 46 41 48 49 44 47 42 45 40 \
+ 43 38 55 40 54 42 53 48 52 46 61 44 60 55 59 54 \
+ 61 58 59 57 57 56 50 50 52 51 54 52 61 53 59 54 \
+ 57 55 50 60 52 58 54 56 61 51 59 53 57 55 56 60 \
+ 51 58 53 56 55 51 60 53 58 55 56 60 51 58 53 50 \
+ 55 52 67 54 -1 61 -1 59 -1 57 -1 50 -1 52 -1 54 \
+ \
+ 45 43 44 44 43 49 38 42 -2 41 -2 40 -2 39 -2 38 \
+ -2 39 -2 40 47 41 -2 42 -2 49 -2 48 39 47 41 46 \
+ 49 44 47 38 45 40 43 42 39 48 40 46 42 44 61 38 \
+ 60 41 59 61 55 60 54 59 53 58 52 57 51 56 58 50 \
+ 57 52 56 53 50 54 51 55 52 51 54 50 61 56 59 57 \
+ 57 58 50 59 52 60 54 61 61 55 59 53 57 55 50 60 \
+ 52 58 61 56 60 51 59 53 57 55 50 60 -1 58 -1 56 \
+ -1 51 -1 53 -1 58 -1 56 67 51 -1 62 -1 62 -1 63 \
+ \
+ 43 31 42 30 -1 29 40 28 39 28 38 29 44 30 45 31 \
+ 46 32 47 37 48 36 49 35 49 34 48 33 47 33 46 34 \
+ 45 35 44 36 38 37 39 32 40 31 41 30 42 29 43 28 \
+ 49 40 48 41 47 42 46 43 45 39 44 38 38 44 47 45 \
+ 39 46 40 47 41 48 42 49 49 43 48 42 46 41 45 40 \
+ 44 39 43 38 37 44 38 45 45 46 46 47 47 48 48 49 \
+ -1 43 -1 42 -1 41 -1 40 -1 39 -1 38 -1 44 -1 45 \
+ -1 46 -1 47 -1 48 -1 49 -1 43 -1 42 -1 41 -1 40 \
+ \
+ 58 39 57 38 56 56 50 57 51 58 52 59 53 60 54 61 \
+ 55 50 61 51 60 52 59 53 58 54 57 55 56 61 50 60 \
+ 51 59 52 55 53 54 54 53 55 52 61 51 60 50 59 56 \
+ 58 57 57 58 56 59 50 60 51 -1 52 55 53 54 54 53 \
+ 70 52 71 51 72 50 73 50 65 51 64 52 63 53 62 54 \
+ 68 55 69 61 70 60 69 59 68 58 67 57 66 56 62 71 \
+ 63 72 64 73 69 67 68 66 67 65 73 64 72 63 71 62 \
+ 70 68 66 69 65 70 68 71 69 72 70 73 71 67 72 66 \
+ \
+ 72 65 73 64 70 66 65 67 63 73 67 72 70 71 72 70 \
+ 66 69 65 63 64 64 62 65 69 66 71 67 73 66 66 64 \
+ 64 62 62 69 69 71 71 73 73 68 66 63 75 68 74 70 \
+ 74 72 77 67 76 65 75 63 76 68 77 70 78 72 82 81 \
+ 84 80 83 79 82 74 81 75 80 76 79 77 78 78 82 79 \
+ 83 80 84 81 85 85 81 84 80 83 79 82 78 -1 77 -1 \
+ 76 76 75 77 74 78 74 79 75 80 76 85 86 84 87 83 \
+ 88 85 89 87 90 88 91 -2 92 87 91 88 90 89 89 90 \
+ \
+ 64 65 66 64 73 63 71 62 69 62 65 63 67 65 70 67 \
+ 69 72 79 70 78 68 77 66 76 68 77 73 78 72 79 75 \
+ 80 74 81 79 82 78 84 74 84 83 83 75 82 80 81 81 \
+ 80 82 74 76 75 83 76 84 -1 85 77 77 85 79 78 84 \
+ 83 83 85 82 78 81 76 85 74 80 75 80 77 81 78 82 \
+ 84 84 81 79 79 77 75 75 76 74 77 76 88 82 87 85 \
+ 86 83 86 80 87 74 88 88 89 87 89 86 90 89 -1 86 \
+ 91 87 92 88 92 89 90 90 -1 91 92 90 -1 92 91 91 \
+ \
+ 60 61 58 60 56 59 51 58 53 57 55 56 63 50 62 51 \
+ 68 52 69 53 70 54 71 55 72 61 73 59 67 57 66 50 \
+ 67 52 65 54 64 61 63 64 62 65 68 66 69 73 70 72 \
+ 71 71 72 70 73 69 67 68 66 71 65 62 64 63 63 64 \
+ 82 62 81 68 80 69 80 70 81 71 82 72 83 78 84 79 \
+ 85 85 79 84 77 83 89 78 75 76 74 86 80 87 81 88 \
+ 88 89 87 -1 86 -1 86 -1 87 -1 -1 -1 89 -1 90 -1 \
+ 90 -1 91 -1 92 -1 91 -1 92 -1 91 -1 92 -1 93 -1 \
+ \
+ -2 60 46 -2 -2 58 44 57 58 56 57 50 56 51 52 52 \
+ 54 53 61 54 59 51 57 53 56 55 58 60 60 58 61 56 \
+ 54 57 52 59 50 55 73 53 72 51 67 -2 66 67 65 66 \
+ 64 65 63 64 62 63 68 62 69 68 70 69 71 70 72 71 \
+ 66 73 64 67 62 66 69 65 71 64 73 63 66 62 64 68 \
+ 63 69 62 70 63 71 65 72 67 73 72 67 71 65 70 63 \
+ 69 68 68 70 82 72 83 67 84 65 85 62 78 64 76 66 \
+ 74 73 81 -1 83 77 -1 75 78 80 77 82 76 84 75 79 \
+ \
+ 80 81 -1 82 94 83 94 84 85 79 78 77 76 75 74 74 \
+ 75 76 77 78 79 85 84 83 82 81 80 74 75 76 77 78 \
+ 79 85 84 83 82 81 80 74 94 75 92 76 91 77 91 78 \
+ 90 79 94 85 90 84 89 83 89 82 88 81 87 80 86 74 \
+ 94 93 94 92 95 93 95 91 95 92 96 90 93 89 94 88 \
+ 96 87 96 93 95 85 96 91 95 90 96 89 95 92 96 88 \
+ 95 87 95 86 96 90 95 91 94 92 93 93 93 91 93 89 \
+ 93 88 92 87 92 86 92 90 91 89 91 88 90 87 94 86 \
+ \
+ 94 93 -1 90 -1 89 -1 89 -1 88 -1 88 -1 87 -1 88 \
+ -1 93 -1 86 -1 87 -1 86 -1 87 -1 86 -1 -1 -1 -1 \
+ -1 -1 -1 -1 94 -1 95 -1 96 -1 96 -1 94 93 95 -1 \
+ 96 -1 95 -1 96 -1 96 -1 95 -1 96 -1 95 -1 94 93 \
+ -1 -1 -1 -1 -1 -1 -1 93 -1 -1 94 -1 -1 -1 -1 -1 \
+ -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 \
+ -1 93 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 95 -1 96 94 \
+ 96 96 95 95 96 96 95 96 94 92 94 93 93 91 94 92 \
+ \
+ 73 64 -2 61 -2 -2 -2 68 -2 69 -2 70 -2 -1 -2 72 \
+ 75 71 75 -2 74 66 76 65 77 64 78 63 79 62 80 -1 \
+ 81 78 85 80 84 83 85 83 82 82 75 75 -1 76 -2 78 \
+ -2 80 -2 85 -2 83 -2 82 -2 75 -2 74 -2 -1 -2 -1 \
+ -2 -1 92 -1 93 -1 91 -1 90 -1 89 -1 88 -1 87 -1 \
+ 86 -1 86 -1 87 -1 88 -1 89 -1 90 91 91 90 92 89 \
+ 94 88 93 87 94 86 91 92 89 93 87 90 86 88 87 94 \
+ 89 86 91 88 95 90 94 93 95 95 95 95 96 96 96 96
+ChannelsToPadXB: Int_t \
+ 8 6 8 6 8 6 9 6 9 6 9 6 9 7 9 7 \
+ 9 7 10 7 10 7 10 7 10 8 11 8 11 8 11 5 \
+ 12 10 12 10 12 11 13 11 13 11 13 12 14 12 14 12 \
+ 14 13 14 13 15 13 15 14 15 14 15 5 15 5 16 5 \
+ 16 5 16 5 16 4 16 4 17 4 17 4 17 4 17 4 \
+ 17 3 18 3 18 3 18 3 18 3 18 3 19 2 19 2 \
+ 19 2 19 2 19 2 20 2 20 1 20 1 20 1 20 1 \
+ 21 1 21 1 21 0 21 0 21 0 22 0 22 0 21 0 \
+ \
+17 14 17 15 18 15 18 15 18 15 18 15 18 15 18 16 \
+ 19 16 19 16 19 16 19 16 19 16 19 17 20 17 20 17 \
+ 20 17 20 14 20 14 20 14 21 14 21 14 21 13 21 13 \
+ 21 13 21 13 22 13 22 13 22 12 22 12 22 12 22 12 \
+ 22 12 22 12 22 12 13 12 13 12 13 11 13 11 13 11 \
+ 14 11 14 11 14 10 14 14 15 15 15 15 15 16 16 16 \
+ 16 16 17 17 17 17 17 18 18 18 18 18 19 19 19 19 \
+ 19 20 20 20 20 20 21 21 21 21 21 22 22 22 22 22 \
+ \
+ 7 0 7 0 7 6 8 6 8 6 8 6 8 5 8 5 \
+ 8 5 9 5 9 5 9 5 9 4 9 4 9 4 10 4 \
+ 10 4 10 3 10 3 10 3 10 3 11 3 11 2 11 2 \
+ 11 2 11 2 11 2 12 1 12 1 7 1 7 1 7 1 \
+ 6 6 6 6 6 0 6 0 7 0 7 4 7 4 7 3 \
+ 7 3 6 3 8 3 8 3 8 2 8 2 8 2 9 2 \
+ 9 2 9 2 9 1 9 1 10 1 10 1 10 1 10 1 \
+ 5 0 5 0 5 0 5 0 5 0 4 0 4 0 4 0 \
+ \
+ 5 4 5 4 5 4 5 5 6 5 6 4 6 4 6 4 \
+ 6 4 6 3 7 3 7 3 7 3 7 3 7 3 7 3 \
+ 8 2 8 2 8 2 8 2 8 2 8 2 9 1 9 1 \
+ 9 1 9 1 9 1 9 1 10 0 10 0 10 0 10 0 \
+ 4 2 4 2 4 2 4 2 3 2 3 2 3 2 3 2 \
+ 3 2 3 1 3 1 3 1 3 1 4 1 4 1 4 1 \
+ 4 1 4 1 5 0 5 0 5 0 5 0 5 0 5 0 \
+ 5 0 5 0 5 0 6 0 6 0 6 0 6 0 6 0 \
+ \
+ 11 8 10 8 10 8 11 8 11 8 11 7 11 7 11 8 \
+ 12 8 12 7 12 8 12 8 12 7 12 7 13 7 13 7 \
+ 13 7 13 7 13 6 13 6 14 6 14 6 14 9 14 9 \
+ 14 9 14 9 15 9 15 9 15 9 15 9 15 9 15 10 \
+ 12 11 12 11 12 11 12 11 13 11 13 11 13 12 13 12 \
+ 13 12 13 12 17 12 10 13 10 13 10 13 10 8 10 8 \
+ 10 8 10 8 10 8 10 8 10 8 10 11 11 11 11 11 \
+ 11 10 11 10 11 10 11 10 11 10 11 10 11 10 11 10 \
+ \
+ 16 14 15 14 16 14 16 14 16 14 16 14 16 14 16 14 \
+ 17 14 17 15 17 15 17 15 17 15 17 15 17 15 18 15 \
+ 18 15 18 15 18 16 18 16 18 16 18 16 19 16 19 16 \
+ 19 16 19 16 19 16 19 17 19 17 20 17 20 17 20 17 \
+ 20 17 20 17 20 17 20 18 21 18 21 18 21 18 21 18 \
+ 21 18 21 18 22 18 22 18 22 19 22 19 22 19 22 19 \
+ 22 19 22 19 22 19 22 19 21 19 21 20 21 20 21 20 \
+ 21 20 21 20 21 20 22 20 22 20 22 20 22 21 22 21 \
+ \
+ 20 19 20 19 20 20 20 20 20 20 21 20 21 16 21 16 \
+ 21 16 21 16 21 16 21 16 21 16 21 16 21 16 21 16 \
+ 22 16 22 16 22 17 22 17 22 17 22 17 22 17 22 17 \
+ 22 17 22 17 22 17 22 17 22 17 22 17 22 18 22 18 \
+ 19 19 19 19 19 19 20 19 20 19 20 19 20 19 20 20 \
+ 20 20 21 20 21 20 21 20 21 20 21 21 21 21 21 21 \
+ 22 21 22 21 22 22 22 22 22 14 22 14 22 14 22 14 \
+ 22 14 22 14 22 14 22 14 22 14 22 14 22 13 22 14 \
+ \
+16 15 16 15 16 15 16 15 17 15 17 15 17 15 17 15 \
+ 17 15 17 16 17 16 17 16 18 16 18 16 18 16 18 16 \
+ 18 17 18 17 19 17 19 18 19 18 19 18 19 18 20 18 \
+ 12 19 12 19 13 19 13 19 13 20 13 11 13 11 13 11 \
+ 13 11 14 11 14 11 14 12 14 11 14 13 14 12 14 12 \
+ 14 12 14 12 14 12 14 12 15 12 15 12 15 12 15 13 \
+ 15 13 15 13 15 14 15 15 18 15 18 15 18 15 18 10 \
+ 18 10 18 10 18 10 18 10 18 10 18 10 19 10 19 10 \
+ \
+ 12 8 12 8 12 8 12 7 12 7 12 7 12 7 12 7 \
+ 12 7 12 7 13 7 13 7 13 6 13 6 13 8 13 6 \
+ 13 6 13 9 13 9 13 9 13 9 14 9 14 9 14 9 \
+ 14 9 14 9 14 9 14 9 14 11 14 12 15 14 15 14 \
+ 8 7 8 7 9 7 9 7 9 8 9 8 9 8 9 8 \
+ 9 8 9 7 9 7 9 8 9 8 9 8 10 8 10 8 \
+ 10 7 10 7 10 7 10 7 10 7 10 7 10 6 10 6 \
+ 10 6 10 6 11 6 13 6 11 6 11 6 11 6 11 6 \
+ \
+ 4 3 4 3 5 3 4 3 4 3 4 3 4 3 4 3 \
+ 4 3 4 3 4 2 4 2 3 2 5 2 5 2 5 2 \
+ 5 2 5 2 5 2 5 2 5 2 5 1 5 1 6 1 \
+ 7 1 6 1 7 1 6 1 6 1 6 1 6 1 6 1 \
+ 0 0 0 0 0 0 4 3 4 3 4 3 4 3 4 3 \
+ 4 3 4 3 4 3 3 3 3 2 3 2 4 2 4 2 \
+ 4 2 4 2 5 2 5 2 5 2 5 2 5 2 5 2 \
+ 5 1 5 1 5 1 5 1 5 1 5 1 6 1 6 1 \
+ \
+ 1 0 1 0 1 0 1 0 0 0 3 0 3 2 3 2 \
+ 3 2 3 2 2 2 2 2 2 1 2 1 2 1 2 1 \
+ 2 1 3 1 3 1 3 1 3 1 3 1 3 1 3 1 \
+ 3 1 4 0 4 0 4 0 4 0 4 0 4 0 4 0 \
+ 0 0 0 0 0 2 0 2 0 2 0 2 0 2 0 1 \
+ 0 1 0 1 0 1 0 1 3 1 2 1 2 1 2 1 \
+ 3 1 3 1 3 0 3 0 3 0 3 0 3 0 3 0 \
+ 3 0 3 0 4 0 4 0 4 0 4 2 4 2 4 2 \
+ \
+ 7 6 7 6 7 6 7 6 7 6 7 6 7 6 7 6 \
+ 7 6 7 6 7 6 7 6 7 6 8 5 8 5 8 5 \
+ 8 5 8 5 8 5 8 5 8 5 8 5 8 5 8 5 \
+ 8 5 8 5 9 4 9 4 9 4 9 4 9 4 9 4 \
+ 6 5 6 5 6 5 6 5 6 5 7 5 7 5 7 5 \
+ 7 5 7 5 7 5 7 4 7 4 7 4 7 4 7 4 \
+ 7 6 7 6 7 6 7 6 7 6 7 6 7 6 7 6 \
+ 7 6 7 6 8 6 8 6 8 6 8 6 8 6 7 6 \
+ \
+ 11 9 11 9 11 9 11 9 11 9 12 9 12 9 12 10 \
+ 12 10 12 10 12 10 12 11 12 11 12 11 12 11 12 11 \
+ 12 11 12 11 13 11 13 9 13 9 13 9 13 9 13 9 \
+ 13 9 13 9 13 9 13 9 14 9 13 9 13 10 14 10 \
+ 11 8 11 8 11 8 11 8 11 8 11 8 11 8 11 8 \
+ 11 8 11 8 11 8 12 10 12 10 12 10 12 10 12 10 \
+ 12 10 12 10 12 10 10 10 10 9 10 9 10 9 10 9 \
+ 10 9 10 9 10 9 10 9 10 9 11 9 11 10 11 9 \
+ \
+ 16 15 16 15 16 15 16 15 16 15 16 15 16 15 16 15 \
+ 16 15 17 15 17 15 17 15 17 15 17 16 17 16 18 16 \
+ 18 16 18 17 18 17 18 17 18 17 18 17 19 17 19 17 \
+ 19 18 19 18 14 18 14 18 14 18 14 18 14 19 14 19 \
+ 14 12 14 12 14 12 14 13 14 13 14 13 15 13 14 13 \
+ 13 13 13 13 13 13 13 13 14 13 14 13 14 13 14 13 \
+ 14 12 14 12 14 12 14 12 16 12 15 12 15 12 14 12 \
+ 15 12 15 12 15 12 15 13 15 13 15 13 15 13 16 13 \
+ \
+ 20 19 20 19 20 19 20 19 20 19 20 19 20 19 20 20 \
+ 20 20 20 20 21 16 21 16 21 16 21 16 21 16 21 16 \
+ 21 16 21 16 21 16 21 16 21 16 21 17 21 17 22 17 \
+ 22 17 22 17 22 17 22 17 22 17 22 17 22 17 22 17 \
+ 18 15 18 15 18 15 18 15 18 15 19 15 19 15 19 15 \
+ 19 15 19 15 19 15 19 18 20 18 20 18 20 18 20 18 \
+ 20 18 20 19 21 19 21 19 21 19 21 20 21 20 21 20 \
+ 21 20 22 20 22 21 22 21 22 21 22 21 22 17 22 17 \
+ \
+ 22 22 21 21 21 22 22 22 21 22 22 22 22 22 22 22 \
+ 22 22 22 22 21 21 21 21 21 21 21 21 21 21 21 20 \
+ 18 20 19 20 20 20 20 20 21 20 21 20 21 20 20 20 \
+ 20 20 20 20 17 19 19 19 18 19 19 19 19 19 19 19 \
+ 16 19 -1 19 20 19 -1 19 -1 19 -1 18 19 18 18 18 \
+ -1 18 17 15 15 18 16 18 -1 18 20 18 19 18 18 18 \
+ 17 18 16 17 -1 17 15 17 15 17 -1 17 14 14 14 17 \
+ 14 17 13 17 13 16 17 16 16 16 16 16 16 16 15 16 \
+ \
+ 12 18 12 17 13 18 12 17 -1 12 -1 16 -1 18 -1 17 \
+ -1 -1 11 11 -1 -1 -1 -1 -1 -1 -1 -1 -1 15 -1 -1 \
+ 11 14 -1 16 -1 15 -1 16 -1 13 -1 12 -1 10 -1 -1 \
+ -1 11 -1 11 -1 11 -1 11 -1 11 16 11 17 11 13 11 \
+ -1 9 -1 11 -1 -1 -1 -1 -1 8 -1 -1 -1 -1 -1 -1 \
+ -1 -1 -1 -1 -1 -1 -1 -1 -1 14 -1 13 -1 15 -1 14 \
+ -1 13 -1 7 -1 12 -1 10 -1 9 10 11 11 12 9 10 \
+ 8 11 7 8 -1 9 10 10 10 9 9 8 9 7 6 6 \
+ \
+ 3 0 4 0 4 0 4 0 4 0 4 0 4 0 4 0 \
+ 4 3 4 3 4 3 4 3 5 3 5 3 5 3 5 3 \
+ 5 3 5 3 5 2 5 2 5 2 5 2 5 2 6 2 \
+ 8 2 8 2 8 2 8 2 8 2 8 1 1 1 1 1 \
+ 5 1 5 1 4 1 5 1 3 1 6 0 6 0 4 0 \
+ 7 1 8 0 7 0 6 1 5 0 3 1 4 0 5 1 \
+ 8 2 8 3 7 3 7 3 6 4 6 4 5 2 5 3 \
+ 4 2 3 2 4 2 -1 2 -1 2 -1 1 -1 1 1 1
+ChannelsToPadYB: Int_t \
+ 0 2 1 3 2 4 5 5 4 1 3 0 2 5 1 4 \
+ 0 3 5 2 4 1 2 0 0 5 4 4 2 3 0 0 \
+ 4 3 2 1 0 5 4 3 3 1 1 5 5 3 3 1 \
+ 2 5 0 2 4 0 3 4 2 1 1 1 0 2 4 3 \
+ 3 4 2 5 1 0 0 1 4 2 3 3 2 4 1 5 \
+ 0 0 4 1 3 2 2 3 1 4 0 5 4 0 3 1 \
+ 2 2 1 3 0 4 4 5 3 0 2 1 1 2 0 3 \
+ 5 4 4 5 3 0 2 1 1 3 1 5 0 7 0 6 \
+ \
+ 9 11 10 5 5 6 6 7 7 8 8 9 9 10 10 5 \
+ 5 6 6 7 7 8 8 9 9 10 10 5 5 6 6 7 \
+ 7 8 8 10 9 9 10 8 6 7 7 6 8 11 9 10 \
+ 10 9 11 8 2 7 3 6 4 11 5 10 6 9 7 8 \
+ 10 12 9 13 8 14 16 15 15 16 14 12 13 13 12 14 \
+ 16 15 15 16 14 12 12 13 14 15 12 13 11 15 14 13 \
+ 12 11 15 14 13 12 11 15 14 13 12 11 15 14 13 12 \
+ 11 15 14 13 12 11 16 15 14 13 12 14 15 13 12 11 \
+ \
+ 9 4 10 2 11 6 6 7 7 8 8 9 9 11 10 10 \
+ 11 9 6 8 7 7 8 6 9 10 10 9 11 8 6 7 \
+ 7 6 8 10 9 9 10 8 11 7 6 6 7 10 8 9 \
+ 9 8 10 7 11 6 6 10 7 9 8 8 7 7 6 6 \
+ 15 10 14 11 13 10 12 9 16 8 15 14 14 15 13 11 \
+ 12 12 16 13 16 14 15 15 14 11 13 12 12 13 16 14 \
+ 15 15 14 16 13 11 12 12 16 13 15 14 14 15 13 16 \
+ 12 14 13 15 14 16 15 17 16 18 11 13 12 12 13 11 \
+ \
+ 22 20 21 21 20 22 19 17 17 18 18 19 19 18 20 17 \
+ 21 16 22 22 17 21 18 20 19 19 20 18 21 17 22 16 \
+ 17 22 18 21 19 20 20 19 21 18 22 17 17 22 18 21 \
+ 19 20 20 19 21 18 22 17 17 22 18 21 19 20 20 19 \
+ 28 23 29 24 30 25 31 26 23 27 24 28 25 29 26 30 \
+ 27 31 28 23 29 24 30 25 31 26 27 27 26 28 25 29 \
+ 24 30 23 31 31 23 30 24 29 25 28 26 27 27 26 28 \
+ 25 29 24 30 23 31 26 36 31 35 30 34 29 33 28 32 \
+ \
+ 17 23 22 24 21 25 18 26 19 27 20 23 21 24 22 28 \
+ 17 29 18 25 19 30 20 31 21 26 22 28 17 29 18 27 \
+ 19 30 20 31 21 23 22 24 17 25 18 27 19 31 20 30 \
+ 21 29 22 28 16 27 17 26 18 25 19 24 20 23 21 31 \
+ 29 28 27 27 25 26 23 25 30 24 28 23 26 31 25 30 \
+ 24 28 23 26 26 24 34 31 35 29 33 27 32 42 36 41 \
+ 37 40 38 39 39 38 40 37 41 36 42 29 32 30 33 31 \
+ 34 23 35 24 36 25 37 26 38 27 39 28 40 29 41 30 \
+ \
+ 16 28 22 27 17 26 18 29 19 30 20 31 21 25 22 24 \
+ 16 23 17 31 18 30 19 29 20 28 21 27 22 26 16 25 \
+ 17 24 18 23 19 31 20 30 21 29 22 28 16 27 17 26 \
+ 18 25 19 24 20 23 21 31 22 30 16 29 17 28 18 27 \
+ 19 25 20 24 21 23 22 31 17 30 18 29 19 28 20 27 \
+ 21 26 22 25 16 24 17 23 18 31 19 30 20 29 21 28 \
+ 25 27 24 26 23 25 22 24 29 23 28 31 27 30 26 29 \
+ 25 28 24 27 23 26 30 25 29 24 28 23 27 31 26 30 \
+ \
+ 38 40 39 42 40 33 41 35 42 36 32 37 33 49 34 48 \
+ 35 47 36 46 37 45 38 44 39 43 40 50 41 51 42 52 \
+ 32 53 33 54 34 43 35 44 36 45 37 46 38 47 39 48 \
+ 40 49 41 50 42 51 46 52 45 53 44 54 43 43 31 44 \
+ 51 44 52 45 53 46 43 48 45 49 47 50 49 54 51 44 \
+ 53 46 43 48 45 50 46 52 47 54 49 44 51 48 53 50 \
+ 47 52 49 54 51 48 52 50 53 65 54 63 55 67 56 62 \
+ 57 61 58 60 64 59 63 58 62 57 61 56 60 55 59 55 \
+ \
+ 39 34 40 35 41 36 42 37 32 38 33 39 34 40 35 41 \
+ 36 42 38 32 39 33 41 34 32 35 34 36 36 37 38 38 \
+ 40 37 42 40 32 42 34 33 36 35 38 37 39 39 32 41 \
+ 51 33 52 35 45 37 46 41 47 34 48 54 49 53 50 51 \
+ 54 50 43 49 44 48 45 43 46 52 47 44 48 44 49 45 \
+ 50 46 51 47 52 48 54 49 44 50 46 53 48 54 49 51 \
+ 50 52 52 53 53 53 54 43 45 45 46 47 47 51 48 65 \
+ 49 66 50 67 51 64 52 63 53 62 54 61 43 60 47 59 \
+ \
+ 33 35 34 34 35 33 36 40 37 42 38 41 39 39 40 34 \
+ 41 35 42 36 32 37 33 38 34 37 35 38 36 32 37 42 \
+ 38 41 39 32 40 33 41 34 42 35 32 36 33 37 34 38 \
+ 35 39 36 40 37 41 39 42 41 42 42 32 32 38 33 40 \
+ 53 52 54 51 43 53 44 54 45 43 46 44 47 45 48 46 \
+ 49 47 50 50 51 49 52 48 53 49 54 50 43 51 44 52 \
+ 45 48 46 47 47 46 48 45 49 44 50 43 51 54 52 53 \
+ 53 52 54 51 43 50 43 49 44 48 45 43 46 44 47 46 \
+ \
+ 42 41 41 40 32 39 40 38 39 37 38 36 37 35 36 34 \
+ 35 33 34 32 33 42 32 41 42 40 33 39 34 38 35 37 \
+ 36 36 37 35 38 34 39 33 40 32 41 42 42 41 32 40 \
+ 32 39 33 38 33 37 34 36 35 35 36 34 39 33 40 32 \
+ 42 39 41 38 40 37 50 51 49 50 48 49 47 48 46 47 \
+ 45 46 44 45 43 44 54 43 53 54 52 53 51 52 52 51 \
+ 53 50 54 49 43 48 44 47 45 46 46 45 47 44 48 43 \
+ 49 54 50 53 51 52 52 51 53 50 54 49 45 48 47 47 \
+ \
+ 46 48 45 47 44 46 43 45 49 44 63 43 64 62 65 63 \
+ 66 64 67 65 55 66 56 67 57 55 58 56 59 57 60 58 \
+ 61 59 62 60 61 61 60 62 59 63 58 64 57 65 56 66 \
+ 55 67 67 58 66 59 65 60 64 61 63 62 62 63 61 64 \
+ 65 51 66 50 67 76 68 77 69 78 70 75 57 74 56 68 \
+ 55 69 54 70 53 71 52 72 78 73 68 74 69 75 70 76 \
+ 77 77 76 78 75 71 74 72 73 73 72 74 71 75 70 76 \
+ 69 77 68 78 78 79 77 80 76 81 75 73 74 72 73 71 \
+ \
+ 67 66 66 67 65 65 64 64 63 63 62 62 61 61 60 60 \
+ 59 59 58 58 57 57 56 56 55 55 67 55 66 56 65 57 \
+ 64 58 63 59 62 60 61 61 60 62 59 63 58 64 57 65 \
+ 56 66 55 67 67 55 66 56 65 57 64 58 63 59 62 60 \
+ 72 69 71 70 70 71 69 68 68 72 78 73 77 74 76 75 \
+ 75 76 74 77 73 78 72 68 71 69 70 70 69 71 68 72 \
+ 83 78 84 77 82 76 81 75 80 74 79 73 85 83 86 82 \
+ 87 81 88 80 79 79 80 84 81 85 82 86 83 87 89 88 \
+ \
+ 59 56 58 57 57 58 56 59 55 60 67 61 66 55 65 58 \
+ 64 57 63 56 62 55 61 67 60 66 59 65 58 64 57 63 \
+ 56 62 55 61 67 60 66 72 65 71 64 70 63 69 62 68 \
+ 61 73 60 74 59 75 58 76 64 77 57 78 56 78 66 77 \
+ 78 68 77 69 76 70 75 71 74 72 73 73 72 74 71 75 \
+ 70 76 69 77 68 78 78 76 77 75 76 74 75 73 74 72 \
+ 73 71 72 70 71 69 86 68 87 87 88 86 89 85 85 84 \
+ 84 83 82 82 83 81 81 80 80 79 79 88 80 79 81 89 \
+ \
+ 66 67 65 66 64 65 63 64 62 63 60 62 58 61 56 60 \
+ 55 59 66 58 64 57 62 56 60 55 58 67 56 61 67 59 \
+ 65 57 63 67 61 65 60 63 58 61 56 59 67 57 65 55 \
+ 63 66 62 64 78 62 77 59 76 57 75 55 74 66 73 64 \
+ 72 70 71 69 70 68 69 78 68 77 79 76 79 75 80 74 \
+ 88 73 87 72 86 71 85 70 81 69 82 68 83 80 84 79 \
+ 85 89 86 88 88 87 87 86 79 85 80 84 81 83 89 82 \
+ 82 81 83 80 84 79 85 81 86 82 87 83 88 84 80 89 \
+ \
+ 64 61 63 60 62 59 61 58 60 57 59 56 58 55 57 67 \
+ 56 66 55 65 67 74 66 75 65 73 64 72 63 71 62 70 \
+ 61 76 60 77 59 78 58 69 57 68 56 78 55 77 73 76 \
+ 72 75 71 74 70 73 69 72 68 71 67 70 66 69 65 68 \
+ 72 71 71 72 70 73 69 74 68 75 78 76 77 77 76 78 \
+ 75 70 73 69 71 68 69 78 78 77 76 76 74 75 72 74 \
+ 70 73 68 74 77 72 75 70 74 68 72 77 70 75 69 73 \
+ 68 71 80 69 79 78 78 76 77 73 76 71 75 79 74 80 \
+ \
+ 96 93 96 93 95 92 95 91 94 90 94 89 88 87 86 85 \
+ 84 83 82 81 89 88 87 86 85 84 83 82 81 80 79 89 \
+ 94 88 94 87 94 86 93 85 92 84 91 83 90 82 92 81 \
+ 91 80 90 79 94 89 93 88 93 87 92 86 91 85 90 84 \
+ 94 83 -1 82 95 81 -1 80 -1 79 -1 89 95 88 95 87 \
+ -1 86 95 93 94 85 95 84 -1 83 96 82 96 81 96 80 \
+ 96 79 96 89 -1 88 95 87 96 86 -1 85 94 93 95 83 \
+ 96 82 96 81 94 88 84 86 89 85 87 83 84 82 89 81 \
+ \
+ 94 91 95 92 95 92 96 91 -1 93 -1 92 -1 90 -1 90 \
+ -1 -1 94 93 -1 -1 -1 -1 -1 -1 -1 -1 -1 92 -1 -1 \
+ 95 92 -1 91 -1 91 -1 90 -1 92 -1 92 -1 93 -1 -1 \
+ -1 89 -1 88 -1 87 -1 86 -1 85 93 84 93 83 93 82 \
+ -1 93 -1 92 -1 -1 -1 -1 -1 93 -1 -1 -1 -1 -1 -1 \
+ -1 -1 -1 -1 -1 -1 -1 -1 -1 91 -1 91 -1 90 -1 90 \
+ -1 90 -1 93 -1 91 -1 92 -1 92 94 91 96 90 94 91 \
+ 94 90 94 92 -1 91 95 90 96 90 95 91 96 92 94 93 \
+ \
+ 89 85 79 86 80 84 81 83 82 82 83 87 84 88 85 89 \
+ 86 79 87 80 88 81 89 82 79 83 80 84 81 85 82 86 \
+ 83 87 84 88 85 89 86 88 87 87 88 86 89 85 89 84 \
+ 84 83 85 82 86 81 87 80 88 79 89 89 86 88 85 87 \
+ 90 84 91 83 92 82 92 81 93 80 90 90 91 91 93 92 \
+ 90 93 90 93 91 94 92 94 93 95 94 95 94 96 94 96 \
+ 95 94 96 90 95 91 96 92 95 90 96 91 95 95 96 96 \
+ 95 96 95 93 96 92 -1 91 -1 90 -1 90 -1 91 79 92
+ChannelsToPadXRpc: Int_t \
+ 0 0 0 0 0 0 0 0 1 0 1 1 1 0 1 1 \
+ 1 1 1 1 2 0 2 2 2 2 2 2 2 2 1 3 \
+ 3 2 3 3 3 3 3 3 0 3 2 3 4 4 4 1 \
+ 4 4 4 4 4 4 3 4 2 5 1 0 5 4 5 5 \
+ 5 5 5 5 5 5 3 5 3 2 2 4 0 1 -3 -3 \
+ 6 6 6 6 6 6 6 6 6 6 4 6 3 5 1 2 \
+ 7 0 7 7 7 7 4 7 7 5 7 7 7 7 5 6 \
+ 0 1 3 2 8 6 8 8 8 8 3 8 8 4 8 8 \
+ \
+ 8 8 5 0 6 1 7 7 9 9 9 0 4 2 9 9 \
+ 9 3 9 9 9 9 4 9 3 4 0 2 5 1 5 5 \
+ 6 6 6 6 7 7 7 7 8 7 8 8 8 8 2 1 \
+ 4 3 6 5 0 7 8 1 9 9 9 2 3 8 5 4 \
+ 0 6 1 8 7 2 9 9 4 3 6 5 8 7 6 7 \
+ 4 5 2 3 0 1 8 8 9 9 9 9 4 9 6 5 \
+ 8 7 9 9 8 9 2 3 7 1 6 4 9 5 8 0 \
+ \
+ 0 0 2 1 8 9 6 7 0 5 1 4 1 1 9 3 \
+ 2 7 5 6 2 3 2 3 0 1 8 4 9 4 7 8 \
+ 5 6 8 9 7 1 2 6 4 5 9 3 7 8 5 6 \
+ 3 4 9 2 0 7 3 5 1 -3 4 -3 6 -3 8 -3 \
+ 9 2 3 8 5 4 7 6 8 9 6 7 4 5 1 3 \
+ 9 2 0 8 6 7 4 5 2 3 8 9 5 6 4 1 \
+ 7 3 2 9 -3 8 -3 7 -3 6 -3 5 -3 4 -3 3 \
+ -3 1 -3 2 8 9 -3 7 1 6 0 0 -3 -3 0 -3 \
+ 6 7 1 5 3 2 3 4 1 2 0 0 0 0 -3 -3 \
+ \
+ 0 0 1 0 3 2 5 4 8 9 7 0 5 6 4 1 \
+ 2 3 8 9 0 1 6 7 4 5 2 3 8 9 1 7 \
+ 5 6 3 4 9 2 7 8 1 6 2 5 4 3 8 9 \
+ 6 7 4 5 2 3 9 1 7 8 5 6 3 4 9 2 \
+ 7 8 5 6 3 4 2 0 9 1 7 8 5 6 3 4 \
+ 1 2 8 9 6 7 4 5 2 3 1 1 -3 -3 -3 -3 \
+ -3 -3 -3 -3 -3 -3 0 0 0 -3 0 -3 -3 -3 -3 -3 \
+ -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3
+ChannelsToPadYRpc: Int_t \
+ 1 0 3 2 5 4 7 6 0 8 2 1 3 9 5 4 \
+ 7 6 9 8 0 10 2 1 4 3 6 5 8 7 10 0 \
+ 1 9 3 2 5 4 7 6 11 8 10 9 1 0 2 11 \
+ 4 3 6 5 8 7 10 9 11 0 12 12 1 10 3 2 \
+ 5 4 7 6 9 8 11 10 12 12 13 11 13 13 -3 -3 \
+ 1 0 3 2 5 4 7 6 9 8 12 10 13 11 14 14 \
+ 0 14 2 1 4 3 13 5 6 12 8 7 10 9 13 11 \
+ 15 15 14 15 0 12 2 1 4 3 15 5 6 14 8 7 \
+ \
+ 10 9 14 16 13 16 11 12 1 0 3 2 17 15 16 4 \
+ 5 16 7 6 9 8 16 10 17 17 18 17 15 17 17 16 \
+ 15 14 17 16 15 14 17 16 12 13 14 13 16 15 18 18 \
+ 18 18 18 18 19 18 11 19 12 11 13 19 19 17 19 19 \
+ 20 19 20 18 19 20 15 14 20 20 20 20 20 20 21 21 \
+ 21 21 21 21 21 21 19 21 17 16 19 18 22 20 22 22 \
+ 22 22 22 21 23 23 22 22 23 22 23 23 24 23 24 22 \
+ 24 24 23 24 23 23 24 24 24 24 23 24 26 25 -3 -3 \
+ \
+ 27 28 25 25 25 25 25 25 29 25 26 25 28 27 26 25 \
+ 26 26 26 26 27 26 28 27 30 29 26 26 27 27 27 27 \
+ 27 27 28 28 28 30 29 28 28 28 29 28 29 29 29 29 \
+ 29 29 30 30 31 30 30 30 31 -3 30 -3 30 -3 30 -3 \
+ 31 31 31 31 31 31 31 31 32 32 32 32 32 32 32 32 \
+ 33 32 32 33 33 33 33 33 33 33 34 34 34 34 34 33 \
+ 34 34 34 35 -3 35 -3 35 -3 35 -3 35 -3 35 -3 35 \
+ -3 34 -3 35 36 36 -3 36 35 36 35 34 -3 -3 33 -3 \
+ \
+ 37 36 36 38 36 36 36 36 37 37 37 39 37 37 37 37 \
+ 37 37 37 37 40 38 38 38 38 38 38 38 39 39 39 39 \
+ 39 39 39 39 40 39 40 40 40 40 40 40 40 40 41 41 \
+ 41 41 41 41 41 41 42 41 42 42 42 42 42 42 43 42 \
+ 43 43 43 43 43 43 43 41 44 42 44 44 44 44 44 44 \
+ 43 44 45 45 45 45 45 45 45 45 44 45 -3 -3 -3 -3 \
+ -3 -3 -3 -3 -3 -3 44 45 43 -3 42 -3 -3 -3 -3 -3 \
+ -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3
+RealX: Double_t \
+ 184.654 187.877 191.156 194.492 197.887 201.34 204.855 208.43 212.068 215.769 219.535 223.366 227.265 231.231 235.267 239.373 243.551 247.802 252.127 256.527 261.005 265.56 270.195 274.911 279.709 284.591 28
+9.558 294.612 299.753 304.985 \
+ 310.308 315.724 321.234 326.841 332.546 338.35 344.255 350.263 356.376 362.596 368.925 375.364 381.915 388.581 395.363 402.263 409.284 416.427 423.695 431.09 438.614 446.27 454.058 461.983 470.046 478.25 48
+6.597 495.09 503.731 512.523 \
+ 521.468 530.569 539.829 549.251 558.837 568.591 578.515 588.612 598.885 609.338 619.973 630.793 641.802 653.004 664.401 675.997 687.795 699.8 712.014 724.441 737.084 749.949 763.038 776.356 789.906 803.692 81
+7.719 831.991 846.512 861.286 \
+ 876.319 891.613 907.175 923.008 939.118 955.508 972.185 184.626 187.848 191.127 194.463 197.857 201.31 204.823 208.398 212.035 215.736 219.501 223.332 227.23 231.196 235.231 239.337 243.514 247.764 252.089 25
+6.488 260.965 265.52 270.154 \
+ 274.869 279.666 284.547 289.514 294.567 299.708 304.939 310.261 315.676 321.186 326.791 332.495 338.298 344.202 350.21 356.322 362.541 368.869 375.307 381.857 388.522 395.303 402.202 409.222 416.364 423.631 43
+1.025 438.548 446.202 453.989 \
+ 461.913 469.975 478.177 486.523 495.015 503.654 512.445 521.389 530.488 539.747 549.168 558.752 568.504 578.427 588.522 598.794 609.245 619.878 630.697 641.705 652.905 664.3 675.894 687.691 699.693 711.905 7
+24.33 736.972 749.835 762.922 \
+ 776.237 789.785 803.57 817.595 831.864 846.383 861.155 876.185 891.478 907.037 922.867 938.975 955.363 972.037 184.542 187.762 191.039 194.374 197.766 201.218 204.73 208.303 211.939 215.638 219.401 223.23 22
+7.127 231.091 235.124 239.228 \
+ 243.403 247.651 251.973 256.371 260.846 265.398 270.03 274.743 279.539 284.417 289.381 294.432 299.571 304.799 310.119 315.532 321.039 326.642 332.343 338.143 344.045 350.05 356.159 362.376 368.7 375.135 38
+1.683 388.344 395.122 402.018 \
+ 409.035 416.174 423.437 430.828 438.347 445.998 453.782 461.702 469.76 477.959 486.301 494.788 503.424 512.211 521.15 530.246 539.501 548.917 558.497 568.245 578.162 588.253 598.52 608.966 619.595 630.409 64
+1.412 652.606 663.996 675.585 \
+ 687.376 699.373 711.58 723.999 736.635 749.492 762.573 775.883 789.424 803.202 817.221 831.484 845.996 860.762 875.785 891.07 906.622 922.446 938.546 954.926 971.593 184.401 187.619 190.894 194.226 197.616 20
+1.065 204.574 208.144 211.777 \
+ 215.473 219.234 223.06 226.954 230.915 234.945 239.045 243.218 247.462 251.781 256.176 260.647 265.196 269.825 274.534 279.326 284.201 289.161 294.208 299.343 304.567 309.883 315.291 320.794 326.393 332.09 33
+7.886 343.783 349.783 355.888 \
+ 362.099 368.419 374.849 381.392 388.048 394.821 401.712 408.723 415.857 423.115 430.5 438.013 445.658 453.436 461.35 469.402 477.595 485.93 494.412 503.041 511.82 520.753 529.842 539.09 548.499 558.072 56
+7.812 577.722 587.805 598.064 \
+ 608.502 619.123 629.929 640.923 652.109 663.491 675.071 686.853 698.841 711.038 723.448 736.074 748.921 761.992 775.292 788.823 802.591 816.598 830.851 845.352 860.106 875.118 890.391 905.932 921.743 937.831 95
+4.199 970.853 184.204 187.419 \
+ 190.69 194.018 197.405 200.85 204.356 207.922 211.551 215.243 219 222.822 226.711 230.668 234.694 238.79 242.958 247.198 251.513 255.903 260.369 264.913 269.537 274.241 279.028 283.897 288.852 293.894 29
+9.023 304.242 309.552 314.955 \
+ 320.452 326.045 331.735 337.525 343.416 349.41 355.508 361.713 368.026 374.449 380.985 387.634 394.4 401.283 408.287 415.413 422.663 430.04 437.546 445.182 452.952 460.858 468.901 477.085 485.412 493.884 50
+2.504 511.274 520.198 529.277 \
+ 538.514 547.913 557.476 567.206 577.106 587.178 597.426 607.853 618.462 629.256 640.239 651.413 662.783 674.35 686.12 698.095 710.279 722.676 735.289 748.122 761.179 774.464 787.981 801.734 815.727 829.964 8
+44.45 859.188 874.184 889.441 \
+ 904.965 920.76 936.83 953.181 969.817 183.951 187.162 190.428 193.752 197.134 200.574 204.075 207.637 211.261 214.948 218.699 222.517 226.4 230.352 234.372 238.463 242.624 246.859 251.168 255.551 260.012 2
+64.55 269.167 273.865 278.645 \
+ 283.508 288.456 293.49 298.613 303.825 309.127 314.523 320.012 325.597 331.28 337.062 342.945 348.93 355.02 361.217 367.521 373.935 380.462 387.102 393.858 400.733 407.727 414.843 422.083 429.45 436.945 44
+4.571 452.331 460.225 468.258 \
+ 476.43 484.746 493.206 501.814 510.572 519.484 528.55 537.775 547.161 556.711 566.427 576.313 586.372 596.606 607.019 617.613 628.393 639.36 650.519 661.873 673.425 685.178 697.137 709.304 721.684 734.28 74
+7.095 760.134 773.401 786.9 \
+ 800.634 814.607 828.825 843.291 858.009 872.984 888.22 903.723 919.496 935.544 951.872 968.486 183.642 186.848 190.109 193.427 196.803 200.238 203.732 207.288 210.906 214.587 218.332 222.143 226.02 229.965 23
+3.978 238.062 242.217 246.445 \
+ 250.746 255.122 259.575 264.105 268.715 273.405 278.177 283.032 287.972 292.998 298.111 303.314 308.608 313.994 319.475 325.051 330.724 336.496 342.369 348.344 354.424 360.61 366.904 373.308 379.823 386.452 39
+3.197 400.06 407.042 414.146 \
+ 421.374 428.729 436.212 443.825 451.571 459.452 467.471 475.63 483.932 492.378 500.971 509.715 518.611 527.663 536.872 546.242 555.776 565.476 575.346 585.387 595.604 606 616.576 627.338 638.287 649.427 66
+0.761 672.294 684.028 695.966 \
+ 708.113 720.472 733.047 745.841 758.858 772.103 785.578 799.289 813.24 827.433 841.875 856.568 871.518 886.729 902.205 917.952 933.973 950.274 966.859 183.278 186.476 189.731 193.042 196.412 199.84 203.328 20
+6.876 210.487 214.161 217.899 \
+ 221.702 225.571 229.508 233.514 237.589 241.736 245.955 250.248 254.615 259.059 263.581 268.181 272.862 277.624 282.469 287.399 292.416 297.519 302.712 307.995 313.371 318.84 324.405 330.067 335.828 341.689 34
+7.652 353.72 359.894 366.175 \
+ 372.566 379.068 385.684 392.416 399.265 406.233 413.323 420.537 427.877 435.345 442.943 450.674 458.54 466.543 474.685 482.97 491.4 499.976 508.702 517.581 526.614 535.806 545.157 554.672 564.353 574.203 58
+4.224 594.421 604.796 615.351 \
+ 626.091 637.019 648.137 659.449 670.958 682.669 694.584 706.706 719.041 731.59 744.359 757.35 770.569 784.018 797.701 811.624 825.789 840.202 854.866 869.787 884.967 900.413 916.128 932.118 948.386 964.939 18
+2.857 186.048 189.296 192.599 \
+ 195.961 199.381 202.861 206.402 210.004 213.669 217.398 221.193 225.053 228.981 232.978 237.044 241.181 245.39 249.673 254.031 258.465 262.976 267.565 272.235 276.987 281.821 286.74 291.744 296.836 302.017 30
+7.288 312.651 318.108 323.66 \
+ 329.309 335.057 340.905 346.854 352.908 359.068 365.335 371.711 378.198 384.799 391.515 398.348 405.301 412.375 419.572 426.895 434.346 441.927 449.64 457.487 465.472 473.596 481.862 490.272 498.829 507.535 51
+6.393 525.406 534.576 543.906 \
+ 553.399 563.058 572.885 582.883 593.057 603.408 613.939 624.654 635.556 646.649 657.935 669.418 681.102 692.989 705.084 717.39 729.911 742.651 755.612 768.8 782.218 795.871 809.761 823.894 838.274 852.904 8
+67.79 882.936 898.346 914.025 \
+ 929.978 946.209 962.724 182.381 185.564 188.802 192.098 195.45 198.862 202.332 205.864 209.457 213.113 216.832 220.616 224.467 228.385 232.371 236.426 240.553 244.751 249.023 253.369 257.791 262.291 266.868 27
+1.526 276.265 281.087 285.993 \
+ 290.984 296.063 301.23 306.488 311.837 317.28 322.817 328.451 334.184 340.016 345.951 351.989 358.132 364.383 370.743 377.213 383.797 390.495 397.311 404.245 411.301 418.479 425.783 433.214 440.775 448.468 45
+6.295 464.259 472.362 480.606 \
+ 488.995 497.529 506.213 515.048 524.037 533.183 542.489 551.957 561.591 571.392 581.365 591.512 601.836 612.34 623.027 633.901 644.964 656.221 667.674 679.328 691.184 703.248 715.521 728.01 740.716 753.644 76
+6.797 780.181 793.797 807.652 \
+ 821.748 836.09 850.683 865.53 880.636 896.006 911.644 927.556 943.744 960.216 181.849 185.023 188.252 191.537 194.88 198.282 201.742 205.263 208.846 212.491 216.2 219.973 223.812 227.719 231.693 235.737 23
+9.851 244.037 248.297 252.63 \
+ 257.039 261.526 266.09 270.734 275.459 280.267 285.159 290.136 295.2 300.352 305.594 310.927 316.354 321.876 327.493 333.209 339.025 344.942 350.962 357.088 363.32 369.661 376.113 382.677 389.356 396.152 40
+3.066 410.101 417.259 424.541 \
+ 431.951 439.49 447.16 454.965 462.905 470.985 479.205 487.568 496.078 504.736 513.546 522.509 531.628 540.907 550.347 559.953 569.726 579.669 589.787 600.08 610.554 621.21 632.052 643.083 654.307 665.727 67
+7.346 689.168 701.196 713.435 \
+ 725.886 738.556 751.446 764.561 777.905 791.482 805.296 819.351 833.652 848.201 863.005 878.068 893.393 908.986 924.85 940.992 957.415 181.261 184.425 187.644 190.919 194.251 197.641 201.091 204.6 208.171 21
+1.805 215.501 219.263 223.089 \
+ 226.983 230.945 234.975 239.077 243.249 247.495 251.814 256.209 260.681 265.231 269.86 274.57 279.362 284.238 289.199 294.246 299.382 304.607 309.923 315.332 320.836 326.436 332.133 337.93 343.828 349.829 35
+5.934 362.147 368.467 374.898 \
+ 381.442 388.099 394.873 401.764 408.777 415.911 423.17 430.556 438.07 445.716 453.495 461.41 469.463 477.657 485.994 494.476 503.106 511.887 520.821 529.911 539.16 548.57 558.144 567.886 577.797 587.882 59
+8.142 608.582 619.204 630.011 \
+ 641.007 652.194 663.577 675.159 686.942 698.932 711.131 723.542 736.17 749.019 762.092 775.393 788.926 802.695 816.705 830.959 845.462 860.218 875.232 890.508 906.05 921.863 937.953 954.323 180.619 183.771 18
+6.979 190.242 193.562 196.941 \
+ 200.378 203.875 207.434 211.054 214.738 218.485 222.299 226.179 230.126 234.143 238.229 242.387 246.617 250.922 255.301 259.757 264.291 268.903 273.597 278.372 283.23 288.174 293.203 298.32 303.527 308.825 31
+4.215 319.699 325.279 330.956 \
+ 336.732 342.609 348.589 354.673 360.863 367.161 373.569 380.089 386.723 393.473 400.34 407.328 414.437 421.67 429.03 436.518 444.136 451.888 459.775 467.799 475.964 484.271 492.723 501.323 510.073 518.975 52
+8.033 537.249 546.626 556.166 \
+ 565.873 575.749 585.798 596.022 606.425 617.009 627.778 638.734 649.882 661.225 672.765 684.507 696.454 708.61 720.977 733.561 746.364 759.39 772.644 786.129 799.85 813.81 828.014 842.465 857.169 872.129 88
+7.351 902.838 918.596 934.628 \
+ 950.94 179.921 183.062 186.257 189.507 192.815 196.18 199.604 203.088 206.632 210.239 213.908 217.642 221.44 225.305 229.237 233.238 237.309 241.451 245.665 249.953 254.315 258.754 263.27 267.865 272.54 27
+7.297 282.136 287.061 292.071 \
+ 297.168 302.355 307.632 313.001 318.464 324.022 329.678 335.432 341.286 347.243 353.303 359.469 365.743 372.127 378.622 385.23 391.953 398.794 405.754 412.836 420.042 427.373 434.832 442.421 450.143 457.999 46
+5.993 474.126 482.401 490.82 \
+ 499.387 508.103 516.971 525.994 535.174 544.515 554.018 563.688 573.526 583.536 593.72 604.083 614.626 625.353 636.268 647.373 658.671 670.167 681.864 693.765 705.873 718.193 730.728 743.481 756.458 769.66 78
+3.093 796.761 810.667 824.816 \
+ 839.212 853.859 868.761 883.924 899.351 915.048 931.019 947.268 179.169 182.296 185.478 188.715 192.009 195.36 198.769 202.239 205.768 209.36 213.014 216.732 220.514 224.363 228.279 232.263 236.317 240.441 24
+4.638 248.908 253.252 257.672 \
+ 262.169 266.745 271.4 276.137 280.957 285.86 290.849 295.926 301.091 306.346 311.692 317.132 322.667 328.299 334.029 339.859 345.791 351.826 357.966 364.214 370.571 377.038 383.619 390.314 397.127 404.058 4
+11.11 418.285 425.586 433.013 \
+ 440.571 448.26 456.084 464.044 472.143 480.384 488.768 497.299 505.978 514.809 523.794 532.936 542.238 551.701 561.33 571.128 581.096 591.238 601.557 612.056 622.738 633.607 644.666 655.917 667.365 679.013 69
+0.864 702.922 715.19 727.672 \
+ 740.373 753.295 766.442 779.819 793.429 807.277 821.367 835.702 850.288 865.129 880.228 895.591 911.222 927.126 943.307 178.362 181.475 184.642 187.865 191.144 194.48 197.874 201.328 204.842 208.417 212.054 21
+5.755 219.521 223.352 227.251 \
+ 231.217 235.252 239.358 243.536 247.787 252.111 256.511 260.988 265.543 270.178 274.894 279.691 284.573 289.54 294.593 299.735 304.966 310.289 315.704 321.214 326.821 332.525 338.328 344.233 350.241 356.354 36
+2.574 368.902 375.34 381.891 \
+ 388.556 395.338 402.238 409.258 416.401 423.669 431.063 438.587 446.242 454.03 461.954 470.017 478.22 486.567 495.059 503.699 512.491 521.435 530.536 539.796 549.217 558.802 568.555 578.479 588.575 598.847 60
+9.299 619.934 630.753 641.762 \
+ 652.963 664.359 675.955 687.752 699.756 711.969 724.395 737.038 749.902 762.99 776.307 789.856 803.642 817.668 831.939 846.459 861.232 876.264 891.557 907.118 922.95 939.059 177.501 180.599 183.751 186.958 19
+0.221 193.541 196.919 200.356 \
+ 203.853 207.41 211.03 214.714 218.461 222.274 226.153 230.101 234.117 238.203 242.36 246.59 250.894 255.273 259.728 264.261 268.873 273.566 278.341 283.199 288.142 293.171 298.287 303.493 308.79 314.18 31
+9.663 325.242 330.919 336.695 \
+ 342.571 348.55 354.633 360.823 367.12 373.528 380.047 386.68 393.429 400.296 407.282 414.391 421.623 428.982 436.469 444.087 451.838 459.724 467.747 475.911 484.217 492.668 501.267 510.016 518.917 527.974 53
+7.189 546.565 556.104 565.81 \
+ 575.685 585.733 595.956 606.357 616.94 627.708 638.663 649.81 661.151 672.691 684.431 696.377 708.531 720.897 733.479 746.281 759.306 772.558 786.042 799.761 813.72 827.922 842.372 857.074 872.032 887.252 90
+2.738 918.494 934.524 176.585 \
+ 179.667 182.803 185.994 189.24 192.543 195.903 199.323 202.801 206.341 209.942 213.606 217.335 221.128 224.987 228.914 232.909 236.974 241.11 245.318 249.6 253.956 258.389 262.899 267.487 272.155 276.905 28
+1.738 286.656 291.659 296.749 \
+ 301.928 307.198 312.56 318.015 323.565 329.213 334.958 340.805 346.753 352.805 358.962 365.227 371.602 378.087 384.686 391.4 398.232 405.182 412.254 419.449 426.77 434.218 441.797 449.508 457.353 465.335 47
+3.457 481.72 490.128 498.682 \
+ 507.386 516.241 525.252 534.419 543.746 553.236 562.892 572.717 582.712 592.883 603.23 613.759 624.471 635.37 646.459 657.742 669.222 680.902 692.786 704.877 717.18 729.697 742.433 755.39 768.574 781.989 79
+5.637 809.523 823.652 838.028 \
+ 852.654 867.536 882.677 898.083 913.757 929.705 175.616 178.681 181.8 184.973 188.201 191.486 194.828 198.229 201.688 205.209 208.79 212.434 216.142 219.914 223.752 227.658 231.631 235.674 239.787 243.972 2
+48.23 252.563 256.971 261.456 \
+ 266.019 270.662 275.386 280.192 285.082 290.058 295.121 300.271 305.512 310.844 316.27 321.79 327.406 333.12 338.934 344.85 350.868 356.992 363.223 369.562 376.012 382.575 389.252 396.046 402.958 409.991 41
+7.147 424.428 431.835 439.372 \
+ 447.041 454.843 462.782 470.859 479.077 487.438 495.945 504.601 513.408 522.369 531.486 540.762 550.2 559.803 569.574 579.514 589.629 599.92 610.39 621.044 631.883 642.911 654.132 665.549 677.165 688.984 70
+1.009 713.244 725.692 738.358 \
+ 751.245 764.357 777.697 791.27 805.081 819.132 833.429 847.975 862.775 877.833 893.154 908.742 924.603 174.594 177.641 180.741 183.896 187.106 190.371 193.694 197.074 200.514 204.014 207.574 211.197 214.883 21
+8.634 222.45 226.332 230.282 \
+ 234.301 238.391 242.551 246.785 251.092 255.474 259.933 264.47 269.086 273.782 278.561 283.422 288.369 293.402 298.523 303.733 309.034 314.428 319.916 325.499 331.18 336.961 342.842 348.825 354.914 361.108 3
+67.41 373.823 380.347 386.986 \
+ 393.74 400.612 407.604 414.718 421.956 429.321 436.814 444.438 452.194 460.087 468.117 476.287 484.6 493.058 501.663 510.419 519.327 528.391 537.613 546.997 556.543 566.257 576.14 586.196 596.427 606.836 61
+7.427 628.204 639.168 650.323 \
+ 661.674 673.222 684.972 696.927 709.091 721.467 734.059 746.87 759.906 773.169 786.663 800.393 814.362 828.576 843.037 857.751 872.721 887.953 903.451 919.219 173.518 176.546 179.628 182.763 185.953 189.198
+192.5 195.86 199.278 202.757 \
+ 206.295 209.896 213.559 217.287 221.079 224.937 228.863 232.858 236.922 241.057 245.264 249.545 253.9 258.332 262.84 267.428 272.095 276.844 281.676 286.592 291.594 296.684 301.862 307.13 312.491 317.945 32
+3.494 329.14 334.884 340.729 \
+ 346.676 352.727 358.883 365.147 371.52 378.004 384.601 391.314 398.144 405.092 412.163 419.356 426.675 434.122 441.699 449.408 457.252 465.232 473.352 481.614 490.02 498.572 507.274 516.127 525.135 534.301 54
+3.626 553.114 562.768 572.59 \
+ 582.584 592.752 603.097 613.623 624.333 635.229 646.316 657.597 669.074 680.751 692.633 704.722 717.021 729.536 742.268 755.223 768.405 781.816 795.461 809.344 823.47 837.842 852.466 867.344 882.482 897.884 91
+3.555 172.389 175.398 178.459 \
+ 181.574 184.743 187.968 191.248 194.586 197.982 201.438 204.953 208.531 212.17 215.873 219.641 223.474 227.375 231.343 235.381 239.489 243.669 247.922 252.249 256.651 261.131 265.688 270.325 275.044 279.844 28
+4.728 289.698 294.754 299.898 \
+ 305.132 310.458 315.876 321.39 326.999 332.706 338.513 344.421 350.432 356.549 362.772 369.103 375.545 382.1 388.769 395.554 402.458 409.482 416.629 423.9 431.299 438.826 446.485 454.278 462.206 470.273 47
+8.481 486.832 495.329 503.974 \
+ 512.77 521.72 530.825 540.09 549.516 559.107 568.866 578.794 588.896 599.174 609.632 620.272 631.098 642.112 653.319 664.722 676.324 688.128 700.138 712.357 724.79 737.44 750.311 763.407 776.73 790.287 8
+04.08 818.114 832.393 846.921 \
+ 861.702 876.742 892.044 907.613 171.208 174.196 177.237 180.33 183.477 186.68 189.938 193.253 196.626 200.058 203.549 207.102 210.716 214.394 218.136 221.943 225.817 229.758 233.768 237.848 241.999 246.223 2
+50.52 254.893 259.342 263.868 \
+ 268.473 273.159 277.927 282.777 287.713 292.734 297.843 303.042 308.331 313.712 319.188 324.758 330.427 336.194 342.061 348.031 354.106 360.286 366.574 372.972 379.482 386.105 392.844 399.7 406.676 413.774 42
+0.996 428.343 435.819 443.426 \
+ 451.165 459.039 467.051 475.203 483.497 491.935 500.521 509.257 518.145 527.188 536.39 545.751 555.277 564.968 574.828 584.861 595.069 605.455 616.022 626.774 637.713 648.843 660.167 671.69 683.413 695.341 70
+7.477 719.824 732.388 745.17 \
+ 758.176 771.409 784.872 798.571 812.509 826.689 841.118 855.798 870.735 885.932 901.394
+PadSize: Double_t \
+ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 \
+ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 \
+ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 \
+ 0 0 0 0 0 0 0 3.22266 3.2789 3.33613 3.39436 3.4536 3.51388 3.5752 3.6376 3.70109 3.76569 3.83141 3.89828 3.96632 4.03555 4.10598 4.17764 4.25056 4.32474 4.40022 4.
+47702 4.55516 4.63466 4.71555 \
+ 4.79785 4.88159 4.96679 5.05348 5.14168 5.23142 5.32272 5.41562 5.51014 5.60631 5.70416 5.80372 5.90501 6.00808 6.11294 6.21963 6.32818 6.43863 6.551 6.66534 6.78167 6.90003 7.02046 7.14299 7.26766 7.39451 7.
+52356 7.65488 7.78848 7.92441 \
+ 8.06272 8.20344 8.34662 8.49229 8.64051 8.79132 8.94475 9.10087 9.25971 9.42132 9.58576 9.75306 9.92328 10.0965 10.2727 10.452 10.6344 10.82 11.0089 11.201 11.3965 11.5954 11.7978 12.0037 12.2132 12.4264 12
+.6432 12.8639 13.0884 13.3169 \
+ 13.5493 13.7858 14.0264 14.2712 14.5202 14.7737 15.0315 15.2939 15.5608 15.8324 16.1087 16.3899 16.6759 16.967 6.44433 6.55681 6.67124 6.78768 6.90615 7.02668 7.14932 7.2741 7.40106 7.53023 7.66166 7.79538 7.
+93143 8.06986 8.21071 8.35401 \
+ 8.49982 8.64817 8.79911 8.95268 9.10893 9.26791 9.42967 9.59425 9.7617 9.93207 10.1054 10.2818 10.4612 10.6438 10.8296 11.0186 11.2109 11.4066 11.6057 11.8082 12.0143 12.224 12.4374 12.6544 12.8753 13.1 13
+.3286 13.5613 13.798 14.0388 \
+ 14.2838 14.5331 14.7868 15.0448 15.3074 15.5746 15.8464 16.123 16.4044 16.6907 16.982 17.2784 17.58 17.8868 18.199 18.5166 18.8398 19.1686 19.5031 19.8435 20.1899 20.5423 20.9008 21.2656 21.6367 22.0144 22
+.3986 22.7895 23.1873 23.592 \
+ 24.0037 24.4227 24.8489 25.2826 25.7239 26.1728 26.6296 27.0944 27.5673 28.0484 28.538 29.0361 29.5428 30.0585 30.5831 31.1169 31.6599 32.2125 32.7747 33.3468 33.9288 9.66404 9.83271 10.0043 10.1789 10.3566 10
+.5373 10.7213 10.9084 11.0988 \
+ 11.2925 11.4896 11.6901 11.8941 12.1017 12.3129 12.5278 12.7465 12.969 13.1953 13.4256 13.6599 13.8983 14.1409 14.3877 14.6388 14.8943 15.1543 15.4188 15.6879 15.9617 16.2403 16.5237 16.8121 17.1055 17.4041 17
+.7078 18.0169 18.3314 18.6513 \
+ 18.9768 19.308 19.645 19.9879 20.3368 20.6917 21.0528 21.4203 21.7941 22.1745 22.5615 22.9553 23.3559 23.7636 24.1783 24.6003 25.0297 25.4665 25.911 26.3632 26.8234 27.2915 27.7679 28.2525 28.7456 29.2473 29
+.7578 30.2771 30.8056 31.3432 \
+ 31.8903 32.4469 33.0132 33.5893 34.1756 34.7721 35.379 35.9964 36.6247 37.2639 37.9143 38.576 39.2493 39.9343 40.6313 41.3405 42.062 42.7961 43.543 44.303 45.0762 45.863 46.6634 47.4779 48.3065 49.1496 50
+.0074 50.8802 12.8808 13.1056 \
+ 13.3344 13.5671 13.8039 14.0448 14.2899 14.5393 14.7931 15.0513 15.314 15.5813 15.8532 16.1299 16.4114 16.6978 16.9893 17.2858 17.5875 17.8945 18.2068 18.5245 18.8479 19.1768 19.5115 19.852 20.1985 20.5511 20
+.9097 21.2747 21.646 22.0238 \
+ 22.4082 22.7993 23.1972 23.6021 24.014 24.4331 24.8596 25.2934 25.7349 26.1841 26.6411 27.106 27.5791 28.0605 28.5502 29.0485 29.5555 30.0713 30.5962 31.1302 31.6735 32.2263 32.7888 33.3611 33.9433 34.5357 35
+.1385 35.7518 36.3758 37.0106 \
+ 37.6566 38.3138 38.9825 39.6629 40.3552 41.0595 41.7761 42.5052 43.2471 44.0019 44.7699 45.5513 46.3463 47.1552 47.9782 48.8156 49.6676 50.5344 51.4164 52.3138 53.2268 54.1558 55.101 56.0627 57.0412 58.0368 59
+.0497 60.0803 61.1289 62.1958 \
+ 63.2813 64.3858 65.5095 66.6529 67.8162 16.0937 16.3745 16.6603 16.9511 17.247 17.548 17.8543 18.1659 18.4829 18.8055 19.1337 19.4677 19.8074 20.1532 20.5049 20.8628 21.2269 21.5974 21.9743 22.3578 22.7481 23
+.1451 23.549 23.9601 24.3782 \
+ 24.8037 25.2366 25.6771 26.1252 26.5812 27.0451 27.5172 27.9974 28.4861 28.9833 29.4891 30.0038 30.5275 31.0603 31.6024 32.1539 32.7151 33.2861 33.8671 34.4581 35.0596 35.6715 36.294 36.9275 37.572 38.2278 3
+8.895 39.5738 40.2645 40.9672 \
+ 41.6823 42.4097 43.1499 43.903 44.6693 45.4489 46.2422 47.0492 47.8704 48.7059 49.556 50.4209 51.3009 52.1963 53.1073 54.0342 54.9772 55.9368 56.9131 57.9064 58.917 59.9453 60.9916 62.0561 63.1392 64.2411 65
+.3624 66.5031 67.6638 68.8448 \
+ 70.0464 71.2689 72.5128 73.7784 75.0661 76.3762 77.7092 79.0655 80.4455 81.8495 83.278 84.7315 19.3016 19.6385 19.9812 20.33 20.6848 21.0458 21.4131 21.7869 22.1671 22.554 22.9476 23.3482 23.7557 24.1703 24
+.5921 25.0213 25.458 25.9024 \
+ 26.3544 26.8144 27.2824 27.7586 28.2431 28.736 29.2375 29.7478 30.267 30.7953 31.3328 31.8796 32.436 33.0021 33.5781 34.1642 34.7605 35.3672 35.9844 36.6125 37.2515 37.9016 38.5632 39.2362 39.921 40.6178 41
+.3267 42.048 42.7818 43.5285 \
+ 44.2882 45.0612 45.8477 46.6479 47.462 48.2904 49.1332 49.9908 50.8633 51.751 52.6542 53.5732 54.5082 55.4596 56.4275 57.4124 58.4144 59.4339 60.4713 61.5267 62.6005 63.6931 64.8048 65.9358 67.0866 68.2575 69
+.4488 70.6609 71.8942 73.149 \
+ 74.4257 75.7247 77.0463 78.391 79.7592 81.1513 82.5676 84.0087 85.4749 86.9667 88.4846 90.0289 91.6002 93.199 94.8256 96.4806 98.1645 99.8778 101.621 22.5037 22.8964 23.296 23.7026 24.1163 24.5372 24.9655 25
+.4012 25.8446 26.2956 26.7546 \
+ 27.2215 27.6966 28.18 28.6719 29.1723 29.6814 30.1995 30.7266 31.2628 31.8085 32.3636 32.9285 33.5032 34.0879 34.6829 35.2882 35.9041 36.5308 37.1683 37.817 38.4771 39.1486 39.8319 40.5271 41.2344 41.9541 42
+.6863 43.4314 44.1894 44.9606 \
+ 45.7453 46.5438 47.3561 48.1826 49.0236 49.8792 50.7497 51.6355 52.5367 53.4536 54.3866 55.3358 56.3016 57.2842 58.284 59.3013 60.3363 61.3894 62.4608 63.551 64.6601 65.7887 66.9369 68.1052 69.2938 70.5032 71
+.7337 72.9857 74.2596 75.5556 \
+ 76.8743 78.216 79.5812 80.9701 82.3833 83.8212 85.2841 86.7726 88.2871 89.828 91.3958 92.9909 94.6139 96.2653 97.9454 99.6549 101.394 103.164 104.964 106.796 108.66 110.557 112.486 114.45 116.447 118.48 25
+.6989 26.1474 26.6038 27.0681 \
+ 27.5405 28.0212 28.5102 29.0078 29.5141 30.0292 30.5534 31.0866 31.6292 32.1812 32.7429 33.3143 33.8958 34.4874 35.0893 35.7017 36.3248 36.9588 37.6039 38.2602 38.928 39.6074 40.2987 41.002 41.7176 42.4457 43
+.1865 43.9403 44.7072 45.4875 \
+ 46.2814 47.0892 47.911 48.7472 49.598 50.4637 51.3444 52.2406 53.1523 54.08 55.0239 55.9842 56.9613 57.9555 58.967 59.9962 61.0433 62.1087 63.1927 64.2956 65.4178 66.5596 67.7213 68.9032 70.1058 71.3294 72
+.5743 73.841 75.1297 76.441 \
+ 77.7751 79.1326 80.5137 81.9189 83.3487 84.8034 86.2835 87.7894 89.3216 90.8806 92.4668 94.0806 95.7226 97.3933 99.0931 100.823 102.582 104.373 106.194 108.048 109.934 111.852 113.804 115.791 117.812 119.868 1
+21.96 124.089 126.254 128.458 \
+ 130.7 132.981 135.302 28.8862 29.3904 29.9034 30.4253 30.9563 31.4966 32.0463 32.6056 33.1747 33.7537 34.3428 34.9422 35.5521 36.1726 36.8039 37.4463 38.0998 38.7648 39.4414 40.1297 40.8301 41.5428 42.2678 43
+.0055 43.7561 44.5198 45.2968 \
+ 46.0874 46.8918 47.7102 48.5429 49.3901 50.2521 51.1292 52.0216 52.9295 53.8533 54.7932 55.7496 56.7226 57.7126 58.7198 59.7447 60.7874 61.8484 62.9278 64.0261 65.1436 66.2806 67.4374 68.6144 69.8119 71.0304 72
+.2701 73.5315 74.8148 76.1206 \
+ 77.4491 78.8009 80.1762 81.5756 82.9993 84.4479 85.9218 87.4214 88.9472 90.4997 92.0792 93.6863 95.3214 96.9851 98.6778 100.4 102.152 103.935 105.749 107.595 109.473 111.383 113.327 115.305 117.318 119.365 12
+1.449 123.568 125.725 127.919 \
+ 130.152 132.424 134.735 137.086 139.479 141.913 144.39 146.91 149.474 152.083 32.0648 32.6245 33.1939 33.7732 34.3627 34.9624 35.5726 36.1935 36.8252 37.4679 38.1218 38.7872 39.4641 40.1529 40.8537 41.5668 42
+.2922 43.0304 43.7814 44.5455 \
+ 45.323 46.114 46.9189 47.7378 48.5709 49.4187 50.2812 51.1588 52.0516 52.9601 53.8844 54.8249 55.7818 56.7554 57.7459 58.7538 59.7792 60.8226 61.8841 62.9642 64.0631 65.1813 66.3189 67.4764 68.654 69.8523 71
+.0714 72.3119 73.5739 74.8581 \
+ 76.1646 77.4939 78.8464 80.2226 81.6227 83.0473 84.4967 85.9715 87.472 88.9986 90.552 92.1324 93.7404 95.3765 97.0411 98.7348 100.458 102.211 103.995 105.81 107.657 109.536 111.448 113.393 115.372 117.386 11
+9.434 121.519 123.64 125.798 \
+ 127.993 130.227 132.5 134.813 137.166 139.56 141.995 144.474 146.995 149.561 152.171 154.827 157.529 160.279 163.076 165.922 168.818 35.2336 35.8486 36.4743 37.1109 37.7586 38.4176 39.0881 39.7703 40.4644 41
+.1707 41.8892 42.6203 43.3642 \
+ 44.121 44.8911 45.6746 46.4718 47.2829 48.1081 48.9477 49.802 50.6713 51.5556 52.4555 53.371 54.3025 55.2502 56.2145 57.1957 58.1939 59.2096 60.243 61.2944 62.3642 63.4527 64.5601 65.6869 66.8334 67.9998 69
+.1867 70.3942 71.6228 72.8729 \
+ 74.1447 75.4388 76.7555 78.0951 79.4581 80.8449 82.2559 83.6916 85.1522 86.6384 88.1506 89.6891 91.2544 92.8471 94.4676 96.1164 97.7939 99.5008 101.237 103.004 104.802 106.631 108.492 110.386 112.312 114.273 11
+6.267 118.296 120.361 122.462 \
+ 124.599 126.774 128.986 131.238 133.528 135.859 138.23 140.642 143.097 145.595 148.136 150.721 153.352 156.028 158.751 161.522 164.341 167.209 170.128 173.097 176.118 179.192 182.32 185.502 38.3917 39.0618 39
+.7435 40.4372 41.143 41.861 \
+ 42.5917 43.335 44.0914 44.8609 45.6439 46.4405 47.251 48.0757 48.9148 49.7685 50.6372 51.5209 52.4202 53.3351 54.2659 55.2131 56.1767 57.1572 58.1547 59.1697 60.2024 61.2532 62.3222 63.41 64.5167 65.6427 66
+.7884 67.9541 69.1401 70.3468 \
+ 71.5746 72.8238 74.0948 75.388 76.7038 78.0425 79.4046 80.7905 82.2006 83.6352 85.0949 86.5801 88.0912 89.6287 91.193 92.7847 94.4041 96.0517 97.7281 99.4338 101.169 102.935 104.732 106.559 108.419 110.312 11
+2.237 114.196 116.189 118.217 \
+ 120.28 122.379 124.515 126.688 128.9 131.149 133.438 135.767 138.137 140.548 143.001 145.497 148.036 150.62 153.248 155.923 158.645 161.413 164.231 167.097 170.013 172.981 176 179.072 182.197 185.377 18
+8.612 191.904 195.254 198.661 \
+ 202.129 41.5381 42.2631 43.0007 43.7512 44.5148 45.2918 46.0822 46.8865 47.7049 48.5375 49.3846 50.2465 51.1235 52.0158 52.9236 53.8473 54.7871 55.7433 56.7162 57.7061 58.7133 59.738 60.7807 61.8415 62.9208 6
+4.019 65.1363 66.2732 67.4299 \
+ 68.6067 69.8041 71.0225 72.262 73.5232 74.8065 76.1121 77.4405 78.7921 80.1673 81.5664 82.99 84.4385 85.9122 87.4117 88.9373 90.4895 92.0689 93.6758 95.3107 96.9742 98.6668 100.389 102.141 103.924 105.737 10
+7.583 109.461 111.371 113.315 \
+ 115.293 117.305 119.352 121.435 123.555 125.711 127.905 130.138 132.409 134.72 137.071 139.463 141.898 144.374 146.894 149.458 152.066 154.72 157.421 160.168 162.964 165.808 168.702 171.646 174.642 177.69 18
+0.791 183.947 187.157 190.424 \
+ 193.747 197.129 200.569 204.07 207.632 211.255 214.943 218.694 44.6718 45.4515 46.2448 47.0519 47.8731 48.7087 49.5588 50.4238 51.3038 52.1992 53.1103 54.0372 54.9804 55.94 56.9163 57.9097 58.9204 59.9487 6
+0.995 62.0596 63.1428 64.2448 \
+ 65.3661 66.5069 67.6677 68.8487 70.0504 71.273 72.5169 73.7826 75.0703 76.3806 77.7136 79.07 80.45 81.8542 83.2828 84.7363 86.2153 87.72 89.251 90.8087 92.3936 94.0062 95.6469 97.3163 99.0148 100.743 10
+2.501 104.29 106.11 107.962 \
+ 109.847 111.764 113.715 115.699 117.719 119.773 121.864 123.99 126.155 128.356 130.597 132.876 135.195 137.555 139.955 142.398 144.883 147.412 149.985 152.603 155.266 157.976 160.733 163.539 166.393 169.297 17
+2.252 175.258 178.317 181.429 \
+ 184.596 187.817 191.095 194.431 197.824 201.277 204.79 208.364 212.001 215.701 219.465 223.296 227.193 231.158 235.193 47.792 48.6261 49.4748 50.3383 51.2169 52.1108 53.0203 53.9456 54.8872 55.8451 56.8198 57
+.8115 58.8205 59.8471 60.8916 \
+ 61.9544 63.0357 64.1359 65.2553 66.3942 67.553 68.732 69.9316 71.1521 72.394 73.6575 74.9431 76.2511 77.5819 78.936 80.3137 81.7154 83.1416 84.5927 86.0691 87.5713 89.0997 90.6548 92.237 93.8469 95.4848 97
+.1513 98.8469 100.572 102.327 \
+ 104.113 105.931 107.779 109.66 111.574 113.522 115.503 117.519 119.57 121.657 123.78 125.941 128.139 130.375 132.651 134.966 137.321 139.718 142.157 144.638 147.162 149.731 152.344 155.003 157.708 160.461 16
+3.261 166.111 169.01 171.96 \
+ 174.961 178.015 181.122 184.283 187.499 190.771 194.101 197.489 200.936 204.443 208.011 211.641 215.335 219.093 222.917 226.808 230.767 234.794 238.892 243.062 247.304 251.62 50.8975 51.7859 52.6897 53.6093 5
+4.545 55.497 56.4656 57.4511 \
+ 58.4538 59.474 60.512 61.5681 62.6427 63.736 64.8484 65.9803 67.1318 68.3035 69.4956 70.7086 71.9426 73.1983 74.4758 75.7757 77.0982 78.4438 79.8129 81.2059 82.6233 84.0653 85.5325 87.0253 88.5442 90.0896 9
+1.662 93.2618 94.8895 96.5456 \
+ 98.2307 99.9451 101.689 103.464 105.27 107.107 108.977 110.879 112.814 114.783 116.786 118.825 120.898 123.009 125.155 127.34 129.562 131.824 134.124 136.465 138.847 141.27 143.736 146.245 148.797 151.394 15
+4.036 156.725 159.46 162.243 \
+ 165.075 167.956 170.888 173.87 176.905 179.992 183.134 186.33 189.582 192.891 196.258 199.683 203.168 206.714 210.322 213.993 217.727 221.528 225.394 229.328 233.33 237.403 241.546 245.762 250.051 254.416 25
+8.856 263.374 267.971 53.9876 \
+ 54.9299 55.8886 56.864 57.8565 58.8663 59.8937 60.939 62.0026 63.0848 64.1858 65.306 66.4458 67.6055 68.7855 69.986 71.2075 72.4503 73.7148 75.0014 76.3104 77.6423 78.9974 80.3761 81.779 83.2063 84.6585 86
+.1361 87.6394 89.169 90.7253 \
+ 92.3088 93.9199 95.5591 97.2269 98.9238 100.65 102.407 104.194 106.013 107.863 109.746 111.661 113.61 115.593 117.61 119.663 121.752 123.877 126.039 128.238 130.477 132.754 135.071 137.428 139.827 142.267 1
+44.75 147.277 149.847 152.462 \
+ 155.123 157.831 160.586 163.388 166.24 169.141 172.093 175.097 178.153 181.262 184.426 187.645 190.92 194.252 197.642 201.092 204.602 208.173 211.806 215.503 219.264 223.091 226.984 230.946 234.977 239.078 24
+3.251 247.496 251.816 256.211 \
+ 260.683 265.232 269.861 274.571 279.364 284.239 57.0612 58.0571 59.0704 60.1014 61.1504 62.2176 63.3035 64.4084 65.5325 66.6763 67.84 69.024 70.2287 71.4545 72.7016 73.9705 75.2615 76.575 77.9115 79.2713 80
+.6549 82.0626 83.4948 84.9521 \
+ 86.4348 87.9434 89.4783 91.04 92.6289 94.2456 95.8905 97.5641 99.2669 100.999 102.762 104.556 106.381 108.237 110.126 112.048 114.004 115.994 118.018 120.078 122.174 124.306 126.476 128.683 130.929 133.214 13
+5.539 137.905 140.312 142.761 \
+ 145.252 147.787 150.367 152.991 155.661 158.378 161.142 163.955 166.816 169.728 172.69 175.704 178.771 181.891 185.066 188.296 191.582 194.926 198.328 201.789 205.311 208.895 212.54 216.25 220.024 223.864 22
+7.772 231.747 235.792 239.907 \
+ 244.094 248.354 252.689 257.099 261.587 266.152 270.797 275.524 280.332 285.225 290.203 295.268 300.422 60.1175 61.1667 62.2343 63.3205 64.4256 65.5501 66.6941 67.8581 69.0425 70.2475 71.4736 72.721 73.9902 75
+.2816 76.5955 77.9324 79.2925 \
+ 80.6765 82.0845 83.5172 84.9748 86.4579 87.9669 89.5022 91.0643 92.6537 94.2708 95.9161 97.5902 99.2935 101.026 102.79 104.584 106.409 108.266 110.156 112.078 114.035 116.025 118.05 120.11 122.207 124.339 1
+26.51 128.718 130.964 133.25 \
+ 135.576 137.942 140.349 142.799 145.291 147.827 150.407 153.032 155.703 158.421 161.186 163.999 166.861 169.773 172.736 175.751 178.819 181.94 185.115 188.346 191.633 194.978 198.381 201.843 205.366 208.95 21
+2.597 216.308 220.083 223.924 \
+ 227.833 231.809 235.855 239.971 244.16 248.421 252.757 257.168 261.657 266.223 270.87 275.597 280.407 285.301 290.281 295.347 300.502 305.747 311.083 316.512 63.1554 64.2577 65.3792 66.5202 67.6812 68.8625 70
+.0644 71.2872 72.5314 73.7973 \
+ 75.0854 76.3958 77.7292 79.0858 80.4661 81.8705 83.2994 84.7533 86.2325 87.7376 89.2689 90.8269 92.4121 94.025 95.6661 97.3358 99.0346 100.763 102.522 104.311 106.132 107.984 109.869 111.786 113.737 115.722 11
+7.742 119.797 121.888 124.015 \
+ 126.18 128.382 130.623 132.902 135.222 137.582 139.983 142.427 144.912 147.442 150.015 152.633 155.297 158.008 160.765 163.571 166.426 169.331 172.286 175.293 178.353 181.465 184.633 187.855 191.134 194.47 19
+7.864 201.317 204.831 208.406 \
+ 212.043 215.744 219.509 223.341 227.239 231.205 235.24 239.346 243.523 247.773 252.098 256.498 260.974 265.529 270.164 274.879 279.676 284.558 289.524 294.577 299.719 304.95 310.272 315.687 321.197 326.803 33
+2.507 66.1741 67.329 68.5041 \
+ 69.6998 70.9163 72.154 73.4133 74.6946 75.9983 77.3247 78.6743 80.0474 81.4445 82.866 84.3122 85.7838 87.281 88.8043 90.3542 91.9312 93.5357 95.1682 96.8292 98.5192 100.239 101.988 103.768 105.579 107.422 10
+9.297 111.204 113.145 115.12 \
+ 117.129 119.174 121.254 123.37 125.523 127.714 129.943 132.211 134.518 136.866 139.255 141.685 144.158 146.674 149.234 151.839 154.489 157.185 159.929 162.72 165.56 168.45 171.39 174.381 177.424 180.521 18
+3.672 186.877 190.139 193.458 \
+ 196.834 200.269 203.765 207.321 210.94 214.621 218.367 222.178 226.056 230.001 234.016 238.1 242.256 246.484 250.786 255.163 259.616 264.147 268.758 273.448 278.221 283.077 288.017 293.044 298.159 303.363 30
+8.657 314.045 319.526 325.102 \
+ 330.777 336.55 342.424 348.4 69.1726 70.3799 71.6083 72.8581 74.1297 75.4235 76.7399 78.0792 79.442 80.8285 82.2392 83.6746 85.135 86.6208 88.1327 89.6709 91.2359 92.8283 94.4484 96.0969 97.7741 99.4806 10
+1.217 102.983 104.781 106.61 \
+ 108.47 110.363 112.29 114.249 116.243 118.272 120.337 122.437 124.574 126.748 128.96 131.211 133.501 135.831 138.202 140.614 143.068 145.565 148.106 150.69 153.321 155.996 158.719 161.489 164.308 167.176 17
+0.093 173.062 176.082 179.156 \
+ 182.283 185.464 188.701 191.994 195.345 198.755 202.224 205.753 209.344 212.998 216.715 220.498 224.346 228.262 232.246 236.299 240.423 244.62 248.889 253.233 257.653 262.15 266.725 271.38 276.117 280.936 28
+5.839 290.828 295.904 301.068 \
+ 306.323 311.669 317.109 322.644 328.275 334.004 339.834 345.765 351.8 357.94 364.187
+######################
diff --git a/macro/run/much_run_unpack_tsa.C b/macro/run/much_run_unpack_tsa.C
new file mode 100644
index 0000000000000000000000000000000000000000..f2800834849c685aa5aa6ab99e60ffa76f0de37a
--- /dev/null
+++ b/macro/run/much_run_unpack_tsa.C
@@ -0,0 +1,336 @@
+/* Copyright (C) 2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Jan de Cuveland, Volker Friese, Pierre-Alain Loizeau, Pascal Raisig [committer], Dominik Smith, Adrian A. Weber */
+
+/** @file run_unpack_tsa.C
+ ** @author Volker Friese <v.friese@gsi.de>
+ ** @since May 2021
+ **/
+
+
+// --- Includes needed for IDE
+#include <RtypesCore.h>
+
+#include <cstdint>
+#include <memory>
+#include <string>
+#include <vector>
+#if !defined(__CLING__)
+
+#include <FairLogger.h>
+#include <FairRootFileSink.h>
+#include <FairRunOnline.h>
+#include <Logger.h>
+
+#include <TStopwatch.h>
+#include <TSystem.h>
+#endif
+
+std::shared_ptr<CbmMuchUnpackMonitor> GetMuchMonitor(std::string treefilename, bool bDebugMode = false);
+const char* defaultSetupName = "mcbm_beam_2021_07_surveyed";
+
+void much_run_unpack_tsa(std::vector<std::string> infile = {"test.tsa"}, UInt_t runid = 0,
+ const char* setupName = defaultSetupName, std::int32_t nevents = -1,
+ std::string outpath = "data/withTof/")
+{
+
+ // ========================================================================
+ // Adjust this part according to your requirements
+
+ // --- Logger settings ----------------------------------------------------
+ TString logLevel = "INFO";
+ //TString logLevel = "DEBUG";
+ TString logVerbosity = "LOW";
+ // ------------------------------------------------------------------------
+
+ // ----- Environment --------------------------------------------------
+ TString myName = "run_unpack_tsa"; // this macro's name for screen output
+ TString srcDir = gSystem->Getenv("VMCWORKDIR"); // top source directory
+ // ------------------------------------------------------------------------
+
+
+ // ----- Output filename ----------------------------------------------
+ std::string outfilename = infile[0];
+ auto filenamepos = infile[0].find_last_of("/");
+ filenamepos++;
+ std::string filename = infile[0].substr(filenamepos);
+ if (filename.find("*") != infile[0].npos) filename = std::to_string(runid) + ".tsa";
+ if (filename.find(";") != infile[0].npos) filename = std::to_string(runid) + "_merged" + ".tsa";
+ if (outpath.empty()) { outpath = infile[0].substr(0, filenamepos); }
+ outfilename = outpath + filename;
+ outfilename.replace(outfilename.find(".tsa"), 4, ".digi.root");
+ std::cout << "-I- " << myName << ": Output file will be " << outfilename << std::endl;
+ // ------------------------------------------------------------------------
+
+
+ // ----- Performance profiling ----------------------------------------
+ // Set to true if you want some minimal performance profiling output
+ bool doPerfProfiling = true;
+ // Define if you want a special path and name for the performance profiling output file
+ std::string perfProfFileName = outpath + filename;
+ perfProfFileName.replace(perfProfFileName.find(".tsa"), 4, ".perf.root");
+ // ------------------------------------------------------------------------
+
+
+ // ----- CbmSetup -----------------------------------------------------
+ auto cbmsetup = CbmSetup::Instance();
+ cbmsetup->LoadSetup(setupName);
+ // ------------------------------------------------------------------------
+
+ // ----- UnpackerConfigs ----------------------------------------------
+
+ /*
+ // ---- STS ----
+ std::shared_ptr<CbmStsUnpackConfig> stsconfig = nullptr;
+
+ stsconfig = std::make_shared<CbmStsUnpackConfig>(std::string(setupName), runid);
+ if (stsconfig) {
+ // stsconfig->SetDebugState();
+ stsconfig->SetDoWriteOutput();
+ stsconfig->SetDoWriteOptOutA("StsDigiPulser");
+ std::string parfilesbasepathSts = Form("%s/macro/beamtime/mcbm2021/", srcDir.Data());
+ stsconfig->SetParFilesBasePath(parfilesbasepathSts);
+ /// Enable duplicates rejection, Ignores the ADC for duplicates check
+ stsconfig->SetDuplicatesRejection(true, true);
+ /// Enable Monitor plots
+ //stsconfig->SetMonitor(GetStsMonitor(outfilename, true));
+ stsconfig->SetSystemTimeOffset(-2221); // [ns] value to be updated
+
+ stsconfig->SetMinAdcCut(1, 1);
+ stsconfig->SetMinAdcCut(2, 1);
+ stsconfig->SetMinAdcCut(3, 1);
+ stsconfig->SetMinAdcCut(4, 1);
+
+ stsconfig->MaskNoisyChannel(3, 56);
+ stsconfig->MaskNoisyChannel(3, 75);
+ stsconfig->MaskNoisyChannel(3, 79);
+ stsconfig->MaskNoisyChannel(3, 85);
+ stsconfig->MaskNoisyChannel(7, 123);
+ stsconfig->MaskNoisyChannel(7, 124);
+ stsconfig->MaskNoisyChannel(7, 125);
+ stsconfig->MaskNoisyChannel(7, 158);
+ stsconfig->MaskNoisyChannel(7, 159);
+ stsconfig->MaskNoisyChannel(7, 162);
+ stsconfig->MaskNoisyChannel(7, 715);
+ stsconfig->MaskNoisyChannel(9, 709);
+ stsconfig->MaskNoisyChannel(12, 119);
+
+ // Time Walk correction
+ std::map<uint32_t, CbmStsParModule> walkMap;
+ auto parAsic = new CbmStsParAsic(128, 31, 31., 1., 5., 800., 1000., 3.9789e-3);
+
+ // Module params: number of channels, number of channels per ASIC
+ auto parMod = new CbmStsParModule(2048, 128);
+
+ // default
+ double p0 = 0, p1 = 0, p2 = 0, p3 = 0;
+ parAsic->SetWalkCoef({p0, p1, p2, p3});
+ parMod->SetAllAsics(*parAsic);
+
+ walkMap[0x10107C02] = CbmStsParModule(*parMod); // Make a copy for storage
+ walkMap[0x101FFC02] = CbmStsParModule(*parMod); // Make a copy for storage
+
+ /// To be replaced by a storage in a new parameter class later
+ int sensor, asic;
+ std::ifstream asicTimeWalk_par(Form("%s/mStsAsicTimeWalk.par", parfilesbasepathSts.data()));
+ while (asicTimeWalk_par >> std::hex >> sensor >> std::dec >> asic >> p0 >> p1 >> p2 >> p3) {
+ std::cout << Form("Setting time-walk parametersfor: module %x, ASIC %u\n", sensor, asic);
+ parAsic->SetWalkCoef({p0, p1, p2, p3});
+
+ if (walkMap.find(sensor) == walkMap.end()) { walkMap[sensor] = CbmStsParModule(*parMod); }
+ walkMap[sensor].SetAsic(asic, *parAsic);
+ }
+
+ stsconfig->SetWalkMap(walkMap);
+ }
+ // -------------
+ */
+ // ---- MUCH ----
+ std::shared_ptr<CbmMuchUnpackConfig> muchconfig = nullptr;
+
+ muchconfig = std::make_shared<CbmMuchUnpackConfig>(std::string(setupName), runid);
+ if (muchconfig) {
+ // muchconfig->SetDebugState();
+ muchconfig->SetDoWriteOutput();
+ muchconfig->SetDoWriteOptOutA("MuchDigiPulser");
+ std::string parfilesbasepathMuch = Form("%s/macro/beamtime/mcbm2022/", srcDir.Data());
+ muchconfig->SetParFilesBasePath(parfilesbasepathMuch);
+ if (2060 <= runid && runid <= 2162) {
+ /// Starting to use CRI Based MUCH setup with 2GEM and 1 RPC since 09/03/2022 Carbon run
+ muchconfig->SetParFileName("mMuchParUpto26032022.par");
+ }
+ /// Enable duplicates rejection, Ignores the ADC for duplicates check
+ muchconfig->SetDuplicatesRejection(true, true);
+ /// Enable Monitor plots
+ //muchconfig->SetMonitor(GetMuchMonitor(outfilename, true));
+ muchconfig->SetSystemTimeOffset(-2221); // [ns] value to be updated
+
+ // muchconfig->SetMinAdcCut(1, 1);
+
+ // muchconfig->MaskNoisyChannel(3, 56);
+ }
+ // -------------
+
+
+ // ---- TOF ----
+ std::shared_ptr<CbmTofUnpackConfig> tofconfig = nullptr;
+
+ tofconfig = std::make_shared<CbmTofUnpackConfig>("", runid);
+ if (tofconfig) {
+ // tofconfig->SetDebugState();
+ tofconfig->SetDoWriteOutput();
+ // tofconfig->SetDoWriteOptOutA("CbmTofErrors");
+ std::string parfilesbasepathTof = Form("%s/macro/beamtime/mcbm2021/", srcDir.Data());
+ tofconfig->SetParFilesBasePath(parfilesbasepathTof);
+ tofconfig->SetSystemTimeOffset(-1220); // [ns] value to be updated
+ if (runid <= 1659) {
+ /// Switch ON the -4 offset in epoch count (hack for Spring-Summer 2021)
+ tofconfig->SetFlagEpochCountHack2021();
+ }
+ }
+ // -------------
+
+
+ // ------------------------------------------------------------------------
+
+ // In general, the following parts need not be touched
+ // ========================================================================
+
+ // ----- Timer --------------------------------------------------------
+ TStopwatch timer;
+ timer.Start();
+ // ------------------------------------------------------------------------
+
+ // ----- CbmSourceTsArchive -------------------------------------------
+ auto source = new CbmSourceTsArchive(infile);
+ auto unpack = source->GetRecoUnpack();
+ unpack->SetDoPerfProfiling(doPerfProfiling);
+ unpack->SetOutputFilename(perfProfFileName);
+ // Enable full time sorting instead sorting per FLIM link
+ unpack->SetTimeSorting(true);
+
+ //if (stsconfig) unpack->SetUnpackConfig(stsconfig);
+ if (muchconfig) unpack->SetUnpackConfig(muchconfig);
+ if (tofconfig) unpack->SetUnpackConfig(tofconfig);
+ // ------------------------------------------------------------------------
+
+ // ----- FairRunAna ---------------------------------------------------
+ auto run = new FairRunOnline(source);
+ auto sink = new FairRootFileSink(outfilename.data());
+ run->SetSink(sink);
+ auto eventheader = new CbmTsEventHeader();
+ run->SetRunId(runid);
+ run->SetEventHeader(eventheader);
+ // ------------------------------------------------------------------------
+
+
+ // ----- Logger settings ----------------------------------------------
+ FairLogger::GetLogger()->SetLogScreenLevel(logLevel.Data());
+ FairLogger::GetLogger()->SetLogVerbosityLevel(logVerbosity.Data());
+ // ------------------------------------------------------------------------
+
+
+ // ----- Run initialisation -------------------------------------------
+ std::cout << std::endl;
+ std::cout << "-I- " << myName << ": Initialise run" << std::endl;
+ run->Init();
+ // ------------------------------------------------------------------------
+
+
+ // ----- Start run ----------------------------------------------------
+ std::cout << std::endl << std::endl;
+ std::cout << "-I- " << myName << ": Starting run" << std::endl;
+ if (nevents < 0) run->Run(-1, 0);
+ else
+ run->Run(0, nevents);
+ // ------------------------------------------------------------------------
+
+
+ // ----- Finish -------------------------------------------------------
+ timer.Stop();
+ std::cout << "Macro finished successfully." << std::endl;
+ std::cout << "After CpuTime = " << timer.CpuTime() << " s RealTime = " << timer.RealTime() << " s." << std::endl;
+ // ------------------------------------------------------------------------
+
+} // End of main macro function
+
+
+/**
+ * @brief Get the Sts Monitor. Extra function to keep default macro part more silent.
+ * @return std::shared_ptr<CbmStsUnpackMonitor>
+*/
+std::shared_ptr<CbmStsUnpackMonitor> GetStsMonitor(std::string treefilename, bool bDebugMode = false)
+{
+ // ----- Output filename and path -------------------------------------
+ std::string outpath = "";
+ std::string filename = "";
+ auto filenamepos = treefilename.find_last_of("/");
+ if (filenamepos != treefilename.npos) {
+ outpath = treefilename.substr(0, filenamepos);
+ filename = treefilename.substr(filenamepos++);
+ }
+ if (outpath.empty()) outpath = gSystem->GetWorkingDirectory();
+ //std::string mydir = "/qa";
+ //outpath += mydir;
+
+ auto currentdir = gSystem->GetWorkingDirectory();
+
+ if (!gSystem->cd(outpath.data())) gSystem->MakeDirectory(outpath.data());
+ else
+ gSystem->cd(currentdir.data());
+
+ std::string outfilename = outpath + filename;
+ auto filetypepos = outfilename.find(".digi.root");
+ if (filetypepos != outfilename.npos) outfilename.replace(filetypepos, 10, ".mon.sts.root");
+ else
+ outfilename += ".mon.sts.root";
+ // ------------------------------------------------------------------------
+
+ auto monitor = std::make_shared<CbmStsUnpackMonitor>();
+ monitor->SetHistoFileName(outfilename);
+ monitor->SetDebugMode(bDebugMode);
+ return monitor;
+}
+
+/**
+ * @brief Get the Much Monitor. Extra function to keep default macro part more silent.
+ * @return std::shared_ptr<CbmMuchUnpackMonitor>
+*/
+std::shared_ptr<CbmMuchUnpackMonitor> GetMuchMonitor(std::string treefilename, bool bDebugMode = false)
+{
+ // ----- Output filename and path -------------------------------------
+ std::string outpath = "";
+ std::string filename = "";
+ auto filenamepos = treefilename.find_last_of("/");
+ if (filenamepos != treefilename.npos) {
+ outpath = treefilename.substr(0, filenamepos);
+ filename = treefilename.substr(filenamepos++);
+ }
+ if (outpath.empty()) outpath = gSystem->GetWorkingDirectory();
+
+ auto currentdir = gSystem->GetWorkingDirectory();
+
+ if (!gSystem->cd(outpath.data())) gSystem->MakeDirectory(outpath.data());
+ else
+ gSystem->cd(currentdir.data());
+
+ std::string outfilename = outpath + filename;
+ auto filetypepos = outfilename.find(".digi.root");
+ if (filetypepos != outfilename.npos) outfilename.replace(filetypepos, 10, ".mon.much.root");
+ else
+ outfilename += ".mon.much.root";
+ // ------------------------------------------------------------------------
+
+ auto monitor = std::make_shared<CbmMuchUnpackMonitor>();
+ monitor->SetHistoFileName(outfilename);
+ monitor->SetDebugMode(bDebugMode);
+ return monitor;
+}
+
+
+void much_run_unpack_tsa(std::string infile = "test.tsa", UInt_t runid = 0, const char* setupName = defaultSetupName,
+ std::int32_t nevents = -1, std::string outpath = "data/")
+{
+ std::vector<std::string> vInFile = {infile};
+ return much_run_unpack_tsa(vInFile, runid, setupName, nevents, outpath);
+}
diff --git a/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx b/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx
index 02a363a557433ee1a10641e944c5139c43ad856b..e0b34cfa1b6efa6b28dbe77811d71e3e5fcbb06e 100644
--- a/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx
+++ b/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx
@@ -466,11 +466,6 @@ void CbmMuchUnpackAlgo::processHitInfo(const stsxyter::Message& mess)
fuSkippMessNoValidElink++;
return;
}
- //uFebIdx += ( fuCurrDpbIdx * fNrCrobPerDpb + ____CrobIndexCalculationIfNeeded___() ) * fuNbFebsPerCrob;
- //uncomment below line if kuNbCrobsPerDpb > 1
- //uFebIdx += (fuCurrDpbIdx * fNrCrobPerDpb) * fNrFebsPerCrob; //no Crob index calculation for now
-
-
const uint16_t usChan = mess.GetHitChannel();
const uint16_t usRawAdc = mess.GetHitAdc();
const uint16_t usRawTs = mess.GetHitTimeBinning();
@@ -479,6 +474,11 @@ void CbmMuchUnpackAlgo::processHitInfo(const stsxyter::Message& mess)
// Compute the Full time stamp
const uint64_t ulHitTime = getFullTimeStamp(usRawTs);
+ //uFebIdx += ( fuCurrDpbIdx * fNrCrobPerDpb + ____CrobIndexCalculationIfNeeded___() ) * fuNbFebsPerCrob;
+ //uncomment below line if kuNbCrobsPerDpb > 1
+ //uFebIdx += (fuCurrDpbIdx * fNrCrobPerDpb) * fNrFebsPerCrob; //no Crob index calculation for now
+
+
if (fMonitor) fMonitor->CountRawHit(uFebIdx, uChanInFeb);
/// Store hit !!!
@@ -534,6 +534,12 @@ void CbmMuchUnpackAlgo::processHitInfo(const stsxyter::Message& mess)
<< uFebIdx << " Channnel Id " << uChanInFeb;
if (address) {
/// Prevent Clang from onlining
+ // Storing Pulser FEB Data in another output stream (pulser feb gives -7 )
+ if (-7 == uFebIdx) {
+ //presently for Pulser Output CbmMuchDigi address is 0
+ fOptOutAVec->emplace_back(CbmMuchDigi(address, usRawAdc, ulTimeInNs));
+ return;
+ } // if (-7 == uFebIdx)
fOutputVec.emplace_back(CbmMuchDigi(address, usRawAdc, ulTimeInNs));
}
else {
@@ -541,7 +547,7 @@ void CbmMuchUnpackAlgo::processHitInfo(const stsxyter::Message& mess)
fuSkippMessNoValidAddress++;
}
/// Catch the pulser digis and either save them to their own output or drop them
- ///Commented below 6 lines as we do not use FebPulser Vector
+ //Commented below 6 lines as we do not use FebPulser Vector
//if (fvbFebPulser[uFebIdx]) {
// if (fOptOutAVec) {
// fOptOutAVec->emplace_back(CbmMuchDigi(address, usRawAdc, ulTimeInNs));
@@ -580,13 +586,18 @@ void CbmMuchUnpackAlgo::processHitInfo(const stsxyter::Message& mess)
//---------------CreateMuchDigi -------
//To do address need to be checked carefully
-uint32_t CbmMuchUnpackAlgo::CreateMuchAddress(uint32_t dpbidx, uint32_t iFebId, uint32_t usChan)
+uint32_t CbmMuchUnpackAlgo::CreateMuchAddress(uint32_t dpbidx, int32_t iFebId, uint32_t usChan)
{
// For generating Station number (GEM1 station = 0, GEM2 station = 1 and RPC station = 2)
/// Below FebID is according to FEB Position in Module GEM A or Module GEM B (Carefully write MUCH Par file)
int32_t station = -1;
int32_t layer = -1;
- if (dpbidx == 0 || dpbidx == 1 || dpbidx == 2) //First 3 DPBs are for GEM-1
+ if (iFebId == -7) //Pulser FEB
+ {
+ station = 6; // for Pulser
+ layer = 0; //
+ }
+ else if (dpbidx == 0 || dpbidx == 1 || dpbidx == 2) //First 3 DPBs are for GEM-1
{
station = 0; // for mCBM setup
layer = 0; // Station 0 for GEM-A and station 1 for Module GEM-B
diff --git a/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx.org b/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx.org
deleted file mode 100644
index 194e44a2f31302c22eabfbed21480eae314eeacd..0000000000000000000000000000000000000000
--- a/reco/detectors/much/unpack/CbmMuchUnpackAlgo.cxx.org
+++ /dev/null
@@ -1,817 +0,0 @@
-/* Copyright (C) 2022 Fair GmbH, Darmstadt
- SPDX-License-Identifier: GPL-3.0-only
- Authors: Pierre-Alain Loizeau [committer] */
-
-#include "CbmMuchUnpackAlgo.h"
-
-#include "CbmMuchDigi.h"
-
-#include <FairParGenericSet.h>
-#include <FairTask.h>
-#include <Logger.h>
-
-#include <Rtypes.h>
-#include <RtypesCore.h>
-
-#include <cstdint>
-#include <iomanip>
-
-CbmMuchUnpackAlgo::CbmMuchUnpackAlgo() : CbmRecoUnpackAlgo("CbmMuchUnpackAlgo") {}
-
-CbmMuchUnpackAlgo::CbmMuchUnpackAlgo(std::string name) : CbmRecoUnpackAlgo(name) {}
-
-CbmMuchUnpackAlgo::~CbmMuchUnpackAlgo() {}
-
-// ---- GetParContainerRequest ----
-std::vector<std::pair<std::string, std::shared_ptr<FairParGenericSet>>>*
- CbmMuchUnpackAlgo::GetParContainerRequest(std::string /*geoTag*/, std::uint32_t /*runId*/)
-{
- // Basepath for default Much parameter sets (those connected to a geoTag)
- std::string basepath = Form("%s", fParFilesBasePath.data());
- std::string temppath = "";
-
- // // Get parameter container
- temppath = basepath + "mMuchPar.par";
- fParContVec.emplace_back(std::make_pair(temppath, std::make_shared<CbmMuchUnpackPar>()));
-
- return &fParContVec;
-}
-
-// ---- getAsicIndex ----
-uint32_t CbmMuchUnpackAlgo::getAsicIndex(uint32_t dpbidx, uint32_t crobidx, uint16_t elinkidx)
-{
-
- uint32_t asicidx = 0;
- const int32_t uFebIdx = fElinkIdxToFebIdxVec.at(elinkidx);
- const uint32_t febtype = fviFebType[dpbidx][crobidx][uFebIdx];
- // Feb type a
- if (febtype == 0) asicidx = fElinkIdxToAsicIdxVec.at(elinkidx).first;
- // Feb type b
- if (febtype == 1) asicidx = fElinkIdxToAsicIdxVec.at(elinkidx).second;
- // else would be inactive feb, this was not handled in the previous implementation, this I expect it should not happen
-
- const uint32_t uAsicIdx = (dpbidx * fNrCrobPerDpb + crobidx) * fNrAsicsPerCrob + asicidx;
- return uAsicIdx;
-}
-
-// ---- getFullTimeStamp ----
-uint64_t CbmMuchUnpackAlgo::getFullTimeStamp(const uint16_t usRawTs)
-{
- // Use TS w/o overlap bits as they will anyway come from the TS_MSB
- const uint64_t ulTime =
- usRawTs + fulTsMsbIndexInTs[fuCurrDpbIdx] * static_cast<uint64_t>(stsxyter::kuHitNbTsBinsBinning);
- /*
- + static_cast<uint64_t>(stsxyter::kuHitNbTsBinsBinning) * static_cast<uint64_t>(fvulCurrentTsMsb[fuCurrDpbIdx])
- + static_cast<uint64_t>(stsxyter::kulTsCycleNbBinsBinning)
- * static_cast<uint64_t>(fvuCurrentTsMsbCycle[fuCurrDpbIdx]);
-*/
- return ulTime;
-}
-
-// ---- init
-Bool_t CbmMuchUnpackAlgo::init() { return kTRUE; }
-
-// ---- initDpbIdIndexMap ----
-void CbmMuchUnpackAlgo::initDpbIdIndexMap(CbmMuchUnpackPar* parset)
-{
- fDpbIdIndexMap.clear();
- for (uint32_t uDpb = 0; uDpb < parset->GetNrOfDpbs(); ++uDpb) {
- fDpbIdIndexMap[parset->GetDpbId(uDpb)] = uDpb;
- LOG(debug) << "Eq. ID for DPB #" << std::setw(2) << uDpb << " = 0x" << std::setw(4) << std::hex
- << parset->GetDpbId(uDpb) << std::dec << " => " << fDpbIdIndexMap[parset->GetDpbId(uDpb)];
- }
-}
-
-// ---- initParSet(FairParGenericSet* parset) ----
-Bool_t CbmMuchUnpackAlgo::initParSet(FairParGenericSet* parset)
-{
- LOG(info) << fName << "::initParSet - for container " << parset->ClassName();
- if (parset->IsA() == CbmMuchUnpackPar::Class()) return initParSet(static_cast<CbmMuchUnpackPar*>(parset));
-
- // If we do not know the derived ParSet class we return false
- LOG(error)
- << fName << "::initParSet - for container " << parset->ClassName()
- << " failed, since CbmMuchUnpackAlgo::initParSet() does not know the derived ParSet and what to do with it!";
- return kFALSE;
-}
-
-// ---- initParSet(CbmMuchUnpackPar* parset) ----
-Bool_t CbmMuchUnpackAlgo::initParSet(CbmMuchUnpackPar* parset)
-{
- LOG(debug) << fName << "::initParSetAsic - ";
-
- //Type of each module: 0 for connectors on the right, 1 for connectors on the left
- std::vector<int32_t> viModuleType;
-
- // STS address for the first strip of each module
- std::vector<int32_t> viModAddress;
-
- // Idx of the STS module for each FEB, [ NbDpb ][ NbCrobPerDpb ][ NbFebsPerCrob ], -1 if inactive
- std::vector<std::vector<std::vector<int32_t>>> viFebModuleIdx;
-
- // Array to hold the active flag for all CROBs, [ NbDpb ][ NbCrobPerDpb ]
- std::vector<std::vector<bool>> vbCrobActiveFlag;
-
- //STS module side for each FEB, [ NbDpb ][ NbCrobPerDpb ][ NbFebsPerCrob ], 0 = P, 1 = N, -1 if inactive
- std::vector<std::vector<std::vector<int32_t>>> viFebModuleSide;
-/*
- // Total number of STS modules in the setup
- const uint32_t uNbModules = parset->GetNbOfModules();
- LOG(debug) << "Nr. of STS Modules: " << uNbModules;
-
- // Total number of STS DPBs in system
- const uint32_t uNbOfDpbs = parset->GetNrOfDpbs();
- LOG(debug) << "Nr. of STS DPBs: " << uNbOfDpbs;
-
- // Get Nr of Febs
- fuNbFebs = parset->GetNrOfFebs();
- LOG(debug) << "Nr. of FEBs: " << fuNbFebs;
-
- // Get Nr of eLinks per CROB
- fNrElinksPerCrob = parset->GetNbElinkPerCrob();
-
- // Get Nr of ASICs per CROB
- fNrAsicsPerCrob = parset->GetNbAsicsPerCrob();
-
- // Get Nr of CROBs per DPB
- fNrCrobPerDpb = parset->GetNbCrobsPerDpb();
-
- // Get Number of ASICs per FEB
- fNrAsicsPerFeb = parset->GetNbAsicsPerFeb();
-
- // Get Number of Channels per Asic
- fNrChsPerAsic = parset->GetNbChanPerAsic();
-
- // Get Number of Channels per FEB
- fNrChsPerFeb = parset->GetNbChanPerFeb();
-
- // Get Number of FEBs per CROB
- fNrFebsPerCrob = parset->GetNbFebsPerCrob();
-
- for (size_t ielink = 0; ielink < fNrElinksPerCrob; ++ielink) {
- fElinkIdxToFebIdxVec.emplace_back(parset->ElinkIdxToFebIdx(ielink));
- fElinkIdxToAsicIdxVec.emplace_back(
- std::make_pair(parset->ElinkIdxToAsicIdxFebA(ielink), parset->ElinkIdxToAsicIdxFebB(ielink)));
- }
-
- // Get Nr of Asics
- const uint32_t uNbStsXyters = parset->GetNrOfAsics();
- LOG(debug) << "Nr. of StsXyter ASICs: " << uNbStsXyters;
-
- //Initialize temporary "per Feb" fields
- initTempVectors(parset, &viModuleType, &viModAddress, &viFebModuleIdx, &vbCrobActiveFlag, &viFebModuleSide);
-
- // Initialize ADC cuts for FEBs
- fvbFebAdcCut.resize(fuNbFebs, fdAdcCut);
- for (auto cut : fdAdcCut_perFeb) {
- fvbFebAdcCut[cut.first] = cut.second;
- }
-
- // Read dpb index map from parameter container
- initDpbIdIndexMap(parset);
-
- if (fvdTimeOffsetNsAsics.size() < uNbStsXyters) { fvdTimeOffsetNsAsics.resize(uNbStsXyters, 0.0); }
-
- //Initialize class-wide "per Feb" fields
- fviFebType.resize(uNbOfDpbs);
-
- for (uint32_t uDpb = 0; uDpb < uNbOfDpbs; ++uDpb) {
- fviFebType[uDpb].resize(fNrCrobPerDpb);
- for (uint32_t uCrobIdx = 0; uCrobIdx < fNrCrobPerDpb; ++uCrobIdx) {
- fviFebType[uDpb][uCrobIdx].resize(parset->GetNbFebsPerCrob(), -1);
- for (uint32_t uFebIdx = 0; uFebIdx < parset->GetNbFebsPerCrob(); ++uFebIdx) {
- fvbFebPulser.push_back(parset->IsFebPulser(uDpb, uCrobIdx, uFebIdx));
- fvdFebAdcGain.push_back(parset->GetFebAdcGain(uDpb, uCrobIdx, uFebIdx));
- fvdFebAdcOffs.push_back(parset->GetFebAdcOffset(uDpb, uCrobIdx, uFebIdx));
-
- if (0 <= viFebModuleIdx[uDpb][uCrobIdx][uFebIdx]
- && static_cast<uint32_t>(viFebModuleIdx[uDpb][uCrobIdx][uFebIdx]) < uNbModules
- && 0 <= viFebModuleSide[uDpb][uCrobIdx][uFebIdx] && viFebModuleSide[uDpb][uCrobIdx][uFebIdx] < 2) {
- switch (viModuleType[viFebModuleIdx[uDpb][uCrobIdx][uFebIdx]]) {
- case 0: // FEB-8-1 with ZIF connector on the right
- {
- // P side (0) has type A (0)
- // N side (1) has type B (1)
- fviFebType[uDpb][uCrobIdx][uFebIdx] = viFebModuleSide[uDpb][uCrobIdx][uFebIdx];
-
- ///! FIXME: 1) Geometry is using front/back while we are using P/N !!!!
- ///! => Assuming that front facing modules have connectors on right side
- ///! +> Volker warns that the front side should be electrons one so N
- ///! 2) No accessor/setter to change only the side field of an STS address
- ///! => hardcode the shift
- ///! +> The bit is unused in the current scheme: the side is encoded in the Digi channel
- fviFebAddress.push_back(viModAddress[viFebModuleIdx[uDpb][uCrobIdx][uFebIdx]]);
- // + (viFebModuleSide[uDpb][uCrobIdx][uFebIdx] << 25));
- fviFebSide.push_back(viFebModuleSide[uDpb][uCrobIdx][uFebIdx]);
- break;
- } // case 0: // FEB-8-1 with ZIF connector on the right
- case 1: // FEB-8-1 with ZIF connector on the left
- {
- // P side (0) has type B (1)
- // N side (1) has type A (0)
- fviFebType[uDpb][uCrobIdx][uFebIdx] = !(viFebModuleSide[uDpb][uCrobIdx][uFebIdx]);
-
- ///! FIXME: 1) Geometry is using front/back while we are using P/N !!!!
- ///! => Assuming that front facing modules have connectors on right side
- ///! +> Volker warns that the front side should be electrons one so N
- ///! 2) No accessor/setter to change only the side field of an STS address
- ///! => hardcode the shift
- ///! +> The bit is unused in the current scheme: the side is encoded in the Digi channel
- fviFebAddress.push_back(viModAddress[viFebModuleIdx[uDpb][uCrobIdx][uFebIdx]]);
- // + ((!viFebModuleSide[uDpb][uCrobIdx][uFebIdx]) << 25));
- fviFebSide.push_back(viFebModuleSide[uDpb][uCrobIdx][uFebIdx]);
- break;
- } // case 1: // FEB-8-1 with ZIF connector on the left
- default:
- LOG(fatal) << Form("Bad module type for DPB #%02u CROB #%u FEB %02u: %d", uDpb, uCrobIdx, uFebIdx,
- viModuleType[viFebModuleIdx[uDpb][uCrobIdx][uFebIdx]]);
- break;
- }
- } // FEB active and module index OK
- else if (-1 == viFebModuleIdx[uDpb][uCrobIdx][uFebIdx] || -1 == viFebModuleSide[uDpb][uCrobIdx][uFebIdx]) {
- fviFebAddress.push_back(0);
- fviFebSide.push_back(-1);
- } // Module index or type is set to inactive
- else {
- LOG(fatal) << Form("Bad module Index and/or Side for DPB #%02u CROB "
- "#%u FEB %02u: %d %d",
- uDpb, uCrobIdx, uFebIdx, viFebModuleIdx[uDpb][uCrobIdx][uFebIdx],
- viFebModuleSide[uDpb][uCrobIdx][uFebIdx]);
- } // Bad module index or type for this FEB
- }
- }
- }
-*/
- printActiveCrobs(parset, vbCrobActiveFlag);
- printAddressMaps(parset, viFebModuleIdx, viFebModuleSide);
-
- LOG(debug) << "Unpacking data in bin sorter FW mode";
- initInternalStatus(parset);
-
- return kTRUE;
-}
-
-// ---- initTempVectors ----
-void CbmMuchUnpackAlgo::initTempVectors(CbmMuchUnpackPar* parset, std::vector<int32_t>* viModuleType,
- std::vector<int32_t>* viModAddress,
- std::vector<std::vector<std::vector<int32_t>>>* viFebModuleIdx,
- std::vector<std::vector<bool>>* vbCrobActiveFlag,
- std::vector<std::vector<std::vector<int32_t>>>* viFebModuleSide)
-{
-/*
- const uint32_t uNbModules = parset->GetNbOfModules();
- const uint32_t uNbOfDpbs = parset->GetNrOfDpbs();
-
- viModuleType->resize(uNbModules);
- viModAddress->resize(uNbModules);
- for (uint32_t uModIdx = 0; uModIdx < uNbModules; ++uModIdx) {
- (*viModuleType)[uModIdx] = parset->GetModuleType(uModIdx);
- (*viModAddress)[uModIdx] = parset->GetModuleAddress(uModIdx);
- LOG(debug) << "Module #" << std::setw(2) << uModIdx << " Type " << std::setw(4) << (*viModuleType)[uModIdx]
- << " Address 0x" << std::setw(8) << std::hex << (*viModAddress)[uModIdx] << std::dec;
- }
- vbCrobActiveFlag->resize(uNbOfDpbs);
- viFebModuleIdx->resize(uNbOfDpbs);
- viFebModuleSide->resize(uNbOfDpbs);
-
- for (uint32_t uDpb = 0; uDpb < uNbOfDpbs; ++uDpb) {
- (*vbCrobActiveFlag)[uDpb].resize(fNrCrobPerDpb);
- (*viFebModuleIdx)[uDpb].resize(fNrCrobPerDpb);
- (*viFebModuleSide)[uDpb].resize(fNrCrobPerDpb);
- for (uint32_t uCrobIdx = 0; uCrobIdx < fNrCrobPerDpb; ++uCrobIdx) {
- (*vbCrobActiveFlag)[uDpb][uCrobIdx] = parset->IsCrobActive(uDpb, uCrobIdx);
- (*viFebModuleIdx)[uDpb][uCrobIdx].resize(parset->GetNbFebsPerCrob());
- (*viFebModuleSide)[uDpb][uCrobIdx].resize(parset->GetNbFebsPerCrob());
- for (uint32_t uFebIdx = 0; uFebIdx < parset->GetNbFebsPerCrob(); ++uFebIdx) {
- (*viFebModuleIdx)[uDpb][uCrobIdx][uFebIdx] = parset->GetFebModuleIdx(uDpb, uCrobIdx, uFebIdx);
- (*viFebModuleSide)[uDpb][uCrobIdx][uFebIdx] = parset->GetFebModuleSide(uDpb, uCrobIdx, uFebIdx);
- }
- }
- }
-*/
-}
-
-// ---- initInternalStatus ----
-void CbmMuchUnpackAlgo::initInternalStatus(CbmMuchUnpackPar* parset)
-{
-/*
- const uint32_t uNbOfDpbs = parset->GetNrOfDpbs();
- const uint32_t uNbStsXyters = parset->GetNrOfAsics();
-
- fvulCurrentTsMsb.resize(uNbOfDpbs);
- fvuCurrentTsMsbCycle.resize(uNbOfDpbs);
- fulTsMsbIndexInTs.resize(uNbOfDpbs);
- for (uint32_t uDpb = 0; uDpb < uNbOfDpbs; ++uDpb) {
- fvulCurrentTsMsb[uDpb] = 0;
- fvuCurrentTsMsbCycle[uDpb] = 0;
- fulTsMsbIndexInTs[uDpb] = 0;
- }
-
- fvvusLastTsChan.resize(uNbStsXyters);
- fvvusLastAdcChan.resize(uNbStsXyters);
- fvvulLastTsMsbChan.resize(uNbStsXyters);
- for (uint32_t uAsicIdx = 0; uAsicIdx < uNbStsXyters; ++uAsicIdx) {
- fvvusLastTsChan[uAsicIdx].resize(fNrChsPerAsic, 0);
- fvvusLastAdcChan[uAsicIdx].resize(fNrChsPerAsic, 0);
- fvvulLastTsMsbChan[uAsicIdx].resize(fNrChsPerAsic, 0);
- }
-*/
-}
-
-// ---- loopMsMessages ----
-void CbmMuchUnpackAlgo::loopMsMessages(const uint8_t* msContent, const uint32_t uSize, const size_t uMsIdx)
-{
- // If not integer number of message in input buffer, print warning/error
- if (0 != (uSize % sizeof(stsxyter::Message))) {
- LOG(error) << "The input microslice buffer does NOT "
- << "contain only complete sDPB messages!";
- }
- // Compute the number of complete messages in the input microslice buffer
- const uint32_t uNbMessages = (uSize - (uSize % sizeof(stsxyter::Message))) / sizeof(stsxyter::Message);
-
- // Prepare variables for the loop on contents
- const stsxyter::Message* pMess = reinterpret_cast<const stsxyter::Message*>(msContent);
-
- for (uint32_t uIdx = 0; uIdx < uNbMessages; uIdx++) {
- /// Get message type
- const stsxyter::MessType typeMess = pMess[uIdx].GetMessType();
-
- LOG(debug2) << " Msg Idx " << std::setw(6) << uIdx << " Type " << stsxyter::Message::PrintMessType(typeMess);
-
- if (fMonitor)
- if (fMonitor->GetDebugMode()) { fMonitor->ProcessDebugInfo(pMess[uIdx], fuCurrDpbIdx); }
- switch (typeMess) {
- case stsxyter::MessType::Hit: {
- processHitInfo(pMess[uIdx]);
- break;
- }
- case stsxyter::MessType::TsMsb: {
- processTsMsbInfo(pMess[uIdx], uIdx, uMsIdx);
- break;
- }
- case stsxyter::MessType::Epoch: {
- processEpochInfo(pMess[uIdx]);
- if (0 < uIdx) {
- LOG(warning) << "CbmMuchUnpackAlgo::loopMsMessages => "
- << "EPOCH message at unexpected position in MS: message " << uIdx << " VS message 0 expected!";
- }
- break;
- }
- case stsxyter::MessType::Status: {
- processStatusInfo(pMess[uIdx], uIdx);
- break;
- }
- case stsxyter::MessType::Empty: {
- break;
- }
- case stsxyter::MessType::EndOfMs: {
- processErrorInfo(pMess[uIdx]);
- break;
- }
- case stsxyter::MessType::Dummy: {
- break;
- }
- default: {
- LOG(fatal) << "CbmMuchUnpackAlgo::loopMsMessages => "
- << "Unknown message type, should never happen, stopping "
- "here! Type found was: "
- << static_cast<int>(typeMess);
- }
- }
- }
-}
-
-// ---- MaskNoisyChannel ----
-void CbmMuchUnpackAlgo::MaskNoisyChannel(const uint32_t uFeb, const uint32_t uChan, const bool bMasked)
-{
- if (false == fbUseChannelMask) {
- fbUseChannelMask = true;
- fvvbMaskedChannels.resize(fuNbFebs);
- for (uint32_t uFebIdx = 0; uFebIdx < fuNbFebs; ++uFebIdx) {
- fvvbMaskedChannels[uFebIdx].resize(fNrChsPerFeb, false);
- }
- }
- if (uFeb < fuNbFebs && uChan < fNrChsPerFeb) fvvbMaskedChannels[uFeb][uChan] = bMasked;
- else
- LOG(fatal) << "CbmMuchUnpackAlgo::MaskNoisyChannel => Invalid FEB "
- "and/or CHAN index:"
- << Form(" %u vs %u and %u vs %u", uFeb, fuNbFebs, uChan, fNrChsPerFeb);
-}
-
-
-// ---- printActiveCrobs ----
-void CbmMuchUnpackAlgo::printActiveCrobs(CbmMuchUnpackPar* parset,
- const std::vector<std::vector<bool>>& vbCrobActiveFlag)
-{
- for (uint32_t uDpb = 0; uDpb < parset->GetNrOfDpbs(); ++uDpb) {
- TString sPrintoutLine = Form("DPB #%02u CROB Active ?: ", uDpb);
- for (uint32_t uCrobIdx = 0; uCrobIdx < fNrCrobPerDpb; ++uCrobIdx) {
- sPrintoutLine += Form("%1u", (vbCrobActiveFlag[uDpb][uCrobIdx] == true));
- }
- LOG(debug) << sPrintoutLine;
- }
-}
-
-// ---- printAddressMaps ----
-void CbmMuchUnpackAlgo::printAddressMaps(CbmMuchUnpackPar* parset,
- const std::vector<std::vector<std::vector<int32_t>>>& viFebModuleIdx,
- const std::vector<std::vector<std::vector<int32_t>>>& viFebModuleSide)
-{
-/*
- uint32_t uGlobalFebIdx = 0;
- for (uint32_t uDpb = 0; uDpb < parset->GetNrOfDpbs(); ++uDpb) {
- for (uint32_t uCrobIdx = 0; uCrobIdx < fNrCrobPerDpb; ++uCrobIdx) {
- LOG(debug) << Form("DPB #%02u CROB #%u: ", uDpb, uCrobIdx);
- for (uint32_t uFebIdx = 0; uFebIdx < parset->GetNbFebsPerCrob(); ++uFebIdx) {
- if (0 <= viFebModuleIdx[uDpb][uCrobIdx][uFebIdx])
- LOG(debug) << Form(" FEB #%02u (%02u): Mod. Idx = %03d Side %c (%2d) Type %c "
- "(%2d) (Addr. 0x%08x) ADC gain %4.0f e- ADC Offs %5.0f e-",
- uFebIdx, uGlobalFebIdx, viFebModuleIdx[uDpb][uCrobIdx][uFebIdx],
- 1 == viFebModuleSide[uDpb][uCrobIdx][uFebIdx] ? 'N' : 'P',
- viFebModuleSide[uDpb][uCrobIdx][uFebIdx],
- 1 == fviFebType[uDpb][uCrobIdx][uFebIdx] ? 'B' : 'A', fviFebType[uDpb][uCrobIdx][uFebIdx],
- fviFebAddress[uGlobalFebIdx], fvdFebAdcGain[uGlobalFebIdx], fvdFebAdcOffs[uGlobalFebIdx]);
- else
- LOG(debug) << Form("Disabled FEB #%02u (%02u): Mod. Idx = %03d Side %c (%2d) Type %c "
- "(%2d) (Addr. 0x%08x) ADC gain %4.0f e- ADC Offs %5.0f e-",
- uFebIdx, uGlobalFebIdx, viFebModuleIdx[uDpb][uCrobIdx][uFebIdx],
- 1 == viFebModuleSide[uDpb][uCrobIdx][uFebIdx] ? 'N' : 'P',
- viFebModuleSide[uDpb][uCrobIdx][uFebIdx],
- 1 == fviFebType[uDpb][uCrobIdx][uFebIdx] ? 'B' : 'A', fviFebType[uDpb][uCrobIdx][uFebIdx],
- fviFebAddress[uGlobalFebIdx], fvdFebAdcGain[uGlobalFebIdx], fvdFebAdcOffs[uGlobalFebIdx]);
- uGlobalFebIdx++;
- }
- }
- }
-*/
-}
-
-// -------------------------------------------------------------------------
-void CbmMuchUnpackAlgo::processErrorInfo(const stsxyter::Message& mess)
-{
- if (mess.IsMsErrorFlagOn()) {
- // I do pass here the Ts start time instead of the ms time, since, we removed the ms time as member for the time being
- if (fMonitor) { fMonitor->FillMsErrorsEvo(fMsStartTime, mess.GetMsErrorType()); }
- if (fOptOutBVec)
- fOptOutBVec->emplace_back(
- CbmErrorMessage(ECbmModuleId::kMuch, fMsStartTime, fuCurrDpbIdx, 0x20, mess.GetMsErrorType()));
- }
-}
-
-// ---- processHitInfo ----
-void CbmMuchUnpackAlgo::processHitInfo(const stsxyter::Message& mess)
-{
- const uint16_t usElinkIdx = mess.GetLinkIndexHitBinning();
- const uint32_t uCrobIdx = usElinkIdx / fNrElinksPerCrob;
- int32_t uFebIdx = fElinkIdxToFebIdxVec.at(usElinkIdx);
- if (-1 == uFebIdx) {
- LOG(warning) << "CbmMuchUnpackAlgo::processHitInfo => "
- << "Wrong elink Idx! Elink raw " << Form("%d remap %d", usElinkIdx, uFebIdx);
- return;
- }
- //uFebIdx += ( fuCurrDpbIdx * fNrCrobPerDpb + ____CrobIndexCalculationIfNeeded___() ) * fuNbFebsPerCrob;
- uFebIdx += (fuCurrDpbIdx * fNrCrobPerDpb) * fNrFebsPerCrob; //no Crob index calculation for now
-
- // Get the asic index
- uint32_t uAsicIdx = getAsicIndex(fuCurrDpbIdx, uCrobIdx, usElinkIdx);
-
- const uint16_t usChan = mess.GetHitChannel();
- const uint16_t usRawAdc = mess.GetHitAdc();
- const uint16_t usRawTs = mess.GetHitTimeBinning();
- const uint32_t uChanInFeb = usChan + fNrChsPerAsic * (uAsicIdx % fNrAsicsPerFeb);
-
- // Compute the Full time stamp
- const uint64_t ulHitTime = getFullTimeStamp(usRawTs);
-
- if (fMonitor) fMonitor->CountRawHit(uFebIdx, uChanInFeb);
-/*
- /// Store hit for output only if it is mapped to a module!!!
- if (fviFebAddress[uFebIdx] != 0 && fvbFebAdcCut[uFebIdx] < usRawAdc) {
- /// Store only if masking is disabled or if channeld is not masked
- /// 2D vector is safe as always right size if masking enabled
- if (false == fbUseChannelMask || false == fvvbMaskedChannels[uFebIdx][uChanInFeb]) {
-
- /// Duplicate hits rejection
- if (fbRejectDuplicateDigis) {
- if (usRawTs == fvvusLastTsChan[uAsicIdx][usChan]
- && (fbDupliWithoutAdc || usRawAdc == fvvusLastAdcChan[uAsicIdx][usChan])
- && fulTsMsbIndexInTs[fuCurrDpbIdx] == fvvulLastTsMsbChan[uAsicIdx][usChan]) {
- /// FIXME: add plots to check what is done in this rejection
- LOG(debug1) << "CbmMuchUnpackAlgo::processHitInfo => "
- << Form("Rejecting duplicate on Asic %3u channel %3u, TS %3u, ADC %2u", uAsicIdx, usChan, usRawTs,
- usRawAdc);
-
- if (fMonitor) fMonitor->FillDuplicateHitsAdc(uFebIdx, uChanInFeb, usRawAdc);
- return;
- } // if same TS, (ADC,) TS MSB, TS MSB cycle, reject
- fvvusLastTsChan[uAsicIdx][usChan] = usRawTs;
- fvvusLastAdcChan[uAsicIdx][usChan] = usRawAdc;
- fvvulLastTsMsbChan[uAsicIdx][usChan] = fulTsMsbIndexInTs[fuCurrDpbIdx];
- } // if (fbRejectDuplicateDigis)
-
- uint32_t uChanInMod = usChan + fNrChsPerAsic * (uAsicIdx % fNrAsicsPerFeb);
-
- /// FIXME: see issue #1549
- /// N side: 0-1023 => 0-1023
- /// P side: 0-1023 => 2047-1024
-
- if (0 == fviFebSide[uFebIdx])
- uChanInMod = fNrChsPerFeb - uChanInMod - 1 // Invert channel order
- + fNrChsPerFeb; // Offset for P (back) side
-
- const double_t tsreltime = static_cast<double_t>(ulHitTime) * stsxyter::kdClockCycleNs;
- double dTimeInNs = tsreltime - fSystemTimeOffset;
- if (uAsicIdx < fvdTimeOffsetNsAsics.size()) dTimeInNs -= fvdTimeOffsetNsAsics[uAsicIdx];
-
- const uint64_t ulTimeInNs = static_cast<uint64_t>(dTimeInNs);
- const double dCalAdc = fvdFebAdcOffs[uFebIdx] + (usRawAdc - 1) * fvdFebAdcGain[uFebIdx];
-
- if (0 == fviFebAddress[uFebIdx] || -1 == fviFebSide[uFebIdx]) {
- LOG(error) << Form("Digi on disabled FEB %02u has address 0x%08x and side %d", uFebIdx, fviFebAddress[uFebIdx],
- fviFebSide[uFebIdx]);
- }
-
- if (fMonitor) fMonitor->CountDigi(uFebIdx, uChanInFeb);
-
- /// Catch the pulser digis and either save them to their own output or drop them
- if (fvbFebPulser[uFebIdx]) {
- if (fOptOutAVec) {
- fOptOutAVec->emplace_back(CbmMuchDigi(fviFebAddress[uFebIdx], uChanInMod, ulTimeInNs, dCalAdc));
- } // if (fOptOutAVec)
- } // if (fvbFebPulser[uFebIdx])
- else {
- fOutputVec.emplace_back(CbmMuchDigi(fviFebAddress[uFebIdx], uChanInMod, ulTimeInNs, dCalAdc));
- } // else of if (fvbFebPulser[uFebIdx])
-
- /// If EM flag ON, store a corresponding error message with the next flag after all other possible status flags set
- if (mess.IsHitMissedEvts())
- if (fOptOutBVec)
- fOptOutBVec->emplace_back(
- CbmErrorMessage(ECbmModuleId::kMuch, dTimeInNs, uAsicIdx, 1 << stsxyter::kusLenStatStatus, usChan));
- }
- }
-*/
- // Convert the Hit time in bins to Hit time in ns
- const double dHitTimeNs = ulHitTime * stsxyter::kdClockCycleNs;
-
- if (fMonitor) {
- // Check Starting point of histos with time as X axis
- if (-1 == fdStartTime) { fdStartTime = dHitTimeNs; }
- if (fMonitor->GetDebugMode()) {
- fMonitor->FillHitDebugMonitoringHistos(uAsicIdx, usChan, usRawAdc, usRawTs, mess.IsHitMissedEvts());
- }
- const uint32_t uAsicInFeb = uAsicIdx % fNrAsicsPerFeb;
- const double dTimeSinceStartSec = (dHitTimeNs - fdStartTime) * 1e-9;
- const double dCalAdc = fvdFebAdcOffs[uFebIdx] + (usRawAdc - 1) * fvdFebAdcGain[uFebIdx];
- fMonitor->FillHitMonitoringHistos(uFebIdx, usChan, uChanInFeb, usRawAdc, dCalAdc, usRawTs, mess.IsHitMissedEvts());
- fMonitor->FillHitEvoMonitoringHistos(uFebIdx, uAsicIdx, uAsicInFeb, uChanInFeb, dTimeSinceStartSec,
- mess.IsHitMissedEvts());
- }
-}
-
-// ---- processStatusInfo ----
-void CbmMuchUnpackAlgo::processStatusInfo(const stsxyter::Message& mess, uint32_t uIdx)
-{
- // again fMonitor settings used for debugging printouts, I would propose to separate this
- if (fMonitor)
- if (fMonitor->GetDebugMode()) {
- std::cout << Form("DPB %2u TS %12lu mess %5u ", fuCurrDpbIdx, fTsIndex, uIdx);
- mess.PrintMess(std::cout, stsxyter::MessagePrintMask::msg_print_Human);
- }
-
- const uint16_t usElinkIdx = mess.GetStatusLink();
- const uint32_t uCrobIdx = usElinkIdx / fNrElinksPerCrob;
- const int32_t uFebIdx = fElinkIdxToFebIdxVec.at(usElinkIdx);
- if (-1 == uFebIdx) {
- LOG(warning) << "CbmMuchUnpackAlgo::processStatusInfo => "
- << "Wrong elink Idx! Elink raw " << Form("%d remap %d", usElinkIdx, uFebIdx);
- return;
- }
- const uint32_t uAsicIdx = getAsicIndex(fuCurrDpbIdx, uCrobIdx, usElinkIdx);
- if (fMonitor) {
- const uint16_t usStatusField = mess.GetStatusStatus();
- fMonitor->FillMuchStatusMessType(uAsicIdx, usStatusField);
- }
- /// Compute the Full time stamp
- const int64_t ulTime = getFullTimeStamp(0);
-
- /// Convert the time in bins to Hit time in ns
- const double dTimeNs = ulTime * stsxyter::kdClockCycleNs;
- if (fOptOutBVec)
- fOptOutBVec->emplace_back(
- CbmErrorMessage(ECbmModuleId::kMuch, dTimeNs, uAsicIdx, mess.GetStatusStatus(), mess.GetData()));
-}
-
-
-// ---- processTsMsbInfo ----
-void CbmMuchUnpackAlgo::processTsMsbInfo(const stsxyter::Message& mess, uint32_t uMessIdx, uint32_t uMsIdx)
-{
- const uint32_t uVal = mess.GetTsMsbValBinning();
-
- // Update Status counters
- if (uVal < fvulCurrentTsMsb[fuCurrDpbIdx]) {
-
- LOG(debug) << " TS " << std::setw(12) << fTsIndex << " MS Idx " << std::setw(4) << uMsIdx << " Msg Idx "
- << std::setw(5) << uMessIdx << " DPB " << std::setw(2) << fuCurrDpbIdx << " Old TsMsb " << std::setw(5)
- << fvulCurrentTsMsb[fuCurrDpbIdx] << " Old MsbCy " << std::setw(5) << fvuCurrentTsMsbCycle[fuCurrDpbIdx]
- << " new TsMsb " << std::setw(5) << uVal;
-
- fvuCurrentTsMsbCycle[fuCurrDpbIdx]++;
- }
- if (uVal != fvulCurrentTsMsb[fuCurrDpbIdx] + 1
- /// Case where we reach a normal cycle edge
- && !(0 == uVal && stsxyter::kuTsMsbNbTsBinsBinning == fvulCurrentTsMsb[fuCurrDpbIdx])
- /// First TS_MSB in MS may jump if TS dropped by DAQ
- && !(1 == uMessIdx && 0 == uMsIdx)
- /// Msg 1 and 2 will be same TS_MSB if data in 1st bin
- && !(uVal == fvulCurrentTsMsb[fuCurrDpbIdx] && 2 == uMessIdx)
- /// New FW introduced TS_MSB suppression + large TS_MSB => warning only if value not increasing
- && uVal < fvulCurrentTsMsb[fuCurrDpbIdx]) {
- LOG(debug) << "TS MSB Jump in "
- << " TS " << std::setw(12) << fTsIndex << " MS Idx " << std::setw(4) << uMsIdx << " Msg Idx "
- << std::setw(5) << uMessIdx << " DPB " << std::setw(2) << fuCurrDpbIdx << " => Old TsMsb "
- << std::setw(5) << fvulCurrentTsMsb[fuCurrDpbIdx] << " new TsMsb " << std::setw(5) << uVal;
- }
-
- fvulCurrentTsMsb[fuCurrDpbIdx] = uVal;
-
- LOG(debug1) << " TS " << std::setw(12) << fTsIndex << " MS Idx " << std::setw(4) << uMsIdx << " Msg Idx "
- << std::setw(5) << uMessIdx << " DPB " << std::setw(2) << fuCurrDpbIdx << " TsMsb " << std::setw(5)
- << fvulCurrentTsMsb[fuCurrDpbIdx] << " MsbCy " << std::setw(5) << fvuCurrentTsMsbCycle[fuCurrDpbIdx];
-
- fulTsMsbIndexInTs[fuCurrDpbIdx] =
- fvulCurrentTsMsb[fuCurrDpbIdx]
- + (fvuCurrentTsMsbCycle[fuCurrDpbIdx] * static_cast<uint64_t>(1 << stsxyter::kusLenTsMsbValBinning));
- if (fulTsMsbIndexInTs[fuCurrDpbIdx] < fulTsStartInTsMsb) {
- LOG(fatal) << "CbmMuchUnpackAlgo::processTsMsbInfo => "
- << "Value computed from TS_MSB and TS_MSB cycle smaller than Timeslice start in TS_MSB, "
- << "would lead to a negative value so it cannot be recovered!!!!"
- << std::endl
- /// Values Printout
- << "TS_MSB: " << fvulCurrentTsMsb[fuCurrDpbIdx] << " Cycle: " << fvuCurrentTsMsbCycle[fuCurrDpbIdx]
- << " Full TS_MSB: " << fulTsMsbIndexInTs[fuCurrDpbIdx] << " TS Start offset: " << fulTsStartInTsMsb;
- }
- fulTsMsbIndexInTs[fuCurrDpbIdx] -= fulTsStartInTsMsb;
-
- if (fMonitor)
- if (fMonitor->GetDebugMode()) { //also if( 1 < uMessIdx )?
- fMonitor->FillMuchDpbRawTsMsb(fuCurrDpbIdx, fvulCurrentTsMsb[fuCurrDpbIdx]);
- fMonitor->FillMuchDpbRawTsMsbSx(fuCurrDpbIdx, fvulCurrentTsMsb[fuCurrDpbIdx]);
- fMonitor->FillMuchDpbRawTsMsbDpb(fuCurrDpbIdx, fvulCurrentTsMsb[fuCurrDpbIdx]);
- }
-}
-
-// ---- refreshTsMsbFields ----
-void CbmMuchUnpackAlgo::refreshTsMsbFields(const uint32_t imslice, const size_t mstime)
-{
- const uint32_t uTsMsbCycleHeader =
- std::floor(mstime / (stsxyter::kulTsCycleNbBinsBinning * stsxyter::kdClockCycleNs));
- const uint32_t uTsMsbHeader =
- std::floor((mstime - uTsMsbCycleHeader * (stsxyter::kulTsCycleNbBinsBinning * stsxyter::kdClockCycleNs))
- / (stsxyter::kuHitNbTsBinsBinning * stsxyter::kdClockCycleNs));
-
- LOG(debug1) << " TS " << std::setw(12) << fTsIndex << " MS " << std::setw(12) << mstime << " MS Idx " << std::setw(4)
- << imslice << " Msg Idx " << std::setw(5) << 0 << " DPB " << std::setw(2) << fuCurrDpbIdx << " Old TsMsb "
- << std::setw(5) << fvulCurrentTsMsb[fuCurrDpbIdx] << " Old MsbCy " << std::setw(5)
- << fvuCurrentTsMsbCycle[fuCurrDpbIdx] << " header TsMsb " << uTsMsbHeader << " New MsbCy "
- << uTsMsbCycleHeader;
-
- if (0 == imslice) {
- if (uTsMsbCycleHeader != fvuCurrentTsMsbCycle[fuCurrDpbIdx])
- LOG(debug) << " TS " << std::setw(12) << fTsIndex << " MS " << std::setw(12) << mstime << " MS Idx "
- << std::setw(4) << imslice << " Msg Idx " << std::setw(5) << 0 << " DPB " << std::setw(2)
- << fuCurrDpbIdx << " Old TsMsb " << std::setw(5) << fvulCurrentTsMsb[fuCurrDpbIdx] << " Old MsbCy "
- << std::setw(5) << fvuCurrentTsMsbCycle[fuCurrDpbIdx] << " New TsMsb " << uTsMsbHeader << " New MsbCy "
- << uTsMsbCycleHeader;
- fvuCurrentTsMsbCycle[fuCurrDpbIdx] = uTsMsbCycleHeader;
- fvulCurrentTsMsb[fuCurrDpbIdx] = uTsMsbHeader;
-
- fulTsMsbIndexInTs[fuCurrDpbIdx] =
- fvulCurrentTsMsb[fuCurrDpbIdx]
- + (fvuCurrentTsMsbCycle[fuCurrDpbIdx] * static_cast<uint64_t>(1 << stsxyter::kusLenTsMsbValBinning));
- if (fulTsMsbIndexInTs[fuCurrDpbIdx] < fulTsStartInTsMsb) {
- LOG(fatal) << "CbmMuchUnpackAlgo::refreshTsMsbFields => "
- << "Value computed from TS_MSB and TS_MSB cycle smaller than Timeslice start in TS_MSB, "
- << "would lead to a negative value so it cannot be recovered!!!!"
- << std::endl
- /// Values Printout
- << "TS_MSB: " << fvulCurrentTsMsb[fuCurrDpbIdx] << " Cycle: " << fvuCurrentTsMsbCycle[fuCurrDpbIdx]
- << " Full TS_MSB: " << fulTsMsbIndexInTs[fuCurrDpbIdx] << " TS Start offset: " << fulTsStartInTsMsb;
- }
- fulTsMsbIndexInTs[fuCurrDpbIdx] -= fulTsStartInTsMsb;
- }
- else if (uTsMsbCycleHeader != fvuCurrentTsMsbCycle[fuCurrDpbIdx]) {
- if ((stsxyter::kuTsMsbNbTsBinsBinning - 1) == fvulCurrentTsMsb[fuCurrDpbIdx]) {
- /// Transition at the edge of two MS, the first TS_MSB message (idx 1) will make the cycle
- LOG(debug) << " TS " << std::setw(12) << fTsIndex << " MS " << std::setw(12) << mstime << " MS Idx "
- << std::setw(4) << imslice << " Msg Idx " << std::setw(5) << 0 << " DPB " << std::setw(2)
- << fuCurrDpbIdx << " Old TsMsb " << std::setw(5) << fvulCurrentTsMsb[fuCurrDpbIdx] << " Old MsbCy "
- << std::setw(5) << fvuCurrentTsMsbCycle[fuCurrDpbIdx] << " New MsbCy " << uTsMsbCycleHeader;
- }
- else {
- LOG(warning) << "TS MSB cycle from MS header does not match current cycle from data "
- << "for TS " << std::setw(12) << fTsIndex << " MS " << std::setw(12) << mstime << " MsInTs "
- << std::setw(3) << imslice << " ====> " << fvuCurrentTsMsbCycle[fuCurrDpbIdx] << " (cnt) VS "
- << uTsMsbCycleHeader << " (header)";
- fvuCurrentTsMsbCycle[fuCurrDpbIdx] = uTsMsbCycleHeader;
- }
- }
-}
-
-// ---- unpack ----
-bool CbmMuchUnpackAlgo::unpack(const fles::Timeslice* ts, std::uint16_t icomp, UInt_t imslice)
-{
- auto msDescriptor = ts->descriptor(icomp, imslice);
-
- /// If monitoring enable, fill here the "Per Timeslice" histograms with info from the previous TS
- /// as only way to detect TS transitions due to the multiple calls in case of multiple components.
- /// Time set to tlast MS in previous TS. Last TS will be missed.
- if (fMonitor && 0 == imslice && 0 < fMsStartTime && msDescriptor.idx == ts->start_time()
- && fMsStartTime < msDescriptor.idx) {
- fMonitor->FillPerTimesliceCountersHistos(static_cast<double_t>(fMsStartTime) / 1e9);
- }
-
- //Current equipment ID, tells from which DPB the current MS is originating
- const uint32_t uCurrentEquipmentId = msDescriptor.eq_id;
- const uint8_t* msContent = reinterpret_cast<const uint8_t*>(ts->content(icomp, imslice));
- const uint32_t uSize = msDescriptor.size;
-
- fMsStartTime = msDescriptor.idx;
- LOG(debug) << "Microslice: " << fMsStartTime << " from EqId " << std::hex << uCurrentEquipmentId << std::dec
- << " has size: " << uSize << " (index " << imslice << ")";
-
- if (0 == fvbMaskedComponents.size()) fvbMaskedComponents.resize(ts->num_components(), false);
-
- //Temp holder until current equipment ID is properly filled in MS
- const uint32_t uCurrDpbId = static_cast<uint32_t>(uCurrentEquipmentId & 0xFFFF);
-
- if (fMonitor)
- if (fMonitor->GetDebugMode()) {
- const double dMsTime = (1e-9) * static_cast<double>(fMsStartTime);
- if (icomp < fMonitor->GetMaxNbFlibLinks()) {
- if (fdStartTimeMsSz < 0) fdStartTimeMsSz = dMsTime;
- fMonitor->CreateMsComponentSizeHistos(icomp);
- fMonitor->FillMsSize(icomp, uSize);
- fMonitor->FillMsSizeTime(icomp, dMsTime - fdStartTimeMsSz, uSize);
- }
- }
-
- /// Check if this sDPB ID was declared in parameter file and stop there if not
- auto it = fDpbIdIndexMap.find(uCurrDpbId);
- if (it == fDpbIdIndexMap.end()) {
- if (false == fvbMaskedComponents[icomp]) {
- // LOG(debug) << "---------------------------------------------------------------";
- // Had to remove this line otherwise we would get circle dependencies in the current stage of cbmroot, since we still have Unpackers in the fles folders, which require the reco folders
- // LOG(debug) << FormatMsHeaderPrintout(msDescriptor);
- LOG(warning) << fName << "::unpack(...)::Could not find the sDPB index for AFCK id 0x" << std::hex << uCurrDpbId
- << std::dec << " in timeslice " << fNrProcessedTs << " in microslice " << imslice << " component "
- << icomp << "\n"
- << "If valid this index has to be added in the STS "
- "parameter file in the DbpIdArray field";
- fvbMaskedComponents[icomp] = true;
-
- /// If first TS being analyzed, we are probably detecting STS/MUCH boards with same sysid
- /// => Do not report the MS as bad, just ignore it
- if (1 == fNrProcessedTs) return true;
- }
- else
- return true;
-
- // REMARK please check experts, but I think this point can never be reached PR 072021
- return false;
- }
- else
- fuCurrDpbIdx = fDpbIdIndexMap[uCurrDpbId];
-
- if (fMonitor) fMonitor->FillMsCntEvo(fMsStartTime);
-
- if (0 == imslice) {
- /// Extract the time base only on MS 0, assuming that we get all TS of a component in order
- fulTsStartInTsMsb =
- static_cast<uint64_t>(fTsStartTime / (stsxyter::kuHitNbTsBinsBinning * stsxyter::kdClockCycleNs));
- }
-
- // Check the current TS_MSb cycle and correct it if wrong
- refreshTsMsbFields(imslice, fMsStartTime);
-
- //Reset internal monitor variables for debugging info
- if (fMonitor) {
- if (fMonitor->GetDebugMode()) { fMonitor->ResetDebugInfo(); }
- }
-
- //Main processing of MS content
- loopMsMessages(msContent, uSize, imslice);
-
- //Output debugging info
- if (fMonitor) {
- if (fMonitor->GetDebugMode()) { fMonitor->PrintDebugInfo(fMsStartTime, fNrProcessedTs, msDescriptor.flags, uSize); }
- for (auto itHit = fOutputVec.begin(); itHit != fOutputVec.end(); ++itHit) {
- fMonitor->FillDigisTimeInRun(itHit->GetTime());
- }
- fMonitor->FillVectorSize(ts->index(), fOutputVec.size());
- //fMonitor->DrawCanvases();
- }
- return true;
-}
-
-ClassImp(CbmMuchUnpackAlgo)
diff --git a/reco/detectors/much/unpack/CbmMuchUnpackAlgo.h b/reco/detectors/much/unpack/CbmMuchUnpackAlgo.h
index 2c8918b1730396029486f72cbf1d6d459d677613..27317fab5343a58818e792937215ccd72b578aa0 100644
--- a/reco/detectors/much/unpack/CbmMuchUnpackAlgo.h
+++ b/reco/detectors/much/unpack/CbmMuchUnpackAlgo.h
@@ -19,7 +19,6 @@
#include "CbmErrorMessage.h" // REMARK see remark in CbmMuchUnpackConfig
-#include "CbmMcbm2018MuchPar.h"
#include "CbmMuchDigi.h"
#include "CbmMuchUnpackMonitor.h"
#include "CbmMuchUnpackPar.h"
@@ -95,6 +94,14 @@ public:
/** @brief Set a predefined monitor @param monitor predefined unpacking monitor */
void SetMonitor(std::shared_ptr<CbmMuchUnpackMonitor> monitor) { fMonitor = monitor; }
+ /**
+ * @brief Sets the name of the parameter file to be used.
+ *
+ * @param[in] std:string, path should not be included as set in the Config class
+ */
+ void SetParFileName(std::string sNewName) { fParFileName = sNewName; }
+
+
protected:
/**
* @brief Get the Asic Index
@@ -190,7 +197,7 @@ protected:
/** @brief experts please add description here */
void refreshTsMsbFields(const uint32_t imslice, const size_t mstime);
- uint32_t CreateMuchAddress(uint32_t dpbidx, uint32_t febId, uint32_t usChan);
+ uint32_t CreateMuchAddress(uint32_t dpbidx, int32_t febId, uint32_t usChan);
/**
* @brief Set the Derived Ts Parameters
@@ -335,7 +342,9 @@ private:
uint64_t fuSkippMessNoValidAddress = 0;
uint64_t fuSkippMessNoValidFebPosition = 0;
uint64_t fuSkippMessNoValidElink = 0;
- ClassDef(CbmMuchUnpackAlgo, 1)
+ /// Settings from parameter file
+ std::string fParFileName = "mMuchPar.par";
+ ClassDef(CbmMuchUnpackAlgo, 2)
};
#endif // CbmMuchUnpackAlgo_H
diff --git a/reco/detectors/much/unpack/CbmMuchUnpackConfig.h b/reco/detectors/much/unpack/CbmMuchUnpackConfig.h
index 7ff72c2455d389727f958d9690d6e17ad668bb92..2f21646a58bda40b9950c3aeb33cf2cdfe29e72c 100644
--- a/reco/detectors/much/unpack/CbmMuchUnpackConfig.h
+++ b/reco/detectors/much/unpack/CbmMuchUnpackConfig.h
@@ -105,6 +105,14 @@ public:
/** @brief Add a monitor to the unpacker. @param value CbmMuchUnpackMonitor */
void SetMonitor(std::shared_ptr<CbmMuchUnpackMonitor> value) { fMonitor = value; }
+ /**
+ * @brief Sets the name of the parameter file to be used.
+ *
+ * @param[in] std:string, path should not be included as set in the Config class
+ */
+ void SetParFileName(std::string sNewName) { fsParFileName = sNewName; }
+ void LoadParFileName() { fAlgo->SetParFileName(fsParFileName); }
+
protected:
/**
@@ -135,8 +143,11 @@ protected:
/// Temporary storage of user parameters
std::vector<FebChanMaskReco> fvChanMasks = {};
+ /// Parameter file name
+ std::string fsParFileName = "mTofCriPar.par";
+
private:
- ClassDef(CbmMuchUnpackConfig, 1)
+ ClassDef(CbmMuchUnpackConfig, 2)
};
#endif // CbmMuchUnpackConfig_H