FCALHitsTestTool Class Reference

#include <FCALHitsTestTool.h>

Inheritance diagram for FCALHitsTestTool:

Collaboration diagram for FCALHitsTestTool:

Public Member Functions
	FCALHitsTestTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~FCALHitsTestTool ()
StatusCode	initialize ()
virtual StatusCode	processEvent () override

Protected Attributes
std::string	m_detname
double	m_edep_cut
TH1 *	m_lar_eta
TH1 *	m_lar_phi
TH2 *	m_lar_zr
TH2 *	m_lar_etaphi
TH2 *	m_lar_edep_zr
TH1 *	m_eta
TH1 *	m_phi
TH2 *	m_zr
TH2 *	m_etaphi
TH1 *	m_time
TH1 *	m_edep
TH1 *	m_log_edep
TH2 *	m_edep_zr
TH1 *	m_eta_cut1
TH1 *	m_phi_cut1
TH2 *	m_zr_cut1
TH2 *	m_etaphi_cut1
TH1 *	m_time_cut1
TH1 *	m_edep_cut1
TH2 *	m_edep_zr_cut1
TH1 *	m_etot
TH1 *	m_edep_eta
TH1 *	m_edep_phi
TH1 *	m_edep_z
TH1 *	m_edep_r
TH1 *	m_etot_eta
TH1 *	m_etot_phi
std::string	m_path {"/truth/"}
ServiceHandle< ITHistSvc >	m_histSvc {"THistSvc", "SimTestHisto"}

Private Attributes
SG::ReadCondHandleKey< CaloDetDescrManager >	m_caloMgrKey

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

Detailed Description

Definition at line 10 of file FCALHitsTestTool.h.

Constructor & Destructor Documentation

FCALHitsTestTool()

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

Definition at line 11 of file FCALHitsTestTool.cxx.

  : LArHitsTestTool(type, name, parent)
{
}

~FCALHitsTestTool()

virtual FCALHitsTestTool::~FCALHitsTestTool ( )

inlinevirtual

Definition at line 14 of file FCALHitsTestTool.h.

{}

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 FCALHitsTestTool::initialize ( )

virtual

Reimplemented from LArHitsTestTool.

Definition at line 16 of file FCALHitsTestTool.cxx.

{ 
  std::string origpath=m_path;
 
  m_path += "LAr/FCAL/";
 
  // variables specific to sub detector (root histo names have to be unique!)
  //_TH1D(m_eta,"FCAL_eta",25,-5.,5.);
  //_TH1D_WEIGHTED(m_time,"FCAL_time",100,0,25);
  _TH1D(m_etot,"FCAL_etot",100,0.,10.);
  //_TH1D_WEIGHTED(m_edep_eta,"FCAL_edep_eta",25,-5.,5.);
  //_TH1D_WEIGHTED(m_edep_z,"FCAL_edep_z",100,-5200.,5200.); 
  _TH1D_WEIGHTED(m_edep_r,"FCAL_edep_r",100,0.,500.);
  _TPROFILE(m_etot_eta,"FCAL_etot_eta",25,-3.,3.); 
 
  m_path = origpath;
 
  return LArHitsTestTool::initialize();
}

processEvent()

StatusCode LArHitsTestTool::processEvent ( )

overridevirtualinherited

