Geant4 Cross Reference |
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 27 // Author: Ivana Hrivnacova, 10/08/2022 (ivana@ipno.in2p3.fr) 28 29 #include "G4THnToolsManager.hh" 30 31 #include "tools/histo/p1d" 32 33 using namespace G4Analysis; 34 35 #include <fstream> 36 37 // Specialization for P1 type 38 39 40 //_____________________________________________________________________________ 41 template <> 42 tools::histo::p1d* G4THnToolsManager<kDim2, tools::histo::p1d>::CreateToolsHT( 43 const G4String& title, 44 const std::array<G4HnDimension, kDim2>& bins, 45 const std::array<G4HnDimensionInformation, kDim2>& hnInfo) 46 { 47 // Apply hn information to bins 48 auto newXBins(bins[kX]); 49 Update(newXBins, hnInfo[kX]); 50 auto newYBins(bins[kY]); 51 UpdateValues(newYBins, hnInfo[kY]); 52 53 if (hnInfo[kX].fBinScheme == G4BinScheme::kLinear) { 54 if ( newYBins.fMinValue == 0. && newYBins.fMaxValue == 0.) { 55 return new tools::histo::p1d( 56 title, newXBins.fNBins, newXBins.fMinValue, newXBins.fMaxValue); 57 } 58 return new tools::histo::p1d( 59 title, newXBins.fNBins, newXBins.fMinValue, newXBins.fMaxValue, 60 newYBins.fMinValue, newYBins.fMaxValue); 61 } 62 63 if ( newYBins.fMinValue == 0. && newYBins.fMaxValue == 0.) { 64 return new tools::histo::p1d(title, newXBins.fEdges); 65 } 66 67 return new tools::histo::p1d(title, newXBins.fEdges, 68 newYBins.fMinValue, newYBins.fMaxValue); 69 } 70 71 //_____________________________________________________________________________ 72 template <> 73 void G4THnToolsManager<kDim2, tools::histo::p1d>::ConfigureToolsHT( 74 tools::histo::p1d* ht, 75 const std::array<G4HnDimension, kDim2>& bins, 76 const std::array<G4HnDimensionInformation, kDim2>& hnInfo) 77 { 78 // Apply hn information to bins 79 auto newXBins(bins[kX]); 80 Update(newXBins, hnInfo[kX]); 81 auto newYBins(bins[kY]); 82 UpdateValues(newYBins, hnInfo[kY]); 83 84 if (hnInfo[kX].fBinScheme == G4BinScheme::kLinear) { 85 if ( newYBins.fMinValue == 0. && newYBins.fMaxValue == 0.) { 86 ht->configure( 87 newXBins.fNBins, newXBins.fMinValue, newXBins.fMaxValue); 88 return; 89 } 90 ht->configure( 91 newXBins.fNBins, newXBins.fMinValue, newXBins.fMaxValue, 92 newYBins.fMinValue, newYBins.fMaxValue); 93 return; 94 } 95 96 if ( newYBins.fMinValue == 0. && newYBins.fMaxValue == 0.) { 97 ht->configure(newXBins.fEdges); 98 return; 99 } 100 ht->configure(newXBins.fEdges, newYBins.fMinValue, newYBins.fMaxValue); 101 return; 102 } 103 104 //_____________________________________________________________________________ 105 template <> 106 G4bool G4THnToolsManager<kDim2, tools::histo::p1d>::FillHT( 107 tools::histo::p1d* ht, const G4HnInformation& hnInformation, 108 std::array<G4double, kDim2>& value, G4double weight) 109 { 110 const auto& xInfo = hnInformation.GetHnDimensionInformation(kX); 111 const auto& yInfo = hnInformation.GetHnDimensionInformation(kY); 112 113 // Apply hn information to value 114 Update(value[kX], xInfo); 115 Update(value[kY], yInfo); 116 117 // Fill updated value 118 ht->fill(value[kX], value[kY], weight); 119 120 return true; 121 } 122 123 //_____________________________________________________________________________ 124 template <> 125 G4bool G4THnToolsManager<kDim2, tools::histo::p1d>::WriteOnAscii( 126 std::ofstream& output) 127 { 128 // Write selected objects on ASCII file 129 130 // Do nothing if no histograms are selected 131 if ( ! GetHnManager()->IsAscii() ) return true; 132 133 // Write p1 histograms 134 auto id = GetHnManager()->GetFirstId(); 135 for (const auto& [p1, info] : *GetTHnVector()) { 136 137 if ( (p1 == nullptr) || (! info->GetAscii()) ) { 138 // skip writing 139 // if p1 was deleted or writing ascii is not selected 140 id++; 141 continue; 142 } 143 144 Message(kVL3, "write on ascii", "p1d", info->GetName()); 145 146 output << "\n 1D profile " << id++ << ": " << p1->title() 147 << "\n \n \t \t X \t\t MeanY" << G4endl; 148 149 for (G4int j=0; j< G4int(p1->axis().bins()); ++j) { 150 auto sw = p1->bin_Sw(j); 151 auto svw = p1->bin_Svw(j); 152 auto mean = ( sw != 0. ) ? (svw / sw) : 0.; 153 output << " " << j << "\t" 154 << p1->axis().bin_center(j) << "\t" 155 << mean << G4endl; 156 } 157 } 158 159 return output.good(); 160 } 161