MdcROOTGeo Class Reference

#include <MdcROOTGeo.h>

Inheritance diagram for MdcROOTGeo:

Public Member Functions
	MdcROOTGeo ()
	Constructor.
	~MdcROOTGeo ()
	Destructor.
void	InitFromGdml (const char *gdmlFile, const char *setupName)
	Initialize the instance of ROOTGeo.
void	SetNode ()
	Set the pointers to theirs nodes;.
int	CorrectReplica (int layer, int replica)
	Correct some axial layer id to copyNo;.
void	SetDefaultVis ()
	Set default visual attributes;.
void	SetAllVisible ()
	Set all visible;.
void	SetQuarterVisible ()
	Set quater visible;.
void	SetPhysicalNode ()
	Set the pointers to the physical nodes;.
int	GetSegmentNb ()
	Get number of segments;.
int	GetLayerNb ()
	Get number of layers;.
int	GetReplicaNb (int layer)
	Get number of replicas on each layer;.
TGeoVolume *	GetVolumeMdc ()
	Get Mdc volume;.
TGeoVolume *	GetVolumeSegment (int segment)
	Get segment volume;.
TGeoVolume *	GetVolumeLayer (int layer)
	Get layer volume;
TGeoVolume *	GetVolumeReplica (int layer)
	Get replica volume;.
TGeoNode *	GetSegment (int segment, int no)
	Get segment node;.
TGeoNode *	GetLayer (int layer)
	Get layer node;
TGeoNode *	GetReplica (int layer, int replica)
	Get replica node;.
TGeoPhysicalNode *	GetPhysicalReplica (int layer, int replica)
	Get replica physical node;.
	MdcROOTGeo (Bool_t cgem_flag=0)
	Constructor.
	~MdcROOTGeo ()
	Destructor.
void	InitFromGDML (const char *gdmlFile, const char *setupName)
	Initialize ROOTGeo from GDML.
void	InitFromROOT (TGeoVolume *vol)
	Initialize ROOTGeo from TGeoVolume logicalMdc.
void	Init2DGeometry ()
	Initialize 2D Geometry.
void	SetNode ()
	Set the pointers to theirs nodes;.
int	CorrectReplica (int layer, int replica)
	Correct some axial layer id to copyNo;.
void	SetVolumeDefaultVis ()
	Set default visual attributes;.
void	SetAllVisible ()
	Set all visible;.
void	SetQuarterVisible ()
	Set quater visible;.
void	SetHalfVisible ()
	Set half visible;.
void	SetPhysicalNode ()
	Set the pointers to the physical nodes;.
void	SetPhysicalDefaultVis ()
	Set default physical node attributes;.
void	SetDetector ()
	Draw Detecor (what is detector depends on you)
void	ClearHits ()
	Set all physicalNodes corresponding to digiCol;.
void	SetHits ()
void	SetTFire (Bool_t input)
void	SetQFire (Bool_t input)
void	SetQNotOverflow (Bool_t input)
void	SetColorfulWire (Bool_t input)
void	SetMdcTimeSubEvTime (Bool_t input)
void	SetVisMdcDetector ()
	Set Mdc default detector visibility;.
void	SetVisMdcHits ()
	Set Mdc hits visibility;.
int	GetSegmentNb ()
	Get number of segments;.
int	GetLayerNb ()
	Get number of layers;.
int	GetReplicaNb (int layer)
	Get number of replicas on each layer;.
void	SetVolumeMdc (TGeoVolume *vol)
	Set Mdc volume, while initializing from ROOT;.
TGeoVolume *	GetVolumeMdc ()
	Get Mdc volume;.
TGeoVolume *	GetVolumeSegment (int segment)
	Get segment volume;.
TGeoVolume *	GetVolumeLayer (int layer)
	Get layer volume;
TGeoVolume *	GetVolumeReplica (int layer)
	Get replica volume;.
TGeoNode *	GetSegment (int segment, int no)
	Get segment node;.
TGeoNode *	GetLayer (int layer)
	Get layer node;
TGeoNode *	GetReplica (int layer, int replica)
	Get replica node;.
TGeoPhysicalNode *	GetPhysicalSegment (int segment)
	Get segment physical node;.
TGeoPhysicalNode *	GetPhysicalReplica (int layer, int replica)
	Get replica physical node;.
Mdc2DWire *	Get2DWire (Int_t layer, Int_t replica)
	Get Mdc2DWire;.
Int_t	GetSimuLayer (Int_t trueLayer)
	Trans trueLayer to simuLayer (0~42)->(0~49)
Bool_t	IsHit (Int_t layer, Int_t wire)
	Judge whether the digiCol contains (layer, wire)
void	Draw (Option_t *option)
	Draw function.
void	DrawHits (Option_t *option)
	Draw 2D hits.
Public Member Functions inherited from SubDetectorROOTGeo
	SubDetectorROOTGeo ()
	Constructor.
virtual	~SubDetectorROOTGeo ()
	Destructor.
void	ReadGdml (const char *gdmlFile, const char *setupName)
	Initialize the instance of ROOTGeo.
int	IsInitialized ()
	If the ROOT geometry of this subdetctor is initialized;.
void	SetChildNo (int childNo)
TGeoVolume *	GetLogicalVolume (const std::string &vn)
	Get a logical volume by name;.
TGeoVolume *	GetTopVolume ()
	Get the top(world) volume;.
TGeoVolumeAssembly *	GetAssemblyVolume (const std::string &an)
	Get an assembly by name;.
TGeoNode *	GetNode (const std::string &nn)
	Get a node(physical volume) by name;.
	SubDetectorROOTGeo ()
	Constructor.
virtual	~SubDetectorROOTGeo ()
	Destructor.
virtual void	ReadGdml (const char *gdmlFile, const char *setupName)
	Initialize the instance of ROOTGeo.
virtual int	IsInitialized ()
	If the ROOT geometry of this subdetctor is initialized;.
virtual void	SetChildNo (int childNo)
virtual TGeoVolume *	GetLogicalVolume (const std::string &vn)
	Get a logical volume by name;.
virtual TGeoVolume *	GetTopVolume ()
	Get the top(world) volume;.
virtual TGeoVolumeAssembly *	GetAssemblyVolume (const std::string &an)
	Get an assembly by name;.
virtual TGeoNode *	GetNode (const std::string &nn)
	Get a node(physical volume) by name;.
virtual int	GetHitsNum ()
	Get number of hits in HitsArray;.
virtual TGeoPhysicalNode *	GetHit (int i)
	Get ith hit in HitsArray;.
virtual void	SetDetectorOn ()
	Set all physicalNodes in m_DeteorsArray visible;.

Additional Inherited Members
Protected Attributes inherited from SubDetectorROOTGeo
int	m_ROOTGeoInit
int	m_childNo
TGeoVolume *	m_TopVolume
int	m_2DGeoInit
TObjArray *	m_DetectorsArray
TObjArray *	m_HitsArray
TObjArray *	m_2DHitsArray

Detailed Description

Class MdcGeo contains all of the objects necessary to describe the mdc geometry.

Author

Zhengyun You \URL{youzy.nosp@m.@hep.nosp@m..pku..nosp@m.cn}

Definition at line 29 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/ROOTGeo/MdcROOTGeo.h.

Constructor & Destructor Documentation

MdcROOTGeo() [1/2]

MdcROOTGeo::MdcROOTGeo

(

)

Constructor.

Definition at line 36 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{ 
  // Default constructor.
  for (int layer = 0; layer < m_kLayer; layer++) {
    for (int replica = 0; replica < m_kReplicaMax; replica++) {
      m_NodeReplica[layer][replica] = 0;
      m_PhysicalReplica[layer][replica] = 0;
    }
  }
 
  for (int i = 0; i < m_kCorrectLayer; i++) {
    m_CorrectMap[m_kiCorrectLayer[i]] = m_kiCorrectReplica[i];
  }
}

