KVElasticCountRates.h

 
#ifndef __KVElasticCountRates_H
#define __KVElasticCountRates_H
 
#include "TVector3.h"
#include "KVNameValueList.h"
#include "KVPosition.h"
#include "TH2.h"
#include "KVNucleus.h"
#include "KV2Body.h"
#include "KVTarget.h"
#include <map>
 
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 {
 
    TH1F* fDepth;
    TH1F* fTheta;
    TH2F* fGlobalMap;
    Int_t fBinE;
 
    Double_t fEnergy;
    KVNucleus fProj;
 
    KVPosition fAngularRange;
 
    KV2Body* fKinematics;
 
    KVTarget* fTarget;
    Double_t fAtomicDensity;
    Double_t fVolume;
   Double_t xsec;
   Double_t theta, phi;
   Double_t fNtirages;
    Bool_t fMultiLayer;
 
    TVector3 fBeamDirection;
    Int_t fIntLayer;
 
    Double_t fExx;
 
    KVHashList fHistos;
 
   std::map<std::string, KVElasticCountRate> fRates;
 
public:
 
   KVElasticCountRates(Double_t theta_min = 0, Double_t theta_max = 180, Double_t phi_min = 0, Double_t phi_max = 360, const TVector3& beam_dir = {0, 0, 1});
   virtual ~ KVElasticCountRates();
 
   void SetRun(Int_t run);
   void SetProjectile(const Char_t* nuc, Double_t e_sur_a);
   void SetTargetScatteringLayer(const Char_t* name);
   void SetEbinning(Int_t nbins = 500);
 
   void CalculateScattering(Int_t N = 10000);
   void SetTargetExcitedState(Double_t ex)
   {
      fExx = ex;
   }
 
   TH1F* GetDepth()
   {
 
      return fDepth;
   }
   TH1F* GetTheta()
   {
 
      return fTheta;
   }
   Int_t GetEbinning(void)
   {
 
      return fBinE;
   }
 
   void FillHistograms(KVNameValueList*);
   const KVHashList& GetHistos() const
   {
      return fHistos;
   };
 
   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];
   }
 
   ClassDef(KVElasticCountRates, 1)//Calculate elastic scattering count rates in multidetector arrays
};
 
#endif