PFCalcRadialEnergyProfiles Class Reference

#include <PFCalcRadialEnergyProfiles.h>

Inheritance diagram for PFCalcRadialEnergyProfiles:

Collaboration diagram for PFCalcRadialEnergyProfiles:

Public Member Functions
	PFCalcRadialEnergyProfiles ()
	~PFCalcRadialEnergyProfiles ()=default
void	calculate (const PFData &data) const
void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
std::string	m_nm
	Message source name. More...
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Detailed Description

Definition at line 12 of file PFCalcRadialEnergyProfiles.h.

Constructor & Destructor Documentation

PFCalcRadialEnergyProfiles()

PFCalcRadialEnergyProfiles::PFCalcRadialEnergyProfiles ( )

inline

Definition at line 16 of file PFCalcRadialEnergyProfiles.h.

    : asg::AsgMessaging("PFCalcRadialEnergyProfiles")
  {
  }

~PFCalcRadialEnergyProfiles()

PFCalcRadialEnergyProfiles::~PFCalcRadialEnergyProfiles ( )

default

Member Function Documentation

calculate()

void PFCalcRadialEnergyProfiles::calculate

(

const PFData &

data

)

const

Definition at line 11 of file PFCalcRadialEnergyProfiles.cxx.

                                                                  {
 
  ATH_MSG_DEBUG("Accessed radial energy profile function");
 
  for (auto thisEflowCaloObject : *data.caloObjects){
 
    //If there are no clusters available then we cannot calculate any calorimeter shower profiles
    if (thisEflowCaloObject->nClusters() < 1 ) continue;
 
    const std::vector<std::pair<eflowTrackClusterLink*,std::pair<float,float> > > matchedTrackList = thisEflowCaloObject->efRecLink();
 
    for( auto track: matchedTrackList){
 
      eflowRecTrack* efRecTrack = (track.first)->getTrack();
          
      std::vector<eflowRecCluster*> matchedClusters;
      matchedClusters.clear();
      std::vector<eflowTrackClusterLink*> links = efRecTrack->getClusterMatches();
      for (auto *thisEFlowTrackClusterLink : links) matchedClusters.push_back(thisEFlowTrackClusterLink->getCluster());
 
      std::vector<std::pair<xAOD::CaloCluster*, bool> > clusterSubtractionList;
      clusterSubtractionList.reserve(matchedClusters.size());
for (auto *thisEFlowRecCluster : matchedClusters) clusterSubtractionList.emplace_back(thisEFlowRecCluster->getCluster(),false);
 
      eflowCellList calorimeterCellList;
      eflowSubtract::Subtractor::makeOrderedCellList(efRecTrack->getTrackCaloPoints(),clusterSubtractionList,calorimeterCellList);
      
      std::vector<int> layerToStoreVector;
      std::vector<float> radiusToStoreVector;
      std::vector<float> avgEdensityToStoreVector;
 
      //Loop over calorimeter layers and in each layer we will calculate radial energy profiles.
      for (int i=0; i < eflowCalo::nRegions ;i++){
    
        eflowCaloENUM layer = (eflowCaloENUM)i;
        ATH_MSG_DEBUG("layer is "<<layer);
        double ringThickness = eflowRingThicknesses::ringThickness((eflowCaloENUM)i);
        ATH_MSG_DEBUG("ring thickness is "<<ringThickness);
    
        double eta_extr = calorimeterCellList.etaFF(layer);
        ATH_MSG_DEBUG("extrapolated eta ["<<layer<<"] is "<<eta_extr);
        double phi_extr = calorimeterCellList.phiFF(layer);
        ATH_MSG_DEBUG("extrapolated phi ["<<layer<<"] is "<<phi_extr);
    
        if (eta_extr == -999.0){
          continue;
        }
                
        //Loop over rings of calorimeter cells going out in increasing radius from the track impact point in this calorimeter layer
        for (unsigned  int indexOfRing = 0; indexOfRing < 100; indexOfRing++){
 
          CellIt beginRing = calorimeterCellList.getLowerBound((eflowCaloENUM)i, ringThickness*(indexOfRing));
          if(beginRing == calorimeterCellList.end()) break;
 
          int totalCellsinRing = 0;
          double totalEnergyPerRing = 0;
          double energyDensityPerRing = 0;
          double averageEnergyDensityPerRing = 0;
 
          //Get the list of calorimeter cells in this cell ring
          std::vector<std::pair<const CaloCell*,int> > tempVector = (*beginRing).second;
 
          //Loop over the calorimeter cells in this cell ring and calculate total and average energy densities in this cell ring.
          for (auto thisPair : tempVector){
            const CaloDetDescrElement* DDE = (thisPair.first)->caloDDE();
              CaloCell_ID::CaloSample sampling = DDE->getSampling();
 
            if(eflowCalo::translateSampl(sampling)!=(eflowCaloENUM)i) break;
 
            ATH_MSG_DEBUG(" cell eta and phi are " << (thisPair.first)->eta() << " and " << (thisPair.first)->phi() << " with index " << thisPair.second << " and sampling of " << sampling);
              ATH_MSG_DEBUG(" cell energy is " << (thisPair.first)->energy());
            
            totalCellsinRing += 1;
        
            totalEnergyPerRing += (thisPair.first)->energy();
            double totalEnergyCell = (thisPair.first)->energy();
            ATH_MSG_DEBUG(" Total E per Cell is " << totalEnergyCell);
            ATH_MSG_DEBUG(" Total E per Ring is " << totalEnergyPerRing);
        
            double cellVolume = DDE->volume();
            ATH_MSG_DEBUG(" cell volume is " << cellVolume/1000.);
        
            double energyDensityCell = totalEnergyCell/(cellVolume/1000.);
            ATH_MSG_DEBUG(" E density per Cell is " << energyDensityCell);
            ATH_MSG_DEBUG(" Initial added E density per Cell is " << energyDensityPerRing);
            energyDensityPerRing += energyDensityCell;
            ATH_MSG_DEBUG(" Final added E density per Cell is " << energyDensityPerRing);
            averageEnergyDensityPerRing = energyDensityPerRing/((totalCellsinRing)*(efRecTrack->getTrack()->e()/1000.));
 
          }
 
          ATH_MSG_DEBUG(" track eta is " << efRecTrack->getTrack()->eta());
          ATH_MSG_DEBUG(" track E is " << efRecTrack->getTrack()->e()/1000.);
          ATH_MSG_DEBUG(" Average E density per Ring is " << averageEnergyDensityPerRing);
      
          //Fill up the vectors with energy density, layer and ring radii
          if (averageEnergyDensityPerRing != 0){
            avgEdensityToStoreVector.push_back(averageEnergyDensityPerRing);
            int layerToStore = (eflowCaloENUM)i;
            ATH_MSG_DEBUG("layerToStore is "<< layerToStore);
            layerToStoreVector.push_back(layerToStore);
            double radiusToStore = (indexOfRing)*ringThickness;
            ATH_MSG_DEBUG("radiusToStore is "<< radiusToStore);
            radiusToStoreVector.push_back(radiusToStore);
          }
          else {ATH_MSG_DEBUG("averageEnergyDensityPerRing = 0");}
 
        }//loop over rings of calorimeter cells
 
      }//loop over the calorimeter layers
 
      //Add the vectors with energy density, layer and ring radii to this eflowRecTrack
      efRecTrack->setLayerCellOrderVector(layerToStoreVector);
      efRecTrack->setRadiusCellOrderVector(radiusToStoreVector);
      efRecTrack->setAvgEDensityCellOrderVector(avgEdensityToStoreVector);
 
    }//loop over tracks matched to clusters
 
  }//loop on eflowCaloObjectContainer
 
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  m_lvl = m_imsg ?
    static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
    MSG::INFO;
}