~MdcROOTGeo() [1/2]

MdcROOTGeo::~MdcROOTGeo

(

)

Destructor.

Definition at line 51 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{ }

MdcROOTGeo() [2/2]

MdcROOTGeo::MdcROOTGeo

(

Bool_t

cgem_flag = 0

)

Constructor.

Definition at line 51 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

        : SubDetectorROOTGeo()
{
    m_MdcColor         = 3;
    m_segmentColor     = 0;
    m_hypeColor        = 3;
    m_tubeColor        = 4;
    //m_twistedTubsColor = kBlue;
    m_twistedTubsColor = 63;  //Long Peixun's update: Blue -> Green, as same as XY View
    //m_replicaColor     = kCyan;
    m_replicaColor     = 73;  //Long Peixun's update: Grey -> Blue, as same as XY View
 
    //Long Peixun's update: According to CGEM flag, set start wire layer
    if (cgem_flag) m_StartLayer = 8;
    else m_StartLayer = 0;
 
    k_TFire = kTRUE;
    k_QFire = kTRUE;
    k_QNotOverflow = kFALSE;
 
    // Default constructor.
    for (int layer = 0; layer < m_kLayer; layer++) {
        m_NodeLayer[layer] = 0;
        for (int replica = 0; replica < m_kReplicaMax; replica++) {
            m_NodeReplica[layer][replica] = 0;
            //Long Peixun's update: TGeoPhysicalNode* -> MdcPhysicalPointer
            m_PhysicalReplica[layer][replica].set(0);
        }
    }
 
    for (int i = 0; i < m_kCorrectLayer; i++) {
        m_CorrectMap[m_kiCorrectLayer[i]] = m_kiCorrectReplica[i];
    }
 
    for (int layer = 0; layer < m_kTrueLayer; layer++) {
        for (int replica = 0; replica < m_kReplicaMax; replica++) {
            m_Mdc2DWire[layer][replica] = 0;
        }
    }
}

~MdcROOTGeo() [2/2]

MdcROOTGeo::~MdcROOTGeo

(

)

Destructor.

Member Function Documentation

ClearHits()

void MdcROOTGeo::ClearHits

(

)

Set all physicalNodes corresponding to digiCol;.

Definition at line 677 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    //Long Peixun's update: Annotate cout
    // cout << "m_HitsArray size(): "   << m_HitsArray->GetEntries()   << endl;
    // cout << "m_2DHitsArray size(): " << m_2DHitsArray->GetEntries() << endl;
 
    //Long Peixun's update: the element type of m_HitsArray is replaced by MdcPhysicalPointer
    for (int i = 0; i < m_HitsArray->GetEntries(); i++)
    {
        MdcPhysicalPointer* PphyNode = (MdcPhysicalPointer*)m_HitsArray->At(i);
        (*PphyNode)->SetVisibility(0);
    }
    m_HitsArray->Clear("C");
 
    // set previous event 2D hits info to default
    for (int i = 0; i < m_2DHitsArray->GetEntries(); i++) {
        Mdc2DWire *aWire = (Mdc2DWire*)m_2DHitsArray->At(i);
        aWire->ClearInfo();
        aWire->AddInfo(aWire->GetTitle());
        aWire->CloseInfo();
    }
    m_2DHitsArray->Clear("C");
}

Referenced by BesEvent::ClearHits(), and SetHits().

CorrectReplica() [1/2]

int MdcROOTGeo::CorrectReplica	(	int	layer,
		int	replica )

Correct some axial layer id to copyNo;.

Definition at line 103 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  int rep = replica;
  for (intMap::iterator iter = m_CorrectMap.begin(); iter != m_CorrectMap.end(); iter++) {
    if (layer == (*iter).first ||
    (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut &&
     layer - m_kAxialLayerOut == (*iter).first)) {
      rep -= (*iter).second;
      if (rep < 0) rep += m_kReplica[layer];
    }
  }
 
  //std::cout << "layer" << layer << " " << replica << "->" << rep << std::endl; 
  return rep;
}

Referenced by SetNode().

CorrectReplica() [2/2]

int MdcROOTGeo::CorrectReplica	(	int	layer,
		int	replica )

Correct some axial layer id to copyNo;.

Draw()

void MdcROOTGeo::Draw

(

Option_t *

option

)

Draw function.

Definition at line 1074 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    TString opt = option;
    opt.ToUpper();
 
    if (!m_2DGeoInit) cout << "MdcROOTGeo::Draw2D(), 2D Geometry not initialized!" << endl;
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
    if (!view) cout << "MdcROOTGeo::Draw(), BesView not found" << endl;
 
    if (view->GetVisMdcGlobal()) { // MdcVisGlobal
        if (opt.Contains("XY")) {
            m_MdcXY->Draw("");
            //Long Peixun's update: Annotate
            // for (Int_t i = 0; i < 4; i++) {
            //     m_MdcXYSuper[i]->Draw("");
            // }
        }
 
        if (opt.Contains("ZR")) {
            for (Int_t i = 0; i < 2; i++) {
                //m_MdcZR[i]->SetRotatable(true);   //Long Peixun's update: Move to initializer
                m_MdcZR[i]->Draw("");
            }
        }
 
        TString wireOpt;
        if (view->GetVisMdcTubes()) wireOpt += "TUBE";
        if (view->GetVisMdcWires()) wireOpt += ",WIRE";
 
        Int_t replicaDraw = 0; // change draw sequence in a circle, (-1: 0->Max; 1: Max->0)
        //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
        for (Int_t layer = m_StartLayer; layer < m_kTrueLayer; layer++) {
            Int_t simuLayer = GetSimuLayer(layer);
            for (int replica = 0; replica < m_kReplica[simuLayer]; replica++) {
                if (m_kStereoDir[layer] == -1) replicaDraw = replica;
                else if (m_kStereoDir[layer] == 1) replicaDraw = m_kReplica[simuLayer]-1 - replica;
                else replicaDraw = replica;
 
                if (m_Mdc2DWire[layer][replicaDraw]) {
                    //cout << layer << " " << replicaDraw << endl;
                    //m_Mdc2DWire[layer][replicaDraw]->SetHighlighted(false);
                    m_Mdc2DWire[layer][replicaDraw]->SetFired(false);
                    m_Mdc2DWire[layer][replicaDraw]->Draw(wireOpt);
                }
            }
        }
    }
}

Referenced by BesGeometry::Draw2DXY(), and BesGeometry::Draw2DZR().

DrawHits()

void MdcROOTGeo::DrawHits

(

Option_t *

option

)

Draw 2D hits.

