Photospp::Photos Class Reference

#include <Photos.h>

Public Member Functions
Photos &	getInstance ()
Photos &	getInstance ()

Static Public Member Functions
static void	initialize ()
static void	iniInfo ()
static void	processParticle (PhotosParticle *p)
static void	processBranch (PhotosParticle *p)
static void	suppressBremForDecay (int count, int motherID,...)
static void	suppressBremForBranch (int count, int motherID,...)
static void	suppressAll ()
static void	forceBremForDecay (int count, int motherID,...)
static void	forceBremForBranch (int count, int motherID,...)
static void	createHistoryEntries (bool flag, int status)
static void	ignoreParticlesOfStatus (int status)
static void	deIgnoreParticlesOfStatus (int status)
static bool	isStatusCodeIgnored (int status)
static void	setRandomGenerator (double(*gen)())
static void	setSeed (int iseed1, int iseed2)
static void	maxWtInterference (double interference)
static void	setInfraredCutOff (double cut_off)
static void	setAlphaQED (double alpha)
static void	setInterference (bool interference)
static void	setDoubleBrem (bool doub)
static void	setQuatroBrem (bool quatroBrem)
static void	setCorrectionWtForW (bool corr)
static void	setExponentiation (bool expo)
static void	setMeCorrectionWtForScalar (bool corr)
static void	setMeCorrectionWtForW (bool corr)
static void	setMeCorrectionWtForZ (bool corr)
static void	setTopProcessRadiation (bool top)
static void	setStopAtCriticalError (bool stop)
static void	initializeKinematicCorrections (int flag)
static void	forceMassFrom4Vector (bool flag)
static void	forceMassFromEventRecord (int pdgid)
static void	forceMass (int pdgid, double mass)
static void	setMomentumConservationThreshold (double threshold)
static void	initialize ()
static void	iniInfo ()
static void	processParticle (PhotosParticle *p)
static void	processBranch (PhotosParticle *p)
static void	suppressBremForDecay (int count, int motherID,...)
static void	suppressBremForBranch (int count, int motherID,...)
static void	suppressAll ()
static void	forceBremForDecay (int count, int motherID,...)
static void	forceBremForBranch (int count, int motherID,...)
static void	createHistoryEntries (bool flag, int status)
static void	ignoreParticlesOfStatus (int status)
static void	deIgnoreParticlesOfStatus (int status)
static bool	isStatusCodeIgnored (int status)
static void	setRandomGenerator (double(*gen)())
static void	setSeed (int iseed1, int iseed2)
static void	maxWtInterference (double interference)
static void	setInfraredCutOff (double cut_off)
static void	setAlphaQED (double alpha)
static void	setInterference (bool interference)
static void	setDoubleBrem (bool doub)
static void	setQuatroBrem (bool quatroBrem)
static void	setCorrectionWtForW (bool corr)
static void	setExponentiation (bool expo)
static void	setMeCorrectionWtForScalar (bool corr)
static void	setMeCorrectionWtForW (bool corr)
static void	setMeCorrectionWtForZ (bool corr)
static void	setTopProcessRadiation (bool top)
static void	setStopAtCriticalError (bool stop)
static void	initializeKinematicCorrections (int flag)
static void	forceMassFrom4Vector (bool flag)
static void	forceMassFromEventRecord (int pdgid)
static void	forceMass (int pdgid, double mass)
static void	setMomentumConservationThreshold (double threshold)

Static Public Attributes
static const int	VER_MAJOR =3
static const int	VER_MINOR =52
static const int	DAT_DAY =4
static const int	DAT_MONTH =11
static const int	DAT_YEAR =12
static bool	isSuppressed =false
static bool	massFrom4Vector =true
static vector< vector< int > * > *	supBremList = 0
static vector< vector< int > * > *	forceBremList = 0
static vector< pair< int, double > * > *	forceMassList = 0
static vector< int > *	ignoreStatusCodeList = 0
static double	momentum_conservation_threshold = 0.1
static bool	meCorrectionWtForScalar =false
static bool	meCorrectionWtForZ =false
static bool	meCorrectionWtForW =false
static bool	isCreateHistoryEntries =false
static int	historyEntriesStatus = 3
static double(*	randomDouble )() = PhotosRandom::randomReal

Private Member Functions
	Photos ()
	~Photos ()
	Photos (const Photos &)
Photos &	operator= (const Photos &)
	Photos ()
	~Photos ()
	Photos (const Photos &)
Photos &	operator= (const Photos &)

Static Private Attributes
static Photos	_instance

Detailed Description

Definition at line 28 of file Photos.h.

Constructor & Destructor Documentation

Photos::Photos ( )

private

Definition at line 33 of file Photos.cxx.

References phokey_, setAlphaQED(), setCorrectionWtForW(), setDoubleBrem(), setInfraredCutOff(), setInterference(), setQuatroBrem(), and setTopProcessRadiation().

