ALFA_SimHitsTestTool Class Reference

#include <ALFA_SimHitsTestTool.h>

Inheritance diagram for ALFA_SimHitsTestTool:

Collaboration diagram for ALFA_SimHitsTestTool:

Public Member Functions
	ALFA_SimHitsTestTool (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
TH2 *	m_hit_map_h [8]
TH2 *	m_E_layer_sum_h [8]
TH1 *	m_E_full_sum_h [8]
TH1 *	m_E_gvs {}

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

Detailed Description

Definition at line 11 of file ALFA_SimHitsTestTool.h.

Constructor & Destructor Documentation

ALFA_SimHitsTestTool()

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

Definition at line 20 of file ALFA_SimHitsTestTool.cxx.

  : SimTestToolBase(type, name, parent)
{
  for (int i(0); i<8; ++i) {
    m_E_full_sum_h[i]=nullptr;
    m_E_layer_sum_h[i]=nullptr;
    m_hit_map_h[i]=nullptr;
  }
}

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

virtual

Reimplemented from SimTestToolBase.

Definition at line 30 of file ALFA_SimHitsTestTool.cxx.

{
  m_path+="ALFA/";
 
  _TH1D(m_E_gvs,"edep_full_in_gvs", 8, -0.5, 7.5);
  _SET_TITLE(m_E_gvs, "Energy deposit in GVS","Station","E (MeV)");
 
  std::stringstream s;
  for (unsigned int j=0; j<8; j++) {
    s.str("");
    s << j+1;
 
    _TH1D(m_E_full_sum_h[j],("edep_full_in_det_no."+s.str()).c_str(), 100, 0., 5.);
    _SET_TITLE(m_E_full_sum_h[j], "Energy deposit in full detector","E (MeV)","N");
 
    _TH2D(m_E_layer_sum_h[j],("edep_per_layer_detector_no."+s.str()).c_str(), 100, 0., 0.5, 22, 0., 21.);
    _SET_TITLE(m_E_layer_sum_h[j], "Energy deposit in layers","E (MeV)","layer");
 
    _TH2D(m_hit_map_h[j],("hitmap_det_no."+s.str()).c_str(), 22, 0., 21., 66, 0., 65.);
    _SET_TITLE(m_hit_map_h[j], "hit map - layer vs fiber","layer","fiber");
  }
 
  return StatusCode::SUCCESS;
}

processEvent()

StatusCode ALFA_SimHitsTestTool::processEvent

(

)

Definition at line 55 of file ALFA_SimHitsTestTool.cxx.

                                              {
 
  int fiber, plate, sign, station;
 
  double E_fiber_sum[8][10][64][2];
  double E_full_sum[8];
  double E_layer_sum[8][20];
 
  // cleaning
  const int station_max{8};
  const int plate_max{10};
  const int fiber_max{64};
 
  for ( int l = 0; l < station_max; l++ ) {
    E_full_sum[l] = 0.;
    for ( int i = 0; i < plate_max; i++ ) {
      E_layer_sum[l][i] = 0.;
      E_layer_sum[l][i+plate_max] = 0.;
      for ( int j = 0; j < fiber_max;  j++ ) {
        for ( int k = 0; k < 2; k++ ) {
          E_fiber_sum[l][i][j][k] = 0.;
        }
      }
    }
  }
 
 
  const ALFA_HitCollection* coll_handle = nullptr;
  CHECK( evtStore()->retrieve(coll_handle,"ALFA_HitCollection") );
 
 
  for (const ALFA_Hit& hit : *coll_handle) {
 
    station = hit.GetStationNumber();
    plate   = hit.GetPlateNumber();
    fiber   = hit.GetFiberNumber();
    sign    = hit.GetSignFiber();
    if (plate==100) {
      // Treat hits in GVS separately.
      m_E_gvs->Fill(station,(hit.GetEnergyDeposit()));
      continue;
    }
    E_fiber_sum[station-1][plate-1][fiber-1][(1-sign)/2] += (hit.GetEnergyDeposit());
  }
 
  for ( int l = 0; l < station_max; l++ ) {
      for ( int i = 0; i < plate_max; i++ ){
          for ( int j = 0; j < fiber_max;  j++ ) {
              for ( int k = 0; k < 2; k++ ) {
                E_full_sum[l] += E_fiber_sum[l][i][j][k];
                E_layer_sum[l][2*i+k] += E_fiber_sum[l][i][j][k];
                if (E_fiber_sum[l][i][j][k] > 0.) {
                  m_hit_map_h[l]->Fill(2*i+k+1,j+1);
                }
              }
          }
      }
  }
 
  for ( int l = 0; l < station_max; l++ ) {
    m_E_full_sum_h[l]->Fill(E_full_sum[l]);
    for ( int i = 0; i < 2*plate_max; i++ ) {
      m_E_layer_sum_h[l]->Fill(E_layer_sum[l][i],i+1);
    }
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_E_full_sum_h

TH1* ALFA_SimHitsTestTool::m_E_full_sum_h[8]

private

Definition at line 26 of file ALFA_SimHitsTestTool.h.

m_E_gvs

TH1* ALFA_SimHitsTestTool::m_E_gvs {}

private

Definition at line 27 of file ALFA_SimHitsTestTool.h.

m_E_layer_sum_h

TH2* ALFA_SimHitsTestTool::m_E_layer_sum_h[8]

private

Definition at line 25 of file ALFA_SimHitsTestTool.h.

m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

m_hit_map_h

TH2* ALFA_SimHitsTestTool::m_hit_map_h[8]

private

Definition at line 24 of file ALFA_SimHitsTestTool.h.

m_key

std::string SimTestToolBase::m_key

protectedinherited

The MC truth key.

Definition at line 34 of file SimTestToolBase.h.

m_path

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

protectedinherited

Definition at line 34 of file SimTestHisto.h.

---

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

• ALFA_SimHitsTestTool.h
• ALFA_SimHitsTestTool.cxx