TileHitsTestTool Class Reference

#include <TileHitsTestTool.h>

Inheritance diagram for TileHitsTestTool:

Collaboration diagram for TileHitsTestTool:

Public Member Functions
	TileHitsTestTool (const std::string &type, const std::string &name, const IInterface *parent)
StatusCode	initialize ()
StatusCode	processEvent ()

Protected Attributes
std::string	m_path {"/truth/"}
ServiceHandle< ITHistSvc >	m_histSvc {"THistSvc", "SimTestHisto"}

Private Attributes
const TileID *	m_tileID
const TileTBID *	m_tileTBID
const TileDetDescrManager *	m_tileMgr
TH1 *	m_tile_eta
TH1 *	m_tile_phi
TH1 *	m_tile_energy
TH1 *	m_tile_log_energy
TH1 *	m_tile_time
TH2 *	m_tile_rz
TH2 *	m_tile_etaphi
TH2 *	m_tile_energyphi
TH2 *	m_tile_energyeta
TH1 *	m_mbts_side
TH1 *	m_mbts_eta
TH1 *	m_mbts_phi
TH2 *	m_mbts_sidetaphi
TH1 *	m_etot
bool	m_testMBTS

structors and AlgTool implementation
std::string	m_key
	The MC truth key.
HepMC::ConstGenParticlePtr	getPrimary ()

Detailed Description

Definition at line 15 of file TileHitsTestTool.h.

Constructor & Destructor Documentation

TileHitsTestTool()

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

Definition at line 18 of file TileHitsTestTool.cxx.

  : SimTestToolBase(type, name, parent),
    m_tileID(0), m_tileTBID(0), m_tileMgr(0),
    m_tile_eta(0), m_tile_phi(0), m_tile_energy(0), m_tile_log_energy(0), m_tile_time(0),
    m_tile_rz(0), m_tile_etaphi(0), m_tile_energyphi(0), m_tile_energyeta(0),
    m_mbts_side(0), m_mbts_eta(0), m_mbts_phi(0),
    m_mbts_sidetaphi(0),
    m_etot(0),
    m_testMBTS(true)
{
  declareProperty("TestMBTS",m_testMBTS,"");
  }

Member Function Documentation

getPrimary()

HepMC::ConstGenParticlePtr SimTestToolBase::getPrimary ( )

protectedinherited

Definition at line 20 of file SimTestToolBase.cxx.

{
  const McEventCollection* mcCollection;
  if (evtStore()->retrieve(mcCollection,m_key).isSuccess()) {
    DataVector<HepMC::GenEvent>::const_iterator e;
    for (e=mcCollection->begin();e!=mcCollection->end(); ++e) {
      for (auto p : (**e)) {
    if (!HepMC::is_simulation_particle(p)) {
      return p;
    }
      }
    }
  }
  return 0;
}

initialize()

StatusCode TileHitsTestTool::initialize ( )

virtual

Reimplemented from SimTestToolBase.

Definition at line 32 of file TileHitsTestTool.cxx.

                                        {
  CHECK(detStore()->retrieve(m_tileMgr));
  CHECK(detStore()->retrieve(m_tileID));
  CHECK(detStore()->retrieve(m_tileTBID));
 
  std::string detName("Tile"); //could configure this in the future?
  m_path+=detName+"/";
 
  // TileCal histograms
  _TH1D(m_tile_eta,"tile_eta",17,-1.65,1.65);
  _TH1D(m_tile_phi,"tile_phi",16,-3.2,3.2);
  _TH1D(m_tile_energy,"tile_energy",50,0,100);
  _TH1D(m_tile_log_energy,"tile_log_energy",100,-14.,8.);
  _SET_TITLE(m_tile_log_energy,"logarithm of energy deposits","log(energy [MeV])","dN/dlog(energy [MeV])");
  _TH1D(m_tile_time,"tile_time",100,0,250);
 
  _TH2D(m_tile_rz,"tile_rz",16,-6200.,6200.,25,1500.,4000.);
  _TH2D(m_tile_etaphi,"tile_etaphi",65,-1.625,1.625,129,-3.175,3.175);
  _TH2D(m_tile_energyphi,"tile_energyphi",100,0.,10.,129,-3.175,3.175);
  _TH2D(m_tile_energyeta,"tile_energyeta",100,0.,10.,129,-3.175,3.175);
 
  // MBTS histograms
  _TH1D(m_mbts_side,"mbts_side",3,-1.5,1.5);
  _TH1D(m_mbts_eta,"mbts_eta",2,-0.5,1.5);
  _TH1D(m_mbts_phi,"mbts_phi",8,-0.5,7.5);
  _TH2D(m_mbts_sidetaphi,"mbts_sidetaphi",5,-2.5,2.5,7,-0.5,7.5);
 
  _TH1D(m_etot,detName+"_etot",100,0.,500.);
 
  return StatusCode::SUCCESS;
}

