From 1f86361bf8c673369c0d008b6b2213371bb8f0b7 Mon Sep 17 00:00:00 2001
From: bnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Fri, 17 Nov 2006 07:34:05 +0000
Subject: [PATCH] macro for PMD reconstruction

---
 PMD/AliPMDClusterFinder.cxx |  3 +-
 PMD/AliPMDRec.C             | 23 +++++++++++++++
 PMD/AliPMDRootDataRead.C    | 58 +++++++++++++++++++++++++++++++++++++
 3 files changed, 82 insertions(+), 2 deletions(-)
 create mode 100644 PMD/AliPMDRec.C
 create mode 100644 PMD/AliPMDRootDataRead.C

diff --git a/PMD/AliPMDClusterFinder.cxx b/PMD/AliPMDClusterFinder.cxx
index cfea1a90754..75b1a38c6a0 100644
--- a/PMD/AliPMDClusterFinder.cxx
+++ b/PMD/AliPMDClusterFinder.cxx
@@ -571,8 +571,7 @@ void AliPMDClusterFinder::Digits2RecPoints(Int_t ievt, AliRawReader *rawReader)
 
 	  Float_t sig1 = (Float_t) sig;
 	  // CALIBRATION
-	  //Float_t gain = fCalibData->GetGainFact(det,smn,row,col);
-	  Float_t gain = 1.0;
+	  Float_t gain = fCalibData->GetGainFact(det,smn,row,col);
 
 	  //printf("sig = %d gain = %f\n",sig,gain);
 	  sig = (Int_t) (sig1*gain);
diff --git a/PMD/AliPMDRec.C b/PMD/AliPMDRec.C
new file mode 100644
index 00000000000..ec565116845
--- /dev/null
+++ b/PMD/AliPMDRec.C
@@ -0,0 +1,23 @@
+void AliPMDRec()
+{
+  // This macro for the full reconstruction chain. Only PMD is ON.
+  //
+
+  Int_t firstEvent = 0;
+  Int_t lastEvent = 1;
+
+  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
+  AliTracker::SetFieldMap(field,1);
+
+  AliReconstruction rec;
+  rec.SetRunReconstruction("PMD");
+  rec.SetRunVertexFinder(kFALSE);
+  rec.SetFillESD("PMD");
+  rec.SetFillTriggerESD(kFALSE);
+
+  //  rec.Run("./");
+  //  rec.Run("raw.date");
+  rec.Run("raw_6b_61.root",firstEvent,lastEvent);
+
+}
+
diff --git a/PMD/AliPMDRootDataRead.C b/PMD/AliPMDRootDataRead.C
new file mode 100644
index 00000000000..2205a600323
--- /dev/null
+++ b/PMD/AliPMDRootDataRead.C
@@ -0,0 +1,58 @@
+//
+void AliPMDRootDataRead()
+{
+  // To read PMD raw root data and fetch the adc value for each cell
+
+  TObjArray pmdddlcont;
+
+  Int_t ievt = 2;
+
+  Bool_t junk;
+
+  AliRawReaderRoot reader("raw_6b_61.root",ievt);
+  //reader.NextEvent();
+  //reader.NextEvent();
+
+
+
+  /*
+  reader.ReadHeader();
+  cout << "LDC ID =       " << reader.GetLDCId()       << endl;
+  cout << "Equipment ID = " << reader.GetEquipmentId() << endl;
+  cout << "Data Size =    " << reader.GetDataSize()    << endl;
+  */
+
+  AliPMDRawStream stream(&reader);
+
+  Int_t indexDDL = 0;
+
+
+  for (Int_t iddl = 0; iddl < 6; iddl++)
+    {
+      
+      reader.Select("PMD", iddl, iddl);
+
+      junk = stream.DdlData(&pmdddlcont);
+
+      Int_t ientries = pmdddlcont.GetEntries();
+      for (Int_t ient = 0; ient < ientries; ient++)
+	{
+	  AliPMDddldata *pmdddl = (AliPMDddldata*)pmdddlcont.UncheckedAt(ient);
+	  
+	  Int_t det = pmdddl->GetDetector();
+	  Int_t smn = pmdddl->GetSMN();
+	  //Int_t mcm = pmdddl->GetMCM();
+	  //Int_t chno = pmdddl->GetChannel();
+	  Int_t row = pmdddl->GetRow();
+	  Int_t col = pmdddl->GetColumn();
+	  Int_t sig = pmdddl->GetSignal();
+
+	  // cout << row << " " << col << " " << sig << endl;
+
+	}
+      pmdddlcont.Clear();
+
+    }
+
+
+}
-- 
2.43.0