Functions
gnumi2simple.C File Reference
#include "Numerical/RandomGen.h"
#include "FluxDrivers/GNuMIFlux.h"
#include "FluxDrivers/GSimpleNtpFlux.h"
#include "Utils/UnitUtils.h"
#include "TSystem.h"
#include "TStopwatch.h"
#include "TLorentzVector.h"
#include "TNtuple.h"
#include "TFile.h"
#include <iostream>
#include <iomanip>
#include <string>
#include <set>

Go to the source code of this file.

Functions

void fill_simple (string fname="fluxntgenie.root", string cfg="MINOS-Near", int flxset=0, long int nentries=0, double pots=1.0e30, double enumin=0.0, bool doaux=false)
 
void gnumi2simple (string fname, string det, string flux, long int nentries=5000, double pots=1.0e30, double enumin=0.0, bool doaux=false)
 
void make_simple_grid (string fname, string det, string flux, long int nentries=5000, double pots=1.0e30, double enumin=0.0, bool doaux=false)
 

Function Documentation

void fill_simple ( string  fname = "fluxntgenie.root",
string  cfg = "MINOS-Near",
int  flxset = 0,
long int  nentries = 0,
double  pots = 1.0e30,
double  enumin = 0.0,
bool  doaux = false 
)

Definition at line 65 of file gnumi2simple.C.

