macros

diff --git a/PMD/AliPMDDigits2Recpoints.C b/PMD/AliPMDDigits2Recpoints.C
new file mode 100644 (file)
index 0000000..366a26c
--- /dev/null
+++ b/PMD/AliPMDDigits2Recpoints.C
@@ -0,0 +1,42 @@
+// ----------------------------------------------------//
+//                                                     //
+//    This macro does Digits to Reconstructed Points   //
+//                                                     //
+// ----------------------------------------------------//
+
+#include "Riostream.h"
+#include "TROOT.h"
+#include "TFile.h"
+#include "TNetFile.h"
+#include "TRandom.h"
+#include "TTree.h"
+#include "TBranch.h"
+#include "TClonesArray.h"
+#include "TStopwatch.h"
+#include <stdlib.h>
+
+void AliPMDDigits2Recpoints(Int_t nevt=1) 
+{
+  TStopwatch timer;
+  timer.Start();
+
+  // Input file name
+  Char_t *alifile = "galice.root"; 
+
+  // Create the PMD Cluster Finder 
+  AliPMDClusterFinder *clus = new AliPMDClusterFinder();
+  
+  // Open the AliRoot file
+  clus->OpengAliceFile(alifile,"DR");
+  
+  Int_t ievt;
+  
+  for (ievt = 0; ievt < nevt; ievt++)
+    {
+      clus->Digits2RecPoints(ievt);
+    }
+  clus->UnLoad("R");
+  timer.Stop();
+  timer.Print();
+}
+

diff --git a/PMD/AliPMDHits2Digits.C b/PMD/AliPMDHits2Digits.C
new file mode 100644 (file)
index 0000000..11ecf92
--- /dev/null
+++ b/PMD/AliPMDHits2Digits.C
@@ -0,0 +1,36 @@
+// ----------------------------------------------------//
+//                                                     //
+//       This macro does Hits to Digits                //
+//                                                     //
+// ----------------------------------------------------//
+
+#include <stdlib.h>
+#include "Riostream.h"
+#include "TROOT.h"
+#include "TFile.h"
+#include "TNetFile.h"
+#include "TRandom.h"
+#include "TTree.h"
+#include "TBranch.h"
+#include "TClonesArray.h"
+#include "TStopwatch.h"
+
+void AliPMDHits2Digits(Int_t nevt=1) 
+{
+  // Input file name
+  Char_t *alifile = "galice.root"; 
+
+  // Create the PMD digitzer 
+  AliPMDDigitizer *digitizer = new AliPMDDigitizer();
+                                                  
+  // Open the AliRoot file
+  digitizer->OpengAliceFile(alifile,"HS");
+
+  // Create the digits
+  for (Int_t ievt=0; ievt<nevt; ievt++)
+    {
+      digitizer->Hits2Digits(ievt);
+    }
+  digitizer->UnLoad("S");
+}
+

diff --git a/PMD/AliPMDHits2SDigits.C b/PMD/AliPMDHits2SDigits.C
new file mode 100644 (file)
index 0000000..81d7bf2
--- /dev/null
+++ b/PMD/AliPMDHits2SDigits.C
@@ -0,0 +1,47 @@
+// ----------------------------------------------------//
+//                                                     //
+//       This macro does Hits to SDigits               //
+//                                                     //
+// ----------------------------------------------------//
+
+#include <stdlib.h>
+#include "Riostream.h"
+#include "TROOT.h"
+#include "TFile.h"
+#include "TNetFile.h"
+#include "TRandom.h"
+#include "TTree.h"
+#include "TBranch.h"
+#include "TClonesArray.h"
+#include "TStopwatch.h"
+
+void AliPMDHits2SDigits(Int_t nevt=1) 
+{
+  TStopwatch timer;
+  timer.Start();
+  Float_t zpos = 365.0;
+
+  // Input (and output) file name
+  Char_t *alifile = "galice.root"; 
+
+  // Create the PMD digitzer 
+  AliPMDDigitizer *digitizer = new AliPMDDigitizer();
+                                                  
+  // Open the AliRoot file
+  digitizer->OpengAliceFile(alifile,"HS");
+
+  digitizer->SetZPosition(zpos);
+  
+  // Create the sdigits
+
+  for (Int_t ievt=0; ievt<nevt; ievt++)
+    {
+      digitizer->Hits2SDigits(ievt);
+    }
+
+  digitizer->UnLoad("S");
+  timer.Stop();
+  timer.Print();
+  
+}
+

diff --git a/PMD/AliPMDSDigits2Digits.C b/PMD/AliPMDSDigits2Digits.C
new file mode 100644 (file)
index 0000000..f4edbe1
--- /dev/null
+++ b/PMD/AliPMDSDigits2Digits.C
@@ -0,0 +1,34 @@
+#include <stdlib.h>
+
+#include "Riostream.h"
+#include "TROOT.h"
+#include "TFile.h"
+#include "TNetFile.h"
+#include "TRandom.h"
+#include "TTree.h"
+#include "TBranch.h"
+#include "TClonesArray.h"
+#include "TStopwatch.h"
+
+
+void AliPMDSDigits2Digits(Int_t nevt=1) 
+{
+  // Input (and output) file name
+  Char_t *alifile = "galice.root"; 
+
+  // Create the PMD digitzer 
+  AliPMDDigitizer *digitizer = new AliPMDDigitizer();
+                                                  
+  // Open the AliRoot file
+  digitizer->OpengAliceFile(alifile,"SD");
+
+  Int_t ievt;
+
+  for (ievt = 0; ievt < nevt; ievt++)
+    {
+      digitizer->SDigits2Digits(ievt);
+    }
+  digitizer->UnLoad("S");
+
+}
+