Calibrated BBC simulation (#718)

BBC simulation has been calibrated with a 2017 p+p run. There is added the calibrated factors for all BBC channels. It was also discussed with Akio. Thanks --------- Co-authored-by: SEUNGHWAN LEE <[email protected]>
star-bnl · Dec 13, 2024 · c68927c · c68927c
1 parent 640288b
commit c68927c
Showing 1 changed file with 30 additions and 7 deletions.
diff --git a/StRoot/StBbcSimulationMaker/StBbcSimulationMaker.cxx b/StRoot/StBbcSimulationMaker/StBbcSimulationMaker.cxx
@@ -52,6 +52,31 @@ const float SinglePhotoElectronResolution = 0.3; // according to Les Bland
 /* Numbering: the real PMT numbering (used in the map) starts from 1.
 ALL OTHERS start from 0 !
  */
+
+
+// Note: BBC new factors are calibrated by 2017 p+p runs
+const short ADCBin = 4096;
+
+const float LightFactor[NPMT2] = {1., 1., 0.95, 0.92, 1., 1., 1., 1., 0.9, 0.96, 0.95, 1., 0.95, 0.97, 1., 1.,    // East small 
+                                  0.33, 0.298, 0.25, 0.26, 0.22, 0.4, 0.24, 0.26,                                 // East large
+                                  1.2, 0.9, 1., 0.6, 1., 1., 1.2, 1.2, 1.25, 0.4, 1., 1., 0.95, 0.95, 0.94, 0.95, // West small
+                                  0., 0.243, 0., 0.2, 0.24, 0.18, 0.19, 0.2};                                     // West large
+
+const float TileResolution[NPMT2] = {0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, // East small
+                                     0.6, 0.6, 0.4, 0.5, 0.5, 0.6, 0.4, 0.5,                                         // East large
+                                     0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.5, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, // West small
+                                     0., 0.5, 0., 0.5, 0.6, 0.6, 0.4, 0.65};                                         // West large
+
+const float pC_per_ADCBins[NPMT2] = {0.03, 0.03, 0.033, 0.033, 0.03, 0.03, 0.022, 0.03, 0.021, 0.022, 0.024, 0.03, 0.024, 0.024, 0.025, 0.027, // East small
+                                     0.0226, 0.028, 0.021, 0.023, 0.014,  0.0125, 0.012, 0.0125,                                               // East large
+                                     0.031, 0.03, 0.03, 0.018, 0.03, 0.03, 0.023, 0.024, 0.028, 0.005, 0.03, 0.03, 0.028, 0.023, 0.022, 0.028, // West small
+                                     0., 0.027, 0., 0.06, 0.013, 0.014, 0.0129, 0.02};                                                         // West large
+
+const int shift_ADC0[NPMT2] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,   // East small
+                               0, 0, 0, 0, -15, 0, -15, 0,                       // East large
+                               0, 0, 0, 0, 0, 0, 0, 0, 0, -60, 0, 0, 0, 0, 0, 0, // West small
+                               0, 0, 0, 0, -10, -7, -15, -10};                   // West large
+
 //____________________________________________________________________________
 
 bool IsSmall(short iPMT)
@@ -118,27 +143,25 @@ the response based on that + noise.
   ~BbcDE(){};
   void AddDE(uint16_t ipmt, float de)
   {
-    if (!IsSmall(ipmt)) {de *= OuterFactor;}
+    de *= LightFactor[ipmt];
     dE[ipmt] += de;
   }
   float GetDE(uint16_t ipmt)
   {
     /// returns DE in pC of PMT signal
-
     float PoissonMean = dE[ipmt]/dE_1MIP*NPhotoelectrons_1MIP;
     short NPhotoelectrons = BbcRndm.Poisson(PoissonMean);
     float Q = pC_per_Photoelectron*
-      (1+BbcRndm.Gaus(0.,SinglePhotoElectronResolution))*
-       NPhotoelectrons;
+      (1+BbcRndm.Gaus(0.,TileResolution[ipmt]))*NPhotoelectrons;
     return Q;
   }
   short GetADC(uint16_t ipmt)
   {
     /// returns digitized (ADC) amplitude
     float A = this->GetDE(ipmt);
-    if (A<ADC0) {return 0;}
-    short N = (short)((A-ADC0)/pC_perADCbin);
-    if (N>=NADCbins) {return NADCbins-1;}
+    short N = (short)(A/pC_per_ADCBins[ipmt] + shift_ADC0[ipmt]);
+    if (N>=ADCBin) {return ADCBin-1;}
+    if (N<0){return 0;}
     return N;
   }
 };