67 {
68  string fluxpatt = "$FLUXPATH";
69 
70  int flxver = flxset % 10;
71  switch ( flxver ) {
72  case 0:
73  fluxpatt += "/v19/fluka05_le010z185i/root/fluka05_le010z185i_";
74  break;
75  case 1:
76  fluxpatt += "/v19/fluka05_le010z185i_lowcut/root/fluka05_le010z185i_";
77  break;
78  case 2:
79  fluxpatt += "/flugg/flugg_le010z185i_run1/root/flugg_le010z185i_run1_";
80  break;
81  case 3:
82  fluxpatt += "/flugg_lowth/flugg_le010z185i_run1_lowth/root/flugg_le010z185i_run1_lowth_";
83  break;
84  default:
85  cout << " not a legal flux set " << endl;
86  return;
87  }
88 
89  if ( flxset/10 == 1 ) {
90  if ( flxver == 2 || flxver == 3 ) fluxpatt += "001";
91  else fluxpatt += "1";
92  } else if ( flxset/10 == 2 ) {
93  if ( flxver == 2 || flxver == 3 ) fluxpatt += "00[12]";
94  else fluxpatt += "[12]";
95  } else fluxpatt += "*";
96  fluxpatt += ".root";
97 
98  string fluxfname(gSystem->ExpandPathName(fluxpatt.c_str()));
99  flux::GNuMIFlux* gnumi = new GNuMIFlux();
100  gnumi->LoadBeamSimData(fluxfname,cfg);
101  gnumi->SetEntryReuse(1); // don't reuse entries when reformatting
102 
103  //gnumi->PrintConfig();
104  //cout << " change to cm:" << endl;
105  //double cm_unit = genie::utils::units::UnitFromString("cm");
106  //gnumi->SetLengthUnits(cm_unit);
107  //gnumi->PrintConfig();
108  //cout << " change to m:" << endl;
109  ///////RWH: the following 2 lines should work ...
110  //double m_unit = genie::utils::units::UnitFromString("m");
111  //gnumi->SetLengthUnits(m_unit);
112  gnumi->PrintConfig();
113  //double scaleusr = 100.;
114 
115  gnumi->SetEntryReuse(1); // reuse entries
116  gnumi->SetUpstreamZ(-3e38); // leave ray on flux window
117 
118  GFluxI* fdriver = dynamic_cast<GFluxI*>(gnumi);
119 
120  if (nentries == 0) nentries = 100000*20;
121 
122  // so as not to include scan time generate 1 nu
123  fdriver->GenerateNext();
124 
125  TFile* file = TFile::Open(fname.c_str(),"RECREATE");
126  TTree* fluxntp = new TTree("flux","a simple flux n-tuple");
127  TTree* metantp = new TTree("meta","metadata for flux n-tuple");
128  genie::flux::GSimpleNtpEntry* fentry = new genie::flux::GSimpleNtpEntry;
129  genie::flux::GSimpleNtpNuMI* fnumi = new genie::flux::GSimpleNtpNuMI;
130  genie::flux::GSimpleNtpAux* faux = new genie::flux::GSimpleNtpAux;
131  genie::flux::GSimpleNtpMeta* fmeta = new genie::flux::GSimpleNtpMeta;
132  fluxntp->Branch("entry",&fentry);
133  fluxntp->Branch("numi",&fnumi);
134  if (doaux) fluxntp->Branch("aux",&faux);
135  metantp->Branch("meta",&fmeta);
136 
137  TLorentzVector p4u, x4u;
138  //TLorentzVector p4b, x4b;
139  long int ngen = 0;
140 
141  set<int> pdglist;
142  double maxe = 0;
143  double minwgt = +1.0e10;
144  double maxwgt = -1.0e10;
145 
146  UInt_t metakey = TString::Hash(fname.c_str(),strlen(fname.c_str()));
147  cout << "metakey " << metakey << endl;
148  // ensure that we don't get smashed by UInt_t vs Int_t
149  metakey &= 0x7FFFFFFF;
150  cout << "metakey " << metakey << " after 0x7FFFFFFF" << endl;
151  cout << "=========================== Start " << endl;
152 
153  TStopwatch sw;
154  sw.Start();
155  int nok = 0, nwrite = 0;
156  while ( nok < nentries && gnumi->UsedPOTs() < pots ) {
157  //for (int i=0; i<nentries; ++i) {
158  fdriver->GenerateNext();
159  ngen++;
160  fentry->Reset();
161  fnumi->Reset();
162  faux->Reset();
163 
164  fentry->metakey = metakey;
165  fentry->pdg = fdriver->PdgCode();
166  fentry->wgt = gnumi->Weight();
167  x4u = fdriver->Position();
168  fentry->vtxx = x4u.X();
169  fentry->vtxy = x4u.Y();
170  fentry->vtxz = x4u.Z();
171  fentry->dist = gnumi->GetDecayDist();
172  p4u = fdriver->Momentum();
173  fentry->px = p4u.Px();
174  fentry->py = p4u.Py();
175  fentry->pz = p4u.Pz();
176  fentry->E = p4u.E();
177 
178  fnumi->run = gnumi->PassThroughInfo().run;
179  fnumi->evtno = gnumi->PassThroughInfo().evtno;
180  fnumi->entryno = gnumi->GetEntryNumber();
181 
182  fnumi->tpx = gnumi->PassThroughInfo().tpx;
183  fnumi->tpy = gnumi->PassThroughInfo().tpy;
184  fnumi->tpz = gnumi->PassThroughInfo().tpz;
185  fnumi->vx = gnumi->PassThroughInfo().vx;
186  fnumi->vy = gnumi->PassThroughInfo().vy;
187  fnumi->vz = gnumi->PassThroughInfo().vz;
188 
189  fnumi->ndecay = gnumi->PassThroughInfo().ndecay;
190  fnumi->ppmedium = gnumi->PassThroughInfo().ppmedium;
191  fnumi->tptype = gnumi->PassThroughInfo().tptype;
192 
193 
194  if ( doaux ) {
195  faux->auxint.push_back(gnumi->PassThroughInfo().tptype);
196  faux->auxdbl.push_back(gnumi->PassThroughInfo().tpx);
197  faux->auxdbl.push_back(gnumi->PassThroughInfo().tpy);
198  faux->auxdbl.push_back(gnumi->PassThroughInfo().tpz);
199  }
200 
201  if ( fentry->E > enumin ) ++nok;
202  fluxntp->Fill();
203  ++nwrite;
204 
205  // accumulate meta data
206  pdglist.insert(fentry->pdg);
207  minwgt = TMath::Min(minwgt,fentry->wgt);
208  maxwgt = TMath::Max(maxwgt,fentry->wgt);
209  maxe = TMath::Max(maxe,fentry->E);
210 
211  }
212  cout << "=========================== Complete " << endl;
213 
214  //fmeta->nflavors = pdglist.size();
215  //if ( fmeta->nflavors > 0 ) {
216  // fmeta->flavor = new Int_t[fmeta->nflavors];
217  // set<int>::const_iterator itr = pdglist.begin();
218  // int indx = 0;
219  // for ( ; itr != pdglist.end(); ++itr, ++indx ) fmeta->flavor[indx] = *itr;
220  //}
221  fmeta->pdglist.clear();
222  set<int>::const_iterator setitr = pdglist.begin();
223  for ( ; setitr != pdglist.end(); ++setitr) fmeta->pdglist.push_back(*setitr);
224 
225  fmeta->maxEnergy = maxe;
226  fmeta->minWgt = minwgt;
227  fmeta->maxWgt = maxwgt;
228  fmeta->protons = gnumi->UsedPOTs();
229  TVector3 p0, p1, p2;
230  gnumi->GetFluxWindow(p0,p1,p2);
231  TVector3 d1 = p1 - p0;
232  TVector3 d2 = p2 - p0;
233  fmeta->windowBase[0] = p0.X();
234  fmeta->windowBase[1] = p0.Y();
235  fmeta->windowBase[2] = p0.Z();
236  fmeta->windowDir1[0] = d1.X();
237  fmeta->windowDir1[1] = d1.Y();
238  fmeta->windowDir1[2] = d1.Z();
239  fmeta->windowDir2[0] = d2.X();
240  fmeta->windowDir2[1] = d2.Y();
241  fmeta->windowDir2[2] = d2.Z();
242  if ( doaux ) {
243  fmeta->auxintname.push_back("tptype");
244  fmeta->auxdblname.push_back("tpx");
245  fmeta->auxdblname.push_back("tpy");
246  fmeta->auxdblname.push_back("tpz");
247  }
248  string fname = "flux_pattern=" + fluxfname;
249  fmeta->infiles.push_back(fluxfname); // should get this expanded from gnumi
250  vector<string> flist = gnumi->GetFileList();
251  for (size_t i = 0; i < flist.size(); ++i) {
252  fname = "flux_infile=" + flist[i];
253  fmeta->infiles.push_back(fname);
254  }
255 
256  fmeta->seed = RandomGen::Instance()->GetSeed();
257  fmeta->metakey = metakey;
258 
259  metantp->Fill();
260 
261  sw.Stop();
262  cout << "Generated " << nwrite << " (" << nok << " w/ E >" << enumin << " ) "
263  << " ( request " << nentries << " ) "
264  << endl
265  << gnumi->UsedPOTs() << " POTs " << " ( request " << pots << " )"
266  << endl
267  << " pulled NFluxNeutrinos " << gnumi->NFluxNeutrinos()
268  << endl
269  << "Time to generate: " << endl;
270  sw.Print();
271 
272  cout << "===================================================" << endl;
273 
274  cout << "Last GSimpleNtpEntry: " << endl
275  << *fentry
276  << endl
277  << "Last GSimpleNtpMeta: " << endl
278  << *fmeta
279  << endl;
280 
281  file->cd();
282 
283  // write meta as simple object at top of file
284  //fmeta->Write();
285 
286  // write ntuples out
287  fluxntp->Write();
288  metantp->Write();
289  file->Close();
290 
291  cout << endl << endl;
292  gnumi->PrintConfig();
293 }
genie::flux::GSimpleNtpEntry::E
Double_t E
energy in lab frame
Definition: GSimpleNtpFlux.h:80
genie::flux::GNuMIFlux::GetFluxWindow
void GetFluxWindow(TVector3 &p1, TVector3 &p2, TVector3 &p3) const
3 points define a plane in beam coordinate
Definition: GNuMIFlux.cxx:1060
genie::flux::GNuMIFlux::SetEntryReuse
void SetEntryReuse(long int nuse=1)
of times to use entry before moving to next
Definition: GNuMIFlux.cxx:924
genie::flux::GSimpleNtpEntry::px
Double_t px
x momentum in lab frame
Definition: GSimpleNtpFlux.h:77
i
size_t i(0)
genie::flux::GNuMIFlux::GetEntryNumber
Long64_t GetEntryNumber()
index in chain
Definition: GNuMIFlux.h:253
genie::flux::GSimpleNtpEntry::vtxy
Double_t vtxy
y position in lab frame
Definition: GSimpleNtpFlux.h:73
genie::flux::GSimpleNtpNuMI::vx
Double_t vx
vertex position of hadron/muon decay
Definition: GSimpleNtpFlux.h:106
genie::flux::GSimpleNtpMeta::maxEnergy
Double_t maxEnergy
maximum energy
Definition: GSimpleNtpFlux.h:173
genie::GFluxI::Position
virtual const TLorentzVector & Position(void)=0
returns the flux neutrino 4-position (note: expect SI rather than physical units) ...
genie::flux::GSimpleNtpMeta::auxintname
std::vector< std::string > auxintname
tagname of aux ints associated w/ entry
Definition: GSimpleNtpFlux.h:182
genie::flux::GSimpleNtpNuMI::tptype
Int_t tptype
parent particle type at target exit
Definition: GSimpleNtpFlux.h:119
genie::flux::GSimpleNtpAux::auxdbl
std::vector< Double_t > auxdbl
additional doubles associated w/ entry
Definition: GSimpleNtpFlux.h:144
dbjson.cfg
cfg
Definition: dbjson.py:29
genie::flux::GSimpleNtpMeta::windowDir1
Double_t windowDir1[3]
dx,dy,dz of window direction 1
Definition: GSimpleNtpFlux.h:179
genie::flux::GSimpleNtpAux::auxint
std::vector< Int_t > auxint
additional ints associated w/ entry
Definition: GSimpleNtpFlux.h:143
cout
the ParameterSet object passed in for the configuration of a destination should be the only source that can affect the behavior of that destination This is to eliminate the dependencies of configuring a destination from multiple mostly from the defaults It suppresses possible glitches about changing the configuration file somewhere outside of a destination segament might still affect the behavior of that destination In the previous configuration for a specific the value of a certain e may come from following and have been suppressed It the configuring ParameterSet object for each destination will be required to carry a parameter list as complete as possible If a parameter still cannot be found in the ParameterSet the configuration code will go look for a hardwired default directly The model is a great simplicity comparing with the previous especially when looking for default values Another great advantage is most of the parameters now have very limited places that allows to appear Usually they can only appear at one certain level in a configuration file For in the old configuring model or in a default ParameterSet object inside of a or in a category or in a severity object This layout of multiple sources for a single parameter brings great confusion in both writing a configuration and in processing the configuration file Under the new the only allowed place for the parameter limit to appear is inside of a category which is inside of a destination object Other improvements simplify the meaning of a destination name In the old a destination name has multiple folds of meanings the e cout and cerr have the special meaning of logging messages to standard output or standard error the name also serves as the output filename if the destination is a file these names are also references to look up for detailed configurations in configuring the MessageFacility The multi purpose of the destination name might cause some unwanted behavior in either writing or parsing the configuration file To amend in the new model the destination name is now merely a name for a which might represent the literal purpose of this or just an id All other meanings of the destinations names now go into the destination ParameterSet as individual such as the type parameter and filename parameter Following is the deatiled rule for the new configuring Everything that is related with MessageFacility configuration must be wrapped in a single ParameterSet object with the name MessageFacility The MessageFacility ParameterSet object contains a series of top level parameters These parameters can be chosen a vector of string listing the name of debug enabled models Or use *to enable debug messages in all modules a vector of string a vector of string a vector of string a ParameterSet object containing the list of all destinations The destinations ParameterSet object is a combination of ParameterSet objects for individual destinations There are two types of destinations that you can insert in the destinations ParameterSet ordinary including cout
Definition: MessageFacilityAsStandAlonePackage.txt:82
genie::flux::GSimpleNtpMeta::protons
Double_t protons
represented number of protons-on-target
Definition: GSimpleNtpFlux.h:176
genie::flux::GSimpleNtpMeta::seed
Int_t seed
random seed used in generation
Definition: GSimpleNtpFlux.h:186
genie::flux::GFluxFileConfigI::SetUpstreamZ
virtual void SetUpstreamZ(double z0)
Definition: GFluxFileConfigI.cxx:74
genie::flux::GSimpleNtpMeta::windowDir2
Double_t windowDir2[3]
dx,dy,dz of window direction 2
Definition: GSimpleNtpFlux.h:180
genie::flux::GNuMIFlux
A GENIE flux driver encapsulating the NuMI neutrino flux. It reads-in the official GNUMI neutrino flu...
Definition: GNuMIFlux.h:217
genie::flux::GNuMIFlux::Weight
double Weight(void)
returns the flux neutrino weight (if any)
Definition: GNuMIFlux.h:234
genie::flux::GNuMIFlux::LoadBeamSimData
virtual void LoadBeamSimData(const std::vector< std::string > &filenames, const std::string &det_loc)
Definition: GNuMIFlux.cxx:640
genie::flux::GSimpleNtpMeta::pdglist
std::vector< Int_t > pdglist
list of neutrino flavors
Definition: GSimpleNtpFlux.h:171
genie::flux::GSimpleNtpMeta::windowBase
Double_t windowBase[3]
x,y,z position of window base point
Definition: GSimpleNtpFlux.h:178
genie::flux::GNuMIFlux::GetFileList
std::vector< std::string > GetFileList()
list of files currently part of chain
Definition: GNuMIFlux.cxx:2592
const_iterator
intermediate_table::const_iterator const_iterator
Definition: intermediate_table.cc:25
genie::flux::GNuMIFlux::PassThroughInfo
const GNuMIFluxPassThroughInfo & PassThroughInfo(void)
GNuMIFluxPassThroughInfo.
Definition: GNuMIFlux.h:252
genie::flux::GSimpleNtpEntry::vtxz
Double_t vtxz
z position in lab frame
Definition: GSimpleNtpFlux.h:74
genie::flux::GSimpleNtpMeta::minWgt
Double_t minWgt
minimum weight
Definition: GSimpleNtpFlux.h:174
genie::flux::GSimpleNtpNuMI::ppmedium
Int_t ppmedium
tracking medium where parent was produced
Definition: GSimpleNtpFlux.h:118
genie::flux::GSimpleNtpMeta::infiles
std::vector< std::string > infiles
list of input files
Definition: GSimpleNtpFlux.h:184
genie::flux::GSimpleNtpNuMI::tpx
Double_t tpx
parent particle px at target exit
Definition: GSimpleNtpFlux.h:103
genie::GFluxI::GenerateNext
virtual bool GenerateNext(void)=0
generate the next flux neutrino (return false in err)
genie::flux::GNuMIFlux::GetDecayDist
double GetDecayDist() const
dist (user units) from dk to current pos
Definition: GNuMIFlux.cxx:551
genie::flux::GSimpleNtpMeta::maxWgt
Double_t maxWgt
maximum weight
Definition: GSimpleNtpFlux.h:175
genie::flux::GSimpleNtpMeta::metakey
UInt_t metakey
index key to tie to individual entries
Definition: GSimpleNtpFlux.h:187
genie::flux::GSimpleNtpEntry::vtxx
Double_t vtxx
x position in lab frame
Definition: GSimpleNtpFlux.h:72
genie::flux::GSimpleNtpEntry::pz
Double_t pz
z momentum in lab frame
Definition: GSimpleNtpFlux.h:79
genie::GFluxI::PdgCode
virtual int PdgCode(void)=0
returns the flux neutrino pdg code
genie::GFluxI::Momentum
virtual const TLorentzVector & Momentum(void)=0
returns the flux neutrino 4-momentum
genie::flux::GNuMIFlux::PrintConfig
void PrintConfig()
print the current configuration
Definition: GNuMIFlux.cxx:2492
genie::flux::GSimpleNtpEntry::dist
Double_t dist
distance from hadron decay
Definition: GSimpleNtpFlux.h:75
file
< separator(=)> module_type Type Source location< separator(-)> DummyAnalyzer analyzer< path > DummyAnalyzer_module cc DummyFilter filter< path > DummyFilter_module cc *DummyProducer producer< path > DummyProducer_module cc *DummyProducer producer< path > DummyProducer_module cc< separator(=)> The modules marked *above are degenerate i e specifying the short module_type value leads to an ambiguity In order to use a degenerate in your configuration file
Definition: BasicOptions_03-ref.txt:15
genie::flux::GNuMIFlux::NFluxNeutrinos
virtual long int NFluxNeutrinos() const
of rays generated
Definition: GNuMIFlux.cxx:288
genie::flux::GSimpleNtpEntry::metakey
UInt_t metakey
key to meta data
Definition: GSimpleNtpFlux.h:83
genie::flux::GSimpleNtpEntry::py
Double_t py
y momentum in lab frame
Definition: GSimpleNtpFlux.h:78
genie::flux::GSimpleNtpMeta::auxdblname
std::vector< std::string > auxdblname
tagname of aux doubles associated w/ entry
Definition: GSimpleNtpFlux.h:183
genie::GFluxI
GENIE Interface for user-defined flux classes.
Definition: GFluxI.h:37
genie::flux::GNuMIFlux::UsedPOTs
double UsedPOTs(void) const
of protons-on-target used
Definition: GNuMIFlux.cxx:620
void gnumi2simple ( string  fname,
string  det,
string  flux,
long int  nentries = 5000,
double  pots = 1.0e30,
double  enumin = 0.0,
bool  doaux = false 
)

