---

VisCSCWireDigiTwig.cc

Go to the documentation of this file.

//<<<<<< INCLUDES                                                       >>>>>>

#include "VisReco/VisMuonCSC/interface/VisCSCWireDigiTwig.h"
#include "VisReco/VisMuonCSC/interface/VisCSCDigiTools.h"
#include "VisReco/VisMuonCSC/interface/VisCSCChamberDrawer.h"
#include "VisFramework/VisFrameworkBase/interface/VisTwigFactroyService.h"
#include "VisFramework/VisFrameworkBase/interface/VisEventSelector.h"
#include "VisFramework/VisFrameworkBase/interface/debug.h"
#include "DataFormats/CSCDigi/interface/CSCWireDigi.h"
#include "DataFormats/CSCDigi/interface/CSCWireDigiCollection.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.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 "Iguana/Models/interface/IgTextRep.h"
#include "Iguana/GLModels/interface/Ig3DRep.h"
#include "Iguana/GLModels/interface/IgLegoRep.h"
#include "Iguana/GLModels/interface/IgRPhiRep.h"
#include "Iguana/GLModels/interface/IgRZRep.h"
#include "Iguana/GLModels/interface/IgSpareRep.h"
#include "Iguana/Inventor/interface/IgSoTower.h"
#include "Iguana/Studio/interface/IgQtLock.h"
#include <Inventor/nodes/SoDrawStyle.h>
#include <Inventor/nodes/SoMaterial.h>
#include <Inventor/nodes/SoSeparator.h>
#include <Inventor/nodes/SoCube.h>
#include <Inventor/nodes/SoLineSet.h>
#include <Inventor/nodes/SoTranslation.h>
#include <Inventor/nodes/SoPointSet.h>
#include <Inventor/nodes/SoTransform.h>
#include <Inventor/nodes/SoRotation.h>
#include <Inventor/nodes/SoText3.h>
#include <Inventor/nodes/SoFont.h>
#include <classlib/utils/DebugAids.h>
#include <qstring.h>
#include <iostream>
#include <iomanip>

//<<<<<< PRIVATE DEFINES                                                >>>>>>
//<<<<<< PRIVATE CONSTANTS                                              >>>>>>
//<<<<<< PRIVATE TYPES                                                  >>>>>>
//<<<<<< PRIVATE VARIABLE DEFINITIONS                                   >>>>>>
//<<<<<< PUBLIC VARIABLE DEFINITIONS                                    >>>>>>
//<<<<<< CLASS STRUCTURE INITIALIZATION                                 >>>>>>
//<<<<<< PRIVATE FUNCTION DEFINITIONS                                   >>>>>>

namespace 
{
    VisQueuedTwig *
    createThisTwig (IgState *state, IgTwig *parent,
                    const std::string &name,
                    const std::string &friendlyName,
                    const std::string &modLabel,
                    const std::string &instanceName,
                    const std::string &processName)
    {
        return new VisCSCWireDigiTwig (state, parent, "Wire digis (" + name + ")", 
                                       friendlyName, modLabel, instanceName, processName);
    }
}

//<<<<<< PUBLIC FUNCTION DEFINITIONS                                    >>>>>>
//<<<<<< MEMBER FUNCTION DEFINITIONS                                    >>>>>>

VisCSCWireDigiTwig::VisCSCWireDigiTwig (IgState *state, IgTwig *parent,
                                        const std::string &name /* = "" */,
                                        const std::string &friendlyName /* = "" */,
                                        const std::string &moduleLabel /* = "" */,
                                        const std::string &instanceName /* = "" */,
                                        const std::string &processName /* = "" */)
    : VisQueuedTwig (state, parent, name),
      m_text (name),
      m_friendlyName (friendlyName),
      m_moduleLabel (moduleLabel),
      m_instanceName (instanceName),
      m_processName (processName)
{
    VisTwigFactroyService *tfService = VisTwigFactroyService::get (state);
    if (! tfService)
    {
        tfService = new VisTwigFactroyService (state);
    }   
    edm::TypeID caloJetsID (typeid (CSCWireDigiCollection));
    tfService->registerTwig (caloJetsID.friendlyClassName (), &createThisTwig);
}

