---

VisBasicClusterCollectionTwig Class Reference

#include <VisReco/VisEcal/interface/VisBasicClusterCollectionTwig.h>

Inheritance diagram for VisBasicClusterCollectionTwig:

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)
	VisBasicClusterCollectionTwig (IgState *state, IgTwig *parent, const std::string &name="", const std::string &friendlyName="", const std::string &moduleLabel="", const std::string &instanceName="", const std::string &processName="", unsigned color=0xff009c00)
Private Attributes
VisCaloAnnotation	m_annotation
VisBinningStrategy	m_binning
std::vector< EcalRecHit >	m_clusterRecHits
std::vector< reco::BasicCluster >	m_clusters
VisEnergyCut	m_cut
VisCaloEnergyScale	m_escale
const std::string	m_friendlyName
const std::string	m_instanceName
const std::string	m_moduleLabel
std::string	m_name
const std::string	m_processName
unsigned	m_rgba
VisLegoScale	m_scale
std::string	m_text

---

Detailed Description

Definition at line 28 of file VisBasicClusterCollectionTwig.h.

---

Constructor & Destructor Documentation

VisBasicClusterCollectionTwig::VisBasicClusterCollectionTwig	(	IgState *	state,
		IgTwig *	parent,
		const std::string &	name = "",
		const std::string &	friendlyName = "",
		const std::string &	moduleLabel = "",
		const std::string &	instanceName = "",
		const std::string &	processName = "",
		unsigned	color = 0xff009c00
	)

Definition at line 71 of file VisBasicClusterCollectionTwig.cc.

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

    : VisQueuedTwig (state, parent, name),
      m_name (name),
      m_text (name),
      m_friendlyName (friendlyName),
      m_moduleLabel (moduleLabel),
      m_instanceName (instanceName),
      m_processName (processName),
      m_binning (state, lat::CreateCallback (this, &VisBasicClusterCollectionTwig::twigChanged)),
      m_scale (state, lat::CreateCallback (this, &VisBasicClusterCollectionTwig::twigChanged)),
      m_cut (state, lat::CreateCallback (this, &VisBasicClusterCollectionTwig::twigChanged)),
      m_escale (state, lat::CreateCallback (this, &VisBasicClusterCollectionTwig::twigChanged)),
      m_annotation (state, lat::CreateCallback (this, &VisBasicClusterCollectionTwig::twigChanged)),
      m_rgba (colour)
{    
    VisTwigFactroyService *tfService = VisTwigFactroyService::get (state);
    if (! tfService)
    {
        tfService = new VisTwigFactroyService (state);
    }   
    edm::TypeID basicClusterID (typeid (reco::BasicClusterCollection));
    tfService->registerTwig (basicClusterID.friendlyClassName (), &createThisTwig);
}

---

Member Function Documentation

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

Reimplemented from VisQueuedTwig.

Definition at line 106 of file VisBasicClusterCollectionTwig.cc.

