[G12] GFlux
Michael C. Kunkel
mkunkel at jlab.org
Fri Mar 28 15:55:45 EDT 2014
Greetings,
I am wondering where we are at with the flux issue?
Meaning, what is going on with the flux in the 3.5 GeV region?
I feel that the flux is in error for a few reasons.
First when I bin flux in 15 MeV, I see 0 flux at 3.515 - 3.53.
Second, when I bin flux in 50 MeV binning I see a drop in flux at 3.5 MeV
Here is the output
//
3.45 3.5 1.06825209003e+12
3.5 3.55 470202243457.0
3.55 3.6 1.69234003548e+12
//
Third reason is that if you look at my attached pdf, I show the pi0
cross section for various energies, if you were to go to compare the
cross sections around 3.5 GeV to the neighboring Egamma bins, you will
notice a sizable discrepancy. See G12_XSection_ALL.pdf
Fourth reason, if I were to compare the pi0 cross section in unit of
dsigma/dt and compare it other experiments at 4 and 5 GeV, our cross
section is off by a factor of 2, see t_vs_dsigdt_Compare.pdf
In order for any g12 papers to be published, we need verification of
cross sections, which as of this moment I do not believe we have.
BR
MK
On 3/25/14 10:16 PM, Eugene Pasyuk wrote:
> This is not quite correct, especially when bins are small and
> comparable with the width of E-bin.
> Just imagine, You energy bin is 3.000-3.015 GeV. E-bin let's say
> covers 3.0121-3.0181 GeV. The center of E-bin will be 3.01505 GeV.
> When you enter this into the histogram it will fall to 3.015-3.030 bin
> then. However, this E-bin should be split between two adjacent bins
> in new binning scheme. There are two possible ways of doing this.
> Either randomize energy within the width of the E-bin before you enter
> it to th ehistogram, or split it between adjacent bins using weight
> which depends on how the bin boundary splits E-bin and have two
> entries to the histogram with sum of the weights equal one. You can
> imagine, if the bin width is chosen to be less than the widths of
> E-bin, then you will get completely empty bins in resulting histogram.
>
> -Eugene
>
> ------------------------------------------------------------------------
>
> *From: *"Johann Goetz" <theodore.goetz at gmail.com>
> *To: *"Eugene Pasyuk" <pasyuk at jlab.org>
> *Sent: *Tuesday, March 25, 2014 9:12:49 PM
> *Subject: *Re: [G12] GFlux
>
> I consider the beam energy to be the center of the energy paddle
> and fill the histogram at that energy with the number if counts in
> the gflux file. Is that wrong?
>
> On Mar 25, 2014 5:59 PM, "Eugene Pasyuk" <pasyuk at jlab.org
> <mailto:pasyuk at jlab.org>> wrote:
>
> I see four numbers in each row in this file. I figured the
> first one is E-bin number. What are the other three? I guess
> the second and the fourth are some energies. What is the third
> number which is integer and varies a lot.
> When you say you "histogram" does it mean you take one fixed
> number for each E-bin (beam centroid?) and put it in the
> histogram, or something else?
>
> -Eugene
>
> ------------------------------------------------------------------------
>
> *From: *"Johann Goetz" <theodore.goetz at gmail.com
> <mailto:theodore.goetz at gmail.com>>
> *To: *"g12" <g12 at jlab.org <mailto:g12 at jlab.org>>
> *Sent: *Tuesday, March 25, 2014 5:40:42 PM
> *Subject: *Re: [G12] GFlux
>
> The output of gflux includes a file that has the energies
> of the individual energy paddles and the counts. These
> files are then used to fill a histogram in the g12-gflux
> script where the range and bins are specified on the
> command line. The energies are slightly "corrected" to the
> ones found in the following file:
> /group/clas/parms/pcor/g12/tagger_energies.dat
> But, the correction is very small and averaged over the
> whole run period (i.e. not one a run-by-run basis).
>
> And before you ask, the the summing of the flux is valid
> to at least 4 figures - probably 5.
>
> -- Johann <https://sites.google.com/site/theodoregoetz/>
>
>
> On Tue, Mar 25, 2014 at 4:28 PM, Michael C. Kunkel
> <mkunkel at jlab.org <mailto:mkunkel at jlab.org>> wrote:
>
> Greetings,
>
> If the data shows a hole at beam energies, how does
> the gflux for this
> energy have events?
> Below is a table of Gflux for 15MeV binning, which is
> to fine for
> analysis but great for illustrating this question.
>
> Also, notice the 0 entries at 3.5 GeV, and low events
> in neighbors.
> FORMAT: Emin:Emax:Flux
>
>
> 1.115 1.13 0.0
> 1.13 1.145 25237172620.6
> 1.145 1.16 37307196694.5
> 1.16 1.175 596831493113.0
> 1.175 1.19 673624146169.0
> 1.19 1.205 1.04119981246e+12
> 1.205 1.22 1.13209512976e+12
> 1.22 1.235 1.31100861162e+12
> 1.235 1.25 1.41003737583e+12
> 1.25 1.265 708211172683.0
> 1.265 1.28 1.00917704833e+12
> 1.28 1.295 1.73677052381e+12
> 1.295 1.31 1.03214604848e+12
> 1.31 1.325 885037422087.0
> 1.325 1.34 1.24724378455e+12
> 1.34 1.355 857127526644.0
> 1.355 1.37 278577324131.0
> 1.37 1.385 881394713415.0
> 1.385 1.4 1.00989383319e+12
> 1.4 1.415 1.17772439858e+12
> 1.415 1.43 1.49431916406e+12
> 1.43 1.445 676795968575.0
> 1.445 1.46 1.39265456336e+12
> 1.46 1.475 1.37504480556e+12
> 1.475 1.49 737560406053.0
> 1.49 1.505 1.30281481429e+12
> 1.505 1.52 1.31266851454e+12
> 1.52 1.535 1.14281511307e+12
> 1.535 1.55 610057791731.0
> 1.55 1.565 1.19234909002e+12
> 1.565 1.58 498029161281.0
> 1.58 1.595 774221479808.0
> 1.595 1.61 1.19361035315e+12
> 1.61 1.625 1.06766392605e+12
> 1.625 1.64 933519396135.0
> 1.64 1.655 585102241894.0
> 1.655 1.67 800861230220.0
> 1.67 1.685 968140194800.0
> 1.685 1.7 616308054000.0
> 1.7 1.715 912309428155.0
> 1.715 1.73 853024554819.0
> 1.73 1.745 710609203700.0
> 1.745 1.76 1.08815959342e+12
> 1.76 1.775 1.02918324853e+12
> 1.775 1.79 752183611352.0
> 1.79 1.805 1.03490939367e+12
> 1.805 1.82 730065288937.0
> 1.82 1.835 1.01433906239e+12
> 1.835 1.85 1.02979831726e+12
> 1.85 1.865 615893688685.0
> 1.865 1.88 874150396134.0
> 1.88 1.895 896177857404.0
> 1.895 1.91 571292987771.0
> 1.91 1.925 871214929019.0
> 1.925 1.94 822118426450.0
> 1.94 1.955 611615944721.0
> 1.955 1.97 965927104595.0
> 1.97 1.985 822267766243.0
> 1.985 2.0 492626746790.0
> 2.0 2.015 862753090828.0
> 2.015 2.03 941649256469.0
> 2.03 2.045 622082896902.0
> 2.045 2.06 820144299324.0
> 2.06 2.075 884652021494.0
> 2.075 2.09 661345117113.0
> 2.09 2.105 957598861367.0
> 2.105 2.12 708490042103.0
> 2.12 2.135 566369924275.0
> 2.135 2.15 474546075393.0
> 2.15 2.165 692894638389.0
> 2.165 2.18 633696565449.0
> 2.18 2.195 1.25916127707e+12
> 2.195 2.21 452535985730.0
> 2.21 2.225 888034787498.0
> 2.225 2.24 927439628701.0
> 2.24 2.255 765262297794.0
> 2.255 2.27 604926963834.0
> 2.27 2.285 629335652375.0
> 2.285 2.3 514910298238.0
> 2.3 2.315 812968923663.0
> 2.315 2.33 540255690719.0
> 2.33 2.345 710261806621.0
> 2.345 2.36 768397230819.0
> 2.36 2.375 582062396011.0
> 2.375 2.39 630704170058.0
> 2.39 2.405 705645287962.0
> 2.405 2.42 583508371550.0
> 2.42 2.435 719849272183.0
> 2.435 2.45 596094143707.0
> 2.45 2.465 564090037414.0
> 2.465 2.48 674316483127.0
> 2.48 2.495 495228291920.0
> 2.495 2.51 796123221870.0
> 2.51 2.525 771594160133.0
> 2.525 2.54 492411586554.0
> 2.54 2.555 587323426393.0
> 2.555 2.57 395800237494.0
> 2.57 2.585 497141999612.0
> 2.585 2.6 569909994439.0
> 2.6 2.615 471441993088.0
> 2.615 2.63 471790834414.0
> 2.63 2.645 577739104258.0
> 2.645 2.66 446897506369.0
> 2.66 2.675 425924604219.0
> 2.675 2.69 657284720888.0
> 2.69 2.705 509839235972.0
> 2.705 2.72 475203266123.0
> 2.72 2.735 412635840375.0
> 2.735 2.75 674519870397.0
> 2.75 2.765 389624544425.0
> 2.765 2.78 609893184411.0
> 2.78 2.795 598744224578.0
> 2.795 2.81 422552580860.0
> 2.81 2.825 678610683392.0
> 2.825 2.84 409695314563.0
> 2.84 2.855 627296617807.0
> 2.855 2.87 573039567654.0
> 2.87 2.885 435133697461.0
> 2.885 2.9 649347542502.0
> 2.9 2.915 438132822181.0
> 2.915 2.93 596045706663.0
> 2.93 2.945 612238952007.0
> 2.945 2.96 422719982862.0
> 2.96 2.975 588371819982.0
> 2.975 2.99 616191013826.0
> 2.99 3.005 540926171982.0
> 3.005 3.02 318397775104.0
> 3.02 3.035 364624375976.0
> 3.035 3.05 369843650705.0
> 3.05 3.065 275693588687.0
> 3.065 3.08 365604500894.0
> 3.08 3.095 250216265266.0
> 3.095 3.11 375675031024.0
> 3.11 3.125 476810476117.0
> 3.125 3.14 251390414115.0
> 3.14 3.155 419690996436.0
> 3.155 3.17 286013104573.0
> 3.17 3.185 380995241526.0
> 3.185 3.2 327030646036.0
> 3.2 3.215 274517324906.0
> 3.215 3.23 416726244951.0
> 3.23 3.245 559740826159.0
> 3.245 3.26 387075736248.0
> 3.26 3.275 343494698289.0
> 3.275 3.29 467326628572.0
> 3.29 3.305 494266990137.0
> 3.305 3.32 472838022140.0
> 3.32 3.335 249310915013.0
> 3.335 3.35 283933122192.0
> 3.35 3.365 597494155644.0
> 3.365 3.38 347709473851.0
> 3.38 3.395 469145628168.0
> 3.395 3.41 235898812493.0
> 3.41 3.425 362199412402.0
> 3.425 3.44 532701421569.0
> 3.44 3.455 467292197773.0
> 3.455 3.47 453447282168.0
> 3.47 3.485 361436533945.0
> 3.485 3.5 449182174859.0
> 3.5 3.515 179551703309.0
> 3.515 3.53 0.0
> 3.53 3.545 350280790151.0
> 3.545 3.56 637140555692.0
> 3.56 3.575 588057645110.0
> 3.575 3.59 541814354100.0
> 3.59 3.605 767261690036.0
> 3.605 3.62 332154566963.0
> 3.62 3.635 366426679394.0
> 3.635 3.65 541240533201.0
> 3.65 3.665 361641638898.0
> 3.665 3.68 223165320203.0
> 3.68 3.695 541460269417.0
> 3.695 3.71 282057181004.0
> 3.71 3.725 318811531530.0
> 3.725 3.74 452785379096.0
> 3.74 3.755 274496168574.0
> 3.755 3.77 422404313178.0
> 3.77 3.785 434777507704.0
> 3.785 3.8 282152873690.0
> 3.8 3.815 388083191100.0
> 3.815 3.83 399265887124.0
> 3.83 3.845 322338229290.0
> 3.845 3.86 338018634439.0
> 3.86 3.875 219276967417.0
> 3.875 3.89 465804890065.0
> 3.89 3.905 377518940783.0
> 3.905 3.92 494456690177.0
> 3.92 3.935 433453516978.0
> 3.935 3.95 252274468259.0
> 3.95 3.965 475373426165.0
> 3.965 3.98 282811468571.0
> 3.98 3.995 406045132039.0
> 3.995 4.01 312730315238.0
> 4.01 4.025 470594492917.0
> 4.025 4.04 349257729547.0
> 4.04 4.055 279260364376.0
> 4.055 4.07 434362699962.0
> 4.07 4.085 238282355243.0
> 4.085 4.1 424665409149.0
> 4.1 4.115 280380792467.0
> 4.115 4.13 447525411083.0
> 4.13 4.145 389218745036.0
> 4.145 4.16 272773744561.0
> 4.16 4.175 416668774571.0
> 4.175 4.19 457741149259.0
> 4.19 4.205 270474957581.0
> 4.205 4.22 462790713555.0
> 4.22 4.235 332735623646.0
> 4.235 4.25 275662853810.0
> 4.25 4.265 305513621376.0
> 4.265 4.28 410076213126.0
> 4.28 4.295 387557655839.0
> 4.295 4.31 233523813972.0
> 4.31 4.325 400879966345.0
> 4.325 4.34 396154678062.0
> 4.34 4.355 214336871573.0
> 4.355 4.37 207054972335.0
> 4.37 4.385 354432043454.0
> 4.385 4.4 268714617035.0
> 4.4 4.415 174979042786.0
> 4.415 4.43 173375313599.0
> 4.43 4.445 179570800942.0
> 4.445 4.46 361918729954.0
> 4.46 4.475 429970441172.0
> 4.475 4.49 297198972337.0
> 4.49 4.505 459418373115.0
> 4.505 4.52 345874555789.0
> 4.52 4.535 426103517141.0
> 4.535 4.55 309209988159.0
> 4.55 4.565 304951470359.0
> 4.565 4.58 450559434174.0
> 4.58 4.595 277916355849.0
> 4.595 4.61 200280173919.0
> 4.61 4.625 420977338325.0
> 4.625 4.64 346795076479.0
> 4.64 4.655 303571900157.0
> 4.655 4.67 349237466285.0
> 4.67 4.685 273303785583.0
> 4.685 4.7 345523192246.0
> 4.7 4.715 362458201636.0
> 4.715 4.73 252864810085.0
> 4.73 4.745 287431010095.0
> 4.745 4.76 407384848522.0
> 4.76 4.775 223841624305.0
> 4.775 4.79 203527451994.0
> 4.79 4.805 186064192378.0
> 4.805 4.82 414281456782.0
> 4.82 4.835 376032383951.0
> 4.835 4.85 200364115883.0
> 4.85 4.865 401230883097.0
> 4.865 4.88 391468496185.0
> 4.88 4.895 354667712464.0
> 4.895 4.91 279441107884.0
> 4.91 4.925 290521222455.0
> 4.925 4.94 427805095943.0
> 4.94 4.955 238396835185.0
> 4.955 4.97 237087398075.0
> 4.97 4.985 238677728278.0
> 4.985 5.0 466393033680.0
> 5.0 5.015 201556281112.0
> 5.015 5.03 299745410214.0
> 5.03 5.045 341961220657.0
> 5.045 5.06 305926541114.0
> 5.06 5.075 385083874020.0
> 5.075 5.09 238658682950.0
> 5.09 5.105 338588394694.0
> 5.105 5.12 422223883425.0
> 5.12 5.135 298956563429.0
> 5.135 5.15 387297275757.0
> 5.15 5.165 280783183169.0
> 5.165 5.18 257466628751.0
> 5.18 5.195 393576308924.0
> 5.195 5.21 314127373912.0
> 5.21 5.225 357980308412.0
> 5.225 5.24 245073968465.0
> 5.24 5.255 196708502651.0
> 5.255 5.27 281339784057.0
> 5.27 5.285 181803381088.0
> 5.285 5.3 325929092225.0
> 5.3 5.315 219474051464.0
> 5.315 5.33 449669312291.0
> 5.33 5.345 187979242858.0
> 5.345 5.36 424496739655.0
> 5.36 5.375 195913478448.0
> 5.375 5.39 407973783769.0
> 5.39 5.405 390351891472.0
> 5.405 5.42 191462081973.0
> 5.42 5.435 425907962446.0
> 5.435 5.45 242697081405.0
> 5.45 5.465 0.0
> 5.465 5.48 0.0
> 5.48 5.495 0.0
>
> _______________________________________________
> G12 mailing list
> G12 at jlab.org <mailto:G12 at jlab.org>
> https://mailman.jlab.org/mailman/listinfo/g12
>
>
>
> _______________________________________________
> G12 mailing list
> G12 at jlab.org <mailto:G12 at jlab.org>
> https://mailman.jlab.org/mailman/listinfo/g12
>
>
>
>
>
> _______________________________________________
> G12 mailing list
> G12 at jlab.org
> https://mailman.jlab.org/mailman/listinfo/g12
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://mailman.jlab.org/pipermail/g12/attachments/20140328/96e70376/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: G12_XSection_ALL.pdf
Type: application/pdf
Size: 287344 bytes
Desc: not available
Url : https://mailman.jlab.org/pipermail/g12/attachments/20140328/96e70376/attachment-0002.pdf
-------------- next part --------------
A non-text attachment was scrubbed...
Name: t_vs_dsigdt_Comare.pdf
Type: application/pdf
Size: 20872 bytes
Desc: not available
Url : https://mailman.jlab.org/pipermail/g12/attachments/20140328/96e70376/attachment-0003.pdf
More information about the G12
mailing list