---

ReactorTree.cpp

Go to the documentation of this file.

#include <iostream.h>
#include <math.h>
#include <string.h>
#include <cstdio>
#include "ReactorEvent.hh"
#include "ReactorDetector.hh"
#include "ReactorFortran.hh"
#include "ReactorTree.hh"
#include "RMCEvent.hh"
#include "TTree.h"
#include "TBranch.h"

ReactorTree::ReactorTree(){
}

ReactorTree::ReactorTree(const char* TreeName){

  // define Tree branches here
  //
  tree = new TTree(TreeName,"ReactorTree");
  tree->SetAutoSave(1000000);
  int bufsize =64000000;
  int split = 2;
  if (split) bufsize /=4 ;
  int branchstyle=1;
  TTree::SetBranchStyle(branchstyle);
  rmce = new RMCEvent();
  TBranch *branch=tree->Branch("RMCEvent","RMCEvent",&rmce,bufsize,split);
//  tree->Show();
  branch->SetAutoDelete(kFALSE); 
}

ReactorTree::~ReactorTree(){
  ;
}

void ReactorTree::Fill(ReactorEvent& Event,ReactorDetector& Detector){

    static const int Ndim=7;
// Harvest information from ReactorEvent 
    std::string EventType = Event.eventType;
    int Nevent = Event.GetNevent();
    rmce->Build(Event,Detector,EventType,Nevent);
/*    if (EventType == "reactor"){   
    }
    else if(EventType == "cf252"){   
            rmce->SetEventType(1);
    }
    else if(EventType == "pmtrad"){   
            rmce->SetEventType(2);
    }
    else if(EventType == "muon"){   
            rmce->SetEventType(3);
    }
    else if(EventType == "all"){   
            rmce->SetEventType(4);
    }
    */
    tree->Fill();

}

void ReactorTree::CloseOut(){

  tree->Print();

}

---