Definition at line 1124 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    // cout << "Draw Hits" << endl;
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
    if (!view) cout << "MdcROOTGeo::DrawHits(), BesView not found" << endl;
 
    //cout << "VisMdcHitsGlobal " << view->GetVisMdcHitsGlobal() << endl;
    //cout << "VisMdcHits       " << view->GetVisMdcHits() << endl;
 
 
    if (view->GetVisMdcHitsGlobal()) {
        TString wireOpt("");
        if (view->GetVisMdcTubes()) wireOpt += ",TUBE";
        if (view->GetVisMdcHits())  wireOpt += ",WIRE";
 
        //Long Peixun's update: Use m_2DHitsArray which contains ready hits rather than rebuild hits
        for (int i = 0; i < m_2DHitsArray->GetEntries(); ++i)
        {
            Mdc2DWire* aWire = (Mdc2DWire*)m_2DHitsArray->At(i);
            if (aWire)
            {
                if ((k_TFire) && (aWire->GetTimeChannel() == 0x7FFFFFFF))
                    continue;
                if ((k_QFire) && (aWire->GetChargeChannel() == 0x7FFFFFFF))
                    continue;
                if ((k_QNotOverflow) && aWire->GetQOverflow() )
                    continue;
                aWire->SetFired(true);
                aWire->SetColorfulWire(k_ColorfulWire, k_MdcTimeSubEvTime);
                aWire->Draw(wireOpt);
            }
        }
 
/*         if (m_MdcDigiCol) {
            for (int i = 0; i < m_MdcDigiCol->GetEntries(); i++) {
                TMdcDigi *aMdcDigi = (TMdcDigi*)m_MdcDigiCol->At(i);
 
                Identifier aMdcID( aMdcDigi->getIntId() );
                int layer  = MdcID::layer( aMdcID );
                int wire   = MdcID::wire( aMdcID );
                unsigned int overflow = aMdcDigi->getOverflow();
 
                Mdc2DWire *aWire = 0;
                aWire = m_Mdc2DWire[layer][wire];
                if (aWire){
                    if ((k_TFire) && (aMdcDigi->getTimeChannel() == 0x7FFFFFFF)){
                        continue;
                    }
                    if ((k_QFire) && (aMdcDigi->getChargeChannel() == 0x7FFFFFFF)){
                        continue;
                    }
                    //if ((!k_TOverflow) && ((overflow&1)>0) ){
                    //    continue;
                    //}
                    if ((k_QNotOverflow) && ((overflow&2)>0) ){
                        continue;
                    }
                    aWire->SetFired(true);
                    aWire->SetColorfulWire(k_ColorfulWire,k_MdcTimeSubEvTime);
                    aWire->Draw(wireOpt);
                }
            }
        } */
    }
}

Referenced by BesEvent::DrawHits().

Get2DWire()

Mdc2DWire * MdcROOTGeo::Get2DWire	(	Int_t	layer,
		Int_t	replica )

Get Mdc2DWire;.

Definition at line 1044 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    if (m_Mdc2DWire[layer][wire]) return m_Mdc2DWire[layer][wire];
    else return 0;
}

Referenced by BesEvent::ConstructMdcTrackFromRec().

GetLayer() [1/2]

TGeoNode * MdcROOTGeo::GetLayer

(

int

layer

)

Get layer node;

Definition at line 334 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  std::stringstream osname;
  if (layer >= 0 && layer < m_kStereoLayerIn) {
    osname << "pv_" << "logical" << "Mdc" << "Stereo" << "Layer" << layer << "_" << layer; 
  }
  else if (layer >= m_kStereoLayerIn && layer < m_kStereoLayerIn + m_kAxialLayerIn) {
    osname << "pv_" << "logical" << "Mdc" << "Axial" << "Layer" << layer << "_" << layer; 
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
       layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) {
    osname << "pv_" << "logical" << "Mdc" << "Stereo" << "Layer" << layer << "_" << layer;
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut &&
       layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut) {
    osname << "pv_" << "logical" << "Mdc" << "Axial" << "Layer" << layer << "_0" << "_" 
       << 2*layer - m_kStereoLayerIn - m_kAxialLayerIn - m_kStereoLayerOut; 
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut &&
       layer <  m_kLayer) {
    osname << "pv_" << "logical" << "Mdc" << "Axial" << "Layer" << layer - m_kAxialLayerOut << "_1" << "_" 
       << 2*(layer - m_kAxialLayerOut) - m_kStereoLayerIn - m_kAxialLayerIn - m_kStereoLayerOut + 1;
  }
  
  //std::cout << osname.str() << std::endl;
  return GetNode( osname.str() );
}

Referenced by Init2DGeometry(), SetAllVisible(), SetDefaultVis(), SetHalfVisible(), SetPhysicalNode(), SetQuarterVisible(), and SetVolumeDefaultVis().

GetLayer() [2/2]

TGeoNode * MdcROOTGeo::GetLayer

(

int

layer

)

Get layer node;

GetLayerNb() [1/2]

int MdcROOTGeo::GetLayerNb

(

)

Get number of layers;.

GetLayerNb() [2/2]

int MdcROOTGeo::GetLayerNb ( )

inline

Get number of layers;.

Definition at line 99 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/BesVisLib/MdcROOTGeo.h.

{ return m_kLayer; }

GetPhysicalReplica() [1/2]

TGeoPhysicalNode * MdcROOTGeo::GetPhysicalReplica	(	int	layer,
		int	replica )

Get replica physical node;.

Definition at line 376 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  if (m_PhysicalReplica[layer][replica] != 0) {
    return m_PhysicalReplica[layer][replica];
  }
  else {
    std::cout << "PhysicalNode: " << "Layer" << layer << "Replica" << replica << " not found" << std::endl;
    return 0;
  } 
}

Referenced by BesEvent::ConstructMdcTrackFromRec(), BesGeometry::GetPhysicalMdcReplica(), Init2DGeometry(), BesVisDisplay::SetMdcFiredCell(), and SetPhysicalDefaultVis().

GetPhysicalReplica() [2/2]

TGeoPhysicalNode * MdcROOTGeo::GetPhysicalReplica	(	int	layer,
		int	replica )

Get replica physical node;.

GetPhysicalSegment()

TGeoPhysicalNode * MdcROOTGeo::GetPhysicalSegment

(

int

segment

)

Get segment physical node;.

Definition at line 1019 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    if (m_PhysicalSegment[segment] != 0) {
        return m_PhysicalSegment[segment];
    }
    else {
        //std::cout << "PhysicalNode: " << "Segment" << segment << " not found" << std::endl;
        return 0;
    }
}

GetReplica() [1/2]

TGeoNode * MdcROOTGeo::GetReplica	(	int	layer,
		int	replica )

Get replica node;.

Definition at line 363 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  if (m_NodeReplica[layer][replica] != 0) {
    //std::cout << " replica " << layer <<  "  " << replica << " found " << std::endl;
    return m_NodeReplica[layer][replica];
  }
  else {
    std::cout << "Node: " << "Layer" << layer << "Replica" << replica << " not found" << std::endl;
    return 0;
  } 
}

Referenced by SetAllVisible(), SetDefaultVis(), SetHalfVisible(), SetPhysicalNode(), SetQuarterVisible(), and SetVolumeDefaultVis().

GetReplica() [2/2]

TGeoNode * MdcROOTGeo::GetReplica	(	int	layer,
		int	replica )

Get replica node;.

GetReplicaNb() [1/2]

int MdcROOTGeo::GetReplicaNb ( int layer )

inline

Get number of replicas on each layer;.

Definition at line 67 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/ROOTGeo/MdcROOTGeo.h.

{ return m_kReplica[layer]; }

GetReplicaNb() [2/2]

int MdcROOTGeo::GetReplicaNb ( int layer )

inline

Get number of replicas on each layer;.

Definition at line 102 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/BesVisLib/MdcROOTGeo.h.

{ return m_kReplica[layer]; }

GetSegment() [1/2]

TGeoNode * MdcROOTGeo::GetSegment	(	int	segment,
		int	no )

Get segment node;.

Definition at line 320 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  std::stringstream osname;
  if (segment == 1 || segment == 2) {
    osname << "pv_" << "logical" << "Mdc" << "Segment" << segment << "_" << m_kLayer + segment - 1; 
  }
  else {
    osname << "pv_" << "logical" << "Mdc" << "Segment" << segment << "_" << m_kLayer + 2 + (segment - 3)*2 + no; 
  }    
 
  return GetNode( osname.str() );
}

Referenced by SetAllVisible(), SetDefaultVis(), SetHalfVisible(), SetQuarterVisible(), and SetVolumeDefaultVis().