{
        setAlphaQED           (0.00729735039);
        setInfraredCutOff     (0.01);
        setInterference       (true);
        setDoubleBrem         (true);
        setQuatroBrem         (false);
        setTopProcessRadiation(true);
        setCorrectionWtForW   (true);

        // setExponentiation(true) moved to initialize() due to communication
        // problems with Fortran under MacOS.
        phokey_.iexp = 1;
}

Photospp::Photos::~Photos ( )

inlineprivate

Definition at line 198 of file Photos.h.

{}

Photospp::Photos::Photos ( const Photos &  )

private

Photospp::Photos::Photos ( )

private

Photospp::Photos::~Photos ( )

inlineprivate

Definition at line 198 of file 3.52/src/photosCInterfaces/Photos.h.

{}

Photospp::Photos::Photos ( const Photos &  )

private

Member Function Documentation

void Photos::createHistoryEntries ( bool  flag, int  status  )

static

If event record allows it, create history entries of particles before Photos processing

Definition at line 289 of file Photos.cxx.

References historyEntriesStatus, ignoreParticlesOfStatus(), isCreateHistoryEntries, status, and Photospp::Log::Warning().

{
  if(status<3)
  {
    Log::Warning()<<"Photos::createHistoryEntries: status must be >=3"<<endl;
    return;
  }

  isCreateHistoryEntries = flag;
  historyEntriesStatus   = status;
  ignoreParticlesOfStatus(status);
}

static void Photospp::Photos::createHistoryEntries ( bool  flag, int  status  )

static

If event record allows it, create history entries of particles before Photos processing

static void Photospp::Photos::deIgnoreParticlesOfStatus ( int  status )

static

Remove 'status' from the list of ignored status codes

void Photos::deIgnoreParticlesOfStatus ( int  status )

static

Remove 'status' from the list of ignored status codes

Definition at line 319 of file Photos.cxx.

References i, and ignoreStatusCodeList.

{
  if(!ignoreStatusCodeList) return;

  for(unsigned int i=0;i<ignoreStatusCodeList->size();i++)
  {
    if( status==ignoreStatusCodeList->at(i) )
    {
      ignoreStatusCodeList->erase(ignoreStatusCodeList->begin()+i);
      return;
    }
  }
}

void Photos::forceBremForBranch ( int  count, int  motherID,   ...  )

static

Force processing of a whole decay branch

Definition at line 273 of file Photos.cxx.

References count, forceBremList, i, and v.

{
        va_list arg;
        va_start(arg, motherID);
        vector<int> *v = new vector<int>();
        v->push_back(motherID);
        for(int i = 0;i<count;i++)
        {
                v->push_back(va_arg(arg,int));
        }
        va_end(arg);
        v->push_back(1);
        if(!forceBremList) forceBremList = new vector< vector<int>* >();
        forceBremList->push_back(v);
}

static void Photospp::Photos::forceBremForBranch ( int  count, int  motherID,   ...  )

static

Force processing of a whole decay branch

static void Photospp::Photos::forceBremForDecay ( int  count, int  motherID,   ...  )

static

Force processing of a single decay

void Photos::forceBremForDecay ( int  count, int  motherID,   ...  )

static

Force processing of a single decay

Definition at line 257 of file Photos.cxx.

References count, forceBremList, i, and v.

{
        va_list arg;
        va_start(arg, motherID);
        vector<int> *v = new vector<int>();
        v->push_back(motherID);
        for(int i = 0;i<count;i++)
        {
                v->push_back(va_arg(arg,int));
        }
        va_end(arg);
        v->push_back(0);
        if(!forceBremList) forceBremList = new vector< vector<int>* >();
        forceBremList->push_back(v);
}

void Photos::forceMass ( int  pdgid, double  mass  )

static

When particles with PDGID and -PDGID will be processed by Photos, their mass value will be given by user instead of being calculated from 4-vector.

This works only if 'forceMassFrom4Vector' is set to 'true' (default) This routine may be executed several times with different PDGID values.

Definition at line 394 of file Photos.cxx.

References forceMassList, and Photospp::Log::Warning().

{
  if(mass<0.0)
  {
    Log::Warning()<<"Photos::forceMass: Mass must be > 0.0"<<endl;
    return;
  }
  
  if(!forceMassList) forceMassList = new vector<pair<int,double>* >();
  forceMassList->push_back( new pair<int,double>(pdgid, mass) );
}

static void Photospp::Photos::forceMass ( int  pdgid, double  mass  )

static

When particles with PDGID and -PDGID will be processed by Photos, their mass value will be given by user instead of being calculated from 4-vector.

This works only if 'forceMassFrom4Vector' is set to 'true' (default) This routine may be executed several times with different PDGID values.

static void Photospp::Photos::forceMassFrom4Vector ( bool  flag )

inlinestatic

