ValidationTree.cc

Go to the documentation of this file.

#include "ValidationTree.h"

using namespace std;

ValidationTree::ValidationTree()
{}

ValidationTree::~ValidationTree()
{}

int ValidationTree::create()
{
  m_file=new TFile("monitor.root","recreate");
  m_tree=new TTree("t","Analysis Tree");    

  m_tree->Branch("nhits",&SimHitsEntries,"SimHitsEntries/I");
  m_tree->Branch("hitTime",hitTime,"hitTime[SimHitsEntries]/D");
  m_tree->Branch("hitx",hitx,"hitx[SimHitsEntries]/D");
  m_tree->Branch("hity",hity,"hity[SimHitsEntries]/D");
  m_tree->Branch("hitz",hitz,"hitz[SimHitsEntries]/D");
  m_tree->Branch("sensID",sensID,"sensID[SimHitsEntries]/I");
  m_tree->Branch("ancestor",ancestorPdg,"ancestor[SimHitsEntries]/I");
  m_tree->Branch("wavelength",wavelength,"wavelength[SimHitsEntries]/D");
  m_tree->Branch("polx",polx,"polx[SimHitsEntries]/D");
  m_tree->Branch("poly",poly,"poly[SimHitsEntries]/D");
  m_tree->Branch("polz",polz,"polz[SimHitsEntries]/D");
  m_tree->Branch("px",px,"px[SimHitsEntries]/D");
  m_tree->Branch("py",py,"py[SimHitsEntries]/D");
  m_tree->Branch("pz",pz,"pz[SimHitsEntries]/D");
  m_tree->Branch("weight",weight,"weight[SimHitsEntries]/F");

  //Muon vertex Info
  m_tree->Branch("MuonE",&MuonE,"MuonE/D");
  m_tree->Branch("MuonX",&MuonX,"MuonX/D");
  m_tree->Branch("MuonY",&MuonY,"MuonY/D");
  m_tree->Branch("MuonZ",&MuonZ,"MuonZ/D");
  m_tree->Branch("MuonPx",&MuonPx,"MuonPx/D");
  m_tree->Branch("MuonPy",&MuonPy,"MuonPy/D");
  m_tree->Branch("MuonPz",&MuonPz,"MuonPz/D");
  m_tree->Branch("MuonTime",&MuonTime,"MuonTime/D");
  m_tree->Branch("StopFlag",&StopFlag,"StopFlag/I");


  m_tree->Branch("MuonTrkLengthInOws",&MuonTrkLengthInOws,"MuonTrkLengthInOws/D");
  m_tree->Branch("MuonTrkLengthInIws",&MuonTrkLengthInIws,"MuonTrkLengthInIws/D");
  m_tree->Branch("MuonTrkLengthInLS",&MuonTrkLengthInLS,"MuonTrkLengthInLS/D");
  m_tree->Branch("MuonTrkLengthInGdLS",&MuonTrkLengthInGdLS,"MuonTrkLengthInGdLS/D");
  m_tree->Branch("dEInn",&dEInn,"dEInn/D");
  m_tree->Branch("dEOut",&dEOut,"dEOut/D");

  m_tree->Branch("trk1_pdgId",  &trk1_pdgId, "trk1_pdgId/I");  
                                                                              
  m_tree->Branch("trk1_t",   &trk1_t,  "trk1_t/D");   
  m_tree->Branch("trk1_x",   &trk1_x,  "trk1_x/D");   
  m_tree->Branch("trk1_y",   &trk1_y,  "trk1_y/D");   
  m_tree->Branch("trk1_z",   &trk1_z,  "trk1_z/D");   
                                                                              
  m_tree->Branch("trk1_e",   &trk1_e,  "trk1_e/D");   
  m_tree->Branch("trk1_px",  &trk1_px, "trk1_px/D");  
  m_tree->Branch("trk1_py",  &trk1_py, "trk1_py/D");  
  m_tree->Branch("trk1_pz",  &trk1_pz, "trk1_pz/D");  
                                                                              
  m_tree->Branch("trk1_ke",  &trk1_ke, "trk1_ke/D");   // kinematic energy
  m_tree->Branch("trk1_ve",  &trk1_ve, "trk1_ve/D");   // visible energy
                                                                              
  m_tree->Branch("trk1_TrkLength_GD", &trk1_TrkLength_GD, "trk1_TrkLength_GD/D");  // pathlength in GDLS
  m_tree->Branch("trk1_TrkLength_iAV",&trk1_TrkLength_iAV,"trk1_TrkLength_iAV/D"); // pathlength in inner AV
  m_tree->Branch("trk1_TrkLength_LS", &trk1_TrkLength_LS, "trk1_TrkLength_LS/D");  // pathlength in LS
  m_tree->Branch("trk1_TrkLength_oAV",&trk1_TrkLength_oAV,"trk1_TrkLength_oAV/D"); // pathlength in outer AV
  m_tree->Branch("trk1_TrkLength_Oil",&trk1_TrkLength_Oil,"trk1_TrkLength_Oil/D"); // pathlength in oil

  m_tree->Branch("trk2_pdgId",  &trk2_pdgId,  "trk2_pdgId/I");

  m_tree->Branch("trk2_t",   &trk2_t,   "trk2_t/D");
  m_tree->Branch("trk2_x",   &trk2_x,   "trk2_x/D");
  m_tree->Branch("trk2_y",   &trk2_y,   "trk2_y/D");
  m_tree->Branch("trk2_z",   &trk2_z,   "trk2_z/D");

  m_tree->Branch("trk2_e",   &trk2_e,   "trk2_e/D");
  m_tree->Branch("trk2_px",  &trk2_px,  "trk2_px/D");
  m_tree->Branch("trk2_py",  &trk2_py,  "trk2_py/D");
  m_tree->Branch("trk2_pz",  &trk2_pz,  "trk2_pz/D");

  m_tree->Branch("trk2_ke",  &trk2_ke,  "trk2_ke/D");   // kinematic energy
  m_tree->Branch("trk2_ve",  &trk2_ve,  "trk2_ve/D");   // visible energy

  m_tree->Branch("trk2_TrkLength_GD", &trk2_TrkLength_GD, "trk2_TrkLength_GD/D");  // pathlength in GDLS
  m_tree->Branch("trk2_TrkLength_iAV",&trk2_TrkLength_iAV,"trk2_TrkLength_iAV/D"); // pathlength in inner AV
  m_tree->Branch("trk2_TrkLength_LS", &trk2_TrkLength_LS, "trk2_TrkLength_LS/D");  // pathlength in LS
  m_tree->Branch("trk2_TrkLength_oAV",&trk2_TrkLength_oAV,"trk2_TrkLength_oAV/D"); // pathlength in outer AV
  m_tree->Branch("trk2_TrkLength_Oil",&trk2_TrkLength_Oil,"trk2_TrkLength_Oil/D"); // pathlength in oil

  // initial optical photons
  m_tree->Branch("number_opti_pho", &number_opti_pho, "number_opti_pho/I");
  m_tree->Branch("wavelength_opti_pho",wavelength_opti_pho,",wavelength_opti_pho[number_opti_pho]/D");

  //
  m_tree->Branch("n",       &n,       "n/I");
  m_tree->Branch("execKine",&execKine,"execKine/I");   // execution number in kinematic
  m_tree->Branch("execDets",&execDets,"execDets/I");   // in DetSim                    
  m_tree->Branch("execElec",&execElec,"execElec/I");   // in ElecSim                   
  m_tree->Branch("execTrig",&execTrig,"execTrig/I");   // in TrigRead                  
  m_tree->Branch("execSing",&execSing,"execSing/I");   // in SingleLoader              

  m_tree->Branch("time",&time,"time/D"); // arrival time in second

  // readout 
  m_tree->Branch("adc",adc,"adc[8][24]/I");   // adc
  m_tree->Branch("tdc",tdc,"tdc[8][24]/I");   // tdc

  m_tree->Branch("adc_sum",&adc_sum,"adc_sum/D");   // adc_sum


  return  1; // 1 for SUCCESS;
}

