]>
Commit | Line | Data |
---|---|---|
1 | void lumi () | |
2 | { | |
3 | // Dynamically link some shared libs | |
4 | gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libLHC"); | |
5 | ||
6 | // | |
7 | // LHC | |
8 | // | |
9 | AliLHC* lhc = new AliLHC(); | |
10 | lhc->SetRadius(2665887./2./TMath::Pi()); | |
11 | lhc->SetAverageBeta(5050.); | |
12 | lhc->SetAverageDisp(129.); | |
13 | lhc->SetSetUpTime(1./3. * 3600.); | |
14 | lhc->SetFillingTime(3. * 3600.); | |
15 | // | |
16 | // The Interaction regions | |
17 | // | |
18 | // | |
19 | ||
20 | Float_t bstar = 50.; | |
21 | // | |
22 | AliLhcIRegion* alice = new AliLhcIRegion(lhc, "ALICE", "ALICE"); | |
23 | alice->SetBetaStar(bstar); | |
24 | AliLhcIRegion* cms = new AliLhcIRegion(lhc, "CMS" , "CMS" ); | |
25 | cms->SetBetaStar(bstar); | |
26 | AliLhcIRegion* atlas = new AliLhcIRegion(lhc, "ATLAS", "ATLAS"); | |
27 | atlas->SetBetaStar(bstar); | |
28 | ||
29 | lhc->AddIRegion(alice); | |
30 | // lhc->AddIRegion(cms); | |
31 | // lhc->AddIRegion(atlas); | |
32 | // | |
33 | // The beams | |
34 | // | |
35 | Float_t n = 6.8e7; | |
36 | Float_t epsH = 1.5e-4; | |
37 | Float_t epsL = 2.5e-9; | |
38 | Float_t de = 1.14e-4; | |
39 | Float_t energy = 7000; | |
40 | Int_t a1 = 208; | |
41 | Int_t z1 = 82; | |
42 | Int_t a2 = 208; | |
43 | Int_t z2 = 82; | |
44 | ||
45 | ||
46 | // | |
47 | AliLhcBeam* beam1 = new AliLhcBeam(lhc); | |
48 | beam1->SetN(n); | |
49 | beam1->SetNEmittance(epsH); // (cm) | |
50 | beam1->SetLongEmittance(epsL); // (GeV s) | |
51 | beam1->SetEnergy(energy); // (GeV) | |
52 | beam1->SetParticle(a1,z1); | |
53 | beam1->SetEnergySpread(de); | |
54 | ||
55 | AliLhcBeam* beam2 = new AliLhcBeam(lhc); | |
56 | beam2->SetN(n); | |
57 | beam2->SetNEmittance(epsH); // (cm) | |
58 | beam2->SetLongEmittance(epsL); // (GeV s) | |
59 | beam2->SetEnergy(energy); // (GeV) | |
60 | beam2->SetParticle(a2,z2); | |
61 | beam2->SetEnergySpread(de); | |
62 | ||
63 | ||
64 | lhc->SetBeams(beam1, beam2); | |
65 | // | |
66 | // The Processes | |
67 | // | |
68 | AliLhcProcessBB* bb = new AliLhcProcessBB(lhc, "BB", "Beam-Beam Losses"); | |
69 | bb->SetCrossSection(505.); | |
70 | AliLhcProcessIBS* ibs = new AliLhcProcessIBS(lhc, "IBS", "Intra Beam Scattering"); | |
71 | AliLhcProcessBT* bt = new AliLhcProcessBT(lhc, "BT", "Beta* Tuning"); | |
72 | bt->SetBetaMin(50.); | |
73 | ||
74 | lhc->AddProcess(ibs); | |
75 | lhc->AddProcess(bb); | |
76 | // lhc->AddProcess(bt); | |
77 | // | |
78 | // Run the collider | |
79 | // | |
80 | lhc->SetTime(100., 20.*3600.); | |
81 | lhc->Init(); | |
82 | lhc->EvolveTime(); | |
83 | lhc->Evaluate(); | |
84 | } | |
85 | ||
86 | ||
87 | ||
88 | ||
89 | ||
90 | ||
91 | ||
92 | ||
93 | ||
94 |