#include <CgemROOTGeo.h>

Inheritance diagram for CgemROOTGeo:

Public Member Functions
	CgemROOTGeo ()

virtual	~CgemROOTGeo ()

void	InitFromGDML (const char gdmlFile, const char setupName)

void	InitFromROOT (TGeoVolume *vol)

void	InitGeometry ()

void	SetNode ()

void	SetVolumeDefaultVis ()

void	SetPhysicalDefaultVis ()

void	SetAllVisible ()

void	SetQuarterVisible ()

void	SetHalfVisible ()

void	SetPhysicalNode ()

void	SetDetector ()

void	SetHits ()

void	ClearHits ()

void	SetVisCgemDetector ()

void	SetVisCgemHits ()

void	SetVolumeCgem (TGeoVolume *vol)

TGeoVolume *	GetVolumeCgem () const

Cgem2DStripBase *	Get2DStrip (bool is_x, int layer, int sheet, int strip) const

TGeoNode *	GetCgemAnode (int layer)

TGeoNode *	GetCgemCathode (int layer)

TGeoNode *	GetCgemFoil (int layer, int foil)

void	Draw (Option_t *option)

void	DrawUF (int layer, Option_t *option="")

void	DrawHits (Option_t *option)

void	Draw3DHits (Option_t *option)

double	calcZFromLayerSheetVIDPhi (int layer, int sheet, double v, double phi) const

double	getRFromLayerSheet (int layer, int sheet) const

double	getSheetWidth (int layer, int sheet) const

double	isPointInSheet (int layer, int sheet, double phi, double z) const

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

Definition at line 23 of file CgemROOTGeo.h.

Constructor & Destructor Documentation

◆ CgemROOTGeo()

CgemROOTGeo::CgemROOTGeo ( )

Definition at line 35 of file CgemROOTGeo.cxx.

                         : SubDetectorROOTGeo()
{
    m_CgemColor = kSpring;
    m_XstripColor = 828;
    m_VstripColor = kTeal;
 
    for (int layer = 0; layer < m_kCgemLayer; ++layer)
    {
        m_NodeCgemAnode[layer] = 0;
        m_NodeCgemCathode[layer] = 0;
        for (int foil = 0; foil < m_kCgemFoil; ++foil)
            m_NodeCgemFoil[layer][foil] = 0;
 
        for (int sheet = 0; sheet < 2; ++sheet)
        {
            NXStrip[layer][sheet] = 0;
            NVStrip[layer][sheet] = 0;
            m_CgemXStrip[layer][sheet] = 0;
            m_CgemVStrip[layer][sheet] = 0;
            m_CgemSheet[layer][sheet] = 0;
        }
    }
}

◆ ~CgemROOTGeo()

CgemROOTGeo::~CgemROOTGeo ( )

virtual

Definition at line 60 of file CgemROOTGeo.cxx.

{
    cout << "delete old CgemROOTGeo" << endl;
    for (int layer = 0; layer < m_kCgemLayer; ++layer)
    {
        for (int sheet = 0; sheet < 2; ++sheet)
        {
            if (layer == 0 && sheet == 1) continue;
            delete m_CgemSheet[layer][sheet];
 
            for (int strip = 0; strip < NXStrip[layer][sheet]; ++strip)
                delete m_CgemXStrip[layer][sheet][strip];
            delete m_CgemXStrip[layer][sheet];
            for (int strip = 0; strip < NVStrip[layer][sheet]; ++strip)
                delete m_CgemVStrip[layer][sheet][strip];
            delete m_CgemVStrip[layer][sheet];
        }
    }
}

Member Function Documentation

◆ calcZFromLayerSheetVIDPhi()

double CgemROOTGeo::calcZFromLayerSheetVIDPhi	(	int	layer,
		int	sheet,
		double	v,
		double	phi ) const

Definition at line 907 of file CgemROOTGeo.cxx.

{
    if (layer < 0 || layer >= m_kCgemLayer || sheet < 0 || sheet >= 2) return 0.0;
    double ruf = (RV[layer][sheet] + RX[layer][sheet]) / 2;
    double len = SheetLength[layer][sheet];
    double wid = SheetWidth[layer][sheet];
    double ang = StereoAngle[layer][sheet];
    double xmin = XMin[layer][sheet];
    double zmin = ZMin[layer][sheet];
    double pit = VPitch[layer][sheet];
    double vuf = v / RV[layer][sheet] * ruf;
    CgemVStripCalculator calc(ruf, len, wid, ang, xmin, zmin, pit);
    return calc.getZFromVPhi(v, phi);
}

◆ ClearHits()

void CgemROOTGeo::ClearHits ( )

Definition at line 479 of file CgemROOTGeo.cxx.

{
    for (int i = 0; i < m_2DHitsArray->GetEntries(); ++i)
    {
        Cgem2DStripBase* strip = (Cgem2DStripBase*)m_2DHitsArray->At(i);
        strip->ClearInfo();
        strip->AddInfo(strip->GetTitle());
        strip->CloseInfo();
    }
    m_2DHitsArray->Clear("C");
}

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

◆ Draw()

void CgemROOTGeo::Draw ( Option_t * option )

Definition at line 596 of file CgemROOTGeo.cxx.

