AliEn API fixed so we can again use TGrid::OpenDir().
[u/mrichter/AliRoot.git] / MUON / README
CommitLineData
50837721 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/* $Id$ */
17
88cb7938 18==========================================================
19Please add to this README file all information concerning
20config files, simulation, digitalization, clusterization,
21reconstruction and macro analysis
22
23==========================================================
7985603c 24 How to run a MUON generation
88cb7938 25==========================================================
26aliroot
27root [0] gAlice->Run(10,"$ALICE_ROOT/MUON/Config_MUON_test.C");
28
291 single muon of 7 GeV/c in the MUON spectrometer
30acceptance will be simulated using geant3.
31Hit information will be store in the root file in the
32execution directory.
88cb7938 33If you want to change the option or to define a new directory
34for hits, you have to do the following before:
35root [0] gAlice->SetConfigFunction("Config( \"/home/martinez/aliroot/work_NewIO/test/\" , \"box\" );");
36
37============================================================
38 How to run MUON Digitalization
39============================================================
40To process the digitalization of hits the AliRunDigitizer
41framework is used:
42
2058bf30 43AliRunDigitizer * manager = new AliRunDigitizer(1,1);
44manager->SetInputStream(0,"galice.root");
e6dd99fc 45AliMUONDigitizerv1* dMUON = new AliMUONDigitizerv1(manager)
2058bf30 46manager->AddDigitizer(dMUON);
47manager->Exec("deb");
88cb7938 48
49============================================================
50 How to run MUON RecPoint clusterization
51============================================================
52aliroot -b << EOF
278a86d5 53.includepath $ALICE_ROOT/STEER
54.includepath $ALICE_ROOT/MUON
55.L $ALICE_ROOT/MUON/MUONrawclusters.C++
bf17dbfd 56MUONrawclusters("galice.root",0,9999);
57to Load
58gSystem->Load("$ALICE_ROOT/MUON/MUONrawclusters_C.so")
88cb7938 59
88cb7938 60============================================================
61 How to run MUON Trigger
62============================================================
ba9436c6 63.includepath $ALICE_ROOT/STEER
64.includepath $ALICE_ROOT/MUON
65.L $ALICE_ROOT/MUON/MUONtrigger.C++
bf17dbfd 66MUONtrigger("galice.root",0,9999);
50837721 67To Load
68gSystem->Load("$ALICE_ROOT/MUON/MUONtrigger_C.so")
88cb7938 69============================================================
dcd2690d 70 How to run MUON Tracking reconstruction (old output)
fb76e306 71===========================================================
2058bf30 72.includepath $ALICE_ROOT/STEER
73.includepath $ALICE_ROOT/MUON
74.L $ALICE_ROOT/MUON/MUONrecoNtuple.C++
50837721 75MUONrecoNtuple(0,9999,0,"galice.root")
76To Load
77gSystem->Load("$ALICE_ROOT/MUON/MUONrecoNtuple_C.so")
88cb7938 78First event and last event are given by hand, but this will change.
79
dcd2690d 80===========================================================
81 How to run MUON Tracking reconstruction (new output)
82===========================================================
83.includepath $ALICE_ROOT/STEER
84.includepath $ALICE_ROOT/MUON
85.L $ALICE_ROOT/MUON/MUONTracker.C++
bf17dbfd 86MUONTracker("galice.root",0,9999);
50837721 87To Load
88gSystem->Load("$ALICE_ROOT/MUON/MUONTracker_C.so")
dcd2690d 89Output in MUON.Tracks.root using the new IO
90as a TClonesArray of AliMUONTrack
91
374ebd7d 92Note1.In the actual version (dec-03) the following data
93members of ALiMUONTrack: fgFitter, fEventReconstructor,
94and fTrackHitsPtr are not saved in the MUON.Track.root file
95
a3448039 96===========================================================
97 How to run MUON Trigger Track reconstruction (new output)
98===========================================================
99.includepath $ALICE_ROOT/STEER
100.includepath $ALICE_ROOT/MUON
101.L $ALICE_ROOT/MUON/MUONTriggerTracker.C++
102MUONTriggerTracker("galice.root",0,10);
103To Load
104gSystem->Load("$ALICE_ROOT/MUON/MUONTriggerTracker_C.so")
105Output in MUON.Tracks.root using the new IO
106as a TClonesArray of AliMUONTriggerTrack
107
108
88cb7938 109============================================================
110 How to run MUONCheck macro
111============================================================
112To check the content of a root data file, the MUONCheck
113provides a ascii output on screen.
114
115To compile MUONCheck.C
2058bf30 116.includepath $ALICE_ROOT/STEER
117.includepath $ALICE_ROOT/MUON
118.L $ALICE_ROOT/MUON/MUONCheck.C++
bf17dbfd 119To Load
120gSystem->Load("$ALICE_ROOT/MUON/MUONCheck_C.so")
2058bf30 121
278a86d5 122To print Kine : (default file is galice.root )
123MUONkine()
88cb7938 124
125To print hits : (default file is galice.root if not MUONhits("toto.root""); )
2058bf30 126MUONhits()
88cb7938 127
128To print digits : (default file is galice.root)
2058bf30 129MUONdigits()
88cb7938 130
131To print rawcluster : (default file is galice.root)
2058bf30 132MUONrecpoints()
88cb7938 133
134To print trigger : (default file is galice.root)
2058bf30 135MUONTestTrigger()
88cb7938 136
02d8f072 137
138============================================================
139 How to check the Geometry
140============================================================
141gAlice->Init("$ALICE_ROOT/MUON/Config_MUON_test.C");
142.L $ALICE/geant3/TGeant3/G3GUI.C
143G3GUI()
144
145
7985603c 146============================================================
147 How to run MUONdisplay
148============================================================
1eccde20 149First you need to perform a full simulation:
7985603c 150generation, digitalisation and clusterisation
151.L $ALICE_ROOT/MUON/MUONdisplay.C
152MUONdisplay(0,"galice.root")
153
02d8f072 154============================================================
155 Tracking parameters, cuts, energy loss and physics processes
156============================================================
157Tracking parameters in MUON are automatically defined by GEANT
158MUON takes the default values of CUTs and physics processes
159defined by the Config files, except for the gas mixture medium
160of the tracking chambers. The CUT's and physics processes of
161the gas mixture medium is then defined in the galice.cuts file
162in the data directory. In particular ILOSS parameter MUST be
163equal unity (1) in order simulate a realistic energy loss
164distribution (mean value and fluctuations) in the active gas.
a88eb0d0 165
166============================================================
167 Tracking of particle in the magnetic field
168============================================================
169GEANT has two ways for tracking charged particles in the
170magnetic field: HELIX et RKUTA.
171HELIX is faster and works well if the gradient of magnetic
172field is small.
173For MUON, HELIX is a not a good approximation and we must
174use RKUTA to get the optimal mass resolution of the
175spectrometer. The choice of HELIX or RKUTA is done in the
176config file when the magnetic field is defined:
177 AliMagFMaps* field = new AliMagFMaps("Maps","Maps", TRACKING, FACTOR, MAXB, AliMagFMaps::k4kG);
178 gAlice->SetField(field);
179TRACKING must be 1 for RKUTA and 2 for HELIX (the default value for aliroot is 2 (HELIX))
180FACTOR allows you to set the magnetic field to 0, just putting FACTOR=0. Default value is 1.
181MAXB is the maximum magnetic field which is 10.T
f4f795ed 182===========================================================
183 MUON cocktail for physics ..............
184===========================================================
185There is a MUON cocktail generator of the muon sources in the
186EVGEN directory. This class derives from AliGenCocktail.
187In the init of this class I have filled the cocktail with
188the muon sources: J/Psi, Upsilon, Open Charm, Open Beauty,
189Pion, Kaons. The code needs only the production cross section
190at 4pi (for the moment this values are in the code since I
191prefere them do not be modified), and the code calculates the
192rate of particles in the acceptance, making the scaling based
193on the number of collisions for the hard probes and on the
194number of participants for soft sources: Pions and Kaons.
195
196In the Genereate of this class all entries in the cocktail
197are called and we define a "primordial trigger" with requires
198a minimum number of muons above a Pt cut in the required acceptance.
199In order to normalized to the real number of simulated events,
200there are 2 data members in the class fNsuceeded adn fNGenerate
201which tell us what is the biais source.
202
203Enclose an example to use this generator:
204AliGenMUONCocktail * gener = new AliGenMUONCocktail();
205gener->SetPtRange(1.,100.); // Transverse momentum range
206gener->SetPhiRange(0.,360.); // Azimuthal angle range
207gener->SetYRange(-4.0,-2.4);
208gener->SetMuonPtCut(1.);
209gener->SetMuonThetaCut(171.,178.);
210gener->SetMuonMultiplicity(2);
211gener->SetNumberOfCollisions(1950.);
212gener->SetNumberOfParticipants(400.);
213gener->SetVertexSmear(kPerTrack);
214gener->SetOrigin(0,0,0); // Vertex position
215gener->SetSigma(0,0,0.0); // Sigma in (X,Y,Z) (cm) on IP position
216gener->Init();
217
a88eb0d0 218
88cb7938 219===========================================================
220 Still working ..............
221===========================================================