Definition at line 31 of file gnumi2simple.C.

34 {
35 
36  cout << "output file: " << fname << endl;
37  cout << "request " << nentries << " entries, or " << pots << " POTs" << endl;
38 
39  string cfg = "MINOS-NearDet";
40  if ( det == "minos" ) cfg = "MINOS-NearDet";
41  if ( det == "novaipnd" ) cfg = "NOvA-IPNDShed";
42  if ( det == "uboone" ) cfg = "microboone-numi";
43 
44  cout << "det \"" << det << "\" ==> cfg " << cfg << endl;
45 
46  int flxset = 0;
47  if ( flux == "old" ) flxset = 0;
48  if ( flux == "old1" ) flxset = 10;
49  if ( flux == "old2" ) flxset = 20;
50  if ( flux == "lowcut" ) flxset = 1;
51  if ( flux == "lowcut1" ) flxset = 11;
52  if ( flux == "lowcut2" ) flxset = 21;
53  if ( flux == "flugg" ) flxset = 2;
54  if ( flux == "flugg1" ) flxset = 12;
55  if ( flux == "flugg2" ) flxset = 22;
56  if ( flux == "flugglowth" ) flxset = 3;
57  if ( flux == "flugglowth1" ) flxset = 13;
58  if ( flux == "flugglowth2" ) flxset = 23;
59 
60  cout << "flux \"" << flux << "\" ==> flxset " << flxset << endl;
61 
62  fill_simple(fname,cfg,flxset,nentries,pots,enumin,doaux);
63 }
fill_simple
void fill_simple(string fname="fluxntgenie.root", string cfg="MINOS-Near", int flxset=0, long int nentries=0, double pots=1.0e30, double enumin=0.0, bool doaux=false)
Definition: gnumi2simple.C:65
dbjson.cfg
cfg
Definition: dbjson.py:29
cout
the ParameterSet object passed in for the configuration of a destination should be the only source that can affect the behavior of that destination This is to eliminate the dependencies of configuring a destination from multiple mostly from the defaults It suppresses possible glitches about changing the configuration file somewhere outside of a destination segament might still affect the behavior of that destination In the previous configuration for a specific the value of a certain e may come from following and have been suppressed It the configuring ParameterSet object for each destination will be required to carry a parameter list as complete as possible If a parameter still cannot be found in the ParameterSet the configuration code will go look for a hardwired default directly The model is a great simplicity comparing with the previous especially when looking for default values Another great advantage is most of the parameters now have very limited places that allows to appear Usually they can only appear at one certain level in a configuration file For in the old configuring model or in a default ParameterSet object inside of a or in a category or in a severity object This layout of multiple sources for a single parameter brings great confusion in both writing a configuration and in processing the configuration file Under the new the only allowed place for the parameter limit to appear is inside of a category which is inside of a destination object Other improvements simplify the meaning of a destination name In the old a destination name has multiple folds of meanings the e cout and cerr have the special meaning of logging messages to standard output or standard error the name also serves as the output filename if the destination is a file these names are also references to look up for detailed configurations in configuring the MessageFacility The multi purpose of the destination name might cause some unwanted behavior in either writing or parsing the configuration file To amend in the new model the destination name is now merely a name for a which might represent the literal purpose of this or just an id All other meanings of the destinations names now go into the destination ParameterSet as individual such as the type parameter and filename parameter Following is the deatiled rule for the new configuring Everything that is related with MessageFacility configuration must be wrapped in a single ParameterSet object with the name MessageFacility The MessageFacility ParameterSet object contains a series of top level parameters These parameters can be chosen a vector of string listing the name of debug enabled models Or use *to enable debug messages in all modules a vector of string a vector of string a vector of string a ParameterSet object containing the list of all destinations The destinations ParameterSet object is a combination of ParameterSet objects for individual destinations There are two types of destinations that you can insert in the destinations ParameterSet ordinary including cout
Definition: MessageFacilityAsStandAlonePackage.txt:82
void make_simple_grid ( string  fname,
string  det,
string  flux,
long int  nentries = 5000,
double  pots = 1.0e30,
double  enumin = 0.0,
bool  doaux = false 
)

Definition at line 296 of file gnumi2simple.C.

299 {
300  gnumi2simple(fname,det,flux,nentries,pots,enumin,doaux);
301 }
gnumi2simple
void gnumi2simple(string fname, string det, string flux, long int nentries=5000, double pots=1.0e30, double enumin=0.0, bool doaux=false)
Definition: gnumi2simple.C:31
Generated by   doxygen 1.8.11