D3PD::SCFillerTool Class Reference

Block filler tool for EM samplings from a CaloCluster. More...

#include <SCFillerTool.h>

Inheritance diagram for D3PD::SCFillerTool:

Public Member Functions
	SCFillerTool (const std::string &type, const std::string &name, const IInterface *parent)
	Standard Gaudi tool constructor.
virtual StatusCode	book ()
	Book variables for this block.
virtual StatusCode	initialize ()
virtual StatusCode	fill (const CaloCellContainer &p)
	Fill one block — type-safe version.
virtual void	fillHashTables (const CaloCellContainer &p)
virtual void	dumpHashTables (const CaloCellContainer &p)
virtual void	resetEnergies ()
	BlockFillerToolImpl (const std::string &type, const std::string &name, const IInterface *parent)
	Standard Gaudi tool constructor.
virtual StatusCode	configureD3PD (IAddVariable *tree, const std::type_info &ti)
	Configure during initialization: type-check.
StatusCode	configureD3PD (IAddVariable *tree)
	Configure the parent tree.
virtual StatusCode	fillUntyped (const void *p, bool again=false)
	Fill one block.
virtual StatusCode	fillAgain (const CaloCellContainer &p)
	Fill one block, after AGAIN has been returned (type-safe).
virtual StatusCode	addVariable (const std::string &name, const std::type_info &ti, void *&ptr, const std::string &docstring="", const void *defval=0)
	Add a variable to the tuple.
StatusCode	addVariable (const std::string &name, CaloCellContainer *&ptr, const std::string &docstring="")
	Make the template implementation from IAddVariable visible.
StatusCode	addVariable (const std::string &name, CaloCellContainer *&ptr, const std::string &docstring, const U &defval)
	Make the template implementation from IAddVariable visible.
virtual StatusCode	addDimensionedVariable (const std::string &name, const std::type_info &ti, void *&ptr, const std::string &dim, const std::string &docstring="", const void *defval=0)
	Add a variable to the tuple.
StatusCode	addDimensionedVariable (const std::string &name, CaloCellContainer *&ptr, const std::string &dim, const std::string &docstring="")
	Add a variable to the tuple.
StatusCode	addDimensionedVariable (const std::string &name, CaloCellContainer *&ptr, const std::string &dim, const std::string &docstring, const U &defval)
	Add a variable to the tuple.

Protected Member Functions
StatusCode	configureImpl (IAddVariable *tree, const std::type_info &ti, const std::type_info &fill_ti)
	Configure during initialization: type-check.
StatusCode	configureImpl (IAddVariable *tree, const std::type_info &ti, const std::vector< const std::type_info * > &tis, size_t &which)
	Configure during initialization: type-check.
StatusCode	convert (void const *&p) const
	Do type conversion.

Protected Attributes
ObjectMetadata	m_metadata
	Metadata about the variables created by this object.

Private Attributes
const LArEM_ID *	m_emid
const LArFCAL_ID *	m_fcalid
const LArHEC_ID *	m_hecid
const TileID *	m_tileid
const LArOnlineID *	m_onlineHelper
SG::ReadCondHandleKey< LArOnOffIdMapping >	m_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}
int *	m_nSC
	parameters
std::vector< double > *	m_E
std::vector< double > *	m_eta
std::vector< double > *	m_phi
std::vector< double > *	m_Et
std::vector< double > *	m_time
std::vector< double > *	m_quality
std::vector< int > *	m_calo
std::vector< int > *	m_region
std::vector< int > *	m_sampling
std::vector< int > *	m_ieta
std::vector< int > *	m_jphi
std::vector< int > *	m_hashid
bool	m_tableFilled
std::map< const Identifier, int >	m_sc_hashTable
std::vector< int >	m_sc_list
std::map< int, double >	m_sc_energy
std::map< int, double >	m_sc_eta
std::map< int, double >	m_sc_phi
std::map< int, std::vector< const CaloCell * > >	m_sc_compId
std::vector< float >	m_etaCut
std::vector< float >	m_phiCut
std::vector< unsigned int >	m_caloNums
std::vector< unsigned int >	m_caloLayers
bool	m_caloEtaSelection
bool	m_caloPhiSelection
bool	m_caloLayerSelection
bool	m_caloSelection
bool	m_tileDLayerOption
std::string	m_prefix
	Property: the variable prefix for this block.
std::string	m_suffix
	Property: the variable prefix for this block.
std::string	m_blockName
	Property: the name of this block.
TypeConverter	m_converter
	Helper to do type conversions.
IAddVariable *	m_tree
	Reference to the block name.

Static Private Attributes
static const std::string	s_emptyString
	Used to provide a default value for blockName.

Detailed Description

Block filler tool for EM samplings from a CaloCluster.

Definition at line 44 of file SCFillerTool.h.

Constructor & Destructor Documentation

SCFillerTool()

D3PD::SCFillerTool::SCFillerTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Standard Gaudi tool constructor.

Parameters

type	The name of the tool type.
name	The tool name.
parent	The tool's Gaudi parent.

