LDMX-Software / G4DarkBreM

Dark Bremmstrahlung in Geant4 with MadGraph event libraries.
https://ldmx-software.github.io/G4DarkBreM/
Apache License 2.0
0 stars 1 forks source link

kinematic limit within /differential/ xsec is MA not 2MA #18

Closed tomeichlersmith closed 1 year ago

tomeichlersmith commented 1 year ago

this was the issue causing the higher effective minimum. x is already limited from above by xmax which includes one of the values of MA, the minimum value of x within the differential xsec is MA because it must be enough energy for the A' to exist

To Do

tomeichlersmith commented 1 year ago

Alright, I'm not going to fully remake all of the plots, but I can quickly do a simpler comparison that shows this update is effective.

I ran the same command on both main and this branch and then compared the resulting CSV files. It shows we recover a small slice of previously excluded phase space, bringing up all values of the cross section and lifting the range between $2m_A$ and $3m_A$ away from zero.

Command: g4db-xsec-calc --ap-mass 1 --method fullww

Difference:

2,22c2,22
< 183.84,74,4000,1820.77
< 183.84,74,3900,1567.81
< 183.84,74,3800,1332.14
< 183.84,74,3700,1114.76
< 183.84,74,3600,912.299
< 183.84,74,3500,726.956
< 183.84,74,3400,555.584
< 183.84,74,3300,398.905
< 183.84,74,3200,255.426
< 183.84,74,3100,123.297
< 183.84,74,3000,0
< 183.84,74,2900,0
< 183.84,74,2800,0
< 183.84,74,2700,0
< 183.84,74,2600,0
< 183.84,74,2500,0
< 183.84,74,2400,0
< 183.84,74,2300,0
< 183.84,74,2200,0
< 183.84,74,2100,0
< 183.84,74,2000,0
---
> 183.84,74,4000,1964.53
> 183.84,74,3900,1720.96
> 183.84,74,3800,1496.3
> 183.84,74,3700,1290.09
> 183.84,74,3600,1101.88
> 183.84,74,3500,931.169
> 183.84,74,3400,777.474
> 183.84,74,3300,640.27
> 183.84,74,3200,518.995
> 183.84,74,3100,413.04
> 183.84,74,3000,321.736
> 183.84,74,2900,244.335
> 183.84,74,2800,179.999
> 183.84,74,2700,127.783
> 183.84,74,2600,86.6259
> 183.84,74,2500,55.343
> 183.84,74,2400,32.637
> 183.84,74,2300,17.1221
> 183.84,74,2200,7.37628
> 183.84,74,2100,2.02896
> 183.84,74,2000,0.000575518