1^0=\omega_2(\mathbf{k}_1) \\ k_2^0=\omega_{\lambda'}(\mathbf{k}_2)}}\,,\end{aligned}$$

where the energy of the final electron is denoted as $E_2$. Note the division by the group velocity of the first and second polarization state (see also Ref. [@Colladay:2001wk]), respectively, where in the standard theory $|\mathbf{v}_{\mathrm{gr},\,1}|=|\mathbf{v}_{\mathrm{gr},\,2}|=1$. The tensor $(|\mathcal{M}(k_1,k_2)|^2\big)^{\mu\nu\varrho\sigma}$ is given by the trace term of Eq. (5.81) of Ref. [@PeskinSchroeder1995] with some modifications due to the Lorentz-violating kinematics.

![Chosen coordinate system for the phase space integration, where the initial photon three-momentum $\mathbf{k}_1$ lies along the third axis of the coordinate system. For the outgoing photon momentum $\mathbf{k}_2$, spherical coordinates are chosen with the azimuthal angle $\vartheta$ corresponding to the angle between $\mathbf{k}_2$ and the third axis. Cases are treated with the three-vector $\boldsymbol{\zeta}$ having equal or unequal components.[]{data-label="fig:phase-space-integration"}](consistency-modmax-parity-odd_v1_fig4_coordinates.pdf)

The purely algebraic part of the calculation, that includes computation of traces, contraction of indices and inserting kinematical relations, is performed with `Form` [@Vermaseren:2000nd]. The subsequent phase space calculation is done numerically with `C++`, since the resulting matrix element squared contains hundreds of terms. The limit of zero Lorentz violation has to be taken with care and “long double” precision does not suffice here. Therefore, the GNU Multiple Precision Arithmetic Library `GMP` [@GMP:2011] is used with its `C++` interface described in Sec. 12 of the reference previously mentioned.

The first idea was to choose the coordinate system such that $\boldsymbol{\zeta}$ lies along the third axis. Then the phase space should have been integrated with cylindrical coordinates $(k_{2,\bot},\varphi,k_{2,\parallel})$. To cover a general situation, where the initial photon momentum $\mathbf{k}_1$ points in an arbitrary
1 ^ 0=\omega_2(\mathbf{k}_1) \\ k_2 ^ 0=\omega_{\lambda'}(\mathbf{k}_2)}}\,,\end{aligned}$$ 

 where the energy of the final electron is denoted as $ E_2$. Note the part by the group speed of the first and second polarization state (visualize besides Ref.   [ @Colladay:2001wk ]), respectively, where in the standard theory $ |\mathbf{v}_{\mathrm{gr},\,1}|=|\mathbf{v}_{\mathrm{gr},\,2}|=1$. The tensor $ (|\mathcal{M}(k_1,k_2)|^2\big)^{\mu\nu\varrho\sigma}$ is given by the trace condition of Eq. (5.81) of Ref.   [ @PeskinSchroeder1995 ] with some modifications due to the Lorentz - violating kinematics. 

! [ choose coordinate system for the phase space consolidation, where the initial photon three - momentum $ \mathbf{k}_1 $ lies along the third axis of the coordinate system. For the outgoing photon momentum $ \mathbf{k}_2 $, ball-shaped coordinates are chosen with the azimuthal slant $ \vartheta$ corresponding to the slant between $ \mathbf{k}_2 $ and the third bloc. Cases are treated with the three - vector $ \boldsymbol{\zeta}$ having equal or inadequate components.[]{data - label="fig: phase - space - integration"}](consistency - modmax - parity - odd_v1_fig4_coordinates.pdf) 

 The strictly algebraic part of the calculation, that include computation of tracing, contraction of indices and inserting kinematical sexual intercourse, is performed with 'Form '[ @Vermaseren:2000nd ]. The subsequent phase space calculation is done numerically with 'C++ `, since the resulting matrix chemical element squared contains hundreds of terms. The limit of zero Lorentz violation has to be taken with care and “ long double ” preciseness does not suffice here. consequently, the GNU Multiple Precision Arithmetic Library 'GMP '[ @GMP:2011 ] is used with its 'C++ 'interface described in Sec. 12 of the reference previously mention. 

 The first idea was to choose the coordinate system such that $ \boldsymbol{\zeta}$ lie down along the third axis. Then the phase space should have been integrated with cylindrical coordinates $ (k_{2,\bot},\varphi, k_{2,\parallel})$. To brood a general situation, where the initial photon momentum $ \mathbf{k}_1 $ points in an arbitrary
1^0=\omeha_2(\mathbf{k}_1) \\ k_2^0=\omega_{\lambda'}(\mxthbf{k}_2)}}\,,\end{aligneb}$$

qhere vhe enedgy of tfe final electron is denoted aw $E_2$. Nite the division by thd group vvlocity od tht first and seconv polaridction dtatz (wee also Ref. [@Cokladay:2001wk]), raspectively, whard nn the standard theory $|\mathbf{v}_{\mathrm{dr},\,1}|=|\mathbg{v}_{\lathrm{gr},\,2}|=1$. The tgnsor $(|\iathdal{M}(k_1,k_2)|^2\big)^{\mu\nu\varrho\sigma}$ is given gy the urace term of Eq. (5.81) pf Ref. [@PeskinSchroeder1995] with somf modifications duf to the Lotsnts-ciolating kivematics.

![Chosen coordinzte system for the phase space kntegxation, wherg fhf initial phiton nhree-momentum $\mathbf{k}_1$ lies akong the third axms od the coordinate systxm. For the outgoing khoton momancum $\mathbf{k}_2$, spherical ciorditatev ard chuseh xitg the wzijuthal angme $\vartheta$ corresponding to tne qngle between $\mathbs{k}_2$ and the third axis. Cases are treated wpth fhe three-vector $\boldsymvol{\zeta}$ having equal lr unequaj components.[]{data-label="fig:phase-space-integration"}](conshstenry-oodnaw-parkry-ldd_v1_fig4_coordinates.pdf)

The purely algebraic pare og nhe calculation, tmat includes compuyahipg of traces, cuntraccjoh of indices and ijsertind kinwmatical welayions, is performed with `Forn` [@Vermaseren:2000ud]. Rhe subsequent phade space cauculstion is done numerically wich `C++`, sjnce the redulting mzgrix element squxrec wontains hundreds of terms. The limiv of vero Lordntz violaeion has tl be bdken with care and “long dmuble” preclsion does not suffice here. Thecxfore, the GNU Mgltpple Precnsion Srithmetic Lifrary `GMP` [@GMP:2011] ns used cith igs `C++` interrace devcribed in Fec. 12 of the raverence prevmously megtiobed.

Tye firsg idea was to cnoose the coordinatw system such that $\boussymbol{\zeta}$ liew aoong the third sxir. Trej vhe prdse space shmuld haxr beev integrateb dith cylindrical coordindtes $(k_{2,\bot},\varphi,k_{2,\parallrl})$. To covet a generwl situation, ehere the initial khoton mooentuk $\mwthbf{k}_1$ points in an arbitrary
1^0=\omega_2(\mathbf{k}_1) \\ k_2^0=\omega_{\lambda'}(\mathbf{k}_2)}}\,,\end{aligned}$$ where the energy of electron denoted as Note the division the and second polarization (see also Ref. respectively, where in the standard theory The tensor $(|\mathcal{M}(k_1,k_2)|^2\big)^{\mu\nu\varrho\sigma}$ is given by the trace term of Eq. (5.81) of [@PeskinSchroeder1995] with some modifications due to the Lorentz-violating kinematics. ![Chosen coordinate system for phase integration, the photon three-momentum $\mathbf{k}_1$ lies along the third axis of the coordinate system. For the outgoing photon $\mathbf{k}_2$, spherical coordinates are chosen with the azimuthal $\vartheta$ corresponding to the between $\mathbf{k}_2$ and the third Cases treated with three-vector having or unequal components.[]{data-label="fig:phase-space-integration"}](consistency-modmax-parity-odd_v1_fig4_coordinates.pdf) purely algebraic part of the calculation, that includes computation of traces, contraction of indices and inserting kinematical is performed [@Vermaseren:2000nd]. The phase calculation done numerically with the resulting matrix element squared contains The limit of zero Lorentz violation has to taken with and “long double” precision does not here. Therefore, the GNU Multiple Precision Arithmetic Library [@GMP:2011] is used with its `C++` interface described in Sec. 12 of the reference previously first idea was to the coordinate system that lies the axis. Then phase space should have been integrated with cylindrical coordinates $(k_{2,\bot},\varphi,k_{2,\parallel})$. To a general situation, where the initial photon momentum $\mathbf{k}_1$ points arbitrary
1^0=\omega_2(\mathbf{k}_1) \\ k_2^0=\omega_{\lambda'}(\mAthbf{k}_2)}}\,,\end{aLigneD}$$

whEre ThE eneRgy oF the final electROn is Denoted as $E_2$. Note the divisIon by ThE GrouP VeLocitY of the fIRsT ANd sEcOnD poLaRIzAtion StaTe (see alSo Ref. [@CollaDay:2001Wk]), Respectively, WHeRe in the staNdaRd theory $|\mathBf{v}_{\Mathrm{Gr},\,1}|=|\MatHBf{v}_{\maThrM{gr},\,2}|=1$. ThE tensoR $(|\MathcaL{M}(k_1,k_2)|^2\big)^{\mu\Nu\VArrho\sIGma}$ is giVEN bY the Trace term of Eq. (5.81) of ReF. [@peSKinSchroeder1995] wiTh some MoDIfICAtiOns Due to the LoReNtz-viOLating kINeMATIcs.

![cHosen coordinaTe system for THe pHase spAcE inTEgratiOn, wheRe THe iNitial photoN thrEe-momentuM $\mathbF{K}_1$ lies alONg the thIrd axiS of The CoorDInAtE syStEM. FoR ThE ouTGoiNg photon MoMeNtum $\mAthbF{K}_2$, SPHeriCal CoorDinatEs are chosen wiTh tHe azIMutHal anGle $\vaRtheTa$ CorreSpondiNg to tHe Angle between $\matHbf{k}_2$ And the thiRd aXiS. CaSeS are tREated wIth The Three-veCtor $\bolDSymBoL{\ZETa}$ Having equal or unequAl COMpOnents.[]{daTa-labeL="FiG:pHAse-space-InTegRatiON"}](ConsiStenCY-mOdmax-parIty-odd_V1_FiG4_cOordinaTeS.pdf)

ThE pUreLy aLgebrAIc paRt of thE calculaTion, tHAt includes compUTation of traceS, CoNTRaCTion Of iNdices and inSertINg kiNemaTIcAl rELatioNs, is pErFOrMEd with `Form` [@VermasereN:2000nD]. The suBsequEnt phase space CalculatioN IS Done numeRicaLLy WIth `C++`, since the reSultiNg matrix elEMent squaRed coNtains huNdreds of tERMs. The limIt oF zeRo LOreNTZ vIolation has to BE TakeN wIth care And “Long douBle” PreCisIon DoEs not suffIce here. THeReFoRe, The gNU MuLTiple PreCiSioN ARitHmetiC librarY `GMP` [@GmP:2011] is UsEd WIth Its `C++` intERfACE desCrIbEd in sec. 12 Of The reFereNCe pReviousLy mentionEd.

THE firSt IdEa was to Choose the coorDiNate system SuCh tHat $\bolDSYmbol{\zetA}$ lies along the third axis. THEn the phAse Space ShouLd have beeN inTegratEd wITh cyliNdricaL coorDiNatES $(K_{2,\bot},\vARPhI,k_{2,\pArAllel})$. To covER A geNeral SiTuatIon, wherE the initial photon mOMenTum $\mathbf{k}_1$ poiNts In an ARBiTraRY
1^0=\omega_2(\mathbf{k}_1)  \\ k_2^0= \omeg a_{ \la mb da'} (\ma thbf{k}_2)}}\, , \end {aligned}$$

where the  ener gy of t h efinal  electr o ni s  de no te d a s$ E_ 2$. N ote  the di vision bythe  g roup velocit y  o f the firs t a nd second po lar izatio nsta t e (se e a lso R ef. [@ C ollada y:2001wk] ), respec t ively,w h er e in  the standard the o ry $|\mathbf{v}_{ \mathr m{ g r} , \ ,1} |=| \mathbf{v} _{ \math r m{gr},\ , 2} | = 1 $.T he tensor $(| \mathcal{M} ( k_1 ,k_2)| ^2 \bi g )^{\mu \nu\v ar r ho\ sigma}$ isgive n by thetracet erm ofE q. (5.8 1) ofRef . [ @Pes k in Sc hro ed e r19 9 5]  wi t h s ome modi fi ca tions  due t o theLor entz -viol ating kinemat ics .

! [ Cho sen c oordi nate  s ystem  for t he ph as e space integra tion , where t hein iti al  phot o n thre e-m ome ntum $\ mathbf{ k }_1 $l i e salong the third ax is o fthe coor dinate sy st e m. For t he  ou tgoi n g  phot on m o me ntum $\m athbf{ k }_ 2$ , spher ic al coo rd ina tes  arec hose n with  the azi mutha l  angle $\varth e ta$ correspon d in g to theang le between$\ma t hbf{ k}_2 $  a ndt he th ird a xi s .C ases are treated wi th  the t hree- vector $\bold symbol{\ze t a } $ having  equ a lo r unequal comp onent s.[]{data- l abel="fi g:pha se-space -integrat i o n"}](con sis ten cy- mod m a x- parity-odd_v1 _ f ig4_ co ordinat es. pdf)

T hepur ely  al ge braic par t of the  c al cu la tio n, th a t includ es  co mp uta tiono f trac es, c ontr ac ti o n o f indic e sa n d in se rt ingkin em atica l re l ati ons, is  performe d w i th ` Fo rm ` [@Ver maseren:2000n d] . The subs eq uen t phas e space ca lculation is done numer i cally w ith  `C++ `, s ince theres ulting  ma t rix el ementsquar ed  co n t ainsh u nd red sof terms.T h e l imitof  zer o Loren tz violation has t o  be  taken with c are  and “ lo ngd ou b le”  p r eci s i on does not suf fice here.  T h er efore, the GNU  M ultiple  Precis ion A r ithmeti c Library  `GMP` [@ GM P:20 1 1 ] i s used wit h its `C ++` inter f ace d e sc ribed  in  Sec.12  of  therefere n ceprevi ouslyme ntione d.

T he  first i dea was to choose the c oordin ate s yst em such t hat $\b oldsymbol {\ze ta}$ liesalo ngthe t hir d  axis . Th e nthe phase  spa c e shouldh av e b e e nintegratedw i t h c ylind ric a l coor dina tes $(k_{2,\bot}, \ varphi,k_{2,\p aral l e l}) $.T o co ve r a general si tua ti o n , whereth e initial p hoton mo me n tum $ \mathb f{k}_1 $ point s in an arb itra ry
1^0=\omega_2(\mathbf{k}_1) \\_k_2^0=\omega_{\lambda'}(\mathbf{k}_2)}}\,,\end{aligned}$$

where the_energy of the final_electron is_denoted_as $E_2$._Note_the division by_the group velocity_of the first and_second polarization state_(see_also Ref. [@Colladay:2001wk]), respectively, where in the standard theory $|\mathbf{v}_{\mathrm{gr},\,1}|=|\mathbf{v}_{\mathrm{gr},\,2}|=1$. The tensor $(|\mathcal{M}(k_1,k_2)|^2\big)^{\mu\nu\varrho\sigma}$ is given_by_the trace_term_of_Eq. (5.81) of Ref. [@PeskinSchroeder1995] with_some modifications due to the_Lorentz-violating kinematics.

![Chosen_coordinate system for the phase space integration, where_the_initial photon three-momentum_$\mathbf{k}_1$ lies along the third axis of the coordinate_system. For the outgoing photon momentum_$\mathbf{k}_2$, spherical coordinates_are_chosen_with the azimuthal angle_$\vartheta$ corresponding to the angle between_$\mathbf{k}_2$ and the third axis. Cases_are treated with the three-vector $\boldsymbol{\zeta}$ having_equal or unequal components.[]{data-label="fig:phase-space-integration"}](consistency-modmax-parity-odd_v1_fig4_coordinates.pdf)

The purely algebraic_part of the calculation, that_includes computation_of traces, contraction of indices_and inserting kinematical_relations, is_performed with `Form`_[@Vermaseren:2000nd]. The subsequent phase space calculation_is done numerically_with `C++`, since the resulting matrix_element_squared contains hundreds_of_terms._The limit_of zero Lorentz_violation_has to_be_taken with care and “long double”_precision_does not suffice here. Therefore, the GNU_Multiple Precision Arithmetic Library_`GMP`_[@GMP:2011] is used with_its `C++` interface described in_Sec. 12 of the reference previously_mentioned.

The first_idea was_to choose the coordinate system such that $\boldsymbol{\zeta}$ lies along the_third axis. Then the phase space_should have been integrated_with cylindrical_coordinates_$(k_{2,\bot},\varphi,k_{2,\parallel})$. To cover_a_general situation,_where the initial photon momentum $\mathbf{k}_1$ points_in an_arbitrary