RecCgemCluster.cxx

Go to the documentation of this file.

#include "CgemRecEvent/RecCgemCluster.h"
 
RecCgemCluster::RecCgemCluster():
 
  m_clusterid(-1),
  m_trkid(-1),
  m_layerid(-1), //cluster in which CgemLayer//
  m_sheetid(-1), //cluster in which sheet//
  m_flag(-1),    //v cluster or x cluster or mixed cluster//
  m_energydeposit(-99.0), //energy deposit//
  m_recphi(-99.0),   //position for reconstruction//
  m_recv(-99.0),
  m_recZ(-99.0),
  m_recphi_CC(-99.0),   //position for reconstruction with CC//
  m_recv_CC(-99.0),
  m_recZ_CC(-99.0),
  m_recphi_mTPC(-99.0),   //position for reconstruction with uTPC// 
  m_recv_mTPC(-99.0),
  m_recZ_mTPC(-99.0),
  m_slope_mTPC(-9999.0),
  m_inter_mTPC(-9999.0),
  m_clusterflag(-1,-1){
}
 
 
 
RecCgemCluster::RecCgemCluster(const RecCgemCluster& reccluster):
  m_clusterid(reccluster.m_clusterid),
  m_trkid(reccluster.m_trkid),
  m_layerid(reccluster.m_layerid), //cluster in which CgemLayer//
  m_sheetid(reccluster.m_sheetid), //cluster in which sheet//
  m_flag(reccluster.m_flag),    //v cluster or x cluster or mixed cluster//
  m_energydeposit(reccluster.m_energydeposit), //energy deposit//
  m_recphi(reccluster.m_recphi),   //position for reconstruction//
  m_recv(reccluster.m_recv),
  m_recphi_CC(reccluster.m_recphi_CC),   //position for reconstruction with CC// 
  m_recv_CC(reccluster.m_recv_CC),
  m_recZ_CC(reccluster.m_recZ_CC),
  m_recphi_mTPC(reccluster.m_recphi_mTPC),   //position for reconstruction with uTPC// 
  m_recv_mTPC(reccluster.m_recv_mTPC),
  m_recZ_mTPC(reccluster.m_recZ_mTPC),
  m_slope_mTPC(reccluster.m_slope_mTPC),
  m_inter_mTPC(reccluster.m_inter_mTPC),
  m_clusterflag(reccluster.m_clusterflag){
} 
 
RecCgemCluster::~RecCgemCluster(){
}
 
double RecCgemCluster::get_merge_phi(){
    int size = getclusterflage() - getclusterflagb() + 1;
    double CC = getrecphi_CC();
    double TPC = getrecphi_mTPC();
    double delta = CC-TPC;
    double merge = CC;
    if(size>2 && abs(delta)<0.1) merge = (0.5 * CC + 0.5 * TPC);
    return merge;
}
 
double RecCgemCluster::get_merge_v(){
    int size = getclusterflage() - getclusterflagb() + 1;
    double CC = getrecv_CC();
    double TPC = getrecv_mTPC();
    double delta = CC-TPC;
    double merge = CC;
    if(size>2 && abs(delta)<0.1) merge = (0.5 * CC + 0.5 * TPC);
    return merge;
}