Диссертация (1149195), страница 19
Текст из файла (страница 19)
.944.15 Отношение средней множественности к квадрату характерногопоперечного импульса от одного источника в EPEM-модели.4.16 Коэффициент pt-pt корреляции в модели THERMINATOR 2.. . .95. .971244.17 Зависимость коэффициента pt-pt корреляций от центральности вмодели слияния струн (слева) и модели расталкивающихся струн(справа) в Pb-Pb столкновениях. . . .
. . . . . . . . . . . . . . . .98125Список таблиц2.1Распределения ядерной плотности и параметры, использованныев данной работе [86]. . . . . . . . . . . . . . . . . . . . . . . . . . .372.2Результаты фиксации параметров в монте-карловской модели. .443.1Параметры фитирования коэффициентов корреляции . . . . . . .57126Приложение AЗаголовочый файлмонте-карловской моделиDsfmGen.h#include <v e c t o r >2 #include <i o s t r e a m >13#include " T P a r t i c l e . h"5 #include " TParticlePDG .
h"6 #include "TTree . h"7 #include "TH1F . h"8 #include "TF1 . h"9 #include <math . h>10 #include "TMath . h"11 #include "Math/ I F u n c t i o n . h"12 #include "TDatabasePDG . h"4// F i n a l1314*/ #i n c l u d e< f l o a t . h>15 #i n c l u d e<math . h>16 #i n c l u d e<memory . h>17 #i n c l u d e< a s s e r t . h>18 #i n c l u d e< s t d l i b . h>19 * /202122using namespace//quarktypestd ;defenition#define QUARK_VALENCE 024 #define QUARK_DIQUARK 123Particlesformat127#define QUARK_SEA 226 #define QUARK_ANTISEA 32527282930structQuark {Double_tx, y;//quark31Double_tz;//momentum32Int_t type ;//quark33Int_t used ;//flagpositionfractiontypeistrueifquarkisalreadyincludedinthestringformation34 } ;353637typedef structNucleon {Double_t x , y ;38UInt_t baseNumber ;39UInt_t count ;40Bool_t wounded ;41Int_t type ;42//Parton////coordinatesnumbernumberof//ofthebeginninginquarkarrayquarksis// 0itparticipant ?for// ! ! !proton ,// I t1forneutronDOES NOT w o r k !43 } Nucleon ;44454647typedef structDipole {Quark * q1 , * q2 ;// QUARK_SEA QUARK_ANTISEA o r QUARK_VALENCE QUARK_DIQUARKBool_t used ;//−use48UInt_t n p a r e n t ;49Double_t GetLength2 ( ) {if50515253}55565758isifusedthisforisatipcollisionofstringmatrixconstructingcout<<"ERROR: D i p o l e Length i s not d e f i n e d "<<e n d l ;return0;( ( q1−>x−q2−>x ) * ( q1−>x−q2−>x ) + ( q1−>y−q2−>y ) * ( q1−>y−q2−>y ) ) ;}}59Double_t GetP02 ( ) {6061// i tflag( q1==0 | | q2==0) {else {return54the−}return3 .
9 0 0 6 2 5 e −02/GetLength2 ( ) ;62D i p o l e ( ) { q1 =0; q2 = 0 ; } ;63~ D i p o l e ( ) { q1 =0; q2 = 0 ; } ;64 } D i p o l e ;65666768typedef structInt_t s i d e ;Edge {//1//Rapidityor−1edgeforstringfusionalgorithm12869Double_t y ;70UInt_t g r i d c e l l ;71Double_t p2 ;// r a p i d i t y// f r o m// w h a twhichpt2gridastringtipscellthisedgeisformedcarries72 } Edge ;737475typedef structString {Quark * q1 , * q2 ;//two76Double_t x , y ;//string77Double_t p02 ;78Double_t ymin , ymax ;79S t r i n g * next ;80Char_t noform ;//////next//ofcoordinatestransverserapiditystringitstringisof(asingle )stringbordersinequalmomentumthe1clusterifthestringcannotbeformed81 } S t r i n g ;8283class Nucleus {85public :84868788Nucleus ( ) ;void// c o n s t r u c t o rSetA ( UInt_t xA) ;Nucleus ( UInt_t xA) ;// s e t// c o n s t r u c t o r89v e c t o r <Nucleon> Nucleons ;90UInt_t A, Z ,N;91Double_t x ;//92Double_t y ;//93TF1 * NucDist ;94UInt_t maxcount ;95Double_t lambda ;9697voidvoid////////AonlyGenerateNucleons ( ) ;numberofintonucleonscoordinatesofofsea// g e n e r a t eMoveNucleons ( ) ; // t a k i n gofworks ! !DistributionnumberAarrayTransversemaximum// meanAwithnucleonspairsinoverar−S a x o n )( Woodsnucleonpairs ;nucleonaccountpositionsA, x , yactualizev e c t o r <N u c l e o n >Nucleons98UInt_t m;// t o t a lnumberofdipoles \ pairs99 } ;100101class Glauber {103public :102104Glauber ( ) ;105Glauber ( Double_t cmsigma ) ;106Double_t msigma ;107UInt_t NpartA ;// NumberofparticipantinnucleusA108UInt_t NpartB ;// NumberofparticipantinnucleusB109UInt_t NpartAB ; // Numberofparticipantinboth110UInt_t N c o l l ;ofbinary111UInt_t AA,AB;112void// c o n s t r u c t o r// s i g m a// Number// c o n s t r u c t o rwithsigmamultipliermultipliernucleicollisionsProceed ( v e c t o r <Nucleon> & A, v e c t o r <Nucleon> & B) ;129113v e c t o r <UInt_t> nSumA ;114v e c t o r <UInt_t> nSumB ;115v e c t o r <UInt_t> I ;116Double_t sigmaGlauber ;117 } ;118119120121structPrimaryParticle {Double_ty;122Double_tphi ;// p h i123Double_tpt ;// p t124Double_tkstring ;125Char_tcharge ;127UInt_tgridcell ;128Int_tpdg ;129Double_t m;// mass// r a p i d i t y// t e n s i o n// c h a r g eofthestring ,producedtheparticlesign126// f r o mwhatgridcelltheparticleiscoming//PDG c o d e130 } ;131132class DsfmGen{134public :133135TDatabasePDG *136Int_t * EnabledParticlesPDG ;137Int_t E n a b l e d P a r t i c l e s N u m b e r ;138TH1F** E n a b l e d P a r t i c l e s H ;139Glauber PrevGlaub ;140Glauber Glaub ;141Double_t lambda ;142Double_t a ;143Double_t b_min , b_max , b ;144Double_t r_0 ;145Double_t rmax ;146Double_t a_s ;147Double_t r _ s t r ;148Double_t En ;149Double_t r2mean ;150151152153boolboolvoidboolgDatabasePDG ;// mean// s q u a r eanumber*ato// i m p a c t−s q u a r e d// c o n f i n e m e n t// e f f e c t i v eEnabledHard ;nucleitransverseradiuspositions ,ofandfmranges ,radiusofproton ,fmfmproton ,fmconstantstringcollision ,// e n a b l epairsparametercoupling// t r a n s v e r s eofseagenerate// mean// e n e r g yofmomentumradiusGeV( s q r t { s })transferinthedipoles .EnabledPrevGlauber ;GenerateEvent ( ) ;Event ;154DsfmGen ( ) ;// c o n s t r u c t o r155~DsfmGen ( ) ;// d e s t r u c t o r156Nucleus A, B ;containing// C o l l i d i n gtheonlyNuclei .nucleonThatcoluldbealsoaproton ,anucleus130157voidG e n e r a t e N u c l e o n s P r o p e r t i e s ( Nucleus & A, UInt_t AA) ;// g e n e r a t enucleonsm,type158v e c t o r <Quark> qA , qB ;159v e c t o r <Dipole > dipA , dipB ;160v e c t o r <Double_t> v ;161void// a r r a yofallquarksf r o m A and B// D i p o l e s ,formedof// t e m p o r a r yvectorfor// Nuc=Aquarks .quarkAlsoorB;q=qAgenerateorqB .Generate// 2Darray163Double_t prob ;havecollision166bool>168169170171172boolboolboolboolvoidofato// g e n e r a t edipCollMatrix ;AngularConserv ;state167// p r o b a b i l i t yGenerateProbMatrix ( ) ;vector<boolf i l l// 2D// w h e t h e rmatrixofUseGlauber ;theangularallmomentum// w h e t e rtheglauberUseNonGlauber ;// w h e t e r// w h e t e rGenerateCollMatrix () ;momenumisthefractionscalculationPrev .glauber−g l a u b e rtheNon// g e n e r a t etostaffbooleancollisioninnucleusA173UInt_t NpartB ;// NumberofparticipantinnucleusB174UInt_t NpartAB ; // Numberofparticipantinboth175UInt_t N c o l l ;ofbinary176UInt_t GlauberNpartA ;178179aascenarioofUInt_t G l a u b e r N c o l l ;agivenv e c t o r <UInt_t> nSumB ;182UInt_t d i p N c o l l ;v e c t o r <S t r i n g > S t r i n g s ;186UInt_t N s t r i n g s ;188189190191192193voidvoidvoidvoidvoidperformedshouldbebeperformedperformedmatrixnucleicollisionsparticipantinnucleusAaccordingtoGlauberofparticipantinnucleusBaccordingtoGlauberofparticipantinbothgivennucleiaccordingtocollisionofbinarycollisionsaccordingtoGlauberscenario//Calculate−NpartAGlaiberNcollnucCollMatrix ;185187aCalcGeometryStats ( ) ;bool>becollision181vector<initialcollision// Numberv e c t o r <UInt_t> nSumA ;voidthecollision// Numbergiven180184ofUInt_t GlauberNpartAB ; // Numberof183// NumbergivenUInt_t GlauberNpartB ;Glauberinvalueshouldparticipant// Numberequalcalculationofofofprobabilitiesconservedshould// Numberscenariopropertiesarray// P u tUsePrevGlauber ;ofthenucleonEqualFraction ;scenarioALLinteractionUInt_t NpartA ;177calculationdipoles .v e c t o r <Double_t> dipProbMatrix ;165momentumand162voidquarksGenerateQuarks ( Nucleus & A, UInt_t AA, v e c t o r <Quark>& q , v e c t o r <Dipole>&dip ) ;164nuclei////numberArrayofofavailablestringsFillStringParams ( String & string ) ;GenerateStrings () ;//stringsGenerate// f i l lstringsymin ,accordingymax ,tox,yofstringcollisionmatrixFuseStrings () ;DontFuseStrings ( ) ;GenerateGridAndEdges ( ) ;v e c t o r <S t r i n g *> g r i d ;// 2D//GeneratearrayUInt_t grid_w , grid_h , g r i d _ s i z e ;ofstringspointers// S i z eoftogridaccordingfirsttostringcollisionineachmatrixcell131194Double_t d g r i d ;195v e c t o r <Edge> e d g e s ;196v e c t o r <UChar_t> G r i d L o c a l ;interval .Used// s t r i n ginedges// 2DPrimaries( 1Darrayarray )ofstringnumberingivenrapidityModification197v e c t o r <P r i m a r y P a r t i c l e > perv ;198P r i m a r y P a r t i c l e G e n e r a t e P r i m a r y P a r t i c l e ( Double_t ymin , Double_t ymax , Double_t// p r i m a r yparticleseta , Double_t pt4 , UInt_t g r i d c e l l ) ;, ymax199voidwitheta−numberwithoneDouble_t p0 , mu0 ;201ULong_t N p a r t i c l e s ;202TF1 ** T r a n s v e r s e D i s t ;204GeneratePrimaries () ;ModifyPrimaries ( ) ;accord .// g e n e r a t e// p r o c e e d// v e c t o r <T P a r t i c l e >206TClonesArray * I n t e r m e d i a t e ;208209210211212 } ;voidvoidFinal ;// f i n a lwithinyminstrings ;primarywith//particlesanisotrophic// g e n e r a t eGenerateFinalParticles () ;TTree * t r e e ;// f i n a l// O u t p u tFillTree () ;// F i l lofmodificationparticlesevent .array ,afterparticles( quenching )decayarray ,andafterResonancesarray ,Intermediate// R e s o n a n c eparticlestheparticles// I n t e r m e d i a t eGenerateIntermediates () ;TClonesArray * F i n a l ;voidparticlet o M.















