G4GIDI_target Class Reference

#include <G4GIDI_target.hh>

Collaboration diagram for G4GIDI_target:

Public Member Functions
void	init (const char *fileName)
	G4GIDI_target (const char *fileName)
	G4GIDI_target (std::string const &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)
channelID	getChannelsID (int channelIndex)
std::vector< channelID > *	getChannelIDs (void)
std::vector< channelID > *	getProductionChannelIDs (void)
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)
double	getReactionsThreshold (int index)
double	getReactionsDomain (int index, double *EMin, double *EMax)

Public Attributes
std::string	equalProbableBinSampleMethod
int	nElasticIndices
int	nCaptureIndices
int	nFissionIndices
int	nOthersIndices
int *	elasticIndices
int *	captureIndices
int *	fissionIndices
int *	othersIndices
GIDI::statusMessageReporting	smr
int	projectilesPOPID
std::string	name
std::string	sourceFilename
double	mass
GIDI::MCGIDI_target *	target

Detailed Description

Definition at line 58 of file G4GIDI_target.hh.

Constructor & Destructor Documentation

G4GIDI_target::G4GIDI_target

(

const char *

fileName

)

Definition at line 44 of file G4GIDI_target.cc.

                                                   {

    init( fileName );
}

G4GIDI_target::G4GIDI_target

(

std::string const &

fileName

)

Definition at line 51 of file G4GIDI_target.cc.

                                                     {

    init( fileName.c_str( ) );
}

G4GIDI_target::~G4GIDI_target

(

)

Definition at line 134 of file G4GIDI_target.cc.

                               {

    MCGIDI_target_free( &smr, target );
    smr_freeMemory( (void **) &elasticIndices );
    smr_release( &smr );
}

Here is the call graph for this function:

Member Function Documentation

int G4GIDI_target::getA

(

void

)

Definition at line 158 of file G4GIDI_target.cc.

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

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

