/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N S D U V I D E : */
/* */
/* */
/* Author of '$xrq/vide.L6$I' : */
/* */
/* Jean-Francois Colonna (LACTAMME, 1992??????????). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U V I D E : */
/* */
/* */
/* Version 1 du modele : */
/* */
/* Dans ce modele, les particules REELles sont */
/* mobiles, ce qui est garanti au cours du */
/* temps grace a la conjonction de deux phenomenes : */
/* */
/* 1-'CONSERVATION_DU_MODULE_DE_LA_VITESSE_D_UN_QUARK_OU_D_UN_ANTI_QUARK(...)' */
/* qui fait que lors d'un evenement du type 'STRONG_ABSORBTION_QxG___Q(...)' ou */
/* 'STRONG_EMISSION_Q___QxG(...)' un quark ou un anti-quark conserve le module */
/* de sa vitesse. Ainsi, si il etait nul (et donc si il etait immobile), il le */
/* restera tout au long de la simulation. */
/* */
/* 2-les particules REELles ne sont pas soumises aux evenements du type */
/* 'STRONG_TRANSFORMATION_QxAQ___G(...)', seuls le sont les particules */
/* VIRTUELles. Les particules REELles ne peuvent donc pas disparaitre... */
/* */
/* */
/* Version 2 du modele : */
/* */
/* Dans ce modele, les particules REELles sont */
/* mobiles, mais ce, independamment des processus */
/* physiques de l'interaction forte. Cela est implante */
/* grace a deux phenomenes : */
/* */
/* 1-les particules dits REELles ne sont plus visualises */
/* implicitement, */
/* */
/* 2-'ENSEMBLE_DES_QUARKS(...)' et 'ENSEMBLE_DES_ANTI_QUARKS(...)' */
/* sont remplaces respectivement par 'ENSEMBLE_DES_QUARKS_VIRTUELS(...)' */
/* et 'ENSEMBLE_DES_ANTI_QUARKS_VIRTUELS(...)' dans les modules de */
/* deplacement ('DEPLACEMENT_DE_L_ENSEMBLE_DES_QUARKS' et */
/* 'DEPLACEMENT_DE_L_ENSEMBLE_DES_ANTI_QUARK' respectivement), */
/* ainsi que dans les modules d'evenements conditionnels */
/* ('EVENEMENTS_CONDITIONNELS_POUR_CERTAINS_QUARKS' et */
/* 'EVENEMENTS_CONDITIONNELS_POUR_CERTAINS_ANTI_QUARKS' */
/* respectivement). */
/* */
/* */
/*************************************************************************************************************************************/
#nodefine __VERSION__LES_PARTICULES_REELLES_SONT_SOUMISES_A_L_INTERACTION_FORTE
/* Choix du mode de gestion des particules REELLEs (voir la note ci-dessus...). */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U V I D E P A R U N C H A M P U N I F O R M E : */
/* */
/*************************************************************************************************************************************/
#define PONDERATION_g0 \
FRA10(FRA10(FU))
DEFV(Local,DEFV(Float,INIT(ponderation_g0,PONDERATION_g0)));
/* Definition l'amplitude du champ du vide. */
#define MAXIMUM_DU_CHAMP_GLOBAL \
FU \
/* Valeur maximale du champ global qui sera atteinte partout. ATTENTION : cette definition */ \
/* doit preceder cellede 'CHAMP_DU_VIDE(...)', car en effet, elles n'appartiennent pas a la */ \
/* meme passe... */
=define CHAMP_DU_VIDE(cX,cY,cZ) \
MUL2(ponderation_g0,MAXIMUM_DU_CHAMP_GLOBAL)
/* Definition du champ du vide. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S V A L E U R S D E L A C H A R G E D E C O U L E U R : */
/* */
/* */
/* ATTENTION : */
/* */
/* Pour des raisons liees au pre-processing */
/* les definitions qui suivent doivent se trouver */
/* apres les definitions de 'CHAMP_DE_VALENCE_DU_QUARK_REEL(...)', */
/* et 'CHAMP_DE_VALENCE_DE_L_ANTI_QUARK_REEL(...)'... */
/* */
/* */
/*************************************************************************************************************************************/
%define ROUGE_max \
COULEUR_max
/* Definition du maximum de la composante ROUGE de la charge de couleur. */
%define VERTE_max \
COULEUR_max
/* Definition du maximum de la composante VERTE de la charge de couleur. */
%define BLEUE_max \
COULEUR_max
/* Definition du maximum de la composante BLEUE de la charge de couleur. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E L A G E O M E T R I E D U C H A M P " D E F L U X " : */
/* */
/* */
/* Definition du champ de flux : */
/* */
/* FLUX(x,y,z) = 0 */
/* */
/* il est donc inutile... */
/* */
/* */
/*************************************************************************************************************************************/
#define CALCULER_LE_CHAMP_DE_FLUX \
LUNDEF
/* Par compatibilte avec '$xrq/nucleon.LI$I' et '$xrq/nucleon.LV$I'... */
#define DISTANCE_CARACTERISTIQUE_g1 \
FLOT__UNDEF
DEFV(Local,DEFV(Float,INIT(distance_caracteristique_g1,FLOT__UNDEF)));
/* Nota : cette definition est uniquement destinee a la compatibilite avec les deux fichiers */
/* '$xrq/nucleon.LW.1$I' et '$xrq/nucleon.LW.2$I'... */
DEFV(Local,DEFV(Float,INIT(cutoff_du_champ_global,FLOT__UNDEF)));
/* Nota : cette definition est uniquement destinee a la compatibilite avec les deux fichiers */
/* '$xrq/nucleon.LI$I' et '$xrq/nucleon.LV$I'... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E L A G E O M E T R I E D U C H A M P " D E C O N F I N E M E N T " */
/* D U Q U A R K E T D E L ' A N T I - Q U A R K : */
/* */
/* */
/* Definition du champ de confinement : */
/* */
/* CONFINEMENT(x,y,z) = 0 */
/* */
/* il est donc inutile... */
/* */
/* */
/*************************************************************************************************************************************/
#define NE_PAS_CALCULER_LE_CHAMP_DE_CONFINEMENT \
LUNDEF
/* Par compatibilte avec '$xrq/nucleon.LI$I' et '$xrq/nucleon.LV$I'... */
#define DISTANCE_CARACTERISTIQUE_g2 \
FLOT__UNDEF
DEFV(Local,DEFV(Float,INIT(distance_caracteristique_g2,FLOT__UNDEF)));
/* Nota : cette definition est uniquement destinee a la compatibilite avec les deux fichiers */
/* '$xrq/nucleon.LW.1$I' et '$xrq/nucleon.LW.2$I'... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U C H A M P G L O B A L D ' I N T E R A C T I O N D U V I D E : */
/* */
/*************************************************************************************************************************************/
#define MINIMUM_DU_CHAMP_GLOBAL \
FZERO \
/* Valeur minimale (relativement arbitraire) du champ global... */
DEFV(Local,DEFV(Float,INIT(champ_global,FLOT__UNDEF)));
/* Intensite du champ resultant (valence + flux + confinement). */
#define CALCUL_DU_CHAMP_GLOBAL(cX,cY,cZ,calculer_le_champ_de_flux,calculer_le_champ_de_confinement,cutoff) \
Bblock \
EGAL(champ_global \
,CHAMP_DU_VIDE(cX,cY,cZ) \
); \
/* Generation du champ du vide... */ \
Eblock \
/* Intensite du champ resultant... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* G E S T I O N D E S P A R T I C U L E S R E E L L E S : */
/* */
/*************************************************************************************************************************************/
#define CHANGEMENT_EVENTUEL_DE_QUARK_REEL(quarkD,quarkA) \
Bblock \
Eblock \
/* Passage de l'etat REEL d'un quark qui le possedait (et qui va etre detruit) a un qui */ \
/* va l'acquerir (et qui vient d'etre cree). */
#define CHANGEMENT_EVENTUEL_D_ANTI_QUARK_REEL(anti_quarkD,anti_quarkA) \
Bblock \
Eblock \
/* Passage de l'etat REEL d'un anti-quark qui le possedait (et qui va etre detruit) a un qui */ \
/* va l'acquerir (et qui vient d'etre cree). */