Definition at line 31 of file SCFillerTool.cxx.

      : BlockFillerTool<CaloCellContainer> (type, name, parent), 
    m_tableFilled(false),
    m_caloEtaSelection( false ),
    m_caloPhiSelection( false ),
    m_caloLayerSelection( false ),
    m_caloSelection( false ), 
    m_tileDLayerOption( false ) {
 
  declareProperty( "CaloEtaCut", m_etaCut );
  declareProperty( "CaloPhiCut", m_phiCut );
  declareProperty( "CaloLayers", m_caloLayers );
  // 1=EMB; 2=EMEC Outer Wheel; 3=EMEC Inner Wheel; 4=HEC; 5=FCAL; 6=Tile LB; 7=Tile EB
  declareProperty( "CaloDetectors", m_caloNums ); 
  declareProperty( "TileDLayerOption", m_tileDLayerOption);
 
  SCFillerTool::book().ignore(); // Avoid coverity warnings.
}

Member Function Documentation

addDimensionedVariable() [1/3]

StatusCode D3PD::IAddVariable::addDimensionedVariable ( const std::string & name, CaloCellContainer *& ptr, const std::string & dim, const std::string & docstring, const U & defval )

inherited

Add a variable to the tuple.

Parameters

name	The name of the variable.
ptr	Pointer to the type of the variable. The pointer need not be initialized; the D3PD software will set the pointer prior to calling fill().
dim	Dimension for the variable. (Presently unimplemented!)
docstring	Documentation string for this variable.
defval	Pointer to the default value to use for this variable. Only works for basic types.

If called from the constructor, the only effect is to clear ptr.

addDimensionedVariable() [2/3]

StatusCode D3PD::IAddVariable::addDimensionedVariable ( const std::string & name, CaloCellContainer *& ptr, const std::string & dim, const std::string & docstring = "" )

inherited

Add a variable to the tuple.

Parameters

name	The name of the variable.
ptr	Pointer to the type of the variable. The pointer need not be initialized; the D3PD software will set the pointer prior to calling fill().
dim	Dimension for the variable. (Presently unimplemented!)
docstring	Documentation string for this variable.

If called from the constructor, the only effect is to clear ptr.

addDimensionedVariable() [3/3]

StatusCode D3PD::AddVariable::addDimensionedVariable ( const std::string & name, const std::type_info & ti, void *& ptr, const std::string & dim, const std::string & docstring = "", const void * defval = 0 )

virtualinherited

Add a variable to the tuple.

Parameters

name	The name of the variable.
type	The type of the variable.
ptr	Pointer to the type of the variable. The pointer need not be initialized; the D3PD software will set the pointer prior to calling fill().
dim	Dimension for the variable. (Presently unimplemented!)
docstring	Documentation string for this variable.
defval	Pointer to the default value to use for this variable. Null for no default (generally means to fill with zeros). Of the type given by ti. Only works for basic types.

If called from the constructor, the only effect is to clear ptr.

Definition at line 119 of file AddVariable.cxx.

{
  if (!m_tree) {
    // Called from constructor --- just initialize pointer.
    ptr = 0;
    return StatusCode::SUCCESS;
  }
 
  CHECK( m_metadata.addDimensionedVariable( m_prefix + name + m_suffix,
                                            ti, ptr, dim,
                                            docstring, defval ) );
  return m_tree->addDimensionedVariable (m_prefix + name + m_suffix,
                                         ti, ptr, dim,
                                         docstring, defval);
}

addVariable() [1/3]

StatusCode D3PD::IAddVariable::addVariable ( const std::string & name, CaloCellContainer *& ptr, const std::string & docstring, const U & defval )

inherited

Make the template implementation from IAddVariable visible.

addVariable() [2/3]

StatusCode D3PD::IAddVariable::addVariable ( const std::string & name, CaloCellContainer *& ptr, const std::string & docstring = "" )

inherited

Make the template implementation from IAddVariable visible.

addVariable() [3/3]

StatusCode D3PD::AddVariable::addVariable ( const std::string & name, const std::type_info & ti, void *& ptr, const std::string & docstring = "", const void * defval = 0 )

virtualinherited

Add a variable to the tuple.

Parameters

name	The name of the variable.
type	The type of the variable.
ptr	Pointer to the type of the variable. The pointer need not be initialized; the D3PD software will set the pointer prior to calling fill().
docstring	Documentation string for this variable.
defval	Pointer to the default value to use for this variable. Null for no default (generally means to fill with zeros). Of the type given by ti. Only works for basic types.

If called from the constructor, the only effect is to clear ptr.

Definition at line 93 of file AddVariable.cxx.

{
  if (!m_tree) {
    // Called from constructor --- just initialize pointer.
    ptr = 0;
    return StatusCode::SUCCESS;
  }
 
  CHECK( m_metadata.addVariable( m_prefix + name + m_suffix,
                                 ti, ptr, docstring, defval ) );
  return m_tree->addVariable (m_prefix + name + m_suffix,
                              ti, ptr, docstring, defval);
}

BlockFillerToolImpl()

( const std::string & type, const std::string & name, const IInterface * parent )

inherited

Standard Gaudi tool constructor.

Parameters

type	The name of the tool type.
name	The tool name.
parent	The tool's Gaudi parent.

book()

StatusCode D3PD::SCFillerTool::book ( )

virtual

Book variables for this block.

Implements D3PD::BlockFillerTool< CaloCellContainer >.

