SimAna.h

Go to the documentation of this file.

//Constructed by Josh Brown.
//brown.ja@berkeley.edu
#ifndef _SIM_ANA_H_
#define _SIM_ANA_H_
 
#include "PostProcBase.h"
#include "SimulationDataStructures.h"
 
#include "TH1.h"
 
#include <string> 
namespace NSDG4
{
class SimAna : public PostProcBase<FullEvent> 
{
public:
  void printTrackInfo(int a_entry);
 
  std::ostream& printEvents(int a_num =1,
                            std::vector<std::string> a_partList = {},
                            std::ostream& a_stream = std::cout
                            );
 
 
  void accumulateFirstStepInfo(std::string a_lvName,
                               std::string a_particleName =""
                               );
  void vertexStepInfoTree(std::string a_lvName,
                          std::string a_particleName =""
                         );
  void createParticleTTree(std::string a_lvName,
                           std::vector<std::string> a_ignoreList
                           );
  void createEnergyDepositionTree(std::string a_lvName,
                                  std::vector<std::string> a_particleNames
                                 );                                 
  void createTotalEnergyDepositionTree(std::string a_lvName);
 
  //MCGUIRE
  void investigateMultipleScatters(std::string a_lvName);
  //MCGUIRE
  
// #ifdef OSLY_PRESENT
  void createLightOutputTree(std::string a_lvName,
                                std::vector<std::string> a_particleNames,
                                std::string a_scintName = "EJ309"
                               );
// #endif
  void createPartTimeEDTree(std::string a_lvName,
                            double a_tCutWindow
                            );
  void createFluxTree(std::string a_lvName,
                      std::string a_particleName,
                      std::vector<double> a_fluxBoundaryList,
                      int a_dim,
                      std::string a_outfileName= "fluxBoundaryTree.root"
                      );
 
  TH1* getTotalEnergyHistory();
  FullEvent getFullEvent(int a_eventNumber);
 
  virtual ~SimAna();
private:
  int m_lastEvent;
  void setTreeName();
  void setBranchName();
 
};
 
 
}
#endif