int ValidationTree::reset()
{
  SimHitsEntries=0;

  //muon vertex information.
  MuonE=0;
  MuonX=0;
  MuonY=0;
  MuonZ=0;
  MuonPx=0;
  MuonPy=0;
  MuonPz=0;
  MuonTime=0;
  StopFlag=0;
  
  MuonTrkLengthInOws=0;
  MuonTrkLengthInIws=0;
  MuonTrkLengthInLS=0;
  MuonTrkLengthInGdLS=0;
  dEInn=0;
  dEOut=0;
  

  trk1_pdgId=0;

  trk1_t=0;
  trk1_x=0;
  trk1_y=0;
  trk1_z=0;

  trk1_e=0;
  trk1_px=0;
  trk1_py=0;
  trk1_pz=0;

  trk1_ke=0;  // kinematic energy
  trk1_ve=0;  // visible energy

  trk1_TrkLength_GD=0;  // pathlength in GDLS
  trk1_TrkLength_iAV=0; //            in inner AV
  trk1_TrkLength_LS=0;  //            in LS
  trk1_TrkLength_oAV=0; //            in outer AV
  trk1_TrkLength_Oil=0; //            in oil

  // MC truth for the second primary track ///////////
  trk2_pdgId=0;

  trk2_t=0;
  trk2_x=0;
  trk2_y=0;
  trk2_z=0;

  trk2_e=0;
  trk2_px=0;
  trk2_py=0;
  trk2_pz=0;

  trk2_ke=0;  // kinematic energy
  trk2_ve=0;  // visible energy

  trk2_TrkLength_GD=0;  // pathlength in GDLS
  trk2_TrkLength_iAV=0; //            in inner AV
  trk2_TrkLength_LS=0;  //            in LS
  trk2_TrkLength_oAV=0; //            in outer AV
  trk2_TrkLength_Oil=0; //            in oil


  number_opti_pho=0;

  n=0;
  execKine=0;
  execDets=0;
  execElec=0;
  execTrig=0;
  execSing=0;

  time =0;

  adc_sum=0;

  return 1;  // 1 for SUCCESS;
}

int ValidationTree::fill()
{
  m_tree->Fill();

  return  1; // 1 for SUCCESS;
}

int ValidationTree::close()
{
  m_file->Write();
  return  1; // 1 for SUCCESS;
}