program example call ltdemo end *CMZ : 03/02/95 12.26.53 by A.Rozanov *-- Author : A.Rozanov 01/02/95 SUBROUTINE LTDEMO * * *----------------------------------------------------------------------* * * * Name : LTDEMO * * (module) * * * * Description : * * template program to demonstrate LEPTOP package * * * * * * * *----------------------------------------------------------------------* * * * Author : A.Rozanov Date : 01/02/95 * * * * * * Last modifications : * * [name] Date : [date] * * * * Keywords : * * LTINIT, LTGET, LTPUT, LTFLAG * *----------------------------------------------------------------------* * * write(6,*)'start ltdemo' * set all LEPTOP print flags off CALL LTFLAG('PRNT',0) * initialisation CALL LTINIT(0) * get some constants CALL LTGET('MZ',AMZ) CALL LTGET('GFERMI',GFERMI) CALL LTGET('ALBAR',ALBAR) * write(6,'(' MZ= ',F8.4,'GFu= ',E16.8,'alphabar= 1./',F6.2)') * 1 AMZ,GFERMI,1./ALBAR * get more constants CALL LTGET('MELE',amele) CALL LTGET('MMUO',ammuo) CALL LTGET('MTAU',amtau) CALL LTGET('MS',amstr) CALL LTGET('MC',amchrm) CALL LTGET('MB',ambot) * write(6,10010)amele,ammuo,amtau,amstr,amchrm,ambot *10010 format(' me= ',f10.6,' mmu= ',f6.3,' mtau= ',f6.3, * 1 ' ms= ',f6.3,' mc= ',f6.3,' mb= ',f6.3) * modify some constants CALL LTPUT('MT',175.) CALL LTPUT('MH',300.) CALL LTPUT('ALSHAT',0.125) * get some physical quantities CALL LTGET('GV',gv) CALL LTGET('GA',ga) CALL LTGET('MW',amw) CALL LTGET('GNU',gnu) CALL LTGET('GE',ge) CALL LTGET('GMUON',gmuon) CALL LTGET('GTAU',gtau) write(6,*)'+LTDEMO: gv, ga, amw', 1 gv, ga, amw write(6,*)'+LTDEMO: gnu,ge,gmuon,gtau', 1 gnu,ge,gmuon,gtau * test fit of mtop and alsbar CALL LTFIT2('MT,ALS',amt,als,eamt,eals,rho,chi2,ndf) write(6,*)'+LTDEMO: mt= ',amt,' +- ',eamt, 1 ' als= ',als,' +- ',eals write(6,*)'+LTDEMO: rho,chi2= ',rho,chi2 *-----get data from Lepton Photon 99 conference CALL LTFPUT('LP99','ALL',0.) write(6,*)'+LTDEMO: read Lepton-Photon 99 data ' *----- modify value of Z total width CALL LTFPUT('VALUE','GZ',2.4974) write(6,*)'+LTDEMO: put GZ=2.4974 ' *-----test fit of mtop and mh CALL LTFIT2('MT,MH',amt,amh,eamt,eamh,rho,chi2,ndf) write(6,*)'+LTDEMO: mt= ',amt,' +- ',eamt, 1 ' amh= ',amh,' +- ',eamh write(6,*)'+LTDEMO: rho,chi2= ',rho,chi2 *-----test fit to mt CALL LTFIT1('MT',amt,eamt,chi2,ndf) write(6,*)'+LTDEMO: mt= ',amt,' +- ',eamt write(6,*)'+LTDEMO: chi2= ',chi2 *-----test fit to mh CALL LTFIT1('MH',amh,eamh,chi2,ndf) write(6,*)'+LTDEMO:', 1 ' mh= ',amh,' +- ',eamh write(6,*)'+LTDEMO: chi2= ',chi2 *-----test fit to als CALL LTFIT1('ALS',als,eals,chi2,ndf) write(6,*)'+LTDEMO: als= ',als,' +- ',eals write(6,*)'+LTDEMO: chi2= ',chi2 END ! LTDEMO