Force mass value to be sqrt(e^2-p^2) for all particle momenta taken from event record. May be important for numerical stability. May lead to faulty results due to rounding errors for hiper-relativistic electron, for example.

Definition at line 131 of file 3.52/src/photosCInterfaces/Photos.h.

References massFrom4Vector.

{ massFrom4Vector=flag; }

static void Photospp::Photos::forceMassFrom4Vector ( bool  flag )

inlinestatic

Force mass value to be sqrt(e^2-p^2) for all particle momenta taken from event record. May be important for numerical stability. May lead to faulty results due to rounding errors for hiper-relativistic electron, for example.

Definition at line 131 of file Photos.h.

References massFrom4Vector.

{ massFrom4Vector=flag; }

static void Photospp::Photos::forceMassFromEventRecord ( int  pdgid )

static

When particles with PDGID and -PDGID will be processed by Photos, their mass value will be taken from event record instead of being calculated from 4-vector.

This works only if 'forceMassFrom4Vector' is set to 'true' (default) This routine may be executed several times with different PDGID values.

void Photos::forceMassFromEventRecord ( int  pdgid )

static

When particles with PDGID and -PDGID will be processed by Photos, their mass value will be taken from event record instead of being calculated from 4-vector.

This works only if 'forceMassFrom4Vector' is set to 'true' (default) This routine may be executed several times with different PDGID values.

Definition at line 388 of file Photos.cxx.

References forceMassList.

{
  if(!forceMassList) forceMassList = new vector<pair<int,double>* >();
  forceMassList->push_back( new pair<int,double>(pdgid, -1.0) );
}

Photos& Photospp::Photos::getInstance ( )

inline

Get instance of Photos

Definition at line 193 of file 3.52/src/photosCInterfaces/Photos.h.

References _instance.

{ return _instance; }

Photos& Photospp::Photos::getInstance ( )

inline

Get instance of Photos

Definition at line 193 of file Photos.h.

References _instance.

{ return _instance; }

void Photos::ignoreParticlesOfStatus ( int  status )

static

Ignore particles with given status code

Definition at line 302 of file Photos.cxx.

References i, ignoreStatusCodeList, and Photospp::Log::Warning().

Referenced by createHistoryEntries().

{
  if(status<3)
  {
    Log::Warning()<<"Photos::ignoreParticlesOfStatus: status must be >=3"<<endl;
    return;
  }
  
  if(!ignoreStatusCodeList) ignoreStatusCodeList = new vector<int>();

  // Do not add duplicate entries to the list
  for(unsigned int i=0;i<ignoreStatusCodeList->size();i++)
    if( status==ignoreStatusCodeList->at(i) ) return;
  
  ignoreStatusCodeList->push_back(status);
}

static void Photospp::Photos::ignoreParticlesOfStatus ( int  status )

static

Ignore particles with given status code

void Photos::iniInfo ( )

static

Prints info on Photos initialization (reinitialization) status

Definition at line 170 of file Photos.cxx.

References meCorrectionWtForScalar, meCorrectionWtForW, meCorrectionWtForZ, phocop_, and phokey_.

{
// prints infomation like Photos::initialize; may be called after reinitializations.

/*******************************************************************************
  Once  parameter setters are called after PHOINI one may want to know and write
  into output info including all reinitializations.
*******************************************************************************/


//------------------------------------------------------------------------------
// Print PHOTOS header again
//------------------------------------------------------------------------------
        int                coutPrec = cout.precision(6);
        ios_base::fmtflags flags    = cout.setf(ios_base::floatfield);
        cout<<endl;
        cout<<"********************************************************************************"<<endl<<endl;
        cout<<"                            ========================================="<<endl;
        cout<<"                            PHOTOS, information routine"<<endl;
        cout<<"                            Input parameters after reinitialization: "<<endl<<endl;
        cout<<"                            ========================================="<<endl<<endl;
        cout<<"********************************************************************************"<<endl<<endl;
        cout<<"                    INTERF= "<<phokey_.interf<<" ISEC= " <<phokey_.isec <<" ITRE= "<<phokey_.itre
                             <<" IEXP= "  <<phokey_.iexp  <<" IFTOP= "<<phokey_.iftop<<" IFW= " <<phokey_.ifw <<endl;
        cout<<"                    ALPHA_QED= "<<phocop_.alpha<<" XPHCUT= "<<phocop_.xphcut<<endl<<endl;

        if(phokey_.interf) cout<<"                    Option with interference is active"<<endl;
        if(phokey_.isec)   cout<<"                    Option with double photons is active"<<endl;
        if(phokey_.itre)   cout<<"                    Option with triple/quatric photons is active"<<endl;
        if(phokey_.iexp)   cout<<"                    Option with exponentiation is active EPSEXP="<<phokey_.expeps<<endl;
        if(phokey_.iftop)  cout<<"                    Emision in t tbar production is active"<<endl;
        if(phokey_.ifw)    cout<<"                    Correction wt in decay of W is active"<<endl;
        if(meCorrectionWtForZ)    cout<<"                    ME in decay of Z is active"<<endl;
        if(meCorrectionWtForW)    cout<<"                    ME in decay of W is active"<<endl;
        if(meCorrectionWtForScalar)    cout<<"                    ME in decay of Scalar is active"<<endl;

        cout<<endl<<"          WARNING:  /HEPEVT/ is not anymore used."<<endl<<endl;
        // Revert output stream flags and precision
        cout.precision(coutPrec);
        cout.flags    (flags);
}

