---

VisPSimHitTwig Class Reference

#include <VisSimulation/VisSimHits/interface/VisPSimHitTwig.h>

Inheritance diagram for VisPSimHitTwig:

List of all members.

Public Member Functions
virtual void	onNewEvent (const edm::Event &event, const edm::EventSetup &eventSetup)
virtual void	twigChanged (void)
virtual void	update (IgRZRep *rep)
virtual void	update (IgRPhiRep *rep)
virtual void	update (IgLegoRep *rep)
virtual void	update (Ig3DRep *rep)
virtual void	update (IgTextRep *rep)
	VisPSimHitTwig (IgState *state, IgTwig *parent, const std::string &name="", const std::string &friendlyName="", const std::string &moduleLabel="", const std::string &instanceName="", const std::string &processName="")
Private Attributes
VisSimHitConnect	m_connect
std::string	m_friendlyName
std::vector< PSimHit >	m_hits
std::string	m_instanceName
std::string	m_moduleLabel
std::string	m_processName
std::string	m_text

---

Detailed Description

Definition at line 20 of file VisPSimHitTwig.h.

---

Constructor & Destructor Documentation

VisPSimHitTwig::VisPSimHitTwig	(	IgState *	state,
		IgTwig *	parent,
		const std::string &	name = "",
		const std::string &	friendlyName = "",
		const std::string &	moduleLabel = "",
		const std::string &	instanceName = "",
		const std::string &	processName = ""
	)

Definition at line 64 of file VisPSimHitTwig.cc.

References createThisTwig(), edm::TypeID::friendlyClassName(), DBSPlugin::get(), and VisTwigFactroyService::registerTwig().

    : VisQueuedTwig (state, parent, name),
      m_text (name),
      m_friendlyName (friendlyName),
      m_moduleLabel (moduleLabel),
      m_instanceName (instanceName),
      m_processName (processName),
      m_connect (state, lat::CreateCallback (this, &VisPSimHitTwig::twigChanged))
{
    VisTwigFactroyService *tfService = VisTwigFactroyService::get (state);
    if (! tfService)
    {
        tfService = new VisTwigFactroyService (state);
    }   
    edm::TypeID caloJetsID (typeid (edm::PSimHitContainer));
    tfService->registerTwig (caloJetsID.friendlyClassName (), &createThisTwig);
}

---

Member Function Documentation

void VisPSimHitTwig::onNewEvent	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)			[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 91 of file VisPSimHitTwig.cc.

References cms::Exception::append(), arg, edm::Event::bunchCrossing(), c, e, exception, i, edm::Event::id(), IgRepSet::invalidate(), edm::Event::luminosityBlock(), m_friendlyName, m_hits, m_instanceName, m_moduleLabel, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_processName, m_text, IgSimpleTwig::name(), VisQueuedTwig::onNewEvent(), edm::Event::orbitNumber(), IgTwig::SELF_MASK, and IgTwig::STRUCTURE_MASK.

{
    // Get debugging dump.
    VisQueuedTwig::onNewEvent (event, eventSetup);

    m_hits.clear ();
    m_text = (QString ("Run %1, Event %2, LS %3, Orbit %4, BX %5")
              .arg (event.id ().run ())
              .arg (event.id ().event ())
              .arg (event.luminosityBlock ())
              .arg (event.orbitNumber ())
              .arg (event.bunchCrossing ())
              .latin1 ());

    std::vector<edm::Handle<edm::PSimHitContainer> > hitCollections;
    
    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, hitCollections);
        }
    }
    catch (cms::Exception& e)
    {
        e.append (" from VisPSimHitTwig::onNewEvent ");
        e.append (this->name ());    

        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 (! hitCollections.empty ())
    {
        for (std::vector<edm::Handle<edm::PSimHitContainer> >::iterator i = hitCollections.begin (), iEnd = hitCollections.end (); i != iEnd; ++i) 
        {
            const edm::PSimHitContainer& c = *(*i);
            QString sizeStr = (QString ("%1").arg (c.size ()));
            QString nameStr = QString (this->name ());
            int ib = nameStr.find ("[");
            int ie = nameStr.find ("]");
            nameStr.replace (ib + 1, ie - 1, sizeStr);
            
            this->name (nameStr);

            for (std::vector<PSimHit>::const_iterator isim = c.begin (), isimEnd = c.end ();
                 isim != isimEnd; ++isim)
            {
                m_hits.push_back ((*isim));
            }    
        }
    }

    IgRepSet::invalidate (this, IgTwig::SELF_MASK | IgTwig::STRUCTURE_MASK);
}

void VisPSimHitTwig::twigChanged

(

void

)

[virtual]

Definition at line 87 of file VisPSimHitTwig.cc.

