ISF::FastCaloSimV2ParamSvc Class Reference

#include <FastCaloSimV2ParamSvc.h>

Inheritance diagram for ISF::FastCaloSimV2ParamSvc:

Collaboration diagram for ISF::FastCaloSimV2ParamSvc:

Public Member Functions
	FastCaloSimV2ParamSvc (const std::string &name, ISvcLocator *pSvcLocator)
	Constructor with parameters. More...
virtual	~FastCaloSimV2ParamSvc ()=default
	Destructor. More...
virtual StatusCode	initialize () override final
	Athena algorithm's interface methods. More...
virtual StatusCode	simulate (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const override final
	Simulation Call. More...

Private Member Functions
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...

Private Attributes
ServiceHandle< StoreGateSvc >	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::string	m_paramsFilename {"TFCSparam.root"}
std::string	m_paramsObject {"SelPDGID"}
TFCSParametrizationBase *	m_param {}
std::unique_ptr< CaloGeometryFromCaloDDM >	m_caloGeo {}
bool	m_printParametrization {false}
bool	m_CompressMemory {true}
bool	m_runOnGPU {false}

Detailed Description

Definition at line 22 of file FastCaloSimV2ParamSvc.h.

Constructor & Destructor Documentation

FastCaloSimV2ParamSvc()

ISF::FastCaloSimV2ParamSvc::FastCaloSimV2ParamSvc	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Constructor with parameters.

Constructor.

Definition at line 44 of file FastCaloSimV2ParamSvc.cxx.

  : base_class(name, svc)
  , m_detStore( "StoreGateSvc/DetectorStore", name )
 
{
  declareProperty( "DetStore",
                   m_detStore,
                   "Handle to a StoreGateSvc/DetectorStore instance: it will be used to "
                   "retrieve data during the course of the job" );
  declareProperty("ParamsInputFilename"            ,       m_paramsFilename);
  declareProperty("ParamsInputObject"              ,       m_paramsObject);
  declareProperty("PrintParametrization"           ,       m_printParametrization);
  declareProperty("CompressMemory"                 ,       m_CompressMemory);
  declareProperty("RunOnGPU"                       ,       m_runOnGPU);
}

~FastCaloSimV2ParamSvc()

virtual ISF::FastCaloSimV2ParamSvc::~FastCaloSimV2ParamSvc ( )

virtualdefault

Destructor.

Member Function Documentation

detStore()

const ServiceHandle< StoreGateSvc > & ISF::FastCaloSimV2ParamSvc::detStore ( ) const

inlineprivate

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

Definition at line 74 of file FastCaloSimV2ParamSvc.h.

{ return m_detStore; }

initialize()

StatusCode ISF::FastCaloSimV2ParamSvc::initialize ( )

finaloverridevirtual

Athena algorithm's interface methods.

framework methods

does not take ownership

Definition at line 61 of file FastCaloSimV2ParamSvc.cxx.

{
  const FCALDetectorManager *fcalManager{};
  ATH_CHECK(detStore()->retrieve(fcalManager));
 
  m_caloGeo = std::make_unique<CaloGeometryFromCaloDDM>();
 
  const CaloDetDescrManager* caloMgr = detStore()->tryConstRetrieve<CaloDetDescrManager>(caloMgrStaticKey);
  if(caloMgr) {
    m_caloGeo->LoadGeometryFromCaloDDM(caloMgr);
  }
  else {
    std::unique_ptr<CaloDetDescrManager> caloMgrPtr = buildCaloDetDescrNoAlign(serviceLocator()
                                           , Athena::getMessageSvc());
    m_caloGeo->LoadGeometryFromCaloDDM(caloMgrPtr.get());
    ATH_CHECK(detStore()->record(std::move(caloMgrPtr), caloMgrStaticKey));
  }
 
  if (!m_caloGeo->LoadFCalChannelMapFromFCalDDM(fcalManager)) {
    ATH_MSG_FATAL("Found inconsistency between FCal_Channel map and GEO file. Please, check if they are configured properly.");
    return StatusCode::FAILURE;
  }
 
  const std::string fileName = m_paramsFilename; //?
  std::string inputFile=PathResolverFindCalibFile(fileName);
  std::unique_ptr<TFile> paramsFile(TFile::Open( inputFile.c_str(), "READ" ));
  if (paramsFile == nullptr) {
    ATH_MSG_ERROR("file = "<<m_paramsFilename<< " not found");
    return StatusCode::FAILURE;
  }
  ATH_MSG_INFO("Opened parametrization file = "<<m_paramsFilename);
  paramsFile->ls();
  m_param=static_cast<TFCSParametrizationBase*>(paramsFile->Get(m_paramsObject.c_str()));
  if (!m_param) {
    ATH_MSG_WARNING("file = "<<m_paramsFilename<< ", object "<< m_paramsObject<<" not found");
    return StatusCode::FAILURE;
  }
 
  paramsFile->Close();
 
  if(m_CompressMemory) m_param->RemoveDuplicates();
  m_param->set_geometry(m_caloGeo.get()); 
  m_param->setLevel(msg().level());
  if (m_printParametrization) {
    m_param->Print("short");
  }
  if(m_CompressMemory) m_param->RemoveNameTitle();
 
#ifdef USE_GPU
 
  // if no GPU is found then run on CPU
  if( !m_runOnGPU )
    return StatusCode::SUCCESS;
  if( AthCUDA::Info::instance().nDevices() == 0 ) {
    ATH_MSG_WARNING( "no GPU device is found, will run everything on the CPU." );
    m_runOnGPU=false;
    return StatusCode::SUCCESS;
  }
 
  // FCS-GPU
  // copy the geometry and parameterization files to GPU
  ATH_MSG_INFO("start copying all parametrization files to GPU");
  m_param->Copy2GPU();
  //random seed for curand in GPU
  unsigned long long seed=23424352;
  m_rd4h = CaloGpuGeneral::Rand4Hits_init( MAXHITS, MAXBINS, seed, true );
  //construct and load Geometry in GPU
  ATH_MSG_INFO("start loading Geometry in GPU");
  m_gl = new GeoLoadGpu();
  CaloGeometry* calogeo=m_caloGeo.get();
  m_gl->set_ncells( calogeo->get_cells()->size() );
 
  if(calogeo->get_cells()->size() ==0 ) return StatusCode::SUCCESS;
 
  m_gl->set_max_sample( CaloGeometry::MAX_SAMPLING);
  int nrgns= calogeo->get_tot_regions();
  m_gl->set_nregions( nrgns);
 
 
  if(!convert_cellmap(calogeo->get_cells(),&m_cellmap_gpu)){
    ATH_MSG_ERROR("can not convert the CaloDetDescrElement to CaloDetDescrElement_GPU ");
  }
  m_gl->set_cellmap( &m_cellmap_gpu);
  GeoRegion* GR_ptr = (GeoRegion*)malloc( nrgns *sizeof( GeoRegion ));
  m_gl->set_regions( GR_ptr );
 
  Rg_Sample_Index* si = (Rg_Sample_Index*)malloc(CaloGeometry::MAX_SAMPLING * sizeof( Rg_Sample_Index ) )    ;
 
  m_gl->set_sample_index_h( si );
  //set regions
  int i = 0;
  for ( int is = 0; is < CaloGeometry::MAX_SAMPLING; ++is ) {
    si[is].index = i;
    int nr       = calogeo->get_n_regions( is );
    si[is].size  = nr;
    for ( int ir = 0; ir < nr; ++ir ) region_data_cpy( calogeo->get_region( is, ir), &GR_ptr[i++] );
  }
 
  if(m_gl->LoadGpu())
    ATH_MSG_INFO("calo geometry loaded successfully in GPU");
  else
    ATH_MSG_ERROR("can not load calo geoemtry in GPU");
#endif
 
  return StatusCode::SUCCESS;
}

simulate()

StatusCode ISF::FastCaloSimV2ParamSvc::simulate ( TFCSSimulationState &  simulstate, const TFCSTruthState *  truth, const TFCSExtrapolationState *  extrapol  ) const

finaloverridevirtual

Simulation Call.

Definition at line 169 of file FastCaloSimV2ParamSvc.cxx.

{
#ifdef USE_GPU
  if( m_runOnGPU ) {
    simulstate.set_gpu_rand (m_rd4h );
    simulstate.set_geold( m_gl );
  }
  else
    simulstate.set_geold( nullptr );
#endif
 
  FCSReturnCode status = m_param->simulate(simulstate, truth, extrapol);
  if (status != FCSSuccess) {
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_caloGeo

std::unique_ptr<CaloGeometryFromCaloDDM> ISF::FastCaloSimV2ParamSvc::m_caloGeo {}

private

Definition at line 61 of file FastCaloSimV2ParamSvc.h.

m_CompressMemory

bool ISF::FastCaloSimV2ParamSvc::m_CompressMemory {true}

private

Definition at line 64 of file FastCaloSimV2ParamSvc.h.

m_detStore

ServiceHandle<StoreGateSvc> ISF::FastCaloSimV2ParamSvc::m_detStore

private

Pointer to StoreGate (detector store by default)

Definition at line 54 of file FastCaloSimV2ParamSvc.h.

m_param

TFCSParametrizationBase* ISF::FastCaloSimV2ParamSvc::m_param {}

private

Definition at line 59 of file FastCaloSimV2ParamSvc.h.

m_paramsFilename

std::string ISF::FastCaloSimV2ParamSvc::m_paramsFilename {"TFCSparam.root"}

private

Definition at line 56 of file FastCaloSimV2ParamSvc.h.

m_paramsObject

std::string ISF::FastCaloSimV2ParamSvc::m_paramsObject {"SelPDGID"}

private

Definition at line 57 of file FastCaloSimV2ParamSvc.h.

m_printParametrization

bool ISF::FastCaloSimV2ParamSvc::m_printParametrization {false}

private

Definition at line 63 of file FastCaloSimV2ParamSvc.h.

m_runOnGPU

bool ISF::FastCaloSimV2ParamSvc::m_runOnGPU {false}

private

Definition at line 65 of file FastCaloSimV2ParamSvc.h.

---

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

• FastCaloSimV2ParamSvc.h
• FastCaloSimV2ParamSvc.cxx