void
VisCSCWireDigiTwig::onNewEvent (const edm::Event &event,
                                const edm::EventSetup &eventSetup)
{
    // Get debugging dump.
    VisQueuedTwig::onNewEvent (event, eventSetup);

    m_digis.clear ();
    m_text = (QString ("Run # %1, event # %2")
              .arg (event.id ().run ())
              .arg (event.id ().event ()).latin1 ());

    std::vector<edm::Handle<CSCWireDigiCollection> > digiCollections;

    try
    {
        if ((! m_friendlyName.empty ()) || (! m_moduleLabel.empty ()) || (! m_instanceName.empty ()) || (! m_processName.empty ()))
        {
            VisEventSelector visSel (m_friendlyName, m_moduleLabel, m_instanceName, m_processName);
            event.getMany (visSel, digiCollections);
        }
        else
        {
            event.getManyByType (digiCollections);
        }
        
        eventSetup.get<MuonGeometryRecord> ().get (m_pDD);
    }
    catch (cms::Exception& e)
    {
        if (this->m_onCmsException)
            this->m_onCmsException (&e);
    }
    catch (lat::Error &e) 
    {
        if (this->m_onError)
            this->m_onError (&e);
    }
    catch (std::exception &e) 
    {
        if (this->m_onException)
            this->m_onException (&e);
    }
    catch (...) 
    {
        if (this->m_onUnhandledException)
            this->m_onUnhandledException ();
    }

    if (! digiCollections.empty ())
    {
        std::vector<edm::Handle<CSCWireDigiCollection> >::iterator i;
        std::vector<edm::Handle<CSCWireDigiCollection> >::iterator iEnd;
        for (i = digiCollections.begin (), iEnd = digiCollections.end (); i != iEnd; ++i) 
        {
            const CSCWireDigiCollection& c = *(*i);
        
            // Loop over the DetUnits with digis
            CSCWireDigiCollection::DigiRangeIterator detUnitIt;
            CSCWireDigiCollection::DigiRangeIterator detUnitItEnd;
            for (detUnitIt = c.begin (), detUnitItEnd = c.end (); detUnitIt != detUnitItEnd; ++detUnitIt)
            {
                const CSCDetId &id = (*detUnitIt).first;
                const CSCWireDigiCollection::Range &range = (*detUnitIt).second;

                // Loop over the digis of this DetUnit
                for (CSCWireDigiCollection::const_iterator digiIt = range.first;
                     digiIt != range.second;
                     ++digiIt)
                {
                    std::pair< CSCDetId, CSCWireDigi > p;
                    p.first = id;
                    p.second = (*digiIt);
                    m_digis.push_back (p);
//                  std::cout << "Hit: Id/Digi: " << id <<", " << *digiIt << std::endl; //DEBUG
                } // for digis in layer
            } // for layers
        }
    }
    
    VisQueuedTwig::onBaseInvalidate ();
}

