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)

Public Attributes
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() [1/2]

G4GIDI_target::G4GIDI_target

(

const char *

fileName

)

Definition at line 75 of file G4GIDI_target.cc.

                                                   {
 
    init( fileName );
}

G4GIDI_target() [2/2]

G4GIDI_target::G4GIDI_target

(

std::string &

fileName

)

Definition at line 82 of file G4GIDI_target.cc.

                                               {
 
    init( fileName.c_str( ) );
}

~G4GIDI_target()

G4GIDI_target::~G4GIDI_target

(

)

Definition at line 150 of file G4GIDI_target.cc.

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

Member Function Documentation

getA()

int G4GIDI_target::getA

(

void

)

Definition at line 174 of file G4GIDI_target.cc.

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

getCaptureCrossSectionAtE()

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 ) );
}

getCaptureFinalState()

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.

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

Referenced by G4LENDCapture::ApplyYourself().

getChannelIDs()

vector< channelID > * G4GIDI_target::getChannelIDs

(

void

)

Definition at line 244 of file G4GIDI_target.cc.

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

getChannelIDs2()

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

Definition at line 258 of file G4GIDI_target.cc.

                                                                                 {
 
    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 );
}

Referenced by getChannelIDs(), and getProductionChannelIDs().

getElasticCrossSectionAtE()

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 ) );
}

getElasticFinalState()

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

Definition at line 349 of file G4GIDI_target.cc.

                                                                                                          {
 
    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 );
}

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

getEnergyGridAtTIndex()

vector< double > * G4GIDI_target::getEnergyGridAtTIndex

(

int

index

)

Definition at line 270 of file G4GIDI_target.cc.

                                                                {
 
    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 );
}

getEqualProbableBinSampleMethod()

string G4GIDI_target::getEqualProbableBinSampleMethod

(

void

)

Definition at line 209 of file G4GIDI_target.cc.

                                                            {
 
    return( equalProbableBinSampleMethod );
}

getFilename()

string * G4GIDI_target::getFilename

(

void

)

Definition at line 163 of file G4GIDI_target.cc.

{ return( &sourceFilename ); }

getFinalState()

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.

                                                                                                                                                           {
 
#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
}

Referenced by getCaptureFinalState(), getFissionFinalState(), and getOthersFinalState().

getFissionCrossSectionAtE()

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 ) );
}

getFissionFinalState()

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.

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

Referenced by G4LENDFission::ApplyYourself().

getM()

int G4GIDI_target::getM

(

void

)

Definition at line 181 of file G4GIDI_target.cc.

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

getMass()

double G4GIDI_target::getMass

(

void

)

Definition at line 188 of file G4GIDI_target.cc.

                                    {
 
    return( mass );
}

getName()

string * G4GIDI_target::getName

(

void

)

Definition at line 159 of file G4GIDI_target.cc.

{ return( &name ); }

getNumberOfChannels()

int G4GIDI_target::getNumberOfChannels

(

void

)

Definition at line 230 of file G4GIDI_target.cc.

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

getNumberOfProductionChannels()

int G4GIDI_target::getNumberOfProductionChannels

(

void

)

Definition at line 237 of file G4GIDI_target.cc.

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

getOthersCrossSectionAtE()

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 ) );
}

getOthersFinalState()

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.

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

Referenced by G4LENDInelastic::ApplyYourself().

getProductionChannelIDs()

vector< channelID > * G4GIDI_target::getProductionChannelIDs

(

void

)

Definition at line 251 of file G4GIDI_target.cc.

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

getTemperatures()

int G4GIDI_target::getTemperatures

(

double *

temperatures

)

Definition at line 195 of file G4GIDI_target.cc.

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

getTotalCrossSectionAtE()

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

Definition at line 286 of file G4GIDI_target.cc.

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

getZ()

int G4GIDI_target::getZ

(

void

)

Definition at line 167 of file G4GIDI_target.cc.

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

init()

void G4GIDI_target::init

(

const char *

fileName

)

Definition at line 89 of file G4GIDI_target.cc.

                                               {
 
    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++;
        }
    }
}

Referenced by G4GIDI_target().

readTemperature()

