In This Package:

MuonTree.cc

Go to the documentation of this file.

#include "MuonTree.h"

#include <sstream>

MuonTree::MuonTree()
{}

MuonTree::~MuonTree()
{}

int MuonTree::create(const char* filename)
{
  cout<<"TTree Created"<<endl;
  
  m_file=new TFile(filename,"recreate");
  m_tree=new TTree("t","Analysis Tree");    

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

  m_tree->Branch("nhits_RPC",&nhits_RPC,"nhits_RPC/I");
  m_tree->Branch("nhits_AD1",&nhits_AD1,"nhits_AD1/I");
  m_tree->Branch("nhits_AD2",&nhits_AD2,"nhits_AD2/I");
  m_tree->Branch("nhits_IWS",&nhits_IWS,"nhits_IWS/I");
  m_tree->Branch("nhits_OWS",&nhits_OWS,"nhits_OWS/I");

  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_r",   &trk1_r,  "trk1_r/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_phi", &trk1_phi,"trk1_phi/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("npe",npe,"npe[8][24]/D");   // npe                                                                                 
  m_tree->Branch("tpe",tpe,"tpe[8][24]/D");   // tpe 

  m_tree->Branch("TotDeInOil1", &TotDeInOil1, "TotDeInOil1/D");
  m_tree->Branch("TotDeInOil2", &TotDeInOil2, "TotDeInOil2/D");
  m_tree->Branch("TotDeInLso1", &TotDeInLso1, "TotDeInLso1/D");
  m_tree->Branch("TotDeInLso2", &TotDeInLso2, "TotDeInLso2/D");

  m_tree->Branch("MuDeInOil1",  &MuDeInOil1,  "MuDeInOil1/D"); 
  m_tree->Branch("MuDeInOil2",  &MuDeInOil2,  "MuDeInOil2/D"); 
  m_tree->Branch("MuDeInLso1",  &MuDeInLso1,  "MuDeInLso1/D"); 
  m_tree->Branch("MuDeInLso2",  &MuDeInLso2,  "MuDeInLso2/D"); 

  m_tree->Branch("MuDxInOil1",  &MuDxInOil1,  "MuDxInOil1/D"); 
  m_tree->Branch("MuDxInOil2",  &MuDxInOil2,  "MuDxInOil2/D"); 
  m_tree->Branch("MuDxInLso1",  &MuDxInLso1,  "MuDxInLso1/D"); 
  m_tree->Branch("MuDxInLso2",  &MuDxInLso2,  "MuDxInLso2/D"); 

  m_tree->Branch("his_t1",   &his_t1,  "his_t1/D");
  m_tree->Branch("his_x1",   &his_x1,  "his_x1/D");
  m_tree->Branch("his_y1",   &his_y1,  "his_y1/D");
  m_tree->Branch("his_z1",   &his_z1,  "his_z1/D");

  m_tree->Branch("oil1_p_x",  &oil1_p_x, "oil1_p_x/D");
  m_tree->Branch("oil1_p_y",  &oil1_p_y, "oil1_p_y/D");
  m_tree->Branch("oil1_p_z",  &oil1_p_z, "oil1_p_z/D");
  m_tree->Branch("oil1_e_x",  &oil1_e_x, "oil1_e_x/D");
  m_tree->Branch("oil1_e_y",  &oil1_e_y, "oil1_e_y/D");
  m_tree->Branch("oil1_e_z",  &oil1_e_z, "oil1_e_z/D");
  m_tree->Branch("oil1_p_px",  &oil1_p_px, "oil1_p_px/D");
  m_tree->Branch("oil1_p_py",  &oil1_p_py, "oil1_p_py/D");
  m_tree->Branch("oil1_p_pz",  &oil1_p_pz, "oil1_p_pz/D");

  m_tree->Branch("oil2_p_x",  &oil2_p_x, "oil2_p_x/D");
  m_tree->Branch("oil2_p_y",  &oil2_p_y, "oil2_p_y/D");
  m_tree->Branch("oil2_p_z",  &oil2_p_z, "oil2_p_z/D");
  m_tree->Branch("oil2_e_x",  &oil2_e_x, "oil2_e_x/D");
  m_tree->Branch("oil2_e_y",  &oil2_e_y, "oil2_e_y/D");
  m_tree->Branch("oil2_e_z",  &oil2_e_z, "oil2_e_z/D");
  m_tree->Branch("oil2_p_px",  &oil2_p_px, "oil2_p_px/D");
  m_tree->Branch("oil2_p_py",  &oil2_p_py, "oil2_p_py/D");
  m_tree->Branch("oil2_p_pz",  &oil2_p_pz, "oil2_p_pz/D");


  m_tree->Branch("phiRec",&phiRec,"phiRec/D");   // recontructed phi
  //
  m_tree->Branch("oil1_p_x_rec",  &oil1_p_x_rec, "oil1_p_x_rec/D");
  m_tree->Branch("oil1_p_y_rec",  &oil1_p_y_rec, "oil1_p_y_rec/D");
  m_tree->Branch("oil1_p_z_rec",  &oil1_p_z_rec, "oil1_p_z_rec/D");
  //
  m_tree->Branch("oil2_p_x_rec",  &oil2_p_x_rec, "oil2_p_x_rec/D");
  m_tree->Branch("oil2_p_y_rec",  &oil2_p_y_rec, "oil2_p_y_rec/D");
  m_tree->Branch("oil2_p_z_rec",  &oil2_p_z_rec, "oil2_p_z_rec/D");

  // Number of neutrons
  m_tree->Branch("nNeutron",      &nNeutron,     "nNeutron/I");
  m_tree->Branch("nInOws",      &nInOws,     "nInOws/I");
  m_tree->Branch("nInOil1",      &nInOil1,     "nInOil1/I");
  m_tree->Branch("nInOil2",      &nInOil2,     "nInOil2/I");
  m_tree->Branch("nInLso1",      &nInLso1,     "nInLso1/I");
  m_tree->Branch("nInLso2",      &nInLso2,     "nInLso2/I");


  return  1; // 1 for SUCCESS;
}

