Page MenuHomeHEPForge
Files F11221343

HdecayWrap.f
No OneTemporary
Actions

HdecayWrap.f
View Options

subroutine hdecaycall2hdm(ABRB_out,ABRL_out,ABRM_out,ABRS_out,
&ABRC_out,
&ABRT_out,ABRG_out,ABRGA_out, ABRZGA_out,ABRZ_out,AWDTH_out,
&HLBRB_out,HLBRL_out,HLBRM_out,HLBRS_out,HLBRC_out,HLBRT_out,
&HLBRG_out,HLBRGA_out,HLBRZGA_out,HLBRW_out,HLBRZ_out,HLBRA_out,
&HLBRAZ_out,HLBRHW_out,HLWDTH_out,HHBRB_out,HHBRL_out,HHBRM_out,
&HHBRS_out,HHBRC_out,HHBRT_out,HHBRG_out,HHBRGA_out,HHBRZGA_out,
&HHBRW_out,HHBRZ_out,HHBRH_out,HHBRA_out,HHBRAZ_out,HHBRHW_out,
&HHWDTH_out,HCBRB_out,HCBRL_out,HCBRM_out,HCBRBU_out,HCBRS_out,
&HCBRC_out,HCBRT_out,HCBRW_out,HCBRA_out,HCWDTH_out,hcbrcd_out,
&hcbrts_out,hcbrtd_out,hcbrwhh_out,hhbrchch_out,hlbrchch_out,
&abrhhaz_out,abrhawphm_out,gamt0_out,gamt1_out)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
PARAMETER(K=6,NI=87,NSA=85,NSB=86,NLA=88,NLB=89,NHA=90,NHB=91,
. NHC=92,NAA=93,NAB=94,NCA=95,NCB=96,NCC=50,NRA=97,NRB=98,
. NSUSYL=81,NSUSYA=82,NSUSYH=83,NSUSYC=84,NPAR=80,
. NSUSYLA=79,NSUSYLB=78,NSUSYLC=77,NSUSYLD=76,NSUSYLE=75,
. NSUSYLF=59,NSUSYHF=58,
. NSUSYHA=74,NSUSYHB=73,NSUSYHC=72,NSUSYHD=71,NSUSYHE=70,
. NSUSYAA=69,NSUSYAB=68,NSUSYAC=67,NSUSYAD=66,NSUSYAE=65,
. NSUSYCA=64,NSUSYCB=63,NSUSYCC=62,NSUSYCD=61,NSUSYCE=60,
. ninlha=22)
double precision minval(1:20),smval(1:30),massval(1:50),
. nmixval(4,4),umixval(2,2),vmixval(2,2),
. stopmixval(2,2),sbotmixval(2,2),staumixval(2,2),
. hmixval(1:10),gaugeval(1:3),msoftval(1:100),
. auval(3,3),adval(3,3),aeval(3,3),yuval(3,3),
. ydval(3,3),yeval(3,3),qvalue(1:20),extval(0:100),
. m_softval(1:100)
double precision slhaneut(1:4),slhaxneut(1:4),slhachar(1:2),
. slhaxchar(1:2),
. slhau(2,2),slhav(2,2),slhaz(4,4),
. slhast(2),slhasb(2),slhasu(2),slhasd(2),
. slhase(2),slhasl(2),slhasn(2),slhasnl(2),
. warning(1:10)
double precision vckmval(4)
integer imod(1:2)
integer check(1:22)
double precision mbmsbar,mbl,mbu
character spinfo1*100,spinfo2*100,modselval*100,mincom(1:20)*20,
. extcom(0:100)*20,softcom(1:100)*20,hmixcom(1:10)*20,
. m_softcom(1:100)*20
DIMENSION GMN(4),XMN(4),GMC(2),GMST(2),GMSB(2),GMSL(2),
. GMSU(2),GMSD(2),GMSE(2),GMSN(2),GMSN1(2)
DIMENSION HLBRSC(2,2),HLBRSN(4,4),HHBRSC(2,2),HHBRSN(4,4),
. HABRSC(2,2),HABRSN(4,4),HCBRSU(2,4),
. HHBRST(2,2),HHBRSB(2,2),HCBRSTB(2,2)
DIMENSION AC1(2,2),AC2(2,2),AC3(2,2),
. AN1(4,4),AN2(4,4),AN3(4,4),
. ACNL(2,4),ACNR(2,4)
DIMENSION GLTT(2,2),GLBB(2,2),GHTT(2,2),GHBB(2,2),GCTB(2,2),
. GLEE(2,2),GHEE(2,2),GCEN(2,2)
DIMENSION AGDL(4),AGDA(4),AGDH(4),AGDC(2)
c -------------- common block given by read_leshouches ------------ c
COMMON/SLHA_leshouches1_HDEC/spinfo1,spinfo2,modselval,mincom,
. extcom,softcom,hmixcom
COMMON/SLHA_leshouches2_HDEC/minval,extval,smval,massval,nmixval,
. umixval,vmixval,stopmixval,sbotmixval,
. staumixval,hmixval,gaugeval,msoftval,auval,
. adval,aeval,yuval,ydval,yeval,alphaval,
. qvalue,imod
COMMON/SLHA_leshouches3_HDEC/vckmval
COMMON/SD_scaleofewsb/scaleofewsb
c -------------- common blocks needed in HDECAY subroutines ---------- c
COMMON/SLHA_vals_HDEC/islhai,islhao
COMMON/SLHA_m1_HDEC/am1
COMMON/SLHA_gaug_HDEC/slhaneut,slhaxneut,slhachar,slhau,slhav,
. slhaz,slhaxchar
COMMON/SLHA_sfer_HDEC/slhast,slhasb,slhasu,slhasd,slhase,slhasl,
. slhasn,slhasnl,slhacot,slhasit,slhacob,slhasib,
. slhacol,slhasil
COMMON/SLHA_hmass_HDEC/slhaml,slhamh,slhamc,slha_alpha
COMMON/SLHAVAL_HDEC/g1ew,g2ew
COMMON/SLHA_checkval_HDEC/check
COMMON/MASSES_HDEC/AMS,AMC,AMB,AMT
COMMON/STRANGE_HDEC/AMSB
COMMON/PARAM_HDEC/GF,ALPH,AMTAU,AMMUON,AMZ,AMW
COMMON/CKMPAR_HDEC/VTB,VTS,VTD,VCB,VCS,VCD,VUB,VUS,VUD
COMMON/HMASS_HDEC/AMSM,AMA,AML,AMH,AMCH,AMAR
COMMON/BREAKSCALE_HDEC/SUSYSCALE
COMMON/BREAK_HDEC/AMEL,AMER,AMSQ,AMUR,AMDR,AL,AU,AD,AMU,AM2
COMMON/BREAKGLU_HDEC/AMGLU
COMMON/SFER1ST_HDEC/AMQL1,AMUR1,AMDR1,AMEL1,AMER1
COMMON/GLUINO_HDEC/AMGLUINO,XMSB1,XMSB2,STHB,CTHB,
. XLBB(2,2),XHBB(2,2),XABB(2,2),
. XMST1,XMST2,STHT,CTHT,
. XLTT(2,2),XHTT(2,2),XATT(2,2)
COMMON/WZWDTH_HDEC/GAMC0,GAMT0,GAMT1,GAMW,GAMZ
COMMON/COUP_HDEC/GAT,GAB,GLT,GLB,GHT,GHB,GZAH,GZAL,
. GHHH,GLLL,GHLL,GLHH,GHAA,GLAA,GLVV,GHVV,
. GLPM,GHPM,B,A
COMMON/ALS_HDEC/XLAMBDA,AMC0,AMB0,AMT0,N0
COMMON/FLAG_HDEC/IHIGGS,NNLO,IPOLE
COMMON/SM4_HDEC/AMTP,AMBP,AMNUP,AMEP,ISM4,IGGELW
COMMON/FERMIOPHOBIC_HDEC/IFERMPHOB
COMMON/MODEL_HDEC/IMODEL
COMMON/ONSHELL_HDEC/IONSH,IONWZ,IOFSUSY
COMMON/OLDFASH_HDEC/NFGG
COMMON/WIDTHSM_HDEC/SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT,SMBRG,
. SMBRGA,SMBRZGA,SMBRW,SMBRZ,SMWDTH
COMMON/WIDTHA_HDEC/ABRB,ABRL,ABRM,ABRS,ABRC,ABRT,ABRG,ABRGA,
. ABRZGA,ABRZ,AWDTH
COMMON/WIDTHHL_HDEC/HLBRB,HLBRL,HLBRM,HLBRS,HLBRC,HLBRT,HLBRG,
. HLBRGA,HLBRZGA,HLBRW,HLBRZ,HLBRA,HLBRAZ,HLBRHW,
. HLWDTH
COMMON/WIDTHHH_HDEC/HHBRB,HHBRL,HHBRM,HHBRS,HHBRC,HHBRT,HHBRG,
. HHBRGA,HHBRZGA,HHBRW,HHBRZ,HHBRH,HHBRA,HHBRAZ,
. HHBRHW,HHWDTH
COMMON/WIDTHHC_HDEC/HCBRB,HCBRL,HCBRM,HCBRBU,HCBRS,HCBRC,HCBRT,
. HCBRW,HCBRA,HCWDTH
COMMON/WISUSY_HDEC/HLBRSC,HLBRSN,HHBRSC,HHBRSN,HABRSC,HABRSN,
. HCBRSU,HLBRCHT,HHBRCHT,HABRCHT,HLBRNET,HHBRNET,
. HABRNET,HCBRCNT,HLBRSL,HHBRSL,HCBRSL,HABRSL,HABRST,
. HABRSB,HHBRSQ,HHBRST,HHBRSB,HHBRSQT,HCBRSQ,HCBRSTB,
. HCBRSQT,HLBRSQ,HLBRSQT
COMMON/WISFER_HDEC/BHLSLNL,BHLSLEL,BHLSLER,BHLSQUL,BHLSQUR,
. BHLSQDL,BHLSQDR,BHLST(2,2),BHLSB(2,2),BHLSTAU(2,2),
. BHHSLNL,BHHSLEL,BHHSLER,BHHSQUL,BHHSQUR,BHHSQDL,
. BHHSQDR,BHHST(2,2),BHHSB(2,2),BHHSTAU(2,2),
. BHASTAU,BHASB,BHAST,
. BHCSL00,BHCSL11,BHCSL21,BHCSQ,BHCSTB(2,2)
COMMON/SMASS_HDEC/GMN,XMN,GMC,GMST,GMSB,GMSL,GMSU,GMSD,GMSE,GMSN
. ,GMSN1
COMMON/GOLDST_HDEC/AXMPL,AXMGD,IGOLD
COMMON/WIGOLD_HDEC/HLBRGD,HABRGD,HHBRGD,HCBRGD
COMMON/FLAGS_HDEC/INDIDEC
COMMON/CPSM_HDEC/CPW,CPZ,CPTAU,CPMU,CPT,CPB,CPC,CPS,
. CPGAGA,CPGG,CPZGA,ICOUPELW
c MMM changed 21/8/13
COMMON/THDM_HDEC/TGBET2HDM,ALPH2HDM,AMHL2HDM,AMHH2HDM,
. AMHA2HDM,AMHC2HDM,AM12SQ,ITYPE2HDM,I2HDM
COMMON/WIDTH_HC_ADD/hcbrcd,hcbrts,hcbrtd
COMMON/WIDTH_2HDM/hcbrwhh,hhbrchch,hlbrchch,abrhhaz,abrhawphm
c end MMM changed 21/8/13
EXTERNAL READ_HDEC,HEAD_HDEC,HDEC,CLOSE_HDEC
c double precision totalwidth
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c COMMON/DAVID/QSUSY1,QSUSY2,LOOP
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
CALL READ_HDEC(TGBET,AMABEG,AMAEND,NMA)
if(islhao.ne.1) then
CALL HEAD_HDEC(TGBET,AMABEG)
endif
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c write(6,*)'Loop, Factor_QCD, Factor_elw = ?'
c read(5,*)LOOP,QSUSY1,QSUSY2
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c write(*,*)TGBET2HDM,"----------"
DO 9999 II=1,NMA
IF(NMA.NE.1)THEN
AMAR = AMABEG + (AMAEND-AMABEG)/(NMA-1D0)*(II-1D0)
ELSE
AMAR = AMABEG
ENDIF
AMSM = AMAR
AMA = AMAR
CALL HDEC(TGBET)
c CALL WRITE_HDEC(TGBET)
c totalwidth=HLWDTH
9999 CONTINUE
ABRB_out=ABRB
ABRL_out=ABRL
ABRM_out=ABRM
ABRS_out=ABRS
ABRC_out=ABRC
ABRT_out=ABRT
ABRG_out=ABRG
ABRGA_out=ABRGA
ABRZGA_out=ABRZGA
ABRZ_out=ABRZ
AWDTH_out=AWDTH
HLBRB_out=HLBRB
HLBRL_out=HLBRL
HLBRM_out=HLBRBM
HLBRS_out=HLBRA
HLBRC_out=HLBRC
HLBRT_out=HLBRT
HLBRG_out=HLBRG
HLBRGA_out=HLBRGA
HLBRZGA_out=HLBRZGA
HLBRW_out=HLBRW
HLBRZ_out=HLBRZ
HLBRA_out=HLBRA
HLBRAZ_out=HLBRAZ
HLBRHW_out=HLBRHW
HLWDTH_out=HLWDTH
HHBRB_out=HHBRB
HHBRL_out=HHBRL
HHBRM_out=HHBRM
HHBRS_out=HHBRS
HHBRC_out=HHBRC
HHBRT_out=HHBRT
HHBRG_out=HHBRG
HHBRGA_out=HHBRGA
HHBRZGA_out=HHBRZGA
HHBRW_out=HHBRW
HHBRZ_out=HHBRZ
HHBRH_out=HHBRH
HHBRA_out=HHBRA
HHBRAZ_out=HHBRAZ
HHBRHW_out=HHBRHW
HHWDTH_out=HHWDTH
HCBRB_out=HCBRB
HCBRL_out=HCBRL
HCBRM_out=HCBRBM
HCBRBU_out=HCBRBU
HCBRS_out=HCBRS
HCBRC_out=HCBRC
HCBRT_out=HCBRT
HCBRW_out=HCBRW
HCBRA_out=HCBRA
HCWDTH_out=HCWDTH
hcbrcd_out=hcbrcd
hcbrts_out=hcbrts
hcbrtd_out=hcbrtd
hcbrwhh_out=hcbrwhh
hhbrchch_out=hhbrchch
hlbrchch_out=hlbrchch
abrhhaz_out=abrhhaz
abrhawphm_out=abrhawphm
gamt0_out=gamt0
gamt1_out=gamt1
CALL CLOSE_HDEC
END
subroutine hdecaycallsm(SMBRB_out,SMBRL_out,SMBRM_out,SMBRS_out,
&SMBRC_out,SMBRT_out,SMBRG_out,SMBRGA_out,SMBRZGA_out,SMBRW_out,
&SMBRZ_out,SMWDTH_out)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
PARAMETER(K=6,NI=87,NSA=85,NSB=86,NLA=88,NLB=89,NHA=90,NHB=91,
. NHC=92,NAA=93,NAB=94,NCA=95,NCB=96,NCC=50,NRA=97,NRB=98,
. NSUSYL=81,NSUSYA=82,NSUSYH=83,NSUSYC=84,NPAR=80,
. NSUSYLA=79,NSUSYLB=78,NSUSYLC=77,NSUSYLD=76,NSUSYLE=75,
. NSUSYLF=59,NSUSYHF=58,
. NSUSYHA=74,NSUSYHB=73,NSUSYHC=72,NSUSYHD=71,NSUSYHE=70,
. NSUSYAA=69,NSUSYAB=68,NSUSYAC=67,NSUSYAD=66,NSUSYAE=65,
. NSUSYCA=64,NSUSYCB=63,NSUSYCC=62,NSUSYCD=61,NSUSYCE=60,
. ninlha=22)
double precision minval(1:20),smval(1:30),massval(1:50),
. nmixval(4,4),umixval(2,2),vmixval(2,2),
. stopmixval(2,2),sbotmixval(2,2),staumixval(2,2),
. hmixval(1:10),gaugeval(1:3),msoftval(1:100),
. auval(3,3),adval(3,3),aeval(3,3),yuval(3,3),
. ydval(3,3),yeval(3,3),qvalue(1:20),extval(0:100),
. m_softval(1:100)
double precision slhaneut(1:4),slhaxneut(1:4),slhachar(1:2),
. slhaxchar(1:2),
. slhau(2,2),slhav(2,2),slhaz(4,4),
. slhast(2),slhasb(2),slhasu(2),slhasd(2),
. slhase(2),slhasl(2),slhasn(2),slhasnl(2),
. warning(1:10)
double precision vckmval(4)
integer imod(1:2)
integer check(1:22)
double precision mbmsbar,mbl,mbu
character spinfo1*100,spinfo2*100,modselval*100,mincom(1:20)*20,
. extcom(0:100)*20,softcom(1:100)*20,hmixcom(1:10)*20,
. m_softcom(1:100)*20
DIMENSION GMN(4),XMN(4),GMC(2),GMST(2),GMSB(2),GMSL(2),
. GMSU(2),GMSD(2),GMSE(2),GMSN(2),GMSN1(2)
DIMENSION HLBRSC(2,2),HLBRSN(4,4),HHBRSC(2,2),HHBRSN(4,4),
. HABRSC(2,2),HABRSN(4,4),HCBRSU(2,4),
. HHBRST(2,2),HHBRSB(2,2),HCBRSTB(2,2)
DIMENSION AC1(2,2),AC2(2,2),AC3(2,2),
. AN1(4,4),AN2(4,4),AN3(4,4),
. ACNL(2,4),ACNR(2,4)
DIMENSION GLTT(2,2),GLBB(2,2),GHTT(2,2),GHBB(2,2),GCTB(2,2),
. GLEE(2,2),GHEE(2,2),GCEN(2,2)
DIMENSION AGDL(4),AGDA(4),AGDH(4),AGDC(2)
c -------------- common block given by read_leshouches ------------ c
COMMON/SLHA_leshouches1_HDEC/spinfo1,spinfo2,modselval,mincom,
. extcom,softcom,hmixcom
COMMON/SLHA_leshouches2_HDEC/minval,extval,smval,massval,nmixval,
. umixval,vmixval,stopmixval,sbotmixval,
. staumixval,hmixval,gaugeval,msoftval,auval,
. adval,aeval,yuval,ydval,yeval,alphaval,
. qvalue,imod
COMMON/SLHA_leshouches3_HDEC/vckmval
COMMON/SD_scaleofewsb/scaleofewsb
c -------------- common blocks needed in HDECAY subroutines ---------- c
COMMON/SLHA_vals_HDEC/islhai,islhao
COMMON/SLHA_m1_HDEC/am1
COMMON/SLHA_gaug_HDEC/slhaneut,slhaxneut,slhachar,slhau,slhav,
. slhaz,slhaxchar
COMMON/SLHA_sfer_HDEC/slhast,slhasb,slhasu,slhasd,slhase,slhasl,
. slhasn,slhasnl,slhacot,slhasit,slhacob,slhasib,
. slhacol,slhasil
COMMON/SLHA_hmass_HDEC/slhaml,slhamh,slhamc,slha_alpha
COMMON/SLHAVAL_HDEC/g1ew,g2ew
COMMON/SLHA_checkval_HDEC/check
COMMON/MASSES_HDEC/AMS,AMC,AMB,AMT
COMMON/STRANGE_HDEC/AMSB
COMMON/PARAM_HDEC/GF,ALPH,AMTAU,AMMUON,AMZ,AMW
COMMON/CKMPAR_HDEC/VTB,VTS,VTD,VCB,VCS,VCD,VUB,VUS,VUD
COMMON/HMASS_HDEC/AMSM,AMA,AML,AMH,AMCH,AMAR
COMMON/BREAKSCALE_HDEC/SUSYSCALE
COMMON/BREAK_HDEC/AMEL,AMER,AMSQ,AMUR,AMDR,AL,AU,AD,AMU,AM2
COMMON/BREAKGLU_HDEC/AMGLU
COMMON/SFER1ST_HDEC/AMQL1,AMUR1,AMDR1,AMEL1,AMER1
COMMON/GLUINO_HDEC/AMGLUINO,XMSB1,XMSB2,STHB,CTHB,
. XLBB(2,2),XHBB(2,2),XABB(2,2),
. XMST1,XMST2,STHT,CTHT,
. XLTT(2,2),XHTT(2,2),XATT(2,2)
COMMON/WZWDTH_HDEC/GAMC0,GAMT0,GAMT1,GAMW,GAMZ
COMMON/COUP_HDEC/GAT,GAB,GLT,GLB,GHT,GHB,GZAH,GZAL,
. GHHH,GLLL,GHLL,GLHH,GHAA,GLAA,GLVV,GHVV,
. GLPM,GHPM,B,A
COMMON/ALS_HDEC/XLAMBDA,AMC0,AMB0,AMT0,N0
COMMON/FLAG_HDEC/IHIGGS,NNLO,IPOLE
COMMON/SM4_HDEC/AMTP,AMBP,AMNUP,AMEP,ISM4,IGGELW
COMMON/FERMIOPHOBIC_HDEC/IFERMPHOB
COMMON/MODEL_HDEC/IMODEL
COMMON/ONSHELL_HDEC/IONSH,IONWZ,IOFSUSY
COMMON/OLDFASH_HDEC/NFGG
COMMON/WIDTHSM_HDEC/SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT,SMBRG,
. SMBRGA,SMBRZGA,SMBRW,SMBRZ,SMWDTH
COMMON/WIDTHA_HDEC/ABRB,ABRL,ABRM,ABRS,ABRC,ABRT,ABRG,ABRGA,
. ABRZGA,ABRZ,AWDTH
COMMON/WIDTHHL_HDEC/HLBRB,HLBRL,HLBRM,HLBRS,HLBRC,HLBRT,HLBRG,
. HLBRGA,HLBRZGA,HLBRW,HLBRZ,HLBRA,HLBRAZ,HLBRHW,
. HLWDTH
COMMON/WIDTHHH_HDEC/HHBRB,HHBRL,HHBRM,HHBRS,HHBRC,HHBRT,HHBRG,
. HHBRGA,HHBRZGA,HHBRW,HHBRZ,HHBRH,HHBRA,HHBRAZ,
. HHBRHW,HHWDTH
COMMON/WIDTHHC_HDEC/HCBRB,HCBRL,HCBRM,HCBRBU,HCBRS,HCBRC,HCBRT,
. HCBRW,HCBRA,HCWDTH
COMMON/WISUSY_HDEC/HLBRSC,HLBRSN,HHBRSC,HHBRSN,HABRSC,HABRSN,
. HCBRSU,HLBRCHT,HHBRCHT,HABRCHT,HLBRNET,HHBRNET,
. HABRNET,HCBRCNT,HLBRSL,HHBRSL,HCBRSL,HABRSL,HABRST,
. HABRSB,HHBRSQ,HHBRST,HHBRSB,HHBRSQT,HCBRSQ,HCBRSTB,
. HCBRSQT,HLBRSQ,HLBRSQT
COMMON/WISFER_HDEC/BHLSLNL,BHLSLEL,BHLSLER,BHLSQUL,BHLSQUR,
. BHLSQDL,BHLSQDR,BHLST(2,2),BHLSB(2,2),BHLSTAU(2,2),
. BHHSLNL,BHHSLEL,BHHSLER,BHHSQUL,BHHSQUR,BHHSQDL,
. BHHSQDR,BHHST(2,2),BHHSB(2,2),BHHSTAU(2,2),
. BHASTAU,BHASB,BHAST,
. BHCSL00,BHCSL11,BHCSL21,BHCSQ,BHCSTB(2,2)
COMMON/SMASS_HDEC/GMN,XMN,GMC,GMST,GMSB,GMSL,GMSU,GMSD,GMSE,GMSN
. ,GMSN1
COMMON/GOLDST_HDEC/AXMPL,AXMGD,IGOLD
COMMON/WIGOLD_HDEC/HLBRGD,HABRGD,HHBRGD,HCBRGD
COMMON/FLAGS_HDEC/INDIDEC
COMMON/CPSM_HDEC/CPW,CPZ,CPTAU,CPMU,CPT,CPB,CPC,CPS,
. CPGAGA,CPGG,CPZGA,ICOUPELW
c MMM changed 21/8/13
COMMON/THDM_HDEC/TGBET2HDM,ALPH2HDM,AMHL2HDM,AMHH2HDM,
. AMHA2HDM,AMHC2HDM,AM12SQ,ITYPE2HDM,I2HDM
COMMON/WIDTH_HC_ADD/hcbrcd,hcbrts,hcbrtd
COMMON/WIDTH_2HDM/hcbrwhh,hhbrchch,hlbrchch,abrhhaz,abrhawphm
c end MMM changed 21/8/13
EXTERNAL READ_HDEC,HEAD_HDEC,HDEC,CLOSE_HDEC
c double precision totalwidth
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c COMMON/DAVID/QSUSY1,QSUSY2,LOOP
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
CALL READ_HDEC(TGBET,AMABEG,AMAEND,NMA)
if(islhao.ne.1) then
CALL HEAD_HDEC(TGBET,AMABEG)
endif
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c write(6,*)'Loop, Factor_QCD, Factor_elw = ?'
c read(5,*)LOOP,QSUSY1,QSUSY2
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c write(*,*)TGBET2HDM,"----------"
DO 9999 II=1,NMA
IF(NMA.NE.1)THEN
AMAR = AMABEG + (AMAEND-AMABEG)/(NMA-1D0)*(II-1D0)
ELSE
AMAR = AMABEG
ENDIF
AMSM = AMAR
AMA = AMAR
CALL HDEC(TGBET)
c CALL WRITE_HDEC(TGBET)
c totalwidth=HLWDTH
9999 CONTINUE
SMBRB_out=SMBRB
SMBRL_out=SMBRL
SMBRM_out=SMBRM
SMBRS_out=SMBR
SMBRC_out=SMBRC
SMBRT_out=SMBRT
SMBRG_out=SMBRG
SMBRGA_out=SMBRGA
SMBRZGA_out=SMBRZGA
SMBRW_out=SMBRW
SMBRZ_out=SMBRZ
SMWDTH_out=SMWDTH
CALL CLOSE_HDEC
END
subroutine hdecaycallplain(dummy)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
EXTERNAL READ_HDEC,HEAD_HDEC,HDEC,CLOSE_HDEC
c double precision totalwidth
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c COMMON/DAVID/QSUSY1,QSUSY2,LOOP
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
CALL READ_HDEC(TGBET,AMABEG,AMAEND,NMA)
if(islhao.ne.1) then
CALL HEAD_HDEC(TGBET,AMABEG)
endif
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c write(6,*)'Loop, Factor_QCD, Factor_elw = ?'
c read(5,*)LOOP,QSUSY1,QSUSY2
c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
c write(*,*)TGBET2HDM,"----------"
DO 9999 II=1,NMA
IF(NMA.NE.1)THEN
AMAR = AMABEG + (AMAEND-AMABEG)/(NMA-1D0)*(II-1D0)
ELSE
AMAR = AMABEG
ENDIF
AMSM = AMAR
AMA = AMAR
CALL HDEC(TGBET)
c CALL WRITE_HDEC(TGBET)
c totalwidth=HLWDTH
9999 CONTINUE
CALL CLOSE_HDEC
END

File Metadata

Mime Type
text/x-fortran
Expires
Wed, May 14, 10:17 AM (1 d, 11 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5111088
Default Alt Text
HdecayWrap.f (18 KB)

Event Timeline