### 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

typehistogram2Dreal :: 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 typehistogram2D

## 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