Update calibration db to provide nonlinearity params

Author: asanchezcastillo

sbndcode/Calibration/PDSDatabaseInterface/PMTCalibrationDatabase.h

@@ -32,6 +32,8 @@ namespace sbndDB {
     virtual double getSPEAmplitudeStd(unsigned int channelID) const = 0;
     virtual double getGaussFilterPower(unsigned int channelID) const = 0;
     virtual double getGaussFilterWC(unsigned int channelID) const = 0;
+    virtual double getNonLineatiryPESat(unsigned int channelID) const = 0;
+    virtual double getNonLineatiryAlpha(unsigned int channelID) const = 0;
     virtual std::vector<double> getSER(unsigned int channelID) const = 0;
   }; // end class
 

sbndcode/Calibration/PDSDatabaseInterface/PMTCalibrationDatabaseProvider.cxx

@@ -178,6 +178,21 @@ void sbndDB::PMTCalibrationDatabaseProvider::ReadPMTCalibration(uint32_t run)
         << dbname << "\n";
     fPMTCalibrationData[channel].gauss_wc = _gauss_wc;
 
+    double _nonlinearity_pesat = 0.;
+    error = db.GetNamedChannelData(channel, "nonlinearity_pesat", _nonlinearity_pesat);
+    if (error)
+      throw cet::exception("PMTTimingCorrectionsProvider")
+        << "Encountered error (code " << error
+        << ") while trying to access 'nonlinearity_pesat' on table " << dbname << "\n";
+    fPMTCalibrationData[channel].nonlinearity_pesat = _nonlinearity_pesat;
+    double _nonlinearity_alpha = 0.;
+    error = db.GetNamedChannelData(channel, "nonlinearity_alpha", _nonlinearity_alpha);
+    if (error)
+      throw cet::exception("PMTTimingCorrectionsProvider")
+        << "Encountered error (code " << error
+        << ") while trying to access 'nonlinearity_alpha' on table " << dbname << "\n";
+    fPMTCalibrationData[channel].nonlinearity_alpha = _nonlinearity_alpha;
+    
     // Read SER
     std::vector<double> _ser;
     std::string name_base = "ser_vec_";

sbndcode/Calibration/PDSDatabaseInterface/PMTCalibrationDatabaseProvider.h

@@ -71,6 +71,12 @@ class sbndDB::PMTCalibrationDatabaseProvider : public PMTCalibrationDatabase {
         double getGaussFilterWC( unsigned int channelID ) const override {
             return getChannelCorrOrDefault(channelID).gauss_wc;
         };
+        double getNonLineatiryPESat( unsigned int channelID ) const override {
+            return getChannelCorrOrDefault(channelID).nonlinearity_pesat;
+        };
+        double getNonLineatiryAlpha( unsigned int channelID ) const override {
+            return getChannelCorrOrDefault(channelID).nonlinearity_alpha;
+        };
         std::vector<double> getSER( unsigned int channelID ) const override {
             return getChannelCorrOrDefault(channelID).ser;
         };
@@ -96,11 +102,12 @@ class sbndDB::PMTCalibrationDatabaseProvider : public PMTCalibrationDatabase {
             double spe_amplitude_std=0.;
             double gauss_wc_power=0.;
             double gauss_wc=0.;
+            double nonlinearity_pesat=0.;
+            double nonlinearity_alpha=0.;
             std::vector<double> ser={};
             };
 
-        const PMTCalibrationDB CorrectionDefaults = {0, 0, 0, 0, 0.0, 0.0, 0.0, {}};
-
+        const PMTCalibrationDB CorrectionDefaults = {0, 0, 0, 0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, {}};
 	    /// Map of corrections by channel
         std::map<unsigned int, PMTCalibrationDB> fPMTCalibrationData;