---

VisTkSimTrackTwig Class Reference

#include <VisSimulation/VisSimTracker/interface/VisTkSimTrackTwig.h>

Inheritance diagram for VisTkSimTrackTwig:

List of all members.

Public Member Functions
virtual void	onNewEvent (const edm::Event &event, const edm::EventSetup &eventSetup)
virtual void	update (IgRZRep *rep)
virtual void	update (IgRPhiRep *rep)
virtual void	update (IgTextRep *rep)
virtual void	update (IgLegoRep *rep)
virtual void	update (Ig3DRep *rep)
	VisTkSimTrackTwig (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
unsigned int	m_detType
const std::string	m_friendlyName
const std::string	m_instanceName
const std::string	m_moduleLabel
edm::ESHandle< TrackerGeometry >	m_pDD
const std::string	m_processName
unsigned int	m_subDetType
std::string	m_text
std::vector< PSimHit >	m_trackerHits
std::vector< std::vector < PSimHit > >	m_trackerHitsSorted

---

Detailed Description

Definition at line 29 of file VisTkSimTrackTwig.h.

---

Constructor & Destructor Documentation

VisTkSimTrackTwig::VisTkSimTrackTwig	(	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 75 of file VisTkSimTrackTwig.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)
{
    VisTwigFactroyService *tfService = VisTwigFactroyService::get (state);
    if (! tfService)
    {
        tfService = new VisTwigFactroyService (state);
    }   
    edm::TypeID simtrackID (typeid (edm::PSimHitContainer));
    tfService->registerTwig (simtrackID.friendlyClassName (), &createThisTwig);
}

---

Member Function Documentation

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

Reimplemented from VisQueuedTwig.

Definition at line 98 of file VisTkSimTrackTwig.cc.

References arg, c, IgCompoundTwig::clear(), e, exception, edm::EventSetup::get(), i, edm::Event::id(), edm::ESHandle< T >::isValid(), m_friendlyName, m_instanceName, m_moduleLabel, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_pDD, m_processName, m_text, m_trackerHits, m_trackerHitsSorted, VisQueuedTwig::onBaseInvalidate(), VisQueuedTwig::onNewEvent(), python::multivaluedict::sort(), sort_TOF(), and track.

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

    m_trackerHits.clear ();
  
    m_text = (QString ("Run # %1, event # %2")
              .arg (event.id ().run ())
              .arg (event.id ().event ()).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);
        }
        else
        {
            event.getManyByType (hitCollections);
        }
    }
    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 (! hitCollections.empty ())
    {
        std::vector<edm::Handle<edm::PSimHitContainer> >::iterator i;
        for (i = hitCollections.begin (); i != hitCollections.end (); i++) 
        {
            const edm::PSimHitContainer& c = *(*i);
          
            for (std::vector<PSimHit>::const_iterator isim = c.begin (), isimEnd = c.end ();
                 isim != isimEnd; ++isim)
            {
                m_trackerHits.push_back ((*isim));
            }    
        }
    }
    try 
    { 
        eventSetup.get<TrackerDigiGeometryRecord> ().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 (! m_trackerHits.empty () && m_pDD.isValid ())
    {
        try 
        {         
            m_trackerHitsSorted.clear ();
            unsigned int maxid = 0;
            //get number of tracks track 65535 (FFFF) is exluded. I suppose it is no real track.
            for (std::vector<PSimHit>::iterator isim = m_trackerHits.begin (), isimEnd = m_trackerHits.end (); 
                 isim != isimEnd; ++isim)
            {
                if ((maxid < (*isim).trackId ()) && ((*isim).trackId () < 65535)) maxid = (*isim).trackId ();
            }

            //array of vector  is assigned and initialized
            std::vector<PSimHit> *m_trackerHitsUnsorted = new std::vector<PSimHit> [maxid + 2];
            maxid += 2;
            
            for (unsigned int i = 0; i < maxid; ++i)
            {
                (m_trackerHitsUnsorted [i]).clear ();
            }

            //filling of array 

            for (std::vector<PSimHit>::iterator isim = m_trackerHits.begin (), isimEnd = m_trackerHits.end (); 
                 isim != isimEnd; ++isim)
            {
                if ((*isim).trackId () == 65535)
                {
                    (m_trackerHitsUnsorted [0]).push_back (*isim);
                }
                else 
                {
                    (m_trackerHitsUnsorted[((*isim).trackId ()) + 1]).push_back (*isim);
                }             
            }

            //copying data into theTrackerHitsSorted and deleting theTrackerHitsUnsorted
            for (unsigned int i = 0; i < maxid; ++i) 
                m_trackerHitsSorted.push_back (m_trackerHitsUnsorted [i]);
            delete [] m_trackerHitsUnsorted;
            // calling theTrackerHitsSorted[i][j] gives Event j of TrackID i-1 
 
            //Every Track should be sorted by ToF and IDs which do not exist shoul dbe removed
            for (std::vector<std::vector<PSimHit> >::iterator track = m_trackerHitsSorted.begin (), 
                                                   trackEnd = m_trackerHitsSorted.end (); 
                 track != trackEnd; ++track)
            {
                unsigned int length = (*track).size ();
                if (length == 0)
                {
                    m_trackerHitsSorted.erase (track); // Remove all elements. 
                    track--;
                    continue;
                }

                unsigned int id = ((*track) [0]).trackId ();
                if (id == 65535) continue;
                std::sort ((*track).begin (), (*track).end (), sort_TOF);
            }
        }
        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 ();
        }
    }

    VisQueuedTwig::onBaseInvalidate ();
}

