Add independent variable to store a normalization energy for TFCSPCAEnergyParametrization

Simulation/ISF/ISF_FastCaloSim/ISF_FastCaloSimEvent/ISF_FastCaloSimEvent/TFCSPCAEnergyParametrization.h

@@ -37,6 +37,9 @@ class TFCSPCAEnergyParametrization:public TFCSEnergyParametrization
   void clean();
   
   void Print(Option_t *option = "") const override;
+  
+  float get_total_energy_normalization() const {return m_total_energy_normalization;};
+  void  set_total_energy_normalization(float norm) {m_total_energy_normalization=norm;};
 
   int                       do_rescale;
   
@@ -53,7 +56,9 @@ class TFCSPCAEnergyParametrization:public TFCSEnergyParametrization
   
   int m_numberpcabins;
   
-  ClassDefOverride(TFCSPCAEnergyParametrization,1)  //TFCSPCAEnergyParametrization
+  float m_total_energy_normalization{1};
+  
+  ClassDefOverride(TFCSPCAEnergyParametrization,2)  //TFCSPCAEnergyParametrization
  
 };
 

Simulation/ISF/ISF_FastCaloSim/ISF_FastCaloSimEvent/src/TFCSPCAEnergyParametrization.cxx

@@ -50,7 +50,7 @@ void TFCSPCAEnergyParametrization::Print(Option_t *option) const
   TFCSEnergyParametrization::Print(option);
   
   if(longprint) {
-    ATH_MSG(INFO) << optprint <<"  #bins="<<m_numberpcabins<<", layers=";
+    ATH_MSG(INFO) << optprint <<"  #bins="<<m_numberpcabins<<", Enorm="<<m_total_energy_normalization<<", layers=";
     for(unsigned int i=0;i<m_RelevantLayers.size();i++) {
       if(i>0) msg()<<", ";
       msg()<<m_RelevantLayers[i];
@@ -127,7 +127,7 @@ FCSReturnCode TFCSPCAEnergyParametrization::simulate(TFCSSimulationState& simuls
     simdata[l]*=scalefactor;
    }
    
-   double total_energy=simdata[layerNr.size()]*simulstate.E()/Ekin_nominal();
+   double total_energy=simdata[layerNr.size()]*simulstate.E()/m_total_energy_normalization;
    simulstate.set_E(total_energy);
    ATH_MSG_DEBUG("set E to total_energy="<<total_energy);