Add pmtEdepFromFCS_StepInfo method for use by TileFCSmStepToTileHitVec

The `TileFCSmStepToTileHitVec` used previously public static methods of
`TileGeoG4SDCalc`.  This is no longer possible now that `TileGeoG4SDCalc`
is a service accessed via the `ITileCalculator` interface. As a workaround
the (now private) static method calls have been wrapped up in a new
public member function `pmtEdepFromFCS_StepInfo`.


Former-commit-id: 982d059a

TileCalorimeter/TileG4/TileG4Interfaces/TileG4Interfaces/ITileCalculator.h

@@ -85,6 +85,9 @@ class ITileCalculator : virtual public IService {
   virtual TileMicroHit GetTileMicroHit(const G4Step*, TileHitData& hitData) const = 0;
   ///
   virtual TileGeoG4LookupBuilder* GetLookupBuilder() const = 0;
+  /// Method used by TileFastCaloSim/TileFCSmStepToTileHitVec
+  virtual void pmtEdepFromFCS_StepInfo(TileHitData& hitData, double ene, double yLocal, double halfYLocal, double zLocal, int Ushape) const = 0;
+
 };
 //class ITileCalculator
 

TileCalorimeter/TileG4/TileGeoG4SD/src/TileGeoG4SDCalc.cc

@@ -881,6 +881,42 @@ G4double TileGeoG4SDCalc::BirkLaw(const G4Step* aStep) const
   }
 }
 
+void TileGeoG4SDCalc::pmtEdepFromFCS_StepInfo(TileHitData& hitData, double ene, double yLocal, double halfYLocal, double zLocal, int Ushape) const
+{
+  /// Code moved here from TileFastCaloSim/src/TileFCSmStepToTileHitVec.cxx
+  switch (Ushape) {
+  case 1:
+    hitData.edep_down = ene * this->Tile_1D_profile(hitData.tileSize, yLocal, zLocal, 0, hitData.nDetector, hitData.nSide);
+    hitData.edep_up = ene * this->Tile_1D_profile(hitData.tileSize, yLocal, zLocal, 1, hitData.nDetector, hitData.nSide);
+    break;
+
+  case 14:
+    hitData.edep_down = ene * this->Tile_1D_profileRescaled(hitData.tileSize, yLocal, zLocal, 0, hitData.nDetector, hitData.nSide);
+    hitData.edep_up = ene * this->Tile_1D_profileRescaled(hitData.tileSize, yLocal, zLocal, 1, hitData.nDetector, hitData.nSide);
+    break;
+
+  case 13:
+    hitData.edep_down = ene * this->Tile_1D_profileFunc(hitData.tileSize, yLocal, zLocal, 0, hitData.nDetector, hitData.nSide);
+    hitData.edep_up = ene * this->Tile_1D_profileFunc(hitData.tileSize, yLocal, zLocal, 1, hitData.nDetector, hitData.nSide);
+    break;
+
+  case 12:
+    hitData.edep_down = ene * this->Tile_1D_profileSym(hitData.tileSize, yLocal, zLocal, 0, hitData.nDetector, hitData.nSide);
+    hitData.edep_up = ene * this->Tile_1D_profileSym(hitData.tileSize, yLocal, zLocal, 1, hitData.nDetector, hitData.nSide);
+    break;
+
+  case 11:
+    hitData.edep_down = ene * this->Tile_1D_profileAsym(hitData.tileSize, yLocal, zLocal, 0, hitData.nDetector, hitData.nSide);
+    hitData.edep_up = ene * this->Tile_1D_profileAsym(hitData.tileSize, yLocal, zLocal, 1, hitData.nDetector, hitData.nSide);
+    break;
+
+  default:
+    hitData.edep_down = ene * ( 0.5 - 0.2 * yLocal/halfYLocal ); // should get 0.7 for yLocal = -halfYLocal
+    hitData.edep_up = ene - hitData.edep_down;
+  }
+  return;
+}
+
 G4double TileGeoG4SDCalc::Tile_1D_profileAsym(int row, G4double x, G4double y, int PMT, int nDetector, int) {
 
   const double xlow = -0.0505; //dPhi low [rad]

TileCalorimeter/TileG4/TileGeoG4SD/src/TileGeoG4SDCalc.hh

@@ -69,6 +69,8 @@ public:
   {
     return m_lookup;
   }
+  /// Method used by TileFastCaloSim/TileFCSmStepToTileHitVec
+  virtual void pmtEdepFromFCS_StepInfo(TileHitData& hitData, double ene, double yLocal, double halfYLocal, double zLocal, int Ushape) const override final;
 private:
   void CreateScintHit(int pmt, TileHitData& hitData) const;
   void UpdateScintHit(int pmt, TileHitData& hitData) const;