It is currently Thu Dec 05, 2024 4:11 am

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 25 posts ]  Go to page 1, 2  Next
Author Message
PostPosted: Sun Jun 17, 2007 12:47 pm 
Offline

Joined: Sun Feb 11, 2007 4:13 pm
Posts: 9
Location: Minneapolis, MN
I have taken your default GALDEF file, 599278, and made just a change to the energy limits. The default limits are 1e1 - 1e7 MeV/n.

I have tried 3e1 - 1e9 and 1e1 - 1e8 MeV/n, and gotten really nonsensical results.

I have also tried 1e3 - 1e9 and 1e2 - 1e8 MeV/n, and these produced good output that agrees with data when plotted on the same graph. The default limits also worked as expected.

Any idea why changing the energy limits would cause problems?

Thanks,
Taylor Childers
University of Minnesota


Top
 Profile  
 
 Post subject: energy limits
PostPosted: Mon Jun 18, 2007 12:59 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
I can see no obvious reason, will try it myself.

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject:
PostPosted: Mon Jun 18, 2007 6:46 am 
Offline

Joined: Sun Feb 11, 2007 4:13 pm
Posts: 9
Location: Minneapolis, MN
Here is a plot that I made of the data with the default file 599278 and the same file with a modified energy range 600003.

http://webusers.physics.umn.edu/~childers/BCRatioA.jpg

I've included the galdef file for you to look over.

http://webusers.physics.umn.edu/~childe ... 600003.txt

Thanks for you help,
Taylor


Top
 Profile  
 
 Post subject: B/C
PostPosted: Wed Jun 20, 2007 2:07 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
could you plot the B and C separately, to see which one is causing the problem
in the ratio

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject: B and C spectra
PostPosted: Wed Jun 20, 2007 6:58 am 
Offline

Joined: Sun Feb 11, 2007 4:13 pm
Posts: 9
Location: Minneapolis, MN
Here they are:

http://webusers.physics.umn.edu/~childe ... ctrumA.jpg

http://webusers.physics.umn.edu/~childe ... ctrumA.jpg


Top
 Profile  
 
 Post subject: B/C
PostPosted: Wed Jun 20, 2007 7:09 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
so B is OK
but C is going to zero at below 10 GeV
.. very odd

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject: test C
PostPosted: Wed Jun 20, 2007 9:38 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
I ran 600003 and C looks normal.
which deepens the mystery

The only thing I changed was Z up to 6 only , and no electrons etc, to make it go faster.
Maybe you can try that, just in case.

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject: New B and C plots
PostPosted: Thu Jun 21, 2007 9:21 am 
Offline

Joined: Sun Feb 11, 2007 4:13 pm
Posts: 9
Location: Minneapolis, MN
I did as you said, running 600003 without tracking electrons, positrons, etc. and also setting all "use_Z_*" = 0 for Z > 6. Here are the results, they still don't look good.


http://webusers.physics.umn.edu/~childe ... ctrumB.jpg

http://webusers.physics.umn.edu/~childe ... ctrumB.jpg


Top
 Profile  
 
 Post subject: puzzled
PostPosted: Thu Jun 21, 2007 9:48 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
indeed I did not expect any change, was just a check.
I admit I'm stumped.
Must be something obscure connected with the compiler or CPU. ? Unexpected.

Tomorrow I will work out how to investigate using the log file - maybe switching
on a debug option.

Meanwhile can you also plot some more - p, He etc to see if there's a similar effect?

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject: B/C
PostPosted: Sun Jun 24, 2007 11:26 pm 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
no progress yet.
I asked Igor to run it too.
Meanwhile can you see if it is only C or are any other nuclei affected?

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject:
PostPosted: Wed Jun 27, 2007 9:24 am 
Offline

Joined: Tue Nov 29, 2005 6:07 pm
Posts: 34
Location: Stanford
Sorry, I am very busy preparing for ICRC.
I will try to do it this week.

_________________
Igor Moskalenko
Stanford University


Top
 Profile  
 
 Post subject: Other Nuclei Plots
PostPosted: Thu Jun 28, 2007 4:30 pm 
Offline

Joined: Sun Feb 11, 2007 4:13 pm
Posts: 9
Location: Minneapolis, MN
Here are those other nuclei plots you asked for:

http://webusers.physics.umn.edu/~childe ... eiPlot.jpg

proton is red
helium is green
lithium is blue
beryllium is yellow
boron is cyan
carbon is light blue

Sorry this took so long, I've been getting ready for the ICRC, too. Perhaps I'll see you there.

Thanks,
Taylor


Top
 Profile  
 
 Post subject: B/C diffusion constant
PostPosted: Thu Jun 28, 2007 4:39 pm 
Offline

Joined: Sun Feb 11, 2007 4:13 pm
Posts: 9
Location: Minneapolis, MN
I've also been trying to reproduce your diffusion coefficient work from your 1998 ApJ paper.

If I take 599278 and change its parameters to match those of your run with the index 01000.I don't get the same results

Here is my galdef file

http://webusers.physics.umn.edu/~childe ... 50p_600025

here is a plot of its output:

http://webusers.physics.umn.edu/~childers/BCRatioA.jpg

As you can see, it doesn't reproduce your fitting results of the HEAO-3 data.

Thanks,
Taylor


Top
 Profile  
 
 Post subject:
PostPosted: Thu Jun 28, 2007 6:51 pm 
Offline

Joined: Tue Nov 29, 2005 6:07 pm
Posts: 34
Location: Stanford
The plots look terrible!
Something is definitively wrong.

_________________
Igor Moskalenko
Stanford University


Top
 Profile  
 
 Post subject:
PostPosted: Thu Jun 28, 2007 6:54 pm 
Offline

Joined: Tue Nov 29, 2005 6:07 pm
Posts: 34
Location: Stanford
The paper of 1998 had used old cross sections.
I suggest you to try to use the parameters of Ptuskin et al. (2006) or
from Strong et al. (2004).

I am currently thinking of cancelling my travel to Merida, but if I go
there we can certainly solve the problem.

_________________
Igor Moskalenko
Stanford University


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 25 posts ]  Go to page 1, 2  Next

All times are UTC - 8 hours [ DST ]


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group