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 <micromegas/MicromegasTile.h>
#include <phparameter/PHParameterInterface.h>

#include <fun4all/SubsysReco.h>

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

class CylinderGeomMicromegas;
class PHCompositeNode;
class PHG4Hit;
class TVector3;

class PHG4MicromegasHitReco : public SubsysReco, public PHParameterInterface
{

  public:
  explicit PHG4MicromegasHitReco(
    const std::string &name = "PHG4MicromegasHitReco",
    const std::string &detector = "MICROMEGAS");

  int InitRun(PHCompositeNode*) override;

  int process_event(PHCompositeNode*) override;

  void SetDefaultParameters() override;

  void set_tiles( const MicromegasTile::List& tiles )
  { m_tiles = tiles; }

  private:

  std::string bare_geonodename() const
  { return "CYLINDERGEOM_" + m_detector; }

  std::string full_geonodename() const
  { return "CYLINDERGEOM_" + m_detector + "_FULL"; }

  void setup_tiles(PHCompositeNode*);


  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 TVector3& local_position, double sigma ) const;

  std::string m_detector;

  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;

  bool m_zigzag_strips = true;

  MicromegasTile::List m_tiles;

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


  std::unique_ptr<gsl_rng, Deleter> m_rng;

};

#endif