int MuonTree::reset()
{
  run=0;
  evt=0;

  nhits_RPC=0;
  nhits_AD1=0;
  nhits_AD2=0;
  nhits_IWS=0;
  nhits_OWS=0;

  trk1_pdgId=0;

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

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

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

  TotDeInOil1=0;
  TotDeInOil2=0;
  TotDeInLso1=0;
  TotDeInLso2=0;

  MuDeInOil1=0;
  MuDeInOil2=0;
  MuDeInLso1=0;
  MuDeInLso2=0;

  MuDxInOil1=0;
  MuDxInOil2=0;
  MuDxInLso1=0;
  MuDxInLso2=0;


  oil1_p_x=0;
  oil1_p_y=0;
  oil1_p_z=0;

  oil1_e_x=0;
  oil1_e_y=0;
  oil1_e_z=0;

  oil1_p_px=0;
  oil1_p_py=0;
  oil1_p_pz=0;

  oil2_p_x=0;
  oil2_p_y=0;
  oil2_p_z=0;

  oil2_e_x=0;
  oil2_e_y=0;
  oil2_e_z=0;

  oil2_p_px=0;
  oil2_p_py=0;
  oil2_p_pz=0;

  phiRec=0;

  oil1_p_x_rec=0;
  oil1_p_y_rec=0;
  oil1_p_z_rec=0;

  oil2_p_x_rec=0;
  oil2_p_y_rec=0;
  oil2_p_z_rec=0;

  for(int ring=0;ring<8;ring++) {
    for(int column=0;column<24;column++) {
      npe[ring][column]=0;  // something invalid                                                                                     
      tpe[ring][column]=99999;  // something invalid                                                                                 
    }
  }

  nNeutron = 0;
  nInOws = 0;
  nInOil1 = 0;
  nInOil2 = 0;
  nInLso1 = 0;
  nInLso2 = 0;

  return 1;  // 1 for SUCCESS;
}

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

  return  1; // 1 for SUCCESS;
}

int MuonTree::close()
{
  m_file->Write();

  return  1; // 1 for SUCCESS;
}

string MuonTree::IntToString(int intValue) {

  ostringstream oss (ostringstream::out);

  oss <<intValue;

  string strRetVal=oss.str();

  return(strRetVal);
}

---