Definition at line 138 of file SCFillerTool.cxx.

{
  CHECK( addVariable ("nSC", m_nSC ));
  CHECK( addVariable ("E", m_E));
  CHECK( addVariable ("eta", m_eta));
  CHECK( addVariable ("phi", m_phi));
  CHECK( addVariable ("Et", m_Et));
  CHECK( addVariable ("calo", m_calo));
  CHECK( addVariable ("region", m_region));
  CHECK( addVariable ("sampling", m_sampling));
  CHECK( addVariable ("ieta", m_ieta));
  CHECK( addVariable ("jphi", m_jphi));
  CHECK( addVariable ("hid", m_hashid));
  CHECK( addVariable ("time", m_time));
  CHECK( addVariable ("quality", m_quality));
 
  return StatusCode::SUCCESS;
}

configureD3PD() [1/2]

StatusCode D3PD::AddVariable::configureD3PD ( IAddVariable * tree )

inherited

Configure the parent tree.

Parameters

tree	The parent IAddVariable instance.

Definition at line 69 of file AddVariable.cxx.

{
  m_tree = tree;
  m_metadata.clear();
  return StatusCode::SUCCESS;
}

configureD3PD() [2/2]

virtual StatusCode D3PD::BlockFillerTool< CaloCellContainer >::configureD3PD ( IAddVariable * tree, const std::type_info & ti )

virtualinherited

Configure during initialization: type-check.

Parameters

tree	Our parent for tuple making.
ti	Gives the type of the object being passed to fillUntyped.

configureD3PD should check that the type of the object coming as input is compatible with what it expects, and raise an error otherwise.

configureImpl() [1/2]

StatusCode D3PD::BlockFillerToolImpl::configureImpl ( IAddVariable * tree, const std::type_info & ti, const std::type_info & fill_ti )

protectedinherited

Configure during initialization: type-check.

Parameters

tree	Our parent for tuple making.
ti	Gives the type of the object being passed to fillUntyped.
fill_ti	is the type of object that the user code will expect.

This is the common implementation for the configureD3PD method of IBlockFillerTool. It checks that the type of the object coming as input is compatible with what the user code is expecting, and raises an error otherwise.

Definition at line 70 of file BlockFillerToolImpl.cxx.

{
  CHECK( AddVariable::configureD3PD (tree) );
  CHECK( m_converter.init (ti, fill_ti) );
 
  return StatusCode::SUCCESS;
}

configureImpl() [2/2]

StatusCode D3PD::BlockFillerToolImpl::configureImpl ( IAddVariable * tree, const std::type_info & ti, const std::vector< const std::type_info * > & tis, size_t & which )

protectedinherited

Configure during initialization: type-check.

Parameters

	tree	Our parent for tuple making.
	ti	Gives the type of the object being passed to fillUntyped.
	tis	List of possible input types that we can accept.
[out]	which	Index of the accepted type.

This is the common implementation for the configureD3PD method of IBlockFillerTool. It checks that the type of the object coming as input is compatible with what the user code is expecting, and raises an error otherwise.

This version allows for one of a set of types to match. If successful, WHICH is set to the index (0-based) of the matching type.

Definition at line 91 of file BlockFillerToolImpl.cxx.

{
  CHECK( AddVariable::configureD3PD (tree) );
  CHECK( m_converter.init (ti, tis, which) );
  return StatusCode::SUCCESS;
}

convert()

StatusCode D3PD::BlockFillerToolImpl::convert ( void const *& p ) const

protectedinherited

Do type conversion.

Parameters

p[in,out]

Pointer to convert.

On entry, p is a pointer to an object of the type that was passed into configureImpl as ti. Convert the pointer to point to an object of the type that was given as fill_ti. If the conversion fails, the pointer is set to 0 (and FAILURE is returned).

Definition at line 107 of file BlockFillerToolImpl.cxx.

{
  p = m_converter.convertUntyped (p);
  if (!p) {
    REPORT_MESSAGE (MSG::ERROR)
      << "Pointer conversion from " << m_converter.srcName() << " to "
      << m_converter.dstName() << "failed.";
    return StatusCode::FAILURE;
  }
  return StatusCode::SUCCESS;
}

dumpHashTables()

void D3PD::SCFillerTool::dumpHashTables ( const CaloCellContainer & p )

virtual