{
    TString opt = option;
    opt.ToUpper();
 
    if (!m_2DGeoInit) cout << "CgemROOTGeo::Draw(), Geometry not initialized!" << endl;
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
    if (!view) cout << "CgemROOTGeo::Draw(), BesView not found" << endl;
 
    if (view->GetVisCgemGlobal())
    {
        if (opt.Contains("XY")) m_CgemXY->Draw("");
        if (opt.Contains("ZR"))
        {
            m_CgemZR[0]->Draw("");
            m_CgemZR[1]->Draw("");
        }
 
        if (view->GetVisCgemLayers())
        {
            if (opt.Contains("XY"))
            {
                for (int layer = 0; layer < m_kCgemLayer; ++layer)
                {
                    m_CgemAnodeXY[layer]->Draw("");
                    m_CgemCathodeXY[layer]->Draw("");
                    for (int foil = 0; foil < m_kCgemFoil; ++foil)
                        m_CgemFoilXY[layer][foil]->Draw("");
                }
            }
 
            if (opt.Contains("ZR"))
            {
                for (int layer = 0; layer < m_kCgemLayer; ++layer)
                {
                    m_CgemAnodeZR[layer][0]->Draw("");
                    m_CgemAnodeZR[layer][1]->Draw("");
                    m_CgemCathodeZR[layer][0]->Draw("");
                    m_CgemCathodeZR[layer][1]->Draw("");
                    for (int foil = 0; foil < m_kCgemFoil; ++foil)
                    {
                        m_CgemFoilZR[layer][foil][0]->Draw("");
                        m_CgemFoilZR[layer][foil][1]->Draw("");
                    }
                }
            }
        }
 
        if (opt.Contains("ZR"))     //In ZR View, draw X-strip first
        {
            if (view->GetVisCgemXStrips())
            {
                for (int layer = 0; layer < m_kCgemLayer; ++layer)
                {
                    for (int sheet = 0; sheet < 2; ++sheet)
                    {
                        if (layer == 0 && sheet == 1) continue;
                        for (int strip = 0; strip < NXStrip[layer][sheet]; ++strip)
                        {
                            m_CgemXStrip[layer][sheet][strip]->SetFired(false);
                            m_CgemXStrip[layer][sheet][strip]->Draw(opt);
                        }
                    }
                }
            }
            if (view->GetVisCgemVStrips())
            {
                for (int layer = 0; layer < m_kCgemLayer; ++layer)
                {
                    for (int sheet = 0; sheet < 2; ++sheet)
                    {
                        if (layer == 0 && sheet == 1) continue;
                        for (int strip = 0; strip < NVStrip[layer][sheet]; ++strip)
                        {
                            m_CgemVStrip[layer][sheet][strip]->SetFired(false);
                            m_CgemVStrip[layer][sheet][strip]->Draw(opt);
                        }
                    }
                }
            }
        }
        else    //In XY View, draw V-strip first
        {
            if (view->GetVisCgemVStrips())
            {
                for (int layer = 0; layer < m_kCgemLayer; ++layer)
                {
                    for (int sheet = 0; sheet < 2; ++sheet)
                    {
                        if (layer == 0 && sheet == 1) continue;
                        for (int strip = 0; strip < NVStrip[layer][sheet]; ++strip)
                        {
                            m_CgemVStrip[layer][sheet][strip]->SetFired(false);
                            m_CgemVStrip[layer][sheet][strip]->Draw(opt);
                        }
                    }
                }
            }
            if (view->GetVisCgemXStrips())
            {
                for (int layer = 0; layer < m_kCgemLayer; ++layer)
                {
                    for (int sheet = 0; sheet < 2; ++sheet)
                    {
                        if (layer == 0 && sheet == 1) continue;
                        for (int strip = 0; strip < NXStrip[layer][sheet]; ++strip)
                        {
                            m_CgemXStrip[layer][sheet][strip]->SetFired(false);
                            m_CgemXStrip[layer][sheet][strip]->Draw(opt);
                        }
                    }
                }
            }
        }
    }
}

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

◆ Draw3DHits()

void CgemROOTGeo::Draw3DHits ( Option_t * option )

Definition at line 802 of file CgemROOTGeo.cxx.

