Removed method

[u/mrichter/AliRoot.git] / MUON / MUONreco.C

diff --git a/MUON/MUONreco.C b/MUON/MUONreco.C
index eff1164fe37ca7efc777d81b7740ac4e555464b3..99395f18b088f3e729078b90bcf263e05938bbed 100644 (file)
--- a/MUON/MUONreco.C
+++ b/MUON/MUONreco.C
@@ -1,4 +1,4 @@
-// Macro for testing the new C++ reconstruction code
+// Macro MUONreco.C for testing the C++ reconstruction code
 
 // Arguments:
 //   FirstEvent (default 0)
@@ -7,10 +7,11 @@
 //   FileName (for signal) (default "galice.root")
 //   BkgGeantFileName (for background),
 //      needed only if RecGeantHits = 1 and background to be added
+#include <iostream.h>
 void MUONreco (Int_t FirstEvent = 0, Int_t LastEvent = 0, Int_t RecGeantHits = 0, Text_t *FileName = "galice.root", Text_t *BkgGeantFileName = "")
 {
   //
-  cout << "MUONtest_reco" << endl;
+  cout << "MUONreco" << endl;
   cout << "FirstEvent " << FirstEvent << endl;
   cout << "LastEvent " << LastEvent << endl;
   cout << "RecGeantHits " << RecGeantHits << endl;
@@ -22,7 +23,8 @@ void MUONreco (Int_t FirstEvent = 0, Int_t LastEvent = 0, Int_t RecGeantHits = 0
     loadlibs();
   }
 
-  // Connect the Root Galice file containing Geometry, Kine and Hits
+  // Connect the Root Galice file containing Geometry, Kine, Hits
+  // and eventually RawClusters
   TFile *file = (TFile*) gROOT->GetListOfFiles()->FindObject(FileName);
   if (!file) {
     printf("\n Creating file %s\n", FileName);
@@ -43,19 +45,23 @@ void MUONreco (Int_t FirstEvent = 0, Int_t LastEvent = 0, Int_t RecGeantHits = 0
   // Initializations
   // AliMUON *MUON  = (AliMUON*) gAlice->GetModule("MUON"); // necessary ????
   AliMUONEventReconstructor *Reco = new AliMUONEventReconstructor();
+  //Reco->SetTrackMethod(2); //AZ - use Kalman
   Reco->SetRecGeantHits(RecGeantHits);
 
   // The right place for changing AliMUONEventReconstructor parameters
   // with respect to the default ones
 //   Reco->SetMaxSigma2Distance(100.0);
-  Reco->SetPrintLevel(10);
+  //Reco->SetPrintLevel(10);
+  //Reco->SetPrintLevel(2);
+  //Reco-> SetChamberThicknessInX0(0.05);
+  //Reco->SetEfficiency(1.); //0.95);
   cout << "AliMUONEventReconstructor: actual parameters" << endl;
   Reco->Dump();
 
   // Loop over events
   for (Int_t event = FirstEvent; event <= LastEvent; event++) {
     cout << "Event: " << event << endl;
-    AliMUON *MUON  = (AliMUON*) gAlice->GetModule("MUON"); // necessary ????
+    //    AliMUON *MUON  = (AliMUON*) gAlice->GetModule("MUON"); // necessary ????
     Int_t nparticles = gAlice->GetEvent(event);
     cout << "nparticles: " << nparticles << endl;
     // prepare background file and/or event if necessary
@@ -64,6 +70,8 @@ void MUONreco (Int_t FirstEvent = 0, Int_t LastEvent = 0, Int_t RecGeantHits = 0
       if (Reco->GetBkgGeantFile())Reco->NextBkgGeantEvent();
     }
     Reco->EventReconstruct();
+    // Write this event in a tree
+    Reco->FillEvent();
     // Dump current event
     Reco->EventDump();
   } // Event loop