Update centrality definition for FXT (#584)

Update of the centrality definition for FXT (3, 3.2, 3.5, 3.9, and 7.7 GeV) and Collider (7.7 and 200 GeV) --------- Co-authored-by: Vinh Luong <[email protected]> Co-authored-by: lbavinh <[email protected]>
star-bnl · Nov 13, 2023 · 93381cc · 93381cc
1 parent f1f9389
commit 93381cc
Show file tree

Hide file tree

Showing 8 changed files with 1,029 additions and 210 deletions.
diff --git a/StRoot/StRefMultCorr/BadRun.h b/StRoot/StRefMultCorr/BadRun.h
diff --git a/StRoot/StRefMultCorr/CentralityMaker.cxx b/StRoot/StRefMultCorr/CentralityMaker.cxx
@@ -46,6 +46,7 @@ CentralityMaker::CentralityMaker() {
   fRefMult2Corr = new StRefMultCorr("refmult2") ;
   fRefMult3Corr = new StRefMultCorr("refmult3") ;
   fRefMultCorr_Isobar  = new StRefMultCorr("refmult","Isobar") ;
+  fRefMultCorrFxt = new StRefMultCorr("fxtmult");
  // fTofTrayMultCorr = new StRefMultCorr("toftray") ;
   fgRefMultCorr  = new StRefMultCorr("grefmult") ;
   fgRefMultCorr_Run14_AuAu200_VpdMB5_P16id = new StRefMultCorr("grefmult","Run14_AuAu200_VpdMB5","P16id") ;
@@ -91,6 +92,11 @@ StRefMultCorr* CentralityMaker::getRefMultCorr_Isobar() {
   return fRefMultCorr_Isobar ;
 }
 
+//_________________
+StRefMultCorr* CentralityMaker::getRefMultCorrFxt() {
+  return fRefMultCorrFxt;
+}
+
 /*
 //_________________
 StRefMultCorr* CentralityMaker::getTofTrayMultCorr() {

diff --git a/StRoot/StRefMultCorr/CentralityMaker.h b/StRoot/StRefMultCorr/CentralityMaker.h
@@ -43,7 +43,7 @@
 #ifndef __CentralityMaker_h__
 #define __CentralityMaker_h__
 
-class StRefMultCorr ;
+class StRefMultCorr;
 #include "Rtypes.h"
 
 //_________________
@@ -57,6 +57,7 @@ class CentralityMaker {
     StRefMultCorr* getRefMult2Corr() ; // For refmult2
     StRefMultCorr* getRefMult3Corr() ; // For refmult3
     StRefMultCorr* getRefMultCorr_Isobar()  ; // For refmult
+    StRefMultCorr* getRefMultCorrFxt(); // For fixed-target data
     // StRefMultCorr* getTofTrayMultCorr() ; // For TOF tray multiplicity
     StRefMultCorr* getgRefMultCorr()  ; // For grefmult //Run14 AuAu200GeV
     StRefMultCorr* getgRefMultCorr_Run14_AuAu200_VpdMB5_P16id()  ; 
@@ -74,12 +75,13 @@ class CentralityMaker {
     static CentralityMaker* fInstance ; // Static pointer of CentralityMaker
 
     // Centrality correction classes
-    StRefMultCorr* fRefMultCorr  ; // refmult based centrality
-    StRefMultCorr* fRefMult2Corr ; // refmult2 based centrality
-    StRefMultCorr* fRefMult3Corr ; // refmult3 based centrality
+    StRefMultCorr* fRefMultCorr  ; // refmult-based centrality
+    StRefMultCorr* fRefMult2Corr ; // refmult2-based centrality
+    StRefMultCorr* fRefMult3Corr ; // refmult3-based centrality
     StRefMultCorr* fRefMultCorr_Isobar  ; // refmult based centrality
-    // StRefMultCorr* fTofTrayMultCorr ; // tofTrayMult based centrality
-    StRefMultCorr* fgRefMultCorr  ; // grefmult based centrality
+    StRefMultCorr* fRefMultCorrFxt; // fxtMult-based centrality
+    // StRefMultCorr* fTofTrayMultCorr ; // tofTrayMult-based centrality
+    StRefMultCorr* fgRefMultCorr  ; // grefmult-based centrality
     StRefMultCorr* fgRefMultCorr_Run14_AuAu200_VpdMB5_P16id  ;
     StRefMultCorr* fgRefMultCorr_Run14_AuAu200_VpdMB30_P16id  ;
     StRefMultCorr* fgRefMultCorr_Run14_AuAu200_VpdMBnoVtx_P16id  ;

diff --git a/StRoot/StRefMultCorr/Param.cxx b/StRoot/StRefMultCorr/Param.cxx
@@ -17,6 +17,7 @@ const string getParamX( const int x, const int y1, const int y2 ) {
     case 2:  str = mParamStr_ref2[y1][y2]; break;
     case 3:  str = mParamStr_ref3[y1][y2]; break;
     case 4:  str = mParamStr_ref4[y1][y2]; break;
+    case 5:  str = mParamStr_ref5[y1][y2]; break;
     default: str = "0";
   }
   return str;