TABLE OF CONTENTS
- 1. /hist2Df90
- 1.1. hist2Df90/histogram2D
- 1.2. hist2Df90/AddHist2D
- 1.3. hist2Df90/CreateHist2D
- 1.4. hist2Df90/WriteHist2D_SYSTEM
- 1.5. hist2Df90/RemoveHist2D
- 1.6. hist2Df90/CopyHist2D
- 1.7. hist2Df90/WriteHist2D_Gnuplot
- 1.8. hist2Df90/WriteHist2D_Gnuplot_Bspline
- 1.9. hist2Df90/DivideHist2D
- 1.10. hist2Df90/ReadHist2D_Gnuplot
- 1.11. hist2Df90/DumpHist2D
- 1.12. hist2Df90/FetchHist2D
- 1.13. hist2Df90/IntegrateHist2D
- 1.14. hist2Df90/AverageHist2D
/hist2Df90 [ Modules ]
NAME
module hist2Df90
PURPOSE
Encapsulate all routines and datas for 2D Histograms.
Features of Histograms provided by this module:
- store paramaters of the (x1,x2)-binning
- enable two y-values (y and y2)
- track keeping of under-/over-score the given extreme values of x.
- provide simple-to-understand output routines
- provide simple histogram arithmetic
- many others...
NOTES
Programming is fairly similar to "module histf90"
INPUTS
...(This module needs no input)
hist2Df90/histogram2D [ Types ]
[ Top ] [ hist2Df90 ] [ Types ]
NAME
type histogram2D
PURPOSE
Type definition to store all information for a 2D Histogram.
SOURCE
type histogram2D real :: xMin(2) ! smallest x-value: (x1,x2) real :: xMax(2) ! largest x-value: (x1,x2) real :: xBin(2) ! width of x-Bin: (x1,x2) real :: xExtreme(2,2) ! extremes of used x-values ! (x1,x2), (min,max) character*(NameLength) :: Name ! name to be written real,allocatable :: yVal(:,:,:) ! histogram values: ! (x1),(x2),(yy1,yy2,yy3) end type histogram2D
hist2Df90/AddHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine AddHist2D(H, x,y,y2)
subroutine AddHist2D(H1,H2, x,y,y2)
PURPOSE
Add to the given histogram(s) at the given (x1,x2)-value the weight y (y2). If y2 is given, also the second column is filled.
INPUTS
- type(histogram2D) :: H -- Histogram to be used
or:
- type(histogram2D) :: H1,H2 -- Histograms to be used
- real,dimension(2) :: x -- (x1,x2)-value
- real :: y -- weight to be added
- real :: y2 -- second weight to be added [OPTIONAL]
OUTPUT
H is changed
NOTES
This routine is overloaded: Giving two histograms at the same time is like calling the routine twice.
hist2Df90/CreateHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine CreateHist2D(H, HName,x1,x2,bin,verbose)
PURPOSE
This is the Constructor of a 2D-Histogram! Allocate Memory for the entries and put additional variables to their default.
INPUTS
- type(histogram2D) :: H -- Histogram to be created
- character*(*) :: HName -- Name of Histogram
- real,dimension(2) :: x1,x2,bin -- Minimal/maximal values for x-coordinates to be considered, bin-width
- logical :: verbose -- switch for verbosity [OPTIONAL]
OUTPUT
H is changed
hist2Df90/WriteHist2D_SYSTEM [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine WriteHist2D_SYSTEM(H, iFile)
PURPOSE
Write out System information of the histogram
INPUTS
- type(histogram2D) :: H -- Histogram to be used
- integer :: iFile -- File number output to redirect
OUTPUT
written to file iFile
hist2Df90/RemoveHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine RemoveHist2D(H)
PURPOSE
Free the allocated memory
INPUTS
- type(histogram2D) :: H -- Histogram to be used
OUTPUT
H is changed
hist2Df90/CopyHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine CopyHist2D(H1,H2)
PURPOSE
Copies Histogram H1 into H2 (H2 is overwritten)
INPUTS
- type(histogram2D) :: H1,H2 -- Histograms to be used
OUTPUT
H2 is changed
hist2Df90/WriteHist2D_Gnuplot [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine WriteHist2D_Gnuplot(H,iFile,add,mul, iColumn,DoAve,MaxVal,H2,file,dump,SwapXY)
PURPOSE
Write out the 2D-histogram. Format is suitable as input to gnuplots "splot" command
The entries are multiplied by 'mul' and 'add' is added.
If iColumn is not given, all 3 columns are writen.
INPUTS
- type(histogram2D) :: H -- Histogram to be used
- integer :: iFile -- File number output to redirect [OPTIONAL]
- real :: add -- factor to add [OPTIONAL]
- real :: mul -- factor to multiply [OPTIONAL]
- integer :: iColumn -- Column to write [OPTIONAL]
- logical :: DoAve -- write also "average" (cf Notes) [OPTIONAL]
- real :: maxVal -- value used instead "divide by zero" [OPTIONAL]
- type(histogram) :: H2 -- Histogram to divide by [OPTIONAL]
- character*(*) :: file -- name of file to open and close [OPTIONAL]
- logical :: dump -- if true, also dump it binary to file [OPTIONAL]
- logical :: SwapXY -- write data sorted according y:x:z instead of x:y:z [OPTIONAL]
OUTPUT
write to file number
NOTES
The Histogram Data is not affected!!!
if DoAve is given: Write out the 2D-histogram, but now also divide column 3 by column 1. If column 2 is zero, i.e. no entries, then parameter maxVal is written instead. This column output is NOT divided by bin widths.
hist2Df90/WriteHist2D_Gnuplot_Bspline [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine WriteHist2D_Gnuplot_Bspline(H,iFile,add,mul, iColumn,nPoints,DoAve,MaxVal)
PURPOSE
Write out the 2D-histogram. Format is suitable as input to gnuplots "splot" command
The entries are multiplied by 'mul' and 'add' is added.
If iColumn is not given, all 3 columns are writen.
INPUTS
- type(histogram2D) :: H -- Histogram to be used
- integer :: iFile -- File number output to redirect
- real :: add -- factor to add [OPTIONAL]
- real :: mul -- factor to multiply [OPTIONAL]
- integer :: iColumn -- Column to write [OPTIONAL]
- integer :: nPoints -- number of sub-points [OPTIONAL]
- logical :: DoAve -- write also "average" (cf Notes) [OPTIONAL]
- real :: maxVal -- value used instead "divide by zero" [OPTIONAL]
OUTPUT
write to file number
NOTES
The Histogram Data is not affected!!!
hist2Df90/DivideHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine DivideHist2D(H1,H2, maxVal, MulWithBin)
PURPOSE
calculate H1 = H1/H2 on a bin-by-bin basis
INPUTS
- type(histogram2D) :: H1,H2 -- Histograms to be used
- real :: maxVal -- value used instead "divide by zero" [OPTIONAL]
- logical :: MulWithBin -- flag: multiply with bin-width [OPTIONAL]
Without explicit "MulWithBin" parameter, the default behaviour is as with "MulWithBin=.TRUE."
OUTPUT
H1 is changed
NOTES
This Dividing-Routine ignores all EXTREME values etc. !
The "MulWithBin=.TRUE." feature is provided, because during output the bin entry values are divided by the bin-width (as usual). Therefore if you divide a histogram by itself and writ it out, you would not get the bin-Values "1" but "1/bin-width". This flag here corrects for this.
hist2Df90/ReadHist2D_Gnuplot [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine ReadHist2D_Gnuplot(H,iFile,add,mul, iColumn, DoAve)
PURPOSE
Read in a 2D-histogram from a file, which has been written in a gnuplot format.
hist2Df90/DumpHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine DumpHist2D(H,file,iFile)
PURPOSE
Write all the histogram information unformatted (i.e. binary) to a file
INPUTS
- type(histogram2D):: H -- Histogram to be used
- character*(*) :: file -- name of file to open and close
- integer,OPTIONAL :: iFile -- File number output to redirect [OPTIONAL]
- real :: add -- factor to add [OPTIONAL]
- real :: mul -- factor to multiply [OPTIONAL]
OUTPUT
H is written UNFORMATTED to the given file
hist2Df90/FetchHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine FetchHist2D(H,file,iFile, add,mul, flagOK)
PURPOSE
Read in all the histogram information previously dumped unformatted (i.e. binary) to a file
INPUTS
- character*(*) :: file -- name of file to open and close
- integer,OPTIONAL:: iFile -- File number input to redirect [OPTIONAL]
OUTPUT
- type(histogram2D) :: H -- Histogram to be used
- logical :: flagOK -- flag, if reading was okay [OPTIONAL]
H is read UNFORMATTED from the given file. Sizes are calculated as in CreateHist, also memory is allocated.
NOTES
No checks about input are performed!
hist2Df90/IntegrateHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine IntegrateHist2D(H2D,H,axis)
PURPOSE
Create a 1D histogram, which represents the integration along one axis of the 2D histogram
INPUTS
- type(histogram2D) :: H2D -- 2D histogram to integrate
- integer :: axis -- axis to preserve (integration is done for the other axis)
OUTPUT
hist2Df90/AverageHist2D [ Subroutines ]
[ Top ] [ hist2Df90 ] [ Subroutines ]
NAME
subroutine AverageHist2D(H2D,H,axis)
PURPOSE
Create a 1D histogram, which represents the average along one axis of the 2D histogram
INPUTS
- type(histogram2D) :: H2D -- 2D histogram to integrate
- integer :: axis -- axis to preserve (integration is done for the other axis)
OUTPUT