Traruh Synred

Ok. I will up load...

Answered: Traruh Synred 55

October 12 2023

0 0

@dharr Thanks.

I did find that by converting Vector to Matriix I can output it.

I will up load...

>

Vector[column](%id = 4409128066)

(1)

File output has the four numbers, one per line

>

Error, (in Export) invalid input: ExportMatrix expects its 2nd argument, M, to be of type {Matrix, list(Matrix)}, but received Vector(4, {(1) = 1, (2) = 3, (3) = 5, (4) = 7})

>

(2)

>

(3)

>

(4)

>

Matrix(%id = 4409118722)

(5)

>

(6)

>

Download testfile.mw

I get an error...

Answered: Traruh Synred 55

October 12 2023

0 0

@dharr It seems like it doesn't like Vector (which I thought was a kind of Matrix). Could this be because I only have Maple 18?

op mehtod does not work for me1...

Answered: Traruh Synred 55

October 12 2023

0 0

@gulliet That method doesn't work for me (Maple 18). It gives strings describing the list not concatenated list.

Just adding them with ''+' produces a vector. It's not hard. I could write my own and it's the obvious function to override '+' with.

@Carl Love Here is another file fro...

Answered: Traruh Synred 55

October 11 2023

0 0

@Carl Love Here is another file from an attempt at 400 event fun. It only made it to 217.

DStar3DKinRecoLoopDebugCombo400EventCrashed217.mw

Memory Leak, continued...

Answered: Traruh Synred 55

October 11 2023

0 0

