VP1CC_TileBarEc Class Reference

#include <VP1CaloCells.h>

Inheritance diagram for VP1CC_TileBarEc:

Collaboration diagram for VP1CC_TileBarEc:

Public Member Functions
	VP1CC_TileBarEc (const CaloCell *caloCell, const TileID *tile_id, const VP1CC_SeparatorMap *separators)
virtual	~VP1CC_TileBarEc ()
virtual void	build3DObjects (VP1CC_SoNode2CCMap *node2cc, bool useEt, const QPair< bool, double > &scale, bool outline, const VP1CC_GlobalCuts &)
bool	isInsideClipVolume (const VP1CC_GlobalCuts &globalCuts)
virtual void	remove3DObjects (VP1CC_SoNode2CCMap *node2cc)
std::vector< std::string >	ToString (const CaloCell_ID *calo_id, const std::string &extrainfos="")
int	GetFragChannel (const TileHWID *tile_hw_id, bool up, int &frag, int &channel)
double	energyToTransverse (const double &) const
bool	cutPassed (const VP1CC_GlobalCuts &globalCuts)
const CaloCell *	getCaloCell ()
Identifier	getID ()
void	updateScene (VP1CC_SoNode2CCMap *node2cc, bool useEt, const QPair< bool, double > &scale, bool outline, const VP1CC_GlobalCuts &globalCuts)
double	cellDepth (const QPair< bool, double > &scale, const double &energy)

Protected Member Functions
std::string	id2name (Identifier &id)

Protected Attributes
SoGenericBox *	m_hitUp
SoGenericBox *	m_hitDown
VP1ExtraSepLayerHelper *	m_helperUp
VP1ExtraSepLayerHelper *	m_helperDown
const TileID *	m_tileID
const CaloCell *	m_caloCell

Detailed Description

Definition at line 313 of file VP1CaloCells.h.

Constructor & Destructor Documentation

VP1CC_TileBarEc()

VP1CC_TileBarEc::VP1CC_TileBarEc	(	const CaloCell *	caloCell,
		const TileID *	tile_id,
		const VP1CC_SeparatorMap *	separators )

caloCell was already dereferenced in the initialiser list, so this is redundant: you've already crashed. if(!caloCell) throw std::runtime_error("VP1CC_TileBarEc: 0 pointer to CaloCell");

Definition at line 466 of file VP1CaloCells.cxx.

                                                          :
  VP1CC_Tile(caloCell,tile_id),
  m_hitUp(0),
  m_hitDown(0)
{
  if(!separators)
    throw std::runtime_error("VP1CC_TileBarEc: 0 pointer to VP1CC Separator Map");
 
  // The CaloCell has to be of type TileCell
  const TileCell* tile_cell = dynamic_cast<const TileCell*>(caloCell);
  if(tile_cell==0)
    throw std::runtime_error("VP1CC_TileBarEc: CaloCell is not TileCell");
 
  // define Separator types (Up/Down) depending on Cell energy
  VP1CC_SeparatorTypes mySeparatorTypeUp;
  VP1CC_SeparatorTypes mySeparatorTypeDown;
 
  double energy1 = tile_cell->ene1();
  double energy2 = tile_cell->ene2();
 
  if(caloCell->energy()<0)
  {
    if(energy1>=0)
      mySeparatorTypeDown = VP1CC_SepTileNegativePos;
    else
      mySeparatorTypeDown = VP1CC_SepTileNegativeDown;
 
    if(energy2>=0)
      mySeparatorTypeUp = VP1CC_SepTileNegativePos;
    else
      mySeparatorTypeUp = VP1CC_SepTileNegativeUp;
  }
  else
  {
    if(energy1<=0)
      mySeparatorTypeDown = VP1CC_SepTilePositiveNeg;
    else
      mySeparatorTypeDown = VP1CC_SepTilePositiveDown;
 
    if(energy2<=0)
      mySeparatorTypeUp = VP1CC_SepTilePositiveNeg;
    else
      mySeparatorTypeUp = VP1CC_SepTilePositiveUp;
  }
 
  VP1CC_SeparatorMap::const_iterator it = separators->find(mySeparatorTypeDown);
  if(it ==separators->end())
    throw std::runtime_error("VP1CC_TileBarEc: Missing separator helper in the map");
 
  m_helperDown = it->second;
 
  it = separators->find(mySeparatorTypeUp);
  if(it ==separators->end())
    throw std::runtime_error("VP1CC_TileBarEc: Missing separator helper in the map");
 
  m_helperUp = it->second;
 
  if(!m_helperDown || !m_helperUp)
    throw std::runtime_error("VP1CC_TileBarEc: 0 pointer to one of separator helpers" );
}