static void Photospp::Photos::iniInfo ( )

static

Prints info on Photos initialization (reinitialization) status

void Photos::initialize ( )

static

Initalize Photos with the parameters previously set via the setter methods

Definition at line 48 of file Photos.cxx.

References DAT_DAY, DAT_MONTH, DAT_YEAR, i, Photospp::PhotosRandom::initialize(), initializeKinematicCorrections(), iphekl_(), iphqrk_(), maxWtInterference(), meCorrectionWtForW, meCorrectionWtForZ, phocop_, phokey_, pholun_, phosta_, phpico_, setExponentiation(), suppressBremForDecay(), VER_MAJOR, and VER_MINOR.

{
// Should return if already initialized?

/*******************************************************************************
  All the following parameter setters can be called after PHOINI.
  Initialization of kinematic correction against rounding errors.
  The set values will be used later if called with zero.
  Default parameter is 1 (no correction) optionally 2, 3, 4
  In case of exponentiation new version 5 is needed in most cases.
  Definition given here will be thus overwritten in such a case
  below in routine PHOCIN
*******************************************************************************/
        if(!phokey_.iexp) initializeKinematicCorrections(1);
        else              setExponentiation(true);

        int buf=1;
        iphqrk_(&buf); // Blocks emission from quarks if buf=1 (default); enables if buf=2
                       // Physical treatment will be 3, option 2 is not realistic and for tests only
        buf=2;
        iphekl_(&buf); // Blocks emission in  pi0 to gamma e+ e- if parameter is >1 (enables otherwise)

// Initialize status counter for warning messages
        for(int i=0;i<10;i++) phosta_.status[i]=0;
// elementary security level, should remain 1 but we may want to have a method to change.
        phosta_.ifstop=1; 

        pholun_.phlun=6; // Logical output unit for printing error messages

// Define pi and 2*pi
        phpico_.pi   =3.14159265358979324;
        phpico_.twopi=6.28318530717958648;

// Further initialization done automatically
// see places with - VARIANT A - VARIANT B - all over to switch between options

//----------- SLOWER VARIANT A, but stable ------------------
//--- it is limiting choice for small XPHCUT in fixed orer
//--- modes of operation

// Best choice is if FINT=2**N where N+1 is maximal number
// of charged daughters
// see report on overweihted events
        if(phokey_.interf) maxWtInterference(2.0);
        else               maxWtInterference(1.0);

/*
----------- FASTER VARIANT B  ------------------
-- it is good for tests of fixed order and small XPHCUT
-- but is less promising for more complex cases of interference
-- sometimes fails because of that

        if(phokey_.interf) maxWtInterference(1.8);
        else               maxWtInterference(0.0);
------------END VARIANTS A B -----------------------
*/


//------------------------------------------------------------------------------
// Print PHOTOS header
//------------------------------------------------------------------------------
        int                coutPrec = cout.precision(6);
        ios_base::fmtflags flags    = cout.setf(ios_base::floatfield);
        cout<<endl;
        cout<<"********************************************************************************"<<endl<<endl;

        cout<<"                            ========================="<<endl;
        cout<<"                              PHOTOS, Version:  "<<VER_MAJOR<<"."<<VER_MINOR<<endl;
        cout<<"                              Released at:  "<<DAT_DAY<<"/"<<DAT_MONTH<<"/"<<DAT_YEAR<<endl;
        cout<<"                            ========================="<<endl<<endl;

        cout<<"                     Photos QED corrections in Particle Decays"<<endl<<endl;

        cout<<"           Monte Carlo Program - by E. Barberio, B. van Eijk and Z. Was"<<endl;
        cout<<"           From version 2.09   - by P. Golonka and Z. Was"<<endl;
        cout<<"           From version 3.00   - by N. Davidson, T. Przedzinski and Z. Was"<<endl;

        cout<<"********************************************************************************"<<endl<<endl;

        cout<<"                  Internal (default) input parameters: "<<endl<<endl;
        cout<<"                    INTERF= "<<phokey_.interf<<" ISEC= " <<phokey_.isec <<" ITRE= "<<phokey_.itre
                             <<" IEXP= "  <<phokey_.iexp  <<" IFTOP= "<<phokey_.iftop<<" IFW= " <<phokey_.ifw <<endl;
        cout<<"                    ALPHA_QED= "<<phocop_.alpha<<" XPHCUT= "<<phocop_.xphcut<<endl<<endl;

        if(phokey_.interf) cout<<"                    Option with interference is active"<<endl;
        if(phokey_.isec)   cout<<"                    Option with double photons is active"<<endl;
        if(phokey_.itre)   cout<<"                    Option with triple/quatric photons is active"<<endl;
        if(phokey_.iexp)   cout<<"                    Option with exponentiation is active EPSEXP="<<phokey_.expeps<<endl;
        if(phokey_.iftop)  cout<<"                    Emision in t tbar production is active"<<endl;
        if(phokey_.ifw)    cout<<"                    Correction wt in decay of W is active"<<endl;
        if(meCorrectionWtForZ)    cout<<"                    ME in decay of Z is active"<<endl;
        if(meCorrectionWtForW)    cout<<"                    ME in decay of W is active"<<endl;
        cout<<endl<<"          WARNING:  /HEPEVT/ is not anymore used."<<endl<<endl;
/*
        cout<<endl<<"            WARNING (1): /HEPEVT/ is not anymore the standard common block"<<endl<<endl;

        cout<<"            PHOTOS expects /HEPEVT/ to have REAL*8 variables. To change to"<<endl;
        cout<<"            REAL*4 modify its declaration in subr. PHOTOS_GET PHOTOS_SET:"<<endl;
        cout<<"                 REAL*8  d_h_phep,  d_h_vhep"<<endl;
        cout<<"            WARNING (2): check dims. of /hepevt/ /phoqed/ /ph_hepevt/."<<endl;
        cout<<"            HERE:                     d_h_nmxhep=10000 and  NMXHEP=10000"<<endl<<endl;
*/
        cout<<"********************************************************************************"<<endl;
        // Revert output stream flags and precision
        cout.precision(coutPrec);
        cout.flags    (flags);

        // Add reggeon, pomeron and its diffractive states to the list
        // of decays where bremsstrahlung is suppressed
        Photos::suppressBremForDecay(0,110);
        Photos::suppressBremForDecay(0,990);
        Photos::suppressBremForDecay(0,9902110);
        Photos::suppressBremForDecay(0,9902210);
        Photos::suppressBremForDecay(0,9900110);
        Photos::suppressBremForDecay(0,9900210);
        Photos::suppressBremForDecay(0,9900220);
        Photos::suppressBremForDecay(0,9900330);
        Photos::suppressBremForDecay(0,9900440);

// Initialize Marsaglia and Zaman random number generator
        PhotosRandom::initialize();
}