Definition at line 284 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 371 of file G4GIDI_target.cc.

                                                                                                                                      {

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

Here is the caller graph for this function:

vector< channelID > * G4GIDI_target::getChannelIDs

(

void

)

Definition at line 241 of file G4GIDI_target.cc.

                                                      { 

    int i, n = MCGIDI_target_numberOfReactions( &smr, target );
    MCGIDI_reaction *reaction;
    vector<channelID> *listOfChannels;

    listOfChannels = new vector<channelID>( n );
    for( i = 0; i < n; i++ ) {
        reaction = MCGIDI_target_heated_getReactionAtIndex( target->baseHeatedTarget, i );
        (*listOfChannels)[i] = reaction->outputChannelStr;
    }
    return( listOfChannels );
}

Here is the call graph for this function:

channelID G4GIDI_target::getChannelsID

(

int

channelIndex

)

Definition at line 228 of file G4GIDI_target.cc.

                                                         { 

    MCGIDI_reaction *reaction;

    if( ( reaction = MCGIDI_target_heated_getReactionAtIndex_smr( &smr, target->baseHeatedTarget, channelIndex ) ) == NULL ) {
        smr_print( &smr, 1 );
        throw 1;
    }
    return( string( reaction->outputChannelStr ) );    /* Only works because channelID is defined to be string. */
}

Here is the call graph for this function:

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

Definition at line 277 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 342 of file G4GIDI_target.cc.

                                                                                                                     {

    MCGIDI_decaySamplingInfo decaySamplingInfo;
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex_smr( &smr, target->baseHeatedTarget, elasticIndices[0] );
    MCGIDI_product *product;
    MCGIDI_quantitiesLookupModes mode( projectilesPOPID );

    if( ( product = MCGIDI_outputChannel_getProductAtIndex( &smr, &(reaction->outputChannel), 0 ) ) == NULL ) {
        smr_print( &smr, 1 );
        throw 1;
    }

    mode.setProjectileEnergy( e_in );
    mode.setCrossSectionMode( MCGIDI_quantityLookupMode_pointwise );
    mode.setTemperature( temperature );

    decaySamplingInfo.isVelocity = 0;
    decaySamplingInfo.rng = rng;
    decaySamplingInfo.rngState = rngState;
    if( MCGIDI_product_sampleMu( &smr, product, mode, &decaySamplingInfo ) ) {
        smr_print( &smr, 1 );
        throw 1;
    }

    return( decaySamplingInfo.mu );
}

Here is the call graph for this function:

Here is the caller graph for this function:

std::vector<double>* G4GIDI_target::getEnergyGridAtTIndex

(

int

index

)

string G4GIDI_target::getEqualProbableBinSampleMethod

(

void

)

Definition at line 193 of file G4GIDI_target.cc.

                                                            {

    return( equalProbableBinSampleMethod );
}

string * G4GIDI_target::getFilename

(

void

)

Definition at line 147 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 392 of file G4GIDI_target.cc.

                                              {

    int index = 0, i, n;
    vector<G4GIDI_Product> *products = NULL;
    MCGIDI_decaySamplingInfo decaySamplingInfo;
    MCGIDI_sampledProductsDatas sampledProductsDatas;
    MCGIDI_sampledProductsData *productData;
    MCGIDI_quantitiesLookupModes mode( projectilesPOPID );

    decaySamplingInfo.isVelocity = 0;
    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 );
        }
    }

    MCGIDI_sampledProducts_initialize( &smr, &sampledProductsDatas, 1000 );
    if( !smr_isOk( &smr ) ) {
        smr_print( &smr, 1 );
        throw 1;
    }

    mode.setProjectileEnergy( e_in );
    mode.setCrossSectionMode( MCGIDI_quantityLookupMode_pointwise );
    mode.setTemperature( temperature );

    n = MCGIDI_target_heated_sampleIndexReactionProductsAtE( &smr, target->baseHeatedTarget, index, mode,
            &decaySamplingInfo, &sampledProductsDatas );
    if( !smr_isOk( &smr ) ) {
        smr_print( &smr, 1 );
        throw 1;
    }
    if( n > 0 ) {
        if( ( products = new vector<G4GIDI_Product>( n ) ) != NULL ) {
            for( i = 0; i < n; i++ ) {
                productData = &(sampledProductsDatas.products[i]);
                (*products)[i].A = productData->pop->A;
                (*products)[i].Z = productData->pop->Z;
                (*products)[i].m = productData->pop->m;
                (*products)[i].kineticEnergy = productData->kineticEnergy;
                (*products)[i].px = productData->px_vx;
                (*products)[i].py = productData->py_vy;
                (*products)[i].pz = productData->pz_vz;
            (*products)[i].birthTimeSec = productData->birthTimeSec;
            }
        }
    }
    MCGIDI_sampledProducts_release( &smr, &sampledProductsDatas );

    return( products );
}

Here is the call graph for this function:

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