void
VisCSCWireDigiTwig::update (IgTextRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    // Prepare property description.
    std::ostringstream  text;
    text << "Total " << m_digis.size () << " CSC wire digis: " << m_text << "<br>";

    text << "<table width='100%' border=1>"
         << "<TR align = center>"
         << "<TH>Number</TH>"
         << "<TH>Position</TH>"
         << "<TH>Wire group</TH>"
         << "<TH>Beam-crossing tag</TH>"
         << "<TH>Time bin</TH>"
         << "</TR>";
    text << setiosflags (std::ios::showpoint | std::ios::fixed);
    text.setf (std::ios::right, std::ios::adjustfield);
    
    if ( ! m_digis.empty () ) 
    {
       
        int nDigis = 0;
        try 
        {
            for (std::vector<std::pair<CSCDetId, CSCWireDigi> >::const_iterator i = m_digis.begin (), iEnd = m_digis.end (); i != iEnd; ++i) 
            {
                Surface::GlobalPoint pos( NAN, NAN, NAN );
                if (m_pDD.isValid ()) pos  =  (m_pDD->idToDet ((*i).first))->surface ().position ();
                text << "<TR align = right>"
                     << "<TD>" << std::setw (3) << ++nDigis << "</TD>"
                     << "<TD>" << pos.x () << ", " << pos.y () << ", " << pos.z () << "</TD>"
                     << "<TD>" << std::setw (5) << std::setprecision (3) << (*i).second.getWireGroup () << "</TD>"
                     << "<TD>" << std::setw (5) << std::setprecision (3) << (*i).second.getBeamCrossingTag () << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).second.getTimeBin () << "</TD>"
                     << "</TR>";
            }
        }
        catch (cms::Exception& e)
        {
            if (this->m_onCmsException)
                this->m_onCmsException (&e);
        }
        catch (lat::Error &e) 
        {
            if (this->m_onError)
                this->m_onError (&e);
        }
        catch (std::exception &e) 
        {
            if (this->m_onException)
                this->m_onException (&e);
        }
        catch (...) 
        {
            if (this->m_onUnhandledException)
                this->m_onUnhandledException ();
        }
    }
    text << "</table>";
    
    // Send it over.
    IgQtLock ();
    rep->setText (text.str ());
}

void
VisCSCWireDigiTwig::update (Ig3DRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    //
    IgQtLock ();
    rep->clear ();
    m_dets.clear ();
    //
    if (! m_digis.empty () && m_pDD.isValid ())
    {
        VisCSCChamberDrawer d;
        const float width= 0.02;// wire group visualization 
        const float depth = 0.01;// wire group visualization 
        //    
        //
        SoDrawStyle *drawStyle = new SoDrawStyle;
        drawStyle->pointSize = 5.0; // big point; useful for the centers of the layers
        //
        SoSeparator *sep = new SoSeparator;// global
        sep->addChild (drawStyle);
        //
        try 
        { 
            SoMaterial *matl = new SoMaterial;
            matl->diffuseColor.setValue (0.5, 0.0, 1.0); // purple
            matl->emissiveColor.setValue(0.5, 0.0, 1.0);          
            matl->shininess = 1.;  
            //
            SoDrawStyle *drawStylel = new SoDrawStyle;
            drawStylel->pointSize = 1.;
            //
            SoSeparator *pulses = new SoSeparator;// wire groups
            pulses->addChild (matl);
            pulses->addChild (drawStylel);
            //
            for (std::vector<std::pair<CSCDetId, CSCWireDigi> >::const_iterator 
                     i = m_digis.begin (), iEnd = m_digis.end (); i != iEnd; ++i)
            {
                const CSCDetId& id = (*i).first;
                if (m_dets.find(id) == m_dets.end() )
                {
                    m_dets.insert(id);  
                    // This is a frame around the (active) layer
                    d.getOutline(sep, m_pDD, id);

                }
                //
                const CSCLayer *layer_p = m_pDD->layer (id);
                const CSCLayerGeometry *layerGeom = layer_p->geometry ();
                const int wireGroup_id = ((*i).second).getWireGroup ();
                const LocalPoint centerWireGroup = layerGeom->localCenterOfWireGroup (wireGroup_id);
                const Surface::GlobalPoint pos = (m_pDD->idToDet (id))->surface ().toGlobal (centerWireGroup);
                SbVec3f pt( pos.x(), pos.y(), pos.z() );
                pt *= 0.01F; // cm -> m
                //
                SoTranslation   *inc = new SoTranslation;
                inc->translation = pt;
                //
                SbVec3f axis (0., 0., 1.);      
                float angle = -atan (pt[0] / pt[1]) - M_PI / 2;
                SbRotation r (axis, angle);
                //
                SoTransform *xform = new SoTransform;
                xform ->rotation = r;
                //
                SoCube *hit = new SoCube;
                hit->width  = width;
                hit->height = layerGeom->lengthOfWireGroup (wireGroup_id) / 100.;
                hit->depth  = depth;
                //
                SoSeparator *pulse = new SoSeparator;// wire group
                pulse->addChild (inc);
                pulse->addChild (xform);
                pulse->addChild (hit);
                pulses->addChild (pulse);
                // Visualize wire groups 
                sep->addChild (pulses);// wire groups
            }
        }
        catch (cms::Exception& e)
        {
            if (this->m_onCmsException)
                this->m_onCmsException (&e);
        }
        catch (lat::Error &e) 
        {
            if (this->m_onError)
                this->m_onError (&e);
        }
        catch (std::exception &e) 
        {
            if (this->m_onException)
                this->m_onException (&e);
        }
        catch (...) 
        {
            if (this->m_onUnhandledException)
                this->m_onUnhandledException ();
        }
        rep->node ()->addChild (sep);
    }
}