static void Photospp::Photos::initialize ( )

static

Initalize Photos with the parameters previously set via the setter methods

static void Photospp::Photos::initializeKinematicCorrections ( int  flag )

inlinestatic

Initialize kinematic corrections

Definition at line 125 of file 3.52/src/photosCInterfaces/Photos.h.

References phcork_().

{ phcork_(&flag); }

static void Photospp::Photos::initializeKinematicCorrections ( int  flag )

inlinestatic

Initialize kinematic corrections

Definition at line 125 of file Photos.h.

References phcork_().

Referenced by initialize(), and setExponentiation().

{ phcork_(&flag); }

bool Photos::isStatusCodeIgnored ( int  status )

static

Returns 'true' if status code is ignored

Definition at line 333 of file Photos.cxx.

References i, and ignoreStatusCodeList.

Referenced by Photospp::PhotosHepMCParticle::checkMomentumConservation(), and Photospp::PhotosHepMCParticle::getDaughters().

{
  if(!ignoreStatusCodeList) return false;

  for(unsigned int i=0;i<ignoreStatusCodeList->size();i++)
    if( status==ignoreStatusCodeList->at(i) ) return true;

  return false;
}

static bool Photospp::Photos::isStatusCodeIgnored ( int  status )

static

Returns 'true' if status code is ignored

static void Photospp::Photos::maxWtInterference ( double  interference )

inlinestatic

Maximum interference weight

Definition at line 82 of file Photos.h.

References phokey_.

Referenced by initialize().

{ phokey_.fint=interference; }

static void Photospp::Photos::maxWtInterference ( double  interference )

inlinestatic

Maximum interference weight

Definition at line 82 of file 3.52/src/photosCInterfaces/Photos.h.

References phokey_.

{ phokey_.fint=interference; }

Photos& Photospp::Photos::operator= ( const Photos &  )

private

Photos& Photospp::Photos::operator= ( const Photos &  )

private

void Photos::processBranch ( PhotosParticle *  p )

static

Process decay of whole decay branch starting from given particle

Definition at line 218 of file Photos.cxx.

References Photospp::PhotosBranch::createBranches(), Photospp::PhotosParticle::getDecayTree(), and i.

{
        vector<PhotosParticle *> particles = p->getDecayTree();
        vector<PhotosBranch *>   branches = PhotosBranch::createBranches(particles);
        for(int i=0;i<(int)branches.size();i++) branches.at(i)->process();
}

static void Photospp::Photos::processBranch ( PhotosParticle *  p )

static

Process decay of whole decay branch starting from given particle

static void Photospp::Photos::processParticle ( PhotosParticle *  p )

static

Process decay of single particle

void Photos::processParticle ( PhotosParticle *  p )

static

Process decay of single particle

Definition at line 212 of file Photos.cxx.

References b, Photospp::PhotosBranch::getSuppressionStatus(), and Photospp::PhotosBranch::process().

{
        PhotosBranch b(p);
        if(!b.getSuppressionStatus()) b.process();
}

static void Photospp::Photos::setAlphaQED ( double  alpha )

inlinestatic

Coupling constant alpha QED

Definition at line 88 of file Photos.h.

References alpha, and phocop_.

Referenced by Photos().

{ phocop_.alpha=alpha; }

static void Photospp::Photos::setAlphaQED ( double  alpha )

inlinestatic

Coupling constant alpha QED

Definition at line 88 of file 3.52/src/photosCInterfaces/Photos.h.

References alpha, and phocop_.

{ phocop_.alpha=alpha; }

static void Photospp::Photos::setCorrectionWtForW ( bool  corr )

inlinestatic

Definition at line 100 of file Photos.h.

References phokey_.

Referenced by Photos().

{ phokey_.ifw=(int)corr; }

static void Photospp::Photos::setCorrectionWtForW ( bool  corr )

inlinestatic

Definition at line 100 of file 3.52/src/photosCInterfaces/Photos.h.

References phokey_.

{ phokey_.ifw=(int)corr; }

static void Photospp::Photos::setDoubleBrem ( bool  doub )

inlinestatic

Set double bremsstrahlung generation

Definition at line 94 of file Photos.h.

References phokey_.

Referenced by Photos(), and setExponentiation().

{ phokey_.isec=(int)doub; }

static void Photospp::Photos::setDoubleBrem ( bool  doub )

inlinestatic

Set double bremsstrahlung generation

Definition at line 94 of file 3.52/src/photosCInterfaces/Photos.h.

References phokey_.

{ phokey_.isec=(int)doub; }

void Photos::setExponentiation ( bool  expo )

static

Set exponentiation mode

Definition at line 349 of file Photos.cxx.

References initializeKinematicCorrections(), phokey_, setDoubleBrem(), setInfraredCutOff(), and setQuatroBrem().

Referenced by initialize().

{
        phokey_.iexp = (int) expo;
        if(expo)
        {
                setDoubleBrem(false);
                setQuatroBrem(false);
                setInfraredCutOff(0.0000001);
                initializeKinematicCorrections(5);
                phokey_.expeps=0.0001;
        }
}

static void Photospp::Photos::setExponentiation ( bool  expo )

static

Set exponentiation mode

static void Photospp::Photos::setInfraredCutOff ( double  cut_off )

inlinestatic

Minimal energy (in units of decaying particle mass) for photons to be explicitly generated

Definition at line 85 of file 3.52/src/photosCInterfaces/Photos.h.

References phocop_.

{ phocop_.xphcut=cut_off; }

static void Photospp::Photos::setInfraredCutOff ( double  cut_off )

inlinestatic

Minimal energy (in units of decaying particle mass) for photons to be explicitly generated

Definition at line 85 of file Photos.h.

References phocop_.

Referenced by Photos(), and setExponentiation().

{ phocop_.xphcut=cut_off; }

static void Photospp::Photos::setInterference ( bool  interference )

inlinestatic

Key for interference, matrix element weight

Definition at line 91 of file Photos.h.

References phokey_.

Referenced by Photos().

{ phokey_.interf=(int)interference; }

static void Photospp::Photos::setInterference ( bool  interference )

inlinestatic

Key for interference, matrix element weight

Definition at line 91 of file 3.52/src/photosCInterfaces/Photos.h.

References phokey_.

{ phokey_.interf=(int)interference; }

static void Photospp::Photos::setMeCorrectionWtForScalar ( bool  corr )

static

Switch for complete effects of matrix element (in scalar to 2 scalars decays)

void Photos::setMeCorrectionWtForScalar ( bool  corr )

static

Switch for complete effects of matrix element (in scalar to 2 scalars decays)

Definition at line 371 of file Photos.cxx.

References corr, and meCorrectionWtForScalar.

{
        meCorrectionWtForScalar=corr;
}

void Photos::setMeCorrectionWtForW ( bool  corr )

static

Switch for complete effects of matrix element (in leptonic W decays)

Definition at line 362 of file Photos.cxx.

References corr, and meCorrectionWtForW.

{
        meCorrectionWtForW=corr;
}

static void Photospp::Photos::setMeCorrectionWtForW ( bool  corr )