int G4GIDI_target::readTemperature

(

int

index

)

Definition at line 202 of file G4GIDI_target.cc.

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

sampleChannelCrossSectionAtE()

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.

                                                  {
 
    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] );
}

Referenced by getFinalState().

setEqualProbableBinSampleMethod()

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 );
}

sumChannelCrossSectionAtE()

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

Definition at line 321 of file G4GIDI_target.cc.

                                                                                                             {
 
    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 );
}

Referenced by getCaptureCrossSectionAtE(), getElasticCrossSectionAtE(), getFissionCrossSectionAtE(), getOthersCrossSectionAtE(), and sampleChannelCrossSectionAtE().

Member Data Documentation

captureIndices

int * G4GIDI_target::captureIndices

Definition at line 94 of file G4GIDI_target.hh.

Referenced by getCaptureCrossSectionAtE(), getCaptureFinalState(), and init().

elasticIndices

int* G4GIDI_target::elasticIndices

Definition at line 94 of file G4GIDI_target.hh.

Referenced by getElasticCrossSectionAtE(), getElasticFinalState(), init(), and ~G4GIDI_target().

equalProbableBinSampleMethod

std::string G4GIDI_target::equalProbableBinSampleMethod

Definition at line 92 of file G4GIDI_target.hh.

Referenced by getEqualProbableBinSampleMethod(), init(), and setEqualProbableBinSampleMethod().

fissionIndices

int * G4GIDI_target::fissionIndices

Definition at line 94 of file G4GIDI_target.hh.

Referenced by getFissionCrossSectionAtE(), getFissionFinalState(), and init().

mass

double G4GIDI_target::mass

Definition at line 100 of file G4GIDI_target.hh.

Referenced by getMass(), and init().

name

std::string G4GIDI_target::name

Definition at line 98 of file G4GIDI_target.hh.

Referenced by getName(), and init().

nCaptureIndices

int G4GIDI_target::nCaptureIndices

Definition at line 93 of file G4GIDI_target.hh.

Referenced by getCaptureCrossSectionAtE(), getCaptureFinalState(), and init().

nElasticIndices

int G4GIDI_target::nElasticIndices

Definition at line 93 of file G4GIDI_target.hh.

Referenced by getElasticCrossSectionAtE(), and init().

nFissionIndices

int G4GIDI_target::nFissionIndices

Definition at line 93 of file G4GIDI_target.hh.

Referenced by getFissionCrossSectionAtE(), getFissionFinalState(), and init().

nOthersIndices

int G4GIDI_target::nOthersIndices

Definition at line 93 of file G4GIDI_target.hh.

Referenced by getOthersCrossSectionAtE(), getOthersFinalState(), and init().

othersIndices

int * G4GIDI_target::othersIndices

Definition at line 94 of file G4GIDI_target.hh.

Referenced by getOthersCrossSectionAtE(), getOthersFinalState(), and init().

smr

GIDI::statusMessageReporting G4GIDI_target::smr

Definition at line 97 of file G4GIDI_target.hh.

Referenced by getChannelIDs(), getElasticFinalState(), getEnergyGridAtTIndex(), getFinalState(), getNumberOfChannels(), getNumberOfProductionChannels(), getProductionChannelIDs(), getTemperatures(), init(), readTemperature(), sampleChannelCrossSectionAtE(), sumChannelCrossSectionAtE(), and ~G4GIDI_target().

sourceFilename

std::string G4GIDI_target::sourceFilename

Definition at line 99 of file G4GIDI_target.hh.

Referenced by getFilename(), and init().

target

GIDI::tpia_target* G4GIDI_target::target

Definition at line 101 of file G4GIDI_target.hh.

Referenced by getA(), getChannelIDs(), getElasticFinalState(), getEnergyGridAtTIndex(), getFinalState(), getM(), getNumberOfChannels(), getNumberOfProductionChannels(), getProductionChannelIDs(), getTemperatures(), getTotalCrossSectionAtE(), getZ(), init(), readTemperature(), sampleChannelCrossSectionAtE(), sumChannelCrossSectionAtE(), and ~G4GIDI_target().

---

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

• G4GIDI_target.hh
• G4GIDI_target.cc