Added Siemens Quadra scanner in the scanner.cxx and .h (#1628)

Added Siemens Quadra scanner in the scanner.cxx and .h, using a work-around
to be able to insert a virtual axial crystal between every "vision 600".

---------

Co-authored-by: Z.K Li <54258479+christanmod@users.noreply.github.com>

Author: z-k-li

src/buildblock/Scanner.cxx

@@ -476,6 +476,36 @@ Scanner::Scanner(Type scanner_type)
       );
       break;
 
+    case Siemens_Quadra:
+      // Info from Prenosil et al., J Nucl Med 2022; 63(3):476-484, DOI: 10.2967/jnumed.121.261972
+      set_params(Siemens_Quadra,                                  // type
+                 string_list("Siemens Quadra", "Quadra", "1232"), // names
+                 323,                                             // rings, because of the max ring difference is 322
+                 520,                                             // max n non-arc-corr bins
+                 520,                                             // default n arc-corr bins
+                 (20 + 1) * 38,                                   // num detector per ring
+                 410.0F,                                          // inner ring radius (mm)
+                 7.0F,                                            // unsure on this                          // avg DoI (mm)
+                 3.29114F,                                        // ring spacing (mm)
+                 1.6F,                                            // bin size (mm)
+                 0.0F,                                            // intrinsic tilt
+                 // ONLY used for CTI scanners for normalisation. These values will be ignored here.
+                 4,
+                 1, // n axial/trans blocks per bucket
+                 10 * 8 + 1,
+                 20 + 1, // n axial/trans xtals per block
+                 0,
+                 0, // n axial/trans xtals per singles unit
+                 1, // n detector layers
+                 // energy
+                 0.F,
+                 511.F,
+                 33 * 8,      // max n timing position
+                 143.231 / 8, // timing position bin size
+                 214.F        // timing resolution
+      );
+      break;
+
     case RPT:
 
       set_params(RPT,
@@ -1608,6 +1638,7 @@ Scanner::get_num_virtual_axial_crystals_per_block() const
     {
     case E1080:
     case Siemens_mCT:
+    case Siemens_Quadra:
       return 1;
     default:
       return 0;
@@ -1624,6 +1655,7 @@ Scanner::get_num_virtual_transaxial_crystals_per_block() const
     case E1080:
     case Siemens_mCT:
     case Siemens_Vision_600:
+    case Siemens_Quadra:
     case Siemens_mMR:
     case UPENN_5rings:
     case UPENN_6rings:

src/include/stir/Scanner.h

@@ -152,6 +152,7 @@ class Scanner
     Siemens_mMR,
     Siemens_mCT,
     Siemens_Vision_600,
+    Siemens_Quadra,
     RPT,
     HiDAC,
     Advance,