IAP GITLAB

Commit decf3bad authored by Tanguy Pierog's avatar Tanguy Pierog

introduce the possibility to switch off Breit-Wigner mass broadening (default...

introduce the possibility to switch off Breit-Wigner mass broadening (default is off, but after some test the difference is very small ... excepte for resonances of course !)


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4473 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 2ef89b55
......@@ -187,6 +187,7 @@ c fragmentation and decay options
!also several "1"s may be used with obvious meaning
maxres=99999 !maximum resonance spin
aouni=0. !technical parameter for development
ibreit=0 !Breit-Wigner mass broadening (0=off, 1=on)
c lepton-nucleon and e+e- options
......@@ -4242,6 +4243,7 @@ c printout options
if(linex(ix:jx).eq.'modsho')modsho=nint(val)
if(linex(ix:jx).eq.'idensi')idensi=nint(val)
if(linex(ix:jx).eq.'infragm')infragm=nint(val)
if(linex(ix:jx).eq.'ibreit')ibreit=nint(val)
if(linex(ix:jx).eq.'ishevt')ishevt=nint(val)
if(linex(ix:jx).eq.'iwseed')iwseed=nint(val)
if(linex(ix:jx).eq.'jwseed')jwseed=nint(val)
......
......@@ -332,6 +332,7 @@ c Decay of alpha
goto 111
endif
if(ibreit.eq.1)then
c Introduce mass broadening due to life time (Breit-Wigner distribution)
ctau=tivptl(2,ip)-tivptl(1,ip) !ctau*gamma (fm)
c ctau=ctau*pptl(5,ip)/pptl(4,ip) !ctau (fm)
......@@ -353,6 +354,7 @@ c only if mean mass is used
c write(*,*)'broadening:',ip,idlv1,abs(pptl(5,ip)-amss)/wi,amss,ntry
endif
endif
endif
c select one of the decay channel
ipoint=look(iabs(idlv1))-1
......
......@@ -1328,6 +1328,7 @@ c initializes /crema/
c masses are limit between stable state (so the average between 2 mass states)
c width=hbar(6.582e-25 GeV.s)/tau for 151, 251, 351, 451 arbitrary
c (no data found) !!!!!!!!!!!
c used only in epos-ids.f
c-----------------------------------------------------------------------
parameter (mxindx=1000,mxre=100,mxma=11,mxmx=6)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment