| 1 | -*-*- END -*-*-
|
|---|
| 2 | 2002/08/19: Thomas Bretz
|
|---|
| 3 |
|
|---|
| 4 | * MCascade.[h,cc]:
|
|---|
| 5 | - removed argument of DoInvCompton
|
|---|
| 6 | - delete GRandom only if set
|
|---|
| 7 | - set gRandom to 0 after deletion
|
|---|
| 8 |
|
|---|
| 9 | * MElectron.[h,cc]:
|
|---|
| 10 | - added Sigma_ge
|
|---|
| 11 | - removed argument from DoInvCompton
|
|---|
| 12 | - Implemented Omega_sigmae and RandomThetaE
|
|---|
| 13 | - Changed old wrong dertermination of interaction angle to
|
|---|
| 14 | correct simulation.
|
|---|
| 15 | - Reengenieered some formulas to be more correct in terms of numerics
|
|---|
| 16 |
|
|---|
| 17 |
|
|---|
| 18 |
|
|---|
| 19 | 2002/08/02: Thomas Bretz
|
|---|
| 20 |
|
|---|
| 21 | * MCascade.[h,cc]:
|
|---|
| 22 | - doesn't overwrite existing files anymore
|
|---|
| 23 | - removed redefinition of default argument
|
|---|
| 24 |
|
|---|
| 25 | * MFit.cc:
|
|---|
| 26 | - removed redefinition of default argument
|
|---|
| 27 |
|
|---|
| 28 | * analp.C:
|
|---|
| 29 | - added smearing
|
|---|
| 30 |
|
|---|
| 31 |
|
|---|
| 32 |
|
|---|
| 33 | 2002/07/29: Thomas Bretz
|
|---|
| 34 |
|
|---|
| 35 | * MCascade.[h,cc]:
|
|---|
| 36 | - replaces draw by fIsBatch, using gROOT->IsBatch in addition.
|
|---|
| 37 |
|
|---|
| 38 | * MFit.h:
|
|---|
| 39 | - added default=0 to SetParameter for min and max
|
|---|
| 40 |
|
|---|
| 41 |
|
|---|
| 42 |
|
|---|
| 43 | 2002/07/29: Thomas Bretz
|
|---|
| 44 |
|
|---|
| 45 | * MElectron.[h,cc]:
|
|---|
| 46 | - if the photon energy is to small return a straight line in
|
|---|
| 47 | Compton
|
|---|
| 48 | - changed lolimit of p_e integration -0.5
|
|---|
| 49 | - changed number of calculated points to 50 for fP, fQ
|
|---|
| 50 | - changed usage of TRandom to gRandom
|
|---|
| 51 | - changed some TVector and TMatrix objects to static
|
|---|
| 52 | - added new copy constructors
|
|---|
| 53 |
|
|---|
| 54 | * MPairProduction.cc:
|
|---|
| 55 | - use the new copy constructor of MElectron
|
|---|
| 56 |
|
|---|
| 57 | * MParticle.[h,cc]:
|
|---|
| 58 | - added new copy constructors
|
|---|
| 59 | - Added a TAxis object in Planck
|
|---|
| 60 | - replaced usage of TRandom by gRandom
|
|---|
| 61 | - replaced TVector, TMatrix by TVector3, TRotation
|
|---|
| 62 | - fixed a bug in SetNewPosition (fX was not incremented if the
|
|---|
| 63 | particle reached the observers plain)
|
|---|
| 64 | - incremented version-number
|
|---|
| 65 |
|
|---|
| 66 | * MPhoton.[h,cc]:
|
|---|
| 67 | - added Fill member function
|
|---|
| 68 | - added new copy constructors
|
|---|
| 69 |
|
|---|
| 70 | * MFit.[h,cc], MCascade.[h,cc]:
|
|---|
| 71 | - added
|
|---|
| 72 |
|
|---|
| 73 | * phys.C:
|
|---|
| 74 | - moved code to MCascade
|
|---|
| 75 |
|
|---|
| 76 | * Makefile, PhysLinkDef.h:
|
|---|
| 77 | - added MFit, MCascade
|
|---|
| 78 |
|
|---|
| 79 | * analp.C:
|
|---|
| 80 | - simplified
|
|---|
| 81 | - added fit for cutoff
|
|---|
| 82 |
|
|---|
| 83 |
|
|---|
| 84 |
|
|---|
| 85 | 2002/07/26: Thomas Bretz
|
|---|
| 86 |
|
|---|
| 87 | * MParticle.[h,cc]:
|
|---|
| 88 | - added fX data member
|
|---|
| 89 | - corrected a small bug in the SetNewPosition calculation
|
|---|
| 90 | (dr*cos(fTheta)>R*cos(fTheta) --> x(2)>R*cos(fTheta))
|
|---|
| 91 | - added usage of fX to the =-operator and SetNewPosition
|
|---|
| 92 | - added fSrcR (unsaved)
|
|---|
| 93 |
|
|---|
| 94 | * MElectron.cc:
|
|---|
| 95 | - added usage of fX to the =-operator and SetNewPosition
|
|---|
| 96 | - changed to use the more accurate SetNewPosition in case B==0
|
|---|
| 97 |
|
|---|
| 98 | * MPhoton.cc:
|
|---|
| 99 | - return 1e50 as InteractionLength in case Eg<100GeV
|
|---|
| 100 |
|
|---|
| 101 | * phys.C:
|
|---|
| 102 | - added support for a B-Field-Bubble
|
|---|
| 103 |
|
|---|
| 104 |
|
|---|
| 105 |
|
|---|
| 106 | 2002/07/21: Thomas Bretz
|
|---|
| 107 |
|
|---|
| 108 | * MElectron.cc:
|
|---|
| 109 | - changed upper limit in InteractionLength from 1e-11 to 2e-11
|
|---|
| 110 | - changed upper limit of p_e in EnergyLoss from -10.8 to -10.6
|
|---|
| 111 |
|
|---|
| 112 | * MParticle.[h,cc]:
|
|---|
| 113 | - replaced Planck by tanjas background
|
|---|
| 114 |
|
|---|
| 115 | * background.txt:
|
|---|
| 116 | - added
|
|---|
| 117 |
|
|---|
| 118 | * MPhoton.[h,cc]:
|
|---|
| 119 | - changed lower limit of Int2-Integration for Interaction Length
|
|---|
| 120 | in case Eg<5e4 to fit better its needs
|
|---|
| 121 | - changed max/min of drawing interaction length
|
|---|
| 122 | - changed lower limit of interaction length plot
|
|---|
| 123 |
|
|---|
| 124 | * analp.C:
|
|---|
| 125 | - added plotting 1h1426 data
|
|---|
| 126 | - changed alpha/plot from -2/2 to -1/1 (which means I do the well
|
|---|
| 127 | known E^2 plot of an E^-2 spectrum)
|
|---|
| 128 | - added a cutoff to simulate cutoff spectra
|
|---|
| 129 | - some small corrections for the layout of the plots
|
|---|
| 130 |
|
|---|
| 131 | * phys.C:
|
|---|
| 132 | - added a maximum number of inv. Compton processes
|
|---|
| 133 | - don't throw away photons which were reemitted lower than the lower
|
|---|
| 134 | limit of the plot (so they are still written to the file)
|
|---|
| 135 | - changed the lower limit of the PhotonSpect simulation function
|
|---|
| 136 | to a physically more realistic value. No photon which don't have
|
|---|
| 137 | an corresponding angle are produced anymore.
|
|---|
| 138 | - stop the electron only if its energy is less than a thousand
|
|---|
| 139 | of its primary energy
|
|---|
| 140 |
|
|---|
| 141 |
|
|---|
| 142 |
|
|---|
| 143 | 2002/06/21: Thomas Bretz
|
|---|
| 144 |
|
|---|
| 145 | * analp.C:
|
|---|
| 146 | - clean up the code
|
|---|
| 147 | - added the point-spread histograms from different directions
|
|---|
| 148 |
|
|---|
| 149 | * mphys.C:
|
|---|
| 150 | - removed the angular histograms
|
|---|
| 151 | - added the plot spectral index
|
|---|
| 152 | - added the timer to get a usable ubdated window
|
|---|
| 153 |
|
|---|
| 154 |
|
|---|
| 155 |
|
|---|
| 156 | 2002/06/20: Thomas Bretz
|
|---|
| 157 |
|
|---|
| 158 | * MElectron.cc:
|
|---|
| 159 | - Some changes to the use of beta (simplifications)
|
|---|
| 160 | - changed the default of rho (B==0) to 0
|
|---|
| 161 |
|
|---|
| 162 | * MParticle.cc:
|
|---|
| 163 | - made TRandom static (InitRandom)
|
|---|
| 164 |
|
|---|
| 165 | * mphys.C:
|
|---|
| 166 | - added support for Emis and Esum
|
|---|
| 167 |
|
|---|
| 168 |
|
|---|
| 169 |
|
|---|
| 170 | 2002/06/17: Thomas Bretz
|
|---|
| 171 |
|
|---|
| 172 | * MElectron.[h,cc]:
|
|---|
| 173 | - Added SetNewPositionB
|
|---|
| 174 | - changed the constructor to support positrons
|
|---|
| 175 |
|
|---|
| 176 | * MPairProduction.cc:
|
|---|
| 177 | - use new constructor of MElectron
|
|---|
| 178 |
|
|---|
| 179 | * MParticle.[h,cc]:
|
|---|
| 180 | - added InitRandom
|
|---|
| 181 | - added GetType
|
|---|
| 182 |
|
|---|
| 183 | * anale.C:
|
|---|
| 184 | - added SetDirectory(NULL)
|
|---|
| 185 |
|
|---|
| 186 | * analp.C:
|
|---|
| 187 | - changed to an automatic scaling
|
|---|
| 188 |
|
|---|
| 189 | * phys.C:
|
|---|
| 190 | - changed limit to 60min
|
|---|
| 191 | - added lower limit for electrons 5*lo
|
|---|
| 192 | - Changed to use a B-field
|
|---|
| 193 |
|
|---|
| 194 |
|
|---|
| 195 |
|
|---|
| 196 | 2002/06/17: Thomas Bretz
|
|---|
| 197 |
|
|---|
| 198 | * MElectron.cc:
|
|---|
| 199 | - Added SetNewDirection for the electron itself
|
|---|
| 200 |
|
|---|
| 201 | * MParticle.cc:
|
|---|
| 202 | - changed sqrt() to asin(sqrt()) - no effect (but correct)
|
|---|
| 203 |
|
|---|
| 204 |
|
|---|
| 205 |
|
|---|
| 206 | 2002/06/14: Thomas Bretz
|
|---|
| 207 |
|
|---|
| 208 | * MElectron.[h,cc]:
|
|---|
| 209 | - Moved Planck function to MParticle
|
|---|
| 210 | - Added the DiSum function
|
|---|
| 211 | - changed Li2 to use the DiSum function (speed reasons)
|
|---|
| 212 | - changed the eps of all integrals to 1e-2
|
|---|
| 213 | - changed the p_e function to use the Compton function
|
|---|
| 214 | - made all function which are used for integration only static
|
|---|
| 215 | and changed the function definition from=0, to=0
|
|---|
| 216 |
|
|---|
| 217 | * MPhoton.[h,cc]:
|
|---|
| 218 | - removed the planck function
|
|---|
| 219 | - changed the eps of all integrals to 1e-2
|
|---|
| 220 |
|
|---|
| 221 | * MParticle.[h,cc]:
|
|---|
| 222 | - added the planck function
|
|---|
| 223 |
|
|---|
| 224 | * Makefile:
|
|---|
| 225 | - removed unused source files
|
|---|
| 226 |
|
|---|
| 227 | * analp.C:
|
|---|
| 228 | - did some small fixes
|
|---|
| 229 | - changed the sclae of the radius- and phi-plots
|
|---|
| 230 |
|
|---|
| 231 | * phys.C:
|
|---|
| 232 | - changed the integral eps to 1e-2
|
|---|
| 233 | - small changed to the output
|
|---|
| 234 | - removed integral==0, replaced by Ep==pow(10,0)
|
|---|
| 235 |
|
|---|
| 236 |
|
|---|
| 237 |
|
|---|
| 238 |
|
|---|
| 239 | 2002/06/13: Thomas Bretz
|
|---|
| 240 |
|
|---|
| 241 | * MParticle.h:
|
|---|
| 242 | - Implemented Set/IsPrimary
|
|---|
| 243 | - removed pure abstract method (root can't read this)
|
|---|
| 244 | - added default constructor
|
|---|
| 245 |
|
|---|
| 246 | * MElectron.cc:
|
|---|
| 247 | - changed many variables to const
|
|---|
| 248 | - changed the output
|
|---|
| 249 |
|
|---|
| 250 | * MPairProduction.cc:
|
|---|
| 251 | - changed many variables to const
|
|---|
| 252 |
|
|---|
| 253 | * MPhoton.[h,cc]:
|
|---|
| 254 | - changed the upper and lower limit of the integration in Int2
|
|---|
| 255 | - added upper and lower limit of draw to DrawInteractionLength
|
|---|
| 256 |
|
|---|
| 257 | * energyloss.C:
|
|---|
| 258 | - changed output range for ZofR
|
|---|
| 259 |
|
|---|
| 260 | * phys.C:
|
|---|
| 261 | - changed energy randomizer
|
|---|
| 262 | - added root file output
|
|---|
| 263 | - removed log scale from point spread plots
|
|---|
| 264 |
|
|---|
| 265 | * analp.C:
|
|---|
| 266 | - added macro to analize the photons from the root file
|
|---|
| 267 |
|
|---|
| 268 | * anale.C:
|
|---|
| 269 | - added macro to analize the electrons from the root file
|
|---|
| 270 |
|
|---|
| 271 |
|
|---|
| 272 |
|
|---|
| 273 | 2002/06/12: Thomas Bretz
|
|---|
| 274 |
|
|---|
| 275 | * MElectron.[h,cc]:
|
|---|
| 276 | - added a primitive theta dependancy to DoInvCompton
|
|---|
| 277 | - added DrawInteractionLength
|
|---|
| 278 | - changed algorithm in ZofR and RofZ
|
|---|
| 279 |
|
|---|
| 280 | * MParticle.[h,cc]:
|
|---|
| 281 | - added RofZ and ZofR as static functions
|
|---|
| 282 | - changed from MParContainer to TObject
|
|---|
| 283 | - removed unnecessary or commented inline functions
|
|---|
| 284 |
|
|---|
| 285 | * MPhoton.[h,cc]:
|
|---|
| 286 | - added DrawInteractionLength
|
|---|
| 287 | - fixed the bug causing the 'strange factor': sigma_gg needs
|
|---|
| 288 | the sqrt of s.
|
|---|
| 289 |
|
|---|
| 290 | * energyloss.C:
|
|---|
| 291 | - added
|
|---|
| 292 |
|
|---|
| 293 | * MPairProduction.[h,cc]:
|
|---|
| 294 | - simplified formula
|
|---|
| 295 | - changed arguments to Energy of Photon and interaction angle
|
|---|
| 296 |
|
|---|
| 297 | * phys.C:
|
|---|
| 298 | - added correct photon energy for pair production
|
|---|
| 299 | - added correct angle for pair production
|
|---|
| 300 | - changed output to E^2*Counts
|
|---|
| 301 | - changed gamma production from E^-2 to Uniform by weighting
|
|---|
| 302 | - removed all unecessary functions and code (s.energyloss.C)
|
|---|
| 303 |
|
|---|
| 304 |
|
|---|
| 305 |
|
|---|
| 306 | 2002/06/10: Thomas Bretz
|
|---|
| 307 |
|
|---|
| 308 | * MElectron.[h,cc], MPhoton.[h,cc], MParticle.[h,cc]:
|
|---|
| 309 | - implemented SetNewDirection, SetNewPosition
|
|---|
| 310 | - implemented DoInvCompton
|
|---|
| 311 | - =operator
|
|---|
| 312 |
|
|---|
| 313 | * MPairProduction.cc:
|
|---|
| 314 | - changed to use SetNewDirection Method
|
|---|
| 315 |
|
|---|
| 316 | * phys.C:
|
|---|
| 317 | - changed to new style functions
|
|---|
| 318 |
|
|---|
| 319 | * MGeinIRPhoton.cc, MGenPrimaryParticle.[h,cc]:
|
|---|
| 320 | - make compile (not used)
|
|---|
| 321 | - changed Process to get angle
|
|---|
| 322 |
|
|---|
| 323 |
|
|---|
