LVL1::gSuperCellTowerMapper Class Reference

#include <gSuperCellTowerMapper.h>

Inheritance diagram for LVL1::gSuperCellTowerMapper:

Collaboration diagram for LVL1::gSuperCellTowerMapper:

Public Member Functions
	gSuperCellTowerMapper (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~gSuperCellTowerMapper ()
virtual StatusCode	initialize () override
	standard Athena-Algorithm method
virtual StatusCode	AssignSuperCellsToTowers (std::unique_ptr< gTowerContainer > &my_gTowerContainerRaw) const override
virtual StatusCode	AssignTriggerTowerMapper (std::unique_ptr< gTowerContainer > &my_gTowerContainerRaw) const override
virtual void	reset () override
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
virtual int	FindAndConnectTower (std::unique_ptr< gTowerContainer > &my_gTowerContainerRaw, CaloSampling::CaloSample sample, const int region, const int pos_neg, const int eta_index, const int phi_index, Identifier ID, float et, int prov, bool doPrint) const override
virtual void	ConnectSuperCellToTower (std::unique_ptr< gTowerContainer > &my_gTowerContainerRaw, int iGTower, Identifier ID, float et) const override
virtual int	FindTowerIDForSuperCell (int towereta, int towerphi, int nphi) const override
virtual void	PrintCellSpec (const CaloSampling::CaloSample sample, const int region, const int eta_index, const int phi_index, const int pos_neg, int towereta, int towerphi, int nphi, int towerID_Modifier, int iGTower, int prov, Identifier ID) const override
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
SG::ReadHandleKey< CaloCellContainer >	m_scellsCollectionSGKey {this, "SCell", "SCell", "SCell"}
SG::ReadHandleKey< xAOD::TriggerTowerContainer >	m_triggerTowerCollectionSGKey {this, "xODTriggerTowers", "xAODTriggerTowers", "xAODTriggerTowers"}
Gaudi::Property< bool >	m_apply_masking {this, "SCellMasking", false, "Applies masking. Only use for data"}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 25 of file gSuperCellTowerMapper.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

gSuperCellTowerMapper()

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

Definition at line 22 of file gSuperCellTowerMapper.cxx.

                                                                                                                :
  AthAlgTool(type,name,parent)
{
  declareInterface<IgSuperCellTowerMapper>(this);
 
}

~gSuperCellTowerMapper()

LVL1::gSuperCellTowerMapper::~gSuperCellTowerMapper ( )

virtual

Definition at line 29 of file gSuperCellTowerMapper.cxx.

{
}

Member Function Documentation

AssignSuperCellsToTowers()

StatusCode LVL1::gSuperCellTowerMapper::AssignSuperCellsToTowers ( std::unique_ptr< gTowerContainer > & my_gTowerContainerRaw ) const

overridevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 99 of file gSuperCellTowerMapper.cxx.

{
  bool doPrint = true;
 
  SG::ReadHandle<CaloCellContainer> scellsCollection(m_scellsCollectionSGKey/*,ctx*/);
  if(!scellsCollection.isValid()){
    ATH_MSG_FATAL("Could not retrieve scellsCollection " << m_scellsCollectionSGKey.key() );
    return StatusCode::FAILURE;
  }
 
 
  const CaloCell_Base_ID* idHelper = nullptr;
  ATH_CHECK( detStore()->retrieve (idHelper, "CaloCell_SuperCell_ID") );
 
  for (const CaloCell* cell : * scellsCollection){
 
    const CaloSampling::CaloSample sample = (cell)->caloDDE()->getSampling();
    const Identifier ID = (cell)->ID(); // super cell unique ID
    int region = idHelper->region(ID);
    int pos_neg = idHelper->pos_neg(ID);
    int eta_index = idHelper->eta(ID);
    const int phi_index = idHelper->phi(ID);
    float et = (cell)->energy()/std::cosh((cell)->eta());
    int prov = (cell)->provenance();
 
 
    // std::cout<< ID << " " << (cell)->eta() << " " << (cell)->phi() << " " << sample << " " << region << std::endl;
    // The following is to check if any SuperCells from data are permanently masked, and if so the masking is applied 
    int SCprov = prov&0xFFF;
    bool isMasked = (SCprov&0x80)==0x80;//prov looks like 0000 0000 1000 0000 if the cell is masked
    if (isMasked and m_apply_masking) et = 0;
 
    // I removed the LOCAL TO GLOBAL ETA INDEX PATCH for gFEX
    // Since in any case the SC assignment won't be regular, the eta and phi bins are combined directly in the FindAndConnectTower
    FindAndConnectTower(my_gTowerContainerRaw,sample,region,pos_neg,eta_index,phi_index,ID,et,prov,doPrint);
 
    
  }
 
  for (auto gTower : *my_gTowerContainerRaw)
  {
    gTower->setET();
 
  }
 
 
  return StatusCode::SUCCESS;
 
}

AssignTriggerTowerMapper()

StatusCode LVL1::gSuperCellTowerMapper::AssignTriggerTowerMapper ( std::unique_ptr< gTowerContainer > & my_gTowerContainerRaw ) const

overridevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 42 of file gSuperCellTowerMapper.cxx.

                                                                                                                       {
 
  static constexpr float delta_phi = M_PI/32;
  static constexpr float delta_eta = 0.1;
 
  SG::ReadHandle<xAOD::TriggerTowerContainer> triggerTowerCollection(m_triggerTowerCollectionSGKey/*,ctx*/);
  if(!triggerTowerCollection.isValid()){
    ATH_MSG_FATAL("Could not retrieve triggerTowerCollection " << m_triggerTowerCollectionSGKey.key() );
    return StatusCode::FAILURE;
  }
 
 
  for(auto eachTower : *triggerTowerCollection) {
 
    if(std::abs(eachTower->eta())<1.5 && eachTower->sampling()==1) {
      int tile_iphi = int(eachTower->phi()/delta_phi);
      int tower_iphi = tile_iphi/2;
      int tile_ieta = int(std::abs(eachTower->eta())/delta_eta);
      int tower_ieta = tile_ieta/2;
      int nphi = 32;
      int etaSign{-1};
      int towerID_Modifier{100000};
 
 
      if (eachTower->eta() > 0) {
        etaSign = 1;
        towerID_Modifier = 200000;
      }
 
      if(tower_ieta * etaSign == -7) {
        towerID_Modifier = 300000;
      }
 
      else if (tower_ieta * etaSign == 7) {
        towerID_Modifier = 400000;
      }
      int towerid = FindTowerIDForSuperCell(tower_ieta, tower_iphi, nphi) + towerID_Modifier;
 
      LVL1::gTower * targetTower;
      if((targetTower = my_gTowerContainerRaw->findTower(towerid)))
      {
        // Set the ET to the gTower, with layer 1 to specify this comes from the HAD
        targetTower->addET(static_cast<int>(eachTower->jepET()) * 1000., 1);
      } else
      {
        ATH_MSG_WARNING("\n==== gSuperCellTowerMapper ============ Tower id is officially unknown - it will be ignored. (Needs investigation).  Please report this!");
      }
    }
  }
  return StatusCode::SUCCESS;
}

ConnectSuperCellToTower()

void LVL1::gSuperCellTowerMapper::ConnectSuperCellToTower ( std::unique_ptr< gTowerContainer > & my_gTowerContainerRaw, int iGTower, Identifier ID, float et ) const

overrideprivatevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 150 of file gSuperCellTowerMapper.cxx.

                                                                                                                                                      {
 
  LVL1::gTower * tmpTower = my_gTowerContainerRaw->findTower(iGTower);
 
  if(tmpTower){
    tmpTower->setSCID(ID);
    tmpTower->addET(et, 0); // layer is always 0 (EM) for SuperCells
 
  }
 
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

FindAndConnectTower()

int LVL1::gSuperCellTowerMapper::FindAndConnectTower ( std::unique_ptr< gTowerContainer > & my_gTowerContainerRaw, CaloSampling::CaloSample sample, const int region, const int pos_neg, const int eta_index, const int phi_index, Identifier ID, float et, int prov, bool doPrint ) const

overrideprivatevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 162 of file gSuperCellTowerMapper.cxx.

{
 
  bool validcell = true;
 
  // We tell the gTower which supercell unique ID is in each tower
  int nphi = -1; // to be filled with number of phi bins
  int towereta = -99; // to be filled with tower eta value to help find it
  int towerphi = -99; // to be filled with tower phi value to help find it
  int iGTower = -99; // The tower ID which the supercell will be assigned to
  int towerID_Modifier = -999999999; // adjustable tower ID modifier to fit tower Id scheme
 
 
  //----------------------------------------------------------
 
  // Begin assigning supercells to Towers
 
  //----------------------------------------------------------
 
  switch (sample) {
      // E.M. Barrel (EMB)
    case CaloSampling::PreSamplerB: {
      // The following divides integers by 2 and truncates away the decimals (rounds DOWN to an integer)
      towereta = eta_index / 2; // Layer 0 has super cells which are 0.1 x 0.1 - gTower is 0.2 x 0.2
      towerphi = phi_index / 2; // phi is 32 for gTowers
 
      nphi = 32;
 
      if(eta_index == 14){ // special treatment for the tower where EMB and EMEC overlap (for gFEX: Eta Range 1.4-1.6)
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else{ // standard treatment for non-transition towers
        if(pos_neg < 0){ towerID_Modifier = 100000; }
        else if(pos_neg > 0){ towerID_Modifier = 200000; }
      }
 
      break;
    }
    case CaloSampling::EMB1:
    case CaloSampling::EMB2: {
 
      nphi = 32;
 
      if (region == 0) {
        towereta = eta_index / 8;
        towerphi = phi_index / 2;
        if(pos_neg < 0){ towerID_Modifier = 100000; }
        else if(pos_neg > 0){ towerID_Modifier = 200000; }
      }
      else if (region == 1){
        towereta = 7;
        towerphi = phi_index/2;
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else {
        ATH_MSG_DEBUG("[CaloSampling::EMB1 or CaloSampling::EMB2] -> invalid 'region' value: " << region << " (Under investigation) ");
      }
 
      switch(sample){
      case CaloSampling::EMB1: {
        break;
      }
      case CaloSampling::EMB2: {
        break;
      }
      default: {
        ATH_MSG_DEBUG("CaloSampling::EMBX -> invalid sample for assigning iCell value! " << sample << " (Under investigation) ");
        break;
      }
      }
 
      break;
    }
    case CaloSampling::EMB3: {
 
      nphi = 32;
      towereta = eta_index/2; // Analogous to PreSamplerB
      towerphi = phi_index/2;
 
 
      if(pos_neg < 0){ towerID_Modifier = 100000; }
      else if(pos_neg > 0){ towerID_Modifier = 200000; }
 
      break;
    }
      // E.M. End-cap (EMEC)
    case CaloSampling::PreSamplerE: {
 
      nphi = 32;
      // towereta = 7 + std::ceil(eta_index / 2.); //This is a special case, as there are only 3 cells in eta 1.56 (gFEX bin 7), 1.66 (gFEX bin 8), 1.76 (gFEX bin 8)
      if (eta_index < 1) {towereta = 7;}
        else {towereta = 8;}
      towerphi = phi_index / 2; // phi is 32 for gTowers
 
 
      if (towereta == 7){
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else if (towereta != 7){
        if(pos_neg < 0){ towerID_Modifier = 500000; }
        else if(pos_neg > 0){ towerID_Modifier = 600000; }
      }
 
      break;
    }
    case CaloSampling::EME1: {
 
      nphi = 32;
 
      switch (region) {
      case 0: {
 
        towereta = 7;
        towerphi = phi_index/2;
 
        break;
      }
      case 2: {
 
        if (eta_index < 4) {towereta = 7;}
        else {towereta = 8;}
        towerphi = phi_index/2;
 
        break;
      }
      case 3: {
 
        if (eta_index < 11) {towereta = 8;}
        else {towereta = 9;}
        towerphi = phi_index/2;
        break;
 
      }
      case 4: {
 
        if (eta_index < 8) {towereta = 10;}
        else {towereta = 11;}
        towerphi = phi_index/2;
 
        break;
      }
      case 5: {
 
        towereta = 11 ;
        towerphi = phi_index/2;
 
        break;
      }
      default: {
        ATH_MSG_DEBUG("CaloSampling::EME1 -> invalid 'region' value: " << region << " (Under investigation) ");
        break;
      }
      break;
      }
 
      if (towereta == 7){
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else if (towereta != 7){
        if(pos_neg < 0){ towerID_Modifier = 500000; }
        else if(pos_neg > 0){ towerID_Modifier = 600000; }
      }
 
      break;
    }
    case CaloSampling::EME2: {
 
      nphi = 32;
 
      switch (region) {
      case 0: { 
 
        if(abs(pos_neg) == 2){
          towereta = 7;
          towerphi = phi_index/2;
        }
        else if(abs(pos_neg) == 3){
          if (eta_index == 0){
            towereta = 12;
            towerphi = phi_index;
          }
          else if (eta_index == 1){
            towereta = 13;
            towerphi = phi_index;
          }
          else if (eta_index == 2){
            towereta = 14;
            towerphi = phi_index;
          }
        }
        break;
      }
      case 1: {
        if(abs(pos_neg) == 2){
          towerphi = phi_index/2;
          if (eta_index < 7) {towereta = 7;}
          else if (eta_index < 15) {towereta = 8;}
          else if (eta_index < 23) {towereta = 9;}
          else if (eta_index < 31) {towereta = 10;}
          else if (eta_index < 43) {towereta = 11;}
        }
        else if(abs(pos_neg) == 3){
          towerphi = phi_index;
          towereta = 15;
        }
        break;
      }
 
      default: {
        ATH_MSG_DEBUG("CaloSampling::EME2 -> invalid 'region' value: " << region << " (Under investigation) ");
        break;
      }
        break;
      }
 
      if (towereta == 7){
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else if (towereta == 12 || towereta == 13 || towereta == 14 || towereta == 15){
        if(pos_neg < 0){ towerID_Modifier = 700000; }
        else if(pos_neg > 0){ towerID_Modifier = 800000; }
      }
      else {
        if(pos_neg < 0){ towerID_Modifier = 500000; }
        else if(pos_neg > 0){ towerID_Modifier = 600000; }
      }
 
      break;
    }
    case CaloSampling::EME3: {
 
      nphi = 32;
      switch (region) {
      case 0: {
        if(abs(pos_neg) == 2){
          towerphi = phi_index/2;
          if (eta_index == 0){
            towereta = 7;
          }
          else if(eta_index < 3){
            towereta = 8;
          }
          else if(eta_index < 5){
            towereta = 9;
          }
          else if(eta_index < 7){
            towereta = 10;
          }
          else if(eta_index < 10){
            towereta = 11;
          }
 
        }
        else if(abs(pos_neg) == 3){
          if (eta_index == 0){
            towereta = 12;
            towerphi = phi_index;
          }
          else if (eta_index == 1){
            towereta = 13;
            towerphi = phi_index;
          }
          else if (eta_index == 2){
            towereta = 14;
            towerphi = phi_index;
          }
        }
        break;
      }
      case 1: { // we only have the case: abs(pos_neg) = 3
        towerphi = phi_index;
        towereta = 15;
        break;
      }
      default: {
        ATH_MSG_DEBUG("CaloSampling::EME3 -> invalid 'region' value: " << region << " (Under investigation) ");
        break;
      }
        break;
      }
 
      if (towereta == 7){
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else if (towereta == 12 || towereta == 13 || towereta == 14 || towereta == 15){
        if(pos_neg < 0){ towerID_Modifier = 700000; }
        else if(pos_neg > 0){ towerID_Modifier = 800000; }
      }
      else {
        if(pos_neg < 0){ towerID_Modifier = 500000; }
        else if(pos_neg > 0){ towerID_Modifier = 600000; }
      }
      break;
    }
      // Hadronic End-cap (HEC)
    case CaloSampling::HEC0:{
      nphi = 32;
      switch(region){
      case 0: {
        if(abs(pos_neg) == 2){
          towerphi = phi_index/2;
          if (eta_index == 0){
            towereta = 7;
          }
          else if(eta_index < 3){
            towereta = 8;
          }
          else if(eta_index < 5){
            towereta = 9;
          }
          else if(eta_index < 7){
            towereta = 10;
          }
          else if(eta_index < 10){
            towereta = 11;
          }
        }
        break;
 
      }
 
      case 1: {
        if(abs(pos_neg) == 2){
          if (eta_index == 0){
            towereta = 12;
            towerphi = phi_index;
          }
          else if(eta_index == 1){
            towereta = 13;
            towerphi = phi_index;
          }
          else if(eta_index == 2){
            towereta = 14;
            towerphi = phi_index;
          }
          else if(eta_index == 3){
            towereta = 15;
            towerphi = phi_index;
          }
        }
        break;
 
      }
 
      default: {
        ATH_MSG_DEBUG("CaloSampling::HEC0 -> invalid 'region' value: " << region << " (Under investigation) ");
        break;
      }
 
      }
 
      // Tower connection
      if (towereta == 7){
        nphi = 32;
        if(pos_neg < 0){ towerID_Modifier = 300000; }
        else if(pos_neg > 0){ towerID_Modifier = 400000; }
      }
      else if (towereta == 12 || towereta == 13 || towereta == 14 || towereta == 15){
        if(pos_neg < 0){ towerID_Modifier = 700000; }
        else if(pos_neg > 0){ towerID_Modifier = 800000; }
      }
      else {
        nphi = 32;
        if(pos_neg < 0){ towerID_Modifier = 500000; }
        else if(pos_neg > 0){ towerID_Modifier = 600000; }
      }
      
      break;
    }
    case CaloSampling::TileBar0:
    case CaloSampling::TileBar1:
    case CaloSampling::TileBar2: {
      REPORT_MESSAGE_WITH_CONTEXT (MSG::DEBUG, "gSuperCellTowerMapper") << "Supercell is from Tile Barrel - it will be ignored.";
      validcell = false;
      ATH_MSG_DEBUG("\n==== gSuperCellTowerMapper ============ Supercell is from Tile Barrel - it will be ignored.");
      break;
    }
    case CaloSampling::TileGap1:
    case CaloSampling::TileGap2:
    case CaloSampling::TileGap3: {
      ATH_MSG_DEBUG("\n==== gSuperCellTowerMapper ============ Supercell is from Tile Gap (ITC and scintillator) - it will be ignored.");
      validcell = false;
      break;
    }
    case CaloSampling::TileExt0:
    case CaloSampling::TileExt1:
    case CaloSampling::TileExt2: {
      ATH_MSG_DEBUG("\n==== gSuperCellTowerMapper ============ Supercell is from Tile Extended Barrel - it will be ignored.");
      validcell = false;
      break;
    }
    case CaloSampling::FCAL0:{
      nphi = 16;
      towerphi = phi_index;
      if(eta_index < 5){
        towereta = 16;
      }
      else if(eta_index < 8){
        towereta = 17;
      }
      else if(eta_index < 10){
        towereta = 18;
        }
      else if(eta_index < 12){
        towereta = 19;
      }
 
      if(pos_neg < 0){ towerID_Modifier = 900000; }
        else if(pos_neg > 0){ towerID_Modifier = 1000000; }
 
      break;
    }
    case CaloSampling::FCAL1:{
      nphi = 16;
      towerphi = phi_index;
      if (eta_index < 4){
        towereta = 16;
      }
      else if(eta_index < 6){
        towereta = 17;
      }
      else if(eta_index < 7){
        towereta = 18;
      }
      else if(eta_index < 8){
        towereta = 19;
      }
 
      if(pos_neg < 0){ towerID_Modifier = 900000; }
      else if(pos_neg > 0){ towerID_Modifier = 1000000; }
 
      break;
    }
    case CaloSampling::FCAL2:{
      nphi = 16;
      towerphi = phi_index;
      if (eta_index < 1){
        towereta = 16;
      }
      else if(eta_index < 3){
        towereta = 17;
      }
      else if(eta_index < 4){
        towereta = 18;
      }
 
      if(pos_neg < 0){ towerID_Modifier = 900000; }
      else if(pos_neg > 0){ towerID_Modifier = 1000000; }
      break;
    }
 
 
    case CaloSampling::MINIFCAL0:
    case CaloSampling::MINIFCAL1:
    case CaloSampling::MINIFCAL2:
    case CaloSampling::MINIFCAL3: {
      ATH_MSG_DEBUG("\n==== gSuperCellTowerMapper ============ Supercells is from MiniFCAL - it will be ignored.");
      validcell = false;
      break;
    }
    case CaloSampling::Unknown: {
      ATH_MSG_WARNING("\n==== gSuperCellTowerMapper ============ Supercell sampling is officially unknown - it will be ignored. (Needs investigation).  Please report this!");
      validcell = false;
      break;
    }
    default: {
      ATH_MSG_DEBUG("\n==== gSuperCellTowerMapper ============ Supercell has invalid CaloSampling value: " << sample << " (Needs investigation).  Please report this!");
      validcell = false;
      break;
    }
  }
 
  if(validcell){
    iGTower = FindTowerIDForSuperCell(towereta, towerphi, nphi) + towerID_Modifier;
    if(doPrint){
      PrintCellSpec(sample, region, eta_index, phi_index, pos_neg, towereta, towerphi, nphi, towerID_Modifier, iGTower, prov, ID);
    }
    ConnectSuperCellToTower(my_gTowerContainerRaw, iGTower, ID, et);
  }
 
  // END ITERATING OVER SUPER CELLS+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ +++++++++++++++++++++++++++++++++++++++++++++
 
  return 1;
 
}

FindTowerIDForSuperCell()

int LVL1::gSuperCellTowerMapper::FindTowerIDForSuperCell ( int towereta, int towerphi, int nphi ) const

overrideprivatevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 654 of file gSuperCellTowerMapper.cxx.

{
 
  return (towerphi + (nphi * towereta));
 
}

initialize()

StatusCode LVL1::gSuperCellTowerMapper::initialize ( )

overridevirtual

standard Athena-Algorithm method

Definition at line 33 of file gSuperCellTowerMapper.cxx.

{
 
  ATH_CHECK( m_scellsCollectionSGKey.initialize() );
  ATH_CHECK( m_triggerTowerCollectionSGKey.initialize() );
 
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

interfaceID()

const InterfaceID & LVL1::IgSuperCellTowerMapper::interfaceID ( )

inlinestaticinherited

Definition at line 46 of file IgSuperCellTowerMapper.h.

  {
    return IID_IgSuperCellTowerMapper;
  }

msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

PrintCellSpec()

void LVL1::gSuperCellTowerMapper::PrintCellSpec ( const CaloSampling::CaloSample sample, const int region, const int eta_index, const int phi_index, const int pos_neg, int towereta, int towerphi, int nphi, int towerID_Modifier, int iGTower, int prov, Identifier ID ) const

overrideprivatevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 661 of file gSuperCellTowerMapper.cxx.

{
 
  std::string sampleName = "";
 
  switch (sample) {
  case CaloSampling::PreSamplerB: { sampleName = "PreSamplerB"; break; }
  case CaloSampling::EMB1: { sampleName = "EMB1"; break; }
  case CaloSampling::EMB2: { sampleName = "EMB2"; break; }
  case CaloSampling::EMB3: { sampleName = "EMB3"; break; }
  case CaloSampling::PreSamplerE: { sampleName = "PreSamplerE"; break; }
  case CaloSampling::EME1: { sampleName = "EME1"; break; }
  case CaloSampling::EME2: { sampleName = "EME2"; break; }
  case CaloSampling::EME3: { sampleName = "EME3"; break; }
  case CaloSampling::HEC0: { sampleName = "HEC0"; break; }
  case CaloSampling::HEC1: { sampleName = "HEC1"; break; }
  case CaloSampling::HEC2: { sampleName = "HEC2"; break; }
  case CaloSampling::HEC3: { sampleName = "HEC3"; break; }
  case CaloSampling::TileBar0: { sampleName = "TileBar0"; break; }
  case CaloSampling::TileBar1: { sampleName = "TileBar1"; break; }
  case CaloSampling::TileBar2: { sampleName = "TileBar2"; break; }
  case CaloSampling::TileGap1: { sampleName = "TileGap1"; break; }
  case CaloSampling::TileGap2: { sampleName = "TileGap2"; break; }
  case CaloSampling::TileGap3: { sampleName = "TileGap3"; break; }
  case CaloSampling::TileExt0: { sampleName = "TileExt0"; break; }
  case CaloSampling::TileExt1: { sampleName = "TileExt1"; break; }
  case CaloSampling::TileExt2: { sampleName = "TileExt2"; break; }
  case CaloSampling::FCAL0: { sampleName = "FCAL0"; break; }
  case CaloSampling::FCAL1: { sampleName = "FCAL1"; break; }
  case CaloSampling::FCAL2: { sampleName = "FCAL2"; break; }
 
  default: {
    ATH_MSG_DEBUG("\n==== gSuperCellTowerMapper ============ Supercell has invalid CaloSampling value: " << sample << " (Needs investigation).  Please report this!");
    break;
  }
  }
 
  ATH_MSG_DEBUG("ASSIGNED CELL:::  CASE: " << sampleName
                << "\tSample:           " << sample
                << "\tRegion:           " << region
                << "\tEta_Index:        " << eta_index
                << "\tPhi_Index:        " << phi_index
                << "\tPosNeg:           " << pos_neg
            << "\tgTower Eta:       " << towereta
            << "\tgTower Phi:       " << towerphi
            << "\tgTower nphi:      " << nphi
            << "\tgTower modifier:  " << towerID_Modifier
                << "\tiGTower:          " << iGTower
                << "\tProvenance:       " << prov
                << "\tID:               " << ID
            << " ");
 
 
  return;
}

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

reset()

void LVL1::gSuperCellTowerMapper::reset ( )

overridevirtual

Implements LVL1::IgSuperCellTowerMapper.

Definition at line 94 of file gSuperCellTowerMapper.cxx.

                                 {
  return;
}

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_apply_masking

Gaudi::Property<bool> LVL1::gSuperCellTowerMapper::m_apply_masking {this, "SCellMasking", false, "Applies masking. Only use for data"}

private

Definition at line 44 of file gSuperCellTowerMapper.h.

{this, "SCellMasking", false, "Applies masking. Only use for data"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_scellsCollectionSGKey

SG::ReadHandleKey<CaloCellContainer> LVL1::gSuperCellTowerMapper::m_scellsCollectionSGKey {this, "SCell", "SCell", "SCell"}

private

Definition at line 41 of file gSuperCellTowerMapper.h.

{this, "SCell", "SCell", "SCell"};

m_triggerTowerCollectionSGKey

SG::ReadHandleKey<xAOD::TriggerTowerContainer> LVL1::gSuperCellTowerMapper::m_triggerTowerCollectionSGKey {this, "xODTriggerTowers", "xAODTriggerTowers", "xAODTriggerTowers"}

private

Definition at line 42 of file gSuperCellTowerMapper.h.

{this, "xODTriggerTowers", "xAODTriggerTowers", "xAODTriggerTowers"};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• gSuperCellTowerMapper.h
• gSuperCellTowerMapper.cxx