Definition at line 116 of file LArHitsTestTool.cxx.

                                         {
  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey};
  ATH_CHECK(caloMgrHandle.isValid());
  const CaloDetDescrManager* caloMgr = *caloMgrHandle;
 
  std::string lArkey = "LArHit"+m_detname;
 
  double etot=0;
  const LArHitContainer* hitContainer = evtStore()->tryConstRetrieve<LArHitContainer>(lArkey);
  if(hitContainer) {
    for(const LArHit* larHit : *hitContainer) {
      const CaloDetDescrElement* ddElement = caloMgr->get_element(larHit->cellID());
 
      double eta = ddElement->eta();
      double phi = ddElement->phi();
      double radius = ddElement->r();
      double z  = ddElement->z();
      
      double energy = larHit->energy();
 
      m_lar_eta->Fill(eta);
      m_lar_phi->Fill(phi);
      m_lar_zr->Fill(z,radius);
      m_lar_etaphi->Fill(eta,phi);
      m_lar_edep_zr->Fill(z,radius,energy);
      
 
      m_eta->Fill(eta);
      m_phi->Fill(phi);
      m_time->Fill( larHit->time(),energy);
      m_edep->Fill( energy);
      m_log_edep->Fill( energy > 0 ? log(energy) : -1 );
      
      m_zr->Fill(z,radius);
      m_etaphi->Fill(eta,phi);
      m_edep_zr->Fill(z,radius,energy);
      
      m_edep_eta->Fill(eta,energy);
      m_edep_phi->Fill(phi,energy);
 
      m_edep_z->Fill(z,energy);
      m_edep_r->Fill(radius,energy);
 
      etot+=energy;
 
      if (larHit->energy()>m_edep_cut) {
    m_eta_cut1->Fill(eta);
    m_phi_cut1->Fill(phi);
    m_time_cut1->Fill( larHit->time(), larHit->energy());
    m_edep_cut1->Fill( larHit->energy());
    
    m_zr_cut1->Fill(z,radius);
    m_etaphi_cut1->Fill(eta,phi);
    m_edep_zr_cut1->Fill(z,radius,larHit->energy());
      }
    }
  }
 
  //For FastCaloSim Container with _Fast postfix, just try to retrieve, if exist, collect information from this container, it not just skip it.
  
  LArHitContainer::const_iterator hi_fast;
  const std::string lArkey_fast="LArHit"+m_detname+"_Fast";
  
  const LArHitContainer* iter_fast;
  if( evtStore()->contains<LArHitContainer>(lArkey_fast) &&
      (evtStore()->retrieve(iter_fast,lArkey_fast)).isSuccess())
  {
    ATH_MSG_DEBUG ( "Read hit info from FastCaloSim Container" );
    for(hi_fast=(*iter_fast).begin();hi_fast!=(*iter_fast).end();++hi_fast) 
    {
      const LArHit* larHit = *hi_fast;
      const CaloDetDescrElement* ddElement = caloMgr->get_element(larHit->cellID());
      double eta=ddElement->eta();
      double phi=ddElement->phi();
      double radius=ddElement->r();
      double z=ddElement->z();
      double energy=larHit->energy();
 
      m_lar_eta->Fill(eta);
      m_lar_phi->Fill(phi);
      m_lar_zr->Fill(z,radius);
      m_lar_etaphi->Fill(eta,phi);
      m_lar_edep_zr->Fill(z,radius,energy);
 
      m_eta->Fill(eta);
      m_phi->Fill(phi);
      m_time->Fill((*hi_fast)->time(),energy);
      m_edep->Fill(energy);
      m_log_edep->Fill( log(energy));
 
      m_zr->Fill(z,radius);
      m_etaphi->Fill(eta,phi);
      m_edep_zr->Fill(z,radius,energy);
 
      m_edep_eta->Fill(eta,energy);
      m_edep_phi->Fill(phi,energy);
 
      m_edep_z->Fill(z,energy);
      m_edep_r->Fill(radius,energy);
 
      etot+=energy;
 
      if((*hi_fast)->energy() > m_edep_cut)
      {
        m_eta_cut1->Fill(eta); 
        m_phi_cut1->Fill(phi);
        m_time_cut1->Fill((*hi_fast)->time(),(*hi_fast)->energy());
        m_edep_cut1->Fill((*hi_fast)->energy());
 
        m_zr_cut1->Fill(z,radius);
        m_etaphi_cut1->Fill(eta,phi);
        m_edep_zr_cut1->Fill(z,radius,(*hi_fast)->energy());
      }
    }
  }
  auto primary = getPrimary();
  if (primary) {
    m_etot_eta->Fill(primary->momentum().eta(),etot);
    m_etot_phi->Fill(primary->momentum().phi(),etot);
  }
  m_etot->Fill(etot);
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> LArHitsTestTool::m_caloMgrKey

privateinherited

Initial value:

{ this
      , "CaloDetDescrManager"
      , "CaloDetDescrManager"
      , "SG Key for CaloDetDescrManager in the Condition Store" }

Definition at line 48 of file LArHitsTestTool.h.

m_detname

std::string LArHitsTestTool::m_detname

protectedinherited

Definition at line 24 of file LArHitsTestTool.h.

m_edep

TH1 * LArHitsTestTool::m_edep

protectedinherited

Definition at line 35 of file LArHitsTestTool.h.

m_edep_cut

double LArHitsTestTool::m_edep_cut

