---

VisHcalDetTwig.cc

Go to the documentation of this file.

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

#include "VisReco/VisHcal/interface/VisHcalDetTwig.h"
#include "Iguana/Framework/interface/IgRepSet.h"
#include "Iguana/Models/interface/IgTextRep.h"
#include "Iguana/GLModels/interface/IgLegoRep.h"
#include "Iguana/GLModels/interface/IgRPhiRep.h"
#include "Iguana/GLModels/interface/IgRZRep.h"
#include "Iguana/GLBrowsers/interface/Ig3DShapeManager.h"
#include "Iguana/Studio/interface/IgQtLock.h"
#include <Inventor/nodes/SoIndexedLineSet.h>
#include <Inventor/nodes/SoSeparator.h>
#include <Inventor/nodes/SoText2.h>
#include <Inventor/nodes/SoTranslation.h>
#include <Inventor/SoInput.h>
#include <qfileinfo.h>
#include <qmessagebox.h>
#include <sstream>
#include <iostream>
#include <iomanip>

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

VisHcalDetTwig::VisHcalDetTwig (IgState *state, IgTwig *parent)
    : VisQueuedTwig (state, parent, "HCAL Detector"),
      m_scale (state, lat::CreateCallback (this, &VisHcalDetTwig::scaleChanged))
{
    Ig3DShapeManager *shapeManager = Ig3DShapeManager::get (state);
    ASSERT (shapeManager);
    
    m_det = shapeManager->lookup ("hcalbarrel.iv");
    m_legoDet = buildLegoDet ();
    m_legoDet->ref ();
}

void
VisHcalDetTwig::scaleChanged (void)
{     
    IgRepSet::invalidate (this, SELF_MASK); 
}
    