{
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
    if (!view) cout << "CgemROOTGeo::Draw3DHits(), BesView not found" << endl;
 
    if (view->GetVisCgemHitsGlobal())
    {
        int NDigiCol = 0;
        if (gEvent)
        {
            m_CgemDigiCol = gEvent->GetCgemDigiCol();
            if (m_CgemDigiCol) NDigiCol = m_CgemDigiCol->GetEntries();
            else NDigiCol = 0;
        }
 
        for (int i = 0; i < NDigiCol; ++i)
        {
            TCgemDigi* aCgemDigi = (TCgemDigi*)m_CgemDigiCol->At(i);
            Identifier aCgemID(aCgemDigi->getIntId());
            int layer = CgemID::layer(aCgemID);
            int sheet = CgemID::sheet(aCgemID);
            int strip = CgemID::strip(aCgemID);
            bool is_x = CgemID::is_xstrip(aCgemID);
            TPolyLine3D* pStrip = 0;
 
            if (is_x)
            {
                if (view->GetVisCgemHitsXStrip() && strip < NXStrip[layer][sheet] && strip >= 0)
                {
                    //Construct X-strip 3D hits
                    Double_t P[6];
                    double phi = XStripPhi[layer][sheet][strip];
                    double r = RX[layer][sheet];
                    double zmin = ZMin[layer][sheet];
                    P[0] = P[3] = TMath::Cos(phi) * r;
                    P[1] = P[4] = TMath::Sin(phi) * r;
                    P[2] = -zmin;
                    P[5] = zmin;
 
                    pStrip = new TPolyLine3D(2, P);
                    pStrip->SetBit(kCanDelete);
                    pStrip->SetLineWidth(4);
                    pStrip->SetLineColor(m_XstripColor);
                    pStrip->Draw(option);
                }
            }
            else
            {
                if (view->GetVisCgemHitsVStrip() && strip < NVStrip[layer][sheet] && strip >= 0)
                {
                    //Construct V-strip 3D hits
                    double r = RV[layer][sheet];
                    double len = SheetLength[layer][sheet];
                    double wid = SheetWidth[layer][sheet];
                    double ang = StereoAngle[layer][sheet];
                    double xmin = XMin[layer][sheet];
                    double zmin = ZMin[layer][sheet];
                    double pit = VPitch[layer][sheet];
                    CgemVStripCalculator calc(r, len, wid, ang, xmin, zmin, pit);
 
                    int N = calc.getPointNumberFromV(strip);
                    if (N > 0)
                    {
                        pStrip = new TPolyLine3D(N);
                        pStrip->SetBit(kCanDelete);
                        pStrip->SetLineWidth(3);
                        pStrip->SetLineColor(m_VstripColor);
                        
                        double tphi = 0, tz = 0;
                        for (int i = 0; i < N; ++i)
                        {
                            //New Geom
                            double v = VStripV[layer][sheet][strip];
                            calc.getNthPointFromV(v, i, tphi, tz);
                            double x = r * TMath::Cos(tphi);
                            double y = r * TMath::Sin(tphi);
                            double z = tz;
                            pStrip->SetPoint(i, x, y, z);
                        }
                        
                        pStrip->Draw(option);
                    }
                }
            }
        }
    }
}

Referenced by BesEvent::DrawCgem3DHits().

◆ DrawHits()

void CgemROOTGeo::DrawHits ( Option_t * option )

Definition at line 752 of file CgemROOTGeo.cxx.

{
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
    if (!view) cout << "CgemROOTGeo::DrawHits(), BesView not found" << endl;
 
    if (view->GetVisCgemHitsGlobal())
    {
        for (int i = 0; i < m_2DHitsArray->GetEntries(); ++i)
        {
            Cgem2DStripBase* aStrip = (Cgem2DStripBase*)m_2DHitsArray->At(i);
            bool skip = false;
            if (aStrip->GetType() == CGEM_X_STRIP)
            {
                if (!view->GetVisCgemHitsXStrip()) skip = true;
            }
            else if (aStrip->GetType() == CGEM_V_STRIP)
            {
                if (!view->GetVisCgemHitsVStrip()) skip = true;
            }
 
            if (!skip && aStrip)
            {
                //normal view
                if ((view->GetViewType() & kXYView) || (view->GetViewType() & kZRView))
                {
                    aStrip->SetFired(true);
                    aStrip->Draw(option);
                }
                
                //unfolded view
                else if ((view->GetViewType() & kCgemUFView0) && aStrip->GetLayer() == 0)
                {
                    aStrip->SetFired(true);
                    aStrip->Draw(option);
                }
                else if ((view->GetViewType() & kCgemUFView1) && aStrip->GetLayer() == 1)
                {
                    aStrip->SetFired(true);
                    aStrip->Draw(option);
                }
                else if ((view->GetViewType() & kCgemUFView2) && aStrip->GetLayer() == 2)
                {
                    aStrip->SetFired(true);
                    aStrip->Draw(option);
                }
            }
        }   //for end
    }
}

Referenced by BesEvent::DrawCgemUFHits(), and BesEvent::DrawHits().

◆ DrawUF()

void CgemROOTGeo::DrawUF	(	int	layer,
		Option_t *	option = "" )

Definition at line 714 of file CgemROOTGeo.cxx.

{
    if (layer < 0 || layer >= m_kCgemLayer) return;
    TString opt = option;
    opt += ",UF";
    opt.ToUpper();
    if (!m_2DGeoInit) cout << "CgemROOTGeo::Draw(), Geometry not initialized!" << endl;
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
    if (!view) cout << "CgemROOTGeo::Draw(), BesView not found" << endl;
 
    if (view->GetVisCgemGlobal())
    {
        for (int sheet = 0; sheet < 2; ++sheet)
        {
            if (layer == 0 && sheet == 1) continue;
            if (view->GetVisCgemLayers()) m_CgemSheet[layer][sheet]->Draw(opt);
 
            if (view->GetVisCgemXStrips())
            {
                for (int strip = 0; strip < NXStrip[layer][sheet]; ++strip)
                {
                    m_CgemXStrip[layer][sheet][strip]->SetFired(false);
                    m_CgemXStrip[layer][sheet][strip]->Draw(opt);
                }
            }
            if (view->GetVisCgemVStrips())
            {
                for (int strip = 0; strip < NVStrip[layer][sheet]; ++strip)
                {
                    m_CgemVStrip[layer][sheet][strip]->SetFired(false);
                    m_CgemVStrip[layer][sheet][strip]->Draw(opt);
                }
            }
        }
    }
}