GetSegment() [2/2]

TGeoNode * MdcROOTGeo::GetSegment	(	int	segment,
		int	no )

Get segment node;.

GetSegmentNb() [1/2]

int MdcROOTGeo::GetSegmentNb

(

)

Get number of segments;.

GetSegmentNb() [2/2]

int MdcROOTGeo::GetSegmentNb ( )

inline

Get number of segments;.

Definition at line 96 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/BesVisLib/MdcROOTGeo.h.

{ return m_kSegment; }

GetSimuLayer()

Int_t MdcROOTGeo::GetSimuLayer

(

Int_t

trueLayer

)

Trans trueLayer to simuLayer (0~42)->(0~49)

Definition at line 1051 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    if (trueLayer < 36 || trueLayer >= m_kTrueLayer) return trueLayer;
    if (trueLayer >= 36) return (36 + 2*(trueLayer-36));
    return trueLayer;
}

Referenced by Draw(), Init2DGeometry(), and SetHits().

GetVolumeLayer() [1/2]

TGeoVolume * MdcROOTGeo::GetVolumeLayer

(

int

layer

)

Get layer volume;

Definition at line 268 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  std::stringstream osname;
  if (layer >= 0 && layer < m_kStereoLayerIn) {
    osname << "logical" << "Mdc" << "Stereo" << "Layer" << layer; 
  }
  else if (layer >= m_kStereoLayerIn && layer < m_kStereoLayerIn + m_kAxialLayerIn) {
    osname << "logical" << "Mdc" << "Axial" << "Layer" << layer; 
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
       layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) {
    osname << "logical" << "Mdc" << "Stereo" << "Layer" << layer;
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut &&
       layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut) {
    osname << "logical" << "Mdc" << "Axial" << "Layer" << layer << "_0";
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut &&
       layer <  m_kLayer) {
    osname << "logical" << "Mdc" << "Axial" << "Layer" << layer - m_kAxialLayerOut << "_1";
  }
  
  return GetLogicalVolume( osname.str() );
}

Referenced by BesGeometry::GetVolumeMdcLayer(), SetAllVisible(), SetDefaultVis(), SetHalfVisible(), SetQuarterVisible(), and SetVolumeDefaultVis().

GetVolumeLayer() [2/2]

TGeoVolume * MdcROOTGeo::GetVolumeLayer

(

int

layer

)

Get layer volume;

GetVolumeMdc() [1/2]

TGeoVolume * MdcROOTGeo::GetVolumeMdc ( )

inline

Get Mdc volume;.

Definition at line 70 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/ROOTGeo/MdcROOTGeo.h.

{ return m_Mdc; }

Referenced by BesEvent::ConstructMcTrack(), BesEvent::ConstructMdcTrackFromRec(), BesGeometry::GetVolumeMdc(), Init2DGeometry(), and BesGeometry::InitFromGDML().

GetVolumeMdc() [2/2]

TGeoVolume * MdcROOTGeo::GetVolumeMdc ( )

inline

Get Mdc volume;.

Definition at line 108 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/BesVisLib/MdcROOTGeo.h.

{ return m_Mdc; }

GetVolumeReplica() [1/2]

TGeoVolume * MdcROOTGeo::GetVolumeReplica

(

int

layer

)

Get replica volume;.

Definition at line 294 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  std::stringstream osname;
  if (layer >= 0 && layer < m_kStereoLayerIn) {
    osname << "logical" << "Mdc" << "Stereo" << "Layer" << layer << "TwistedTubs"; 
  }
  else if (layer >= m_kStereoLayerIn && layer < m_kStereoLayerIn + m_kAxialLayerIn) {
    osname << "logical" << "Mdc" << "Axial" << "Layer" << layer << "Replica"; 
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
       layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) {
    osname << "logical" << "Mdc" << "Stereo" << "Layer" << layer << "TwistedTubs";
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut &&
       layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut) {
    osname << "logical" << "Mdc" << "Axial" << "Layer" << layer << "_0" << "Replica";
  }
  else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut &&
       layer <  m_kLayer) {
    osname << "logical" << "Mdc" << "Axial" << "Layer" << layer - m_kAxialLayerOut << "_1" << "Replica";
  }
  
  return GetLogicalVolume( osname.str() );
}

Referenced by BesGeometry::GetVolumeMdcReplica(), SetAllVisible(), SetDefaultVis(), SetHalfVisible(), SetQuarterVisible(), and SetVolumeDefaultVis().

GetVolumeReplica() [2/2]

TGeoVolume * MdcROOTGeo::GetVolumeReplica

(

int

layer

)

Get replica volume;.

GetVolumeSegment() [1/2]

TGeoVolume * MdcROOTGeo::GetVolumeSegment

(

int

segment

)

Get segment volume;.

Definition at line 260 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  std::stringstream osname;
  osname << "logical" << "Mdc" << "Segment" << segment;
  return GetLogicalVolume( osname.str() );
}

Referenced by BesGeometry::GetVolumeMdcSegment(), SetAllVisible(), SetDefaultVis(), SetHalfVisible(), SetQuarterVisible(), and SetVolumeDefaultVis().

GetVolumeSegment() [2/2]

TGeoVolume * MdcROOTGeo::GetVolumeSegment

(

int

segment

)

Get segment volume;.

Init2DGeometry()

void MdcROOTGeo::Init2DGeometry

(

)

Initialize 2D Geometry.

Definition at line 126 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    if (m_ROOTGeoInit != 1 && m_ROOTGeoInit != 2) {
        cout << "MdcROOTGeo::Init2DGeometry, ROOT Geometry not Initialized yet!" << endl;
        return;
    }
 
    m_2DGeoInit = 1;
    Int_t mdcColor = 11; //1002;  // 41, 29
    Int_t mdcLineColor = 15;
    Int_t mdcXYStyle = 1001;//3001;
    Int_t mdcZRStyle = 1001;//3007;
    Int_t mdcStereoSuperColor = 38;
    Int_t mdcAxialSuperColor  = 4;
 
    Double_t local[3] = {0.0, 0.0, 0.0};
    Double_t master[3] = {0.0, 0.0, 0.0};
    //Int_t    nPoints = 4;
    Double_t  P[306] = {0.0};   //Long Peixun's update: 300 -> 306, at most 102 points
    Double_t  center[3] = {0.0, 0.0, 0.0};
    TString  name;
 
    //----------XY-----------
    // Mdc
    TGeoTube *mdcShape = (TGeoTube*)GetVolumeMdc()->GetShape();
    //Long Peixun's update: Consider CGEM detector
    Double_t rmin = mdcShape->GetRmin();
    if (m_StartLayer > 0)   //if contains CGEM
    {
        TGeoTube *mdcShapeOuter = (TGeoTube*)(GetLayer(m_StartLayer)->GetVolume())->GetShape();
        rmin = mdcShapeOuter->GetRmin();
    }
 
    m_MdcXY = new BesCircle2D("Mdc", "Mdc", rmin, mdcShape->GetRmax(), &center[0]);
    m_MdcXY->SetNSegment(360);
    m_MdcXY->SetFillColor(mdcColor);
    m_MdcXY->SetFillStyle(mdcXYStyle);
    m_MdcXY->SetLineColor(mdcLineColor);
 
    //Long Peixun's update: Annotate