void VisTkSimTrackTwig::update

(

IgRZRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 532 of file VisTkSimTrackTwig.cc.

References Ig3DBaseRep::clear(), IgSoSimpleTrajectory::controlPoints, e, exception, SurfaceOrientation::inner, edm::ESHandle< T >::isValid(), VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_pDD, m_trackerHits, m_trackerHitsSorted, IgSoSimpleTrajectory::markerPoints, Ig3DBaseRep::node(), track, VisQueuedTwig::update(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), y, PV3DBase< T, PVType, FrameType >::z(), and z.

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    IgQtLock ();
    rep->clear ();
    if (! m_trackerHits.empty () && m_pDD.isValid ())
    {
        SoSeparator *sep = new SoSeparator; // global
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (1.0, 1.0, 0.0);
        sep->addChild (mat);
        try 
        {
            for (std::vector<std::vector<PSimHit> >::iterator track = m_trackerHitsSorted.begin (),
                                                           trackEnd = m_trackerHitsSorted.end (); 
                 track != trackEnd; ++track)
            {
                unsigned int id = ((*track) [0]).trackId ();
                if (id == 65535) continue;
                IgSoSimpleTrajectory *simulatedTrack = new IgSoSimpleTrajectory;
                int nSimHit = 0;
                for (std::vector<PSimHit>::iterator jsim = (*track).begin (),
                                                 jsimEnd = (*track).end (); 
                     jsim != jsimEnd; ++jsim)
                {
                    GlobalPoint inner = m_pDD->idToDet (DetId ((*jsim).detUnitId ()))->surface ().toGlobal ((*jsim).localPosition ()); 
                    float x = inner.x () / 100;
                    float y = inner.y () / 100;
                    float z = inner.z () / 100;
                    simulatedTrack->controlPoints.set1Value (nSimHit, (SbVec3f (x, y, z)));
                    simulatedTrack->markerPoints.set1Value (nSimHit, (SbVec3f (x, y, z)));            
                    ++nSimHit;
                }
                sep->addChild (simulatedTrack);
                //if (id>30) break;;
                //delete simulatedTrack;
            }
        }
        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 VisTkSimTrackTwig::update

(

IgRPhiRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 468 of file VisTkSimTrackTwig.cc.

References Ig3DBaseRep::clear(), IgSoSimpleTrajectory::controlPoints, e, exception, SurfaceOrientation::inner, edm::ESHandle< T >::isValid(), VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_pDD, m_trackerHits, m_trackerHitsSorted, IgSoSimpleTrajectory::markerPoints, Ig3DBaseRep::node(), track, VisQueuedTwig::update(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), y, PV3DBase< T, PVType, FrameType >::z(), and z.

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

    if (! m_trackerHits.empty () && m_pDD.isValid ())
    {
        SoSeparator *sep = new SoSeparator; // global
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (1.0, 1.0, 0.0);
        sep->addChild (mat);
        try 
        {
            for (std::vector<std::vector<PSimHit> >::iterator track = m_trackerHitsSorted.begin (),
                                                           trackEnd = m_trackerHitsSorted.end (); 
                 track != trackEnd; ++track)
            {
                unsigned int id = ((*track) [0]).trackId ();
                if (id == 65535) continue;
                IgSoSimpleTrajectory *simulatedTrack = new IgSoSimpleTrajectory;
                int nSimHit = 0;
                for (std::vector<PSimHit>::iterator jsim = (*track).begin (), jsimEnd = (*track).end (); 
                     jsim != jsimEnd; ++jsim)
                {
                    GlobalPoint inner = m_pDD->idToDet (DetId ((*jsim).detUnitId ()))->surface ().toGlobal ((*jsim).localPosition ()); 
                    float x = inner.x () / 100;
                    float y = inner.y () / 100;
                    float z = inner.z () / 100;
                    simulatedTrack->controlPoints.set1Value (nSimHit, (SbVec3f (x, y, z)));
                    simulatedTrack->markerPoints.set1Value (nSimHit, (SbVec3f (x, y, z)));            
                    ++nSimHit;
                }
                sep->addChild (simulatedTrack);
                //if (id>30) break;;
                //delete simulatedTrack;
            }
        }
        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 VisTkSimTrackTwig::update

(

IgTextRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 271 of file VisTkSimTrackTwig.cc.

References e, exception, edm::ESHandle< T >::isValid(), VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_pDD, m_text, m_trackerHits, m_trackerHitsSorted, IgTextRep::setText(), and VisQueuedTwig::update().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    // Prepare property description.
    std::ostringstream  text;
    text << "Total " << m_trackerHitsSorted.size () << " Tracker SimTracks from";
    if (m_trackerHitsSorted.size () > 1000) 
    {
        text << " input file: printing only first 1000 " << "<br>"; 
    }
    else 
    {
        text << " input file: " << m_text << "<br>";
    }

    text << "<table width='100%' border=1>"
         << "<TR align = center>"
         << "<TH>TrackID</TH>"
         << "<TH>NumberSimHits</TH>"
         << "</TR>";
    text << setiosflags (std::ios::showpoint | std::ios::fixed);
    text.setf (std::ios::right, std::ios::adjustfield);

    if (! m_trackerHits.empty () && m_pDD.isValid ()) 
    {
        int nSimtracks = 0;
        try 
        {
            for (std::vector<std::vector<PSimHit> >::iterator isim = m_trackerHitsSorted.begin (), 
                                                           isimEnd = m_trackerHitsSorted.end (); 
                 isim != isimEnd; ++isim)
            {
                unsigned int length = (*isim).size ();
                unsigned int id = ((*isim) [0]).trackId ();
        
                text << "<TR align = right>"
                     << "<TD>" << id << "</TD>"
                     << "<TD>" << length << "</TD>";
                ++nSimtracks;

                if (nSimtracks > 1000) break;
            }
        }
        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 VisTkSimTrackTwig::update

(

IgLegoRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 407 of file VisTkSimTrackTwig.cc.

References Ig3DBaseRep::clear(), e, IgSoTower::energy, eta, PV3DBase< T, PVType, FrameType >::eta(), exception, i, iter, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_pDD, m_trackerHitsSorted, Ig3DBaseRep::node(), PV3DBase< T, PVType, FrameType >::phi(), phi, IgSoTower::position, IgSoTower::scaleFactor, and VisQueuedTwig::update().

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

    IgQtLock ();
    rep->clear ();
    
    SoSeparator *sep = new SoSeparator; // global
    
    try 
    {
        for (std::vector<std::vector<PSimHit> >::const_iterator i = m_trackerHitsSorted.begin (),
                                                             iEnd = m_trackerHitsSorted.end (); 
             i != iEnd; ++i) 
        {
            if ((*(*i).begin ()).trackId () != 65535)
            {
                for (std::vector<PSimHit>::const_iterator iter = (*i).begin (), iterEnd = (*i).end (); iter != iterEnd; ++iter) 
                {
                    GlobalPoint SimHitPosition = m_pDD->idToDet (DetId ((*iter).detUnitId ()))->surface ().toGlobal ((*iter).localPosition ());

                    double eta = SimHitPosition.eta ();
                    double phi = SimHitPosition.phi ();
                    if (phi < 0.0) phi = 2 * M_PI + phi; // correction in absence of one convention

                    float EnergyLoss = (*iter).energyLoss(); // the energy deposit in the PSimHit, in ???

                    IgSoTower *tower = new IgSoTower;
                    tower->position = SbVec2f (phi, eta); // eta, phi
                    tower->scaleFactor = 20.0;
                    tower->energy = EnergyLoss;
                    sep->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 VisTkSimTrackTwig::update

(

Ig3DRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 344 of file VisTkSimTrackTwig.cc.

References Ig3DBaseRep::clear(), IgSoSimpleTrajectory::controlPoints, e, exception, SurfaceOrientation::inner, edm::ESHandle< T >::isValid(), VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_pDD, m_trackerHits, m_trackerHitsSorted, IgSoSimpleTrajectory::markerPoints, Ig3DBaseRep::node(), track, VisQueuedTwig::update(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), y, PV3DBase< T, PVType, FrameType >::z(), and z.

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    IgQtLock ();
    rep->clear ();
    if (! m_trackerHits.empty () && m_pDD.isValid ())
    {
        SoSeparator *sep = new SoSeparator; // global
        SoMaterial *mat = new SoMaterial;
        mat->diffuseColor.setValue (1.0, 1.0, 0.0);
        sep->addChild (mat);
        try 
        {
            for (std::vector<std::vector<PSimHit> >::iterator track = m_trackerHitsSorted.begin (),
                                                           trackEnd = m_trackerHitsSorted.end (); 
                 track != trackEnd; ++track)
            {
                unsigned int id = ((*track) [0]).trackId ();
                if (id == 65535) continue;
                IgSoSimpleTrajectory *simulatedTrack = new IgSoSimpleTrajectory;
                int nSimHit = 0;
                for (std::vector<PSimHit>::iterator jsim = (*track).begin (), jsimEnd = (*track).end (); 
                     jsim != jsimEnd; ++jsim)
                {
                    GlobalPoint inner = m_pDD->idToDet (DetId ((*jsim).detUnitId ()))->surface ().toGlobal ((*jsim).localPosition ()); 
                    float x = inner.x () / 100;
                    float y = inner.y () / 100;
                    float z = inner.z () / 100;
                    simulatedTrack->controlPoints.set1Value (nSimHit, (SbVec3f (x, y, z)));
                    simulatedTrack->markerPoints.set1Value (nSimHit, (SbVec3f (x, y, z)));            
                    ++nSimHit;
                }
                sep->addChild (simulatedTrack);
                //if (id>30) break;;
                //delete simulatedTrack;
            }
        }
        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);
    }
}

---

Member Data Documentation

unsigned int VisTkSimTrackTwig::m_detType [private]

Definition at line 57 of file VisTkSimTrackTwig.h.

const std::string VisTkSimTrackTwig::m_friendlyName [private]

Definition at line 52 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent().

const std::string VisTkSimTrackTwig::m_instanceName [private]

Definition at line 54 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent().

const std::string VisTkSimTrackTwig::m_moduleLabel [private]

Definition at line 53 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent().

edm::ESHandle<TrackerGeometry> VisTkSimTrackTwig::m_pDD [private]

Definition at line 56 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent(), and update().

const std::string VisTkSimTrackTwig::m_processName [private]

Definition at line 55 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent().

unsigned int VisTkSimTrackTwig::m_subDetType [private]

Definition at line 58 of file VisTkSimTrackTwig.h.

std::string VisTkSimTrackTwig::m_text [private]

Definition at line 51 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent(), and update().

std::vector<PSimHit> VisTkSimTrackTwig::m_trackerHits [private]

Definition at line 59 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent(), and update().

std::vector<std::vector<PSimHit> > VisTkSimTrackTwig::m_trackerHitsSorted [private]

Definition at line 60 of file VisTkSimTrackTwig.h.

Referenced by onNewEvent(), and update().

---

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

• VisSimulation/VisSimTracker/interface/VisTkSimTrackTwig.h
• VisSimulation/VisSimTracker/src/VisTkSimTrackTwig.cc

---