[G12] GFlux
Eugene Pasyuk
pasyuk at jlab.org
Tue Mar 25 22:16:46 EDT 2014
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
----- Original Message -----
> 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 > 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 >
> >
>
> > > To: "g12" < 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
> >
>
> > > On Tue, Mar 25, 2014 at 4:28 PM, Michael C. Kunkel <
> > > 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
> > >
> >
>
> > > > 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/20140325/877ba040/attachment-0001.html
More information about the G12
mailing list