#include <CosmicsGenerator.h>

Inheritance diagram for CosmicsGenerator:

Collaboration diagram for CosmicsGenerator:

[legend]

Public Member Functions
	CosmicsGenerator ()

	~CosmicsGenerator () override

Bool_t	ReadEvent (FairPrimaryGenerator *) override

virtual Bool_t	Init (Bool_t largeMom)

Bool_t	Init (const char *inFile) override

Bool_t	Init (const char *inFile, int startEvent) override

Public Member Functions inherited from SHiP::Generator
	Generator ()

virtual	~Generator ()

virtual Bool_t	Init (const char *, int)=0

virtual Bool_t	Init (const char *)=0

virtual Bool_t	Init (const std::vector< std::string > &inFiles, int startNumber)

virtual Bool_t	Init (const std::vector< std::string > &inFiles)

virtual void	UseExternalFile (std::string x, Int_t i)

virtual void	UseExternalFile (std::vector< std::string > &inFiles, Int_t i)

Public Attributes
double	z0

double	yBox

double	xBox

double	zBox

double	xdist

double	zdist

double	minE

int	n_EVENTS

Protected Member Functions
void	GenerateDynamics ()

Bool_t	DetectorBox ()

Protected Attributes
double	P

double	px

double	py

double	pz

double	x

double	y

double	z

double	weighttest

double	weight

double	mass

double	FluxIntegral

double	theta

int	PID

int	nInside

int	nTest

Bool_t	high

Protected Attributes inherited from SHiP::Generator
std::optional< std::string >	fextFile

Int_t	firstEvent = 0

Private Attributes
Co3Rng *	fRandomEngine

Detailed Description

Definition at line 47 of file CosmicsGenerator.h.

Constructor & Destructor Documentation

◆ CosmicsGenerator()

CosmicsGenerator::CosmicsGenerator ( )

inline

constructor,destructor

Definition at line 50 of file CosmicsGenerator.h.

50{};

◆ ~CosmicsGenerator()

CosmicsGenerator::~CosmicsGenerator ( )

inlineoverride

Definition at line 51 of file CosmicsGenerator.h.

                               {
    delete fRandomEngine;
    std::cout << nInside << " events have been generated." << std::endl;
    std::cout << "There is a total of " << nInside << "/" << nTest
              << " muons that passed close enough to the detector."
              << std::endl;
    std::cout << "Including the given weight this corresponds to ";
    std::cout << weighttest / FluxIntegral / xdist / zdist * 10000
              << " spills (1 spill = " << xdist * zdist * FluxIntegral / 10000;
    std::cout << " real cosmic muons = " << n_EVENTS << " simulated events)."
              << std::endl;
  };

Member Function Documentation

◆ DetectorBox()

Bool_t CosmicsGenerator::DetectorBox ( )

protected

Definition at line 59 of file CosmicsGenerator.cxx.

                                     {
  // check, if a given staring setup x,y,z,px,py,pz will lead to a
  // close enough to the detector
 
  return (TMath::Abs(x - (y + yBox) * px / py) < xBox &&
          TMath::Abs(z - z0 - (y + yBox) * pz / py) < zBox) ||
         (TMath::Abs(x - (y - yBox) * px / py) < xBox &&
          TMath::Abs(z - z0 - (y - yBox) * pz / py) < zBox) ||
         (TMath::Abs(y - (x + xBox) * py / px) < yBox &&
          TMath::Abs(z - z0 - (x + xBox) * pz / px) < zBox) ||
         (TMath::Abs(y - (x - xBox) * py / px) < yBox &&
          TMath::Abs(z - z0 - (x - xBox) * pz / px) < zBox);
}

◆ GenerateDynamics()

void CosmicsGenerator::GenerateDynamics ( )

protected

Definition at line 73 of file CosmicsGenerator.cxx.

                                        {
  // generate starting conditions for CM until the DetectorBox is hit
  do {
    weighttest += weight;
    nTest++;  // book keeping
    // momentum components
    double phi = fRandomEngine->Uniform(0, 2 * TMath::Pi());
    theta = fRandomEngine->fTheta->GetRandom();
    px = TMath::Sin(phi) * TMath::Sin(theta);
    pz = TMath::Cos(phi) * TMath::Sin(theta);
    py = -TMath::Cos(theta);
    // staring location
    x = fRandomEngine->Uniform(-xdist / 2, xdist / 2);
    z = fRandomEngine->Uniform(z0 - zdist / 2, z0 + zdist / 2);
  } while (!DetectorBox());
  nInside++;
}

