BeamParams Class Reference

#include <BeamParams.h>

List of all members.

Public Member Functions
	BeamParams (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
ITHistSvc *	m_thsvc
int	m_sel_number [15]
int	m_pid
int	m_vertexFind
double	m_vz0Cut
double	m_cosThetaCut
double	m_minDistance
double	m_minPointX
double	m_minPointY
double	m_meanPointX
double	m_meanPointY
double	m_chisqCut
int	m_trackIteration
int	m_vertexIteration
double	m_chi2CutforTrkIter
double	m_chi2CutforSmooth
int	m_trackNumberCut
bool	m_useDedx
bool	m_useTof1
bool	m_useTof2
bool	m_useTofE
bool	m_useTofQ
bool	m_useEmc
bool	m_useMuc
double	m_PidProbCut
double	m_pivotX
double	m_pivotY
double	m_pivotZ
int	m_trkNum
int	m_runNo
int	m_parVer
std::string	m_fileNameDst
std::string	m_fileNameHadron
std::string	m_figsName
TH1D *	m_vertex_x
TH1D *	m_vertex_y
TH1D *	m_vertex_z
TH1D *	m_vertex_x_kal
TH1D *	m_vertex_y_kal
TH1D *	m_vertex_z_kal
HepVector	m_Cpuvc
double	m_cpu [5]
int	m_hadronFile
NTuple::Tuple *	m_tuple1
NTuple::Item< double >	m_xc
NTuple::Item< double >	m_yc
NTuple::Item< double >	m_zc
NTuple::Item< double >	m_mind
NTuple::Tuple *	m_tuple2
NTuple::Item< double >	m_chis
NTuple::Item< double >	m_chif
NTuple::Item< double >	m_probs
NTuple::Item< double >	m_probf
NTuple::Tuple *	m_tuple3
NTuple::Item< double >	m_chik
NTuple::Item< long >	m_ndofk
NTuple::Item< double >	m_probk
NTuple::Item< double >	m_kvx
NTuple::Item< double >	m_kvy
NTuple::Item< double >	m_kvz
NTuple::Item< long >	m_numTrk
NTuple::Tuple *	m_tuple4
NTuple::Item< double >	m_chig
NTuple::Item< long >	m_ndofg
NTuple::Item< double >	m_probg
NTuple::Item< double >	m_gvx
NTuple::Item< double >	m_gvy
NTuple::Item< double >	m_gvz
NTuple::Tuple *	m_tuple5
NTuple::Item< double >	m_pull_drho
NTuple::Item< double >	m_pull_phi
NTuple::Item< double >	m_pull_kapha
NTuple::Item< double >	m_pull_dz
NTuple::Item< double >	m_pull_lamb
NTuple::Item< double >	m_pull_momentum
NTuple::Tuple *	m_tuple6
NTuple::Item< double >	m_mdcTrk_x
NTuple::Item< double >	m_mdcTrk_y
NTuple::Item< double >	m_mdcTrk_z
NTuple::Item< double >	m_mdcTrk_r
NTuple::Item< double >	m_rxy
NTuple::Item< double >	m_mdcTrk_dr
NTuple::Item< double >	m_mdcKalTrk_z
NTuple::Tuple *	m_tuple7
NTuple::Item< double >	m_gpull_drho
NTuple::Item< double >	m_gpull_phi
NTuple::Item< double >	m_gpull_kapha
NTuple::Item< double >	m_gpull_dz
NTuple::Item< double >	m_gpull_lamb

---

Detailed Description

Definition at line 12 of file BeamParams.h.

---

Constructor & Destructor Documentation

BeamParams::BeamParams	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 44 of file BeamParams.cxx.

References m_chi2CutforSmooth, m_chi2CutforTrkIter, m_chisqCut, m_cosThetaCut, m_figsName, m_fileNameDst, m_fileNameHadron, m_hadronFile, m_meanPointX, m_meanPointY, m_minDistance, m_minPointX, m_minPointY, m_parVer, m_pid, m_trackIteration, m_trackNumberCut, m_useDedx, m_useEmc, m_useMuc, m_useTof1, m_useTof2, m_useTofE, m_useTofQ, m_vertexFind, m_vertexIteration, m_vz0Cut, and deljobs::string.

                                                                      :
  Algorithm (name, pSvcLocator)
{
  // Declare the properties
  //declareProperty("RunNo", m_runNo = 9999);

  declareProperty("ParVer", m_parVer = 1);
  declareProperty("TrackNumberCut", m_trackNumberCut = 1);
  declareProperty("Vz0Cut", m_vz0Cut = 50.0);
  declareProperty("CosThetaCut", m_cosThetaCut=0.93);

  // particle identification
  declareProperty("OptionPID", m_pid = 0);
  declareProperty("PidUseDedx", m_useDedx = true);
  declareProperty("PidUseTof1", m_useTof1 = true);
  declareProperty("PidUseTof2", m_useTof2 = true);
  declareProperty("PidUseTofE", m_useTofE = false);
  declareProperty("PidUseTofQ", m_useTofQ = false);
  declareProperty("PidUseEmc", m_useEmc = false);
  declareProperty("PidUseMuc", m_useMuc = false);

  // vertex finding
  declareProperty("OptionVertexFind", m_vertexFind = 0);
  declareProperty("MinDistance", m_minDistance = 100.0);
  declareProperty("MinPointX", m_minPointX = 0.1);
  declareProperty("MinPointY", m_minPointY = 0.1);
  declareProperty("MeanPointX", m_meanPointX = 0.1);
  declareProperty("MeanPointY", m_meanPointY = -0.25);

  // Kalman vertex fitting
  declareProperty("ChisqCut", m_chisqCut = 200.0);
  declareProperty("TrackIteration", m_trackIteration = 5);
  declareProperty("VertexIteration", m_vertexIteration = 5);
  declareProperty("Chi2CutforTrkIter", m_chi2CutforTrkIter = 0.1);
  declareProperty("Chi2CutforSmooth", m_chi2CutforSmooth = 200.0);

  // output file name
  declareProperty("HadronFile", m_hadronFile = 0);
  declareProperty("DstFileName", m_fileNameDst = std::string("dst.txt"));
  declareProperty("HadronFileName", m_fileNameHadron = std::string("hadron.txt"));
  declareProperty("FigsName", m_figsName = std::string("figs.ps"));
}

---

Member Function Documentation

StatusCode BeamParams::execute

(

)

Definition at line 262 of file BeamParams.cxx.

References KalmanVertexFit::addTrack(), TrackPool::AddTrack(), VertexFit::AddVertex(), DstMdcTrack::charge(), KalmanVertexFit::chiF(), KalmanVertexFit::chiS(), KalmanVertexFit::chisq(), VertexFit::chisq(), cos(), cross(), Bes_Common::DEBUG, check_raw_filter::dist, DstMdcKalTrack::err(), EventModel::EvtRec::EvtRecEvent, EventModel::EvtRec::EvtRecTrackCol, Bes_Common::FATAL, KalmanVertexFit::filter(), VertexFit::Fit(), DstMdcKalTrack::helix(), DstMdcTrack::helix(), genRecEmupikp::i, Bes_Common::INFO, KalmanVertexFit::init(), VertexFit::init(), KalmanVertexFit::initVertex(), KalmanVertexFit::instance(), VertexFit::instance(), ParticleID::instance(), ganga-rec::j, m_chi2CutforSmooth, m_chi2CutforTrkIter, m_chif, m_chig, m_chik, m_chis, m_chisqCut, m_cosThetaCut, m_gpull_drho, m_gpull_dz, m_gpull_kapha, m_gpull_lamb, m_gpull_phi, m_gvx, m_gvy, m_gvz, m_kvx, m_kvy, m_kvz, m_mdcKalTrk_z, m_mdcTrk_dr, m_mdcTrk_r, m_mdcTrk_x, m_mdcTrk_y, m_mdcTrk_z, m_mind, m_minDistance, m_ndofg, m_ndofk, m_pid, m_probf, m_probg, m_probk, m_probs, m_pull_drho, m_pull_dz, m_pull_kapha, m_pull_lamb, m_pull_momentum, m_pull_phi, m_runNo, m_sel_number, m_trackIteration, m_trackNumberCut, m_tuple1, m_tuple2, m_tuple3, m_tuple4, m_tuple5, m_tuple6, m_tuple7, m_useDedx, m_useEmc, m_useMuc, m_useTof1, m_useTof2, m_useTofE, m_useTofQ, m_vertex_x, m_vertex_x_kal, m_vertex_y, m_vertex_y_kal, m_vertex_z, m_vertex_z_kal, m_vertexFind, m_vertexIteration, m_vz0Cut, m_xc, m_yc, m_zc, max, HTrackParameter::minDistanceTwoHelix(), msgSvc(), KalmanVertexFit::ndof(), pid, DstMdcKalTrack::pion, KalmanVertexFit::pull(), VertexFit::pull(), KalmanVertexFit::pullmomentum(), DstMdcTrack::r(), KalmanVertexFit::remove(), KalmanVertexFit::setChisqCut(), VertexParameter::setEvx(), DstMdcKalTrack::setPidType(), KalmanVertexFit::setTrackIteration(), KalmanVertexFit::setTrackIterationCut(), KalmanVertexFit::setVertexIteration(), VertexParameter::setVx(), KalmanVertexFit::smooth(), DstMdcTrack::theta(), VertexFit::Vx(), KalmanVertexFit::x(), DstMdcTrack::x(), xmass, DstMdcTrack::y(), and DstMdcTrack::z().

{
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in execute()" << endmsg;

  //  Read REC data
  SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
  SmartDataPtr<EvtRecEvent> recEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);
  SmartDataPtr<EvtRecTrackCol> recTrackCol(eventSvc(), EventModel::EvtRec::EvtRecTrackCol);
  log << MSG::DEBUG << "run and event = " << eventHeader->runNumber() 
                                        << " " << eventHeader->eventNumber() << endreq;
  log << MSG::DEBUG <<"ncharg, nneu, tottks = "
          << recEvent->totalCharged() << " , "
          << recEvent->totalNeutral() << " , "
      << recEvent->totalTracks() <<endreq;
  
  m_runNo = eventHeader->runNumber();
  //FIXME : cout
  if (eventHeader->eventNumber() % 5000 == 0) {
    cout << "Event number = " << eventHeader->eventNumber()<<" run number = "<< m_runNo << endl;
  }
  m_sel_number[0]++;

  // Cut 1 : Number of tracks. For anything sample, at least 3 charged tracks
  if (recEvent->totalCharged() < m_trackNumberCut) return StatusCode::SUCCESS;
  m_sel_number[1]++;

  //  Select good charged tracks in MDC ( option for PID )
  Vint icp, icm, iGood, jGood;
  icp.clear();
  icm.clear();
  iGood.clear();
  jGood.clear();

  map<int, int> ParticleType; //0:e, 1:mu, 2:pi, 3:K, 4:p
  ParticleID* pid = ParticleID::instance();

  for(unsigned int i = 0; i < recEvent->totalCharged(); i++) {
    jGood.push_back(0);
    EvtRecTrackIterator itTrk = recTrackCol->begin() + i;
    if(!(*itTrk)->isMdcTrackValid()) continue;
    if(!(*itTrk)->isMdcKalTrackValid()) continue;
    RecMdcTrack* mdcTrk = (*itTrk)->mdcTrack();
    if (fabs(cos(mdcTrk->theta())) >= m_cosThetaCut) continue;
    if (fabs(mdcTrk->z()) >= m_vz0Cut) continue;
    //iGood.push_back((*itTrk)->trackId());
    iGood.push_back(i);

    if (m_pid == 0) {
      if (mdcTrk->charge() > 0) {
        //icp.push_back((*itTrk)->trackId());
        icp.push_back(i);
      } else if (mdcTrk->charge() < 0) {
        //icm.push_back((*itTrk)->trackId());
        icm.push_back(i);
      }
    } else if (m_pid == 1) {
      ParticleType[i] = 2; // FIXME default value is Pion ?
      // pid info
      pid->init();
      pid->setChiMinCut(4);
      pid->setRecTrack(*itTrk);
      pid->setMethod(pid->methodProbability());
      //pid->setMethod(pid->methodLikelihood());    
      // use PID sub-system
      if(m_useDedx)pid->usePidSys(pid->useDedx());
      if(m_useTof1)pid->usePidSys(pid->useTof1());
      if(m_useTof2)pid->usePidSys(pid->useTof2());
      if(m_useTofE)pid->usePidSys(pid->useTofE());
      if(m_useTofQ)pid->usePidSys(pid->useTofQ());
      if(m_useEmc)pid->usePidSys(pid->useEmc());
      if(m_useMuc)pid->usePidSys(pid->useMuc());
      pid->identify(pid->onlyElectron());
      pid->identify(pid->onlyMuon());
      pid->identify(pid->onlyPion());
      pid->identify(pid->onlyKaon());
      pid->identify(pid->onlyProton());
      pid->calculate();
      if(!pid->IsPidInfoValid()) continue;
      double prob_value[5];
      prob_value[0] = pid->probElectron();
      prob_value[1] = pid->probMuon();
      prob_value[2] = pid->probPion();
      prob_value[3] = pid->probKaon();
      prob_value[4] = pid->probProton();

      int max = 0;
      for (int i = 1; i < 5; i++) {
        if (prob_value[i] > prob_value[max]) {
          max = i;
        }
      }
    
//cout << "PID : n = " << n << endl;
 
      if(max == 0) ParticleType[i] = 0; //m_sel_number[3]++;}
      if(max == 1) ParticleType[i] = 1; //m_sel_number[4]++;}
      if(max == 2) ParticleType[i] = 2; //m_sel_number[5]++;}
      if(max == 3) ParticleType[i] = 3;//m_sel_number[6]++;}
      if(max == 4) ParticleType[i] =4;//m_sel_number[7]++;}
    }

    // fill z position of track parameters 
    RecMdcKalTrack *kalTrk = (*itTrk)->mdcKalTrack();
    RecMdcKalTrack::setPidType(RecMdcKalTrack::pion);
    m_mdcTrk_x = mdcTrk->x();
    m_mdcTrk_y = mdcTrk->y();
    m_mdcTrk_z = mdcTrk->helix(3);
    m_mdcTrk_r = mdcTrk->r();
    m_mdcTrk_dr = mdcTrk->helix(0);
    m_mdcKalTrk_z = kalTrk->helix()[3];
    m_tuple6->write();
  }

  // Cut 2 : for anything sample, at least 3 good charged tracks 
  if (iGood.size() < m_trackNumberCut) return StatusCode::SUCCESS;
  m_sel_number[2]++;

  //   Vertex Finding
  if (m_vertexFind == 1) {
    int nGood = iGood.size();
    for(int i1 = 0; i1 < nGood; i1++) {
      int id_trk1 = iGood[i1];
      EvtRecTrackIterator trk1 = (recTrackCol->begin()+id_trk1);
      RecMdcKalTrack *kalTrk1 = (*trk1)->mdcKalTrack();
      //FIXME default set is ?
      HTrackParameter htrk1, htrk2;
      if (m_pid == 0) {
        RecMdcKalTrack::setPidType(RecMdcKalTrack::pion);
        htrk1 = HTrackParameter(kalTrk1->helix(), kalTrk1->err(), id_trk1, 2);
      } else if (m_pid == 1) {
        RecMdcKalTrack::setPidType((RecMdcKalTrack::PidType)ParticleType[id_trk1]);
        htrk1 = HTrackParameter(kalTrk1->helix(),kalTrk1->err(), id_trk1, ParticleType[id_trk1]);
      }
      for(int i2 = i1+1; i2 < nGood; i2++) {
        int id_trk2 = iGood[i2];
        EvtRecTrackIterator trk2 = (recTrackCol->begin()+id_trk2);
        RecMdcKalTrack *kalTrk2 = (*trk2)->mdcKalTrack();
        //FIXME default set is ?
        if (m_pid == 0) {
          RecMdcKalTrack::setPidType(RecMdcKalTrack::pion);
          htrk2 = HTrackParameter(kalTrk2->helix(), kalTrk2->err(), id_trk2, 2);
        } else if (m_pid == 1) {
          RecMdcKalTrack::setPidType((RecMdcKalTrack::PidType)ParticleType[id_trk2]);
          htrk2 = HTrackParameter(kalTrk2->helix(),kalTrk2->err(),id_trk2,ParticleType[id_trk2]);
        }
        HepPoint3D cross(0., 0., 0.);
        double dist = htrk1.minDistanceTwoHelix(htrk2, cross);
        m_mind = dist;
        m_xc = cross.x();
        m_yc = cross.y();
        m_zc = cross.z();
        m_tuple1->write();
        
         if(sqrt(cross.x()*cross.x()+cross.y()*cross.y())>2)continue;
        
        if(fabs(dist) > m_minDistance) continue;  // Cut condition 2
       // double cross_cut = (cross.x() - m_meanPointX) * (cross.x() - m_meanPointX)
       //                              /m_minPointX/m_minPointX
       //                   + (cross.y() - m_meanPointY) * (cross.y() - m_meanPointY)
       //                              /m_minPointY/m_minPointY;
       //                             //FIXME sigma value from Gaussian fit //0.071/0.059 
        double cross_cut=0.;
        if(cross_cut < 3.5 * 3.5) {   //Cut condition 3
          jGood[id_trk1] = 1;
          jGood[id_trk2] = 1;
        } //end if
      }
    }

    iGood.clear();
    for(int i =0; i < jGood.size(); i++) {
      if(jGood[i] == 1) iGood.push_back(i);
    }

  } //end if vertex finding

  if (iGood.size() >= 2) m_sel_number[3]++;
  if (iGood.size() >= 3) m_sel_number[4]++;

  //  Vertex Fit
  HepPoint3D vx(0., 0., 0.);
  HepSymMatrix Evx(3, 0);
  double bx = 1E+6; 
  double by = 1E+6;
  double bz = 1E+6;
  Evx[0][0] = bx*bx;
  Evx[1][1] = by*by;
  Evx[2][2] = bz*bz;
  VertexParameter vxpar;
  vxpar.setVx(vx);
  vxpar.setEvx(Evx);

  VertexFit* vtxfit = VertexFit::instance();
  // Step 1: using MdcKalTrk parameters
  vtxfit->init();
  Vint tlis;
  tlis.clear();
  for(int i = 0; i < iGood.size(); i++) {
    int trk_id = iGood[i];
    EvtRecTrackIterator itTrk = recTrackCol->begin()+trk_id;
    RecMdcKalTrack *kalTrk = (*itTrk)->mdcKalTrack();
    if (m_pid == 0) {
      RecMdcKalTrack::setPidType(RecMdcKalTrack::pion);
      WTrackParameter wtrk(xmass[2], kalTrk->helix(), kalTrk->err());
      vtxfit->AddTrack(i ,wtrk);
    } else if (m_pid == 1) {
      RecMdcKalTrack::setPidType((RecMdcKalTrack::PidType)ParticleType[trk_id]);
      WTrackParameter wtrk(xmass[ParticleType[trk_id]], kalTrk->helix(), kalTrk->err());
      vtxfit->AddTrack(i ,wtrk);
    }
    tlis.push_back(i);
  }
  //if(iGood.size() >= m_trackNumberCut) { //at least N tracks
    vtxfit->AddVertex(0, vxpar, tlis);
    if(!vtxfit->Fit(0)) return StatusCode::SUCCESS;
    m_chig = vtxfit->chisq(0);
    m_ndofg = 2 * iGood.size() - 3;
    m_probg = TMath::Prob(m_chig, m_ndofg);
    HepVector glvtx = vtxfit->Vx(0);
    m_gvx = glvtx[0];
    m_gvy = glvtx[1];
    m_gvz = glvtx[2];
    m_tuple4->write();
 
    if(!(m_vertex_x->Fill(glvtx[0]))){
            log << MSG::FATAL << "Couldn't Fill x of vertex" << endreq;
            exit(1);
    }
    if(!(m_vertex_y->Fill(glvtx[1]))){
            log << MSG::FATAL << "Couldn't Fill y of vertex" << endreq;
            exit(1);
    }
    if(!(m_vertex_z->Fill(glvtx[2]))){
          log << MSG::FATAL << "Couldn't Fill z of vertex" << endreq;
          exit(1);
    }
  //}
  for (int j = 0; j < iGood.size(); j++) {
    HepVector pull(5, 0);
    bool success = vtxfit->pull(0, j, pull);
    if (success) {
      m_gpull_drho = pull[0];
      m_gpull_phi = pull[1];
      m_gpull_kapha = pull[2];
      m_gpull_dz = pull[3];
      m_gpull_lamb = pull[4];
      m_tuple7->write();
    } 
  }
 
  // Step 2 : Kalman Vertex Fitting
  KalmanVertexFit *kalvtx = KalmanVertexFit::instance();
  kalvtx->init();
  kalvtx->initVertex(vxpar);
  kalvtx->setChisqCut(m_chisqCut);
  kalvtx->setTrackIteration(m_trackIteration);
  kalvtx->setVertexIteration(m_vertexIteration);
  kalvtx->setTrackIterationCut(m_chi2CutforTrkIter);
  for(int i = 0; i < iGood.size(); i++) {
    int trk_id = iGood[i];
    EvtRecTrackIterator itTrk = recTrackCol->begin()+trk_id;
    RecMdcKalTrack *kalTrk = (*itTrk)->mdcKalTrack();
    if (m_pid == 0) {
      RecMdcKalTrack::setPidType(RecMdcKalTrack::pion);
      HTrackParameter htrk(kalTrk->helix(), kalTrk->err(), trk_id, 2);
      kalvtx->addTrack(htrk);
    } else if (m_pid == 1) {
      RecMdcKalTrack::setPidType((RecMdcKalTrack::PidType)ParticleType[trk_id]);
      HTrackParameter htrk(kalTrk->helix(), kalTrk->err(), trk_id, ParticleType[trk_id]);
      kalvtx->addTrack(htrk);
    }
  }
  kalvtx->filter();
  //int freedomCut = -3 + 2*m_trackNumberCut;
  int freedomCut = 1;
  if(kalvtx->ndof() >= freedomCut) {  //Cut condition 5 //add 2 when add every track
    //for(int i = 0; i < kalvtx->numTrack(); i++) {
    for(int i = 0; i < iGood.size(); i++) {
      m_chif = kalvtx->chiF(i);
      m_chis = kalvtx->chiS(i);
      m_probs = TMath::Prob(m_chis, 2);
      m_probf = TMath::Prob(m_chif, 2);
      m_tuple2->write();

      if(kalvtx->chiS(i) > m_chi2CutforSmooth) kalvtx->remove(i);
    }
  }
  if(kalvtx->ndof() >= freedomCut) { //FIXME  to Rhopi is 0 , others are 1
    for(int i = 0; i < iGood.size(); i++) {
      kalvtx->smooth(i);

      HepVector Pull(5, 0);
      Pull = kalvtx->pull(i);
      m_pull_drho = Pull[0];
      m_pull_phi = Pull[1];
      m_pull_kapha = Pull[2];
      m_pull_dz = Pull[3];
      m_pull_lamb = Pull[4];
      m_pull_momentum = kalvtx->pullmomentum(i);
      m_tuple5->write();
    }

    m_chik = kalvtx->chisq();
    m_ndofk = kalvtx->ndof();
    m_probk = TMath::Prob(m_chik, m_ndofk);
    HepVector xp(3, 0);
    xp = kalvtx->x();
    m_kvx = xp[0];
    m_kvy = xp[1];
    m_kvz = xp[2];
    m_tuple3->write();

    if(!(m_vertex_x_kal->Fill(xp[0]))){
            log << MSG::FATAL << "Couldn't Fill kal x of vertex" << endreq;
            exit(1);
    }
    if(!(m_vertex_y_kal->Fill(xp[1]))){
            log << MSG::FATAL << "Couldn't Fill kal y of vertex" << endreq;
            exit(1);
    }
    if(!(m_vertex_z_kal->Fill(xp[2]))){
            log << MSG::FATAL << "Couldn't Fill kal z of vertex" << endreq;
            exit(1);
    }

    m_sel_number[3]++;
  }
  return StatusCode::SUCCESS;
}

