/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/Fireworks/ParticleFlow/plugins/FWPFLegoRecHit.cc

Go to the documentation of this file.

#include "FWPFLegoRecHit.h"
#include "FWPFEcalRecHitLegoProxyBuilder.h"

//______________________________________________________________________________
FWPFLegoRecHit::FWPFLegoRecHit( const std::vector<TEveVector> &corners, TEveElement *comp, FWPFEcalRecHitLegoProxyBuilder *pb,
                                const FWViewContext *vc, float e, float et )
   : m_builder(pb), 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( const std::vector<TEveVector> &corners )
{
   for( size_t i = 0; i < corners.size(); ++i )
      m_tower->SetVertex( i, corners[i] );

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

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

//______________________________________________________________________________
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;

   convertToTower( towerCorners, scale );
   setupEveBox( towerCorners );
}

//______________________________________________________________________________
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 )
{
   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 ) / m_builder->getMaxValLog(caloScale->getPlotEt());
   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);
}