---

ReactorFsim.cpp

Go to the documentation of this file.

#include "ReactorFortran.hh"
#include "ReactorConstants.hh"
#include "ReactorEvent.hh"
#include "ReactorDetector.hh"
#include "ReactorXC.hh"
#include "ReactorNtuple.hh"
#include <iostream.h>
#include <math.h>
#include <string.h>

int main(){

  //  ReactorNeutrons Neutrons;
  //
  //  Set number of events, paw file name, and ntuple name
  //
  int nevents = 100000;
  double rmaxgen=225;
  int fast=0;
  char* fileName="/local/data14/reactorfsim/cf252_pmtradon_100k.ntuple";
  char* ntName="event data";
  //
  // init random number seed
  //
  int luxury = 2;
  int seed = 8796896;
  int k1=0;
  int k2=0;
  rluxgo_(luxury,seed,k1,k2);
  //
  // Paw setup.  
  //
  hbook_init__();
  //
  int iostat;
  int ntLun=99;
  int ntBuf=1024;
  char* dirName = "TEST";
  char* chopt = "N";
  hropen_(ntLun,dirName,fileName,chopt,ntBuf,iostat,
          strlen(dirName),strlen(fileName),strlen(chopt));
  //
  // Set up simple detector models
  // Set up simple Ntuple and histogranms 
  // See ReactorNtuple.hh for ntuple/histogram details.
  //
  ReactorDetector DetectorA;
  DetectorA.SetParam_FastNeutronOption(fast);
  DetectorA.SetParam_R0(175);
  DetectorA.SetParam_R1(200);
  DetectorA.SetParam_R2(275);
  ReactorNtuple ModelA(1,"model A",dirName);
  //
  // event loop
  //
  for (int events = 0;events<nevents;events++){
    //
    // create neutrino event
    //
    ReactorEvent Event(nevents,rmaxgen);
    //
    // create fast detector simulation
    //
    DetectorA.LightsOut(Event);
    ModelA.Fill(Event,DetectorA);
    //
  }
  //
  // close up paw
  //
  hldir_(" "," ",strlen(" "),strlen(" "));
  char dirName2[strlen(dirName)+2];
  *dirName2='\0';
  strcat(dirName2,"//");
  strcat(dirName2,dirName);
  hcdir_(dirName2," ",strlen(dirName2),strlen(" "));
  int outopt = 0;
  int icycle;
  hrout_(outopt,icycle," ",strlen(" "));
  hrend_(dirName,strlen(dirName));
}

---