Geant4_10
G4BinScheme.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id:$
27 
28 // Author: Ivana Hrivnacova, 22/08/2013 (ivana@ipno.in2p3.fr)
29 
30 #include "G4BinScheme.hh"
31 
32 namespace G4Analysis
33 {
34 
35 //_____________________________________________________________________________
36 G4BinScheme GetBinScheme(const G4String& binSchemeName)
37 {
38  G4BinScheme binScheme = kLinearBinScheme;
39  if ( binSchemeName != "linear" ) {
40  if ( binSchemeName == "log" )
41  binScheme = kLogBinScheme;
42  else {
43  // There is no name associated with kUserBinScheme
44  G4ExceptionDescription description;
45  description
46  << " \"" << binScheme << "\" binning scheme is not supported." << G4endl
47  << " " << "Linear binning will be applied.";
48  G4Exception("G4Analysis::GetBinScheme",
49  "Analysis_W013", JustWarning, description);
50  }
51  }
52  return binScheme;
53 }
54 
55 //_____________________________________________________________________________
56 void ComputeEdges(G4int nbins, G4double xmin, G4double xmax,
57  G4Fcn fcn, G4BinScheme binScheme,
58  std::vector<G4double>& edges)
59 {
60 // Compute edges from parameters
61 
62  if ( binScheme == kLinearBinScheme ) {
63  G4double dx = (xmax - xmin ) / nbins;
64  G4double binValue = xmin;
65  while ( G4int(edges.size()) <= nbins ) {
66  edges.push_back(binValue);
67  binValue += dx;
68  }
69  }
70  else if ( binScheme == kLogBinScheme ) {
71  G4double dlog
72  = (std::log10(fcn(xmax)) - std::log10(fcn(xmin)))/ nbins;
73  G4double dx = std::pow(10, dlog);
74  G4double binValue = fcn(xmin);
75  while ( G4int(edges.size()) <= nbins ) {
76  edges.push_back(binValue);
77  binValue *= dx;
78  }
79  }
80  else if ( binScheme == kUserBinScheme ) {
81  // This should never happen, but let's make sure about it
82  // by issuing a warning
83  G4ExceptionDescription description;
84  description
85  << " User binning scheme setting was ignored." << G4endl
86  << " Linear binning will be applied with given (nbins, xmin, xmax) values";
87  G4Exception("G4Analysis::ComputeEdges",
88  "Analysis_W013", JustWarning, description);
89  }
90 }
91 
92 //_____________________________________________________________________________
93 void ComputeEdges(const std::vector<G4double>& edges,
94  G4Fcn fcn,
95  std::vector<G4double>& newBins)
96 {
97 // Apply function to defined edges
98  std::vector<G4double>::const_iterator it;
99  for (it = edges.begin(); it != edges.end(); it++ ) {
100  newBins.push_back(fcn(*it));
101  }
102 }
103 
104 }
std::ostringstream G4ExceptionDescription
Definition: globals.hh:76
G4double(* G4Fcn)(G4double)
Definition: G4Fcn.hh:36
int G4int
Definition: G4Types.hh:78
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
G4BinScheme GetBinScheme(const G4String &binSchemeName)
Definition: G4BinScheme.cc:36
#define G4endl
Definition: G4ios.hh:61
void ComputeEdges(G4int nbins, G4double xmin, G4double xmax, G4Fcn fcn, G4BinScheme, std::vector< G4double > &edges)
Definition: G4BinScheme.cc:56
G4BinScheme
Definition: G4BinScheme.hh:40
double G4double
Definition: G4Types.hh:76