PHG4MicromegasHitReco.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PHG4MicromegasHitReco.h

// Tell emacs that this is a C++ source
// -*- C++ -*-.

#ifndef G4MICROMEGAS_PHG4MICROMEGASHITRECO_H
#define G4MICROMEGAS_PHG4MICROMEGASHITRECO_H

#include <phparameter/PHParameterInterface.h>

#include <fun4all/SubsysReco.h>

#include <gsl/gsl_rng.h>
#include <memory>
#include <string>
#include <utility>
#include <vector>

class ActsGeometry;

class CylinderGeomMicromegas;
class PHCompositeNode;
class PHG4Hit;
class TVector2;

class PHG4MicromegasHitReco : public SubsysReco, public PHParameterInterface
{
 public:
  explicit PHG4MicromegasHitReco(const std::string& name = "PHG4MicromegasHitReco");

  int InitRun(PHCompositeNode*) override;

  int process_event(PHCompositeNode*) override;

  void SetDefaultParameters() override;

 private:
  std::string full_geonodename() const
  {
    return "CYLINDERGEOM_MICROMEGAS_FULL";
  }


  uint get_primary_electrons(PHG4Hit*) const;

  uint get_single_electron_amplification() const;

  using charge_pair_t = std::pair<int, double>;

  using charge_list_t = std::vector<charge_pair_t>;

  charge_list_t distribute_charge(CylinderGeomMicromegas*, uint tileid, const TVector2& local_position, double sigma) const;

  ActsGeometry* m_acts_geometry = nullptr;

  double m_tmin = -20;

  double m_tmax = 800;

  double m_electrons_per_gev = 0;

  double m_gain = 0;

  double m_cloud_sigma = 0.04;

  double m_diffusion_trans = 0.03;


  double m_added_smear_sigma_z = 0;
  double m_added_smear_sigma_rphi = 0;

  class Deleter
  {
   public:
    void operator()(gsl_rng* rng) const { gsl_rng_free(rng); }
  };


  std::unique_ptr<gsl_rng, Deleter> m_rng;
};

#endif