Context Navigation

← Previous Ticket
Next Ticket →

#565 closed defect (fixed)

Quadruple tests broken in decays

Reported by:	Juergen Reuter	Owned by:	kilian
Priority:	P0	Milestone:	v2.2.0
Component:	core	Version:	2.1.1
Severity:	normal	Keywords:
Cc:

Description

I did not just add the reference output for decays_5 as quadruple test output, as the diffs do not look compatible with the double precision output. This has to be resolved:

559,562c559,562
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
604c604
<  L0j =  -7.180906037455E-01  5.910932024561E-01  3.759648216366E+00
---
>  L0j =  -3.098386676966E+00  2.210055629348E+00 -7.180906037455E-01
606c606
<  L1j =   1.103130823038E+00 -8.489169492432E-02 -5.399536115122E-01
---
>  L1j =   4.449844718999E-01 -3.174040362379E-01  1.103130823038E+00
608c608
<  L2j =  -3.174040362379E-01  1.155696962911E+00  2.109019936489E+00
---
>  L2j =  -1.369521383465E+00  1.976869176962E+00 -3.174040362379E-01
610c610
<  L3j =  -4.449844718999E-01 -8.092653483294E-02  3.224195979733E+00
---
>  L3j =  -2.920000000000E+00  1.369521383465E+00 -4.449844718999E-01
633,634c633,634
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
638,639c638,639
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
663,664c663,664
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
669,670c669,670
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
697,698c697,698
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
703,704c703,704
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
731,732c731,732
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
737,738c737,738
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
784,785c784,785
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
790,791c790,791
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
817,818c817,818
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
823,824c823,824
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
850,851c850,851
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
856,857c856,857
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
882,883c882,883
<  3 prt(o:6| 3.8043091E+02;-7.0565788E+01, 2.2833644E+02, 2.9173724E+02| 2.500000000000E+03| 3)
<  4 prt(o:-6| 1.1956909E+02;-1.9195538E+01, 4.7920510E+01, 9.5561092E+01| 2.500000000000E+03| 4)
---
>  3 prt(o:6| 2.9729248E+02;-2.1672719E+01, 1.8043105E+02, 2.2990812E+02| 2.500000000000E+03| 3)
>  4 prt(o:-6| 2.0270752E+02;-6.8088606E+01, 9.5825905E+01, 1.5739021E+02| 2.500000000000E+03| 4)
927,930c927,930
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
972c972
<  L0j =   7.180906037455E-01  5.910932024561E-01  3.759648216366E+00
---
>  L0j =  -3.098386676966E+00 -2.210055629348E+00 -7.180906037455E-01
974c974
<  L1j =   1.103130823038E+00  8.489169492432E-02  5.399536115122E-01
---
>  L1j =  -4.449844718999E-01 -3.174040362379E-01 -1.103130823038E+00
976c976
<  L2j =  -3.174040362379E-01 -1.155696962911E+00 -2.109019936489E+00
---
>  L2j =   1.369521383465E+00  1.976869176962E+00  3.174040362379E-01
978c978
<  L3j =  -4.449844718999E-01  8.092653483294E-02 -3.224195979733E+00
---
>  L3j =   2.920000000000E+00  1.369521383465E+00  4.449844718999E-01
1001,1002c1001,1002
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1006,1007c1006,1007
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1031,1032c1031,1032
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1037,1038c1037,1038
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1065,1066c1065,1066
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1071,1072c1071,1072
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1099,1100c1099,1100
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1105,1106c1105,1106
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1152,1153c1152,1153
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1158,1159c1158,1159
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1185,1186c1185,1186
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1191,1192c1191,1192
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1218,1219c1218,1219
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1224,1225c1224,1225
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1250,1251c1250,1251
<  3 prt(o:6| 2.0163091E+02; 6.5739347E+01,-9.1971578E+01,-1.5929404E+02| 2.500000000000E+03| 3)
<  4 prt(o:-6| 2.9836909E+02; 2.4021978E+01,-1.8428538E+02,-2.2800429E+02| 2.500000000000E+03| 4)
---
>  3 prt(o:6| 2.1926665E+02; 5.5466797E+01,-1.4474577E+02,-1.4679904E+02| 2.500000000000E+03| 3)
>  4 prt(o:-6| 2.8073335E+02; 3.4294529E+01,-1.3151118E+02,-2.4049929E+02| 2.500000000000E+03| 4)
1297,1298c1297,1298
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
1303,1304c1303,1304
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
1351,1352c1351,1352
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
1357,1358c1357,1358
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
1405,1406c1405,1406
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
1411,1412c1411,1412
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
1467,1468c1467,1468
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
1473,1474c1473,1474
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
1479,1480c1479,1480
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1485,1486c1485,1486
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1535,1536c1535,1536
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
1541,1542c1541,1542
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
1547,1548c1547,1548
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1553,1554c1553,1554
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02
1603,1604c1603,1604
<  E =   3.804309144358E+02
<  P =  -7.056578792639E+01  2.283364432251E+02  2.917372429569E+02
---
>  E =   2.972924756725E+02
>  P =  -2.167271921575E+01  1.804310488680E+02  2.299081249927E+02
1609,1610c1609,1610
<  E =   1.195690855642E+02
<  P =  -1.919553754180E+01  4.792051044336E+01  9.556109166386E+01
---
>  E =   2.027075243275E+02
>  P =  -6.808860625244E+01  9.582590480045E+01  1.573902096281E+02
1615,1616c1615,1616
<  E =   2.016309123203E+02
<  P =   6.573934737341E+01 -9.197157837558E+01 -1.592940418884E+02
---
>  E =   2.192666489075E+02
>  P =   5.546679660587E+01 -1.447457717987E+02 -1.467990440834E+02
1621,1622c1621,1622
<  E =   2.983690876797E+02
<  P =   2.402197809478E+01 -1.842853752929E+02 -2.280042927323E+02
---
>  E =   2.807333510925E+02
>  P =   3.429452886232E+01 -1.315111818698E+02 -2.404992905373E+02

Change History (2)

comment:1 Changed 10 years ago by kilian

Status:	new → assigned

Good one. The culprit is the Lorentz transformation rest frame -> lab frame for the decaying particle. It turns out that both solutions are entirely correct, because a particle at rest has no preferred direction of flight ...

So, let's see how to resolve the ambiguity.

comment:2 Changed 10 years ago by kilian

Resolution:	→ fixed
Status:	assigned → closed

OK, an obvious solution is not to rotate the decay anymore, after boosting to the rest frame. Done in r4874.

Note: See TracTickets for help on using tickets.

Download in other formats: