/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/Fireworks/ParticleFlow/src/FWPFLegoRecHit.cc

Go to the documentation of this file.

#include "Fireworks/ParticleFlow/interface/FWPFLegoRecHit.h"
#include "Fireworks/Core/interface/FWProxyBuilderBase.h"

//______________________________________________________________________________
FWPFLegoRecHit::FWPFLegoRecHit( const std::vector<TEveVector> &corners, TEveElement *comp, FWProxyBuilderBase*pb,
                                const FWViewContext *vc, float e, float et )
   : m_energy(e), m_et(et), m_isTallest(false)
{
   buildTower( corners, vc );
   buildLineSet( corners, vc );

   pb->setupAddElement( m_tower, comp );
   pb->setupAddElement( m_ls, comp );
}

//______________________________________________________________________________
void
FWPFLegoRecHit::setupEveBox( std::vector<TEveVector> &corners, float scale )
{
   for( size_t i = 0; i < 4; ++i )
   {
      int j = i + 4;
      corners[i+4].fZ = corners[i].fZ + scale;
      m_tower->SetVertex( i, corners[i] );
      m_tower->SetVertex( j, corners[j] );
   }

   m_tower->SetPickable( true );
   m_tower->SetDrawFrame(false);
   m_tower->SetLineWidth( 1.0 );
   m_tower->SetLineColor( kBlack );
}

//______________________________________________________________________________
void
FWPFLegoRecHit::buildTower( const std::vector<TEveVector> &corners, const FWViewContext *vc )
{
   m_tower = new TEveBox( "EcalRecHitTower" );
   std::vector<TEveVector> towerCorners = corners;
   FWViewEnergyScale *caloScale = vc->getEnergyScale();
   float val = caloScale->getPlotEt() ? m_et : m_energy;
   float scale = caloScale->getScaleFactorLego() * val;

   if( scale < 0 )
      scale *= -1;

   setupEveBox( towerCorners, scale );
}

//______________________________________________________________________________
void
FWPFLegoRecHit::buildLineSet( const std::vector<TEveVector> &corners, const FWViewContext *vc )
{
   m_ls = new TEveStraightLineSet( "EcalRecHitLineSet" );

   // no need to set anything, all is re-set in updateScales()
   // reserve space for square outline
   TEveVector c;
   m_ls->AddLine( c.fX, c.fY, c.fZ, c.fX, c.fY, c.fZ );
   m_ls->AddLine( c.fX, c.fY, c.fZ, c.fX, c.fY, c.fZ );
   m_ls->AddLine( c.fX, c.fY, c.fZ, c.fX, c.fY, c.fZ );
   m_ls->AddLine( c.fX, c.fY, c.fZ, c.fX, c.fY, c.fZ );

   // last line is trick to add a marker in line set
   m_ls->SetMarkerStyle( 1 );
   m_ls->AddLine( c.fX, c.fY, c.fZ, c.fX, c.fY, c.fZ );
   m_ls->AddMarker( 0, 0. );
}

//______________________________________________________________________________
void
FWPFLegoRecHit::updateScale( const FWViewContext *vc, float maxLogVal )
{
   FWViewEnergyScale *caloScale = vc->getEnergyScale();
   float val = caloScale->getPlotEt() ? m_et : m_energy;
   float scale = caloScale->getScaleFactorLego() * val;

   // printf("scale %f %f\n",  caloScale->getValToHeight(), val);

   if( scale < 0 )
      scale *= -1;

   // Reposition top points of tower
   const float *data;
   TEveVector c;
   for( unsigned int i = 0; i < 4; ++i )
   {
      data = m_tower->GetVertex( i );
      c.fX += data[0];
      c.fY += data[1];
      m_tower->SetVertex( i, data[0], data[1], 0 );
      m_tower->SetVertex( i+4,  data[0], data[1], scale);
   }
   c *= 0.25;
   // Scale lineset 
   float s = log( 1 + val ) / maxLogVal;
   float d = 0.5 * ( m_tower->GetVertex(1)[0]  -m_tower->GetVertex(0)[0]);
   d *= s;
   float z =  scale * 1.001;
   setLine(0, c.fX - d, c.fY -d, z, c.fX + d, c.fY -d, z);
   setLine(1, c.fX + d, c.fY -d, z, c.fX + d, c.fY +d, z);
   setLine(2, c.fX + d, c.fY +d, z, c.fX - d, c.fY +d, z);
   setLine(3, c.fX - d, c.fY +d, z, c.fX - d, c.fY -d, z);

   if( m_isTallest )
   {
      // This is the tallest tower and hence two additional lines needs scaling
      setLine( 4, c.fX - d, c.fY - d, z, c.fX + d, c.fY + d, z );
      setLine( 5, c.fX - d, c.fY + d, z, c.fX + d, c.fY - d, z );
   }

   TEveStraightLineSet::Marker_t* m = ((TEveStraightLineSet::Marker_t*)(m_ls->GetMarkerPlex().Atom(0)));
   m->fV[0] = c.fX; m->fV[1] = c.fY; m->fV[2] = z;

   // stamp changed elements
   m_tower->StampTransBBox();
   m_ls->StampTransBBox();
}

//______________________________________________________________________________
void FWPFLegoRecHit::setLine(int idx, float x1, float y1, float z1, float x2, float y2, float z2)
{
   // AMT: this func should go in TEveStraightLineSet class

   TEveStraightLineSet::Line_t* l = ((TEveStraightLineSet::Line_t*)(m_ls->GetLinePlex().Atom(idx)));

   l->fV1[0] = x1;
   l->fV1[1] = y1;
   l->fV1[2] = z1; 

   l->fV2[0] = x2;
   l->fV2[1] = y2;
   l->fV2[2] = z2;
}

//______________________________________________________________________________
void
FWPFLegoRecHit::setIsTallest( bool b )
{
   m_isTallest = b;
   
   if( m_isTallest )
   {
      TEveVector vec;
      addLine( vec, vec );
      addLine( vec, vec );
   }
}

//______________________________________________________________________________
void
FWPFLegoRecHit::addLine( float x1, float y1, float z1, float x2, float y2, float z2 )
{
   m_ls->AddLine( x1, y1, z1, x2, y2, z2 );
}

//______________________________________________________________________________
void
FWPFLegoRecHit::addLine( const TEveVector &v1, const TEveVector &v2 )
{
   m_ls->AddLine(v1.fX, v1.fY, v1.fZ, v2.fX, v2.fY, v2.fZ);
}