Macros for Simulation, Reconstruction and ESD analysis
authorschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 4 Feb 2004 17:57:20 +0000 (17:57 +0000)
committerschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 4 Feb 2004 17:57:20 +0000 (17:57 +0000)
PHOS/AnaESD.C [new file with mode: 0644]
PHOS/Reconstruction.C [new file with mode: 0644]
PHOS/Simulation.C [new file with mode: 0644]

diff --git a/PHOS/AnaESD.C b/PHOS/AnaESD.C
new file mode 100644 (file)
index 0000000..954fa95
--- /dev/null
@@ -0,0 +1,17 @@
+#include "TFile.h"
+#include "AliPHOSGetter.h"
+#include "Riostream.h"
+#include "AliESD.h"
+
+void Ana() 
+{
+  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(); 
+    
+  }
+}
diff --git a/PHOS/Reconstruction.C b/PHOS/Reconstruction.C
new file mode 100644 (file)
index 0000000..ddec37d
--- /dev/null
@@ -0,0 +1,44 @@
+#include "AliReconstruction.h"
+#include "TString.h"
+#include "Riostream.h"
+#include "AliPHOSGetter.h"
+#include "AliEMCALGetter.h"
+
+void reco(TString opt, TString name) 
+{
+  AliReconstruction rec ; 
+  if ( !opt.Contains("T") )
+    rec.SetRunTracking(kFALSE) ;
+  if ( opt.Contains("R") ) 
+    rec.SetRunReconstruction(name.Data()) ; 
+  if ( !opt.Contains("E") )
+    rec.SetFillESD("") ; 
+  else 
+    rec.SetFillESD(name.Data()) ; 
+  rec.Run() ;
+
+  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, "RP") ; 
+      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") ) {
+    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, "RP") ; 
+      cout << "       RecPoints  # " << gime->ECARecPoints()->GetEntries() << endl ; 
+      cout << "   Rec Particles  # " << gime->RecParticles()->GetEntries() << endl ; 
+    }
+ } 
+}   
diff --git a/PHOS/Simulation.C b/PHOS/Simulation.C
new file mode 100644 (file)
index 0000000..2a6932d
--- /dev/null
@@ -0,0 +1,49 @@
+#include "AliSimulation.h"
+#include "TString.h"
+#include "AliPHOSGetter.h"
+#include "AliEMCALGetter.h"
+#include "Riostream.h"
+
+void simu(TString opt, TString name) 
+{
+  AliSimulation sim ; 
+  // Generation and simulation
+  if ( !opt.Contains("G") )
+    sim.SetRunGeneration(kFALSE) ;
+  // Making SDigits 
+  if ( !opt.Contains("S") )
+    sim.SetMakeSDigits("") ; 
+  else 
+    sim.SetMakeSDigits(name.Data()) ;
+  // Making Digits 
+  if ( !opt.Contains("D") )
+    sim.SetMakeDigits("") ; 
+  else 
+    sim.SetMakeDigits(name.Data()) ;    
+  //Merging
+  // to implement 
+  sim.Run() ;  
+  // 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 ; 
+    }
+  }
+}