◆ Init() [1/3]

Bool_t CosmicsGenerator::Init ( Bool_t largeMom )

virtual

Definition at line 91 of file CosmicsGenerator.cxx.

                                             {
  // general
  fRandomEngine = new Co3Rng();
  TDatabasePDG* pdgBase = TDatabasePDG::Instance();
  mass = pdgBase->GetParticle(13)->Mass();  // muons!
  cout << "--------------------------------------------------------------------"
          "--"
       << endl;
  cout << "configuration for the CMBG as defined in "
          "$FAIRSHIP/python/CMBG_conf.py: "
       << endl;
  cout << "x_dist:  " << xdist << endl;
  cout << "z_dist:  " << zdist << endl;
  cout << "x_box:   " << xBox << endl;
  cout << "y_box:   " << yBox << endl;
  cout << "z_box:   " << zBox << endl;
  cout << "n_EVENTS:" << n_EVENTS << endl;
  cout << "minE:    " << minE << endl << endl;
  if (xdist * zdist * n_EVENTS == 0) {
    cout << "check the configuration for unphysical behavior." << endl
         << "We stop the execution." << endl
         << endl;
    return kFALSE;
  }
 
  high = largeMom;
  if (high)
    cout << "Simulation for high momentum" << endl;
  else
    cout << "Simulation for low momentum" << endl;
 
  // calculating weights for this run
  // weight_flux: expected #muons per spill/ #simulated events per spill:
  // FluxIntegral*xdist*zdist/EVENTS;
  //              the respective integrals are calculated from the fluxes
  // weight_DetectorBox: only consider CM hitting the DetectorBox
  //                     this is gained from a MC test of 10xEVENTS events
  double weight_flux, weight_DetectorBox;
  FluxIntegral = 0;
  if (!high) {  // momentum range 1 GeV - 100 GeV
    if (minE > 100) {
      cout << "choose minE < 100 !" << endl;
      return kFALSE;
    }
    double dt = TMath::Pi() / 2 / 100;
    for (double t = dt / 2; t < TMath::Pi() / 2; t += dt) {
      FluxIntegral += fRandomEngine->fSpectrumL(t, minE, 0);
    }
    FluxIntegral = 2 * TMath::Pi() / 3 * FluxIntegral * dt * 10000;
    cout << "LowE CM flux with P < minE = " << minE << " : " << FluxIntegral
         << "m-2s-1" << endl;
  } else {  // momentum range 100 GeV - 1000 GeV
    FluxIntegral =
        2 * TMath::Pi() / 3 * fRandomEngine->fSpectrumH->Integral(100, 1000);
    cout << "High E CM flux: " << FluxIntegral << "m-2s-1" << endl;
  }
  weight_flux = FluxIntegral * xdist * zdist / n_EVENTS / 10000;
  nInside = 0;
  nTest = 0;
  weighttest = 0;  // book keeping
  y = 1900;        // all muons start 19m over beam axis
  for (; nInside < 10 * n_EVENTS;) {
    GenerateDynamics();
  }
  weight_DetectorBox = 0.10 * nTest / n_EVENTS;
  weight = weight_flux / weight_DetectorBox;
  cout << "weight_DetectorBox: " << weight_DetectorBox << ", weight: " << weight
       << endl;
  cout << "--------------------------------------------------------------------"
          "--"
       << endl
       << endl;
  nInside = 0;
  nTest = 0;
  weighttest = 0;  // book keeping
 
  return kTRUE;
}

◆ Init() [2/3]

Bool_t CosmicsGenerator::Init ( const char * inFile )

inlineoverridevirtual

Implements SHiP::Generator.

Definition at line 69 of file CosmicsGenerator.h.

69{ return Init(inFile, 0); };

CosmicsGenerator::Init

virtual Bool_t Init(Bool_t largeMom)

Definition: CosmicsGenerator.cxx:91

◆ Init() [3/3]

Bool_t CosmicsGenerator::Init	(	const char *	inFile,
		int	startEvent
	)

inlineoverridevirtual

Implements SHiP::Generator.