processEvent()

StatusCode TileHitsTestTool::processEvent

(

)

Definition at line 65 of file TileHitsTestTool.cxx.

                                          {
  // make sure we still have a TileMgr (SimFlags.ReleaseGeoModel=False)
  CHECK(detStore()->retrieve(m_tileMgr));
 
  const TileHitVector* hitVec = nullptr;
 
  double etot=0.;
  if (evtStore()->retrieve(hitVec,"TileHitVec") == StatusCode::SUCCESS) {
    for (const TileHit& hit : *hitVec) {
 
      Identifier pmt_id = hit.identify();
      Identifier cell_id = m_tileID->cell_id(pmt_id);
 
      const CaloDetDescrElement *ddElement = m_tileMgr->get_cell_element(cell_id);
      if(ddElement) // skip E4' cells which are not yet described in TileDetDescrManager.
        {
          double eta = ddElement->eta();
          double phi = ddElement->phi();
          double radius = ddElement->r();
          double z  = ddElement->z();
 
          int pmt = m_tileID->pmt(pmt_id);
          if (pmt>0) phi += ddElement->dphi()/2.; // to see clearly two PMTs in a cell
 
          // loop over subhits
          double energy=0;
          //    ATH_MSG(INFO) << "TileHit size :" <<hit.size()<<endmsg;
          for (int i=0; i<hit.size();++i) {
            energy+=hit.energy(i);
            m_tile_energy->Fill(hit.energy(i));
            m_tile_log_energy->Fill( std::log(hit.energy(i)) );
            m_tile_time->Fill(hit.time(i));
          }
          etot+=energy;
 
          m_tile_eta->Fill(eta);
          m_tile_phi->Fill(phi);
          m_tile_rz->Fill(z,radius);
          m_tile_etaphi->Fill(eta,phi);
          m_tile_energyphi->Fill(energy,phi);
          m_tile_energyeta->Fill(energy,eta);
        }
    }
  }
  //Check the Hit container with postfix _Fast which generated with FastCaloSim
  const TileHitVector *hitVec_fast = nullptr;
  if(  evtStore()->contains<TileHitVector>("TileHitVec_Fast") &&
       (evtStore()->retrieve(hitVec_fast,"TileHitVec_Fast")).isSuccess())
    {
      ATH_MSG_DEBUG ( "Retrieve FastCaloSim container TileHitVec_Fast." );
      for (const TileHit& hit : *hitVec_fast)
        {
          Identifier pmt_id=hit.identify();
          Identifier cell_id=m_tileID->cell_id(pmt_id);
          const CaloDetDescrElement *ddElement=m_tileMgr->get_cell_element(cell_id);
          if(ddElement) // skip E4' cells which are not yet described in TileDetDescrManager.
            {
              double eta=ddElement->eta();
              double phi=ddElement->phi();
              double radius=ddElement->r();
              double z=ddElement->z();
 
              int pmt=m_tileID->pmt(pmt_id);
              if(pmt>0) phi+=ddElement->dphi()/2.;
              double energy=0;
              for(int i=0;i<hit.size();++i)
                {
                  energy+=hit.energy(i);
                  m_tile_energy->Fill(hit.energy(i));
                  m_tile_log_energy->Fill( std::log(hit.energy(i))) ;
                  m_tile_time->Fill(hit.time(i));
                }
              etot+=energy;
 
              m_tile_eta->Fill(eta);
              m_tile_phi->Fill(phi);
              m_tile_rz->Fill(z,radius);
              m_tile_etaphi->Fill(eta,phi);
              m_tile_energyphi->Fill(energy,phi);
              m_tile_energyeta->Fill(energy,eta);
            }
        }
    }
 
  if(m_testMBTS) {
    if (evtStore()->retrieve(hitVec,"MBTSHits") == StatusCode::SUCCESS) {
      for (const TileHit& hit : *hitVec) {
        Identifier mbts_id = hit.identify();
        double side = m_tileTBID->side(mbts_id); // -1 or +1
        double ieta = m_tileTBID->eta(mbts_id);  // 0 for inner cell, 1 for outer cell
        double iphi = m_tileTBID->phi(mbts_id);  // 0-7, cell 0 at phi=0
        m_mbts_side->Fill(side);
        m_mbts_eta->Fill(ieta);
        m_mbts_phi->Fill(iphi);
        m_mbts_sidetaphi->Fill((ieta+1)*side,iphi);
      }
    }
  }
  m_etot->Fill(etot);
 
 
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_etot

TH1* TileHitsTestTool::m_etot

private

Definition at line 38 of file TileHitsTestTool.h.

m_histSvc

ServiceHandle<ITHistSvc> SimTestHisto::m_histSvc {"THistSvc", "SimTestHisto"}

protectedinherited

Definition at line 35 of file SimTestHisto.h.

{"THistSvc", "SimTestHisto"};

m_key

std::string SimTestToolBase::m_key

protectedinherited

The MC truth key.

Definition at line 34 of file SimTestToolBase.h.

m_mbts_eta

TH1 * TileHitsTestTool::m_mbts_eta

private

Definition at line 35 of file TileHitsTestTool.h.

m_mbts_phi

TH1 * TileHitsTestTool::m_mbts_phi

private

Definition at line 35 of file TileHitsTestTool.h.

m_mbts_side

TH1* TileHitsTestTool::m_mbts_side

private

Definition at line 35 of file TileHitsTestTool.h.

m_mbts_sidetaphi

TH2* TileHitsTestTool::m_mbts_sidetaphi

private

Definition at line 36 of file TileHitsTestTool.h.

m_path

std::string SimTestHisto::m_path {"/truth/"}

protectedinherited

Definition at line 34 of file SimTestHisto.h.

{"/truth/"};

m_testMBTS

bool TileHitsTestTool::m_testMBTS

private

Definition at line 40 of file TileHitsTestTool.h.

m_tile_energy

TH1 * TileHitsTestTool::m_tile_energy

private

Definition at line 32 of file TileHitsTestTool.h.

m_tile_energyeta

TH2 * TileHitsTestTool::m_tile_energyeta

private

Definition at line 33 of file TileHitsTestTool.h.

m_tile_energyphi

TH2 * TileHitsTestTool::m_tile_energyphi

private

Definition at line 33 of file TileHitsTestTool.h.

m_tile_eta

TH1* TileHitsTestTool::m_tile_eta

private

Definition at line 32 of file TileHitsTestTool.h.

m_tile_etaphi

TH2 * TileHitsTestTool::m_tile_etaphi

private

Definition at line 33 of file TileHitsTestTool.h.

m_tile_log_energy

TH1 * TileHitsTestTool::m_tile_log_energy

private

Definition at line 32 of file TileHitsTestTool.h.

m_tile_phi

TH1 * TileHitsTestTool::m_tile_phi

private

Definition at line 32 of file TileHitsTestTool.h.

m_tile_rz

TH2* TileHitsTestTool::m_tile_rz

private

Definition at line 33 of file TileHitsTestTool.h.

m_tile_time

TH1 * TileHitsTestTool::m_tile_time

private

Definition at line 32 of file TileHitsTestTool.h.

m_tileID

const TileID* TileHitsTestTool::m_tileID

private

Definition at line 28 of file TileHitsTestTool.h.

m_tileMgr

const TileDetDescrManager* TileHitsTestTool::m_tileMgr

private

Definition at line 30 of file TileHitsTestTool.h.

m_tileTBID

const TileTBID* TileHitsTestTool::m_tileTBID

private

Definition at line 29 of file TileHitsTestTool.h.

---

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

• TileHitsTestTool.h
• TileHitsTestTool.cxx