Definition at line 291 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 378 of file G4GIDI_target.cc.

                                                                                                                                      {

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

Here is the caller graph for this function:

int G4GIDI_target::getM

(

void

)

Definition at line 165 of file G4GIDI_target.cc.

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

double G4GIDI_target::getMass

(

void

)

Definition at line 172 of file G4GIDI_target.cc.

                                    {

    return( mass );
}

string * G4GIDI_target::getName

(

void

)

Definition at line 143 of file G4GIDI_target.cc.

{ return( &name ); }

int G4GIDI_target::getNumberOfChannels

(

void

)

Definition at line 214 of file G4GIDI_target.cc.

                                             {

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

Here is the call graph for this function:

int G4GIDI_target::getNumberOfProductionChannels

(

void

)

Definition at line 221 of file G4GIDI_target.cc.

                                                       {

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

Here is the call graph for this function:

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

Definition at line 298 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 385 of file G4GIDI_target.cc.

                                                                                                                                     {

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

Here is the caller graph for this function:

vector< channelID > * G4GIDI_target::getProductionChannelIDs

(

void

)

Definition at line 257 of file G4GIDI_target.cc.

                                                                {

    return( NULL );
}

double G4GIDI_target::getReactionsDomain	(	int	index,
		double *	EMin,
		double *	EMax
	)

Definition at line 461 of file G4GIDI_target.cc.

                                                                                {

    return( MCGIDI_target_heated_getReactionsDomain( &smr, target->baseHeatedTarget, index, EMin, EMax ) );
}

Here is the call graph for this function:

double G4GIDI_target::getReactionsThreshold

(

int

index

)

Definition at line 454 of file G4GIDI_target.cc.

                                                       {

    return( MCGIDI_target_heated_getReactionsThreshold( &smr, target->baseHeatedTarget, index ) );
}

Here is the call graph for this function:

int G4GIDI_target::getTemperatures

(

double *

temperatures

)

Definition at line 179 of file G4GIDI_target.cc.

                                                         {

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

Here is the call graph for this function:

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

Definition at line 264 of file G4GIDI_target.cc.

                                                                               {

    MCGIDI_quantitiesLookupModes mode( projectilesPOPID );

    mode.setProjectileEnergy( e_in );
    mode.setCrossSectionMode( MCGIDI_quantityLookupMode_pointwise );
    mode.setTemperature( temperature );

    return( MCGIDI_target_getTotalCrossSectionAtTAndE( NULL, target, mode, true ) );
}

Here is the call graph for this function:

int G4GIDI_target::getZ

(

void

)

Definition at line 151 of file G4GIDI_target.cc.

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

void G4GIDI_target::init

(

const char *

fileName

)

Definition at line 58 of file G4GIDI_target.cc.

                                               {

    int i, j, n, *p, *q, ir;
    MCGIDI_reaction *reaction;

    smr_initialize( &smr, smr_status_Ok, 1 );
    sourceFilename = fileName;
    target = MCGIDI_target_newRead( &smr, fileName );
    if( !smr_isOk( &smr ) ) {
        smr_print( &smr, 1 );
        throw 1;
    }
    projectilesPOPID = target->projectilePOP->globalPoPsIndex;
    name = target->targetPOP->name;
    mass = G4GIDI_targetMass( target->targetPOP->name );
    equalProbableBinSampleMethod = "constant";
    elasticIndices = NULL;
    nElasticIndices = nCaptureIndices = nFissionIndices = nOthersIndices = 0;

    if( ( n = MCGIDI_target_numberOfReactions( &smr, target ) ) > 0 ) {
        if( ( p = elasticIndices = (int *) smr_malloc2( &smr, n * sizeof( double ), 1, "elasticIndices" ) ) == NULL ) {
            smr_print( &smr, 1 );
            throw 1;
        }
        for( i = 0; i < n; i++ ) {      /* Find elastic channel(s). */
            reaction = MCGIDI_target_heated_getReactionAtIndex( target->baseHeatedTarget, i );
            if( MCGIDI_reaction_getENDF_MTNumber( reaction ) == 2 ) {
                *(p++) = i;
                nElasticIndices++;
            }
        }
        captureIndices = p;
        for( i = 0; i < n; i++ ) {      /* Find capture channel(s). */
            reaction = MCGIDI_target_heated_getReactionAtIndex( target->baseHeatedTarget, i );
            if( MCGIDI_reaction_getENDF_MTNumber( reaction ) == 102 ) {
                *(p++) = i;
                nCaptureIndices++;
            }
        }

        fissionIndices = p;
        for( i = 0; i < n; i++ ) {      /* Find fission channel(s). */
            reaction = MCGIDI_target_heated_getReactionAtIndex( target->baseHeatedTarget, i );
            ir = MCGIDI_reaction_getENDF_MTNumber( reaction );
            if( ( ir != 18 ) && ( ir != 19 ) && ( ir != 20 ) && ( ir != 21 ) && ( ir != 38 ) ) continue;
            *(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++;
        }
#if 0
printf( "elastic %d: ", nElasticIndices );
for( i = 0; i < nElasticIndices; i++ ) printf( " %d", elasticIndices[i] );
printf( "\ncapture %d: ", nCaptureIndices );
for( i = 0; i < nCaptureIndices; i++ ) printf( " %d", captureIndices[i] );
printf( "\nfission %d: ", nFissionIndices );
for( i = 0; i < nFissionIndices; i++ ) printf( " %d", fissionIndices[i] );
printf( "\nothers %d: ", nOthersIndices );
for( i = 0; i < nOthersIndices; i++ ) printf( " %d", othersIndices[i] );
printf( "\n" );
#endif
    }
}

Here is the call graph for this function:

int G4GIDI_target::readTemperature

(

int

index

)

Definition at line 186 of file G4GIDI_target.cc.

                                              {

    return( MCGIDI_target_readHeatedTarget( &smr, target, index ) );
}

Here is the call graph for this function:

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

Definition at line 322 of file G4GIDI_target.cc.

                                                  {

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

    mode.setProjectileEnergy( e_in );
    mode.setCrossSectionMode( MCGIDI_quantityLookupMode_pointwise );
    mode.setTemperature( temperature );

    for( i = 0; i < nIndices - 1; i++ ) {
        xsec += MCGIDI_target_getIndexReactionCrossSectionAtE( &smr, target, indices[i], mode, true );
        if( xsec >= rxsec ) break;
    }
    return( indices[i] );
}

Here is the call graph for this function:

int G4GIDI_target::setEqualProbableBinSampleMethod

(

std::string

method

)

Definition at line 200 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 305 of file G4GIDI_target.cc.

                                                                                                             {

    int i;
    double xsec = 0.;
    MCGIDI_quantitiesLookupModes mode( projectilesPOPID );

    mode.setProjectileEnergy( e_in );
    mode.setCrossSectionMode( MCGIDI_quantityLookupMode_pointwise );
    mode.setTemperature( temperature );

    for( i = 0; i < nIndices; i++ ) 
        xsec += MCGIDI_target_getIndexReactionCrossSectionAtE( &smr, target, indices[i], mode, true );
    return( xsec );
}

Here is the call graph for this function:

Member Data Documentation

int * G4GIDI_target::captureIndices

Definition at line 64 of file G4GIDI_target.hh.

int* G4GIDI_target::elasticIndices

Definition at line 64 of file G4GIDI_target.hh.

std::string G4GIDI_target::equalProbableBinSampleMethod

Definition at line 62 of file G4GIDI_target.hh.

int * G4GIDI_target::fissionIndices

Definition at line 64 of file G4GIDI_target.hh.

double G4GIDI_target::mass

Definition at line 71 of file G4GIDI_target.hh.

std::string G4GIDI_target::name

Definition at line 69 of file G4GIDI_target.hh.

int G4GIDI_target::nCaptureIndices

Definition at line 63 of file G4GIDI_target.hh.

int G4GIDI_target::nElasticIndices

Definition at line 63 of file G4GIDI_target.hh.

int G4GIDI_target::nFissionIndices

Definition at line 63 of file G4GIDI_target.hh.

int G4GIDI_target::nOthersIndices

Definition at line 63 of file G4GIDI_target.hh.

int * G4GIDI_target::othersIndices

Definition at line 64 of file G4GIDI_target.hh.

int G4GIDI_target::projectilesPOPID

Definition at line 68 of file G4GIDI_target.hh.

GIDI::statusMessageReporting G4GIDI_target::smr

Definition at line 67 of file G4GIDI_target.hh.

std::string G4GIDI_target::sourceFilename

Definition at line 70 of file G4GIDI_target.hh.

GIDI::MCGIDI_target* G4GIDI_target::target

Definition at line 72 of file G4GIDI_target.hh.

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/lend/include/G4GIDI_target.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/lend/src/G4GIDI_target.cc