subroutine prinmain(nhist,nbins,iteration,title,outdir, | obgnorms,ncnorms,cohnorms,ccnorms,obg_count) implicit none ! Input/Output: integer nhist integer nbins(3) integer iteration character*55 title(22,3) !--> Plot Titles character*45 outdir(22) !--> Output directories ! Normalizatinos: ! 1st Index is nhist, 2nd is for iteration number ! Third index gives: ! 0 = Central Normalization ! 1 = +1 Sigma ! 2 = -1 Sigma ! 3 = MinChisq ! 4 = NbinsUsed double precision obgnorms(22,25) double precision ncnorms(22,0:25,0:4) double precision cohnorms(22,0:25,0:4) double precision ccnorms(22,0:25,0:4) double precision obg_count(22,15) !(nhist,ntyp) integer ihist !--> Do loop double precision SigmaCoh double precision SigmaNC double precision SigmaCC integer length integer lentit character*80 filename double precision OBGmin(22),OBGmax(22) !00000000000000000000000000000000000000000000000! ! Calclate OBG Upper and Lower Normalizations: do ihist=1,nhist OBGmax(ihist) = obgnorms(ihist,iteration)* | (1.0D0+1.0D0/dsqrt(obg_count(ihist,5))) OBGmin(ihist) = obgnorms(ihist,iteration)* | (1.0D0-1.0D0/dsqrt(obg_count(ihist,5))) enddo !00000000000000000000000000000000000000000000000! !----------------------------------------------------------------! ! Write normlist files: ! CohRho0: open(25,file='output/normlist_coh.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) cohnorms(ihist,iteration,0) enddo close(25) open(25,file='output/normlist_coh_-1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) cohnorms(ihist,iteration,2) enddo close(25) open(25,file='output/normlist_coh_+1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) cohnorms(ihist,iteration,1) enddo close(25) ! NCDIS: open(25,file='output/normlist_ncdis.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) ncnorms(ihist,iteration,0) enddo close(25) open(25,file='output/normlist_ncdis_-1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) ncnorms(ihist,iteration,2) enddo close(25) open(25,file='output/normlist_ncdis_+1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) ncnorms(ihist,iteration,1) enddo close(25) ! CCDIS: open(25,file='output/normlist_ccdis.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) ccnorms(ihist,iteration,0) enddo close(25) open(25,file='output/normlist_ccdis_-1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) ccnorms(ihist,iteration,2) enddo close(25) open(25,file='output/normlist_ccdis_+1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,101) ccnorms(ihist,iteration,1) enddo close(25) ! OBG: open(25,file='output/normlist_obg.txt',status='UNKNOWN') do ihist=1,nhist write(25,102) obgnorms(ihist,iteration) enddo close(25) open(25,file='output/normlist_obg_-1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,102) OBGmin(ihist) enddo close(25) open(25,file='output/normlist_obg_+1sigma.txt',status='UNKNOWN') do ihist=1,nhist write(25,102) OBGmax(ihist) enddo close(25) 101 format(F10.4) 102 format(F12.8) !----------------------------------------------------------------! do ihist=1,nhist !--> Loop over sets of histograms length=LEN_TRIM(outdir(ihist)) !Gets length of dirnames without trailing spaces !=================================================! !0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0! ! Calculate Sigmas: SigmaCoh = dsqrt(cohnorms(ihist,iteration,3)/ | cohnorms(ihist,iteration,4)) SigmaNC = dsqrt(ncnorms(ihist,iteration,3)/ | ncnorms(ihist,iteration,4)) SigmaCC = dsqrt(ccnorms(ihist,iteration,3)/ | ccnorms(ihist,iteration,4)) !0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0o0! !=================================================! 104 format(a37,F8.3,a6) 105 format(a26,F7.3,a5) 107 format(a16,F7.3,a8,F5.1,a9) !{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}! ! CohRho0 Latex File: lentit=LEN_TRIM(title(ihist,3)) filename = 'output/'//outdir(ihist)(1:length)//'/chisq_coh.tex' open(55,file=filename,status='UNKNOWN') write(55,*) '\\begin{table}[h!]\\centering' write(55,*) '{\\small{' write(55,*) '\\begin{tabular}{||l||r||r||}' write(55,*) '\\hline' write(55,*) '\\hline' write(55,104) '\\multicolumn{2}{||c||}{$\\chi^{2}$ Min: ', | cohnorms(ihist,iteration,3),'} & \\\\' write(55,*) '\\multicolumn{2}{||c||}{Number of bins used: ', | cohnorms(ihist,iteration,4),'} & \\\\' write(55,104) '\\multicolumn{2}{||c||}{One $\\sigma$: ', | SigmaCoh,'} & \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) ' & \\cohrz & \\\\' write(55,105) 'Norm at Min $\\chi^{2}$ & ', | cohnorms(ihist,iteration,0),' & \\\\' write(55,107) '$-1$ $\\sigma$ & ', | cohnorms(ihist,iteration,2),' & $($', |((cohnorms(ihist,iteration,0)-cohnorms(ihist,iteration,2))/ | cohnorms(ihist,iteration,0))*100.0,'$\\%)$ \\\\' write(55,107) '$+1$ $\\sigma$ & ', | cohnorms(ihist,iteration,1),' & $($', |((cohnorms(ihist,iteration,1)-cohnorms(ihist,iteration,0))/ | cohnorms(ihist,iteration,0))*100.0,'$\\%)$ \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) '\\end{tabular}' write(55,*) '\\caption{$\\chi^{2}$ for \\cohrz on plot:'// | title(ihist,3)(1:lentit)//'}' write(55,*) '\\label{tab-chicoh}' write(55,*) '}}' write(55,*) '\\end{table}' write(55,*) '\\endinput' close(55) !{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}! ! NCDIS Latex File: lentit=LEN_TRIM(title(ihist,1)) filename = 'output/'//outdir(ihist)(1:length)//'/chisq_nc.tex' open(55,file=filename,status='UNKNOWN') write(55,*) '\\begin{table}[h!]\\centering' write(55,*) '{\\small{' write(55,*) '\\begin{tabular}{||l||r||r||}' write(55,*) '\\hline' write(55,*) '\\hline' write(55,104) '\\multicolumn{2}{||c||}{$\\chi^{2}$ Min: ', | ncnorms(ihist,iteration,3),'} & \\\\' write(55,*) '\\multicolumn{2}{||c||}{Number of bins used: ', | ncnorms(ihist,iteration,4),'} & \\\\' write(55,104) '\\multicolumn{2}{||c||}{One $\\sigma$: ', | SigmaNC,'} & \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) ' & NCDIS & \\\\' write(55,105) 'Norm at Min $\\chi^{2}$ & ', | ncnorms(ihist,iteration,0),' & \\\\' write(55,107) '$-1$ $\\sigma$ & ', | ncnorms(ihist,iteration,2),' & $($', |((ncnorms(ihist,iteration,0)-ncnorms(ihist,iteration,2))/ | ncnorms(ihist,iteration,0))*100.0,'$\\%)$ \\\\' write(55,107) '$+1$ $\\sigma$ & ', | ncnorms(ihist,iteration,1),' & $($', |((ncnorms(ihist,iteration,1)-ncnorms(ihist,iteration,0))/ | ncnorms(ihist,iteration,0))*100.0,'$\\%)$ \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) '\\end{tabular}' write(55,*) '\\caption{$\\chi^{2}$ for NCDIS on plot:'// | title(ihist,1)(1:lentit)//'}' write(55,*) '\\label{tab-chinc}' write(55,*) '}}' write(55,*) '\\end{table}' write(55,*) '\\endinput' close(55) !{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}! ! CCDIS Latex File: lentit=LEN_TRIM(title(ihist,2)) filename = 'output/'//outdir(ihist)(1:length)//'/chisq_cc.tex' open(55,file=filename,status='UNKNOWN') write(55,*) '\\begin{table}[h!]\\centering' write(55,*) '{\\small{' write(55,*) '\\begin{tabular}{||l||r||r||}' write(55,*) '\\hline' write(55,*) '\\hline' write(55,104) '\\multicolumn{2}{||c||}{$\\chi^{2}$ Min: ', | ccnorms(ihist,iteration,3),'} & \\\\' write(55,*) '\\multicolumn{2}{||c||}{Number of bins used: ', | ccnorms(ihist,iteration,4),'} & \\\\' write(55,104) '\\multicolumn{2}{||c||}{One $\\sigma$: ', | SigmaCC,'} & \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) ' & CCDIS & \\\\' write(55,105) 'Norm at Min $\\chi^{2}$ & ', | ccnorms(ihist,iteration,0),' & \\\\' write(55,107) '$-1$ $\\sigma$ & ', | ccnorms(ihist,iteration,2),' & $($', |((ccnorms(ihist,iteration,0)-ccnorms(ihist,iteration,2))/ | ccnorms(ihist,iteration,0))*100.0,'$\\%)$ \\\\' write(55,107) '$+1$ $\\sigma$ & ', | ccnorms(ihist,iteration,1),' & $($', |((ccnorms(ihist,iteration,1)-ccnorms(ihist,iteration,0))/ | ccnorms(ihist,iteration,0))*100.0,'$\\%)$ \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) '\\end{tabular}' write(55,*) '\\caption{$\\chi^{2}$ for CCDIS on plot:'// | title(ihist,2)(1:lentit)//'}' write(55,*) '\\label{tab-chicc}' write(55,*) '}}' write(55,*) '\\end{table}' write(55,*) '\\endinput' close(55) !{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}! ! OBG Latex File: 115 format(a26,F10.6,a5) 117 format(a16,F10.6,a8,F5.1,a9) filename = 'output/'//outdir(ihist)(1:length)//'/chisq_obg.tex' open(55,file=filename,status='UNKNOWN') write(55,*) '\\begin{table}[h!]\\centering' write(55,*) '{\\small{' write(55,*) '\\begin{tabular}{||l||r||r||}' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) ' & OBG & \\\\' write(55,115) 'Norm at Min $\\chi^{2}$ & ', | obgnorms(ihist,iteration),' & \\\\' write(55,117) '$-1$ $\\sigma$ & ', | OBGmin(ihist),' & $($', | ((obgnorms(ihist,iteration)-OBGmin(ihist))/ | obgnorms(ihist,iteration))*100.0,'$\\%)$ \\\\' write(55,117) '$+1$ $\\sigma$ & ', | OBGmax(ihist),' & $($', | ((OBGmax(ihist)-obgnorms(ihist,iteration))/ | obgnorms(ihist,iteration))*100.0,'$\\%)$ \\\\' write(55,*) '\\hline' write(55,*) '\\multicolumn{2}{||c||}{Number of OBG used', | obg_count(ihist,5),'} & \\\\' write(55,*) '\\hline' write(55,*) '\\hline' write(55,*) '\\end{tabular}' write(55,*) '\\caption{Normalization for OBG by number of events | within Kaon mass range}' write(55,*) '\\label{tab-normobg}' write(55,*) '}}' write(55,*) '\\end{table}' write(55,*) '\\endinput' close(55) !{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}{}! enddo !--> Loop over sets of histograms return end