void
VisCSCWireDigiTwig::update (IgLegoRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    IgQtLock ();
    rep->clear ();
    
//     SoSeparator *sep = new SoSeparator;    
//     SoSeparator *posSep = new SoSeparator;    
//     SoSeparator *negSep = new SoSeparator;

//     SoMaterial *mat = new SoMaterial;
//     mat->diffuseColor.setValue (0.0, 0.0, 1.0);
//     posSep->addChild (mat);

//     SoMaterial *negMat = new SoMaterial;
//     negMat->diffuseColor.setValue (1.0, 0.0, 0.0);
//     negSep->addChild (negMat);

//     sep->addChild (posSep);
//     sep->addChild (negSep);

//     try 
//     {
//      for (std::vector<std::pair<CSCDetId, CSCWireDigi> >::const_iterator i = m_digis.begin (); i != m_digis.end (); ++i) 
//      {
//          Surface::GlobalPoint pos =  (m_pDD->idToDet ((*i).first))->surface ().position ();
//          double eta = pos.eta ();
//          double phi = pos.phi ();
//          if (phi < 0.0) phi = 2 * M_PI + phi; // correction in absence of one convention

// //       int wire = (*i).second.getWireGroup ();
// //       int tbin = (*i).second.getBeamCrossingTag ();
            
//          IgSoTower *tower = new IgSoTower;
//          tower->position = SbVec2f (phi, eta); // eta, phi
//          tower->scaleFactor = 20.0;
//          posSep->addChild (tower); 
//      }
//     }
//      catch (cms::Exception& e)
//      {
//          if (this->m_onCmsException)
//              this->m_onCmsException (&e);
//      }
//      catch (lat::Error &e) 
//      {
//          if (this->m_onError)
//              this->m_onError (&e);
//      }
//      catch (std::exception &e) 
//      {
//          if (this->m_onException)
//              this->m_onException (&e);
//      }
//      catch (...) 
//      {
//          if (this->m_onUnhandledException)
//              this->m_onUnhandledException ();
//      }
    
//     rep->node ()->addChild (sep);
}

