StDetectorDbMaker/St_spaceChargeCorC.cxx

#include "St_spaceChargeCorC.h"
#include "StDetectorDbRichScalers.h"
#include "StMemStat.h"

 Double_t St_spaceChargeCorC::getSpaceChargeCoulombs(Double_t scaleFactor)
 {
    StDetectorDbRichScalers* scalers = StDetectorDbRichScalers::instance();
    if (! scalers ) return 0;
    Double_t zf = zeroField(0); // potential validity margin for scalers
    if (zf>0 && zf<1) scalers->setValidityMargin(zf);
    Double_t coulombs = 0;
//    StMemStat::PrintMem("Space charge before GetDate");
    int idate  = StMaker::GetChain()->GetDate();
//    StMemStat::PrintMem("Space charge AFTER int GetDate() conversion");  

    bool use_powers = idate > 20090101;
//    bool use_powers = (StMaker::GetChain()->GetDate() > 20090101);
//    StMemStat::PrintMem("Space charge AFTER comparision");
    for (int row=0;row< (int) getNumRows();row++) {
      Double_t mult = 0;
      switch ((int) getSpaceChargeDetector(row)) {
        case (0) : mult = scalers->getMult(); break;
        case (1) : mult = scalers->getBBCX(); break;
        case (2) : mult = scalers->getZDCX(); break;
        case (3) : mult = scalers->getZDCEast()+scalers->getZDCWest(); break;
        case (4) : mult = scalers->getBBCEast()+scalers->getBBCWest(); break;
        case (5) : mult = scalers->getZDCEast(); break;
        case (6) : mult = scalers->getZDCWest(); break;
        case (7) : mult = scalers->getBBCEast(); break;
        case (8) : mult = scalers->getBBCWest(); break;
        case (9) : mult = scalers->getBBCYellowBkg(); break;
        case (10): mult = scalers->getBBCBlueBkg(); break;
        default  : mult = 0.;
      }
      if (mult < 0) {
        Mark();
        return 0; // Unphysical scaler rates will be uncorrected
      } else UnMark();
      Double_t saturation = getSpaceChargeSatRate(row);
      Double_t correction = getSpaceChargeCorrection(scaleFactor,row);
      Double_t factor     = getSpaceChargeFactor(row);
      Double_t offset     = getSpaceChargeOffset(row);
      Double_t intens = (mult < saturation) ? mult : saturation;
      if (use_powers) coulombs += ::pow(intens-offset,factor) * correction ;
      else coulombs += factor * (intens-offset) * correction ;
    }
    return coulombs;
 }

---