G4GIDI_target Class Reference

#include <G4GIDI_target.hh>

Public Member Functions
void	init (const char *fileName)
	G4GIDI_target (const char *fileName)
	G4GIDI_target (std::string &fileName)
	~G4GIDI_target ()
std::string *	getName (void)
std::string *	getFilename (void)
int	getZ (void)
int	getA (void)
int	getM (void)
double	getMass (void)
int	getTemperatures (double *temperatures)
int	readTemperature (int index)
std::string	getEqualProbableBinSampleMethod (void)
int	setEqualProbableBinSampleMethod (std::string method)
int	getNumberOfChannels (void)
int	getNumberOfProductionChannels (void)
std::vector< channelID > *	getChannelIDs (void)
std::vector< channelID > *	getProductionChannelIDs (void)
std::vector< channelID > *	getChannelIDs2 (GIDI::tpia_channel **channels, int n)
std::vector< double > *	getEnergyGridAtTIndex (int index)
double	getTotalCrossSectionAtE (double e_in, double temperature)
double	getElasticCrossSectionAtE (double e_in, double temperature)
double	getCaptureCrossSectionAtE (double e_in, double temperature)
double	getFissionCrossSectionAtE (double e_in, double temperature)
double	getOthersCrossSectionAtE (double e_in, double temperature)
double	sumChannelCrossSectionAtE (int nIndices, int *indices, double e_in, double temperature)
int	sampleChannelCrossSectionAtE (int nIndices, int *indices, double e_in, double temperature, double(*rng)(void *), void *rngState)
double	getElasticFinalState (double e_in, double temperature, double(*rng)(void *), void *rngState)
std::vector< G4GIDI_Product > *	getCaptureFinalState (double e_in, double temperature, double(*rng)(void *), void *rngState)
std::vector< G4GIDI_Product > *	getFissionFinalState (double e_in, double temperature, double(*rng)(void *), void *rngState)
std::vector< G4GIDI_Product > *	getOthersFinalState (double e_in, double temperature, double(*rng)(void *), void *rngState)
std::vector< G4GIDI_Product > *	getFinalState (int nIndices, int *indices, double e_in, double temperature, double(*rng)(void *), void *rngState)

Data Fields
std::string	equalProbableBinSampleMethod
int	nElasticIndices
int	nCaptureIndices
int	nFissionIndices
int	nOthersIndices
int *	elasticIndices
int *	captureIndices
int *	fissionIndices
int *	othersIndices
GIDI::statusMessageReporting	smr
std::string	name
std::string	sourceFilename
double	mass
GIDI::tpia_target *	target

Detailed Description

Definition at line 88 of file G4GIDI_target.hh.

Constructor & Destructor Documentation

G4GIDI_target::G4GIDI_target

(

const char *

fileName

)

Definition at line 75 of file G4GIDI_target.cc.

                                                   {

    init( fileName );
}

G4GIDI_target::G4GIDI_target

(

std::string &

fileName

)

Definition at line 82 of file G4GIDI_target.cc.

                                               {

    init( fileName.c_str( ) );
}

G4GIDI_target::~G4GIDI_target

(

)

Definition at line 150 of file G4GIDI_target.cc.

References smr_release(), tpia_target_free(), and xData_free().

                               {

    tpia_target_free( &smr, target );
    xData_free( &smr, elasticIndices );
    smr_release( &smr );
}

Member Function Documentation

int G4GIDI_target::getA

(

void

)

Definition at line 174 of file G4GIDI_target.cc.

                              {
   
    return( target->targetID->A );
}

double G4GIDI_target::getCaptureCrossSectionAtE	(	double	e_in,
		double	temperature
	)

Definition at line 300 of file G4GIDI_target.cc.

                                                                                 {

    return( sumChannelCrossSectionAtE( nCaptureIndices, captureIndices, e_in, temperature ) );
}