void
VisCSCWireDigiTwig::update (IgRPhiRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    //
    IgQtLock ();
    rep->clear ();
    //
    if (! m_digis.empty () && m_pDD.isValid ())
    {
        const float width= 0.02;// wire group visualization 
        const float depth = 0.01;// wire group visualization 
        //    
        SbColor grey35;
        float transparency = 0.0f;
        grey35.setPackedValue(static_cast<uint32_t>(0x00595959), transparency);
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (grey35); // grey35
        mat->emissiveColor.setValue(grey35);
        //
        SoDrawStyle *drawStyle = new SoDrawStyle;
        drawStyle->pointSize = 5.0; // big point; useful for the centers of the layers
        //
        SoSeparator *sep = new SoSeparator;// global
        sep->addChild (mat);
        sep->addChild (drawStyle);
        //
        try 
        { 
            SoMaterial *matl = new SoMaterial;
            matl->diffuseColor.setValue (0.5, 0.0, 1.0); // purple
            matl->emissiveColor.setValue(0.5, 0.0, 1.0);          
            matl->shininess = 1.;  
            //
            SoDrawStyle *drawStylel = new SoDrawStyle;
            drawStylel->pointSize = 1.;
            //
            SoSeparator *pulses = new SoSeparator;// wire groups
            pulses->addChild (matl);
            pulses->addChild (drawStylel);
            //
            SoVertexProperty *vertices = new SoVertexProperty;
            int nVrtx = 0;
            for (std::vector<std::pair<CSCDetId, CSCWireDigi> >::const_iterator 
                     i = m_digis.begin (), iEnd = m_digis.end (); i != iEnd; ++i)
            {
                const CSCDetId& id = (*i).first;
                // This is the center of the (active) layer
                Surface::GlobalPoint pos =  (m_pDD->idToDet (id))->surface ().position (); 
                float x = pos.x () / 100.0;  // cm -> m
                float y = pos.y () / 100.0;  // cm -> m
                float z = pos.z () / 100.0;  // cm -> m
                vertices->vertex.set1Value (nVrtx++, SbVec3f (x, y, z));
                //
                const CSCLayer* layer_p = m_pDD->layer(id);
                const CSCLayerGeometry* layerGeom=0;
                if (layer_p)
                {
                    layerGeom = layer_p->geometry ();
                }
                else 
                {
                    // to be filled...
                }
                // This is a frame around the (active) layer
                bool framePlot = true;
                if (framePlot)
                {
                    const std::vector<float> trapezBounds = layerGeom->parameters ();
                    std::vector<Surface::GlobalPoint> CrossPoint;
                    CrossPoint.push_back ((m_pDD->idToDet (id))->surface ().toGlobal
                                          (LocalPoint (-trapezBounds [0], -trapezBounds [3], 0.)));
                    CrossPoint.push_back ((m_pDD->idToDet (id))->surface ().toGlobal
                                          (LocalPoint (trapezBounds [0], -trapezBounds [3], 0.)));
                    CrossPoint.push_back ((m_pDD->idToDet (id))->surface ().toGlobal
                                          (LocalPoint(trapezBounds [1], trapezBounds [3], 0.)));
                    CrossPoint.push_back ((m_pDD->idToDet (id))->surface ().toGlobal
                                          (LocalPoint (-trapezBounds [1], trapezBounds [3], 0.)));
                    //
                    SoVertexProperty *vert = new SoVertexProperty;
                    int nVrx = 0;
                    for (std::vector<Surface::GlobalPoint>::iterator _crossP = CrossPoint.begin (); _crossP != CrossPoint.end (); ++_crossP)
                    {
                        x = _crossP->x () / 100.0;  // cm -> m
                        y = _crossP->y () / 100.0;  // cm -> m
                        z = _crossP->z () / 100.0;  // cm -> m
                        vert->vertex.set1Value (nVrx++, SbVec3f (x, y, z));
                    }
                    std::vector<Surface::GlobalPoint>::reference closeCircle 
                        =  CrossPoint.front ();
                    x = closeCircle.x () / 100.0;  // cm -> m
                    y = closeCircle.y () / 100.0;  // cm -> m
                    z = closeCircle.z () / 100.0;  // cm -> m
                    vert->vertex.set1Value (nVrx++, SbVec3f (x, y, 0.f));
                    // 
                    vert->vertex.setNum (nVrx);
                    SoLineSet *trapez = new SoLineSet;
                    trapez->vertexProperty.setValue (vert);
                    trapez->numVertices [CrossPoint.size () + 1];
                    //
                    SoMaterial *mats = new SoMaterial;
                    mats->diffuseColor.setValue (1.0, 1.0, 1.0); // white
                    mats->emissiveColor.setValue (1.0, 1.0, 1.0);
                    //
                    SoSeparator *shape = new SoSeparator;// layer
                    shape->addChild (mats);
                    shape->addChild (trapez);
                    // Visualize the layer bounds
                    sep->addChild (shape);
                }
                //
                int wireGroup_id = ((*i).second).getWireGroup ();
                LocalPoint centerWireGroup = layerGeom->localCenterOfWireGroup (wireGroup_id);
                pos = (m_pDD->idToDet (id))->surface ().toGlobal (centerWireGroup);
                x = pos.x () / 100.0;  // cm -> m
                y = pos.y () / 100.0;  // cm -> m
                z = pos.z () / 100.0;  // cm -> m
                //
                SoTranslation   *inc = new SoTranslation;
                inc->translation = SbVec3f (x, y, z);
                //
                SbVec3f axis (0., 0., 1.);      
                float angle = -atan (x / y) - M_PI / 2;
                SbRotation r (axis, angle);
                //
                SoTransform *xform = new SoTransform;
                xform ->rotation = r;
                //
                SoCube *hit = new SoCube;
                hit->width  = width;
                hit->height = layerGeom->lengthOfWireGroup (wireGroup_id) / 100.;
                hit->depth  = depth;
                //
                SoSeparator *pulse = new SoSeparator;// wire group
                pulse->addChild (inc);
                pulse->addChild (xform);
                pulse->addChild (hit);
                pulses->addChild (pulse);
                // Visualize wire groups 
                sep->addChild (pulses);// wire groups
            }
            //
            vertices->vertex.setNum (m_digis.size ());      
            SoPointSet *points = new SoPointSet;
            points->vertexProperty.setValue (vertices);
            points->numPoints.setValue (m_digis.size ());
            // Visualize the center of the layer
            sep->addChild (points);
        }
        catch (cms::Exception& e)
        {
            if (this->m_onCmsException)
                this->m_onCmsException (&e);
        }
        catch (lat::Error &e) 
        {
            if (this->m_onError)
                this->m_onError (&e);
        }
        catch (std::exception &e) 
        {
            if (this->m_onException)
                this->m_onException (&e);
        }
        catch (...) 
        {
            if (this->m_onUnhandledException)
                this->m_onUnhandledException ();
        }
        rep->node ()->addChild (sep);
    }
}