/*     // Mdc Supers
    Int_t layerIn[4]  = {0,
                         m_kStereoLayerIn,
                         m_kStereoLayerIn+m_kAxialLayerIn,
                         m_kStereoLayerIn+m_kAxialLayerIn+m_kStereoLayerOut
                        };
    Int_t layerOut[4] = {m_kStereoLayerIn-1,
                         m_kStereoLayerIn+m_kAxialLayerIn-1,
                         m_kStereoLayerIn+m_kAxialLayerIn+m_kStereoLayerOut-1,
                         m_kLayer-1
                        };
    for (Int_t i = 0; i < 4; i++) {
        TGeoTube *mdcSuperShapeIn  = (TGeoTube*)(GetLayer(layerIn[i])->GetVolume())->GetShape();
        TGeoTube *mdcSuperShapeOut = (TGeoTube*)(GetLayer(layerOut[i])->GetVolume())->GetShape();
        name = TString("MdcSuper" + i);
        m_MdcXYSuper[i] = new BesCircle2D(name, name, mdcSuperShapeIn->GetRmin(), mdcSuperShapeOut->GetRmax(), &center[0]);
        if (i%2 == 0) m_MdcXYSuper[i]->SetFillColor(mdcStereoSuperColor);
        else m_MdcXYSuper[i]->SetFillColor(mdcAxialSuperColor);
    } */
 
    // Mdc Wires
    //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
    for (Int_t layer = m_StartLayer; layer < m_kTrueLayer; layer++) {
        Int_t simuLayer = GetSimuLayer(layer);
        for (int replica = 0; replica < m_kReplica[simuLayer]; replica++) {
            //cout << "layer " << layer << " replica " << replica << endl;
            TGeoPhysicalNode *wirePhyNode = GetPhysicalReplica(simuLayer, replica);
            name = TString("Mdc Layer ");
            name += layer;
            name += " Wire ";
            name += replica;
 
            if (m_kStereoDir[layer] != 0) { // stereo Layer
                TGeoArb8 *wireShape = (TGeoArb8*)wirePhyNode->GetShape();
                Double_t *localArb8Point, masterArb8Point[24];
                localArb8Point = wireShape->GetVertices();
                for (Int_t i = 0; i < 8; i++) {
                    local[0] = localArb8Point[2*i];
                    local[1] = localArb8Point[2*i+1];
                    if (i < 4) local[2] = wireShape->GetDz() * (-1.0);
                    else local[2] = wireShape->GetDz();
 
                    wirePhyNode->GetMatrix(-1*wirePhyNode->GetLevel())
                    ->LocalToMaster(local, &master[0]); // transform to top
                    for (Int_t j = 0; j < 3; j++) {
                        masterArb8Point[3*i+j] = master[j];
                    }
                }
 
                m_Mdc2DWire[layer][replica] = new Mdc2DWire(name, name, 8, &masterArb8Point[0]);
        
            }
            else { // axial Layer
 
                TGeoTubeSeg *wireShape = (TGeoTubeSeg*)wirePhyNode->GetShape();
 
                Double_t centerR = 0.5*(wireShape->GetRmin() + wireShape->GetRmax());
                if (layer >= m_kStereoLayerIn+m_kAxialLayerIn+m_kStereoLayerOut)
                    centerR = wireShape->GetRmax();
                Double_t centerPhi =
                    0.5*(wireShape->GetPhi1() + wireShape->GetPhi2()) * TMath::DegToRad();
                local[0] = centerR * cos(centerPhi);
                local[1] = centerR * sin(centerPhi);
                local[2] = 0.0;
                wirePhyNode->GetMatrix(-1*wirePhyNode->GetLevel())
                ->LocalToMaster(local, &master[0]); // transform to top
 
                Double_t rmin = wireShape->GetRmin();
                Double_t rmax = wireShape->GetRmax();
                if (layer >= m_kStereoLayerIn+m_kAxialLayerIn+m_kStereoLayerOut)
                    rmax = rmin + 2*(rmax-rmin);
                m_Mdc2DWire[layer][replica] =
                    new Mdc2DWire( name, name, rmin, rmax, wireShape->GetDz(),centerPhi, &master[0] );
            }
        }
    }
 
    //----------ZR-----------
    // Mdc
    Int_t iPoint = 0;
    TGeoTube *aMdcLayer = 0;
 
    //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
    for (Int_t layer = m_StartLayer; layer < m_kLayer; layer++) {
        aMdcLayer = (TGeoTube*)GetLayer(layer)->GetVolume()->GetShape();
        P[3*iPoint] = 0.0;
        P[3*iPoint+1] = aMdcLayer->GetRmin();
        P[3*iPoint+2] = aMdcLayer->GetDz();
        iPoint++;
    }
 
    P[3*iPoint] = 0.0;
    P[3*iPoint+1] = aMdcLayer->GetRmax();
    P[3*iPoint+2] = aMdcLayer->GetDz();
    iPoint++;
    P[3*iPoint] = 0.0;
    P[3*iPoint+1] = aMdcLayer->GetRmax();
    P[3*iPoint+2] = aMdcLayer->GetDz() * (-1.0);
    iPoint++;
 
    for (Int_t layer = m_kLayer-1; layer >= m_StartLayer; layer--) {
        aMdcLayer = (TGeoTube*)GetLayer(layer)->GetVolume()->GetShape();
        P[3*iPoint] = 0.0;
        P[3*iPoint+1] = aMdcLayer->GetRmin();
        P[3*iPoint+2] = aMdcLayer->GetDz() * (-1.0);
        iPoint++;
    }
 
 
    m_MdcZR[0] = new BesPolygon2D("MdcZRUp", "MdcZRUp", iPoint, &P[0]);
 
    for (Int_t i = 0; i < iPoint; i++)
        P[3*i+1] *= -1.0;
 
    m_MdcZR[1] = new BesPolygon2D("MdcZRDown", "MdcZRDown", iPoint, &P[0]);
 
    for (Int_t i = 0; i < 2; i++) {
        m_MdcZR[i]->SetFillColor(mdcColor);
        m_MdcZR[i]->SetFillStyle(mdcZRStyle);
        m_MdcZR[i]->SetLineColor(mdcLineColor);
        m_MdcZR[i]->SetRotatable(true);     //Long Peixun's update
    }
}

Referenced by BesGeometry::InitGeometry().

InitFromGDML()

void MdcROOTGeo::InitFromGDML	(	const char *	gdmlFile,
		const char *	setupName )

Initialize ROOTGeo from GDML.

Definition at line 110 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    m_ROOTGeoInit = 2;
    ReadGdml(gdmlFile, setupName);
    SetNode();
}

Referenced by BesGeometry::InitFromGDML().

InitFromGdml()

void MdcROOTGeo::InitFromGdml	(	const char *	gdmlFile,
		const char *	setupName )

Initialize the instance of ROOTGeo.

Definition at line 55 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  ReadGdml(gdmlFile, setupName);
  SetNode();
  SetDefaultVis();
 
  m_ROOTGeoInit = 1;  
}

InitFromROOT()

void MdcROOTGeo::InitFromROOT

(

TGeoVolume *

vol

)

Initialize ROOTGeo from TGeoVolume logicalMdc.

Definition at line 118 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    m_ROOTGeoInit = 1;
    SetVolumeMdc(vol);
    SetNode();
}

Referenced by BesGeometry::InitFromROOT().

IsHit()

Bool_t MdcROOTGeo::IsHit	(	Int_t	layer,
		Int_t	wire )

Judge whether the digiCol contains (layer, wire)

Definition at line 1059 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    if (m_MdcDigiCol) {
        for (int i = 0; i < m_MdcDigiCol->GetEntries(); i++) {
            Identifier mdcID( ((TMdcDigi*)m_MdcDigiCol->At(i))->getIntId() );
            if (layer == MdcID::layer( mdcID ) && wire == MdcID::wire( mdcID )) {
                return kTRUE;
            }
        }
    }
 
    return kFALSE;
}

SetAllVisible() [1/2]

void MdcROOTGeo::SetAllVisible

(

)

Set all visible;.