Referenced by BesGeometry::DrawCgemUF().

◆ Get2DStrip()

Cgem2DStripBase * CgemROOTGeo::Get2DStrip	(	bool	is_x,
		int	layer,
		int	sheet,
		int	strip ) const

Definition at line 891 of file CgemROOTGeo.cxx.

{
    if (is_x)
    {
        if (strip >= 0 && strip < NXStrip[layer][sheet])
            return m_CgemXStrip[layer][sheet][strip];
    }
    else
    {
        if (strip >= 0 && strip < NVStrip[layer][sheet])
            return m_CgemVStrip[layer][sheet][strip];
    }
    return 0;
}

◆ GetCgemAnode()

TGeoNode * CgemROOTGeo::GetCgemAnode ( int layer )

Definition at line 580 of file CgemROOTGeo.cxx.

{
    return m_NodeCgemAnode[layer];
}

Referenced by InitGeometry().

◆ GetCgemCathode()

TGeoNode * CgemROOTGeo::GetCgemCathode ( int layer )

Definition at line 585 of file CgemROOTGeo.cxx.

{
    return m_NodeCgemCathode[layer];
}

Referenced by InitGeometry().

◆ GetCgemFoil()

TGeoNode * CgemROOTGeo::GetCgemFoil	(	int	layer,
		int	foil )

Definition at line 590 of file CgemROOTGeo.cxx.

{
    return m_NodeCgemFoil[layer][foil];
}

Referenced by InitGeometry().

◆ getRFromLayerSheet()

double CgemROOTGeo::getRFromLayerSheet	(	int	layer,
		int	sheet ) const

Definition at line 922 of file CgemROOTGeo.cxx.

{
    if (layer < 0 || layer >= m_kCgemLayer || sheet < 0 || sheet >= 2) return 0.0;
    return (RX[layer][sheet] + RV[layer][sheet]) / 2;
}

Referenced by BesEvent::SetCgemClusters().

◆ getSheetWidth()

double CgemROOTGeo::getSheetWidth	(	int	layer,
		int	sheet ) const

Definition at line 928 of file CgemROOTGeo.cxx.

{
    if (layer < 0 || layer >= m_kCgemLayer || sheet < 0 || sheet >= 2) return 0.0;
    return SheetWidth[layer][sheet];
}

Referenced by BesEvent::SetCgemClusters().

◆ GetVolumeCgem()

TGeoVolume * CgemROOTGeo::GetVolumeCgem ( ) const

inline

Definition at line 70 of file CgemROOTGeo.h.

70{ return m_CGEM; }

Referenced by BesGeometry::GetVolumeCgem(), BesGeometry::InitFromGDML(), and InitGeometry().

◆ InitFromGDML()

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

Definition at line 81 of file CgemROOTGeo.cxx.

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

Referenced by BesGeometry::InitFromGDML().

◆ InitFromROOT()

void CgemROOTGeo::InitFromROOT ( TGeoVolume * vol )

Definition at line 89 of file CgemROOTGeo.cxx.

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

Referenced by BesGeometry::InitFromROOT().

◆ InitGeometry()

void CgemROOTGeo::InitGeometry ( )

Definition at line 97 of file CgemROOTGeo.cxx.