void
VisCSCWireDigiTwig::update (IgSpareRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    //
    IgQtLock ();
    rep->clear ();
    //
    if (! m_digis.empty () && m_pDD.isValid ()) 
    {
        SoSeparator *sep = new SoSeparator;    
        //
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (0.0, 0.0, 0.0);
        sep->addChild (mat);
        try 
        {
            bool iterations = false;
            float OurPhi = 9999999.;
            float OurPhiSlice = 9999999.;
            float OurZ = 99999999.;
            for (std::vector<std::pair<CSCDetId, CSCWireDigi> >::const_iterator i = 
                     m_digis.begin (); i != m_digis.end (); ++i) 
            {
                const CSCDetId& id = (*i).first;
                //int EndCap  = id.endcap();
                //int Station = id.station();
                //int Ring    = id.ring();
                //int Chamber = id.chamber();
                const CSCLayer* layer_p = m_pDD->layer(id);
                const CSCLayerGeometry* layerGeom=0;
              
                if(layer_p){
                    layerGeom = layer_p->geometry();
                }
                else{// to be filled...
                }
                //
                Surface::GlobalPoint pos =  (m_pDD->idToDet (id))->surface ().position ();
                float layer_x = pos.x () / 100.0;  // cm -> m
                float layer_y = pos.y () / 100.0;  // cm -> m
                float layer_z = pos.z () / 100.0;  // cm -> m
                float layer_phi = pos.phi ();
                float layer_r = sqrt(layer_x*layer_x + layer_y*layer_y);
                if(!iterations){
                    // Make captions, some init
                    iterations=VisCaptions(layer_z, layer_phi, OurPhi, OurZ, sep);
                }
                const std::vector<float> trapezBounds = layerGeom->parameters();
                OurPhiSlice = 2. * atan(trapezBounds[1]/(layer_r + trapezBounds[3])) * 0.8;
                // Only adjacent chambers 
                if (fabs(layer_phi - OurPhi) < 2. * OurPhiSlice){
                    int wireGroup_id = ((*i).second).getWireGroup();
                    LocalPoint centerWireGroup=layerGeom->localCenterOfWireGroup(wireGroup_id);
                    float wireLength = layerGeom->lengthOfWireGroup(wireGroup_id);
                    // Make frames
                    SoLineSet *trapez = new SoLineSet;
                    SoLineSet *rect_s = new SoLineSet;
                    SoLineSet *rect_w = new SoLineSet;
                    // bool FramesOK = VisFrames(trapezBounds,trapez,rect_w,rect_s);
                    // Translations
                    SoTranslation *trans = new SoTranslation; // trapez frames
                    SoTranslation *trans_caption = new SoTranslation; // captions
                    SoTranslation *trans_s = new SoTranslation; // strips
                    SoTranslation *trans_w = new SoTranslation; // wire groups
                    /* 
                    bool TransOK = MakeTrans(trapezBounds,
                                             id.station(), layer_r, layer_phi, layer_z, 
                                             OurPhi, OurZ,
                                             trans, trans_caption, trans_s, trans_w);
                    */                       
                    // caption
                    SoText3 *TextCaption = new SoText3;
                    char ch[10];
                    sprintf( ch,"%1.1f", layer_phi);
                    TextCaption->string = ch;
                    // font
                    SoFont *font_caption = new SoFont;
                    font_caption->name.setValue(VisFont);
                    font_caption->size.setValue(VisFontSize);
                    // Wire group as a line
                    SoVertexProperty *vert = new SoVertexProperty;
                    int nVrx = 0;
                    float x = (centerWireGroup.x()-wireLength/2.)/100.;
                    float y = centerWireGroup.y()/100.;
                    float z = 0.; 
                    vert->vertex.set1Value (nVrx++,SbVec3f (x, y, z) );
                    //
                    x = (centerWireGroup.x()+wireLength/2.)/100.;
                    y = centerWireGroup.y()/100.;
                    z = 0.;
                    vert->vertex.set1Value (nVrx++,SbVec3f (x, y, z) );
                    //
                    vert->vertex.setNum (nVrx);
                    SoLineSet *WireGroupVis = new SoLineSet;
                    WireGroupVis->vertexProperty.setValue (vert);
                    WireGroupVis->numVertices[2];
                    // Wire group projection - a point
                    SoVertexProperty *wireDots = new SoVertexProperty;
                    nVrx = 0;
                    x =  centerWireGroup.y()/ 100.0;  // cm -> m
                    y =  (trapezBounds[2]*1.1 -  id.layer()*2.*trapezBounds[2]/6.)/4.;
                    z = 0. / 100.0;  // cm -> m
                    //  " ,x:"<<x<<" ,y:"<<y<<std::endl;
                    if(3==id.station()) y=-y;// ...
                    //std::cout<<"new y:"<<y<<std::endl;
                    wireDots->vertex.set1Value (nVrx++, SbVec3f (x, y, z));
                    wireDots->vertex.setNum (nVrx);
                    int Npoints=nVrx;
                    // Colors
                    SoMaterial *matw = new SoMaterial;
                    SoMaterial *matf = new SoMaterial;
                    if(OurPhi == layer_phi){
                        matw->diffuseColor.setValue (0.3, 0.0, 0.3); // purple
                        matf->diffuseColor.setValue (0.0, 0.0, 0.0); // black
                    }
                    else{
                        matw->diffuseColor.setValue (0.6, 0.0, 6.0); // light purple
                        matf->diffuseColor.setValue(0.6, 0.6, 0.6); // grey
                    }
                    //  Build the "picture"
                    SoSeparator *shape = new SoSeparator;// frame
                    shape->addChild (matf); 
                    shape->addChild (trapez);
                    //
                    SoSeparator *caption = new SoSeparator;// caption
                    caption->addChild (matf); 
                    caption->addChild (trans_caption);
                    caption->addChild (font_caption);
                    caption->addChild (TextCaption);
                    //
                    SoSeparator *proj_s = new SoSeparator;// projection strip frame
                    proj_s->addChild (matf); 
                    proj_s->addChild (trans_s); 
                    proj_s->addChild (rect_s);
                    //
                    SoSeparator *proj_w = new SoSeparator;// projection wire frame
                    proj_w->addChild (matf); 
                    proj_w->addChild (trans_w); 
                    proj_w->addChild (rect_w);
                    //
                    SoDrawStyle *drawStylel = new SoDrawStyle;
                    drawStylel->pointSize = 2.0; // bigger 
                    //
                    SoSeparator *dotsWG = new SoSeparator; // projection wire groups
                    SoPointSet *points = new SoPointSet;
                    points->vertexProperty.setValue (wireDots);
                    points->numPoints.setValue (Npoints);
                    dotsWG->addChild (matw); 
                    dotsWG->addChild (drawStylel);
                    dotsWG->addChild (trans_w);
                    dotsWG->addChild (points);
                    //
                    SoSeparator *wiregroup = new SoSeparator;// wiregroup
                    wiregroup->addChild (matw);
                    wiregroup->addChild (WireGroupVis);
                    //
                    SoSeparator *layer = new SoSeparator;// layer
                    layer->addChild (trans);
                    layer->addChild (shape);
                    layer->addChild (wiregroup);
                    // Visualization of all
                    sep->addChild (layer);
                    sep->addChild (caption);
                    sep->addChild (proj_s);
                    sep->addChild (proj_w);
                    sep->addChild (dotsWG);
                }
            }
        }
        catch (cms::Exception& e)
        {
            if (this->m_onCmsException)
                this->m_onCmsException (&e);
        }
        catch (lat::Error &e) 
        {
            if (this->m_onError)
                this->m_onError (&e);
        }
        catch (std::exception &e) 
        {
            if (this->m_onException)
                this->m_onException (&e);
        }
        catch (...) 
        {
            if (this->m_onUnhandledException)
                this->m_onUnhandledException ();
        }
        //
        // FIXME: Move to another rep type. 
        rep->node ()->addChild (sep);
    }
}