~VP1CC_TileBarEc()

VP1CC_TileBarEc::~VP1CC_TileBarEc ( )

virtual

Definition at line 534 of file VP1CaloCells.cxx.

{
}

Member Function Documentation

build3DObjects()

void VP1CC_TileBarEc::build3DObjects ( VP1CC_SoNode2CCMap * node2cc, bool useEt, const QPair< bool, double > & scale, bool outline, const VP1CC_GlobalCuts & globalCuts )

virtual

Implements VP1CaloCell.

Definition at line 598 of file VP1CaloCells.cxx.

{
  bool createNewHit = false;
  if(!m_hitUp){
    m_hitUp = new SoGenericBox();
    m_hitDown = new SoGenericBox();
    createNewHit = true;
  }
 
  const TileCell* tile_cell = dynamic_cast<const TileCell*>(m_caloCell);
  if (not tile_cell) return;
  double energy1 = (useEt ? energyToTransverse(tile_cell->ene1()) : tile_cell->ene1());
  double energy2 = (useEt ? energyToTransverse(tile_cell->ene2()) : tile_cell->ene2());
 
  const CaloDetDescrElement* ddElement = m_caloCell->caloDDE();
 
  double eta = ddElement->eta();
  double deta = ddElement->deta()*0.5;
  double phi = ddElement->phi();
  double dphi = ddElement->dphi()*0.5;
  double radius = ddElement->r();
  double dradius = ddElement->dr()*0.5;
 
  double z = ddElement->z();
  double minR = radius - dradius;
  double radialDistance = sqrt(minR*minR + z*z);
 
 
  double depth1 = cellDepth(scale,energy1);
  depth1 = std::min(depth1, std::fabs(globalCuts.clipRadius - radialDistance));
  double depth2 = cellDepth(scale,energy2);
  depth2 = std::min(depth2, std::fabs(globalCuts.clipRadius - radialDistance));
  // std::cout<<"VP1CC_TileBarEc clipRadius: "<<globalCuts.clipRadius<<"\t radius "<<radius<<"\t depth1 "
 //           <<cellDepth(scale,energy1)<<"\t depth2 "<<cellDepth(scale,energy2)<<"\t cd1 "<<depth1<<"\t cd1 "<<depth1<<std::endl;
 
  //double etaMin, double etaMax,
  // double phiMin, double phiMax,
  // double cellDepth, double cellDistance,
  // double etasqueezefact, double phisqueezefact 
  m_hitUp->setParametersForBarrelEtaPhiCell(eta,eta+deta,phi-dphi,phi+dphi,
                        depth1,radius-dradius,0.9,0.9);
  m_hitUp->drawEdgeLines = outline;
  m_hitDown->setParametersForBarrelEtaPhiCell(eta-deta,eta,phi-dphi,phi+dphi,
                          depth2,radius-dradius,0.9,0.9);
  m_hitDown->drawEdgeLines = outline;
 
  if(createNewHit) {
    (*node2cc)[m_hitUp] = this;
    (*node2cc)[m_hitDown] = this;
    m_helperUp->addNode(m_hitUp);
    m_helperDown->addNode(m_hitDown);
  }
}

cellDepth()

double VP1CaloCell::cellDepth ( const QPair< bool, double > & scale, const double & energy )

inlineinherited

Definition at line 189 of file VP1CaloCells.h.

  { return std::max(1.0*Gaudi::Units::mm,scale.second*(scale.first?log(1+fabs(energy)):fabs(energy))); }

cutPassed()

bool VP1CaloCell::cutPassed ( const VP1CC_GlobalCuts & globalCuts )

