Gas pressure peak in TDI removed.
[u/mrichter/AliRoot.git] / LHC / lumi.C
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