]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/README
Move to AOD/ESD schema - ITS readers not needed any more
[u/mrichter/AliRoot.git] / MUON / README
index fc7d66a2c621e6b8dd7e179814db9ff6d92838f9..782858967f4a04bf04d9d9a2d51b0827edee49ea 100644 (file)
@@ -35,7 +35,30 @@ for hits, you have to do the following before:
 root [0] gAlice->SetConfigFunction("Config( \"/home/martinez/aliroot/work_NewIO/test/\" , \"box\" );"); 
 
 ============================================================
- How to run MUON Digitalization
+ How to run MUON Hits to SDigits
+============================================================
+AliRunDigitizer   * manager = new AliRunDigitizer(1,1);
+manager->SetInputStream(0,"galice.root");
+AliMUONDigitizer * dMUON   = new AliMUONSDigitizerv1(manager);
+manager->AddDigitizer(dMUON);
+manager->Exec("deb");
+
+
+============================================================
+ How to run MUON Digitalization from SDigits
+============================================================
+To process the digitalization of hits the AliRunDigitizer 
+framework is used:
+
+AliRunDigitizer   * manager = new AliRunDigitizer(1,1);
+manager->SetInputStream(0,"galice.root");
+AliMUONDigitizerv2* dMUON   = new AliMUONDigitizerv2(manager)
+manager->AddDigitizer(dMUON);
+manager->Exec("deb");
+
+
+============================================================
+ How to run MUON Digitalization from Hits
 ============================================================
 To process the digitalization of hits the AliRunDigitizer 
 framework is used:
@@ -66,6 +89,7 @@ gSystem->Load("$ALICE_ROOT/MUON/MUONrawclusters_C.so")
 MUONtrigger("galice.root",0,9999); 
 To Load
 gSystem->Load("$ALICE_ROOT/MUON/MUONtrigger_C.so")
+
 ============================================================
  How to run MUON Tracking reconstruction (old output)
 ===========================================================
@@ -89,6 +113,33 @@ gSystem->Load("$ALICE_ROOT/MUON/MUONTracker_C.so")
 Output in MUON.Tracks.root using the new IO 
 as a TClonesArray of AliMUONTrack
 
+Note1.In the actual version (dec-03) the following data 
+members  of ALiMUONTrack:  fgFitter, fEventReconstructor, 
+and fTrackHitsPtr are not saved in the MUON.Track.root file
+
+
+===========================================================
+ How to run MUON Trigger Track reconstruction (new output)
+===========================================================
+.includepath $ALICE_ROOT/STEER
+.includepath $ALICE_ROOT/MUON
+.L $ALICE_ROOT/MUON/MUONTriggerTracker.C++
+MUONTriggerTracker("galice.root",0,10);    
+To Load
+gSystem->Load("$ALICE_ROOT/MUON/MUONTriggerTracker_C.so")
+Output in MUON.Tracks.root using the new IO 
+as a TClonesArray of AliMUONTriggerTrack
+
+==========================================================
+ How to run full MUON  reconstruction 
+(clusterisation + trigger + tracking)
+===========================================================
+.includepath $ALICE_ROOT/STEER
+.includepath $ALICE_ROOT/MUON
+.L $ALICE_ROOT/MUON/MUONChallengeTest.C++
+MUONRecoTest("galice.root");
+
+
 ============================================================
  How to run MUONCheck macro
 ============================================================
@@ -103,19 +154,21 @@ To Load
 gSystem->Load("$ALICE_ROOT/MUON/MUONCheck_C.so")
 
 To print Kine : (default file is galice.root )
