MUON/README

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