void
VisCSCWireDigiTwig::update (IgRZRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    IgQtLock ();
    rep->clear ();

    if (! m_digis.empty () && m_pDD.isValid ()) 
    {
        SoSeparator *sep = new SoSeparator;    
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (0.3, 0.8, 0.0);
        sep->addChild (mat);

        SoDrawStyle *drawStyle = new SoDrawStyle;
        drawStyle->pointSize = 3.0;
        sep->addChild (drawStyle);
    
        try 
        {
            SoVertexProperty *vertices = new SoVertexProperty;
            int nVrx = 0;

            for (std::vector<std::pair<CSCDetId, CSCWireDigi> >::const_iterator i = m_digis.begin (), iEnd = m_digis.end (); i != iEnd; ++i) 
            {
                Surface::GlobalPoint pos =  (m_pDD->idToDet ((*i).first))->surface ().position ();
//              int wire = (*i).second.getWireGroup ();
//              int tbin = (*i).second.getBeamCrossingTag ();
            
                float x = pos.x () / 100.0;  // cm -> m
                float y = pos.y () / 100.0;  // cm -> m
                float z = pos.z () / 100.0;  // cm -> m
                
                vertices->vertex.set1Value (nVrx++, SbVec3f (x, y, z));
            }
            vertices->vertex.setNum (m_digis.size ());

            SoPointSet *points = new SoPointSet;
            points->vertexProperty.setValue (vertices);
            points->numPoints.setValue (m_digis.size ());
            sep->addChild (points);
        }    
        catch (cms::Exception& e)
        {
            if (this->m_onCmsException)
                this->m_onCmsException (&e);
        }
        catch (lat::Error &e) 
        {
            if (this->m_onError)
                this->m_onError (&e);
        }
        catch (std::exception &e) 
        {
            if (this->m_onException)
                this->m_onException (&e);
        }
        catch (...) 
        {
            if (this->m_onUnhandledException)
                this->m_onUnhandledException ();
        }
        rep->node ()->addChild (sep);
    }
}

---