static

Switch for complete effects of matrix element (in leptonic W decays)

void Photos::setMeCorrectionWtForZ ( bool  corr )

static

Switch for complete effects of matrix element (in leptonic Z decays)

Definition at line 367 of file Photos.cxx.

References corr, and meCorrectionWtForZ.

{
        meCorrectionWtForZ=corr;
}

static void Photospp::Photos::setMeCorrectionWtForZ ( bool  corr )

static

Switch for complete effects of matrix element (in leptonic Z decays)

static void Photospp::Photos::setMomentumConservationThreshold ( double  threshold )

inlinestatic

set energy momentum conservation threshold

Definition at line 150 of file Photos.h.

References momentum_conservation_threshold, and threshold.

{momentum_conservation_threshold=threshold; }

static void Photospp::Photos::setMomentumConservationThreshold ( double  threshold )

inlinestatic

set energy momentum conservation threshold

Definition at line 150 of file 3.52/src/photosCInterfaces/Photos.h.

References momentum_conservation_threshold, and threshold.

{momentum_conservation_threshold=threshold; }

static void Photospp::Photos::setQuatroBrem ( bool  quatroBrem )

inlinestatic

Set bremsstrahlung generation up to multiplicity of 4

Definition at line 97 of file Photos.h.

References phokey_.

Referenced by Photos(), and setExponentiation().

{ phokey_.itre=(int)quatroBrem; }

static void Photospp::Photos::setQuatroBrem ( bool  quatroBrem )

inlinestatic

Set bremsstrahlung generation up to multiplicity of 4

Definition at line 97 of file 3.52/src/photosCInterfaces/Photos.h.

References phokey_.

{ phokey_.itre=(int)quatroBrem; }

static void Photospp::Photos::setRandomGenerator ( double(*)()  gen )

static

Substitute build-in generator with external one

void Photos::setRandomGenerator ( double(*)()  gen )

static

Substitute build-in generator with external one

Definition at line 343 of file Photos.cxx.

References randomDouble, and Photospp::PhotosRandom::randomReal().

{
  if(gen==NULL) randomDouble = PhotosRandom::randomReal;
  else          randomDouble = gen;
}

static void Photospp::Photos::setSeed ( int  iseed1, int  iseed2  )

inlinestatic

Seed for RANMAR used by fortran part of the Photos

Definition at line 79 of file 3.52/src/photosCInterfaces/Photos.h.

References Photospp::PhotosRandom::setSeed().

{ PhotosRandom::setSeed(iseed1,iseed2); }

static void Photospp::Photos::setSeed ( int  iseed1, int  iseed2  )

inlinestatic

Seed for RANMAR used by fortran part of the Photos

Definition at line 79 of file Photos.h.

References Photospp::PhotosRandom::setSeed().

{ PhotosRandom::setSeed(iseed1,iseed2); }

static void Photospp::Photos::setStopAtCriticalError ( bool  stop )

static

void Photos::setStopAtCriticalError ( bool  stop )

static

Definition at line 376 of file Photos.cxx.

References Photospp::Log::Info(), and phosta_.

{
        phosta_.ifstop=(int)stop;
        if(!stop)
        {
                Log::Info()<<"PHOTOS production mode. Elementary test of data flow from event record disabled. "<<endl
                           <<"Prior checks of the complete configuration "<<endl
                           <<"(for the particular set of input parameters) must have been done! "<<endl;
        }
}

static void Photospp::Photos::setTopProcessRadiation ( bool  top )

inlinestatic

Set photon emission in top pair production in quark (gluon) pair annihilation

Definition at line 115 of file 3.52/src/photosCInterfaces/Photos.h.

References phokey_.

{ phokey_.iftop=(int)top; }

static void Photospp::Photos::setTopProcessRadiation ( bool  top )

inlinestatic

Set photon emission in top pair production in quark (gluon) pair annihilation

Definition at line 115 of file Photos.h.

References phokey_.

Referenced by Photos().

{ phokey_.iftop=(int)top; }

static void Photospp::Photos::suppressAll ( )

inlinestatic

Suppress all processing. Only forced decays will be processed.

Definition at line 54 of file 3.52/src/photosCInterfaces/Photos.h.

References isSuppressed.

{ isSuppressed=true; }

static void Photospp::Photos::suppressAll ( )

inlinestatic

Suppress all processing. Only forced decays will be processed.

Definition at line 54 of file Photos.h.

References isSuppressed.

{ isSuppressed=true; }

static void Photospp::Photos::suppressBremForBranch ( int  count, int  motherID,   ...  )

static

Suppress processing of whole decay branch

void Photos::suppressBremForBranch ( int  count, int  motherID,   ...  )

static

Suppress processing of whole decay branch

Definition at line 241 of file Photos.cxx.

References count, i, supBremList, and v.