References cms::Exception::append(), arg, c, IgCompoundTwig::children(), IgCompoundTwig::destroy(), e, edm::SortedCollection< T, SORT >::end(), exception, edm::SortedCollection< T, SORT >::find(), DBSPlugin::get(), i, edm::Event::id(), IgRepSet::invalidate(), j, m_clusterRecHits, m_clusters, m_friendlyName, m_instanceName, m_moduleLabel, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_processName, m_text, IgSimpleTwig::name(), VisQueuedTwig::onNewEvent(), IgTwig::SELF_MASK, VisQueuedTwig::state(), and IgTwig::STRUCTURE_MASK.

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

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

    std::vector<edm::Handle<reco::BasicClusterCollection> > clusters;

    try 
    {
        QApplication::setOverrideCursor (Qt::waitCursor);
        if (IgQtAppStatusBarService *sbar = IgQtAppStatusBarService::get (state ()))
            sbar->setMessage ("Retrieving a BasicClusterCollection...");

        if ((! m_friendlyName.empty ()) || (! m_moduleLabel.empty ()) || (! m_instanceName.empty ()) || (! m_processName.empty ()))
        {
            if (children ()) destroy ();
            
            VisEventSelector visSel (m_friendlyName, m_moduleLabel, m_instanceName, m_processName);
            event.getMany (visSel, clusters);
        }

        if (IgQtAppStatusBarService *sbar = IgQtAppStatusBarService::get (state ()))
            sbar->setMessage ("BasicClusterCollection done.");
        QApplication::restoreOverrideCursor ();
    }
    catch (cms::Exception& e)
    {
        e.append (" from VisBasicClusterCollectionTwig::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 ();
    }
    
    std::vector< edm::Handle<EBRecHitCollection> > EBRecHits;
    std::vector< edm::Handle<EERecHitCollection> > EERecHits;
    
    try 
    {
        event.getManyByType (EBRecHits);
        event.getManyByType (EERecHits);
    }
    catch (cms::Exception& e)
    {
        e.append (" from VisBasicClusterCollectionTwig::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 ();
    }
    m_clusters.clear ();
    
    if (! clusters.empty ()) 
    {
        std::vector<edm::Handle<reco::BasicClusterCollection> >::iterator i;
        std::vector<edm::Handle<reco::BasicClusterCollection> >::iterator iEnd;
        for (i = clusters.begin (), iEnd = clusters.end (); i != iEnd; ++i) 
        {
            const reco::BasicClusterCollection& 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 (reco::BasicClusterCollection::const_iterator j = c.begin (), jEnd = c.end (); j != jEnd; ++j) 
            {
                m_clusters.push_back (*j);
                QString name = QString ("Basic Cluster Et = %1 GeV")
                               .arg ((*j).energy (), 0, 'f', 2);
                                
                std::vector<DetId> clusterDetIds = (*j).getHitsByDetId ();
                std::vector<DetId>::iterator posCurrent;
                std::vector<DetId>::iterator posEnd;
                m_clusterRecHits.clear ();

                for (posCurrent = clusterDetIds.begin (), posEnd = clusterDetIds.end (); posCurrent != posEnd; ++posCurrent)
                {
                    int EcalNum = (*posCurrent).subdetId ();
                    if (EcalNum == 1) 
                    {
                        for (std::vector<edm::Handle<EBRecHitCollection> >::iterator ieb = EBRecHits.begin (), iebEnd = EBRecHits.end (); 
                             ieb != iebEnd; ++ieb) 
                        {
                            const EBRecHitCollection& ceb = *(*ieb);
                            EBDetId EcalID = (*posCurrent);
                            EBRecHitCollection::const_iterator theRecHit = ceb.find (EcalID);
                            if (theRecHit != ceb.end ())
                            {               
                                m_clusterRecHits.push_back (*theRecHit);
                            }
                        }
                    }               
                    else if (EcalNum == 2) 
                    {
                        for (std::vector<edm::Handle<EERecHitCollection> >::iterator iee = EERecHits.begin (), ieeEnd = EERecHits.end (); 
                             iee != ieeEnd; ++iee) 
                        {
                            const EERecHitCollection& cee = *(*iee);
                            EEDetId EcalID = (*posCurrent);
                            EERecHitCollection::const_iterator theRecHit = cee.find (EcalID);         
                            if (theRecHit != cee.end ())
                            {
                                m_clusterRecHits.push_back (*theRecHit);
                            }
                        }
                    }
                }               
                
                new VisBasicClusterTwig (state (), this, name.latin1 (), m_clusterRecHits);
            }
        }
    }
    
    IgRepSet::invalidate (this->parent (), IgTwig::SELF_MASK | IgTwig::STRUCTURE_MASK);   

    if (IgQtAppStatusBarService *sbar = IgQtAppStatusBarService::get (state ()))
        sbar->setMessage (QString (m_text));
}

void VisBasicClusterCollectionTwig::twigChanged

(

void

)

[virtual]

Definition at line 102 of file VisBasicClusterCollectionTwig.cc.

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

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

void VisBasicClusterCollectionTwig::update

(

IgRZRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 676 of file VisBasicClusterCollectionTwig.cc.

References cms::Exception::append(), ASSERT, IgSoRectColHist::barrelMaxEta, IgSoRectColHist::beamPipeTheta, TestMuL1L2Filter_cff::cerr, Ig3DBaseRep::clear(), e, IgSoRectColHist::endcapMaxTheta, lat::endl(), IgSoRectColHist::energies, relval_parameters_module::energy, exception, IgSoRectColHist::faceColors, DBSPlugin::get(), VisEventSetupService::getBinNumber(), i, IgSoRectColHist::layer, IgSoRectColHist::logScale, m_binning, m_clusters, m_cut, m_escale, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, Ig3DBaseRep::node(), IgSoRectColHist::numR, IgSoRectColHist::numZ, IgSoRectColHist::radiusR, IgSoRectColHist::radiusZ, IgSoRectColHist::scaleFactor, VisQueuedTwig::state(), VisQueuedTwig::update(), and VisActiveConfigurable< T >::value().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
        
    IgQtLock ();
        
    rep->clear ();
        
    SoSeparator *sep = new SoSeparator;    
    SoMaterial *mat = new SoMaterial;
    mat->ambientColor.setValue (1.0, 0.0, 153.0 / 255.0);
    mat->diffuseColor.setValue (1.0, 0.0, 153.0 / 255.0);
    mat->specularColor.setValue (1.0, 0.0, 153.0 / 255.0);
    mat->emissiveColor.setValue (1.0, 0.0, 153.0 / 255.0);
    sep->addChild (mat);    

    VisEventSetupService *esService = VisEventSetupService::get (state ());
    ASSERT (esService);
        
    unsigned binNumber = 0;
    const unsigned nBinsZ = 85*2; // bins along the Z axis (one side only, 85 = a quarter of bins in Z direction)
    const unsigned nBinsR = 39*2;
    const unsigned nbrOfBins = 2 * (nBinsZ + nBinsR);
    const float radiusR = 1.29f;
    const float radiusZ = 3.17f;
    const float barrelMaxEta = 1.479f;                  // max eta of ECAL Barrel
    const float beamPipeTheta = 5.7f * M_PI / 180.f;    // opening angle of the beam gap
    const float endcapMaxTheta = 25.464f * M_PI / 180.f;// max theta of ECAL endcaps
        
    std::vector<float> energies (nbrOfBins * 2); // times 2 for the other energy (which is not used)    
    try 
    {
        if (! m_clusters.empty ()) 
        {
            for (std::vector<reco::BasicCluster>::const_iterator i = m_clusters.begin (), iEnd = m_clusters.end (); i != iEnd; ++i) 
            {
                float energy = (*i).energy ();
                if (energy > m_cut.value ())
                {
                    std::vector<DetId> clusterDetIds = (*i).getHitsByDetId ();
                    std::vector<DetId>::iterator posCurrent;
                    std::vector<DetId>::iterator posEnd;
                                        
                    for (posCurrent = clusterDetIds.begin (), posEnd = clusterDetIds.end (); posCurrent != posEnd; ++posCurrent)
                    {
                        binNumber = esService->getBinNumber (*posCurrent, nBinsR, nBinsZ, radiusR, radiusZ, barrelMaxEta, beamPipeTheta, endcapMaxTheta);

                        if ((binNumber >= 0) && (binNumber < nbrOfBins)) 
                        {
                            if (m_binning.value () == "project")
                                (energies[binNumber] > energy) ? energies[binNumber] : energies[binNumber] = energy;
                            else
                                energies[binNumber] += energy;
                        }
                        else
                        {
                            std::cerr << "Bin number (" << binNumber << ") exceeds the number of bins (" <<  nbrOfBins << ")"<< std::endl;
                        }
                    }
                }
            }
                        
            if (*max_element (energies.begin (), energies.end ()) > 0.0)
            {
                IgSoRectColHist* hist = new IgSoRectColHist;
                hist->radiusR = radiusR;
                hist->radiusZ = radiusZ;
                hist->numR = nBinsR;
                hist->numZ = nBinsZ;
                hist->energies.setValues (0, energies.size(), &energies [0]);
                hist->logScale = false;
                hist->layer = -5.0;
                hist->scaleFactor = m_escale.value ();
                hist->barrelMaxEta = barrelMaxEta;
                hist->beamPipeTheta = beamPipeTheta;
                hist->endcapMaxTheta = endcapMaxTheta;
                                
                std::vector<SbColor> colors (4);
                colors[0] = SbColor (1.0, 0.0, 153.0 / 255.0); // first energy's positive color 
                colors[1] = SbColor (1.0, 1.0, 153.0 / 255.0); // first energy's negative color
                colors[2] = SbColor (51.0 / 255.0, 1.0, 102.0 / 255.0); // second energy's positive color (not used)
                colors[3] = SbColor (51.0 / 255.0, 0.0, 102.0 / 255.0); // second energy's negative color (not used)
                                
                hist->faceColors.setValues (0, colors.size (), &colors[0]);
                sep->addChild (hist);
            }
        }
    }
    catch (cms::Exception& e)
    {
        e.append (" from VisBasicClusterCollectionTwig::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 VisBasicClusterCollectionTwig::update

(

IgRPhiRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 538 of file VisBasicClusterCollectionTwig.cc.

References cms::Exception::append(), ASSERT, Ig3DBaseRep::clear(), e, IgSoCircularHist::energies, relval_parameters_module::energy, exception, first, DBSPlugin::get(), VisEventSetupService::getCellPosition(), i, int, prof2calltree::last, IgSoCircularHist::layer, IgSoCircularHist::logScale, m_annotation, m_binning, m_clusters, m_cut, m_escale, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, IgSoCircularHist::maxRadius, IgSoCircularHist::minRadius, Ig3DBaseRep::node(), IgSoCircularHist::numberOfBins, PV3DBase< T, PVType, FrameType >::phi(), phi, radius(), IgSoCircularHist::scaleFactor, IgSoCircularHist::showAnnotations, VisQueuedTwig::state(), VisQueuedTwig::update(), and VisActiveConfigurable< T >::value().

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

    IgQtLock ();

    rep->clear ();
    
    SoSeparator *sep = new SoSeparator;    
    SoSeparator *posSep = new SoSeparator;    
    SoSeparator *negSep = new SoSeparator;
    sep->addChild (posSep);
    sep->addChild (negSep);
    
    SoMaterial *mat = new SoMaterial;
    mat->ambientColor.setValue (1.0, 0.0, 153.0 / 255.0);
    mat->diffuseColor.setValue (1.0, 0.0, 153.0 / 255.0);
    mat->specularColor.setValue (1.0, 0.0, 153.0 / 255.0);
    mat->emissiveColor.setValue (1.0, 0.0, 153.0 / 255.0);
    posSep->addChild (mat);
    
    SoMaterial *negMat = new SoMaterial;
    negMat->ambientColor.setValue (51.0 / 255.0, 1.0, 102.0 / 255.0);
    negMat->diffuseColor.setValue (51.0 / 255.0, 1.0, 102.0 / 255.0);
    negMat->specularColor.setValue (51.0 / 255.0, 1.0, 102.0 / 255.0);
    negMat->emissiveColor.setValue (51.0 / 255.0, 1.0, 102.0 / 255.0);
    negSep->addChild (negMat);

    int nbrOfBins = 360;
    int binNumber;

    std::vector<float> bufferPositive (nbrOfBins);
    std::vector<float> bufferNegative (nbrOfBins);
    try 
    {
        if (! m_clusters.empty ()) 
        {
            VisEventSetupService *esService = VisEventSetupService::get (state ());
            ASSERT (esService);

            float radius = 1.29;
            
            for (std::vector<reco::BasicCluster>::const_iterator i = m_clusters.begin (), iEnd = m_clusters.end (); i != iEnd; ++i) 
            {
                float energy = (*i).energy ();
                if (energy > m_cut.value ())
                {
                    std::vector<DetId> clusterDetIds = (*i).getHitsByDetId ();
                    std::vector<DetId>::iterator posCurrent;
                    std::vector<DetId>::iterator posEnd;
                    
                    for (posCurrent = clusterDetIds.begin (), posEnd = clusterDetIds.end (); posCurrent != posEnd; ++posCurrent)
                    {
                        const GlobalPoint& pos = esService->getCellPosition (*posCurrent);              
                        float phi = pos.phi ();

                        (phi < 0) ? phi = 2 * M_PI + phi : phi;
       
                        binNumber = (int) floor (phi / (2 * M_PI / nbrOfBins));
                        if (binNumber >= 0 && binNumber < nbrOfBins) 
                            if (m_binning.value () == "project")
                            {                           
                                if ((energy > 0.0) && (energy > bufferPositive [binNumber]))
                                    bufferPositive [binNumber] = energy;
                                else if ((energy < 0.0) && (energy < bufferNegative [binNumber]))
                                    bufferNegative [binNumber] = energy;
                            }
                            else
                                (energy > 0.0) ?  bufferPositive [binNumber] += energy : bufferNegative [binNumber] += energy;
                    }
                }
            }

            std::vector<float>::const_iterator first = bufferPositive.begin ();
            std::vector<float>::const_iterator last = bufferPositive.end ();
            
            if (*max_element (first, last) > 0.0)
            {
                IgSoCircularHist *recHits = new IgSoCircularHist;
                recHits->minRadius = radius;
                recHits->maxRadius = -1;
                recHits->scaleFactor = m_escale.value ();
                recHits->numberOfBins = nbrOfBins;
                recHits->energies.setValues (0, nbrOfBins, &bufferPositive [0]);
                recHits->logScale = false;
                recHits->showAnnotations = m_annotation.value ();
                recHits->layer = 5.5;
                posSep->addChild (recHits);
            }
            
            std::vector<float>::const_iterator nfirst = bufferNegative.begin ();
            std::vector<float>::const_iterator nlast = bufferNegative.end ();
            
            if (fabs(*max_element (nfirst, nlast)) > 0.0)
            {
                IgSoCircularHist *negRecHits = new IgSoCircularHist;
                negRecHits->minRadius = radius;
                negRecHits->maxRadius = -1;
                negRecHits->scaleFactor = m_escale.value ();
                negRecHits->numberOfBins = nbrOfBins;
                negRecHits->energies.setValues (0, nbrOfBins, &bufferNegative [0]);
                negRecHits->logScale = false;
                negRecHits->showAnnotations = m_annotation.value ();
                negRecHits->layer = 5.0;

                negSep->addChild (negRecHits);
            }
        }
    }
    catch (cms::Exception& e)
    {
        e.append (" from VisBasicClusterCollectionTwig::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 VisBasicClusterCollectionTwig::update

(

IgLegoRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 462 of file VisBasicClusterCollectionTwig.cc.

References cms::Exception::append(), ASSERT, Ig3DBaseRep::clear(), e, IgSoTower::emFraction, IgSoTower::energy, relval_parameters_module::energy, eta, PV3DBase< T, PVType, FrameType >::eta(), IgSoTower::etaWidth, exception, DBSPlugin::get(), VisEventSetupService::getCellPosition(), i, m_clusters, m_cut, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_scale, Ig3DBaseRep::node(), PV3DBase< T, PVType, FrameType >::phi(), phi, IgSoTower::phiWidth, IgSoTower::position, IgSoTower::scaleFactor, VisQueuedTwig::state(), VisQueuedTwig::update(), and VisActiveConfigurable< T >::value().

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

    IgQtLock ();

    rep->clear ();
    
    SoSeparator *sep = new SoSeparator;    

    try 
    {
        if (! m_clusters.empty ()) 
        {
            VisEventSetupService *esService = VisEventSetupService::get (state ());
            ASSERT (esService);

            for (std::vector<reco::BasicCluster>::const_iterator i = m_clusters.begin (), iEnd = m_clusters.end (); i != iEnd; ++i) 
            {
                float energy = (*i).energy ();
                if (fabs (energy) > m_cut.value ()) 
                {
                    std::vector<DetId> clusterDetIds = (*i).getHitsByDetId ();
                    std::vector<DetId>::iterator posCurrent;
                    std::vector<DetId>::iterator posEnd;
                    
                    for (posCurrent = clusterDetIds.begin (), posEnd = clusterDetIds.end (); posCurrent != posEnd; ++posCurrent)
                    {
                        const GlobalPoint& pos = esService->getCellPosition (*posCurrent);
                
                        float eta = pos.eta ();
                        float phi = pos.phi ();
                        (phi < 0.0) ? phi = 2 * M_PI + phi : phi; // correction in absence of one convention

                        IgSoTower *tower = new IgSoTower;
                        tower->position = SbVec2f (phi, eta); // eta, phi
                        tower->scaleFactor = m_scale.value ();
                        tower->etaWidth = 0.0174;
                        tower->phiWidth = 0.0174;
                        tower->energy = energy; 
                        tower->emFraction = 1;
                        sep->addChild (tower);
                    }
                }
            }
        } 
    }
    catch (cms::Exception& e)
    {
        e.append (" from VisBasicClusterCollectionTwig::update(IgLegoRep *) ");
        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 VisBasicClusterCollectionTwig::update

(

Ig3DRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 354 of file VisBasicClusterCollectionTwig.cc.

References cms::Exception::append(), ASSERT, Ig3DBaseRep::clear(), e, Ig3DBaseModel::encode(), relval_parameters_module::energy, exception, DBSPlugin::get(), VisEventSetupService::getCorners(), i, m_clusters, m_cut, m_escale, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, Ig3DBaseRep::node(), EZArrayFL< T >::size(), VisQueuedTwig::state(), IgSbColorMap::unpack(), VisQueuedTwig::update(), VisActiveConfigurable< T >::value(), x, y, and z.

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

    rep->clear ();
    
    SoSeparator *clusterSep = new SoSeparator;

    SoMaterial *clusterMaterial = new SoMaterial;
    clusterMaterial->ambientColor.setValue (0.000, 0.000, 0.000); 
    clusterMaterial->diffuseColor.setValue (1.000, 0.000, 0.000);  
    clusterMaterial->emissiveColor.setValue (1.000, 0.000, 0.000); 
    clusterMaterial->specularColor.setValue (0.000, 0.000, 0.000); 
    clusterMaterial->shininess = 0.0;
    clusterMaterial->transparency = 0.0;
    clusterSep->addChild (clusterMaterial);

    if (! m_clusters.empty ())
    {
        try 
        {
            unsigned colour [10] = { 0xff000000, 0xffff0000, 0x00ff0000, // red yellow green
                                     0xff00ff00, 0x00ffff00, 0xd0dfff00, // magenta cyan
                                     0xff000fff, 0x00ffff00, 0xe00fff00,
                                     0xff00ee00};
            int mi = 0;
            
            for (std::vector<reco::BasicCluster>::const_iterator i = m_clusters.begin (), iEnd = m_clusters.end (); i != iEnd; ++i) 
            {
                float energy = (*i).energy ();
                if (energy > m_cut.value ()) 
                {
                    std::vector<DetId> clusterDetIds = (*i).getHitsByDetId ();
                    std::vector<DetId>::iterator posCurrent;
                    std::vector<DetId>::iterator posEnd;

                    SoSeparator *baseSep = new SoSeparator;
                    clusterSep->addChild (baseSep);
                    SoMaterial *baseMat = new SoMaterial;
                    baseMat->setName (Ig3DBaseModel::encode ("Cluster color"));
                    float rgbcomponents [4];
                    IgSbColorMap::unpack (colour [mi++], rgbcomponents);
                    baseMat->diffuseColor.setValue (rgbcomponents);
                    baseSep->addChild (baseMat);
                    
                    VisEventSetupService *esService = VisEventSetupService::get (state ());
                    ASSERT (esService);
                    
                    for (posCurrent = clusterDetIds.begin (), posEnd = clusterDetIds.end (); posCurrent != posEnd; ++posCurrent)
                    {
                        float cEnergy = 0;
                        const CaloCellGeometry::CornersVec& corners = esService->getCorners (*posCurrent);
                        ASSERT (corners.size () == 8);
                    
                        IgSoCrystalHit *crystalHit = new IgSoCrystalHit;
                        crystalHit->energy.setValue (cEnergy);              
                        crystalHit->scale.setValue (m_escale.value ());
                        crystalHit->relativeWidth.setValue (1.0);
                        crystalHit->drawCrystal.setValue (true);
                        crystalHit->drawHit.setValue (true);
                    
                        crystalHit->front1.setValue (corners [3].x () / 100.0, corners [3].y () / 100.0, corners [3].z () / 100.0);
                        crystalHit->front2.setValue (corners [2].x () / 100.0, corners [2].y () / 100.0, corners [2].z () / 100.0);
                        crystalHit->front3.setValue (corners [1].x () / 100.0, corners [1].y () / 100.0, corners [1].z () / 100.0);
                        crystalHit->front4.setValue (corners [0].x () / 100.0, corners [0].y () / 100.0, corners [0].z () / 100.0);
                    
                        crystalHit->back1.setValue  (corners [7].x () / 100.0, corners [7].y () / 100.0, corners [7].z () / 100.0);
                        crystalHit->back2.setValue  (corners [6].x () / 100.0, corners [6].y () / 100.0, corners [6].z () / 100.0);
                        crystalHit->back3.setValue  (corners [5].x () / 100.0, corners [5].y () / 100.0, corners [5].z () / 100.0);
                        crystalHit->back4.setValue  (corners [4].x () / 100.0, corners [4].y () / 100.0, corners [4].z () / 100.0);
                    
                        baseSep->addChild (crystalHit); 
                    }
                }
            }
        }
        catch (cms::Exception& e)
        {
            e.append (" from VisBasicClusterCollectionTwig::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 (clusterSep);
}

void VisBasicClusterCollectionTwig::update

(

IgTextRep *

rep

)

[virtual]

Reimplemented from VisQueuedTwig.

Definition at line 263 of file VisBasicClusterCollectionTwig.cc.

References cms::Exception::append(), e, relval_parameters_module::energy, exception, i, m_clusters, m_cut, m_friendlyName, m_instanceName, m_moduleLabel, VisQueuedTwig::m_onCmsException, VisQueuedTwig::m_onError, VisQueuedTwig::m_onException, VisQueuedTwig::m_onUnhandledException, m_processName, m_text, n, IgSimpleTwig::name(), IgTextRep::setText(), VisQueuedTwig::update(), and VisActiveConfigurable< T >::value().

{
    // Get debugging dump.
    VisQueuedTwig::update (rep);
    
    // Prepare property description.
    std::ostringstream  text;
    text << this->name () << " from ";
    text << m_text << "<br>";

    if (m_friendlyName.empty () && m_moduleLabel.empty () && m_instanceName.empty () && m_processName.empty ())
    {
        text << "Please, select the branches of this Twig for more information.";       
    }
    else
    {
        text << "Displayed only above  " << m_cut.value () << " GeV:<br>";
    
        text << "<table width='100%' border=1>"
             << "<TR align = center>"
             << "<TH>Number</TH>"
             << "<TH>Energy</TH>"
             << "<TH>Position (x, y, z)</TH>"
             << "<TH>Algorithm</TH>"
             << "<TH>Chi-squared</TH>"
             << "<TH>Eta</TH>"
             << "<TH>Phi</TH>"
             << "</TR>";
        text << setiosflags (std::ios::showpoint | std::ios::fixed);
        text.setf (std::ios::right, std::ios::adjustfield);
    
        int n = 0;
        if (! m_clusters.empty ()) 
        {
            try 
            {
                for (std::vector<reco::BasicCluster>::const_iterator i = m_clusters.begin (), iEnd = m_clusters.end (); i != iEnd; ++i) 
                {
                    float energy = (*i).energy ();
                    if (energy > m_cut.value ())
                    {
                        text << "<TR align = right>"
                             << "<TD>" << std::setw (3) << n++ << "</TD>"
                             << "<TD>" << std::setw (5) << std::setprecision (3) << energy << "</TD>"
                             << "<TD>" << std::setw (5) << std::setprecision (3) << (*i).x () << ", " << (*i).y () << ", " << (*i).z () << "</TD>";

                        ((*i).algo () == 0) ?
                                       (text << "<TD>" << std::setw (5) << std::setprecision (3) << "island" << "</TD>") :
                                       (text << "<TD>" << std::setw (5) << std::setprecision (3) << "hybrid" << "</TD>");
                    
                        text << "<TD>" << std::setw (5) << std::setprecision (3) << (*i).chi2 () << "</TD>"
                             << "<TD>" << std::setw (5) << std::setprecision (3) << (*i).eta () << "</TD>"
                             << "<TD>" << std::setw (5) << std::setprecision (3) << (*i).phi ()<< "</TD>"
                             << "</TR>";
                    }
                }
            }
            catch (cms::Exception& e)
            {
                e.append (" from VisBasicClusterCollectionTwig::update(IgTextRep*) ");
                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 ();
            }
        }
    
        text << "</table>";
    }
    
    // Send it over.
    IgQtLock ();
    rep->setText (text.str ());
}

---

Member Data Documentation

VisCaloAnnotation VisBasicClusterCollectionTwig::m_annotation [private]

Definition at line 61 of file VisBasicClusterCollectionTwig.h.

Referenced by update().

VisBinningStrategy VisBasicClusterCollectionTwig::m_binning [private]

Definition at line 57 of file VisBasicClusterCollectionTwig.h.

Referenced by update().

std::vector<EcalRecHit> VisBasicClusterCollectionTwig::m_clusterRecHits [private]

Definition at line 63 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent().

std::vector<reco::BasicCluster> VisBasicClusterCollectionTwig::m_clusters [private]

Definition at line 62 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent(), and update().

VisEnergyCut VisBasicClusterCollectionTwig::m_cut [private]

Definition at line 59 of file VisBasicClusterCollectionTwig.h.

Referenced by update().

VisCaloEnergyScale VisBasicClusterCollectionTwig::m_escale [private]

Definition at line 60 of file VisBasicClusterCollectionTwig.h.

Referenced by update().

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

Definition at line 53 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent(), and update().

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

Definition at line 55 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent(), and update().

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

Definition at line 54 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent(), and update().

std::string VisBasicClusterCollectionTwig::m_name [private]

Reimplemented from IgSimpleTwig.

Definition at line 51 of file VisBasicClusterCollectionTwig.h.

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

Definition at line 56 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent(), and update().

unsigned VisBasicClusterCollectionTwig::m_rgba [private]

Definition at line 64 of file VisBasicClusterCollectionTwig.h.

VisLegoScale VisBasicClusterCollectionTwig::m_scale [private]

Definition at line 58 of file VisBasicClusterCollectionTwig.h.

Referenced by update().

std::string VisBasicClusterCollectionTwig::m_text [private]

Definition at line 52 of file VisBasicClusterCollectionTwig.h.

Referenced by onNewEvent(), and update().

---

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

• VisReco/VisEcal/interface/VisBasicClusterCollectionTwig.h
• VisReco/VisEcal/src/VisBasicClusterCollectionTwig.cc

---