{
    Int_t CgemLineColor = 15;
    Int_t CgemFillColor = 18;
    Int_t CgemFillStyle = 1001;
    Int_t CgemLayerStyle = 3003;
    Int_t CgemElectrodeColor = 12;
    Int_t CgemFoilColor = 12;
    m_2DGeoInit = 1;
 
    TGeoVolume* CgemVol;
    TGeoTube* CgemTube;
    TString name;
    Double_t center[3] = {0.0, 0.0, 0.0};
    Double_t P[12];
 
    //XY View
    TGeoTube* CgemShape = (TGeoTube*)GetVolumeCgem()->GetShape();
    m_CgemXY = new BesCircle2D("CGEM", "CGEM", CgemShape->GetRmin(), CgemShape->GetRmax(), &center[0]);
    m_CgemXY->SetNSegment(360);
    m_CgemXY->SetFillColor(CgemFillColor);
    m_CgemXY->SetFillStyle(CgemFillStyle);
    m_CgemXY->SetLineColor(CgemLineColor);
 
    for (int layer = 0; layer < m_kCgemLayer; ++layer)
    {
        CgemVol = GetCgemAnode(layer)->GetVolume();
        CgemTube = (TGeoTube*)CgemVol->GetShape();
        name = "CGEM Layer";
        name += layer;
        name += " Anode";
        center[2] = CgemTube->GetDz();
        m_CgemAnodeXY[layer] = new BesCircle2D(name, name, CgemTube->GetRmin(), CgemTube->GetRmax(), center);
        m_CgemAnodeXY[layer]->SetNSegment(360);
        m_CgemAnodeXY[layer]->SetFillStyle(CgemLayerStyle);
        m_CgemAnodeXY[layer]->SetLineColor(CgemElectrodeColor);
        m_CgemAnodeXY[layer]->SetFillColor(CgemElectrodeColor);
 
        CgemVol = GetCgemCathode(layer)->GetVolume();
        CgemTube = (TGeoTube*)CgemVol->GetShape();
        name = "CGEM Layer";
        name += layer;
        name += " Cathode";
        center[2] = CgemTube->GetDz();
        m_CgemCathodeXY[layer] = new BesCircle2D(name, name, CgemTube->GetRmin(), CgemTube->GetRmax(), center);
        m_CgemCathodeXY[layer]->SetNSegment(360);
        m_CgemCathodeXY[layer]->SetFillStyle(CgemLayerStyle);
        m_CgemCathodeXY[layer]->SetLineColor(CgemElectrodeColor);
        m_CgemCathodeXY[layer]->SetFillColor(CgemElectrodeColor);
        
        for (int foil = 0; foil < m_kCgemFoil; ++foil)
        {
            CgemVol = GetCgemFoil(layer, foil)->GetVolume();
            CgemTube = (TGeoTube*)CgemVol->GetShape();
            name = "CGEM Layer";
            name += layer;
            name += " Foil";
            name += foil;
            center[2] = CgemTube->GetDz();
            m_CgemFoilXY[layer][foil] = new BesCircle2D(name, name, CgemTube->GetRmin(), CgemTube->GetRmax(), center);
            m_CgemFoilXY[layer][foil]->SetNSegment(360);
            m_CgemFoilXY[layer][foil]->SetFillStyle(CgemLayerStyle);
            m_CgemFoilXY[layer][foil]->SetLineColor(CgemFoilColor);
            m_CgemFoilXY[layer][foil]->SetFillColor(CgemFoilColor);
        }
    }
 
    //ZR View
    for (int j = 0; j < 2; ++j)
    {
        P[0] = P[3] = P[6] = P[9] = 0;
        P[1] = P[10] = CgemShape->GetRmin() * (1 - 2 * j);
        P[2] = P[5] = -CgemShape->GetDz();
        P[4] = P[7] = CgemShape->GetRmax() * (1 - 2 * j);
        P[8] = P[11] = CgemShape->GetDz();
 
        name = "CgemZR ";
        if (j == 0) name += "Up";
        else name += "Down";
        m_CgemZR[j] = new BesPolygon2D(name, name, 4, P);
        m_CgemZR[j]->SetFillColor(CgemFillColor);
        m_CgemZR[j]->SetFillStyle(CgemFillStyle);
        m_CgemZR[j]->SetLineColor(CgemLineColor);
        m_CgemZR[j]->SetRotatable(true);
    }
 
    for (int layer = 0; layer < m_kCgemLayer; ++layer)
    {
        for (int j = 0; j < 2; ++j)
        {
            CgemVol = GetCgemAnode(layer)->GetVolume();
            CgemTube = (TGeoTube*)CgemVol->GetShape();
            name = "CgemZR Layer";
            name += layer;
            if (j == 0) name += " Anode Up";
            else name += " Anode Down";
 
            P[0] = P[3] = P[6] = P[9] = 0;
            P[1] = P[10] = CgemTube->GetRmin() * (1 - 2 * j);
            P[2] = P[5] = -CgemTube->GetDz();
            P[4] = P[7] = CgemTube->GetRmax() * (1 - 2 * j);
            P[8] = P[11] = CgemTube->GetDz();
 
            m_CgemAnodeZR[layer][j] = new BesPolygon2D(name, name, 4, P);
            m_CgemAnodeZR[layer][j]->SetFillStyle(CgemLayerStyle);
            m_CgemAnodeZR[layer][j]->SetLineColor(CgemElectrodeColor);
            m_CgemAnodeZR[layer][j]->SetFillColor(CgemElectrodeColor);
            m_CgemAnodeZR[layer][j]->SetRotatable(true);
 
            CgemVol = GetCgemCathode(layer)->GetVolume();
            CgemTube = (TGeoTube*)CgemVol->GetShape();
            name = "CgemZR Layer";
            name += layer;
            if (j == 0) name += " Cathode Up";
            else name += " Cathode Down";
 
            P[0] = P[3] = P[6] = P[9] = 0;
            P[1] = P[10] = CgemTube->GetRmin() * (1 - 2 * j);
            P[2] = P[5] = -CgemTube->GetDz();
            P[4] = P[7] = CgemTube->GetRmax() * (1 - 2 * j);
            P[8] = P[11] = CgemTube->GetDz();
 
            m_CgemCathodeZR[layer][j] = new BesPolygon2D(name, name, 4, P);
            m_CgemCathodeZR[layer][j]->SetFillStyle(CgemLayerStyle);
            m_CgemCathodeZR[layer][j]->SetLineColor(CgemElectrodeColor);
            m_CgemCathodeZR[layer][j]->SetFillColor(CgemElectrodeColor);
            m_CgemCathodeZR[layer][j]->SetRotatable(true);
 
            for (int foil = 0; foil < m_kCgemFoil; ++foil)
            {
                CgemVol = GetCgemFoil(layer, foil)->GetVolume();
                CgemTube = (TGeoTube*)CgemVol->GetShape();
                name = "CgemZR Layer";
                name += layer;
                name += " Foil";
                name += foil;
                if (j == 0) name += " Up";
                else name += " Down";
 
                P[0] = P[3] = P[6] = P[9] = 0;
                P[1] = P[10] = CgemTube->GetRmin() * (1 - 2 * j);
                P[2] = P[5] = -CgemTube->GetDz();
                P[4] = P[7] = CgemTube->GetRmax() * (1 - 2 * j);
                P[8] = P[11] = CgemTube->GetDz();
 
                m_CgemFoilZR[layer][foil][j] = new BesPolygon2D(name, name, 4, P);
                m_CgemFoilZR[layer][foil][j]->SetFillStyle(CgemLayerStyle);
                m_CgemFoilZR[layer][foil][j]->SetLineColor(CgemFoilColor);
                m_CgemFoilZR[layer][foil][j]->SetFillColor(CgemFoilColor);
                m_CgemFoilZR[layer][foil][j]->SetRotatable(true);
            }
        }
    }
 
    //Calculate CGEM strips by using CgemGeomSvc package
    CgemGeomSvc cgs("CgemGeomSvc", Gaudi::svcLocator());
    //cgs.initGeom();
    cgs.initialize();
    CgemGeoReadoutPlane* plane = 0;
    CgemGeoLayer* geolayer = 0;
 
    for (int layer = 0; layer < 3; ++layer)
    {
        for (int sheet = 0; sheet < 2; ++sheet)
        {
            //Load data of Readout Plane
            if (layer == 0 && sheet == 1) continue;
            plane = cgs.getReadoutPlane(layer, sheet);
            geolayer = cgs.getCgemLayer(layer);
 
            cout << "layer" << layer << " sheet" << sheet << " " << plane << endl;
            cout << "\tRX:" << plane->getRX() << endl;
            cout << "\tRV:" << plane->getRV() << endl;
            cout << "\tXmin:" << plane->getXmin() << endl;
            cout << "\tWidth:" << plane->getWidth() << endl;
            cout << "\tZmin:" << plane->getZmin() << endl;
            cout << "\tLength:" << plane->getLength() << endl;
            cout << "\tXPitch:" << plane->getXPitch() << endl;
            cout << "\tVPitch:" << plane->getVPitch() << endl;
            cout << "\tStereo Angle:" << plane->getStereoAngle() << endl;
            cout << "\tMid R at Gap:" << plane->getMidRAtGap() << endl;
            cout << "\tVmax:" << plane->getVmax() << endl;
            cout << "\tX-strip number:" << plane->getNXstrips() << endl;
            cout << "\tV-strip number:" << plane->getNVstrips() << endl;
            cout << "\tX-strip width:" << geolayer->getWidthOfStripX() << endl;
            cout << "\tV-strip width:" << geolayer->getWidthOfStripV() << endl;
 
            //Construct sheets unfolded geometry
            //P[1] = P[4] = -plane->getLength() / 2;
            P[1] = P[4] = plane->getZmin();    //New Geom
            P[2] = P[5] = P[8] = P[11] = 0;
            //P[7] = P[10] = plane->getLength() / 2;
            P[7] = P[10] = plane->getZmin() + plane->getLength();    //New Geom
            P[0] = P[9] = plane->getXmin();     //New Geom
            P[3] = P[6] = plane->getXmin() + plane->getWidth();     //New Geom
            /*
            if (layer == 0)
            {
                P[0] = P[9] = -plane->getWidth() / 2;
                P[3] = P[6] = plane->getWidth() / 2;
            }
            else
            {
                P[0] = P[9] = plane->getXmin() - plane->getWidth();
                P[3] = P[6] = plane->getXmin();
            }
            */
            
            name = "CGEM Layer";
            name += layer;
            name += " Sheet";
            name += sheet;
            m_CgemSheet[layer][sheet] = new BesPolygon2D(name, name, 4, P);
            m_CgemSheet[layer][sheet]->SetFillColor(CgemFillColor);
            m_CgemSheet[layer][sheet]->SetFillStyle(CgemFillStyle);
            m_CgemSheet[layer][sheet]->SetLineColor(CgemLineColor);
 
 
            //Variables which will be used in constructing strips geometry
            double r, wid, len, phi, xmin, zmin, ang, pit, thick, size, v, ruf;
            ruf = (plane->getRX() + plane->getRV()) / 2;
 
            //Construct X-strips 2D and 3D geometry
            int n = plane->getNXstrips();
            NXStrip[layer][sheet] = n;
            m_CgemXStrip[layer][sheet] = new Cgem2DXStrip*[n];
            XStripPhi[layer][sheet] = new Double_t[n];
            len = plane->getLength();
            wid = plane->getWidth();
            zmin = plane->getZmin();
            r = plane->getRX();
            thick = 0.1;
            size = geolayer->getWidthOfStripX();
            SheetWidth[layer][sheet] = wid;
            RX[layer][sheet] = r;
            ZMin[layer][sheet] = zmin;
 
            for (int strip = 0; strip < n; ++strip)
            {
                name = "CGEM X-Strip Layer";
                name += layer;
                name += " Sheet";
                name += sheet;
                name += " Strip";
                name += strip;
                phi = plane->getPhiFromXID(strip);
                XStripPhi[layer][sheet][strip] = phi;
                m_CgemXStrip[layer][sheet][strip] = new Cgem2DXStrip(name, name,
                    len, wid, phi, zmin, r, size, thick, layer, sheet, strip, ruf);
            }
 
            phi = plane->getPhiFromXID(n - 1);
            phi -= plane->getPhiFromXID(0);
 
            //Construct V-strips 2D and 3D geometry
            n = plane->getNVstrips();
            NVStrip[layer][sheet] = n;
            m_CgemVStrip[layer][sheet] = new Cgem2DVStrip*[n];
            VStripV[layer][sheet] = new Double_t[n];
            r = plane->getRV();
            len = plane->getLength();
            wid = phi * ruf + geolayer->getWidthOfStripX();     //New Geom
            ang = plane->getStereoAngle();
            xmin = plane->getXmin();
            pit = plane->getVPitch();
            thick = 0.1;
            size = geolayer->getWidthOfStripV();
            RV[layer][sheet] = r;
            SheetLength[layer][sheet] = len;
            StereoAngle[layer][sheet] = ang;
            //New Geom
            XMin[layer][sheet] = plane->getPhiFromXID(0) * ruf - geolayer->getWidthOfStripX() / 2;  //New Geom
            xmin = XMin[layer][sheet];
            VPitch[layer][sheet] = pit;
            XPitch[layer][sheet] = plane->getXPitch();
 
            for (int strip = 0; strip < n; ++strip)
            {
                name = "CGEM V-Strip Layer";
                name += layer;
                name += " Sheet";
                name += sheet;
                name += " Strip";
                name += strip;
                //New Geom
                v = plane->getCentralVFromVID(strip);
                VStripV[layer][sheet][strip] = v;
                m_CgemVStrip[layer][sheet][strip] = new Cgem2DVStrip(name, name, v,
                    r, len, wid, ang, xmin, zmin, pit, size, thick, layer, sheet, strip, ruf);
            }
        }   //sheet
    }   //layer
}