Definition at line 170 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  for (int segment = 1; segment <= m_kSegment; segment++) {
    GetVolumeSegment(segment)->SetVisibility(0);
  }
  
  for (int layer = 0; layer < m_kLayer; layer++) {
    GetVolumeLayer(layer)->SetVisibility(1);
    GetVolumeReplica(layer)->SetVisibility(0);
  }
 
  for (int segment = 1; segment <= m_kSegment; segment++) {
    for (int no = 0; no < 2; no++) {
      GetSegment(segment, no)->SetVisibility(1);
    }
  }
 
  for (int layer = 0; layer < m_kLayer; layer++) {
    GetLayer(layer)->SetVisibility(1);
    for (int replica = 0; replica < m_kReplica[layer]; replica++) {
      GetReplica(layer, replica)->SetVisibility(0);
    }
  }
}

SetAllVisible() [2/2]

void MdcROOTGeo::SetAllVisible

(

)

Set all visible;.

SetColorfulWire()

void MdcROOTGeo::SetColorfulWire

(

Bool_t

input

)

Definition at line 860 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

                                       {
    k_ColorfulWire = input;
}

SetDefaultVis()

void MdcROOTGeo::SetDefaultVis

(

)

Set default visual attributes;.

Definition at line 120 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  //std::cout << "begin of set defaultvis" << std::endl;
  int mdcColor         = 3;
  int segmentColor     = 0;
  int hypeColor        = 3;
  int tubeColor        = 4;
  int twistedTubsColor = 3;
  int replicaColor     = 4;
 
  m_Mdc->SetLineColor(mdcColor);
  m_Mdc->SetVisibility(0);
  
  for (int segment = 1; segment <= m_kSegment; segment++) {
    GetVolumeSegment(segment)->SetLineColor(segmentColor);
    GetVolumeSegment(segment)->SetVisibility(0);
  }
 
  for (int layer =0; layer < m_kLayer; layer++) {
    if ( (layer >= 0 && layer < m_kStereoLayerIn) ||
     (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
      layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) ) {
      GetVolumeLayer(layer)->SetLineColor(hypeColor);
      GetVolumeReplica(layer)->SetLineColor(twistedTubsColor);
    }
    else {
      GetVolumeLayer(layer)->SetLineColor(tubeColor);
      GetVolumeReplica(layer)->SetLineColor(replicaColor);
    }
    GetVolumeLayer(layer)->SetVisibility(0);
    GetVolumeReplica(layer)->SetVisibility(1);
  }
 
  for (int segment = 1; segment <= m_kSegment; segment++) {
    for (int no = 0; no < 2; no++) {
      GetSegment(segment, no)->SetVisibility(0);
    }
  }
 
  for (int layer = 0; layer < m_kLayer; layer++) {
    GetLayer(layer)->SetVisibility(0);
    for (int replica = 0; replica < m_kReplica[layer]; replica++) {
      GetReplica(layer, replica)->SetVisibility(0);
    }
  }
 
  //std::cout << "end of set defaultvis" << std::endl;
}

Referenced by InitFromGdml().

SetDetector()

void MdcROOTGeo::SetDetector

(

)

Draw Detecor (what is detector depends on you)

Definition at line 629 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    BesView *view = 0;
    if (gPad) view = dynamic_cast<BesView*>(gPad->GetView());
    //if (view) cout << "viewVisFull3DMdc " << view->GetVisFull3DMdc() << endl;
 
    m_DetectorsArray->Clear();
    for (int segment = 0; segment < 2*m_kSegment-2; segment++) {
        TGeoPhysicalNode *phyNode = 0;
        phyNode = m_PhysicalSegment[segment];
        //cout << m_PhysicalSegment[segment]->GetName() << endl;
        if (phyNode) {
            phyNode->SetVisibility(0); // set all invisible before set any visible
            if ( (segment >= 2 && segment <= 3) || segment > 59) {
                m_DetectorsArray->Add( phyNode );
            }
            else if (view && view->GetVisFull3DMdc()) {
                m_DetectorsArray->Add( phyNode );
            }
        }
    }
    /*
    for (int part = 0; part < m_kPart; part++) {
      int nPhi = (part == 1 ? m_kPhiBr : m_kPhiEc);
      for (int phi = 0; phi < nPhi; phi++) {
        int nTheta = (part == 1 ? m_kThetaBr : m_kThetaEc);
        for (int theta = 0; theta < nTheta; theta++) {
      // This is you detector, set it on
      if ( (part == 1 && theta != 0 && theta != 43 && ((phi != 0 && phi != 60) || theta%2 == 0)) ||
           (part != 1 && theta != 0) ) {
        m_PhysicalCrystal[part][phi][theta]->SetVisibility(0);
      }
      
      // This is you hit
      // when you are to set fired hits to another color, say red,
      // SetIsVolAtt(kFALSE) and to then SetLineColor(kRed);
      if (phi == 12 && theta == 5) {
        m_PhysicalCrystal[part][phi][theta]->SetVisibility(1);
        m_PhysicalCrystal[part][phi][theta]->SetLineColor(2);
      }
        }
      }
    }
    std::cout << "end of set detector on" << std::endl;
    */
}

Referenced by BesClient::HandleViewOptionMenu().

SetHalfVisible()

void MdcROOTGeo::SetHalfVisible

(

)

Set half visible;.

Definition at line 512 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    for (int segment = 1; segment <= m_kSegment; segment++) {
        GetVolumeSegment(segment)->SetVisibility(0);
    }
 
    //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
    for (int layer = m_StartLayer; layer < m_kLayer; layer++) {
        GetVolumeLayer(layer)->SetVisibility(0);
        GetVolumeReplica(layer)->SetVisibility(1);
    }
 
    for (int segment = 1; segment <= m_kSegment; segment++) {
        for (int no = 0; no < 2; no++) {
            GetSegment(segment, no)->SetVisibility(0);
        }
    }
 
    //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
    for (int layer = m_StartLayer; layer < m_kLayer; layer++) {
        GetLayer(layer)->SetVisibility(0);
        for (int replica = 0; replica < m_kReplica[layer]; replica++) {
            if (replica < m_kReplica[layer]/4 ||
                    replica > m_kReplica[layer]*3/4 ) GetReplica(layer, replica)->SetVisibility(1);
            else                                  GetReplica(layer, replica)->SetVisibility(0);
        }
    }
}

SetHits()

void MdcROOTGeo::SetHits

(

)

