LHC related code. First commit.
[u/mrichter/AliRoot.git] / LHC / AliLhcProcessBB.cxx
CommitLineData
11141716 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/*
17$Log$
18*/
19
20#include "AliLhcProcessBB.h"
21#include "AliLHC.h"
22#include "AliLhcIRegion.h"
23#include "AliLhcBeam.h"
24
25ClassImp(AliLhcProcessBB)
26
27AliLhcProcessBB::AliLhcProcessBB(AliLHC* lhc, const char* name, const char* title)
28 :AliLhcProcess(lhc,name,title)
29{
30// Constructor
31}
32
33
34AliLhcProcessBB::~AliLhcProcessBB()
35{
36// Destructor
37
38}
39
40void AliLhcProcessBB::Init()
41{
42 // Initialization
43 printf("\n Initializing Process %s", GetName());
44 printf("\n ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^");
45 fIRegions = fAccelerator->IRegions();
46 fBeam1 = fAccelerator->Beam(0);
47 fBeam2 = fAccelerator->Beam(1);
48}
49
50void AliLhcProcessBB::Evolve(Float_t dt)
51{
52 printf("\n Here process %s %f:", GetName(), dt);
53 TIter next(fIRegions);
54 AliLhcIRegion *region;
55 //
56 // Loop over generators and initialize
57 while((region = (AliLhcIRegion*)next())) {
58 Float_t rate = region->Luminosity()*fCrossSection;
59 Float_t loss = rate*dt;
60 fBeam1->RemoveParticles(loss);
61 fBeam2->RemoveParticles(loss);
62 }
63}
64
65
66AliLhcProcessBB& AliLhcProcessBB::operator=(const AliLhcProcessBB & rhs)
67{
68// Assignment operator
69 return *this;
70}
71
72