Definition at line 248 of file SCFillerTool.cxx.

                                                              {
//---------------------------------------------------------------------------------
  typedef CaloCellContainer CONTAINER;
  CONTAINER::const_iterator f_cell = p.begin();
  CONTAINER::const_iterator l_cell = p.end();
 
  for( ; f_cell != l_cell; ++f_cell ) {
      const CaloCell* cell = ( *f_cell ); 
      Identifier id = cell->ID();
      int posneg = (cell->eta()>=0.0) ? 1: -1;
      int sc_hash = m_sc_hashTable[ id ];
      //      m_sc_compId[sc_hash];
      if ( posneg > 0 ) {
    if ( m_emid->is_lar_em(id) ) {
      if (m_emid->is_em_barrel(id) ) {
        int ieta = m_emid->eta(id);
        int jphi = m_emid->phi(id);
        int lay = m_emid->sampling(id);
        // int reg = m_emid->region(id);
        std::cout << "SCFillerTool::dumpHashTables() [ONLY EM BARREL] ==> posneg,lay,ieta,jphi= " << posneg << ", " << lay << ", " << ieta << ", " << jphi << " ID: " << id << " sc_hash: " << sc_hash << std::endl;
      }
    }
      }
  }
  std::map< int, std::vector<const CaloCell*> >::iterator it = m_sc_compId.begin();
  for( ; it!=m_sc_compId.end(); ++it) {
    int scID = (*it).first;
    int key = scID;
    std::vector<const CaloCell*> v = (*it).second;
    int nl = v.size();
    int pn = ((key >> 28) & 0x1) ? 1: -1;
    int calo = ((key >> 24) & 0xf)*pn;
    int region =  (key>>20) & 0xf;
    unsigned int lay = (key>>16) & 0xf;
    int ieta = (key>>8) & 0xff;
    int jphi = key & 0xff;
    if (abs(calo)==1 && pn>0) {
      std::cout << "SCFillerTool::dumpHashTables() ===> scID = " << scID << " calo,pn,region,lay,ieta,jphi: " << calo << ", " << pn << ", " << region  << ", " << lay << ", " << ieta << ", " << jphi << std::endl;
      std::cout << "SCFillerTool::dumpHashTables() ===> scID = " << scID << " No. of cells: " << nl << std::endl;
      std::cout << "SCFillerTool::dumpHashTables() ===>        " ;
      for (const CaloCell* cell : v) {
    std::cout << cell << ", " ;
      }
      std::cout << "" << std::endl;
    }
  }
}

fill()

StatusCode D3PD::SCFillerTool::fill ( const CaloCellContainer & p )

virtual

Fill one block — type-safe version.

Parameters

p	The input object.

This is called once per object. The caller is responsible for arranging that all the pointers for booked variables are set appropriately upon entry.

Implements D3PD::BlockFillerTool< CaloCellContainer >.

Definition at line 166 of file SCFillerTool.cxx.

                                                         {
 
   if ( !m_tableFilled ) {
     m_tableFilled = true;
     fillHashTables(p);
     dumpHashTables(p);
     ATH_MSG_INFO("**** after fillHashTables in first event");
   }
   
   resetEnergies();
 
   typedef CaloCellContainer CONTAINER;
   CONTAINER::const_iterator f_cell = p.begin();
   CONTAINER::const_iterator l_cell = p.end();
   for( ; f_cell != l_cell; ++f_cell ) {
     const CaloCell* cell = ( *f_cell ); 
     const Identifier id = cell->ID();
     int scid = m_sc_hashTable[id];
     m_sc_energy[scid] += cell->energy();
   }
 
   int nSC = 0;
   std::vector<int>::iterator it;
   it=m_sc_list.begin();
   for ( ; it!=m_sc_list.end(); ++it ) {
     int key = *it;
     //      int sc_hash = (sc_side << 28) | ( sc_calo << 24 ) | (sc_region << 20 ) | ( sc_layer << 16 ) | ( sc_ieta << 8 ) | ( sc_jphi );
     int pn = ((key >> 28) & 0x1) ? 1: -1;
     int calo = ((key >> 24) & 0xf)*pn;
     int region =  (key>>20) & 0xf;
     unsigned int lay = (key>>16) & 0xf;
     int ieta = (key>>8) & 0xff;
     int jphi = key & 0xff;
 
     float feta = m_sc_eta[key];
     float fphi = m_sc_phi[key];
     int lardet = abs((key >> 24) & 0xf);
 
     if (m_caloSelection) {
       std::vector< unsigned int >::const_iterator theFound =
     std::find( m_caloNums.begin(), m_caloNums.end(), lardet );
       if( theFound == m_caloNums.end() ) continue;
     }
     if ( m_caloEtaSelection ) {
       if (feta<m_etaCut[0] || feta>m_etaCut[1]) continue;
     }
     if ( m_caloPhiSelection ) {
       if (fphi<m_phiCut[0] || fphi>m_phiCut[1]) continue;
     }
     if ( m_caloLayerSelection ) {
       bool tmp=false;
       for( unsigned int j=0; j<m_caloLayers.size(); j++) {
     if (m_caloLayers[j]==lay) tmp=true;
       }
       if (!tmp) continue;
    }
     nSC++;
     m_E->push_back( m_sc_energy[ key ]  );
     m_calo->push_back( calo );
     m_region->push_back( region );
     m_sampling->push_back( lay );
     m_ieta->push_back( ieta );
     m_jphi->push_back( jphi );
     m_hashid->push_back( key );
     m_eta->push_back( feta  );
     m_phi->push_back( m_sc_phi[ key ] );
     m_Et->push_back( m_sc_energy[key]/cosh( feta ) );
   }
  *m_nSC = nSC;
 
  return StatusCode::SUCCESS;
}

fillAgain()

virtual StatusCode D3PD::BlockFillerTool< CaloCellContainer >::fillAgain ( const CaloCellContainer & p )

virtualinherited

Fill one block, after AGAIN has been returned (type-safe).

Parameters

p	The input object.

Once fill returns AGAIN, the parent should call fillAgain with the same arguments. This continues until fillAgain returns something other than AGAIN.

By default, this just calls fill().

