Definition of TargetSi detector

libVAMOS/inc/TargetSi.hh

+#ifndef TARGETSI_H
+#define TARGETSI_H
+
+#include "BaseDetector.hh"
+#include "MesytecParameters.hh"
+
+class TargetSi : public BaseDetector
+{
+public:
+    //! Constructor
+    TargetSi(const Char_t * Name,
+         UShort_t NDetectors,
+                   Bool_t RawData,
+                   Bool_t CalData,
+                   Bool_t DefaultCalibration);
+     ~TargetSi(void); //!< Destructor
+
+    void AllocateComponents(void);
+
+    //! Treat the data
+    Bool_t Treat(void);
+
+    //! SetParameters
+    void SetParameters(Parameters* Par,Map* Map);
+
+
+#ifdef WITH_ROOT
+    void SetOpt(TTree *OutTTree, TTree *InTTree);
+#endif
+
+private:
+
+    BaseDetector* Si_front;
+    BaseDetector* Si_back;
+};
+
+#endif // TARGETSI_H

libVAMOS/src/TargetSi.cc

+#include "TargetSi.hh"
+
+TargetSi::TargetSi(const Char_t * Name,
+           UShort_t NDetectors,
+           Bool_t RawData,
+           Bool_t CalData,
+           Bool_t DefaultCalibration)
+    : BaseDetector(Name,0,false,CalData,true,false,"")
+{
+  START;
+  NumberOfSubDetectors = NDetectors;
+  fRawDataSubDetectors  = RawData;
+  AllocateComponents();
+  Si_front = nullptr;
+  Si_back = nullptr;
+  if(CalData)
+  {
+//
+  }
+  END;
+}
+
+
+TargetSi::~TargetSi()
+{
+  START;
+
+  END;
+}
+
+
+
+void TargetSi::AllocateComponents()
+{
+  START;
+  Char_t Name[200];
+  Bool_t RawData = fRawDataSubDetectors;
+  Bool_t CalData = fCalData;
+  Bool_t DefCal = true;
+  Bool_t PosInfo = false;
+
+  // Front E strips
+  sprintf(Name,"%s_Front_E",DetectorName);
+  Si_front = new BaseDetector(Name,NumberOfSubDetectors,RawData,CalData,DefCal,PosInfo,"",false,true);
+  if(RawData)
+  {
+    for (int i=0; i < NumberOfSubDetectors; i++)
+    {
+      sprintf(Name,"%s_Front_%02d_E",DetectorName,i);
+      Si_front->SetParameterName(Name, i);
+    }
+    Si_front->SetGateRaw(0,65534);
+    Si_front->SetRawHistogramsParams(16384,0, 65533,"");
+  }
+  if(CalData)
+  {
+    Si_front->HasCalCharges(false);
+    Si_front->SetCalHistogramsParams(1000,0,1000,"MeV");
+  }
+  AddComponent(Si_front);
+
+  // Back E strips
+  sprintf(Name,"%s_Back_E",DetectorName);
+  Si_back = new BaseDetector(Name,1,RawData,CalData,DefCal,PosInfo,"",false,true);
+  if(RawData)
+  {
+    Si_back->SetParameterName(Name, 0);
+    
+    Si_back->SetGateRaw(0,65534);
+    Si_back->SetRawHistogramsParams(16384,0, 65533,"");
+  }
+  if(CalData)
+  {
+    Si_back->HasCalCharges(false);
+    Si_back->SetCalHistogramsParams(1000,0,1000,"MeV");
+  }
+  AddComponent(Si_back);
+
+  END;
+}
+void TargetSi::SetParameters(Parameters* Par,Map* Map)
+{
+  START;
+  MesytecParameters *MP = MesytecParameters::getInstance();
+  if(isComposite)
+  {
+    for(UShort_t i=0; i< DetList->size(); i++)
+    {
+      DetList->at(i)->SetParametersMesytec(MP,Map);
+    }
+
+  }
+  END;
+}
+
+Bool_t TargetSi::Treat()
+{
+  START;
+
+  Ctr->at(0)++;
+
+
+
+  if(isComposite)
+  {
+    for(UShort_t i=0; i< DetList->size(); i++)
+    {
+
+        DetList->at(i)->Treat();
+    }
+  }
+
+  return IsPresent();
+  END;
+}
+#ifdef WITH_ROOT
+
+void TargetSi::SetOpt(TTree *OutTTree, TTree *InTTree)
+{
+  START;
+
+  // Set histogram Hierarchy
+  SetMainHistogramFolder("TargetSi");
+
+  SetHistogramsCal1DFolder("");
+  SetHistogramsCal2DFolder("");
+  for(UShort_t i = 0;i<DetList->size();i++)
+  {
+    
+    DetList->at(i)->SetMainHistogramFolder("");
+    DetList->at(i)->SetHistogramsRaw1DFolder("TargetSi_R");
+    DetList->at(i)->SetHistogramsRaw2DFolder("TargetSi_R");
+    DetList->at(i)->SetHistogramsCal1DFolder("TargetSi_C");
+    DetList->at(i)->SetHistogramsCal2DFolder("TargetSi_C");
+  }
+
+  if(fMode == MODE_WATCHER)
+  {
+    if(fRawData)
+    {
+      SetHistogramsRaw(true);
+    }
+    if(fCalData)
+    {
+      SetHistogramsCal(true);
+    }
+    for(UShort_t i = 0;i<DetList->size();i++)
+    {
+      if(DetList->at(i)->HasRawData())
+      {
+
+          DetList->at(i)->SetHistogramsRaw1D(true);
+          DetList->at(i)->SetHistogramsRaw2D(true);
+          DetList->at(i)->SetHistogramsRawSummary(true);
+          DetList->at(i)->SetOutAttachRawV(false);
+          DetList->at(i)->SetOutAttachRawI(false);
+      }
+      if(DetList->at(i)->HasCalData())
+      {
+        DetList->at(i)->SetHistogramsCal1D(true);
+        DetList->at(i)->SetHistogramsCal2D(true);
+        DetList->at(i)->SetHistogramsCalSummary(true);
+      }
+    }
+  }
+  else if(fMode == MODE_D2R)
+  {
+    for(UShort_t i = 0;i<DetList->size();i++)
+    {
+      if(DetList->at(i)->HasRawData())
+      {
+
+          DetList->at(i)->SetHistogramsRaw1D(false);
+          DetList->at(i)->SetHistogramsRaw2D(true);
+          DetList->at(i)->SetHistogramsRawSummary(true);
+          DetList->at(i)->SetOutAttachRawV(true);
+          DetList->at(i)->SetOutAttachRawF(false);
+          DetList->at(i)->SetOutAttachRawI(false);
+	  DetList->at(i)->SetOutAttachTS(true);
+      }
+      DetList->at(i)->OutAttach(OutTTree);
+    }
+    OutAttach(OutTTree);
+
+  }
+  else if(fMode == MODE_D2A)
+  {
+    if(fRawData)
+    {
+      SetHistogramsRaw(false);
+      SetOutAttachRawI(false);
+      SetOutAttachRawV(false);
+    }
+    if(fCalData)
+    {
+      SetHistogramsCal(true);
+      SetOutAttachCalI(false);
+      SetOutAttachCalV(false);
+    }
+
+    for(UShort_t i = 0;i<DetList->size();i++)
+    {
+      if(DetList->at(i)->HasRawData())
+      {
+
+          DetList->at(i)->SetHistogramsRaw1D(false);
+          DetList->at(i)->SetHistogramsRaw2D(true);
+          DetList->at(i)->SetHistogramsRawSummary(true);
+          DetList->at(i)->SetOutAttachRawV(true);
+          DetList->at(i)->SetOutAttachRawI(false);
+          DetList->at(i)->SetOutAttachTS(true);
+
+      }
+      if(DetList->at(i)->HasCalData())
+      {
+
+          DetList->at(i)->SetHistogramsCal1D(false);
+          DetList->at(i)->SetHistogramsCal2D(true);
+          DetList->at(i)->SetHistogramsCalSummary(true);
+          DetList->at(i)->SetOutAttachCalV(true);
+          DetList->at(i)->SetOutAttachCalI(false);
+      }
+      DetList->at(i)->OutAttach(OutTTree);
+    }
+
+    OutAttach(OutTTree);
+  }
+  else if(fMode == MODE_R2A)
+  {
+
+    SetNoOutput();
+
+  }
+  else if(fMode == MODE_RECAL)
+  {
+
+    SetNoOutput();
+
+  }
+  else if (fMode == MODE_CALC)
+  {
+    SetNoOutput();
+  }
+  else
+  {
+    Char_t Message[500];
+    sprintf(Message,"In <%s><%s> Trying to set the detector unknown Mode (%d) !", GetName(), GetName(),fMode );
+    MErr * Er= new MErr(WhoamI,0,0, Message);
+    throw Er;
+  }
+  END;
+}
+
+
+
+#endif