Definition at line 702 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    // set previous event hits to default vis
    //Long Peixun's update: move to ClearHits()
    ClearHits();
 
    // set new hits
    //Long Peixun's update: Consider gEvent == NULL
    int NDigiCol = 0;
    if (gEvent) 
    {
        m_MdcDigiCol = gEvent->GetMdcDigiCol();
        //Long Peixun's update: Consider m_MdcDigiCol == NULL
        if (m_MdcDigiCol) NDigiCol = m_MdcDigiCol->GetEntries();
        else NDigiCol = 0;
    }
 
    for (int i = 0; i < NDigiCol; i++) {
        TMdcDigi *aMdcDigi = (TMdcDigi*)m_MdcDigiCol->At(i);
        Identifier aMdcID( aMdcDigi->getIntId() );
        int layer = MdcID::layer( aMdcID );
        int wire  = MdcID::wire( aMdcID );
 
        Double_t time   = RawDataUtil::MdcTime(aMdcDigi->getTimeChannel());
        Int_t charge = RawDataUtil::MdcCharge(aMdcDigi->getChargeChannel());
        unsigned int overflow = aMdcDigi->getOverflow();
 
        //Long Peixun's update: the element type of m_HitsArray is replaced by MdcPhysicalPointer
        MdcPhysicalPointer* phyNode = 0;
        if (layer < 36)
        {   // number of layers in the three inner
            phyNode = &m_PhysicalReplica[layer][wire];
            // Long Peixun's update: Add additional information into pointer
            phyNode->setTimeChannel(aMdcDigi->getTimeChannel());
            phyNode->setChargeChannel(aMdcDigi->getChargeChannel());
            phyNode->setQOverflow(overflow & 2);
            if (phyNode) m_HitsArray->Add(phyNode);
        }
        else
        {
            phyNode = &m_PhysicalReplica[GetSimuLayer(layer)][wire];
            if (phyNode) m_HitsArray->Add(phyNode);
            phyNode->setTimeChannel(aMdcDigi->getTimeChannel());
            phyNode->setChargeChannel(aMdcDigi->getChargeChannel());
            phyNode->setQOverflow(overflow & 2);
 
            phyNode = &m_PhysicalReplica[GetSimuLayer(layer) + 1][wire];
            if (phyNode) m_HitsArray->Add(phyNode);
            phyNode->setTimeChannel(aMdcDigi->getTimeChannel());
            phyNode->setChargeChannel(aMdcDigi->getChargeChannel());
            phyNode->setQOverflow(overflow & 2);
        }
 
        Mdc2DWire *aWire = 0;
        aWire = m_Mdc2DWire[layer][wire];
        if (aWire) {
            aWire->ClearInfo();
            aWire->AddInfo(aWire->GetTitle());
        double evTime = gEvent->GetHeader().GetEvTime();
 
            char data[200];     //Long Peixun's update: 100 -> 200
        if(time>1000000) {
        sprintf(data, "NO time");
        }else{
        sprintf(data, "time=%-.1f ns", time);
        if(fabs(evTime)>0.0001){ 
            sprintf(data, "%s, Tdrift(Raw)=%-.1f ns", data, time - evTime +230); 
        }
        }
        if(charge>1000000) {
        sprintf(data, "%s, NO charge", data);
        }else{
        sprintf(data, "%s, charge=%d", data,charge);
        }
            aWire->AddInfo( TString(data) );
            sprintf(data, "Fired");
            aWire->AddInfo( TString(data) );
 
        aWire->SetTime(time);//yzhang 
        aWire->SetCharge(charge);//yzhang
        aWire->SetQOverflow(overflow&2);//yzhang
        aWire->SetEvTime(evTime);
 
            //Long Peixun's update: Set raw channel data of aWire
            aWire->SetTimeChannel(aMdcDigi->getTimeChannel());
            aWire->SetChargeChannel(aMdcDigi->getChargeChannel());
 
            aWire->CloseInfo();
 
            m_2DHitsArray->Add(aWire);
        }
    }
}

Referenced by BesEvent::SetHits().

SetMdcTimeSubEvTime()

void MdcROOTGeo::SetMdcTimeSubEvTime

(

Bool_t

input

)

Definition at line 865 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

                                           {
    k_MdcTimeSubEvTime = input;
}

SetNode() [1/2]

void MdcROOTGeo::SetNode

(

)

Set the pointers to theirs nodes;.

Definition at line 65 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  m_Mdc = GetTopVolume();
  if(!m_Mdc) std::cout << "m_Mdc = 0" << std::endl;
  
  for (int layer = 0; layer < m_kLayer; layer++) {
    //std::cout << "Layer " << layer << std::endl;
    for (int replica = 0; replica < m_kReplica[layer]; replica++) {
      std::stringstream osname;
      int rep = replica;
      if (layer >= 0 && layer < m_kStereoLayerIn) {
    osname << "pv_" << "logical" << "Mdc" << "Stereo" << "Layer" << layer << "TwistedTubs" << "_" << rep; 
      }
      else if (layer >= m_kStereoLayerIn && layer < m_kStereoLayerIn + m_kAxialLayerIn) {
    rep = CorrectReplica(layer, rep);
    osname << "pv_" << "logical" << "Mdc" << "Axial" << "Layer" << layer << "Replica" << "_" << rep; 
      }
      else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
           layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) {
    osname << "pv_" << "logical" << "Mdc" << "Stereo" << "Layer" << layer << "TwistedTubs" << "_" << rep;
      }
      else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut &&
           layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut) {
    rep = CorrectReplica(layer, rep);
    osname << "pv_" << "logical" << "Mdc" << "Axial" << "Layer" << layer << "_0" << "Replica" << "_" << rep;
      }
      else if (layer >= m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut + m_kAxialLayerOut &&
           layer <  m_kLayer) {
    rep = CorrectReplica(layer, rep);
    osname << "pv_" << "logical" << "Mdc" << "Axial" << "Layer" << layer - m_kAxialLayerOut << "_1" << "Replica" << "_" << rep;
      }
      
      m_NodeReplica[layer][replica] = GetNode( osname.str() );
    }
  }
}

Referenced by InitFromGDML(), InitFromGdml(), and InitFromROOT().

SetNode() [2/2]

void MdcROOTGeo::SetNode

(

)

Set the pointers to theirs nodes;.

SetPhysicalDefaultVis()

void MdcROOTGeo::SetPhysicalDefaultVis

(

)

Set default physical node attributes;.

Definition at line 592 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    // for (int part = 0; part < m_kPart; part++)
    // {
    //     int nPhi = (part == 1 ? m_kPhiBr : m_kPhiEc);
    //     for (int phi = 0; phi < nPhi; phi++) 
    //     {
    //         int nTheta = (part == 1 ? m_kThetaBr : m_kThetaEc);
    //         for (int theta = 0; theta < nTheta; theta++)
    //         {
    //             m_PhysicalCrystal[part][phi][theta]->SetIsVolAtt(kFALSE);
    //             if (part == 1) m_PhysicalCrystal[part][phi][theta]->SetLineColor(m_brCrystalColor);
    //             else m_PhysicalCrystal[part][phi][theta]->SetLineColor(m_ecCrystalColor);
    //         }
    //     }
    // }
 
    for (int layer = m_StartLayer; layer < m_kLayer; layer++)
    {
        if ( (layer >= 0 && layer < m_kStereoLayerIn) ||
            (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
            layer < m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) )
        {
            for (int replica = 0; replica < m_kReplica[layer]; replica++)
                GetPhysicalReplica(layer, replica)->SetLineColor(m_twistedTubsColor);
        }
        else
        {
            for (int replica = 0; replica < m_kReplica[layer]; replica++)
                GetPhysicalReplica(layer, replica)->SetLineColor(m_replicaColor);
        }
    }
  //std::cout << "end of set defaultvis" << std::endl;
}

Referenced by BesGeometry::Draw3D(), and BesGeometry::SetPhysicalDefaultVis().

SetPhysicalNode() [1/2]

void MdcROOTGeo::SetPhysicalNode

(

)

Set the pointers to the physical nodes;.

Definition at line 223 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  int twistedTubsColor = 3;
  int replicaColor     = 4;
 
  if (gGeoManager == 0) std::cout << "Create gGeoManager first" << std::endl; 
  TGeoNode *bes = gGeoManager->GetTopNode();
  //std::cout << "m_childNo " << m_childNo << std::endl;
  TGeoNode *mdc = bes->GetDaughter(m_childNo);
 
  //if(!m_Mdc) std::cout << "m_Mdc = 0" << std::endl;
  
  for (int layer = 0; layer < m_kLayer; layer++) {
    TGeoNode *nodeLayer = GetLayer(layer);
    //std::cout << "Layer " << layer << std::endl;
    for (int replica = 0; replica < m_kReplica[layer]; replica++) {
      TGeoNode *nodeReplica = GetReplica(layer, replica);
      m_PhysicalReplica[layer][replica] = gGeoManager->MakePhysicalNode( TString("/") + bes->GetName() + 
                                     TString("/") + mdc->GetName() +
                                     TString("/") + nodeLayer->GetName() +
                                     TString("/") + nodeReplica->GetName());
      m_PhysicalReplica[layer][replica]->SetVisibility(0); 
      m_PhysicalReplica[layer][replica]->SetIsVolAtt(kFALSE);
      if ( (layer >= 0 && layer < m_kStereoLayerIn) ||
       (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
        layer < m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) ) {
    m_PhysicalReplica[layer][replica]->SetLineColor(twistedTubsColor);
      }
      else {
    m_PhysicalReplica[layer][replica]->SetLineColor(replicaColor);
      } 
      //if (m_PhysicalReplica[layer][replica]->IsVolAttributes()) std::cout << "yes " << std::endl; 
    }
  }
}