Referenced by BesGeometry::InitGeometry().

◆ isPointInSheet()

double CgemROOTGeo::isPointInSheet	(	int	layer,
		int	sheet,
		double	phi,
		double	z ) const

Definition at line 934 of file CgemROOTGeo.cxx.

{
    if (layer < 0 || layer >= m_kCgemLayer || sheet < 0 || sheet >= 2) return 0.0;
    double ruf = (RV[layer][sheet] + RX[layer][sheet]) / 2;
    double len = SheetLength[layer][sheet];
    double wid = SheetWidth[layer][sheet];
    double ang = StereoAngle[layer][sheet];
    double xmin = XMin[layer][sheet];
    double zmin = ZMin[layer][sheet];
    double pit = VPitch[layer][sheet];
    CgemVStripCalculator calc(ruf, len, wid, ang, xmin, zmin, pit);
    return calc.isPointInSheet(phi, z);
}

◆ SetAllVisible()

void CgemROOTGeo::SetAllVisible ( )

Definition at line 449 of file CgemROOTGeo.cxx.

{
 
}

◆ SetDetector()

void CgemROOTGeo::SetDetector ( )

Definition at line 473 of file CgemROOTGeo.cxx.

{
 
}

Referenced by BesClient::HandleViewOptionMenu().

◆ SetHalfVisible()