inherited

Definition at line 63 of file VP1CaloCells.cxx.

{
  // only do this if it's something significant
  if (globalCuts.clipRadius<10e8){
    if (!isInsideClipVolume(globalCuts)) return false;
    // double radius = sqrt(m_caloCell->x()*m_caloCell->x() + m_caloCell->y()*m_caloCell->y() + m_caloCell->z()*m_caloCell->z());
   
  }
  
  // Check side and Eta cut
  bool passed = (m_caloCell->eta()>=0 && globalCuts.sideA) 
    || (m_caloCell->eta()<0 && globalCuts.sideC);
  passed = passed && globalCuts.allowedEta.contains(m_caloCell->eta());
  if(!passed) return false;
  // Check Phi Cut
  passed = false;
  for(const VP1Interval& i : globalCuts.allowedPhi){
    if(i.contains(m_caloCell->phi()) 
       ||i.contains(m_caloCell->phi()+2*M_PI)
       ||i.contains(m_caloCell->phi()-2*M_PI)) {
      passed = true;
      break;
    }
  }
  return passed;
}

energyToTransverse()

double VP1CaloCell::energyToTransverse ( const double & e ) const

inherited

Definition at line 58 of file VP1CaloCells.cxx.

{
  return m_caloCell->sinTh()*e;
}

getCaloCell()

const CaloCell * VP1CaloCell::getCaloCell ( )

inherited

Definition at line 48 of file VP1CaloCells.cxx.

{
  return m_caloCell;
}

GetFragChannel()

int VP1CC_TileBarEc::GetFragChannel ( const TileHWID * tile_hw_id, bool up, int & frag, int & channel )

virtual

Implements VP1CC_Tile.

Definition at line 574 of file VP1CaloCells.cxx.

{
  HWIdentifier hwId;
  IdContext context = tile_hw_id->channel_context();
 
  int result = 0;
 
  if(up)
    result = tile_hw_id->get_id(m_caloCell->caloDDE()->onl2(),hwId,&context);
  else
    result = tile_hw_id->get_id(m_caloCell->caloDDE()->onl1(),hwId,&context);
 
  if(result!=0)
    return result;
 
  frag = tile_hw_id->frag(hwId);
  channel = tile_hw_id->channel(hwId);
 
  return 0;
}

getID()

Identifier VP1CaloCell::getID ( )

inherited

Definition at line 53 of file VP1CaloCells.cxx.

{
  return m_caloCell->ID();
}

id2name()

std::string VP1CC_TileBarEc::id2name ( Identifier & id )

protected

Definition at line 686 of file VP1CaloCells.cxx.

{
  std::string name("");
 
  // Barrel/EC
  if(m_tileID->is_tile_barrel(id))
    name += std::string("LB");
  else
    name += std::string("EB");
 
  // Pos/Neg
  if(m_tileID->side(id)==1)
    name += std::string("A");
  else
    name += std::string("C");
 
  // Module number
  std::ostringstream modnum;
  modnum << m_tileID->module(id) + 1;
  name += (modnum.str() + std::string(" "));
 
  // Sample
  int sample = m_tileID->sample(id);
  switch(sample)
  {
  case 0:
    {
      name += std::string("A");
      break;
    }
  case 1:
    {
      if(m_tileID->is_tile_barrel(id))
      {
    name += std::string("B");
    if(m_tileID->tower(id)<9)
      name += std::string("C");
      }
      else
      {
    if(m_tileID->section(id)==2)
      name += std::string("B");
    else
      name += std::string("C");
      }
      break;
    }
  case 2:
    {
      name += std::string("D");
      break;
    }
  default:
    break;
  }
 
  // Tower
  std::ostringstream townum;
  if(m_tileID->sample(id)==2)
    townum << m_tileID->tower(id)/2;
  else
    townum << m_tileID->tower(id)+1;
 
  name += townum.str();
  return name;
}

isInsideClipVolume()

bool VP1CC_TileBarEc::isInsideClipVolume ( const VP1CC_GlobalCuts & globalCuts )

virtual

Reimplemented from VP1CaloCell.