References IgRepSet::invalidate(), and IgTwig::SELF_MASK.

{ IgRepSet::invalidate (this, SELF_MASK); }

void VisPSimHitTwig::update

(

IgRZRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 485 of file VisPSimHitTwig.cc.

References cms::Exception::append(), ASSERT, e, Ig3DBaseModel::encode(), end, exception, DBSPlugin::get(), i, label, m_connect, m_hits, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, DetId::Muon, Ig3DBaseRep::node(), VisQueuedTwig::state(), VisEventSetupService::toGlobal(), track, DetId::Tracker, VisQueuedTwig::update(), VisActiveConfigurable< T >::value(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    if (! m_hits.empty ())
    {
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (0.0, 0.0, 1.0);
        mat->emissiveColor.setValue(0.0, 0.0, 1.0);

        SoDrawStyle *drawStyle = new SoDrawStyle;
        drawStyle->pointSize = 6.0;

        SoSeparator *sep = new SoSeparator; // global
        sep->addChild (mat);
        sep->addChild (drawStyle);

        VisEventSetupService *esService = VisEventSetupService::get (state ());
        ASSERT (esService);

        try 
        {
            typedef std::map<int, std::vector <SbVec3f> > VisParticleMap;
            VisParticleMap pmap;

            for (std::vector<PSimHit>::const_iterator i = m_hits.begin (); i != m_hits.end (); ++i) 
            {
                int ptype = (*i).trackId ();
                
                VisParticleMap::iterator mit = pmap.find (ptype);
                
                if (mit == pmap.end ())
                {
                    std::vector<SbVec3f> vec;
                    bool done = pmap.insert (VisParticleMap::value_type (ptype, vec)).second;
                    if (done)
                    {
                        mit = pmap.find (ptype);
                    }
                }

                if (mit == pmap.end ())
                {
                    throw std::runtime_error ("Cannot find particle type in the map!");
                }               
                DetId id = (*i).detUnitId ();
                if (! id.null () && (id.det () == DetId::Tracker || id.det () == DetId::Muon))
                {
                    GlobalPoint pos = esService->toGlobal ((*i).localPosition (), id);
                    (mit->second).push_back (SbVec3f (pos.x () / 100.0, pos.y () / 100.0, pos.z () / 100.0));
                }
            }       
            if (! pmap.empty ())
            {
                for (VisParticleMap::const_iterator mit = pmap.begin (); mit != pmap.end (); ++mit)
                {
                    SoSeparator *trackSep = new SoSeparator; 
                    
                    QString label = QString ("ID = %1")
                                    .arg (mit->first);             

                    trackSep->setName (Ig3DBaseModel::encode (label.prepend ("Track with ")));

                    SoVertexProperty *vertices = new SoVertexProperty;
                    int nVrtx = 0;

                    for (std::vector<SbVec3f>::const_iterator vit = (mit->second).begin (); vit != (mit->second).end (); ++vit) 
                    {   
                        vertices->vertex.set1Value (nVrtx++, (*vit));
                    }
                    
                    vertices->vertex.setNum (nVrtx);
                    SoPointSet *points = new SoPointSet;
                    points->vertexProperty.setValue (vertices);
                    points->numPoints.setValue (nVrtx);
                    trackSep->addChild (points);

                    if (m_connect.value ()) 
                    {
                        SoLineSet *track = new SoLineSet;
                        track->vertexProperty.setValue (vertices);
                        track->numVertices.setValue (nVrtx);
                        trackSep->addChild (track);
                    }
                    
                    sep->addChild (trackSep);               
                }
            }
        }
        catch (cms::Exception& e)
        {
            e.append (" from VisPSimHitTwig::update(IgRZRep*) ");
            e.append (this->name ());    

            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 VisPSimHitTwig::update

(

IgRPhiRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 367 of file VisPSimHitTwig.cc.

References cms::Exception::append(), ASSERT, e, Ig3DBaseModel::encode(), end, exception, DBSPlugin::get(), i, label, m_connect, m_hits, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, DetId::Muon, Ig3DBaseRep::node(), VisQueuedTwig::state(), VisEventSetupService::toGlobal(), track, DetId::Tracker, VisQueuedTwig::update(), VisActiveConfigurable< T >::value(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    if (! m_hits.empty ())
    {
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (0.0, 0.0, 1.0);
        mat->emissiveColor.setValue(0.0, 0.0, 1.0);

        SoDrawStyle *drawStyle = new SoDrawStyle;
        drawStyle->pointSize = 6.0;

        SoSeparator *sep = new SoSeparator; // global
        sep->addChild (mat);
        sep->addChild (drawStyle);

        VisEventSetupService *esService = VisEventSetupService::get (state ());
        ASSERT (esService);

        try 
        {
            typedef std::map<int, std::vector <SbVec3f> > VisParticleMap;
            VisParticleMap pmap;

            for (std::vector<PSimHit>::const_iterator i = m_hits.begin (); i != m_hits.end (); ++i) 
            {
                int ptype = (*i).trackId ();
                
                VisParticleMap::iterator mit = pmap.find (ptype);
                
                if (mit == pmap.end ())
                {
                    std::vector<SbVec3f> vec;
                    bool done = pmap.insert (VisParticleMap::value_type (ptype, vec)).second;
                    if (done)
                    {
                        mit = pmap.find (ptype);
                    }
                }

                if (mit == pmap.end ())
                {
                    throw std::runtime_error ("Cannot find particle type in the map!");
                }               
                DetId id = (*i).detUnitId ();
                if (! id.null () && (id.det () == DetId::Tracker || id.det () == DetId::Muon))
                {
                    GlobalPoint pos = esService->toGlobal ((*i).localPosition (), id);
                    (mit->second).push_back (SbVec3f (pos.x () / 100.0, pos.y () / 100.0, pos.z () / 100.0));
                }
            }       
            if (! pmap.empty ())
            {
                for (VisParticleMap::const_iterator mit = pmap.begin (); mit != pmap.end (); ++mit)
                {
                    SoSeparator *trackSep = new SoSeparator; 
                    
                    QString label = QString ("ID = %1")
                                    .arg (mit->first);             

                    trackSep->setName (Ig3DBaseModel::encode (label.prepend ("Track with ")));

                    SoVertexProperty *vertices = new SoVertexProperty;
                    int nVrtx = 0;

                    for (std::vector<SbVec3f>::const_iterator vit = (mit->second).begin (); vit != (mit->second).end (); ++vit) 
                    {   
                        vertices->vertex.set1Value (nVrtx++, (*vit));
                    }
                    
                    vertices->vertex.setNum (nVrtx);
                    SoPointSet *points = new SoPointSet;
                    points->vertexProperty.setValue (vertices);
                    points->numPoints.setValue (nVrtx);
                    trackSep->addChild (points);

                    if (m_connect.value ()) 
                    {
                        SoLineSet *track = new SoLineSet;
                        track->vertexProperty.setValue (vertices);
                        track->numVertices.setValue (nVrtx);
                        trackSep->addChild (track);
                    }
                    
                    sep->addChild (trackSep);               
                }
            }
        }
        catch (cms::Exception& e)
        {
            e.append (" from VisPSimHitTwig::update(IgRPhiRep*) ");
            e.append (this->name ());    

            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 VisPSimHitTwig::update

(

IgLegoRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 360 of file VisPSimHitTwig.cc.

References VisQueuedTwig::update().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
}

void VisPSimHitTwig::update

(

Ig3DRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 240 of file VisPSimHitTwig.cc.

References cms::Exception::append(), ASSERT, Ig3DBaseRep::clear(), e, Ig3DBaseModel::encode(), end, exception, DBSPlugin::get(), i, label, m_connect, m_hits, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, DetId::Muon, Ig3DBaseRep::node(), VisQueuedTwig::state(), VisEventSetupService::toGlobal(), track, DetId::Tracker, VisQueuedTwig::update(), VisActiveConfigurable< T >::value(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);      
    IgQtLock ();
    rep->clear ();

    if (! m_hits.empty ())
    {
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (0.0, 0.0, 1.0);
        mat->emissiveColor.setValue(0.0, 0.0, 1.0);

        SoDrawStyle *drawStyle = new SoDrawStyle;
        drawStyle->pointSize = 6.0;

        SoSeparator *sep = new SoSeparator; // global
        sep->addChild (mat);
        sep->addChild (drawStyle);

        VisEventSetupService *esService = VisEventSetupService::get (state ());
        ASSERT (esService);

        try 
        {
            typedef std::map<int, std::vector <SbVec3f> > VisParticleMap;
            VisParticleMap pmap;

            for (std::vector<PSimHit>::const_iterator i = m_hits.begin (); i != m_hits.end (); ++i) 
            {
                int ptype = (*i).trackId ();
                
                VisParticleMap::iterator mit = pmap.find (ptype);
                
                if (mit == pmap.end ())
                {
                    std::vector<SbVec3f> vec;
                    bool done = pmap.insert (VisParticleMap::value_type (ptype, vec)).second;
                    if (done)
                    {
                        mit = pmap.find (ptype);
                    }
                }

                if (mit == pmap.end ())
                {
                    throw std::runtime_error ("Cannot find particle type in the map!");
                }               
                DetId id = (*i).detUnitId ();
                if (! id.null () && (id.det () == DetId::Tracker || id.det () == DetId::Muon))
                {
                    GlobalPoint pos = esService->toGlobal ((*i).localPosition (), id);
                    (mit->second).push_back (SbVec3f (pos.x () / 100.0, pos.y () / 100.0, pos.z () / 100.0));
                }
            }       
            if (! pmap.empty ())
            {
                for (VisParticleMap::const_iterator mit = pmap.begin (); mit != pmap.end (); ++mit)
                {
                    SoSeparator *trackSep = new SoSeparator; 
                    
                    QString label = QString ("ID = %1")
                                    .arg (mit->first);             

                    trackSep->setName (Ig3DBaseModel::encode (label.prepend ("Track with ")));

                    SoVertexProperty *vertices = new SoVertexProperty;
                    int nVrtx = 0;

                    for (std::vector<SbVec3f>::const_iterator vit = (mit->second).begin (); vit != (mit->second).end (); ++vit) 
                    {   
                        vertices->vertex.set1Value (nVrtx++, (*vit));
                    }
                    
                    vertices->vertex.setNum (nVrtx);
                    SoPointSet *points = new SoPointSet;
                    points->vertexProperty.setValue (vertices);
                    points->numPoints.setValue (nVrtx);
                    trackSep->addChild (points);

                    if (m_connect.value ())
                    {                   
                        SoLineSet *track = new SoLineSet;
                        track->vertexProperty.setValue (vertices);
                        track->numVertices.setValue (nVrtx);
                        trackSep->addChild (track);
                    }
                    
                    sep->addChild (trackSep);               
                }
            }
        }
        catch (cms::Exception& e)
        {
            e.append (" from VisPSimHitTwig::update(Ig3DRep*) ");
            e.append (this->name ());    

            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 VisPSimHitTwig::update

(

IgTextRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 165 of file VisPSimHitTwig.cc.

References e, exception, i, m_hits, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_text, IgTextRep::setText(), and VisQueuedTwig::update().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);

    // Prepare property description.
    std::ostringstream  text;
    text << "Total " << m_hits.size () << " hits from ";
    text << m_text << "<br>";

    text << "<br>";

    text << "<table width='100%' border=1>"
         << "<TR align = center>"
         << "<TH>Number</TH>"
         << "<TH>Position</TH>"
         << "<TH>Energy loss</TH>"
         << "<TH>DetID</TH>"
         << "<TH>Track ID</TH>"
         << "<TH>Time of flight</TH>"
         << "<TH>Particle type</TH>"
         << "<TH>Process type</TH>"
         << "</TR>";
    text << setiosflags (std::ios::showpoint | std::ios::fixed);
    text.setf (std::ios::right, std::ios::adjustfield);

    if ((! m_hits.empty ())) 
    {
        int nHits = 0;
        try 
        {
            for (std::vector<PSimHit>::const_iterator i = m_hits.begin (); i != m_hits.end (); ++i) 
            {
                text << "<TR align = right>"
                     << "<TD>" << std::setw (3) << nHits++ << "</TD>"
                     << "<TD>" << (*i).localPosition ().x () << ", " 
                               << (*i).localPosition ().y () << ", " 
                     << (*i).localPosition ().z () << "; Phi "<< (*i).localPosition ().phi ()  << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).energyLoss () << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).detUnitId () << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).trackId () << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).timeOfFlight () << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).particleType () << "</TD>"
                     << "<TD>" << std::setw (6) << std::setprecision (3) << (*i).processType () << "</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 ());
}

---

Member Data Documentation

VisSimHitConnect VisPSimHitTwig::m_connect [private]

Definition at line 49 of file VisPSimHitTwig.h.

Referenced by update().

std::string VisPSimHitTwig::m_friendlyName [private]

Definition at line 45 of file VisPSimHitTwig.h.

Referenced by onNewEvent().

std::vector<PSimHit> VisPSimHitTwig::m_hits [private]

Definition at line 50 of file VisPSimHitTwig.h.

Referenced by onNewEvent(), and update().

std::string VisPSimHitTwig::m_instanceName [private]

Definition at line 47 of file VisPSimHitTwig.h.

Referenced by onNewEvent().

std::string VisPSimHitTwig::m_moduleLabel [private]

Definition at line 46 of file VisPSimHitTwig.h.

Referenced by onNewEvent().

std::string VisPSimHitTwig::m_processName [private]

Definition at line 48 of file VisPSimHitTwig.h.

Referenced by onNewEvent().

std::string VisPSimHitTwig::m_text [private]

Definition at line 44 of file VisPSimHitTwig.h.

Referenced by onNewEvent(), and update().

---

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

• VisSimulation/VisSimHits/interface/VisPSimHitTwig.h
• VisSimulation/VisSimHits/src/VisPSimHitTwig.cc

---