macro for PMD reconstruction

diff --git a/PMD/AliPMDClusterFinder.cxx b/PMD/AliPMDClusterFinder.cxx
index cfea1a907546153e410fd9be3b61a89129f8bf83..75b1a38c6a09ce85f1ed438c6cd2c0a7f75cf10b 100644 (file)
--- 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 (file)
index 0000000..ec56511
--- /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 (file)
index 0000000..2205a60
--- /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();
+
+    }
+
+
+}