Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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 {
39  if ( binSchemeName != "linear" ) {
40  if ( binSchemeName == "log" )
41  binScheme = G4BinScheme::kLog;
42  else {
43  // There is no name associated with G4BinScheme::kUser
44  G4ExceptionDescription description;
45  description
46  << " \"" << binSchemeName << "\" 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  G4double unit, G4Fcn fcn, G4BinScheme binScheme,
58  std::vector<G4double>& edges)
59 {
60 // Compute edges from parameters
61 
62  // Apply units
63  auto xumin = xmin/unit;
64  auto xumax = xmax/unit;
65 
66  if ( binScheme == G4BinScheme::kLinear ) {
67  auto dx = (fcn(xumax) - fcn(xumin) ) / nbins;
68  auto binValue = fcn(xumin);
69  while ( G4int(edges.size()) <= nbins ) { // Loop checking, 23.06.2015, I. Hrivnacova
70  edges.push_back(binValue);
71  binValue += dx;
72  }
73  }
74  else if ( binScheme == G4BinScheme::kLog ) {
75  // do not apply fcn
76  auto dlog
77  = (std::log10(xumax) - std::log10(xumin))/ nbins;
78  auto dx = std::pow(10, dlog);
79  auto binValue = xumin;
80  while ( G4int(edges.size()) <= nbins ) { // Loop checking, 23.06.2015, I. Hrivnacova
81  edges.push_back(binValue);
82  binValue *= dx;
83  }
84  }
85  else if ( binScheme == G4BinScheme::kUser ) {
86  // This should never happen, but let's make sure about it
87  // by issuing a warning
88  G4ExceptionDescription description;
89  description
90  << " User binning scheme setting was ignored." << G4endl
91  << " Linear binning will be applied with given (nbins, xmin, xmax) values";
92  G4Exception("G4Analysis::ComputeEdges",
93  "Analysis_W013", JustWarning, description);
94  }
95 }
96 
97 //_____________________________________________________________________________
98 void ComputeEdges(const std::vector<G4double>& edges,
99  G4double unit, G4Fcn fcn,
100  std::vector<G4double>& newBins)
101 {
102 // Apply function to defined edges
103 
104  for (auto element : edges) {
105  newBins.push_back(fcn(element/unit));
106  }
107 }
108 
109 }
G4double(*)(G4double) G4Fcn
Definition: G4Fcn.hh:36
std::ostringstream G4ExceptionDescription
Definition: globals.hh:76
void ComputeEdges(G4int nbins, G4double xmin, G4double xmax, G4double unit, G4Fcn fcn, G4BinScheme, std::vector< G4double > &edges)
Definition: G4BinScheme.cc:56
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
G4BinScheme
Definition: G4BinScheme.hh:40
double G4double
Definition: G4Types.hh:76