TargetCount.h

Go to the documentation of this file.

#pragma once

#include "TGeoManager.h"
#include "TVector3.h"
#include "TH1D.h"

namespace ana
{
  /// Count number of targets within the main part of the ND (vSuperBlockND)
  class TargetCount
  {
  public:
    /// \param Z Count instances of this element. -1 = all elements
    TargetCount(int Z = -1);
      
    /// Count the number of nuclear targets of atomic number Z, and
    /// inside the fiducial bounds specified by min and max
    /// The method works by random-polling and is based on
    /// calculateMassesLong function in Geometry service
    /// \param Z Count instances of this element. -1 = all elements
    /// \param npoints is the number of samples in the random poll
    /// This value should be 1e6 or more.
    TargetCount(const TVector3& min,
                const TVector3& max,
                const int npoints,
                int Z = -1);

    /// Mass in kg
    double Mass() const {return fMass;}
    /// Number of nucleons (mass * avogadro's number)
    double NNucleons() const {return fNNucleons;}
    /// Number of nuclei, optionally restricted to a particular element or isotope
    double NNuclei(int Z=-1, int A=-1) const;
    /// Fractional error on mass or number of nucleons
    double ErrorFraction() const {return fErrorFrac;}
    /// Map of mass for each Z
    std::map<int, double> ZMassMap() const {return fZMassMap;}

  protected:
    int SetupGeoManager();
    TGeoManager* fGeoManager;
    double fMass;
    double fErrorFrac;
    double fNNucleons;

    std::map<int, double> fZMassMap;
    std::map<int, std::map<int, double>>  fZIsotopeMassMap;
    std::map<std::string, double> kMassError;

    double CountMass(int Z, TVector3 min , TVector3 max, int npoints, double& errofrac);
    double CountMass(const TGeoVolume* vol, int Z, double& errorfrac);

  
  };
}