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.

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...