fillHashTables()

void D3PD::SCFillerTool::fillHashTables ( const CaloCellContainer & p )

virtual

Definition at line 297 of file SCFillerTool.cxx.

                                                              {
//---------------------------------------------------------------------------------
  typedef CaloCellContainer CONTAINER;
  CONTAINER::const_iterator f_cell = p.begin();
  CONTAINER::const_iterator l_cell = p.end();
 
  int /*posneg,*/reg,lay,ieta,jphi;
  int sc_side, sc_calo, sc_region, sc_layer, sc_ieta, sc_jphi;
  //float sc_feta, sc_fphi;
 
  SG::ReadCondHandle<LArOnOffIdMapping> cabling (m_cablingKey);
 
  for( ; f_cell != l_cell; ++f_cell ) {
      const CaloCell* cell = ( *f_cell ); 
      Identifier id = cell->ID();
      //posneg = (cell->eta()>=0.0) ? 1: -1;
      sc_side = (cell->eta()>=0.0) ? 1: 0;
      sc_ieta = -999;
      sc_jphi = -999;
      //sc_feta = 0.;
      //sc_fphi = 0.;
      sc_calo = -1;
      sc_region = -1;
      sc_layer = -1;
 
      ieta = -999;
      jphi = -999;
 
      if ( m_emid->is_lar_em(id) ) {
    ieta = m_emid->eta(id);
    jphi = m_emid->phi(id);
    lay = m_emid->sampling(id);
    reg = m_emid->region(id);
    sc_region = reg;
    sc_layer = lay;
    if (m_emid->is_em_barrel(id) ) {
      sc_calo = 1;
      if ( lay == 0 ) {
        sc_ieta = int(ieta/4);
        sc_jphi = jphi;
      } else if ( lay == 1 ) {
        sc_ieta = (reg==0) ? int(ieta/8) : 56+ieta;
        sc_jphi = (reg==0) ? jphi: int(jphi/4);
      } else if ( lay == 2 ) { 
        sc_ieta = (reg==0) ? ieta: 56 + ieta;
        sc_jphi = int(jphi/4);
      } else if ( lay == 3 ) { 
        sc_ieta = int(ieta/2); 
        sc_jphi = int(jphi/4);
      }
    } else if ( m_emid->is_em_endcap_outer(id) ) {
      sc_calo = 2;
      if ( lay == 0 ) {
        sc_ieta = int(ieta/4);
        sc_jphi = jphi;
        sc_region = 1;
      } else if (lay == 1) {
        //      sc_ieta = (reg<=1) ? 0 : ((reg==2) ? 1+int(ieta/8) : ((reg==3) ? 13+int(ieta/8) : ((reg==4) ? 19+int(ieta/4): 35)));
      sc_ieta = (reg<=1) ? 0 : ((reg==2) ? 1+int(ieta/8) : ((reg==3) ? 13+int(ieta/4) : ((reg==4) ? 25+int(ieta/4): 41)));
        sc_jphi = jphi;
        sc_region = (reg<2) ? 0: reg-1;
      } else if (lay==2) {
        sc_ieta = (reg==0) ? 0 : 1+ieta;  
        sc_jphi = int(jphi/4);
        sc_region = (ieta<3) ? 0: ((ieta<15) ? 1 : ((ieta<23) ? 2 : ((ieta<39) ? 3 : 4)));
      } else if (lay==3) {
        sc_ieta = int(ieta/2); 
        sc_jphi = int(jphi/4);
        sc_region = (ieta<6) ? 1: ((ieta<10) ? 2 : ((ieta<18) ? 3: 4));
      } else {
        sc_ieta = -999;
        sc_jphi = -999;
        sc_region = -999;
      }
    } else if ( m_emid->is_em_endcap_inner(id) ) {
      sc_calo = 3;
      sc_ieta =  (ieta<6) ? int(ieta/2) : 47;
      sc_jphi =  int(jphi/2);
      sc_region = (ieta<6) ? 5 : 6;
    } 
      } else if (m_hecid->is_lar_hec(id)) {
    sc_layer  = m_hecid->sampling(id);
    sc_region = m_hecid->region(id);
    sc_calo= 4;
    ieta = m_hecid->eta(id);
    jphi = m_hecid->phi(id);
    sc_ieta = (sc_region==0) ? ieta : 10+ieta; 
    sc_jphi = jphi;
    sc_region = 0;
      } else if (m_fcalid->is_lar_fcal(id)) {
    HWIdentifier hwid = cabling->createSignalChannelID( id );
    sc_calo = 5;
    const int lay = m_fcalid->module(id);
    const int s = m_onlineHelper->slot(hwid);
    const int c  = m_onlineHelper->channel(hwid);
    const int c64 = c%64;
    sc_layer = lay;
    if (lay==3) {
      sc_ieta = 4-int((c%16)/4);
      int jphip = 8*(s-14) + 3 + 4*int(c/64)-int(c64/16);
      int jphin = 8*(s-14) + 7 - int(c/16);
      sc_jphi = (sc_side)>0 ? jphip : jphin;
      sc_region = 0;
    } else if (lay==2) {
      sc_ieta = ((c%32)<16 ? 4 : 8) - int((c%32)/4);
      int jphip = 4*(s-10) + 1 + 2*int(c/64)-int(c64/32);
      int jphin =  ((s>=12) ? 14 -4*(s-12): 6 -4*(s-10) ) + int(c64/32) -2*int(c/64);
      sc_jphi = (sc_side>0) ? jphip : jphin;
      sc_region = 0;
    } else if (lay==1) {
      ieta  = 15-int(c64%16) + 16*int(c64/16);
      sc_ieta = int(ieta/4);
      int jphip = int(c/64) + 2*(s- ((s==9) ? 2 : 1));
      int jphin = ((s<5) ? 7 - 2*(s-1) : 15 - 2*(s-((s==9)?6:5))) - int(c/64);
      sc_jphi = (sc_side>0) ? jphip : jphin;
      if (c64<16) {
        sc_region = 0;
      } else if (c64>=16 && c64<32 ) {
        sc_region = 1;
      } else if (c64>=32 && c64<48 ) {
        sc_region = 2;
      } else {
        sc_region = 3;
      }
    }
      } else if (m_tileid->is_tile(id) ) {
    sc_region = 0;
    sc_layer = m_tileid->sampling(id);
    sc_side = (m_tileid->side(id)>0) ? 1 : 0;
    ieta = m_tileid->tower(id);
    jphi = m_tileid->module(id);
    sc_ieta = ieta;
    sc_jphi = jphi;
    if ( m_tileid->is_tile_barrel(id) )
      sc_calo = 6;
    if ( m_tileid->is_tile_extbarrel(id) )
      sc_calo = 7;
      }
      if( sc_ieta<0 || sc_jphi<0 || sc_calo<=0 || sc_region<0 || sc_layer<0 ) continue;
      // This is a kludge so that for HEC/Tile (sc_calo=4,6,7) the SC hash id don't have any layer information...(11/27/12)
      if (sc_calo==4)
    sc_layer = 0;
      if (sc_calo==6 || sc_calo==7) {
    if ( sc_layer!=2 || !m_tileDLayerOption)  
      sc_layer = 0; 
      }
      int sc_hash = (sc_side << 28) | ( sc_calo << 24 ) | (sc_region << 20 ) | ( sc_layer << 16 ) | ( sc_ieta << 8 ) | ( sc_jphi );
      m_sc_hashTable[ id ] =  sc_hash;
      m_sc_compId[sc_hash].push_back( cell );
  }
 
  std::map< int, std::vector<const CaloCell*> >::iterator it = m_sc_compId.begin();
  for( ; it!=m_sc_compId.end(); ++it) {
    int scID = (*it).first;
    m_sc_list.push_back( scID );
    std::vector<const CaloCell*> v = (*it).second;
    int nl = v.size();
    std::vector<const CaloCell*>::iterator vit;
    /*
    if ( ((scID >> 24) & 0xf) >= 6 ) {
      std::cout << "*** KAKKA TILECAL: THIS IS A TILECAL SCELL hash= " << std::hex << scID << " no. of cells: " << nl << std::endl;
      for( vit=v.begin(); vit!=v.end(); ++vit) {
    std::cout  << " " << (*vit)->ID() << " ";
      }
      std::cout  << "" << std::endl;
      for( vit=v.begin(); vit!=v.end(); ++vit) {
    std::cout  << " " << (*vit)->eta() << " ";
      }
      std::cout  << "" << std::endl;
      for( vit=v.begin(); vit!=v.end(); ++vit) {
    std::cout  << " " << (*vit)->phi() << " ";
      }
      std::cout  << "" << std::endl;
    }
    */
    double feta = 0;
    double fphi = 0;
    for( vit=v.begin(); vit!=v.end(); ++vit) {
      feta += (*vit)->eta();
      fphi += (*vit)->phi();
    }
    if (nl>0) {
      feta /= nl;
      fphi /= nl;
    }
    m_sc_eta[ scID ] = feta;
    m_sc_phi[ scID ] = fphi;
  } // for (; it!=m_sc_compId...)
 
 
}