-MUONkine()
+MUONkine() or MUONkine("galice.root",##) for the event number ##
 
 To print hits : (default file is galice.root if not MUONhits("toto.root""); )
-MUONhits()  
+MUONhits()  or MUONhits("galice.root",##) for the event number ##
 
 To print digits : (default file is galice.root)
-MUONdigits()  
+MUONdigits()  or MUONdigits("galice.root",##) for the event number ##
 
 To print rawcluster : (default file is galice.root)
-MUONrecpoints()
+MUONrecpoints() or MUONrecpoints("galice.root",##) for the event number ##
 
 To print trigger : (default file is galice.root)
-MUONTestTrigger()
+MUONTestTrigger() or MUONTestTrigger("galice.root",##) for the event number ##
+
+....
 
 
 ============================================================
@@ -129,7 +182,7 @@ G3GUI()
 ============================================================
  How to run MUONdisplay
 ============================================================
-Frist you need to perfomr a full simulation: 
+First you need to perform a full simulation: 
 generation, digitalisation and clusterisation
 .L $ALICE_ROOT/MUON/MUONdisplay.C
 MUONdisplay(0,"galice.root")
@@ -163,6 +216,72 @@ TRACKING must be 1 for RKUTA and 2 for HELIX (the default value for aliroot is 2
 FACTOR allows you to set the magnetic field to 0, just putting FACTOR=0. Default value is 1.
 MAXB is the maximum magnetic field which is 10.T
 
+===========================================================
+ MUON cocktail for physics ..............
+===========================================================
+There is a MUON cocktail generator of the muon sources in the
+EVGEN directory. This class derives from AliGenCocktail.
+In the init of this class I have filled the cocktail with 
+the muon sources: J/Psi, Upsilon, Open Charm, Open Beauty, 
+Pion, Kaons. The code needs only the production cross section 
+at 4pi (for the moment this values are in the code since I 
+prefere them do not be modified), and the code calculates the  
+rate of particles in the acceptance, making the scaling based 
+on the number of collisions for the hard probes and on the  
+number of participants for soft sources: Pions and Kaons.
+
+In the Genereate of this class all entries in the cocktail 
+are called and we define a "primordial trigger" with requires 
+a minimum number of muons above a Pt cut in the required acceptance.
+In order to normalized to the real number of simulated events, 
+there are 2 data members in the class fNsuceeded adn fNGenerate 
+which tell us what is the biais source.
+
+Enclose an example to use this generator:   
+AliGenMUONCocktail * gener = new AliGenMUONCocktail();
+gener->SetPtRange(1.,100.);       // Transverse momentum range  
+gener->SetPhiRange(0.,360.);    // Azimuthal angle range 
+gener->SetYRange(-4.0,-2.4);
+gener->SetMuonPtCut(1.);
+gener->SetMuonThetaCut(171.,178.);
+gener->SetMuonMultiplicity(2);
+gener->SetNumberOfCollisions(1950.);   
+gener->SetNumberOfParticipants(400.);     
+gener->SetVertexSmear(kPerTrack);  
+gener->SetOrigin(0,0,0);        // Vertex position
+gener->SetSigma(0,0,0.0);       // Sigma in (X,Y,Z) (cm) on IP position
+gener->Init();
+===========================================================
+ csh Script for the full reconstruction
+===========================================================
+
+
+aliroot -b << EOF  
+AliSimulation MuonSim("YourConfig.C")
+MuonSim.Run(XXX)
+.q
+EOF
+
+aliroot -b << EOF 
+TPluginManager* pluginManager = gROOT->GetPluginManager();
+pluginManager->AddHandler("AliReconstructor", "MUON","AliMUONReconstructor", "MUON","AliMUONReconstructor()")
+AliReconstruction MuonRec("galice.root") 
+MuonRec.SetRunVertexFinder(kFALSE)
+MuonRec.SetRunLocalReconstruction("MUON")
+MuonRec.SetFillESD("MUON")
+MuonRec.Run() 
+.q
+EOF
+
+aliroot -b << EOF  
+.includepath $ALICE_ROOT/STEER
+.includepath $ALICE_ROOT/MUON
+.L $ALICE_ROOT/MUON/MUONmassPlot_ESD.C++
+MUONmassPlot("galice.root",0,99999);
+.q
+EOF
+
 ===========================================================
  Still working ..............
 ===========================================================