-*-*- END -*-*- 2002/08/19: Thomas Bretz * MCascade.[h,cc]: - removed argument of DoInvCompton - delete GRandom only if set - set gRandom to 0 after deletion * MElectron.[h,cc]: - added Sigma_ge - removed argument from DoInvCompton - Implemented Omega_sigmae and RandomThetaE - Changed old wrong dertermination of interaction angle to correct simulation. - Reengenieered some formulas to be more correct in terms of numerics 2002/08/02: Thomas Bretz * MCascade.[h,cc]: - doesn't overwrite existing files anymore - removed redefinition of default argument * MFit.cc: - removed redefinition of default argument * analp.C: - added smearing 2002/07/29: Thomas Bretz * MCascade.[h,cc]: - replaces draw by fIsBatch, using gROOT->IsBatch in addition. * MFit.h: - added default=0 to SetParameter for min and max 2002/07/29: Thomas Bretz * MElectron.[h,cc]: - if the photon energy is to small return a straight line in Compton - changed lolimit of p_e integration -0.5 - changed number of calculated points to 50 for fP, fQ - changed usage of TRandom to gRandom - changed some TVector and TMatrix objects to static - added new copy constructors * MPairProduction.cc: - use the new copy constructor of MElectron * MParticle.[h,cc]: - added new copy constructors - Added a TAxis object in Planck - replaced usage of TRandom by gRandom - replaced TVector, TMatrix by TVector3, TRotation - fixed a bug in SetNewPosition (fX was not incremented if the particle reached the observers plain) - incremented version-number * MPhoton.[h,cc]: - added Fill member function - added new copy constructors * MFit.[h,cc], MCascade.[h,cc]: - added * phys.C: - moved code to MCascade * Makefile, PhysLinkDef.h: - added MFit, MCascade * analp.C: - simplified - added fit for cutoff 2002/07/26: Thomas Bretz * MParticle.[h,cc]: - added fX data member - corrected a small bug in the SetNewPosition calculation (dr*cos(fTheta)>R*cos(fTheta) --> x(2)>R*cos(fTheta)) - added usage of fX to the =-operator and SetNewPosition - added fSrcR (unsaved) * MElectron.cc: - added usage of fX to the =-operator and SetNewPosition - changed to use the more accurate SetNewPosition in case B==0 * MPhoton.cc: - return 1e50 as InteractionLength in case Eg<100GeV * phys.C: - added support for a B-Field-Bubble 2002/07/21: Thomas Bretz * MElectron.cc: - changed upper limit in InteractionLength from 1e-11 to 2e-11 - changed upper limit of p_e in EnergyLoss from -10.8 to -10.6 * MParticle.[h,cc]: - replaced Planck by tanjas background * background.txt: - added * MPhoton.[h,cc]: - changed lower limit of Int2-Integration for Interaction Length in case Eg<5e4 to fit better its needs - changed max/min of drawing interaction length - changed lower limit of interaction length plot * analp.C: - added plotting 1h1426 data - changed alpha/plot from -2/2 to -1/1 (which means I do the well known E^2 plot of an E^-2 spectrum) - added a cutoff to simulate cutoff spectra - some small corrections for the layout of the plots * phys.C: - added a maximum number of inv. Compton processes - don't throw away photons which were reemitted lower than the lower limit of the plot (so they are still written to the file) - changed the lower limit of the PhotonSpect simulation function to a physically more realistic value. No photon which don't have an corresponding angle are produced anymore. - stop the electron only if its energy is less than a thousand of its primary energy 2002/06/21: Thomas Bretz * analp.C: - clean up the code - added the point-spread histograms from different directions * mphys.C: - removed the angular histograms - added the plot spectral index - added the timer to get a usable ubdated window 2002/06/20: Thomas Bretz * MElectron.cc: - Some changes to the use of beta (simplifications) - changed the default of rho (B==0) to 0 * MParticle.cc: - made TRandom static (InitRandom) * mphys.C: - added support for Emis and Esum 2002/06/17: Thomas Bretz * MElectron.[h,cc]: - Added SetNewPositionB - changed the constructor to support positrons * MPairProduction.cc: - use new constructor of MElectron * MParticle.[h,cc]: - added InitRandom - added GetType * anale.C: - added SetDirectory(NULL) * analp.C: - changed to an automatic scaling * phys.C: - changed limit to 60min - added lower limit for electrons 5*lo - Changed to use a B-field 2002/06/17: Thomas Bretz * MElectron.cc: - Added SetNewDirection for the electron itself * MParticle.cc: - changed sqrt() to asin(sqrt()) - no effect (but correct) 2002/06/14: Thomas Bretz * MElectron.[h,cc]: - Moved Planck function to MParticle - Added the DiSum function - changed Li2 to use the DiSum function (speed reasons) - changed the eps of all integrals to 1e-2 - changed the p_e function to use the Compton function - made all function which are used for integration only static and changed the function definition from=0, to=0 * MPhoton.[h,cc]: - removed the planck function - changed the eps of all integrals to 1e-2 * MParticle.[h,cc]: - added the planck function * Makefile: - removed unused source files * analp.C: - did some small fixes - changed the sclae of the radius- and phi-plots * phys.C: - changed the integral eps to 1e-2 - small changed to the output - removed integral==0, replaced by Ep==pow(10,0) 2002/06/13: Thomas Bretz * MParticle.h: - Implemented Set/IsPrimary - removed pure abstract method (root can't read this) - added default constructor * MElectron.cc: - changed many variables to const - changed the output * MPairProduction.cc: - changed many variables to const * MPhoton.[h,cc]: - changed the upper and lower limit of the integration in Int2 - added upper and lower limit of draw to DrawInteractionLength * energyloss.C: - changed output range for ZofR * phys.C: - changed energy randomizer - added root file output - removed log scale from point spread plots * analp.C: - added macro to analize the photons from the root file * anale.C: - added macro to analize the electrons from the root file 2002/06/12: Thomas Bretz * MElectron.[h,cc]: - added a primitive theta dependancy to DoInvCompton - added DrawInteractionLength - changed algorithm in ZofR and RofZ * MParticle.[h,cc]: - added RofZ and ZofR as static functions - changed from MParContainer to TObject - removed unnecessary or commented inline functions * MPhoton.[h,cc]: - added DrawInteractionLength - fixed the bug causing the 'strange factor': sigma_gg needs the sqrt of s. * energyloss.C: - added * MPairProduction.[h,cc]: - simplified formula - changed arguments to Energy of Photon and interaction angle * phys.C: - added correct photon energy for pair production - added correct angle for pair production - changed output to E^2*Counts - changed gamma production from E^-2 to Uniform by weighting - removed all unecessary functions and code (s.energyloss.C) 2002/06/10: Thomas Bretz * MElectron.[h,cc], MPhoton.[h,cc], MParticle.[h,cc]: - implemented SetNewDirection, SetNewPosition - implemented DoInvCompton - =operator * MPairProduction.cc: - changed to use SetNewDirection Method * phys.C: - changed to new style functions * MGeinIRPhoton.cc, MGenPrimaryParticle.[h,cc]: - make compile (not used) - changed Process to get angle