Definition at line 656 of file VP1CaloCells.cxx.

                                                                           {
  const CaloDetDescrElement* ddElement = m_caloCell->caloDDE();
  double z = ddElement->z();
  double radius = ddElement->r();
  double dradius = ddElement->dr()*0.5;
  double minR = radius - dradius;
  
  double radialDistance = sqrt(minR*minR + z*z);
  
  if ( radialDistance > globalCuts.clipRadius ) return false;
  return true;
}

remove3DObjects()

void VP1CC_TileBarEc::remove3DObjects ( VP1CC_SoNode2CCMap * node2cc )

virtual

Implements VP1CaloCell.

Definition at line 670 of file VP1CaloCells.cxx.

{
  if(m_hitUp) {
    // Remove entries from node2cc map
    node2cc->erase(m_hitUp);
    node2cc->erase(m_hitDown);
 
    // Drop the hit from the scene and delete hit object
    m_helperUp->removeNode(m_hitUp);
    m_helperDown->removeNode(m_hitDown);
 
    m_hitUp = 0;
    m_hitDown = 0;
  }
}

ToString()

std::vector< std::string > VP1CC_TileBarEc::ToString ( const CaloCell_ID * calo_id, const std::string & extrainfos = "" )

virtual

Implements VP1CaloCell.

Definition at line 538 of file VP1CaloCells.cxx.

{
  std::vector<std::string> result;
  const TileCell* tile_cell = dynamic_cast<const TileCell*>(m_caloCell);
  if (not tile_cell) return result;
  Identifier cellid = m_caloCell->ID();
 
  double total_energy = m_caloCell->energy();
  double total_time = m_caloCell->time();
  double pmt_energy1 = tile_cell->ene1();
  double pmt_time1 = tile_cell->time1();
  double pmt_energy2 = tile_cell->ene2();
  double pmt_time2 = tile_cell->time2();
 
  std::ostringstream msg, msg1, msg2, msg3, msg4;
  msg  << "Cell: " << id2name(cellid) << " " << extrainfos;
  msg1 << "Energy = " << total_energy << " (Mev)";
  msg2 << "Time = " << total_time;
  msg3 << " PMT1 Energy = " << pmt_energy1 << " (Mev)"
       << " PMT1 Time = " << pmt_time1;
  msg4 << " PMT2 Energy = " << pmt_energy2 << " (Mev)"
       << " PMT2 Time = " << pmt_time2;
 
  std::string stars("***");
  result.push_back(stars);
 
  result.push_back(msg.str());
  result.push_back(msg1.str());
  result.push_back(msg2.str());
  result.push_back(msg3.str());
  result.push_back(msg4.str());
 
  result.push_back(stars);
  return result;
}

updateScene()

void VP1CaloCell::updateScene ( VP1CC_SoNode2CCMap * node2cc, bool useEt, const QPair< bool, double > & scale, bool outline, const VP1CC_GlobalCuts & globalCuts )

inherited

Definition at line 99 of file VP1CaloCells.cxx.

{
  if(cutPassed(globalCuts))
    build3DObjects(node2cc,useEt,scale,outline, globalCuts);
  else
    remove3DObjects(node2cc);
}

Member Data Documentation

m_caloCell

const CaloCell* VP1CaloCell::m_caloCell

protectedinherited

Definition at line 204 of file VP1CaloCells.h.

m_helperDown

VP1ExtraSepLayerHelper* VP1CC_TileBarEc::m_helperDown

protected

Definition at line 342 of file VP1CaloCells.h.

m_helperUp

VP1ExtraSepLayerHelper* VP1CC_TileBarEc::m_helperUp

protected

Definition at line 341 of file VP1CaloCells.h.

m_hitDown

SoGenericBox* VP1CC_TileBarEc::m_hitDown

protected

Definition at line 340 of file VP1CaloCells.h.

m_hitUp

SoGenericBox* VP1CC_TileBarEc::m_hitUp

protected

Definition at line 339 of file VP1CaloCells.h.

m_tileID

const TileID* VP1CC_Tile::m_tileID

protectedinherited

Definition at line 307 of file VP1CaloCells.h.

---

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

• VP1CaloCells.h
• VP1CaloCells.cxx