Creation of hidden photon flux code

[tomasoshea]

Created 2 new files TRestAxionSolarHiddenPhotonFlux.h and TRestAxionSolarHiddenPhotonFlux.cxx to go with the TRestAxionSolarFlux code.

[jgalan]

I think we do not need the mass here? It is the method GetRandomEnergyAndRadius where we define the mass. There should be something inside that method that gets the probability depending on the mass.

Then, inside the processing chain, TRestAxionGeneratorProcess will chose a given mass, and it will use that method to get a random particle from the sun, with given energy and position in the sun (for a given mass).

[tomasoshea]

I've been looking at it and instead of having the mass as an argument to all the functions that need it, could we not have it as a member of the class? So for example when GetRandomEnergyAndRadius is called with double_t mass as an argument we can then set HiddenPhotonMass = mass? So that since the flux depends on the mass we wouldn't have to supply mass as an argument to every function that deals with the flux.

[jgalan]

Yes, this is possible.

We just need to make sure then that we have this mass inside TRestAxionSolarFlux and not within TRestAxionQCDSolarFlux or TRestAxionHiddenPhotonFlux.

So, we define a TRestAxionSolarFlux::fMass.

Then, inside the TRestAxionGeneratorProcess we can call TRestAxionSolarFlux::SetMass(m) before calling GetRandomEnergyAndRadius.

TRestAxionGeneratorProcess will only have access to those methods available to TRestAxionSolarFlux.

[jgalan]

Then, inside TRestAxionSolarFlux we should add:

Double_t fMass = 0; //!

Double_t GetMass() { return fMass; }
void SetMass(const Double_t &m ){ fMass = m; }

[tomasoshea]

Ok, then would we remove the double_t mass from the other functions and just have them use fMass? For example in GetTotalFlux(Double_t mass = 0)