{
        va_list arg;
        va_start(arg, motherID);
        vector<int> *v = new vector<int>();
        v->push_back(motherID);
        for(int i = 0;i<count;i++)
        {
                v->push_back(va_arg(arg,int));
        }
        va_end(arg);
        v->push_back(1);
        if(!supBremList) supBremList = new vector< vector<int>* >();
        supBremList->push_back(v);
}

static void Photospp::Photos::suppressBremForDecay ( int  count, int  motherID,   ...  )

static

Suppress processing of a single decay

void Photos::suppressBremForDecay ( int  count, int  motherID,   ...  )

static

Suppress processing of a single decay

Definition at line 225 of file Photos.cxx.

References count, i, supBremList, and v.

Referenced by initialize().

{
        va_list arg;
        va_start(arg, motherID);
        vector<int> *v = new vector<int>();
        v->push_back(motherID);
        for(int i = 0;i<count;i++)
        {
                v->push_back(va_arg(arg,int));
        }
        va_end(arg);
        v->push_back(0);
        if(!supBremList) supBremList = new vector< vector<int>* >();
        supBremList->push_back(v);
}

Member Data Documentation

static Photos Photos::_instance

staticprivate

Definition at line 201 of file Photos.h.

Referenced by getInstance().

static const int Photos::DAT_DAY =4

static

Definition at line 32 of file Photos.h.

Referenced by initialize().

static const int Photos::DAT_MONTH =11

static

Definition at line 32 of file Photos.h.

Referenced by initialize().

static const int Photos::DAT_YEAR =12

static

Definition at line 32 of file Photos.h.

Referenced by initialize().

static vector< vector< int > * > * Photos::forceBremList = 0

static

List of forced decays

Definition at line 163 of file Photos.h.

Referenced by Photospp::PhotosBranch::checkList(), Photospp::PhotosBranch::createBranches(), forceBremForBranch(), and forceBremForDecay().

static vector< pair< int, double > * > * Photos::forceMassList = 0

static

List of forced mass values

Definition at line 166 of file Photos.h.

Referenced by forceMass(), and forceMassFromEventRecord().

static int Photos::historyEntriesStatus = 3

static

Status of history entries

Definition at line 187 of file Photos.h.

Referenced by createHistoryEntries(), and Photospp::PhotosHepMCParticle::createHistoryEntry().

static vector< int > * Photos::ignoreStatusCodeList = 0

static

List of ignored status codes

Definition at line 169 of file Photos.h.

Referenced by deIgnoreParticlesOfStatus(), ignoreParticlesOfStatus(), and isStatusCodeIgnored().

static bool Photos::isCreateHistoryEntries =false

static

Flag for creating historic entries

Definition at line 184 of file Photos.h.

Referenced by createHistoryEntries().

static bool Photos::isSuppressed =false

static

Is in suppressed mode

Definition at line 154 of file Photos.h.

Referenced by Photospp::PhotosBranch::createBranches(), and suppressAll().

static bool Photos::massFrom4Vector =true

static

Is mass from 4-vector or from event record

Definition at line 157 of file Photos.h.

Referenced by forceMassFrom4Vector().

static bool Photos::meCorrectionWtForScalar =false

static

Flag for complete effects of matrix element (in scalars decays)

Definition at line 175 of file Photos.h.

Referenced by iniInfo(), and setMeCorrectionWtForScalar().

static bool Photos::meCorrectionWtForW =false

static

Flag for complete effects of matrix element (in leptonic W decays)

Definition at line 181 of file Photos.h.

Referenced by iniInfo(), initialize(), and setMeCorrectionWtForW().

static bool Photos::meCorrectionWtForZ =false

static

Flag for complete effects of matrix element (in leptonic Z decays)

Definition at line 178 of file Photos.h.

Referenced by iniInfo(), initialize(), and setMeCorrectionWtForZ().

static double Photos::momentum_conservation_threshold = 0.1

static

Threshold for momentum conservation check

Definition at line 172 of file Photos.h.

Referenced by Photospp::PhotosHepMCParticle::checkMomentumConservation(), and setMomentumConservationThreshold().

static double(* Photos::randomDouble)() = PhotosRandom::randomReal

static

Pointer to random generator function

Definition at line 190 of file Photos.h.

Referenced by Photospp::phoranc_(), and setRandomGenerator().

static vector< vector< int > * > * Photos::supBremList = 0

static

List of suppressed decays

Definition at line 160 of file Photos.h.

Referenced by Photospp::PhotosBranch::checkList(), suppressBremForBranch(), and suppressBremForDecay().

static const int Photos::VER_MAJOR =3

static

Definition at line 31 of file Photos.h.

Referenced by initialize().

static const int Photos::VER_MINOR =52

static

Definition at line 31 of file Photos.h.

Referenced by initialize().

---

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

• Photos.h
• 3.52/src/photosCInterfaces/Photos.h
• Photos.cxx