---

VisCSCChamberDrawer.cc

Go to the documentation of this file.

//<<<<<< INCLUDES                                                       >>>>>>
#include "VisReco/VisMuonCSC/interface/VisCSCChamberDrawer.h"

#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "Geometry/CSCGeometry/interface/CSCLayer.h"
#include "Geometry/CSCGeometry/interface/CSCLayerGeometry.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include <Inventor/nodes/SoDrawStyle.h>
#include <Inventor/nodes/SoMaterial.h>
#include <Inventor/nodes/SoSeparator.h>
#include <Inventor/nodes/SoLineSet.h>
#include <Inventor/nodes/SoCoordinate3.h>

//<<<<<< PRIVATE VARIABLE DEFINITIONS                                   >>>>>>
 
VisCSCChamberDrawer::VisCSCChamberDrawer(void) : m_lineColor(0x59595900), m_lineThickness(2.0F), m_chamberFlag(true)
{
}

bool VisCSCChamberDrawer::getOutline(SoSeparator *sep, const edm::ESHandle<CSCGeometry>  pGeom, const CSCDetId id)
{
        if (!pGeom.isValid() || !m_chamberFlag) return false;
        
        SbColor lineColor;
        float transparency = 0.0f;
        lineColor.setPackedValue(static_cast<uint32_t>(m_lineColor), transparency);
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (lineColor); 
        SoDrawStyle *drawStyle = new SoDrawStyle;
        if (m_lineThickness != 0.F) drawStyle->lineWidth = m_lineThickness;
        
        
        sep->addChild(mat);
        sep->addChild(drawStyle);
        
        const CSCLayer *pLayer = pGeom->layer (id);
        if (!pLayer)
        {
            pLayer = pGeom->layer(4); // pick a default if not present
            //assert (pLayer);
        }
        const CSCLayerGeometry *layerGeom = pLayer->geometry ();
        // params: half bottom edge, half top edge, half thickness, half apothem
        const std::vector<float> trapezBounds = layerGeom->parameters ();
        const GeomDet *det = pGeom->idToDet(id);
        std::vector<Surface::GlobalPoint> crossPoint;
        // get vertices, converting from cm to meters
        crossPoint.push_back (det->surface ().toGlobal
                              (LocalPoint (-trapezBounds [0], -trapezBounds [3], 0.0)));
        crossPoint.push_back (det->surface ().toGlobal
                              (LocalPoint ( trapezBounds [0], -trapezBounds [3], 0.0)));
        crossPoint.push_back (det->surface ().toGlobal
                              (LocalPoint ( trapezBounds [1],  trapezBounds [3], 0.0)));
        crossPoint.push_back (det->surface ().toGlobal
                              (LocalPoint (-trapezBounds [1],  trapezBounds [3], 0.0)));
        crossPoint.push_back(crossPoint.front()); // close the loop
        std::vector<SbVec3f> pts(5);
        for (int i=0; i < 5; i++)
        {
                pts[i] = SbVec3f(crossPoint[i].x(), crossPoint[i].y(),  crossPoint[i].z())*0.01F;
        }             
        SoCoordinate3 *vert = new SoCoordinate3;
        vert->point.setValues(0, 5, &pts[0]);
        
        sep->addChild(vert);
        SoLineSet *trapez = new SoLineSet;
        sep->addChild (trapez);

    return true;
}

---