protectedinherited

Definition at line 25 of file LArHitsTestTool.h.

m_edep_cut1

TH1 * LArHitsTestTool::m_edep_cut1

protectedinherited

Definition at line 39 of file LArHitsTestTool.h.

m_edep_eta

TH1* LArHitsTestTool::m_edep_eta

protectedinherited

Definition at line 43 of file LArHitsTestTool.h.

m_edep_phi

TH1 * LArHitsTestTool::m_edep_phi

protectedinherited

Definition at line 43 of file LArHitsTestTool.h.

m_edep_r

TH1 * LArHitsTestTool::m_edep_r

protectedinherited

Definition at line 44 of file LArHitsTestTool.h.

m_edep_z

TH1* LArHitsTestTool::m_edep_z

protectedinherited

Definition at line 44 of file LArHitsTestTool.h.

m_edep_zr

TH2* LArHitsTestTool::m_edep_zr

protectedinherited

Definition at line 36 of file LArHitsTestTool.h.

m_edep_zr_cut1

TH2* LArHitsTestTool::m_edep_zr_cut1

protectedinherited

Definition at line 40 of file LArHitsTestTool.h.

m_eta

TH1* LArHitsTestTool::m_eta

protectedinherited

Definition at line 33 of file LArHitsTestTool.h.

m_eta_cut1

TH1* LArHitsTestTool::m_eta_cut1

protectedinherited

Definition at line 37 of file LArHitsTestTool.h.

m_etaphi

TH2 * LArHitsTestTool::m_etaphi

protectedinherited

Definition at line 34 of file LArHitsTestTool.h.

m_etaphi_cut1

TH2 * LArHitsTestTool::m_etaphi_cut1

protectedinherited

Definition at line 38 of file LArHitsTestTool.h.

m_etot

TH1* LArHitsTestTool::m_etot

protectedinherited

Definition at line 42 of file LArHitsTestTool.h.

m_etot_eta

TH1* LArHitsTestTool::m_etot_eta

protectedinherited

Definition at line 45 of file LArHitsTestTool.h.

m_etot_phi

TH1 * LArHitsTestTool::m_etot_phi

protectedinherited

Definition at line 45 of file LArHitsTestTool.h.

m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

m_key

std::string SimTestToolBase::m_key

protectedinherited

The MC truth key.

Definition at line 34 of file SimTestToolBase.h.

m_lar_edep_zr

TH2* LArHitsTestTool::m_lar_edep_zr

protectedinherited

Definition at line 30 of file LArHitsTestTool.h.

m_lar_eta

TH1* LArHitsTestTool::m_lar_eta

protectedinherited

Definition at line 28 of file LArHitsTestTool.h.

m_lar_etaphi

TH2 * LArHitsTestTool::m_lar_etaphi

protectedinherited

Definition at line 29 of file LArHitsTestTool.h.

m_lar_phi

TH1 * LArHitsTestTool::m_lar_phi

protectedinherited

Definition at line 28 of file LArHitsTestTool.h.

m_lar_zr

TH2* LArHitsTestTool::m_lar_zr

protectedinherited

Definition at line 29 of file LArHitsTestTool.h.

m_log_edep

TH1 * LArHitsTestTool::m_log_edep

protectedinherited

Definition at line 35 of file LArHitsTestTool.h.

m_path

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

protectedinherited

Definition at line 34 of file SimTestHisto.h.

m_phi

TH1 * LArHitsTestTool::m_phi

protectedinherited

Definition at line 33 of file LArHitsTestTool.h.

m_phi_cut1

TH1 * LArHitsTestTool::m_phi_cut1

protectedinherited

Definition at line 37 of file LArHitsTestTool.h.

m_time

TH1* LArHitsTestTool::m_time

protectedinherited

Definition at line 35 of file LArHitsTestTool.h.

m_time_cut1

TH1* LArHitsTestTool::m_time_cut1

protectedinherited

Definition at line 39 of file LArHitsTestTool.h.

m_zr

TH2* LArHitsTestTool::m_zr

protectedinherited

Definition at line 34 of file LArHitsTestTool.h.

m_zr_cut1

TH2* LArHitsTestTool::m_zr_cut1

protectedinherited

Definition at line 38 of file LArHitsTestTool.h.

---

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

• FCALHitsTestTool.h
• FCALHitsTestTool.cxx