msg() [1/2]

MsgStream & asg::AsgMessaging::msg ( ) const

inherited

The standard message stream.

Returns

A reference to the default message stream of this object.

Definition at line 49 of file AsgMessaging.cxx.

                                      {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg();
#else // not XAOD_STANDALONE
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msg() [2/2]

MsgStream & asg::AsgMessaging::msg ( const MSG::Level  lvl ) const

inherited

The standard message stream.

Parameters

lvl	The message level to set the stream to

Returns

A reference to the default message stream, set to level "lvl"

Definition at line 57 of file AsgMessaging.cxx.

                                                          {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg( lvl );
#else // not XAOD_STANDALONE
      m_msg << lvl;
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msgLvl()

bool asg::AsgMessaging::msgLvl ( const MSG::Level  lvl ) const

inherited

Test the output level of the object.

Parameters

lvl	The message level to test against

Returns

boolean Indicting if messages at given level will be printed

true If messages at level "lvl" will be printed

Definition at line 41 of file AsgMessaging.cxx.

                                                       {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msgLvl( lvl );
#else // not XAOD_STANDALONE
      return m_msg.msgLevel( lvl );
#endif // not XAOD_STANDALONE
   }

setLevel()

void AthMessaging::setLevel ( MSG::Level  lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

---

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

• PFCalcRadialEnergyProfiles.h
• PFCalcRadialEnergyProfiles.cxx