Page MenuHomeHEPForge
Files F11221385

notes_on_comparing_LEP_tables_type1.txt
No OneTemporary
Actions

notes_on_comparing_LEP_tables_type1.txt
View Options

*********************************************************
******** to compare LEP tables type 1 ****************
30.03.09 **NEW METHOD*** (due to changes in structure of Oliver's code)
* Add the following to subroutine check_LEP_bound (just before "comparing input to the predicted LEP results"):
!KW added
open(10,file="debugging_tables_leptype1.txt")
predicted_limit = .TRUE.
n=1
do k=0,nLEP_NSneutCHN-1,1
write(10,'(" (pred) id of table:",1I4)')k+1
runid=n+NHIGGS*k
call check_LEP_bound1(Mh_avg(runid),
& 1.0D0,predratio(runid),k+1,predicted_limit)
enddo
n=1
do k=0,nLEP_NScharCHN-1,1
write(10,'(" (pred) id of table:",1I4)')nLEP_NSneutCHN+k+1
runid=n+NHPLUS*k+nLEP_NSneutCHN*NHIGGS
call check_LEP_bound1(Mh_avg(runid),
& 1.0D0,predratio(runid),nLEP_NSneutCHN+k+1,predicted_limit)
enddo
predicted_limit = .FALSE.
n=1
do k=0,nLEP_NSneutCHN-1,1
write(10,'(" (pred) id of table:",1I4)')k+1
runid=n+NHIGGS*k
call check_LEP_bound1(Mh_avg(runid),
& 1.0D0,predratio(runid),k+1,predicted_limit)
enddo
n=1
do k=0,nLEP_NScharCHN-1,1
write(10,'(" (pred) id of table:",1I4)')nLEP_NSneutCHN+k+1
runid=n+NHPLUS*k+nLEP_NSneutCHN*NHIGGS
call check_LEP_bound1(Mh_avg(runid),
& 1.0D0,predratio(runid),nLEP_NSneutCHN+k+1,predicted_limit)
enddo
close(10)
stop"checking tables only"
!...KW added
* Add the following to the end (just before last endif) of subroutine check_LEP_bound1 (note *not* check_LEP_bound2, which is the last subroutine in the file):
!KW added
do i=1,tabsize
write(10,"(2E17.10)"),Mhmin(Sid)+dble(i-1)*sep(Sid),S95_1d(pred_obs_flag,Sid,i)
enddo
write(*,*)"checking tables only!"
!...KW added
* In my code, uncomment this line in subroutine fill_pr in channels.F90:
call check_tables
* Run both codes:
in the KW HiggsBounds example_programs directory:
./demonstrateHiggsBounds_commandline.bat
in the OB HiggsBounds example_programs directory:
../HiggsBounds LandH effC 3 1 '../example_data/'
Note both programs will 'stop'. (i.e. the programs both end when they come across
'stop' statements in the code, rather than finishing naturally.)
* Then need to use a command that's something like:
diff /home/Karina/Work/HiggsBounds/HBsvn/HiggsBounds/trunk/HiggsBounds_KW/example_programs/debugging_tables_tevtype1.txt /home/Karina/Work/HiggsBounds/OB_code/email_02.04.10/HiggsBounds-f77_ed/example_programs/debugging_tables_tevtype1.txt > temp.txt

File Metadata

Mime Type
text/plain
Expires
Wed, May 14, 10:23 AM (1 d, 11 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5068031
Default Alt Text
notes_on_comparing_LEP_tables_type1.txt (2 KB)

Event Timeline