KVElasticCountRates.h

 
#ifndef __KVElasticCountRates_H
#define __KVElasticCountRates_H
#include "KVElasticScatter.h"
#include "TVector3.h"
#include "KVNameValueList.h"
#include "KVPosition.h"
#include "TH2.h"
#include "KVNucleus.h"
#include "KV2Body.h"
#include "KVTarget.h"
#include <map>
#include "KVDetectionSimulator.h"
 
struct KVElasticCountRate {
   double count_rate;
   double mean_energy;
   double intXsec;      // integrated X-section (barn)
   double fluence;      // ions/sec./cm**2
   double dissipation;  // MeV/sec./cm**2
   KVElasticCountRate(double c = 0, double e = 0, double i = 0, double f = 0, double d = 0)
      : count_rate(c), mean_energy(e), intXsec(i), fluence(f), dissipation(d) {}
   virtual ~KVElasticCountRate() {}
 
   ClassDef(KVElasticCountRate, 0) //Elastic scattering rate information for detector
};
 
class KVElasticCountRates : public KVElasticScatter {
    TH2F* fGlobalMap;
    KVPosition fAngularRange;
   Double_t fVolume = 0;            //volume factor for MC integration
   Double_t theta, phi;
    KVHashList fArrayHistos;
   std::map<std::string, KVElasticCountRate> fRates;
   Double_t fBeamIntensity = 1.e+7;
 
   KVDetectionSimulator det_sim;
   bool initial_checks_and_reset() override;
   void reset_before_new_scattering() override;
   std::pair<double, double> get_random_angles_for_scattering(const KV2Body&) override;
   void detect_particle_fill_histograms(KVNucleus* ejectile, double theta, double phi, double xsec) override;
   void end_of_run() override;
public:
 
   KVElasticCountRates(int run)
      : KVElasticScatter{run}, det_sim{gMultiDetArray}
   {}
 
   void SetAngularRange(Double_t theta_min = 0, Double_t theta_max = 180, Double_t phi_min = 0, Double_t phi_max = 360);
   void SetBeamIntensity(Double_t bi)
   {
      fBeamIntensity = bi;
   }
   void FillHistograms(const KVNameValueList*, double theta, double phi, double xsec);
   const KVHashList& GetHistos() const
   {
      return fArrayHistos;
   }
 
   void PrintResults(Double_t beam_intensity = 1.e+07);
   KVNameValueList PutResultsInList(Double_t beam_intensity = 1.e+07);
 
   KVElasticCountRate GetDetector(const std::string& name)
   {
      return fRates[name];
   }
 
   ClassDefOverride(KVElasticCountRates, 1)//Calculate elastic scattering count rates in multidetector arrays
};
 
#endif