Definition at line 71 of file CosmicsGenerator.h.

                                                           {
    LOG(warning) << "Init with files not implemented for CosmicsGenerator.";
    return false;
  };

◆ ReadEvent()

Bool_t CosmicsGenerator::ReadEvent ( FairPrimaryGenerator * cpg )

override

Definition at line 170 of file CosmicsGenerator.cxx.

                                                            {
  // muon or anti-muon
  PID = 26 * (fRandomEngine->Uniform(0, 1) < 1.0 / 2.278) - 13;
  // starting conditions
  GenerateDynamics();
  // momentum in the two regions, < or > 100 GeV
  if (!high)
    P = fRandomEngine->fSpectrumL(theta, minE);
  else
    P = fRandomEngine->fSpectrumH->GetRandom();
  px = px * P;
  py = py * P;
  pz = pz * P;
  // transfer to Geant4
  cpg->AddTrack(
      PID, px, py, pz, x, y, z, -1, true, TMath::Sqrt(P * P + mass * mass), 0,
      weight);  // -1 = Mother ID, true = tracking, SQRT(x) = Energy, 0 = t
  if (!nInside % 10000) {
    cout << nInside / 10000 << "10k events have been simulated" << endl;
  }
  return kTRUE;
}

Member Data Documentation

◆ FluxIntegral

double CosmicsGenerator::FluxIntegral

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ fRandomEngine

Co3Rng* CosmicsGenerator::fRandomEngine

private

Definition at line 80 of file CosmicsGenerator.h.

◆ high

Bool_t CosmicsGenerator::high

protected

Definition at line 85 of file CosmicsGenerator.h.

◆ mass

double CosmicsGenerator::mass

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ minE

double CosmicsGenerator::minE

Definition at line 76 of file CosmicsGenerator.h.

◆ n_EVENTS

int CosmicsGenerator::n_EVENTS

Definition at line 77 of file CosmicsGenerator.h.

◆ nInside

int CosmicsGenerator::nInside

protected

Definition at line 84 of file CosmicsGenerator.h.

◆ nTest

int CosmicsGenerator::nTest

protected

Definition at line 84 of file CosmicsGenerator.h.

◆ P

double CosmicsGenerator::P

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ PID

int CosmicsGenerator::PID

protected

Definition at line 84 of file CosmicsGenerator.h.

◆ px

double CosmicsGenerator::px

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ py

double CosmicsGenerator::py

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ pz

double CosmicsGenerator::pz

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ theta

double CosmicsGenerator::theta

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ weight

double CosmicsGenerator::weight

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ weighttest

double CosmicsGenerator::weighttest

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ x

double CosmicsGenerator::x

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ xBox

double CosmicsGenerator::xBox

Definition at line 76 of file CosmicsGenerator.h.

◆ xdist

double CosmicsGenerator::xdist

Definition at line 76 of file CosmicsGenerator.h.

◆ y

double CosmicsGenerator::y

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ yBox

double CosmicsGenerator::yBox

Definition at line 76 of file CosmicsGenerator.h.

◆ z

double CosmicsGenerator::z

protected

Definition at line 83 of file CosmicsGenerator.h.

◆ z0

double CosmicsGenerator::z0

Definition at line 76 of file CosmicsGenerator.h.

◆ zBox

double CosmicsGenerator::zBox

Definition at line 76 of file CosmicsGenerator.h.

◆ zdist

double CosmicsGenerator::zdist

Definition at line 76 of file CosmicsGenerator.h.

The documentation for this class was generated from the following files:

shipgen/CosmicsGenerator.h
shipgen/CosmicsGenerator.cxx

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CosmicsGenerator()

◆ ~CosmicsGenerator()

Member Function Documentation

◆ DetectorBox()

◆ GenerateDynamics()

◆ Init() [1/3]

◆ Init() [2/3]

◆ Init() [3/3]

◆ ReadEvent()

Member Data Documentation

◆ FluxIntegral

◆ fRandomEngine

◆ high

◆ mass

◆ minE

◆ n_EVENTS

◆ nInside

◆ nTest

◆ P

◆ PID

◆ px

◆ py

◆ pz

◆ theta

◆ weight

◆ weighttest

◆ x

◆ xBox

◆ xdist

◆ y

◆ yBox

◆ z

◆ z0

◆ zBox

◆ zdist