fillUntyped()

virtual StatusCode D3PD::BlockFillerTool< CaloCellContainer >::fillUntyped ( const void * p, bool again = false )

virtualinherited

Fill one block.

Parameters

p	The input object.
again	Set if this is a subsequent call requested by an AGAIN return

This is called once per object. The type of the object at which p points is given by the ti argument to configureD3PD. The caller is responsible for arranging that all the pointers for booked variables are set appropriately upon entry.

If the return status is the special code AGAIN (defined above), then this filler tool wants to make multiple entries. The parent should set up to capture a new ‘row’ and run through the list of block filler tools again, but for this tool call fillAgainUntyped instead of fillUntyped. This should be repeated as long as fillAgainUntyped returns AGAIN.

Once fillUntyped returns AGAIN, the parent should call fillUntyped with the same p argument and again set to true. This continues until fillUntyped returns something other than AGAIN.

Not all parents may support this. In that case, returning AGAIN will be treated as an error.

initialize()

StatusCode D3PD::SCFillerTool::initialize ( )

virtual

Definition at line 53 of file SCFillerTool.cxx.

                                    {
 
  ServiceHandle<StoreGateSvc> detStore("DetectorStore", name());
  CHECK( detStore.retrieve() ) ;
    
  //  CHECK ( detStore->retrieve(m_dd_man)   );
  CHECK ( detStore->retrieve(m_emid)   );
  CHECK ( detStore->retrieve(m_fcalid) );
  CHECK ( detStore->retrieve(m_hecid)  );
  CHECK ( detStore->retrieve(m_tileid) );
  CHECK ( detStore->retrieve(m_onlineHelper));
 
  const unsigned int nSubCalo = 7; // static_cast< int >( CaloCell_ID::NSUBCALO );
 
  //check calo number specified
  m_caloSelection = false;
  if( m_caloNums.size() == 0 ) {
    ATH_MSG_INFO( " *** SCFillerTool: No calorimeter selection" );
    return StatusCode::SUCCESS;
  } else if( m_caloNums.size() > nSubCalo ) {
    REPORT_MESSAGE( MSG::FATAL )
      << " More than " << nSubCalo << " calo specified. "
      << "Must be wrong. Stop.";
    return StatusCode::FAILURE;
  } else {
    m_caloSelection = true;
    for( unsigned int index = 0; index < m_caloNums.size(); ++index ) {
      if( m_caloNums[ index ]> nSubCalo ) {
    REPORT_MESSAGE( MSG::FATAL )
      << "Invalid calo specification:" 
      << m_caloNums[index] << "Stop.";
    return StatusCode::FAILURE;
      } else {
    ATH_MSG_INFO( "*** SCFillerTool:  Select calorimeter " << m_caloNums[ index ] );
      }
    }
  }
 
  m_caloEtaSelection = false;
  if ( m_etaCut.size() >0 ) {
    if ( m_etaCut.size()== 2 && m_etaCut[0]<m_etaCut[1] ) {
      REPORT_MESSAGE( MSG::INFO ) <<  "*** SCFillerTool: Eta Selection specified: " << m_etaCut[0] << ", " << m_etaCut[1];
      m_caloEtaSelection=true;
    } else {
      REPORT_MESSAGE( MSG::FATAL ) << " *** SCFillerTool: Invalid Eta Cut Range specified: Abort!";
      return StatusCode::FAILURE;
    }
  }
 
  m_caloPhiSelection = false;
  if ( m_phiCut.size() >0 ) {
    if ( m_phiCut.size() == 2 && m_phiCut[0]<m_phiCut[1] ) {
      REPORT_MESSAGE( MSG::INFO ) << "*** SCFillerTool: Phi Selection specified: "  << m_phiCut[0] << ", " << m_phiCut[1];
      m_caloPhiSelection = true;
    } else {
      REPORT_MESSAGE( MSG::FATAL ) << " *** SCFillerTool: Invalid Phi Cut Range specified: Abort!";
      return StatusCode::FAILURE;
    }
  }
 
  m_caloLayerSelection = false;
  if ( m_caloLayers.size()>0 ) {
    if ( m_caloLayers.size()>28 ) {
      REPORT_MESSAGE( MSG::FATAL ) << " *** SCFillerTool: Invalid No. of Calo Layer selection: " << m_caloLayers.size() << " Abort!";
      return StatusCode::FAILURE;
    }
    for( unsigned int j=0; j<m_caloLayers.size(); j++) {
      if (m_caloLayers[j]>27  ) {
    REPORT_MESSAGE( MSG::FATAL ) << " *** SCFillerTool: Invalid Calo Layer selection: " << m_caloLayers[j] << " Abort!";
    return StatusCode::FAILURE;
      }
    }
    m_caloLayerSelection = true;
  }
 
  ATH_CHECK( m_cablingKey.initialize() );
 
  ATH_MSG_INFO( " *** SCFillerTool: completed" );
  return StatusCode::SUCCESS;
}