Referenced by BesGeometry::InitGeometry().

SetPhysicalNode() [2/2]

void MdcROOTGeo::SetPhysicalNode

(

)

Set the pointers to the physical nodes;.

SetQFire()

void MdcROOTGeo::SetQFire

(

Bool_t

input

)

Definition at line 854 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

                                {
    k_QFire = input;
    //cout << "Mdc ROOT Geometry ADC Match: " << k_QMatch << endl;
}

SetQNotOverflow()

void MdcROOTGeo::SetQNotOverflow

(

Bool_t

input

)

Definition at line 848 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

                                       {
    k_QNotOverflow = input;
    //cout << "Mdc ROOT Geometry TDC Match: " << k_TMatch << endl;
}

SetQuarterVisible() [1/2]

void MdcROOTGeo::SetQuarterVisible

(

)

Set quater visible;.

Definition at line 196 of file DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx.

{
  for (int segment = 1; segment <= m_kSegment; segment++) {
    GetVolumeSegment(segment)->SetVisibility(0);
  }
  
  for (int layer = 0; layer < m_kLayer; layer++) {
    GetVolumeLayer(layer)->SetVisibility(0);
    GetVolumeReplica(layer)->SetVisibility(1);
  }
 
  for (int segment = 1; segment <= m_kSegment; segment++) {
    for (int no = 0; no < 2; no++) {
      GetSegment(segment, no)->SetVisibility(0);
    }
  }
 
  for (int layer = 0; layer < m_kLayer; layer++) {
    GetLayer(layer)->SetVisibility(0);
    for (int replica = 0; replica < m_kReplica[layer]; replica++) {
      if (replica < m_kReplica[layer]/4) GetReplica(layer, replica)->SetVisibility(0);
      else                               GetReplica(layer, replica)->SetVisibility(1);
    }
  }
}

SetQuarterVisible() [2/2]

void MdcROOTGeo::SetQuarterVisible

(

)

Set quater visible;.

SetTFire()

void MdcROOTGeo::SetTFire

(

Bool_t

input

)

Definition at line 836 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

                                {
    k_TFire= input;
    //cout << "Mdc ROOT Geometry TDC Match: " << k_TMatch << endl;
}

SetVisMdcDetector()

void MdcROOTGeo::SetVisMdcDetector

(

)

Set Mdc default detector visibility;.

Definition at line 797 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    BesView *view = 0;
    if (gPad) view = dynamic_cast<BesView*>(gPad->GetView());
 
    for (int i = 0; i < m_DetectorsArray->GetEntries(); i++) {
        TGeoPhysicalNode *phyNode = (TGeoPhysicalNode*)m_DetectorsArray->At(i);
        if (view && view->GetVisMdcGlobal()) {
            phyNode->SetVisibility(1);
            continue;
        }
        phyNode->SetVisibility(0);
    }
}

Referenced by BesGeometry::Draw3D(), and BesClient::HandleViewOptionMenu().

SetVisMdcHits()

void MdcROOTGeo::SetVisMdcHits

(

)

Set Mdc hits visibility;.

Definition at line 813 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    BesView *view = 0;
    if (gPad) view = dynamic_cast<BesView*>(gPad->GetView());
 
    for (int i = 0; i < m_HitsArray->GetEntries(); i++)
    {
        MdcPhysicalPointer* PphyNode = (MdcPhysicalPointer*)m_HitsArray->At(i);
        if (view && view->GetVisMdcHitsGlobal())
        {
            //Long Peixun's update: Consider overflow options in 3D view
            if ( (k_TFire && PphyNode->getTimeChannel() == 0x7FFFFFFF)
              || (k_QFire && PphyNode->getChargeChannel() == 0x7FFFFFFF)
              || (k_QNotOverflow && PphyNode->getQOverflow())  )
                (*PphyNode)->SetVisibility(0);
            else
                (*PphyNode)->SetVisibility(1);
        }
        else (*PphyNode)->SetVisibility(0);
    }
}

Referenced by BesGeometry::Draw3D().

SetVolumeDefaultVis()

void MdcROOTGeo::SetVolumeDefaultVis

(

)

Set default visual attributes;.

Definition at line 390 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx.

{
    //std::cout << "begin of set defaultvis" << std::endl;
    m_Mdc->SetLineColor(m_MdcColor);
    m_Mdc->SetVisibility(0);
 
    for (int segment = 1; segment <= m_kSegment; segment++) {
        GetVolumeSegment(segment)->SetLineColor(m_segmentColor);
        //if (segment > 3) GetVolumeSegment(segment)->SetVisibility(0);
        GetVolumeSegment(segment)->SetVisibility(1);
    }
 
    //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
    for (int layer = m_StartLayer; layer < m_kLayer; layer++) {
        if ( (layer >= 0 && layer < m_kStereoLayerIn) ||
                (layer >= m_kStereoLayerIn + m_kAxialLayerIn &&
                 layer <  m_kStereoLayerIn + m_kAxialLayerIn + m_kStereoLayerOut) ) {
            GetVolumeLayer(layer)->SetLineColor(m_hypeColor);
            GetVolumeReplica(layer)->SetLineColor(m_twistedTubsColor);
        }
        else {
            GetVolumeLayer(layer)->SetLineColor(m_tubeColor);
            GetVolumeReplica(layer)->SetLineColor(m_replicaColor);
        }
        GetVolumeLayer(layer)->SetVisibility(0);
        GetVolumeReplica(layer)->SetVisibility(1);
    }
 
    for (int segment = 1; segment <= m_kSegment; segment++) {
        for (int no = 0; no < 2; no++) {
            GetSegment(segment, no)->SetVisibility(0);
        }
    }
 
    //Long Peixun's update: Ignore layer 0~7 if detector contains CGEM
    for (int layer = m_StartLayer; layer < m_kLayer; layer++) {
        GetLayer(layer)->SetVisibility(0);
        for (int replica = 0; replica < m_kReplica[layer]; replica++) {
            GetReplica(layer, replica)->SetVisibility(0);
        }
    }
 
    //std::cout << "end of set defaultvis" << std::endl;
}

SetVolumeMdc()

void MdcROOTGeo::SetVolumeMdc ( TGeoVolume * vol )

inline

Set Mdc volume, while initializing from ROOT;.

Definition at line 105 of file EventDisplay/BesVisLib/BesVisLib-00-06-12/BesVisLib/MdcROOTGeo.h.

{ m_Mdc = vol; }

Referenced by InitFromROOT().

---

The documentation for this class was generated from the following files:

• DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/ROOTGeo/MdcROOTGeo.h
• EventDisplay/BesVisLib/BesVisLib-00-06-12/BesVisLib/MdcROOTGeo.h
• DetectorDescription/ROOTGeo/ROOTGeo-00-00-15/src/MdcROOTGeo.cxx
• EventDisplay/BesVisLib/BesVisLib-00-06-12/src/MdcROOTGeo.cxx