Macro to perform reconstruction in the ZDC

diff --git a/ZDC/ZDCDigits2Reco.C b/ZDC/ZDCDigits2Reco.C
new file mode 100644 (file)
index 0000000..cdc1025
--- /dev/null
+++ b/ZDC/ZDCDigits2Reco.C
@@ -0,0 +1,42 @@
+//             Macro to perform ZDC reconstruction
+void ZDCDigits2Reco(Int_t nev=1) 
+{
+   delete gAlice;
+   gAlice=0;
+// Dynamically link some shared libs
+   if (gClassTable->GetID("AliRun") < 0) {
+      gROOT->LoadMacro("loadlibs.C");
+      loadlibs();
+   }
+      
+// Connect the Root Galice file containing Geometry, Kine, Hits and Digits
+   TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
+    if (!file) {
+       printf("\n Creating galice.root \n");
+       file = new TFile("galice.root","UPDATE");
+    } else {
+       printf("\n galice.root found in file list");
+    }
+
+// Get AliRun object from file or create it if not on file
+   if (!gAlice) {
+      gAlice = (AliRun*)file->Get("gAlice");
+      if (gAlice) printf("AliRun object found on file\n");
+      if (!gAlice) {
+           printf("\n create new gAlice object");
+           gAlice = new AliRun("gAlice","Alice test program");
+       }
+   }
+   
+   AliZDC *ZDC  = (AliZDC*) gAlice->GetModule("ZDC");
+   AliZDCMerger *merger = new AliZDCMerger();
+   merger->SetMode(1);
+   merger->SetBackgroundFileName("galice.root");
+   ZDC->SetMerger(merger);
+   
+//   Loop over events to be reconstructed          
+   for(Int_t iev=0; iev<nev; iev++) {
+       merger->SetBackgroundEventNum(iev);
+       gAlice->Digits2Reco("ZDC");
+   }   // event loop 
+}