resetEnergies()

void D3PD::SCFillerTool::resetEnergies ( )

virtual

Definition at line 239 of file SCFillerTool.cxx.

                                  {
  std::map<int,double>::iterator it;
  it=m_sc_energy.begin();
  for ( ; it!=m_sc_energy.end(); ++it ) {
    (*it).second=0.0;
  }
}

Member Data Documentation

m_blockName

std::string D3PD::BlockFillerToolImpl::m_blockName

privateinherited

Property: the name of this block.

Definition at line 118 of file BlockFillerToolImpl.h.

m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> D3PD::SCFillerTool::m_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}

private

Definition at line 74 of file SCFillerTool.h.

{this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"};

m_calo

std::vector<int>* D3PD::SCFillerTool::m_calo

private

Definition at line 88 of file SCFillerTool.h.

m_caloEtaSelection

bool D3PD::SCFillerTool::m_caloEtaSelection

private

Definition at line 108 of file SCFillerTool.h.

m_caloLayers

std::vector< unsigned int > D3PD::SCFillerTool::m_caloLayers

private

Definition at line 106 of file SCFillerTool.h.

m_caloLayerSelection

bool D3PD::SCFillerTool::m_caloLayerSelection

private

Definition at line 110 of file SCFillerTool.h.

m_caloNums

std::vector< unsigned int > D3PD::SCFillerTool::m_caloNums

private

Definition at line 105 of file SCFillerTool.h.

m_caloPhiSelection

bool D3PD::SCFillerTool::m_caloPhiSelection

private

Definition at line 109 of file SCFillerTool.h.

m_caloSelection

bool D3PD::SCFillerTool::m_caloSelection

private

Definition at line 111 of file SCFillerTool.h.

m_converter

TypeConverter D3PD::BlockFillerToolImpl::m_converter

privateinherited

Helper to do type conversions.

Definition at line 121 of file BlockFillerToolImpl.h.

m_E

std::vector<double>* D3PD::SCFillerTool::m_E

private

Definition at line 81 of file SCFillerTool.h.

m_emid

const LArEM_ID* D3PD::SCFillerTool::m_emid

private

Definition at line 69 of file SCFillerTool.h.

m_Et

std::vector<double>* D3PD::SCFillerTool::m_Et

private

Definition at line 84 of file SCFillerTool.h.

m_eta

std::vector<double>* D3PD::SCFillerTool::m_eta

private

Definition at line 82 of file SCFillerTool.h.

m_etaCut

std::vector<float> D3PD::SCFillerTool::m_etaCut

private

Definition at line 103 of file SCFillerTool.h.

m_fcalid

const LArFCAL_ID* D3PD::SCFillerTool::m_fcalid

private

Definition at line 70 of file SCFillerTool.h.

m_hashid

std::vector<int>* D3PD::SCFillerTool::m_hashid

private

Definition at line 93 of file SCFillerTool.h.

m_hecid

const LArHEC_ID* D3PD::SCFillerTool::m_hecid

private

Definition at line 71 of file SCFillerTool.h.

m_ieta

std::vector<int>* D3PD::SCFillerTool::m_ieta

private

Definition at line 91 of file SCFillerTool.h.

m_jphi

std::vector<int>* D3PD::SCFillerTool::m_jphi

private

Definition at line 92 of file SCFillerTool.h.

m_metadata

ObjectMetadata D3PD::AddVariable::m_metadata

protectedinherited

Metadata about the variables created by this object.

Definition at line 129 of file AddVariable.h.

m_nSC

int* D3PD::SCFillerTool::m_nSC

private

parameters

Definition at line 80 of file SCFillerTool.h.

m_onlineHelper

const LArOnlineID* D3PD::SCFillerTool::m_onlineHelper

private

Definition at line 73 of file SCFillerTool.h.

m_phi

std::vector<double>* D3PD::SCFillerTool::m_phi

private

Definition at line 83 of file SCFillerTool.h.

m_phiCut

std::vector<float> D3PD::SCFillerTool::m_phiCut

private

Definition at line 104 of file SCFillerTool.h.

m_prefix

std::string D3PD::BlockFillerToolImpl::m_prefix

privateinherited

Property: the variable prefix for this block.

Definition at line 112 of file BlockFillerToolImpl.h.

m_quality

std::vector<double>* D3PD::SCFillerTool::m_quality

private

Definition at line 86 of file SCFillerTool.h.

m_region

std::vector<int>* D3PD::SCFillerTool::m_region

private

Definition at line 89 of file SCFillerTool.h.

m_sampling

std::vector<int>* D3PD::SCFillerTool::m_sampling

private

Definition at line 90 of file SCFillerTool.h.

m_sc_compId

std::map<int, std::vector<const CaloCell*> > D3PD::SCFillerTool::m_sc_compId

private

Definition at line 101 of file SCFillerTool.h.

m_sc_energy

std::map<int, double> D3PD::SCFillerTool::m_sc_energy

private

Definition at line 98 of file SCFillerTool.h.

m_sc_eta

std::map<int, double> D3PD::SCFillerTool::m_sc_eta

private

Definition at line 99 of file SCFillerTool.h.

m_sc_hashTable

std::map<const Identifier, int> D3PD::SCFillerTool::m_sc_hashTable

private

Definition at line 96 of file SCFillerTool.h.

m_sc_list

std::vector<int> D3PD::SCFillerTool::m_sc_list

private

Definition at line 97 of file SCFillerTool.h.

m_sc_phi

std::map<int, double> D3PD::SCFillerTool::m_sc_phi

private

Definition at line 100 of file SCFillerTool.h.

m_suffix

std::string D3PD::BlockFillerToolImpl::m_suffix

privateinherited

Property: the variable prefix for this block.

Definition at line 115 of file BlockFillerToolImpl.h.

m_tableFilled

bool D3PD::SCFillerTool::m_tableFilled

private

Definition at line 95 of file SCFillerTool.h.

m_tileDLayerOption

bool D3PD::SCFillerTool::m_tileDLayerOption

private

Definition at line 113 of file SCFillerTool.h.

m_tileid

const TileID* D3PD::SCFillerTool::m_tileid

private

Definition at line 72 of file SCFillerTool.h.

m_time

std::vector<double>* D3PD::SCFillerTool::m_time

private

Definition at line 85 of file SCFillerTool.h.

m_tree

IAddVariable* D3PD::AddVariable::m_tree

privateinherited

Reference to the block name.

The parent IAddVariable instance.

Definition at line 143 of file AddVariable.h.

s_emptyString

const std::string D3PD::AddVariable::s_emptyString

staticprivateinherited

Used to provide a default value for blockName.

Definition at line 41 of file AddVariable.h.

---

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

• SCFillerTool.h
• SCFillerTool.cxx