updating
authorschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 5 Feb 2004 15:51:40 +0000 (15:51 +0000)
committerschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 5 Feb 2004 15:51:40 +0000 (15:51 +0000)
PHOS/Checker.C [new file with mode: 0644]
PHOS/Reconstruction.C
PHOS/Simulation.C

diff --git a/PHOS/Checker.C b/PHOS/Checker.C
new file mode 100644 (file)
index 0000000..e32a417
--- /dev/null
@@ -0,0 +1,84 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+/* $Id$ */
+//_________________________________________________________________________
+// Macros checking what has been produced by Simulation.C and Reconstruction.C 
+// Use Case : 
+//          root> .L Checking.C++
+//          root> check("HSDRE", "PHOS EMCAL") --> checks how many hits, Sdigits, Digits and 
+//                                               Reconstruction objects and ESD have been produced
+// author  : Yves Schutz (CERN/SUBATECH)
+// February 2004
+//_________________________________________________________________________
+
+#include "AliPHOSGetter.h"
+#include "AliEMCALGetter.h"
+#include "AliESD.h" 
+#include "Riostream.h"
+#include "TString.h" 
+
+void check(TString opt="HSDRE", TString name="all") 
+{
+  if ( name.Contains("PHOS") ) {
+    AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ; 
+    cout << ">>>>>>>>>>>> PHOS:" ; 
+    Int_t event ; 
+    for (event = 0; event < gime->MaxEvent(); event++) {
+      gime->Event(event) ;
+      cout << " event # " << event << endl ; 
+      if ( opt.Contains("H") )
+       cout <<  "  Hits          # " << gime->Hits()->GetEntries() << endl ; 
+      if ( opt.Contains("S") )
+       cout <<  "  SDigits       # " << gime->SDigits()->GetEntries() << endl ;     
+      if ( opt.Contains("D") )
+       cout <<  "  Digits        # " << gime->Digits()->GetEntries() << endl ;     
+      if ( opt.Contains("R") ){
+       cout <<  "  EMC Recpoints # " << gime->EmcRecPoints()->GetEntries() << endl ;     
+       cout <<  "  CPV Recpoints # " << gime->CpvRecPoints()->GetEntries() << endl ;     
+       cout <<  "  Track Segments# " << gime->TrackSegments()->GetEntries() << endl ;     
+       cout <<  "  Rec Particles # " << gime->RecParticles()->GetEntries() << endl ;     
+      }
+    }
+  }
+  if ( name.Contains("EMCAL") ) {
+    AliEMCALGetter * gime = AliEMCALGetter::Instance("galice.root") ; 
+    cout <<  ">>>>>>>>>>>> EMCAL:" ; 
+    Int_t event ; 
+    for (event = 0; event < gime->MaxEvent(); event++) {
+      gime->Event(event) ;
+      cout << " event # " << event << endl ; 
+      if ( opt.Contains("H") )
+       cout <<  "  Hits          # " << gime->Hits()->GetEntries() << endl ; 
+      if ( opt.Contains("S") )
+       cout <<  "  SDigits       # " << gime->SDigits()->GetEntries() << endl ;     
+      if ( opt.Contains("D") )
+       cout <<  "  Digits        # " << gime->Digits()->GetEntries() << endl ;     
+      if ( opt.Contains("R") ){
+       cout <<  "  EMC Recpoints # " << gime->ECARecPoints()->GetEntries() << endl ;     
+       cout <<  "  Rec Particles # " << gime->RecParticles()->GetEntries() << endl ;     
+      }
+    }
+  }
+  if ( opt.Contains("E") ) {
+    AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ; 
+    Int_t nEvent = gime->MaxEvent() ;  
+    Int_t event ; 
+    AliESD * esd ;
+    for (event = 0 ; event < nEvent; event++) {
+      esd = gime->ESD(event) ; 
+      esd->Print(); 
+    }
+  }
+}
index 21212ff..6d7ab59 100644 (file)
@@ -38,6 +38,8 @@ void reco(TString opt="TRE", TString name="all")
     rec.SetRunTracking(kFALSE) ;
   if ( opt.Contains("R") ) 
     rec.SetRunReconstruction(name.Data()) ; 
+  else 
+    rec.SetRunReconstruction("") ;
   if ( !opt.Contains("E") )
     rec.SetFillESD("") ; 
   else 
index 6a526a3..429a686 100644 (file)
@@ -49,31 +49,34 @@ void simu(Int_t nevents=1, TString opt="GSD", TString name="all")
   else 
     sim.SetMakeDigits(name.Data()) ;    
   //Merging
-  if ( !opt.Contains("M") )
+  if ( opt.Contains("M") )
     sim.MergeWith("./bgrd/galice.root", nevents) ;  
   // to implement 
   sim.Run(nevents) ;  
-  // Checking result
-  if ( name.Contains("PHOS") ) {
-    cout << ">>>>>>>>>>>> PHOS " << endl ; 
-    AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;
-    Int_t event ; 
-    for (event = 0; event < gime->MaxEvent(); event++) {
-      cout << "event # " << event << endl ; 
-      gime->Event(event, "SD") ; 
-      cout << "  SDigits # " << gime->SDigits()->GetEntries() << endl ; 
-      cout << "   Digits # " << gime->Digits()->GetEntries() << endl ; 
+
+  if ( !opt.Contains("M") ) {
+    // Checking result
+    if ( name.Contains("PHOS") ) {
+      cout << ">>>>>>>>>>>> PHOS " << endl ; 
+      AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;
+      Int_t event ; 
+      for (event = 0; event < gime->MaxEvent(); event++) {
+       cout << "event # " << event << endl ; 
+       gime->Event(event, "SD") ; 
+       cout << "  SDigits # " << gime->SDigits()->GetEntries() << endl ; 
+       cout << "   Digits # " << gime->Digits()->GetEntries() << endl ; 
+      }
     }
-  }
-  if ( name.Contains("EMCAL") ) {
-    cout << ">>>>>>>>>>>> EMCAL " << endl ; 
-    AliEMCALGetter * gime = AliEMCALGetter::Instance("galice.root") ;
-    Int_t event ; 
-    for (event = 0; event < gime->MaxEvent(); event++) {
-      cout << "event # " << event << endl ; 
-      gime->Event(event, "SD") ; 
-      cout << "  SDigits # " << gime->SDigits()->GetEntries() << endl ; 
-      cout << "   Digits # " << gime->Digits()->GetEntries() << endl ; 
+    if ( name.Contains("EMCAL") ) {
+      cout << ">>>>>>>>>>>> EMCAL " << endl ; 
+      AliEMCALGetter * gime = AliEMCALGetter::Instance("galice.root"); 
+      Int_t event ; 
+      for (event = 0; event < gime->MaxEvent(); event++) {
+       cout << "event # " << event << endl ; 
+       gime->Event(event, "SD") ; 
+       cout << "  SDigits # " << gime->SDigits()->GetEntries() << endl ; 
+       cout << "   Digits # " << gime->Digits()->GetEntries() << endl ; 
+      }
     }
   }
 }