vector< G4GIDI_Product > * G4GIDI_target::getCaptureFinalState	(	double	e_in,
		double	temperature,
		double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 369 of file G4GIDI_target.cc.

Referenced by G4LENDCapture::ApplyYourself().

                                                                                                                                      {

    return( getFinalState( nCaptureIndices, captureIndices, e_in, temperature, rng, rngState ) );
}

vector< channelID > * G4GIDI_target::getChannelIDs

(

void

)

Definition at line 244 of file G4GIDI_target.cc.

References tpia_target_numberOfChannels().

                                                      { 

    return( getChannelIDs2( target->baseHeatedTarget->channels, tpia_target_numberOfChannels( &smr, target ) ) );
}

vector< channelID > * G4GIDI_target::getChannelIDs2	(	GIDI::tpia_channel **	channels,
		int	n
	)

Definition at line 258 of file G4GIDI_target.cc.

References n, and tpia_channel_s::outputChannel.

                                                                                 {

    int i = 0;
    vector<channelID> *listOfChannels;

    listOfChannels = new vector<channelID>( n );
    for( i = 0; i < n; i++ ) (*listOfChannels)[i].ID = channels[i]->outputChannel;
    return( listOfChannels );
}

double G4GIDI_target::getElasticCrossSectionAtE	(	double	e_in,
		double	temperature
	)

Definition at line 293 of file G4GIDI_target.cc.

                                                                                 {

    return( sumChannelCrossSectionAtE( nElasticIndices, elasticIndices, e_in, temperature ) );
}

double G4GIDI_target::getElasticFinalState	(	double	e_in,
		double	temperature,
		double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 349 of file G4GIDI_target.cc.

References tpia_product_s::angular, tpia_channel_s::decayChannel, tpia_decaySamplingInfo_s::e_in, tpia_decaySamplingInfo_s::frame, tpia_decaySamplingInfo_s::isVelocity, tpia_decaySamplingInfo_s::mu, tpia_decaySamplingInfo_s::rng, tpia_decaySamplingInfo_s::rngState, tpia_decaySamplingInfo_s::samplingMethods, tpia_angular_SampleMu(), tpia_decayChannel_getFirstProduct(), tpia_frame_setColumn(), tpia_referenceFrame_lab, and tpia_target_heated_getChannelAtIndex_smr().

Referenced by G4LENDModel::ApplyYourself(), and G4LENDElastic::ApplyYourself().

                                                                                                          {

    tpia_decaySamplingInfo decaySamplingInfo;
    tpia_channel *channel = tpia_target_heated_getChannelAtIndex_smr( &smr, target->baseHeatedTarget, elasticIndices[0] );
    tpia_product *product;

    decaySamplingInfo.e_in = e_in;
    decaySamplingInfo.isVelocity = 0;
    tpia_frame_setColumn( &smr, &(decaySamplingInfo.frame), 0, tpia_referenceFrame_lab );
    decaySamplingInfo.samplingMethods = &(target->samplingMethods);
    decaySamplingInfo.rng = rng;
    decaySamplingInfo.rngState = rngState;
    product = tpia_decayChannel_getFirstProduct( &(channel->decayChannel) );
    tpia_angular_SampleMu( &smr, &(product->angular), &decaySamplingInfo );

    return( decaySamplingInfo.mu );
}

vector< double > * G4GIDI_target::getEnergyGridAtTIndex

(

int

index

)

Definition at line 270 of file G4GIDI_target.cc.

References n, and tpia_target_getEnergyGridAtTIndex().

                                                                {

    xData_Int i, n;
    double *dEnergyGrid;
    vector<double> *energyGrid;
    vector<double>::iterator iter;

    n = tpia_target_getEnergyGridAtTIndex( &smr, target, index, &dEnergyGrid );
    if( n < 0 ) return( NULL );
    energyGrid = new vector<double>( n );
    for( i = 0, iter = energyGrid->begin( ); i < n; i++, iter++ ) *iter = dEnergyGrid[i];
    return( energyGrid );
}

string G4GIDI_target::getEqualProbableBinSampleMethod

(

void

)

Definition at line 209 of file G4GIDI_target.cc.

                                                            {

    return( equalProbableBinSampleMethod );
}

string * G4GIDI_target::getFilename

(

void

)

Definition at line 163 of file G4GIDI_target.cc.

{ return( &sourceFilename ); }

vector< G4GIDI_Product > * G4GIDI_target::getFinalState	(	int	nIndices,
		int *	indices,
		double	e_in,
		double	temperature,
		double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 390 of file G4GIDI_target.cc.

References tpia_particle_s::A, tpia_decaySamplingInfo_s::e_in, tpia_decaySamplingInfo_s::frame, tpia_decaySamplingInfo_s::isVelocity, tpia_productOutgoingData_s::kineticEnergy, tpia_particle_s::m, n, nProductsMax, tpia_productOutgoingData_s::productID, tpia_productOutgoingData_s::px_vx, tpia_productOutgoingData_s::py_vy, tpia_productOutgoingData_s::pz_vz, tpia_decaySamplingInfo_s::rng, tpia_decaySamplingInfo_s::rngState, tpia_decaySamplingInfo_s::samplingMethods, tpia_frame_setColumn(), tpia_referenceFrame_lab, tpia_target_heated_sampleIndexChannelProductsAtE(), and tpia_particle_s::Z.

                                                                                                                                                           {

#define nProductsMax 50
    int index = 0, i, n;
    vector<G4GIDI_Product> *products = NULL;
    tpia_decaySamplingInfo decaySamplingInfo;
    tpia_productOutgoingData productDatas[nProductsMax], *productData;

    decaySamplingInfo.e_in = e_in;
    decaySamplingInfo.samplingMethods = &(target->samplingMethods);
    decaySamplingInfo.isVelocity = 0;
    tpia_frame_setColumn( &smr, &(decaySamplingInfo.frame), 0, tpia_referenceFrame_lab );
    decaySamplingInfo.rng = rng;
    decaySamplingInfo.rngState = rngState;

    if( nIndices == 0 ) {
        return( NULL ); }
    else {
        if( nIndices == 1 ) {
            index = indices[0]; }
        else {
            index = sampleChannelCrossSectionAtE( nIndices, indices, e_in, temperature, rng, rngState );
        }
    }
    n = tpia_target_heated_sampleIndexChannelProductsAtE( &smr, target->baseHeatedTarget, index, &decaySamplingInfo, nProductsMax, productDatas );
    if( n > 0 ) {
        if( ( products = new vector<G4GIDI_Product>( n ) ) != NULL ) {
            for( i = 0; i < n; i++ ) {
                productData = &(productDatas[i]);
                (*products)[i].A = productData->productID->A;
                (*products)[i].Z = productData->productID->Z;
                (*products)[i].m = productData->productID->m;
                (*products)[i].kineticEnergy = productData->kineticEnergy;
                (*products)[i].px = productData->px_vx;
                (*products)[i].py = productData->py_vy;
                (*products)[i].pz = productData->pz_vz;
            }
        }
    }

    return( products );
#undef nProductsMax
}

double G4GIDI_target::getFissionCrossSectionAtE	(	double	e_in,
		double	temperature
	)

Definition at line 307 of file G4GIDI_target.cc.

                                                                                 {

    return( sumChannelCrossSectionAtE( nFissionIndices, fissionIndices, e_in, temperature ) );
}

vector< G4GIDI_Product > * G4GIDI_target::getFissionFinalState	(	double	e_in,
		double	temperature,
		double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 376 of file G4GIDI_target.cc.

Referenced by G4LENDFission::ApplyYourself().

                                                                                                                                      {

    return( getFinalState( nFissionIndices, fissionIndices, e_in, temperature, rng, rngState ) );
}

int G4GIDI_target::getM

(

void

)

Definition at line 181 of file G4GIDI_target.cc.

                              {
   
    return( target->targetID->m );
}

double G4GIDI_target::getMass

(

void

)

Definition at line 188 of file G4GIDI_target.cc.

References ZAMass::mass.

                                    {

    return( mass );
}

string * G4GIDI_target::getName

(

void

)

Definition at line 159 of file G4GIDI_target.cc.

{ return( &name ); }

int G4GIDI_target::getNumberOfChannels

(

void

)

Definition at line 230 of file G4GIDI_target.cc.

References tpia_target_numberOfChannels().

                                             {

    return( tpia_target_numberOfChannels( &smr, target ) );
}

int G4GIDI_target::getNumberOfProductionChannels

(

void

)

Definition at line 237 of file G4GIDI_target.cc.

References tpia_target_numberOfProductionChannels().

                                                       {

    return( tpia_target_numberOfProductionChannels( &smr, target ) );
}

double G4GIDI_target::getOthersCrossSectionAtE	(	double	e_in,
		double	temperature
	)

Definition at line 314 of file G4GIDI_target.cc.

                                                                                {

    return( sumChannelCrossSectionAtE( nOthersIndices, othersIndices, e_in, temperature ) );
}

vector< G4GIDI_Product > * G4GIDI_target::getOthersFinalState	(	double	e_in,
		double	temperature,
		double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 383 of file G4GIDI_target.cc.

Referenced by G4LENDInelastic::ApplyYourself().

                                                                                                                                     {

    return( getFinalState( nOthersIndices, othersIndices, e_in, temperature, rng, rngState ) );
}

vector< channelID > * G4GIDI_target::getProductionChannelIDs

(

void

)

Definition at line 251 of file G4GIDI_target.cc.

References tpia_target_numberOfProductionChannels().

                                                                {

    return( getChannelIDs2( target->baseHeatedTarget->productionChannels, tpia_target_numberOfProductionChannels( &smr, target ) ) );
}

int G4GIDI_target::getTemperatures

(

double *

temperatures

)

Definition at line 195 of file G4GIDI_target.cc.

References tpia_target_getTemperatures().

                                                         {

    return( tpia_target_getTemperatures( &smr, target, temperatures ) );
}

double G4GIDI_target::getTotalCrossSectionAtE	(	double	e_in,
		double	temperature
	)

Definition at line 286 of file G4GIDI_target.cc.

References tpia_crossSectionType_pointwise, and tpia_target_getTotalCrossSectionAtTAndE().

                                                                               {

    return( tpia_target_getTotalCrossSectionAtTAndE( NULL, target, temperature, -1, e_in, tpia_crossSectionType_pointwise ) );
}

int G4GIDI_target::getZ

(

void

)

Definition at line 167 of file G4GIDI_target.cc.

                              {
   
    return( target->targetID->Z );
}

void G4GIDI_target::init

(

const char *

fileName

)

Definition at line 89 of file G4GIDI_target.cc.

References tpia_channel_s::ENDL_C, tpia_channel_s::fission, G4GIDI_targetMass(), ZAMass::mass, n, smr_initialize(), smr_isOk(), smr_print(), tpia_target_createRead(), tpia_target_heated_getChannelAtIndex(), tpia_target_numberOfChannels(), and xData_malloc2.

                                               {

    int i, j, n, *p, *q;
    tpia_channel *channel;

    smr_initialize( &smr );
    sourceFilename = fileName;
    target = tpia_target_createRead( &smr, fileName );
    if( !smr_isOk( &smr ) ) {
        smr_print( &smr, stderr, 1 );
        throw 1;
    }
    name = target->targetID->name;
    mass = G4GIDI_targetMass( target->targetID->name );
    equalProbableBinSampleMethod = "constant";
    elasticIndices = NULL;
    nElasticIndices = nCaptureIndices = nFissionIndices = nOthersIndices = 0;

    if( ( n = tpia_target_numberOfChannels( &smr, target ) ) > 0 ) {
        p = elasticIndices = (int *) xData_malloc2( NULL, n * sizeof( double ), 1, "elasticIndices" );
        for( i = 0; i < n; i++ ) {      /* Find elastic channel(s). */
            channel = tpia_target_heated_getChannelAtIndex( target->baseHeatedTarget, i );
            if( channel->ENDL_C == 10 ) {
                *(p++) = i;
                nElasticIndices++;
            }
        }
        captureIndices = p;
        for( i = 0; i < n; i++ ) {      /* Find capture channel(s). */
            channel = tpia_target_heated_getChannelAtIndex( target->baseHeatedTarget, i );
            if( channel->ENDL_C == 46 ) {
                *(p++) = i;
                nCaptureIndices++;
            }
        }

        fissionIndices = p;
        for( i = 0; i < n; i++ ) {      /* Find fission channel(s). */
            channel = tpia_target_heated_getChannelAtIndex( target->baseHeatedTarget, i );
            if( channel->fission != NULL ) {
                *(p++) = i;
                nFissionIndices++;
            }
        }
        othersIndices = p;
        for( i = 0; i < n; i++ ) {      /* Find other channel(s). */
            for( j = 0, q = elasticIndices; j < nElasticIndices; j++, q++ ) if( *q == i ) break;
            if( j < nElasticIndices ) continue;
            for( j = 0, q = captureIndices; j < nCaptureIndices; j++, q++ ) if( *q == i ) break;
            if( j < nCaptureIndices ) continue;
            for( j = 0, q = fissionIndices; j < nFissionIndices; j++, q++ ) if( *q == i ) break;
            if( j < nFissionIndices ) continue;
            *p = i;
            p++;
            nOthersIndices++;
        }
    }
}

int G4GIDI_target::readTemperature

(

int

index

)

Definition at line 202 of file G4GIDI_target.cc.

References tpia_target_readHeatedTarget().

                                              {

    return( tpia_target_readHeatedTarget( &smr, target, index, 0 ) );
}

int G4GIDI_target::sampleChannelCrossSectionAtE	(	int	nIndices,
		int *	indices,
		double	e_in,
		double	temperature,
		double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 333 of file G4GIDI_target.cc.

References tpia_crossSectionType_pointwise, tpia_misc_drng(), and tpia_target_getIndexChannelCrossSectionAtE().

                                                  {

    int i;
    double xsec = 0., rxsec = sumChannelCrossSectionAtE( nIndices, indices, e_in, temperature ) * tpia_misc_drng( rng, rngState );

    for( i = 0; i < nIndices - 1; i++ ) {
        xsec += tpia_target_getIndexChannelCrossSectionAtE( &smr, target, indices[i], temperature, -1, e_in, tpia_crossSectionType_pointwise );
        if( xsec >= rxsec ) break;
    }
    return( indices[i] );
}

int G4GIDI_target::setEqualProbableBinSampleMethod

(

std::string

method

)

Definition at line 216 of file G4GIDI_target.cc.

                                                                  {

    if( method == "constant" ) {
        equalProbableBinSampleMethod = "constant"; }
    if( method == "linear" ) {
        equalProbableBinSampleMethod = "linear"; }
    else {
        return( 1 );
    }
    return( 0 );
}

double G4GIDI_target::sumChannelCrossSectionAtE	(	int	nIndices,
		int *	indices,
		double	e_in,
		double	temperature
	)

Definition at line 321 of file G4GIDI_target.cc.

References tpia_crossSectionType_pointwise, and tpia_target_getIndexChannelCrossSectionAtE().

                                                                                                             {

    int i;
    double xsec = 0.;

    for( i = 0; i < nIndices; i++ ) 
        xsec += tpia_target_getIndexChannelCrossSectionAtE( &smr, target, indices[i], temperature, -1, e_in, tpia_crossSectionType_pointwise );
    return( xsec );
}

Field Documentation

int * G4GIDI_target::captureIndices

Definition at line 94 of file G4GIDI_target.hh.

int* G4GIDI_target::elasticIndices

Definition at line 94 of file G4GIDI_target.hh.

std::string G4GIDI_target::equalProbableBinSampleMethod

Definition at line 92 of file G4GIDI_target.hh.

int * G4GIDI_target::fissionIndices

Definition at line 94 of file G4GIDI_target.hh.

double G4GIDI_target::mass

Definition at line 100 of file G4GIDI_target.hh.

std::string G4GIDI_target::name

Definition at line 98 of file G4GIDI_target.hh.

int G4GIDI_target::nCaptureIndices

Definition at line 93 of file G4GIDI_target.hh.

int G4GIDI_target::nElasticIndices

Definition at line 93 of file G4GIDI_target.hh.

int G4GIDI_target::nFissionIndices

Definition at line 93 of file G4GIDI_target.hh.

int G4GIDI_target::nOthersIndices

Definition at line 93 of file G4GIDI_target.hh.

int * G4GIDI_target::othersIndices

Definition at line 94 of file G4GIDI_target.hh.

GIDI::statusMessageReporting G4GIDI_target::smr

Definition at line 97 of file G4GIDI_target.hh.

std::string G4GIDI_target::sourceFilename

Definition at line 99 of file G4GIDI_target.hh.

GIDI::tpia_target* G4GIDI_target::target

Definition at line 101 of file G4GIDI_target.hh.

---

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

• G4GIDI_target.hh
• G4GIDI_target.cc