void CgemROOTGeo::SetHalfVisible ( )

Definition at line 461 of file CgemROOTGeo.cxx.

{
 
}

◆ SetHits()

void CgemROOTGeo::SetHits ( )

Definition at line 492 of file CgemROOTGeo.cxx.

{
    // set previous event 2D hits info to default
    ClearHits();
 
    int NDigiCol = 0;
    if (gEvent)
    {
        m_CgemDigiCol = gEvent->GetCgemDigiCol();
        if (m_CgemDigiCol) NDigiCol = m_CgemDigiCol->GetEntries();
        else NDigiCol = 0;
    }
 
    for (int i = 0; i < NDigiCol; ++i)
    {
        TCgemDigi* aCgemDigi = (TCgemDigi*)m_CgemDigiCol->At(i);
        Identifier aCgemID(aCgemDigi->getIntId());
        int layer = CgemID::layer(aCgemID);
        int sheet = CgemID::sheet(aCgemID);
        int strip = CgemID::strip(aCgemID);
        bool is_x = CgemID::is_xstrip(aCgemID);
 
        //Double_t time   = RawDataUtil::CgemTime(aCgemDigi->getTimeChannel());
        //Int_t charge = RawDataUtil::CgemCharge(aCgemDigi->getChargeChannel());
    
    Double_t charge = aCgemDigi->getCharge_fc();
    //cout<<"charge="<<charge<<endl;
    Double_t time   = aCgemDigi->getTime_ns();
 
        //Construct 2D hits
        Cgem2DStripBase* aStrip = 0;
        if (is_x)
        {
            if (strip < NXStrip[layer][sheet] && strip >= 0)
                aStrip = m_CgemXStrip[layer][sheet][strip];
        }
        else
        {
            if (strip < NVStrip[layer][sheet] && strip > 0)
                aStrip = m_CgemVStrip[layer][sheet][strip];
        }
 
        //Add 2D hits
        if (aStrip)
        {
            aStrip->CloseInfo();
            //aStrip->AddInfo(aStrip->GetTitle());
            double evTime = gEvent->GetHeader().GetEvTime();
 
            char data[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=%-.2f fC", data,charge);
            
            aStrip->AddInfo(TString(data));
            sprintf(data, "Fired");
            aStrip->AddInfo(TString(data));
            //aStrip->SetTime(time);
            //aStrip->SetCharge(charge);
            //unsigned int overflow = aCgemDigi->getOverflow();
            //aWire->SetQOverflow(overflow & 2);
            //aWire->SetEvTime(evTime);
            aStrip->CloseInfo();
            m_2DHitsArray->Add(aStrip);
        }
    }
}

Referenced by BesEvent::SetHits().

◆ SetNode()

void CgemROOTGeo::SetNode ( )

Definition at line 392 of file CgemROOTGeo.cxx.

{
    if (m_ROOTGeoInit == 2)
    {   //from GDML
        m_CGEM = GetTopVolume();
        if (!m_CGEM) std::cout << "m_CGEM = 0" << std::endl;
        // Int_t nDaughters = m_CGEM->GetNodes()->GetEntries();
        // cout << "CGEM contains " << nDaughters  << " nodes : " << endl;
        // analNodes(m_CGEM->GetNodes(), 0);
 
        for (int layer = 0; layer < m_kCgemLayer; ++layer)
        {
            std::stringstream osanode, oscathode;
            osanode << "pv_Anode_logic" << layer << "_0";
            oscathode << "pv_Cathode_logic" << layer << "_8";
            m_NodeCgemAnode[layer] = GetNode(osanode.str());
            m_NodeCgemCathode[layer] = GetNode(oscathode.str());
            // cout << osanode.str() << "   " << m_NodeCgemAnode[layer] << endl;
            // cout << oscathode.str() << "   " << m_NodeCgemCathode[layer] << endl;
 
            for (int foil = 0; foil < m_kCgemFoil; ++foil)
            {
                std::stringstream osfoil;
                osfoil << "pv_GemFoil_logic" << layer << "foil" << foil << "_" << m_kCgemFoil - foil;
                m_NodeCgemFoil[layer][foil] = GetNode(osfoil.str());
                // cout << osfoil.str() << "   " << m_NodeCgemFoil[layer][foil] << endl;
            }
        }
    }
    else if (m_ROOTGeoInit == 1)
    {   //from ROOT
        for (int layer = 0; layer < m_kCgemLayer; ++layer)
        {
            TGeoNode* cgemLayer = m_CGEM->GetNode(m_kCgemLayer - layer);
            m_NodeCgemAnode[layer] = cgemLayer->GetVolume()->GetNode(0);
            m_NodeCgemCathode[layer] = cgemLayer->GetVolume()->GetNode(8);
            for (int foil = 0; foil < m_kCgemFoil; ++foil)
            {
                m_NodeCgemFoil[layer][foil] = cgemLayer->GetVolume()->GetNode(m_kCgemFoil - foil);
            }
        }
    }
}

Referenced by InitFromGDML(), and InitFromROOT().

◆ SetPhysicalDefaultVis()

void CgemROOTGeo::SetPhysicalDefaultVis ( )

Definition at line 443 of file CgemROOTGeo.cxx.

{
    
}

Referenced by BesGeometry::SetPhysicalDefaultVis().

◆ SetPhysicalNode()

void CgemROOTGeo::SetPhysicalNode ( )

Definition at line 467 of file CgemROOTGeo.cxx.

{
 
}

Referenced by BesGeometry::InitGeometry().

◆ SetQuarterVisible()

void CgemROOTGeo::SetQuarterVisible ( )

Definition at line 455 of file CgemROOTGeo.cxx.

{
 
}

◆ SetVisCgemDetector()

void CgemROOTGeo::SetVisCgemDetector ( )

Definition at line 569 of file CgemROOTGeo.cxx.

{
 
}

Referenced by BesClient::HandleViewOptionMenu().

◆ SetVisCgemHits()

void CgemROOTGeo::SetVisCgemHits ( )

Definition at line 575 of file CgemROOTGeo.cxx.

{
 
}

◆ SetVolumeCgem()

void CgemROOTGeo::SetVolumeCgem ( TGeoVolume * vol )

inline

Definition at line 68 of file CgemROOTGeo.h.

68{ m_CGEM = vol; }

Referenced by InitFromROOT().

◆ SetVolumeDefaultVis()

void CgemROOTGeo::SetVolumeDefaultVis ( )

Definition at line 437 of file CgemROOTGeo.cxx.

{
 
}

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CgemROOTGeo()

◆ ~CgemROOTGeo()

Member Function Documentation

◆ calcZFromLayerSheetVIDPhi()

◆ ClearHits()

◆ Draw()

◆ Draw3DHits()

◆ DrawHits()

◆ DrawUF()

◆ Get2DStrip()

◆ GetCgemAnode()

◆ GetCgemCathode()

◆ GetCgemFoil()

◆ getRFromLayerSheet()

◆ getSheetWidth()

◆ GetVolumeCgem()

◆ InitFromGDML()

◆ InitFromROOT()

◆ InitGeometry()

◆ isPointInSheet()

◆ SetAllVisible()

◆ SetDetector()

◆ SetHalfVisible()

◆ SetHits()

◆ SetNode()

◆ SetPhysicalDefaultVis()

◆ SetPhysicalNode()

◆ SetQuarterVisible()

◆ SetVisCgemDetector()

◆ SetVisCgemHits()

◆ SetVolumeCgem()

◆ SetVolumeDefaultVis()