SoNode*         
VisHcalDetTwig::buildLegoDet (void) 
{
    SoSeparator *gridSep = new SoSeparator;

    // HF
    float etaHFBins [43] = { 0, 0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609,
                             0.696, 0.783, 0.87, 0.957, 1.044, 1.131, 1.218, 1.305,
                             1.392, 1.479, 1.566, 1.653, 1.74, 1.83, 1.93, 2.043,
                             2.172, 2.322, 2.5, 2.65, 3,
                             2.853, 2.964, 3.139, 3.314, // HF
                             3.486, 3.664, 3.839, 4.013, 
                             4.191, 4.363, 4.538, 4.716, 
                             4.889, 5.191};
    
    SoVertexProperty *vtx = new SoVertexProperty;
    SoIndexedLineSet *grid = new SoIndexedLineSet;
    SoMFInt32 coords;
    int vertex = 6;
    int coord = 0;
    int tick = 0;
    
    float x = float (2 * M_PI);
    float z = 5.191;
    
    vtx->vertex.set1Value (0, SbVec3f (0, 0, -z));
    vtx->vertex.set1Value (1, SbVec3f (0, 0,  z));
    vtx->vertex.set1Value (2, SbVec3f (x, 0,  z));
    vtx->vertex.set1Value (3, SbVec3f (x, 0, -z));
    vtx->vertex.set1Value (4, SbVec3f (0, 0, -z));
    vtx->vertex.set1Value (5, SbVec3f (0, 1, -z));
    
    for (coord = 0; coord < vertex; coord++)
        coords.set1Value (coord, coord);
    coords.set1Value (coord++, SO_END_LINE_INDEX);
    
    for (tick = 0; tick < 43; tick++)
    {
        vtx->vertex.set1Value (vertex, SbVec3f (0, 0, etaHFBins [tick]));       
        coords.set1Value (coord++, vertex++);
        vtx->vertex.set1Value (vertex, SbVec3f (x, 0, etaHFBins [tick]));
        coords.set1Value (coord++, vertex++);
        coords.set1Value (coord++, SO_END_LINE_INDEX);
        
        vtx->vertex.set1Value (vertex, SbVec3f (0, 0, -etaHFBins [tick]));
        coords.set1Value (coord++, vertex++);
        vtx->vertex.set1Value (vertex, SbVec3f (x, 0, -etaHFBins [tick]));
        coords.set1Value (coord++, vertex++);
        coords.set1Value (coord++, SO_END_LINE_INDEX);
    }
    
    int ticks = 36;
    x = 0;
    
    for (tick = 0; tick < ticks; tick++)
    {
        x = tick * M_PI / 18.0;
        vtx->vertex.set1Value (vertex, SbVec3f (x, 0, -z));
        coords.set1Value (coord++, vertex++);
        vtx->vertex.set1Value (vertex, SbVec3f (x, 0, z));
        coords.set1Value (coord++, vertex++);
        coords.set1Value (coord++, SO_END_LINE_INDEX);
        
        for (int ttow = 1; ttow < 2; ttow++)
        {
            vtx->vertex.set1Value (vertex, SbVec3f (x + ttow * 0.087, 0, -1.74));
            coords.set1Value (coord++, vertex++);
            vtx->vertex.set1Value (vertex, SbVec3f (x + ttow * 0.087, 0, 1.74));
            coords.set1Value (coord++, vertex++);
            coords.set1Value (coord++, SO_END_LINE_INDEX);
        }
    }
    
    grid->vertexProperty = vtx;
    grid->coordIndex = coords;
    
    gridSep->addChild (grid);
    
    SoSeparator         *phiLabels = new SoSeparator;
    SoText2             *text;
    SoText2             *phiLabelText [8];
    SoTranslation       *phiStart = new SoTranslation;
    SoTranslation       *phiTranslation = new SoTranslation;
    
    z = 5.191;
    
    phiStart->translation = SbVec3f (M_PI / 4, 0, -z - 0.6);
    phiTranslation->translation = SbVec3f (M_PI / 4, 0, 0);
    phiLabels->addChild (phiStart);
    
    for (tick = 0; tick < 8; tick ++) 
    {
        text = phiLabelText [tick] = new SoText2;
        phiLabels->addChild (text);
        phiLabels->addChild (phiTranslation);
    }
    
    phiLabelText [0]->string = "PI/4";
    phiLabelText [1]->string = "PI/2";
    phiLabelText [2]->string = "3PI/4";
    phiLabelText [3]->string = "PI";
    phiLabelText [4]->string = "5PI/4";
    phiLabelText [5]->string = "3PI/2";
    phiLabelText [6]->string = "7PI/4";
    phiLabelText [7]->string = "2PI";
    
    gridSep->addChild (phiLabels);
    
    z = -5.191;
    
    // scale
    SoSeparator         *labelScale = new SoSeparator;
    SoText2             *labelScaleText = new SoText2;
    SoTranslation       *labelScaleOffset = new SoTranslation;
    
    char scaleChars [12] = "1.0 GeV";
    sprintf (scaleChars, "%.2G GeV", m_scale.value ());
    labelScaleText->string = scaleChars;
    labelScaleOffset->translation
        = SbVec3f (-0.6, 1, z - 0.6);
    labelScale->addChild (labelScaleOffset);
    labelScale->addChild (labelScaleText);
    gridSep->addChild (labelScale);
    
    SoSeparator         *etaLabels = new SoSeparator;
    SoTranslation       *etaStart = new SoTranslation;
    SoTranslation       *etaTranslation = new SoTranslation;
    
    etaStart->translation = SbVec3f (-0.6, 0, -5.0);
    etaTranslation->translation = SbVec3f (0, 0, 1);
    etaLabels->addChild (etaStart);
    
    for (tick = 0; tick < 11; tick++, z++)
    {
        text = new SoText2;
        char textLabel [5];
        sprintf (textLabel, "%.0f", z);
        text->string = textLabel;
        etaLabels->addChild (text);
        etaLabels->addChild (etaTranslation);
    }
    gridSep->addChild (etaLabels);

    return gridSep;
}

void
VisHcalDetTwig::update (IgTextRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    
    // Prepare property description.
    std::ostringstream  text;
    
    text << "<table width='100%'><tr><td>HCAL detector</td><td>";
    
    text << "</td></tr></table>";
    
    // Send it over.
    IgQtLock ();

    rep->setText (text.str ());
}

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

    IgQtLock ();

    rep->clear ();    
    m_legoDet = buildLegoDet ();
    rep->node ()->addChild (m_legoDet);
}

void
VisHcalDetTwig::update (IgRPhiRep *rep)
{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    IgQtLock ();

    rep->clear ();
    rep->node ()->addChild (m_det);
}

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

    IgQtLock ();

    rep->clear ();
    rep->node ()->addChild (m_det);
}

---