@Carl Love Besides the more or less small continuous leaking, I have found occasional large 'bumps'. I plot the increase using 'MemUse()' for my event loop -- from start to end for each iteration. plot the log of the number of bytes. (I haven't figured out how to make a 'dataplot' using a log scale). That big spike is 10^7..10MB. At least, mostly these are associated with a similar size increase in the proc 'bDecay'. I also count ('nTryAng') the number of trials it takes 'PickAngles' to pick the angles and this does not seem to be strongly associated with the size of the bump from 'bDecay'.

For a given random seed these bumps do occur at the same event, but they vary in size. However, running the same seed multiple times does not result in a crash at that event.

While the bumps I see are not big enough to cause crashes, the crashes are so hard that I can not examine anything afterward.

Doesn't work!...

Answered: Traruh Synred 55

October 08 2023

0 0

@Joe Riel This [ (* *_ doesn't work for me, even though it shows up in Maple Help on commenting.

Multiple runs, different crac point...

Answered: Traruh Synred 55

October 06 2023

0 0

@Carl Love I've been running the 'MemUse()' disabled version with the same 'seed' multiple times and it crashes at a different event every time. This hints to me that it's not the code as I don't think there is the undefined variable being random in old Fortran occurs n Maple.

Another crash...

Answered: Traruh Synred 55

October 06 2023

0 0

@Carl Love Here's one that crashes with 'MemUse()' disabled (replaced with a version that just returns 1); It crashes earlier.DStar3DKinRecoLoopDebug.mw

'Insert Comment' didn't seem to work so I inserted a link. I'll leave it alone for a while or rename it.

Download DStar3DKinRecoLoopDebug.mw

Another Try for one with crash...

Answered: Traruh Synred 55

October 06 2023

0 0

@Carl Love This is an

B-vector mathcing for double D*LNu decays

>

(1)

>

4.9 -- Fix boostTo to use Gram-Schmit orthogonalization

5.3 -- Fixed sign error in two body decay momentum calc in doa loop

Fixed incorrect calculation of v1 and v2 which summed D* and it's pion instead of D* and the lepton

5.9 Fixed bug in 'PickAngles'
6.0 Remove 5.9..that was a mistake...

6.1 Negative values from 'ProbDLNu'. There was a missing sThetaL from the cross terms

6.2 Type causing none zero neutrion masses found and fixe

Does not stop mememory crashs they seem less common when memory is intergated?

>

estart; with(Physics); with(LinearAlgebra); with(VectorCalculus); with(Optimization); with(Statistics); with(ArrayTools); with(plots); with(Threads); with(MmaTranslator[Mma]); Setup(mathematicalnotation = true); Setup(coordinatesystems = cartesian)

$`Default differentiation variables for d_, D_ and dAlembertian are:`*{X = (x, y, z, t)}$

$`Systems of spacetime Coordinates are:`*{X = (x, y, z, t)}; "_noterminate"$

(2)

>

Here I consider the 0C reconstructing of both B-mesons in the decay BBbar --> (D*LNu)(D*Lnu). The vectors v1 and v2 are the three momenta of the D* and the lepton. The angle bbetween v1 and I is called theta12 (T12in Maple code).The B-vectors B1 and B2 branch (unlabed ugh!) of there respective v (D*+lep) vectors at angle specified. Clearly to mathch (pB1=-pB2) the must lie in the v1-v2 plane. Their angles with w.r.t., v1 and v2 are theta1 (T1) and theta2 (T2). The kinematics of the decay should guarentee that one and only one match. The only ambiguous case would be when v1 and v2 are back-to-back and this could be excluded w/o much loss.The kinematics determins the angles cos(theta1) and cos(theta2).In most cases it will be possible to reconstruct the decay angles of both D*LNu's.

>	There are 4 possible values for the angles between the B-mesons -theta12+/-theta1+/-theta2. One must by 180 deg (pi).

>

Angles to try:

>

(3)

>

(4)

>

(5)

>

(6)

>

(7)

>

(8)

>

dot4 := proc (x, y) local result; result := x[1]*y[1]-x[2]*y[2]-x[3]*y[3]-x[4]*y[4]; return result end proc

Calculate theta from contraints

>

P3 := proc (P4) local result; `#print`("P3 on P4:", P4); result := Vector([P4[2], P4[3], P4[4]]); return result end proc

>

Mass := proc (P4) local result; result := Re(sqrt(dot4(P4, P4))); return result end proc

>

momentum := proc (P4) local result, p3; `#print`(" Momentum input P4:", P4); p3 := P3(P4); `#print`("p3:", p3); result := Norm(p3, Euclidean); return result end proc

>

nHat := proc (P4) local p3, magp3, nhat; p3 := P3(P4); magp3 := Norm(p3, Euclidean); nhat := p3/magp3; return nhat end proc

sums second particle is neutriono, taken to have zero mass

>

TwoBodyDecay := proc (mMother, mDaug) local k, z, phi, zsin, result, eDaug; global rFlat; k := (mMother^2-mDaug^2)/(2.0*mMother); eDaug := sqrt(mDaug^2+k^2); z := -1.0+2.0*rFlat(); zsin := sqrt(1-z^2); phi := 2*rFlat()*Pi; result := Vector([eDaug, k*zsin*cos(phi), k*zsin*sin(phi), k*z]); return result end proc

>

Warning, `Mp4DStar` is implicitly declared local to procedure `bDecay`

Warning, `p4LepNuRest` is implicitly declared local to procedure `bDecay`

Warning, `Mp4LepNu` is implicitly declared local to procedure `bDecay`

Warning, `Mp4Lep` is implicitly declared local to procedure `bDecay`

Warning, `Mp4Nu` is implicitly declared local to procedure `bDecay`

Warning, `p4LepNuSum` is implicitly declared local to procedure `bDecay`

>

(9)

>

(10)

>

$boostTo := proc (p4::Vector, p4To::Vector, sign::float := 1.0) local p3, eNergy, pTrans, pz, zHat, yHat, xHat, betaTo, cTheta, sTheta, rcs, xz, gammaTo, eP, pzP, p3P, mp4, mag, result; `#print`("boostTo starts"); mp4 := Mass(p4); betaTo := momentum(p4To)/p4To[1]; gammaTo := 1/sqrt(1.0-betaTo^2); `##print`("betaTo,gammaTo:", betaTo, "\n", gammaTo); zHat := nHat(p4To); `#print`("zHat:", zHat); p3 := P3(p4); `#print`("p3:", p3); eNergy := p4[1]; pz := zHat.p3; pTrans := p3-pz*zHat; `#print`("pTrans:", pTrans); pzP := gammaTo*(pz-sign*betaTo*eNergy); `#print`("eP:", eP, " pzP:", pzP); `#print`("pTrans:", pTrans); `#print`("pzP*zHat:", pzP*zHat); p3P := pzP*zHat+pTrans; `#print`("p3P:", p3P); mag := Norm(p3P, Euclidean); eP := sqrt(mp4^2+mag^2); result := Vector([eP, p3P]); `#print`("boostTo stops"); return result end proc$

>

flip := proc (p4::Vector) local p3, result; p3 := P3(p4); result := Vector([p4[1], -p3]); return result end proc

defHist := proc (n, low, high, name::string) local h, info, result; h := Sample(0, n); info := Sample(0, 5); info[1] := low; info[2] := high; result := [h, info, name]; return result end proc

>

fillHist := proc (aHist, x::float) local nhist, dx, bin, h, info, low, high; h := aHist[1]; info := aHist[2]; nhist := Size(h, 2); low := info[1]; high := info[2]; dx := (high-low)/nhist; bin := trunc((x-low)/dx)+1; if bin < nhist+1 and 0 < nhist then h[bin] := h[bin]+1; info[3] := info[3]+x; info[4] := info[4]+x^2; info[5] := info[5]+1 end if; return nhist end proc

>

$plotHist := proc (aHist) local h, info, low, high, mean, rms, std, nentry, name, xcenter, xBinCenters, nhist, dx, ibin, thePlot; h := aHist[1]; info := aHist[2]; low := info[1]; high := info[2]; name := aHist[3]; mean := info[3]/info[5]; rms := info[4]/info[5]; nentry := info[5]; std := rms-mean^2; print(" mean:", mean, " std:", std, "name:", name); name := cat(name, "\n entries:", nentry, "\n mean:", mean, "\n std:", std); nhist := Size(h, 2); xBinCenters := Sample(0, 0); dx := (high-low)/nhist; for ibin to nhist do xcenter := low+(1/2)*dx+(ibin-1)*dx; Append(xBinCenters, xcenter); thePlot := plot(xBinCenters, h, style = point, axes = boxed, symbol = solidcircle, symbolsize = 10, title = name) end do; return thePlot end proc$

>

Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

>

(11)

>

(12)

Restart here to gen new event

>

(13)

>

(14)

>

count := 0; fa := 1.0; `#memSum` := 0.; nDo := 1

(15)

>

for i to nDo do `#memBefProb` := MemUse(); Prob := ProbDStarLNu(1.0, 1.0, 1.0, .2*fa, .4*fa, 1.01*fa); count := count+1; fa := 1.0001*fa; `#memAftProb` := MemUse(); `#memProb` := memAftProb-memBefProb; `#memSum` := memSum+memProb end do

>

(16)

>

Vector(4, {(1) = 2161.4540183799612632, (2) = 244.93657976024539873, (3) = 478.46812092888873589, (4) = 951.48808564322715396})

Vector(4, {(1) = 250.86998162003873729, (2) = 95.779097099681417500, (3) = 158.14695992125048399, (4) = 169.56306147364397496})

Vector(4, {(1) = 1891.3133351496805806, (2) = -989.98766795689809863, (3) = -1332.5398224582864055, (4) = -906.27143254721546083})

Vector[column](%id = 4697626498)

(17)

>

(18)

>

(19)

>

(20)

>

(21)

>

(22)

>

(23)

>

(24)

>

(25)

>

recoAngles[3]/Pi

(26)

>

(27)

>

(28)

>

(29)

>

(30)

>

Vector[column](%id = 4697181506)

(31)

>

Vector[column](%id = 4697176322)

(32)

>

Vector[column](%id = 4697182786)

(33)

>

Vector[column](%id = 4697185794)

(34)

>

(35)

>

(36)

>

(37)

>

(38)

>

(39)

>

(40)

>

(41)

>

(42)

>

(43)

>

(44)

>

(45)

>

(46)

>

(47)

>

(48)

>

(49)

>

(50)

>

(51)

>

(52)

>

Vector[row](%id = 4697612866)

(53)

>

Vector[row](%id = 4697613570)

(54)

>

(55)

>

(56)

>

(57)

>

(58)

>

(59)

>

(60)

>

(61)

>

(62)

>

(63)

>

(64)

>

(65)

>

Error, (in Vector) Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

HandCalcZ

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

Error, Maple was unable to allocate enough memory to complete this computation. Please see ?alloc

>

contraint works!

>

transformation!

>

Here need to compute theta1U and theta2U from the kinematics...

>

So solution 1 is far from the orignal input, but solution 2 virtually perfect

>

here I have B-reco BTrue comparison

>

for i from 0 to nPhi-1 do phiU := 2.0*i*Pi/nPhi; VList := [VList[], i+3]; nb1pphi := subs(phi1 = phiU, nb1p); VList[i+3] := nb1pphi end do

>

for i from 0 to nPhi-1 do phiU := 2.0*i*Pi/nPhi; VList := [VList[], i+1+nVlist]; nb2pphi := evalf(subs(phi2 = phiU, nb2p)); VList[i+1+nVList] := nb2pphi end do

>

>	Hp :=1.0; Hm:=1.0; H0:=1.0;

>

"do n:=n+1; #`#doa print`(n); angles:=PickAngles(Hp,Hm,H0); thetal:=angles[1]; thetav:=angles[2]; chi:=angles[3]; cL:=cos(thetal); Al[n]:=cL; AV[n]:=cos(thetav); AChi[n]:=chi; #hist(Al,-1.0,1.0,cL); #print("n:",n,"cL:",cL); until n > nMax-1: "

>

Download DStar3DKinRecoLoopDebugMemCrashA.mw

earlier before I turned on more detailed printing in the event loop. That one did save intact...

That's not one with crash......

Answered: Traruh Synred 55

October 06 2023

0 0

@Traruh Synred @Carl Love

That's not one with a crash. I'll try again. in a bit.

It's sort of big......

Answered: Traruh Synred 55

October 06 2023

0 0

@Carl Love It's sort of big with lots of extraneous stuff from when I was trying out things (like timing some 'proc's on there own. I just made a run in which I had disabled 'MemUse' which crashed earlier than somewhere it was being used properly. That seems to be a persistent pattern. The crash takes longer to develop of I look at memory usage.

Hmm,mm the crash is so hard it's not letting me save

B-vector mathcing for double D*LNu decays

>

(1)

>

4.9 -- Fix boostTo to use Gram-Schmit orthogonalization

5.3 -- Fixed sign error in two body decay momentum calc in doa loop

Fixed incorrect calculation of v1 and v2 which summed D* and it's pion instead of D* and the lepton

5.9 Fixed bug in 'PickAngles'
6.0 Remove 5.9..that was a mistake...

6.1 Negative values from 'ProbDLNu'. There was a missing sThetaL from the cross terms

6.2 Type causing none zero neutrion masses found and fixe

Does not stop mememory crashs they seem less common when memory is intergated?

>

estart; with(Physics); with(LinearAlgebra); with(VectorCalculus); with(Optimization); with(Statistics); with(ArrayTools); with(plots); with(Threads); with(MmaTranslator[Mma]); Setup(mathematicalnotation = true); Setup(coordinatesystems = cartesian)

$`Default differentiation variables for d_, D_ and dAlembertian are:`*{X = (x, y, z, t)}$

(2)

>

Here I consider the 0C reconstructing of both B-mesons in the decay BBbar --> (D*LNu)(D*Lnu). The vectors v1 and v2 are the three momenta of the D* and the lepton. The angle bbetween v1 and I is called theta12 (T12in Maple code).The B-vectors B1 and B2 branch (unlabed ugh!) of there respective v (D*+lep) vectors at angle specified. Clearly to mathch (pB1=-pB2) the must lie in the v1-v2 plane. Their angles with w.r.t., v1 and v2 are theta1 (T1) and theta2 (T2). The kinematics of the decay should guarentee that one and only one match. The only ambiguous case would be when v1 and v2 are back-to-back and this could be excluded w/o much loss.The kinematics determins the angles cos(theta1) and cos(theta2).In most cases it will be possible to reconstruct the decay angles of both D*LNu's.

>	There are 4 possible values for the angles between the B-mesons -theta12+/-theta1+/-theta2. One must by 180 deg (pi).

>

Angles to try:

>

(3)

>

(4)

>

(5)

>

(6)

>

(7)

>

(8)

>

dot4 := proc (x, y) local result; result := x[1]*y[1]-x[2]*y[2]-x[3]*y[3]-x[4]*y[4]; return result end proc

Calculate theta from contraints

>

P3 := proc (P4) local result; `#print`("P3 on P4:", P4); result := Vector([P4[2], P4[3], P4[4]]); return result end proc

>

Mass := proc (P4) local result; result := Re(sqrt(dot4(P4, P4))); return result end proc

>

momentum := proc (P4) local result, p3; `#print`(" Momentum input P4:", P4); p3 := P3(P4); `#print`("p3:", p3); result := Norm(p3, Euclidean); return result end proc

>

nHat := proc (P4) local p3, magp3, nhat; p3 := P3(P4); magp3 := Norm(p3, Euclidean); nhat := p3/magp3; return nhat end proc

sums second particle is neutriono, taken to have zero mass

>

TwoBodyDecay := proc (mMother, mDaug) local k, z, phi, zsin, result, eDaug; global rFlat; k := (mMother^2-mDaug^2)/(2.0*mMother); eDaug := sqrt(mDaug^2+k^2); z := -1.0+2.0*rFlat(); zsin := sqrt(1-z^2); phi := 2*rFlat()*Pi; result := Vector([eDaug, k*zsin*cos(phi), k*zsin*sin(phi), k*z]); return result end proc

>

Warning, `Mp4DStar` is implicitly declared local to procedure `bDecay`

Warning, `p4LepNuRest` is implicitly declared local to procedure `bDecay`

Warning, `Mp4LepNu` is implicitly declared local to procedure `bDecay`

Warning, `Mp4Lep` is implicitly declared local to procedure `bDecay`

Warning, `Mp4Nu` is implicitly declared local to procedure `bDecay`

Warning, `p4LepNuSum` is implicitly declared local to procedure `bDecay`

>

(9)

>

(10)

>

$boostTo := proc (p4::Vector, p4To::Vector, sign::float := 1.0) local p3, eNergy, pTrans, pz, zHat, yHat, xHat, betaTo, cTheta, sTheta, rcs, xz, gammaTo, eP, pzP, p3P, mp4, mag, result; `#print`("boostTo starts"); mp4 := Mass(p4); betaTo := momentum(p4To)/p4To[1]; gammaTo := 1/sqrt(1.0-betaTo^2); `##print`("betaTo,gammaTo:", betaTo, "\n", gammaTo); zHat := nHat(p4To); `#print`("zHat:", zHat); p3 := P3(p4); `#print`("p3:", p3); eNergy := p4[1]; pz := zHat.p3; pTrans := p3-pz*zHat; `#print`("pTrans:", pTrans); pzP := gammaTo*(pz-sign*betaTo*eNergy); `#print`("eP:", eP, " pzP:", pzP); `#print`("pTrans:", pTrans); `#print`("pzP*zHat:", pzP*zHat); p3P := pzP*zHat+pTrans; `#print`("p3P:", p3P); mag := Norm(p3P, Euclidean); eP := sqrt(mp4^2+mag^2); result := Vector([eP, p3P]); `#print`("boostTo stops"); return result end proc$

>

flip := proc (p4::Vector) local p3, result; p3 := P3(p4); result := Vector([p4[1], -p3]); return result end proc

defHist := proc (n, low, high, name::string) local h, info, result; h := Sample(0, n); info := Sample(0, 5); info[1] := low; info[2] := high; result := [h, info, name]; return result end proc

>

fillHist := proc (aHist, x::float) local nhist, dx, bin, h, info, low, high; h := aHist[1]; info := aHist[2]; nhist := Size(h, 2); low := info[1]; high := info[2]; dx := (high-low)/nhist; bin := trunc((x-low)/dx)+1; if bin < nhist+1 and 0 < nhist then h[bin] := h[bin]+1; info[3] := info[3]+x; info[4] := info[4]+x^2; info[5] := info[5]+1 end if; return nhist end proc

>

$plotHist := proc (aHist) local h, info, low, high, mean, rms, std, nentry, name, xcenter, xBinCenters, nhist, dx, ibin, thePlot; h := aHist[1]; info := aHist[2]; low := info[1]; high := info[2]; name := aHist[3]; mean := info[3]/info[5]; rms := info[4]/info[5]; nentry := info[5]; std := rms-mean^2; print(" mean:", mean, " std:", std, "name:", name); name := cat(name, "\n entries:", nentry, "\n mean:", mean, "\n std:", std); nhist := Size(h, 2); xBinCenters := Sample(0, 0); dx := (high-low)/nhist; for ibin to nhist do xcenter := low+(1/2)*dx+(ibin-1)*dx; Append(xBinCenters, xcenter); thePlot := plot(xBinCenters, h, style = point, axes = boxed, symbol = solidcircle, symbolsize = 10, title = name) end do; return thePlot end proc$

>

Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

>

MemUse := proc () return 1.0 end proc

>

(11)

>

(12)

Restart here to gen new event

>

(13)

>

(14)

>

count := 0; fa := 1.0; `#memSum` := 0.; nDo := 1

(15)

>

for i to nDo do `#memBefProb` := MemUse(); Prob := ProbDStarLNu(1.0, 1.0, 1.0, .2*fa, .4*fa, 1.01*fa); count := count+1; fa := 1.0001*fa; `#memAftProb` := MemUse(); `#memProb` := memAftProb-memBefProb; `#memSum` := memSum+memProb end do

>

(16)

>

Vector(4, {(1) = 2161.4540183799612632, (2) = 244.93657976024539873, (3) = 478.46812092888873589, (4) = 951.48808564322715396})

Vector(4, {(1) = 250.86998162003873729, (2) = 95.779097099681417500, (3) = 158.14695992125048399, (4) = 169.56306147364397496})

Vector(4, {(1) = 1891.3133351496805806, (2) = -989.98766795689809863, (3) = -1332.5398224582864055, (4) = -906.27143254721546083})

Vector(4, {(1) = 975.70266485031941945, (2) = 649.27199109697128247, (3) = 695.92474160814718574, (4) = -214.77971456965566785})

(17)

>

(18)

>

(19)

>

(20)

>

(21)

>

(22)

>

(23)

>

(24)

>

(25)

>

recoAngles[3]/Pi

(26)

>

(27)

>

(28)

>

(29)

>

(30)

>

Vector(4, {(1) = 2210.2489853584658006, (2) = -582.98360907141692211, (3) = 933.83871487847686035, (4) = 442.32083132756442472})

(31)

>

Vector(4, {(1) = 202.07501464153419981, (2) = -152.21139748029832032, (3) = 113.93423781061827915, (4) = -68.446855193093973541})

(32)

>

Vector(4, {(1) = 1240.2218105923179287, (2) = 336.03729168026476119, (3) = 325.95749810552585223, (4) = -1148.4687746638735673})

(33)

>

Vector(4, {(1) = 1626.7941894076820712, (2) = 399.15771487145048112, (3) = -1373.7304507946209916, (4) = 774.59479852940311618})

(34)

>

(35)

>

(36)

>

(37)

>

(38)

>

(39)

>

(40)

>

(41)

>

(42)

>

(43)

>

(44)

>

(45)

>

(46)

>

(47)

>

(48)

>

(49)

>

(50)

>

(51)

>

(52)

>

$z12List := Vector(4, {(1) = ` 1 .. 1000 `*Vector[row], (2) = `Data Type: `*sfloat, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})$

(53)

>

$zTS1Or2List := Vector(4, {(1) = ` 1 .. 1000 `*Vector[row], (2) = `Data Type: `*sfloat, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})$

(54)

>

(55)

>

(56)

>

(57)

>

(58)

>

(59)

>

(60)

>

(61)

>

(62)

>

(63)

>

(64)

>

(65)

>

(66)

>

(67)

>

HandCalcZ

>

contraint works!

>

transformation!

>

Here need to compute theta1U and theta2U from the kinematics...

>

So solution 1 is far from the orignal input, but solution 2 virtually perfect

>

here I have B-reco BTrue comparison

>

for i from 0 to nPhi-1 do phiU := 2.0*i*Pi/nPhi; VList := [VList[], i+3]; nb1pphi := subs(phi1 = phiU, nb1p); VList[i+3] := nb1pphi end do

>

for i from 0 to nPhi-1 do phiU := 2.0*i*Pi/nPhi; VList := [VList[], i+1+nVlist]; nb2pphi := evalf(subs(phi2 = phiU, nb2p)); VList[i+1+nVList] := nb2pphi end do

>

>	Hp :=1.0; Hm:=1.0; H0:=1.0;

>

"do n:=n+1; #`#doa print`(n); angles:=PickAngles(Hp,Hm,H0); thetal:=angles[1]; thetav:=angles[2]; chi:=angles[3]; cL:=cos(thetal); Al[n]:=cL; AV[n]:=cos(thetav); AChi[n]:=chi; #hist(Al,-1.0,1.0,cL); #print("n:",n,"cL:",cL); until n > nMax-1: "

>

Download DStar3DKinRecoLoopDebug.mw

B-vector mathcing for double D*LNu decays

>

(1)

>

4.9 -- Fix boostTo to use Gram-Schmit orthogonalization

5.3 -- Fixed sign error in two body decay momentum calc in doa loop

Fixed incorrect calculation of v1 and v2 which summed D* and it's pion instead of D* and the lepton

5.9 Fixed bug in 'PickAngles'
6.0 Remove 5.9..that was a mistake...

6.1 Negative values from 'ProbDLNu'. There was a missing sThetaL from the cross terms

6.2 Type causing none zero neutrion masses found and fixe

Does not stop mememory crashs they seem less common when memory is intergated?

>

estart; with(Physics); with(LinearAlgebra); with(VectorCalculus); with(Optimization); with(Statistics); with(ArrayTools); with(plots); with(Threads); with(MmaTranslator[Mma]); Setup(mathematicalnotation = true); Setup(coordinatesystems = cartesian)

$`Default differentiation variables for d_, D_ and dAlembertian are:`*{X = (x, y, z, t)}$

(2)

>

Here I consider the 0C reconstructing of both B-mesons in the decay BBbar --> (D*LNu)(D*Lnu). The vectors v1 and v2 are the three momenta of the D* and the lepton. The angle bbetween v1 and I is called theta12 (T12in Maple code).The B-vectors B1 and B2 branch (unlabed ugh!) of there respective v (D*+lep) vectors at angle specified. Clearly to mathch (pB1=-pB2) the must lie in the v1-v2 plane. Their angles with w.r.t., v1 and v2 are theta1 (T1) and theta2 (T2). The kinematics of the decay should guarentee that one and only one match. The only ambiguous case would be when v1 and v2 are back-to-back and this could be excluded w/o much loss.The kinematics determins the angles cos(theta1) and cos(theta2).In most cases it will be possible to reconstruct the decay angles of both D*LNu's.

>	There are 4 possible values for the angles between the B-mesons -theta12+/-theta1+/-theta2. One must by 180 deg (pi).

>

Angles to try:

>

(3)

>

(4)

>

(5)

>

(6)

>

(7)

>

(8)

>

dot4 := proc (x, y) local result; result := x[1]*y[1]-x[2]*y[2]-x[3]*y[3]-x[4]*y[4]; return result end proc

Calculate theta from contraints

>

P3 := proc (P4) local result; `#print`("P3 on P4:", P4); result := Vector([P4[2], P4[3], P4[4]]); return result end proc

>

Mass := proc (P4) local result; result := Re(sqrt(dot4(P4, P4))); return result end proc

>

momentum := proc (P4) local result, p3; `#print`(" Momentum input P4:", P4); p3 := P3(P4); `#print`("p3:", p3); result := Norm(p3, Euclidean); return result end proc

>

nHat := proc (P4) local p3, magp3, nhat; p3 := P3(P4); magp3 := Norm(p3, Euclidean); nhat := p3/magp3; return nhat end proc

sums second particle is neutriono, taken to have zero mass

>

TwoBodyDecay := proc (mMother, mDaug) local k, z, phi, zsin, result, eDaug; global rFlat; k := (mMother^2-mDaug^2)/(2.0*mMother); eDaug := sqrt(mDaug^2+k^2); z := -1.0+2.0*rFlat(); zsin := sqrt(1-z^2); phi := 2*rFlat()*Pi; result := Vector([eDaug, k*zsin*cos(phi), k*zsin*sin(phi), k*z]); return result end proc

>

Warning, `Mp4DStar` is implicitly declared local to procedure `bDecay`

Warning, `p4LepNuRest` is implicitly declared local to procedure `bDecay`

Warning, `Mp4LepNu` is implicitly declared local to procedure `bDecay`

Warning, `Mp4Lep` is implicitly declared local to procedure `bDecay`

Warning, `Mp4Nu` is implicitly declared local to procedure `bDecay`

Warning, `p4LepNuSum` is implicitly declared local to procedure `bDecay`

>

(9)

>

(10)

>

$boostTo := proc (p4::Vector, p4To::Vector, sign::float := 1.0) local p3, eNergy, pTrans, pz, zHat, yHat, xHat, betaTo, cTheta, sTheta, rcs, xz, gammaTo, eP, pzP, p3P, mp4, mag, result; `#print`("boostTo starts"); mp4 := Mass(p4); betaTo := momentum(p4To)/p4To[1]; gammaTo := 1/sqrt(1.0-betaTo^2); `##print`("betaTo,gammaTo:", betaTo, "\n", gammaTo); zHat := nHat(p4To); `#print`("zHat:", zHat); p3 := P3(p4); `#print`("p3:", p3); eNergy := p4[1]; pz := zHat.p3; pTrans := p3-pz*zHat; `#print`("pTrans:", pTrans); pzP := gammaTo*(pz-sign*betaTo*eNergy); `#print`("eP:", eP, " pzP:", pzP); `#print`("pTrans:", pTrans); `#print`("pzP*zHat:", pzP*zHat); p3P := pzP*zHat+pTrans; `#print`("p3P:", p3P); mag := Norm(p3P, Euclidean); eP := sqrt(mp4^2+mag^2); result := Vector([eP, p3P]); `#print`("boostTo stops"); return result end proc$

>

flip := proc (p4::Vector) local p3, result; p3 := P3(p4); result := Vector([p4[1], -p3]); return result end proc

defHist := proc (n, low, high, name::string) local h, info, result; h := Sample(0, n); info := Sample(0, 5); info[1] := low; info[2] := high; result := [h, info, name]; return result end proc

>

fillHist := proc (aHist, x::float) local nhist, dx, bin, h, info, low, high; h := aHist[1]; info := aHist[2]; nhist := Size(h, 2); low := info[1]; high := info[2]; dx := (high-low)/nhist; bin := trunc((x-low)/dx)+1; if bin < nhist+1 and 0 < nhist then h[bin] := h[bin]+1; info[3] := info[3]+x; info[4] := info[4]+x^2; info[5] := info[5]+1 end if; return nhist end proc

>

$plotHist := proc (aHist) local h, info, low, high, mean, rms, std, nentry, name, xcenter, xBinCenters, nhist, dx, ibin, thePlot; h := aHist[1]; info := aHist[2]; low := info[1]; high := info[2]; name := aHist[3]; mean := info[3]/info[5]; rms := info[4]/info[5]; nentry := info[5]; std := rms-mean^2; print(" mean:", mean, " std:", std, "name:", name); name := cat(name, "\n entries:", nentry, "\n mean:", mean, "\n std:", std); nhist := Size(h, 2); xBinCenters := Sample(0, 0); dx := (high-low)/nhist; for ibin to nhist do xcenter := low+(1/2)*dx+(ibin-1)*dx; Append(xBinCenters, xcenter); thePlot := plot(xBinCenters, h, style = point, axes = boxed, symbol = solidcircle, symbolsize = 10, title = name) end do; return thePlot end proc$

>

Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

>

MemUse := proc () return 1.0 end proc

>

(11)

>

(12)

Restart here to gen new event

>

(13)

>

(14)

>

count := 0; fa := 1.0; `#memSum` := 0.; nDo := 1

(15)

>

for i to nDo do `#memBefProb` := MemUse(); Prob := ProbDStarLNu(1.0, 1.0, 1.0, .2*fa, .4*fa, 1.01*fa); count := count+1; fa := 1.0001*fa; `#memAftProb` := MemUse(); `#memProb` := memAftProb-memBefProb; `#memSum` := memSum+memProb end do

>

(16)

>

Vector(4, {(1) = 2161.4540183799612632, (2) = 244.93657976024539873, (3) = 478.46812092888873589, (4) = 951.48808564322715396})

Vector(4, {(1) = 250.86998162003873729, (2) = 95.779097099681417500, (3) = 158.14695992125048399, (4) = 169.56306147364397496})

Vector(4, {(1) = 1891.3133351496805806, (2) = -989.98766795689809863, (3) = -1332.5398224582864055, (4) = -906.27143254721546083})

Vector(4, {(1) = 975.70266485031941945, (2) = 649.27199109697128247, (3) = 695.92474160814718574, (4) = -214.77971456965566785})

(17)

>

(18)

>

(19)

>

(20)

>

(21)

>

(22)

>

(23)

>

(24)

>

(25)

>

recoAngles[3]/Pi

(26)

>

(27)

>

(28)

>

(29)

>

(30)

>

Vector(4, {(1) = 2210.2489853584658006, (2) = -582.98360907141692211, (3) = 933.83871487847686035, (4) = 442.32083132756442472})

(31)

>

Vector(4, {(1) = 202.07501464153419981, (2) = -152.21139748029832032, (3) = 113.93423781061827915, (4) = -68.446855193093973541})

(32)

>

Vector(4, {(1) = 1240.2218105923179287, (2) = 336.03729168026476119, (3) = 325.95749810552585223, (4) = -1148.4687746638735673})

(33)

>

Vector(4, {(1) = 1626.7941894076820712, (2) = 399.15771487145048112, (3) = -1373.7304507946209916, (4) = 774.59479852940311618})

(34)

>

(35)

>

(36)

>

(37)

>

(38)

>

(39)

>

(40)

>

(41)

>

(42)

>

(43)

>

(44)

>

(45)

>

(46)

>

(47)

>

(48)

>

(49)

>

(50)

>

(51)

>

(52)

>

$z12List := Vector(4, {(1) = ` 1 .. 1000 `*Vector[row], (2) = `Data Type: `*sfloat, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})$

(53)

>

$zTS1Or2List := Vector(4, {(1) = ` 1 .. 1000 `*Vector[row], (2) = `Data Type: `*sfloat, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})$

(54)

>

(55)

>

(56)

>

(57)

>

(58)

>

(59)

>

(60)

>

(61)

>

(62)

>

(63)

>

(64)

>

(65)

>

(66)

>

(67)

>

HandCalcZ

>

contraint works!

>

transformation!

>

Here need to compute theta1U and theta2U from the kinematics...

>

So solution 1 is far from the orignal input, but solution 2 virtually perfect

>

here I have B-reco BTrue comparison

>

for i from 0 to nPhi-1 do phiU := 2.0*i*Pi/nPhi; VList := [VList[], i+3]; nb1pphi := subs(phi1 = phiU, nb1p); VList[i+3] := nb1pphi end do

>

for i from 0 to nPhi-1 do phiU := 2.0*i*Pi/nPhi; VList := [VList[], i+1+nVlist]; nb2pphi := evalf(subs(phi2 = phiU, nb2p)); VList[i+1+nVList] := nb2pphi end do

>

>	Hp :=1.0; Hm:=1.0; H0:=1.0;

>

"do n:=n+1; #`#doa print`(n); angles:=PickAngles(Hp,Hm,H0); thetal:=angles[1]; thetav:=angles[2]; chi:=angles[3]; cL:=cos(thetal); Al[n]:=cL; AV[n]:=cos(thetav); AChi[n]:=chi; #hist(Al,-1.0,1.0,cL); #print("n:",n,"cL:",cL); until n > nMax-1: "

>

Download DStar3DKinRecoLoopDebug.mw

it. Maybe I'll get a chance to do it when I exit Maple...DStar3DKinRecoLoopDebug.mw

Back to Memory...

Answered: Traruh Synred 55

October 06 2023

0 0

@Carl Love

In my event loop occasionally only very large additions to memory. Typically, a single pass through the loop adds about 30KB according to 'MemUse()'. Event 46 the snapshot below is different.

dMemLoop is memory added in one iteration of the event loop. dMemStartLoop is memory change at the start of iteration (from the end of the previous one). It typically gives back some memory.

Another odd thing si that while 'MemUse()' indicates 20MB being used Task Manger indicates 1000s of MB. Could this be a Windows/Java problem rather than a Maple problem?

Another>MB single iteration bump occurs a few events before the crash. The crash might be caused by such a bump, perhaps a bigger one that doesn't get printed could have caused the crash.

218 events! Then memory crash!...

Answered: Traruh Synred 55

October 06 2023

0 0

@Carl Love Trying to run 1000 events with neutrino mass bug fixed (a typo of course) the progam crashes at event 218. This is going to prevent me from collecting adequate stats for my project. Memory grew to ~1500 MB during the course of the run and toal mem use on my machine was about 60%. Still it might have run into a Maple soft limit. I don't known.

Currently in crased state Maple is taking ~200MB.

I suspect there is an occasion catostrophic mem grab. I was hoping that fixing the neutrio mass bug and the resulting complex number problems of 'NPSolve' would solve the problem, but apparently not.

?

Memory grew to ~1500MB. Below error on crash followed by Task Manger out put just before crash:

I need Four Vectors...

Answered: Traruh Synred 55

October 05 2023

0 0

@ecterrab I need four vectors and a boost function. I'll look to see if tensors as such a thing.

I found part of my problems my boost. It still doing weried thank. Boost a zero mass 4-vector some times yields a four vector with mass.

-Thanks

Well, boost problem was simple and stupi...

Answered: Traruh Synred 55

October 05 2023

0 0

I was simply boosting the wrong way...nothing wrong with 'boostTo'..

There is another probllem that a cosine I calculate doesn't always come out less 1 causing NPSolve to chock on complex numbers. I thought this related to non-zero recoil neutriono masses as the caculation is based on neutrino mass being 0, but apparently that's not the problem as it still occurs with neutriono mass problem fixed. I don't think it's related to mem leak either.

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post

55 Reputation

3 Badges

MaplePrimes Activity

These are replies submitted by Traruh Synred

Ok. I will up load...

I get an error...

op mehtod does not work for me1...

@Carl Love Here is another file fro...

Memory Leak, continued...

Doesn't work!...

Multiple runs, different crac point...

Another crash...

Another Try for one with crash...

That's not one with crash......

It's sort of big......

Back to Memory...

218 events! Then memory crash!...

I need Four Vectors...

Well, boost problem was simple and stupi...

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again