StatusCode BeamParams::finalize

(

)

Definition at line 603 of file BeamParams.cxx.

References genRecEmupikp::i, Bes_Common::INFO, m_sel_number, and msgSvc().

{
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in finalize()" << endmsg;

  /*TF1 *func = new TF1("func", "gaus", -0.6, 0.6);
  m_vertex_x->Fit("func", "NR+");
  Double_t MeanX = func->GetParameter(1);
  Double_t SigmaX = func->GetParameter(2);
  Double_t ErrMeanX = func->GetParError(1);
  Double_t ErrSigmaX = func->GetParError(2);

  TF1 *funcY = new TF1("funcY", "gaus", -0.6, 0.1);
  m_vertex_y->Fit("funcY", "NR+");
  Double_t MeanY = funcY->GetParameter(1);
  Double_t SigmaY = funcY->GetParameter(2);
  Double_t ErrMeanY = funcY->GetParError(1);
  Double_t ErrSigmaY = funcY->GetParError(2);

  TF1 *funcZ = new TF1("funcZ", "gaus", -6, 6);
  m_vertex_z->Fit("funcZ", "NR+");
  Double_t MeanZ = funcZ->GetParameter(1);
  Double_t SigmaZ = funcZ->GetParameter(2);
  Double_t ErrMeanZ = funcZ->GetParError(1);
  Double_t ErrSigmaZ = funcZ->GetParError(2);

  TCanvas *myC = new TCanvas("myC", "myC", 1200, 400);
  TPad *background = (TPad*)gPad;
  background->SetFillColor(10);
  myC->Divide(3,1);
  myC->cd(1);

  m_vertex_x_kal->Fit("func", "R+");
  Double_t MeanXKal = func->GetParameter(1);
  Double_t SigmaXKal = func->GetParameter(2);
  Double_t ErrMeanXKal = func->GetParError(1);
  Double_t ErrSigmaXKal = func->GetParError(2);

  myC->cd(2);
  m_vertex_y_kal->Fit("funcY", "R+");
  Double_t MeanYKal = funcY->GetParameter(1);
  Double_t SigmaYKal = funcY->GetParameter(2);
  Double_t ErrMeanYKal = funcY->GetParError(1);
  Double_t ErrSigmaYKal = funcY->GetParError(2);

  myC->cd(3);
  m_vertex_z_kal->Fit("funcZ", "R+");
  Double_t MeanZKal = funcZ->GetParameter(1);
  Double_t SigmaZKal = funcZ->GetParameter(2); 
  Double_t ErrMeanZKal = funcZ->GetParError(1);
  Double_t ErrSigmaZKal = funcZ->GetParError(2);
  
  cout << "Kal: Mean X, Y, Z = "<<MeanXKal << "  " << MeanYKal << "  " << MeanZKal << endl;
  cout << "Kal: Sigma X, Y, Z = "<<SigmaXKal<<"  " <<SigmaYKal <<"  " << SigmaZKal << endl;
  cout << "Gl: Mean X, Y, Z = " << MeanX << "  " << MeanY << "  " << MeanZ << endl;
  cout << "Gl: Sigma X, Y, Z = " << SigmaX << "  " << SigmaY << "  " << SigmaZ << endl;

  // Output a TXT file and a .ps figure for storing the fit results  
  const char *file_name, *figs_name;
  if (m_hadronFile == 0) {
    file_name = m_fileNameDst.c_str();
  } else if (m_hadronFile == 1) {
    file_name = m_fileNameHadron.c_str();
  }

  figs_name = m_figsName.c_str();
  myC->SaveAs(figs_name);

  ofstream file(file_name);
 
  file << getenv("BES_RELEASE") << " " << m_parVer << " " << m_runNo << endl;  
  file << MeanX << " " << MeanY << " " << MeanZ << " "<< SigmaX << " "<< SigmaY << " " << SigmaZ << endl;
  file << ErrMeanX << " " << ErrMeanY<< " " << ErrMeanZ << " " << ErrSigmaX << " " << ErrSigmaY << " " << ErrSigmaZ << endl;
  file << MeanXKal << " "<< MeanYKal << " "<< MeanZKal << " "<< SigmaXKal << " " << SigmaYKal << " " << SigmaZKal << endl;
  file << ErrMeanXKal << " " << ErrMeanYKal<< " " << ErrMeanZKal << " " << ErrSigmaXKal << " " << ErrSigmaYKal << " " << ErrSigmaZKal << endl;*/

  /*delete m_vertex_x;
  delete m_vertex_y;
  delete m_vertex_z;
  delete m_vertex_x_kal;
  delete m_vertex_y_kal;
  delete m_vertex_z_kal;
  */
  log << MSG::ALWAYS << "==================================" << endreq;
  log << MSG::ALWAYS << "survived event :";
  for(int i = 0; i < 10; i++){
    log << MSG::ALWAYS << m_sel_number[i] << " ";
  }
  log << MSG::ALWAYS << endreq;
  log << MSG::ALWAYS << "==================================" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode BeamParams::initialize

(

)

Definition at line 88 of file BeamParams.cxx.

References Bes_Common::FATAL, genRecEmupikp::i, Bes_Common::INFO, m_chif, m_chig, m_chik, m_chis, m_gpull_drho, m_gpull_dz, m_gpull_kapha, m_gpull_lamb, m_gpull_phi, m_gvx, m_gvy, m_gvz, m_kvx, m_kvy, m_kvz, m_mdcKalTrk_z, m_mdcTrk_dr, m_mdcTrk_r, m_mdcTrk_x, m_mdcTrk_y, m_mdcTrk_z, m_mind, m_ndofg, m_ndofk, m_probf, m_probg, m_probk, m_probs, m_pull_drho, m_pull_dz, m_pull_kapha, m_pull_lamb, m_pull_momentum, m_pull_phi, m_rxy, m_sel_number, m_thsvc, m_tuple1, m_tuple2, m_tuple3, m_tuple4, m_tuple5, m_tuple6, m_tuple7, m_vertex_x, m_vertex_x_kal, m_vertex_y, m_vertex_y_kal, m_vertex_z, m_vertex_z_kal, m_xc, m_yc, m_zc, msgSvc(), and ntupleSvc().

{
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in initialize()" << endmsg;
   
  StatusCode sc=service("THistSvc", m_thsvc);
  if(sc.isFailure()) {
          log << MSG::FATAL << "Couldn't get THistSvc" << endreq;
          exit(1);
  }

  for(int i = 0; i < 15; i++){
    m_sel_number[i] = 0;
  }

  //   Book histogram
  m_vertex_x = new TH1D("x_of_vertex", "x of vertex", 200, -5, 5);
  m_vertex_y = new TH1D("y_of_vertex", "y of vertex", 200, -5, 5);
  m_vertex_z = new TH1D("z_of_vertex", "z of vertex", 200, -10, 10);
  m_vertex_x_kal = new TH1D("x_of_vertex_in_kal", "x of vertex in kal", 200, -5, 5);
  m_vertex_y_kal = new TH1D("y_of_vertex_in_kal", "y of vertex in kal", 200, -5, 5);
  m_vertex_z_kal = new TH1D("z_of_vertex_in_kal", "z of vertex in kal", 200, -10, 10);
  if(m_thsvc->regHist("/DQAHist/zhsVER/x_of_vertex",m_vertex_x).isFailure()){
          log << MSG::FATAL << "Couldn't register x of vertex " << endreq;
          exit(1);
  }
  if(m_thsvc->regHist("/DQAHist/zhsVER/y_of_vertex",m_vertex_y).isFailure()){
          log << MSG::FATAL << "Couldn't register y of vertex" << endreq;
          exit(1);
  }
  if(m_thsvc->regHist("/DQAHist/zhsVER/z_of_vertex",m_vertex_z).isFailure()){
          log << MSG::FATAL << "Couldn't register z of vertex" << endreq;
          exit(1);
  }
  if(m_thsvc->regHist("/DQAHist/zhsVER/x_of_vertex_in_kal",m_vertex_x_kal).isFailure()){
          log << MSG::FATAL << "Couldn't register x of vertex in kal" << endreq;
          exit(1);
  }
  if(m_thsvc->regHist("/DQAHist/zhsVER/y_of_vertex_in_kal",m_vertex_y_kal).isFailure()){
          log << MSG::FATAL << "Couldn't register y of vertex in kal" << endreq;
          exit(1);
  }
  if(m_thsvc->regHist("/DQAHist/zhsVER/z_of_vertex_in_kal",m_vertex_z_kal).isFailure()){
          log << MSG::FATAL << "Couldn't register z of vertex in kal" << endreq;
          exit(1);
  }
  // end book

  StatusCode status;

  NTuplePtr nt1(ntupleSvc(), "FILE1/minid");// minimal distance 
  if(nt1) m_tuple1 = nt1;
  else {
    m_tuple1 = ntupleSvc()->book ("FILE1/minid", CLID_ColumnWiseTuple, "minimal distance");
    if(m_tuple1) {
      status = m_tuple1->addItem("xc", m_xc);
      status = m_tuple1->addItem("yc", m_yc);
      status = m_tuple1->addItem("zc", m_zc);
      status = m_tuple1->addItem("mind", m_mind);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple1) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  NTuplePtr nt2(ntupleSvc(), "FILE1/chisq"); //chisq of Kalman filter 
  if(nt2) m_tuple2 = nt2;
  else {
    m_tuple2 = ntupleSvc()->book ("FILE1/chisq", CLID_ColumnWiseTuple, "chi-square of ");
    if(m_tuple2) {
      status = m_tuple2->addItem("chis", m_chis);
      status = m_tuple2->addItem("probs", m_probs);
      status = m_tuple2->addItem("chif", m_chif);
      status = m_tuple2->addItem("probf", m_probf);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple2) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  NTuplePtr nt3(ntupleSvc(), "FILE1/kalvtx"); 
  if(nt3) m_tuple3 = nt3;
  else {
    m_tuple3 = ntupleSvc()->book ("FILE1/kalvtx", CLID_ColumnWiseTuple, "kalman vertex");
    if(m_tuple3) {
      status = m_tuple3->addItem("kvx", m_kvx);
      status = m_tuple3->addItem("kvy", m_kvy);
      status = m_tuple3->addItem("kvz", m_kvz);
      status = m_tuple3->addItem("chik", m_chik);
      status = m_tuple3->addItem("ndofk", m_ndofk);
      status = m_tuple3->addItem("probk", m_probk);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple3) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  NTuplePtr nt4(ntupleSvc(), "FILE1/glbvtx");
  if(nt4) m_tuple4 = nt4;
  else {
    m_tuple4 = ntupleSvc()->book ("FILE1/glbvtx", CLID_ColumnWiseTuple, "global vertex");
    if(m_tuple4) {
      status = m_tuple4->addItem("gvx", m_gvx);
      status = m_tuple4->addItem("gvy", m_gvy);
      status = m_tuple4->addItem("gvz", m_gvz);
      status = m_tuple4->addItem("chig", m_chig);
      status = m_tuple4->addItem("ndofg", m_ndofg);
      status = m_tuple4->addItem("probg", m_probg);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple4) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  NTuplePtr nt5(ntupleSvc(), "FILE1/Pull");
  if(nt5) m_tuple5 = nt5;
  else {
    m_tuple5 = ntupleSvc()->book ("FILE1/Pull", CLID_ColumnWiseTuple, "Pull");
    if(m_tuple5) {
      status = m_tuple5->addItem("pull_drho", m_pull_drho);
      status = m_tuple5->addItem("pull_phi", m_pull_phi);
      status = m_tuple5->addItem("pull_kapha", m_pull_kapha);
      status = m_tuple5->addItem("pull_dz", m_pull_dz);
      status = m_tuple5->addItem("pull_lamb", m_pull_lamb);
      status = m_tuple5->addItem("pull_momentum", m_pull_momentum);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple5) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  NTuplePtr nt6(ntupleSvc(), "FILE1/MdcTrack");
  if(nt6) m_tuple6 = nt6;
  else {
    m_tuple6 = ntupleSvc()->book ("FILE1/MdcTrack", CLID_ColumnWiseTuple, "MdcTrack");
    if(m_tuple6) {
      status = m_tuple6->addItem("MdcTrkX", m_mdcTrk_x);
      status = m_tuple6->addItem("MdcTrkY", m_mdcTrk_y);
      status = m_tuple6->addItem("MdcTrkZ", m_mdcTrk_z);
      status = m_tuple6->addItem("MdcTrkR", m_mdcTrk_r);
      status = m_tuple6->addItem("MdcTrkDrho", m_mdcTrk_dr);
      status = m_tuple6->addItem("Rxy", m_rxy);
      status = m_tuple6->addItem("MdcKalTrkZ", m_mdcKalTrk_z);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple6) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  NTuplePtr nt7(ntupleSvc(), "FILE1/PullG");
  if(nt7) m_tuple7 = nt7;
  else {
    m_tuple7 = ntupleSvc()->book ("FILE1/PullG", CLID_ColumnWiseTuple, "Pull");
    if(m_tuple7) {
      status = m_tuple7->addItem("gpull_drho", m_gpull_drho);
      status = m_tuple7->addItem("gpull_phi", m_gpull_phi);
      status = m_tuple7->addItem("gpull_kapha", m_gpull_kapha);
      status = m_tuple7->addItem("gpull_dz", m_gpull_dz);
      status = m_tuple7->addItem("gpull_lamb", m_gpull_lamb);
    } else {
      log << MSG::FATAL << "Cannot book N-tuple:" << long(m_tuple7) << endmsg;
      return StatusCode::FAILURE;
    }
  }

  log << MSG::INFO << "successfully return from initialize()" <<endmsg;
  return StatusCode::SUCCESS;
}

---

Member Data Documentation

double BeamParams::m_chi2CutforSmooth [private]

Definition at line 38 of file BeamParams.h.

Referenced by BeamParams(), and execute().

double BeamParams::m_chi2CutforTrkIter [private]

Definition at line 37 of file BeamParams.h.

Referenced by BeamParams(), and execute().

NTuple::Item<double> BeamParams::m_chif [private]

Definition at line 84 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_chig [private]

Definition at line 98 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_chik [private]

Definition at line 89 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_chis [private]

Definition at line 83 of file BeamParams.h.

Referenced by execute(), and initialize().

double BeamParams::m_chisqCut [private]

Definition at line 34 of file BeamParams.h.

Referenced by BeamParams(), and execute().

double BeamParams::m_cosThetaCut [private]

Definition at line 28 of file BeamParams.h.

Referenced by BeamParams(), and execute().

double BeamParams::m_cpu[5] [private]

Definition at line 70 of file BeamParams.h.

HepVector BeamParams::m_Cpuvc [private]

Definition at line 69 of file BeamParams.h.

std::string BeamParams::m_figsName [private]

Definition at line 60 of file BeamParams.h.

Referenced by BeamParams().

std::string BeamParams::m_fileNameDst [private]

Definition at line 58 of file BeamParams.h.

Referenced by BeamParams().

std::string BeamParams::m_fileNameHadron [private]

Definition at line 59 of file BeamParams.h.

Referenced by BeamParams().

NTuple::Item<double> BeamParams::m_gpull_drho [private]

Definition at line 123 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gpull_dz [private]

Definition at line 126 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gpull_kapha [private]

Definition at line 125 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gpull_lamb [private]

Definition at line 127 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gpull_phi [private]

Definition at line 124 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gvx [private]

Definition at line 101 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gvy [private]

Definition at line 102 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_gvz [private]

Definition at line 103 of file BeamParams.h.

Referenced by execute(), and initialize().

int BeamParams::m_hadronFile [private]

Definition at line 72 of file BeamParams.h.

Referenced by BeamParams().

NTuple::Item<double> BeamParams::m_kvx [private]

Definition at line 92 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_kvy [private]

Definition at line 93 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_kvz [private]

Definition at line 94 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_mdcKalTrk_z [private]

Definition at line 120 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_mdcTrk_dr [private]

Definition at line 119 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_mdcTrk_r [private]

Definition at line 117 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_mdcTrk_x [private]

Definition at line 114 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_mdcTrk_y [private]

Definition at line 115 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_mdcTrk_z [private]

Definition at line 116 of file BeamParams.h.

Referenced by execute(), and initialize().

double BeamParams::m_meanPointX [private]

Definition at line 32 of file BeamParams.h.

Referenced by BeamParams().

double BeamParams::m_meanPointY [private]

Definition at line 33 of file BeamParams.h.

Referenced by BeamParams().

NTuple::Item<double> BeamParams::m_mind [private]

Definition at line 80 of file BeamParams.h.

Referenced by execute(), and initialize().

double BeamParams::m_minDistance [private]

Definition at line 29 of file BeamParams.h.

Referenced by BeamParams(), and execute().

double BeamParams::m_minPointX [private]

Definition at line 30 of file BeamParams.h.

Referenced by BeamParams().

double BeamParams::m_minPointY [private]

Definition at line 31 of file BeamParams.h.

Referenced by BeamParams().

NTuple::Item<long> BeamParams::m_ndofg [private]

Definition at line 99 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<long> BeamParams::m_ndofk [private]

Definition at line 90 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<long> BeamParams::m_numTrk [private]

Definition at line 95 of file BeamParams.h.

int BeamParams::m_parVer [private]

Definition at line 57 of file BeamParams.h.

Referenced by BeamParams().

int BeamParams::m_pid [private]

Definition at line 24 of file BeamParams.h.

Referenced by BeamParams(), and execute().

double BeamParams::m_PidProbCut [private]

Definition at line 49 of file BeamParams.h.

double BeamParams::m_pivotX [private]

Definition at line 51 of file BeamParams.h.

double BeamParams::m_pivotY [private]

Definition at line 52 of file BeamParams.h.

double BeamParams::m_pivotZ [private]

Definition at line 53 of file BeamParams.h.

NTuple::Item<double> BeamParams::m_probf [private]

Definition at line 86 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_probg [private]

Definition at line 100 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_probk [private]

Definition at line 91 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_probs [private]

Definition at line 85 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_pull_drho [private]

Definition at line 106 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_pull_dz [private]

Definition at line 109 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_pull_kapha [private]

Definition at line 108 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_pull_lamb [private]

Definition at line 110 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_pull_momentum [private]

Definition at line 111 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_pull_phi [private]

Definition at line 107 of file BeamParams.h.

Referenced by execute(), and initialize().

int BeamParams::m_runNo [private]

Definition at line 56 of file BeamParams.h.

Referenced by execute().

NTuple::Item<double> BeamParams::m_rxy [private]

Definition at line 118 of file BeamParams.h.

Referenced by initialize().

int BeamParams::m_sel_number[15] [private]

Definition at line 22 of file BeamParams.h.

Referenced by execute(), finalize(), and initialize().

ITHistSvc* BeamParams::m_thsvc [private]

Definition at line 21 of file BeamParams.h.

Referenced by initialize().

int BeamParams::m_trackIteration [private]

Definition at line 35 of file BeamParams.h.

Referenced by BeamParams(), and execute().

int BeamParams::m_trackNumberCut [private]

Definition at line 39 of file BeamParams.h.

Referenced by BeamParams(), and execute().

int BeamParams::m_trkNum [private]

Definition at line 55 of file BeamParams.h.

NTuple::Tuple* BeamParams::m_tuple1 [private]

Definition at line 76 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Tuple* BeamParams::m_tuple2 [private]

Definition at line 82 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Tuple* BeamParams::m_tuple3 [private]

Definition at line 88 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Tuple* BeamParams::m_tuple4 [private]

Definition at line 97 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Tuple* BeamParams::m_tuple5 [private]

Definition at line 105 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Tuple* BeamParams::m_tuple6 [private]

Definition at line 113 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Tuple* BeamParams::m_tuple7 [private]

Definition at line 122 of file BeamParams.h.

Referenced by execute(), and initialize().

bool BeamParams::m_useDedx [private]

Definition at line 42 of file BeamParams.h.

Referenced by BeamParams(), and execute().

bool BeamParams::m_useEmc [private]

Definition at line 47 of file BeamParams.h.

Referenced by BeamParams(), and execute().

bool BeamParams::m_useMuc [private]

Definition at line 48 of file BeamParams.h.

Referenced by BeamParams(), and execute().

bool BeamParams::m_useTof1 [private]

Definition at line 43 of file BeamParams.h.

Referenced by BeamParams(), and execute().

bool BeamParams::m_useTof2 [private]

Definition at line 44 of file BeamParams.h.

Referenced by BeamParams(), and execute().

bool BeamParams::m_useTofE [private]

Definition at line 45 of file BeamParams.h.

Referenced by BeamParams(), and execute().

bool BeamParams::m_useTofQ [private]

Definition at line 46 of file BeamParams.h.

Referenced by BeamParams(), and execute().

TH1D* BeamParams::m_vertex_x [private]

Definition at line 62 of file BeamParams.h.

Referenced by execute(), and initialize().

TH1D* BeamParams::m_vertex_x_kal [private]

Definition at line 65 of file BeamParams.h.

Referenced by execute(), and initialize().

TH1D* BeamParams::m_vertex_y [private]

Definition at line 63 of file BeamParams.h.

Referenced by execute(), and initialize().

TH1D* BeamParams::m_vertex_y_kal [private]

Definition at line 66 of file BeamParams.h.

Referenced by execute(), and initialize().

TH1D* BeamParams::m_vertex_z [private]

Definition at line 64 of file BeamParams.h.

Referenced by execute(), and initialize().

TH1D* BeamParams::m_vertex_z_kal [private]

Definition at line 67 of file BeamParams.h.

Referenced by execute(), and initialize().

int BeamParams::m_vertexFind [private]

Definition at line 25 of file BeamParams.h.

Referenced by BeamParams(), and execute().

int BeamParams::m_vertexIteration [private]

Definition at line 36 of file BeamParams.h.

Referenced by BeamParams(), and execute().

double BeamParams::m_vz0Cut [private]

Definition at line 27 of file BeamParams.h.

Referenced by BeamParams(), and execute().

NTuple::Item<double> BeamParams::m_xc [private]

Definition at line 77 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_yc [private]

Definition at line 78 of file BeamParams.h.

Referenced by execute(), and initialize().

NTuple::Item<double> BeamParams::m_zc [private]

Definition at line 79 of file BeamParams.h.

Referenced by execute(), and initialize().

---