KVINDRA.cpp

/***************************************************************************
$Id: KVINDRA.cpp,v 1.68 2009/01/21 10:05:51 franklan Exp $
                          kvindra.cpp  -  description
                             -------------------
    begin                : Mon May 20 2002
    copyright            : (C) 2002 by J.D. Frankland
    email                : frankland@ganil.fr
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
 
#include "KVINDRA.h"
#include "KVMaterial.h"
#include "KVDetector.h"
#include "KVTelescope.h"
#include "KVIDTelescope.h"
#include "KVIDSiCsI.h"
#include "KVIDINDRACsI.h"
#include "KVIDChIoSi.h"
#include "KVIDChIoCsI.h"
#include "KVIDChIoSi75.h"
#include "KVIDSi75SiLi.h"
#include "KVIDPhoswich.h"
#include "KVLayer.h"
#include "KVRing.h"
#include "KVGroup.h"
#include "KVChIo.h"
#include "KVCsI.h"
#include "KVNucleus.h"
#include "KVDetectorEvent.h"
#include "KVNucleusEvent.h"
#include "TEnv.h"
#include "INDRAGeometryBuilder.h"
#include <KVASGroup.h>
#include "KVGeoNavigator.h"
#ifdef WITH_BUILTIN_GRU
#include <KVGANILDataReader.h>
#endif
#include <KVDBChIoPressures.h>
#include <KVGeoImport.h>
#include <KVRangeTableGeoNavigator.h>
#include <KVRawDataReader.h>
#ifdef WITH_MFM
#include "KVMFMDataFileReader.h"
#include "MFMEbyedatFrame.h"
#ifdef WITH_MESYTEC
#include "MFMMesytecMDPPFrame.h"
#include "mesytec_buffer_reader.h"
#endif
#endif
#include "KVINDRAGroupReconstructor.h"
#include <cassert>
 
using namespace std;
 
ClassImp(KVINDRA)
 
//Use this static array to translate EBaseIndra_type
//signal type to a string giving the signal type
 
 
 
Char_t KVINDRA::SignalTypes[][3] = {
   "",                          //dummy for index=0
   "GG", "PG", "T",
   "GG", "PG", "T",
   "R", "L", "T",
   "GG", "PG", "T",
   "GG", "PG", "T"
};
 
 
#define KV_DEBUG 1
 
KVINDRA* gIndra;
 
 
KVINDRA::KVINDRA()
{
   //Default constructor
   //Set up lists of ChIo, Si, CsI, Phoswich
   fChIo = new KVHashList();
   fChIo->SetCleanup(kTRUE);
   fSi = new KVHashList();
   fSi->SetCleanup(kTRUE);
   fCsI = new KVHashList();
   fCsI->SetCleanup(kTRUE);
   fPhoswich = new KVHashList();
   fPhoswich->SetCleanup(kTRUE);
   fTrigger = 0;
   gIndra = this;
   fPHDSet = kFALSE;
   fSelecteur = nullptr;
   // unlike a normal multidetector, INDRA does not own its detectors
   // they are owned by the TELESCOPE structures
   SetOwnsDetectors(kFALSE);
   fMesytecData = kFALSE;
   fEbyedatData = kFALSE;
}
 
 
 
 
KVINDRA::~KVINDRA()
{
   //Delets lists of ChIo, Si, etc.
   //Resets global gIndra pointer.
 
   if (fChIo && fChIo->TestBit(kNotDeleted)) {
      fChIo->Clear();
      delete fChIo;
   }
   fChIo = 0;
   if (fSi && fSi->TestBit(kNotDeleted)) {
      fSi->Clear();
      delete fSi;
   }
   fSi = 0;
   if (fCsI && fCsI->TestBit(kNotDeleted)) {
      fCsI->Clear();
      delete fCsI;
   }
   fCsI = 0;
   if (fPhoswich && fPhoswich->TestBit(kNotDeleted)) {
      fPhoswich->Clear();
      delete fPhoswich;
   }
   fPhoswich = 0;
   if (fSelecteur) delete fSelecteur;
   gIndra = 0;
}
 
 
 
 
 
void KVINDRA::BuildGeometry()
{
   // Construction of INDRA detector array.
   //
   // Uses infos in file
   //        $KVROOT/KVFiles/data/indra_struct.[dataset].env
   // or
   //        $KVROOT/KVFiles/data/indra_struct.env
   //
   // if no dataset-specific file found
   //
   // Alternatively, by defining the variable
   //
   //   [dataset].INDRA.StructureFile: [path to file]
 
 
   TString path = Form("indra-struct.%s.env", fDataSet.Data());
   TString path2 = GetDataSetEnv(fDataSet, "INDRA.StructureFile", "");
   if (path2 == "") {
      if (!SearchKVFile(path.Data(), path2, "data")) {
         path = "indra-struct.env";
         SearchKVFile(path.Data(), path2, "data");
      }
   }
   fStrucInfos.ReadFile(path2, kEnvAll);
 
   KVString lruns = fStrucInfos.GetValue("AddOnForRuns", "");
   //test if additional geometrical specification exists
   if (lruns != "") {
      lruns.Begin(",");
      while (!lruns.End()) {
         KVString sruns = lruns.Next();
         KVNumberList nlr(sruns.Data());
         //the current run needs specific geometry
         if (nlr.Contains(fCurrentRun)) {
            path = fStrucInfos.GetValue(sruns.Data(), "");
            Info("BuildGeometry", "Additional geometry for run=%d in file #%s#", fCurrentRun, path.Data());
            SearchKVFile(path.Data(), path2, "data");
            if (path2 == "") {
               Warning("BuildGeometry", "fichier %s inconnu", path.Data());
            }
            else {
               fStrucInfos.ReadFile(path2, kEnvChange);
            }
         }
      }
   }
 
   SetName(fStrucInfos.GetValue("INDRA.Name", ""));
   SetTitle(fStrucInfos.GetValue("INDRA.Title", ""));
 
   KVString layers = fStrucInfos.GetValue("INDRA.Layers", "");
   layers.Begin(" ");
   while (!layers.End()) BuildLayer(layers.Next());
}
 
 
 
 
void KVINDRA::BuildLayer(const Char_t* name)
{
   // Build layer 'name' with infos in file "$KVROOT/KVFiles/data/indra-struct.[dataset].env"
 
   KVLayer* layer = new KVLayer;
   Int_t number = fStrucInfos.GetValue(Form("INDRA.Layer.%s.Number", name), 0);
   layer->SetName(name);
   layer->SetNumber(number);
   Add(layer);
   Info("BuildLayer", "Building layer %d:%s", number, name);
   KVNumberList rings = fStrucInfos.GetValue(Form("INDRA.Layer.%s.Rings", name), "");
   rings.Begin();
   while (!rings.End()) {
      Int_t ring_num = rings.Next();
      TString prefix = Form("INDRA.Layer.%s.Ring.%d", name, ring_num);
      KVRing* ring = BuildRing(ring_num, prefix);
      layer->Add(ring);
   }
}
 
 
 
 
KVRing* KVINDRA::BuildRing(Int_t number, const Char_t* prefix)
{
   // Build ring with infos in file "$KVROOT/KVFiles/data/indra-struct.[dataset].env"
 
   KVRing* ring = new KVRing;
   Info("BuildRing", "Building ring %d (%s)", number, prefix);
 
   Double_t thmin = fStrucInfos.GetValue(Form("%s.ThMin", prefix), 0.0);
   Double_t thmax = fStrucInfos.GetValue(Form("%s.ThMax", prefix), 0.0);
   Double_t phi_0 = fStrucInfos.GetValue(Form("%s.Phi0", prefix), 0.0);
   Int_t ntel = fStrucInfos.GetValue(Form("%s.Ntel", prefix), 0);
   Int_t step = fStrucInfos.GetValue(Form("%s.Step", prefix), 1);
   Int_t num_first = fStrucInfos.GetValue(Form("%s.Nfirst", prefix), 1);
   Double_t dphi = fStrucInfos.GetValue(Form("%s.Dphi", prefix), -1.0);
   KVNumberList absent = fStrucInfos.GetValue(Form("%s.Absent", prefix), "");
 
   ring->SetPolarMinMax(thmin, thmax);
   KVINDRATelescope* kvt = BuildTelescope(prefix, num_first);
   Double_t phi_min = phi_0 - 0.5 * (kvt->GetAzimuthalWidth());
   phi_min = (phi_min < 0. ? phi_min + 360. : phi_min);
   dphi = (dphi < 0. ? 360. / ntel : dphi);
   Double_t phi_max =
      phi_0 + ((ntel - 1.0) * dphi) + 0.5 * (kvt->GetAzimuthalWidth());
   phi_max = (phi_max > 360. ? phi_max - 360. : phi_max);
   ring->SetAzimuthalMinMax(phi_min, phi_max);
   ring->SetNumber(number);
 
   Double_t phi = phi_0;
   UInt_t counter = num_first;
   for (int i = 0; i < ntel; i++) {
      kvt->SetPolarMinMax(thmin, thmax);
      kvt->SetPhi(phi);
      phi = (phi + dphi > 360. ? (phi + dphi - 360.) : (phi + dphi));
      kvt->SetNumber(counter);
      if (!absent.Contains(counter)) ring->Add(kvt);
      else delete kvt;
      counter += step;
      if (i + 1 < ntel) kvt = BuildTelescope(prefix, counter);
   }
   ring->SetName(Form("RING_%d", number));
   return ring;
}
 
 
 
 
KVINDRATelescope* KVINDRA::BuildTelescope(const Char_t* prefix, Int_t module)
{
   // Build telescope from infos in file "$KVROOT/KVFiles/data/indra-struct.[dataset].env"
 
   //Info("BuildTelescope", "Building telescope %s",name);
   KVINDRATelescope* kvt = new KVINDRATelescope;
   KVString telescopes = fStrucInfos.GetValue(Form("%s.Telescope", prefix), "");
   telescopes.Begin(" ");
   KVString name;
   while (!telescopes.End()) {
      name = telescopes.Next();
      KVNumberList mods = fStrucInfos.GetValue(Form("%s.Telescope.%s.Modules", prefix, name.Data()), "1-100");
      if (mods.Contains(module)) break;
   }
 
   KVString detectors = fStrucInfos.GetValue(Form("INDRA.Telescope.%s.Detectors", name.Data()), "");
   Double_t aziwidth = fStrucInfos.GetValue(Form("INDRA.Telescope.%s.AziWidth", name.Data()), 1.0);
   Double_t dist = fStrucInfos.GetValue(Form("%s.Dist", prefix), 0.0);
   kvt->SetDistance(dist);
   detectors.Begin(" ");
   Int_t idet = 1;
   while (!detectors.End()) {
      KVString detector = detectors.Next();
      detector.Begin("()");
      KVString dettype = detector.Next();
      Double_t detthick = detector.Next().Atof();
      KVDetector* det = KVDetector::MakeDetector(Form("%s.%s", fDataSet.Data(), dettype.Data()), detthick);
      kvt->Add(det);
      Double_t depth = fStrucInfos.GetValue(Form("INDRA.Telescope.%s.Depth.%s", name.Data(), dettype.Data()), 0.);
      kvt->SetDepth(idet, depth);
      idet++;
   }
   kvt->SetAzimuthalWidth(aziwidth);
   return kvt;
}
 
 
 
 
void KVINDRA::FillTrajectoryIDTelescopeLists()
{
   // Kludge to make INDRA ROOT geometry work like any other
   // Normally ID telescopes are deduced from successive detectors on the different trajectories
   // Each trajectory then possesses its list of ID telescopes
   // These lists are then used when reconstruction trajectories are calculated
   //
   // When INDRA ROOT geometry is used, trajectory lists need to be filled
   // by hand before reconstruction trajectories are calculated
 
   TIter it_traj(GetTrajectories());
   KVGeoDNTrajectory* tr;
   while ((tr = (KVGeoDNTrajectory*)it_traj())) {
      // get first detector on trajectory
      KVDetector* det = tr->GetNodeAt(0)->GetDetector();
      // list of 'aligned' id telescopes
      TIter it_idt(det->GetAlignedIDTelescopes());
      KVIDTelescope* idt;
      while ((idt = (KVIDTelescope*)it_idt())) {
         if (tr->ContainsAll(idt->GetDetectors())) { // all detectors on trajectory
            if (idt->GetDetectors()->GetEntries() > 1) {
               // make sure dE and E detector are consecutive on trajectory
               // i.e. dE has to be immediately in front of E on the trajectory
               if (tr->GetNodeInFront(idt->GetDetector(2)->GetNode()) == idt->GetDetector(1)->GetNode())
                  tr->AccessIDTelescopeList()->Add(idt);
            }
            else
               tr->AccessIDTelescopeList()->Add(idt);//single-detector telescope
         }
      }
   }
}
 
 
 
 
 
void KVINDRA::Build(Int_t run)
{
   // Overrides KVASMultiDetArray::Build
   // Correspondance between CsI detectors and pin lasers is set up if known.
 
   if (run != -1) {
      fCurrentRun = run;
   }
   BuildGeometry();
 
   MakeListOfDetectors();
 
   SetGroupsAndIDTelescopes();
 
   SetNamesOfIDTelescopes();
 
   SetDetectorThicknesses();
 
   SetIdentifications();
 
   //set flag to say Build() was called
   SetBit(kIsBuilt);
 
   SetPinLasersForCsI();
 
   SetExpectedDetectorSignalNames();
 
   // link EBYEDAT parameter names to detectors
   TIter next_det(GetDetectors());
   KVDetector* det;
   while ((det = (KVDetector*)next_det())) {
      TString det_name(det->GetName());
      if (det_name.BeginsWith("CI") || det_name.BeginsWith("SI")) {
         fEbyedatParamDetMap.SetValue(det_name + "_PG", det_name);
         fEbyedatParamDetMap.SetValue(det_name + "_GG", det_name);
      }
      else if (det_name.BeginsWith("CSI")) {
         fEbyedatParamDetMap.SetValue(det_name + "_R", det_name);
         fEbyedatParamDetMap.SetValue(det_name + "_L", det_name);
      }
      // handle change of "time marker" parameter type for certain experiments
      TString time_marker_type = gDataSet->GetDataSetEnv(Form("KVACQParam.%s.T", det_name.Data()), "T");
      time_marker_type.Prepend("_");
      fEbyedatParamDetMap.SetValue(det_name + time_marker_type, det_name);
   }
}
 
 
//void KVINDRA::SetArrayACQParams()
//{
//   // Overrides KVASMultiDetArray::SetArrayACQParams() in order to
//   // add the following acquisition parameters which are not associated to a detector:
//   //
//   // STAT_EVE
//   // R_DEC
//   // CONFIG
//   // PILA_01_PG
//   // PILA_01_GG
//   // PILA_02_PG
//   // PILA_02_GG
//   // PILA_03_PG
//   // PILA_03_GG
//   // PILA_04_PG
//   // PILA_04_GG
//   // PILA_05_PG
//   // PILA_05_GG
//   // PILA_06_PG
//   // PILA_06_GG
//   // PILA_07_PG
//   // PILA_07_GG
//   // PILA_08_PG
//   // PILA_08_GG
//   // SI_PIN1_PG
//   // SI_PIN1_GG
//   // SI_PIN2_PG
//   // SI_PIN2_GG
//   //
//   // We also create and initialize the KVINDRATriggerInfo object (fSelecteur) used to read
//   // the status of the DAQ trigger event by event (access through GetTriggerInfo()).
 
//   KVACQParam* ste = new KVACQParam("STAT_EVE");
//   AddArrayACQParam(ste);
//   KVACQParam* dec = new KVACQParam("R_DEC");
//   AddArrayACQParam(dec);
//   KVACQParam* conf = new KVACQParam("CONFIG");
//   AddArrayACQParam(conf);
//   fSelecteur = new KVINDRATriggerInfo;
//   fSelecteur->SetSTAT_EVE_PAR(ste);
//   fSelecteur->SetR_DEC_PAR(dec);
//   fSelecteur->SetVXCONFIG_PAR(conf);
 
//   AddArrayACQParam(new KVACQParam("PILA_01_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_01_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_02_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_02_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_03_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_03_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_04_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_04_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_05_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_05_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_06_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_06_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_07_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_07_GG"));
//   AddArrayACQParam(new KVACQParam("PILA_08_PG"));
//   AddArrayACQParam(new KVACQParam("PILA_08_GG"));
//   AddArrayACQParam(new KVACQParam("SI_PIN1_PG"));
//   AddArrayACQParam(new KVACQParam("SI_PIN1_GG"));
//   AddArrayACQParam(new KVACQParam("SI_PIN2_PG"));
//   AddArrayACQParam(new KVACQParam("SI_PIN2_GG"));
//}
 
 
 
 
void KVINDRA::MakeListOfDetectors()
{
   //Overrides KVASMultiDetArray method to add FillListsOfDetectorsByType()
 
   KVASMultiDetArray::MakeListOfDetectors();
   FillListsOfDetectorsByType();
}
 
 
 
 
void KVINDRA::FillListsOfDetectorsByType()
{
   //Fill lists of ChIo, Si, CsI and phoswich
 
   fChIo->Clear();
   fSi->Clear();
   fCsI->Clear();
   fPhoswich->Clear();
   TIter next_det(GetDetectors());
   KVDetector* kvd;
   while ((kvd = (KVDetector*) next_det())) {
      kvd->SetNameOfArray("INDRA");
      if (kvd->InheritsFrom("KVChIo")) {
         fChIo->Add(kvd);
      }
      if (kvd->InheritsFrom("KVSilicon")) {
         fSi->Add(kvd);
      }
      if (kvd->InheritsFrom("KVCsI")) {
         fCsI->Add(kvd);
      }
      if (kvd->InheritsFrom("KVPhoswich")) {
         fPhoswich->Add(kvd);
      }
   }
}
 
 
 
 
 
void KVINDRA::SetGroupsAndIDTelescopes()
{
   //Find groups of telescopes in angular alignment placed on different layers.
   //List is in fGroups.
   //Also creates all ID telescopes in array and stores them in fIDTelescopes.
   //Any previous groups/idtelescopes are deleted beforehand.
   CalculateGroupsFromGeometry();
 
   //now read list of groups and create list of ID telescopes
   CreateIDTelescopesInGroups();
}
 
 
 
 
KVChIo* KVINDRA::GetChIoOf(const Char_t* detname)
{
 
   //Returns a pointer to the Ionisation Chamber placed directly in front of the
   //detector "detname". If no ChIo is present, a null pointer is returned.
 
   KVINDRADetector* kvd;
   kvd = dynamic_cast<KVINDRADetector*>(GetDetector(detname));
   return (kvd ? (KVChIo*)kvd->GetChIo() : NULL);
}
 
 
 
 
 
KVLayer* KVINDRA::GetChIoLayer(void)
{
   //Return pointer to layer in INDRA structure corresponding to ionisation
   //chambers.
 
   return (KVLayer*)GetStructure("LAYER", "CHIO");
}
 
 
 
 
 
void KVINDRA::SetTrigger(UChar_t trig)
{
   // Define multiplicity trigger used for acquisition and filter.
   // Events with multipicity >= trig are OK.
 
   fTrigger = trig;
}
 
 
 
 
KVINDRADetector* KVINDRA::GetDetectorByType(UInt_t cou, UInt_t mod, UInt_t type) const
{
   //Find a detector based on the old BaseIndra type definitions:
   //
//enum EBaseIndra_type {
//      ChIo_GG=1,
//      ChIo_PG,//=2
//      ChIo_T,//=3
//      Si_GG,//=4
//      Si_PG,//=5
//      Si_T,//=6
//      CsI_R,//=7
//      CsI_L,//=8
//      CsI_T,//=9
//      Si75_GG,//=10
//      Si75_PG,//=11
//      Si75_T,//=12
//      SiLi_GG,//=13
//      SiLi_PG,//=14
//      SiLi_T//=15
//};
//enum EBaseIndra_typePhos {
//      Phos_R=1,
//      Phos_L,//=2
//      Phos_T,//=3
//};
 
   KVINDRADetector* det = 0;
 
   char nom_det[20];
   if (cou == 1) {
      if (type >= Phos_R && type <= Phos_T) {
         sprintf(nom_det, "PHOS_%02d", mod);
         det =
            (KVINDRADetector*) GetListOfPhoswich()->FindObjectByName(nom_det);
      }
   }
   if (det)
      return det;
   else if (type >= ChIo_GG && type <= ChIo_T) {
      sprintf(nom_det, "CI_%02d%02d", cou, mod);
      det = (KVINDRADetector*) GetListOfChIo()->FindObject(nom_det);
      return det;
   }
   else if (type >= Si_GG && type <= Si_T) {
      sprintf(nom_det, "SI_%02d%02d", cou, mod);
      det = (KVINDRADetector*) GetListOfSi()->FindObject(nom_det);
      return det;
   }
   else if (type >= SiLi_GG && type <= SiLi_T) {
      sprintf(nom_det, "SILI_%02d", cou);
      det = (KVINDRADetector*) GetListOfSi()->FindObject(nom_det);
      return det;
   }
   else if (type >= Si75_GG && type <= Si75_T) {
      sprintf(nom_det, "SI75_%02d", cou);
      det = (KVINDRADetector*) GetListOfSi()->FindObject(nom_det);
      return det;
   }
   else if (type >= CsI_R && type <= CsI_T) {
      sprintf(nom_det, "CSI_%02d%02d", cou, mod);
      det = (KVINDRADetector*) GetListOfCsI()->FindObject(nom_det);
      return det;
   }
   return 0;
}
 
 
 
 
 
Int_t KVINDRA::GetIDTelescopes(KVDetector* de, KVDetector* e, TCollection* idtels)
{
   //Override KVASMultiDetArray method for special case of "etalon" modules:
   //we need to add ChIo-CsI identification telescope by hand
 
   Int_t n = KVASMultiDetArray::GetIDTelescopes(de, e, idtels);
 
   if (de->InheritsFrom("KVSiLi") && e->InheritsFrom("KVCsI")) {
      KVChIo* chio = (KVChIo*)((KVINDRADetector*)e)->GetChIo();
      if (chio) {
         n += KVASMultiDetArray::GetIDTelescopes(chio, e, idtels);
      }
   }
   return n;
}
 
 
 
 
void KVINDRA::SetNamesOfIDTelescopes() const
{
   // Change default names of ID telescopes to INDRA standard
   //
   // This method also sets the types of the ID telescopes
 
   TIter it(GetListOfIDTelescopes());
   KVIDTelescope* idt;
   const TString etalon_numbers[] = {"1002", "1102", "1202", "1304", "1403", "1503", "1602", "1702"};
   while ((idt = (KVIDTelescope*)it())) {
      KVString N = idt->GetDetector(1)->GetName();
      N.Begin("_");
      KVString de_type = N.Next();
      KVString de_number = N.Next();
      if (idt->GetSize() == 1) {
         if (de_type == "PHOS") {
            idt->SetName(Form("%s_R_L_%s", de_type.Data(), de_number.Data()));
            idt->SetType(Form("%s_R_L", de_type.Data()));
         }
         else {
            // CsI identification telescopes are called either CSI_R_L_RRMM (before 2021)
            // or CSI_RRMM (for datasets from 2021 onwards)
            TString csi_id_name_fmt = KVBase::GetDataSetEnv(fDataSet, "INDRA.CSI.IDTelescopeNameFormat", "CSI_%s");
            idt->SetName(Form(csi_id_name_fmt, de_number.Data()));
            // similarly, the type of CSI identifications was CSI_R_L before 2021, CSI from 2021 onwards
            if (csi_id_name_fmt.Contains("R_L")) {
               idt->SetType("CSI_R_L");
               KVINDRAGroupReconstructor::CSI_ID_TYPE = "CSI_R_L";
            }
            else idt->SetType("CSI");
         }
      }
      else {
         N = idt->GetDetector(2)->GetName();
         N.Begin("_");
         KVString e_type = N.Next();
         KVString e_number = N.Next();
         if (e_type == "SILI" || e_type == "SI75") {
            e_number = etalon_numbers[dynamic_cast<KVINDRADetector*>(idt->GetDetector(2))->GetRingNumber() - 10];
         }
         idt->SetName(Form("%s_%s_%s", de_type.Data(), e_type.Data(), e_number.Data()));
         idt->SetType(Form("%s_%s", de_type.Data(), e_type.Data()));
      }
   }
   // changed names of all id telescope objects in hashlist: must rehash
   dynamic_cast<KVHashList*>(GetListOfIDTelescopes())->Rehash();
}
 
 
 
 
void KVINDRA::PerformClosedROOTGeometryOperations()
{
   // Finalise the ROOT geometry description by performing operations which can
   // only be done once the geometry is closed
   CreateROOTGeometry();
}
 
 
 
 
void KVINDRA::handle_ebyedat_raw_data_parameter(const char* param_name, uint16_t val)
{
   KVString detname;
   KVDetector* det = nullptr;
   KVString sig_type;
   // use look-up table parname => detectorname
   if (fEbyedatParamDetMap.HasParameter(param_name)) {
      det = GetDetector(fEbyedatParamDetMap.GetStringValue(param_name));
      detname = det->GetName();
      // parameter type given by whatever comes after final '_'
      std::string lab(param_name);
      sig_type = lab.substr(lab.rfind('_') + 1);
   }
   else {   // no detector associated to parameter
      sig_type = param_name;
   }
   add_and_set_detector_signal(det, detname, val, sig_type);
}
 
 
 
 
void KVINDRA::copy_fired_parameters_to_recon_param_list()
{
   // Put values of all fired detector signals into the parameter list which will be copied
   // into the reconstructed event.
   //
   // If Mesytec data is being read, we set the parameter "INDRA.MESYTEC", or if Ebyedat data, "INDRA.EBYEDAT".
   // Note that in old reconstructed data these parameters did not exist, therefore in their
   // absence EBYEDAT is assumed.
 
   if (fMesytecData) fReconParameters.SetValue("INDRA.MESYTEC", kTRUE);
   else if (fEbyedatData) fReconParameters.SetValue("INDRA.EBYEDAT", kTRUE);
   KVASMultiDetArray::copy_fired_parameters_to_recon_param_list();
}
 
 
 
 
void KVINDRA::ReadChIoPressures(KVExpDB* db)
{
   // Read ChIo pressures from a file with the following format:
   //
   //~~~
   //PressureSet: A
   //A.RunList: 1-22
   //A.CI_02: 50.0
   //A.CI_04: 30.0
   //A.CI_06: 20.0
   //
   //+PressureSet: B
   //B.Units: torr
   //B.RunList: 23-50
   //B.CI_08: 33.0
   //B.CI_13: 11.6
   //~~~
   //
   // `Units` can be either `mbar` or `torr`: `mbar` is used by default (e.g. for set `A` in the example).
   //
   // The different ionization chambers are designated by `CI_02`, `CI_04`, `CI_06`, `CI_08`, `CI_13`
   // The same pressure will be applied to all cells of each ChIo which are present in the current geometry.
   //
   // If `RunList` is not given, the pressures will be applied to all runs in the database.
 
   TString fullpath;
   if (!db->FindCalibFile("Pressures", fullpath, "INDRA")) return;
   Info("ReadChIoPressures()", "Reading INDRA ChIo pressures used during runs...");
 
   auto pressure_table = db->AddTable("INDRA.ChIo Pressures", "INDRA ChIo gas pressures");
 
   KVString the_chios[] = {"CI_02", "CI_04", "CI_06", "CI_08", "CI_13"};
 
   TEnv env;
   env.ReadFile(fullpath.Data(), kEnvAll);
   KVString pressure_sets = env.GetValue("PressureSet", "");
   pressure_sets.Begin(" ");
   while (!pressure_sets.End()) {
      auto pressure_set = pressure_sets.Next();
      // set units
      auto units = KVUnits::mbar;
      if (KVString(env.GetValue(Form("%s.Units", pressure_set.Data()), "mbar")) != "mbar")
         units = KVUnits::torr;
      // set runlist if given
      KVNumberList runlist(env.GetValue(Form("%s.RunList", pressure_set.Data()), ""));
      if (runlist.IsEmpty()) runlist = db->GetRunList();
      // get and set pressures
      auto dbp = new KVDBChIoPressures;
      int ichio = 0;
      for (auto& chio : the_chios) {
         Float_t P = 0;
         if ((P = env.GetValue(Form("%s.%s", pressure_set.Data(), chio.Data()), 0.0)) > 0.0) {
            dbp->SetPressure(ichio, P * units);
         }
         ++ichio;
      }
      pressure_table->AddRecord(dbp);
      db->LinkRecordToRunRange(dbp, runlist);
   }
}
 
 
 
 
void KVINDRA::SetChIoPressures(KVDBRun* kvrun)
{
   //Update ChIo pressures for this run with values in database (if any)
   //
   //Any ChIo with zero pressure has KVDetector::IsDetecting() set to false
   //(and hence KVDetector::IsOK() will also return false).
 
   KVRList* param_list = kvrun->GetLinks("INDRA.ChIo Pressures");
   if (!param_list) {
      return;
   }
   KVChIo* kvd;
   KVDBChIoPressures* kvps;
   TIter next_ps(param_list);
   KVSeqCollection* chios = GetListOfChIo();
   if (!chios) {
      Error("SetChIoPressures",
            "GetListOfChIo() returns null list pointer");
      return;
   }
   cout << "--> Setting ChIo pressures" << endl;
   TIter next_chio(chios);
   while ((kvps = (KVDBChIoPressures*) next_ps())) {
      if (kvps->GetPressure(CHIO_2_3))
         cout << "       Ring 2/3: " << kvps->
              GetPressure(CHIO_2_3) / KVUnits::mbar << " mbar" << endl;
      if (kvps->GetPressure(CHIO_4_5))
         cout << "       Ring 4/5: " << kvps->
              GetPressure(CHIO_4_5) / KVUnits::mbar << " mbar" << endl;
      if (kvps->GetPressure(CHIO_6_7))
         cout << "       Ring 6/7: " << kvps->
              GetPressure(CHIO_6_7) / KVUnits::mbar << " mbar" << endl;
      if (kvps->GetPressure(CHIO_8_12))
         cout << "      Ring 8/12: " << kvps->
              GetPressure(CHIO_8_12) / KVUnits::mbar << " mbar" << endl;
      if (kvps->GetPressure(CHIO_13_17))
         cout << "     Ring 13/17: " << kvps->
              GetPressure(CHIO_13_17) / KVUnits::mbar << " mbar" << endl;
      while ((kvd = (KVChIo*) next_chio())) {
         if (!strcmp(kvd->GetType(), "CI")) {
            //check detector type: ="CI" for standard INDRA chio
            if (kvd->GetRingNumber() == 2)
               kvd->SetPressure(kvps->GetPressure(CHIO_2_3) / KVUnits::mbar);
            if (kvd->GetRingNumber() == 4)
               kvd->SetPressure(kvps->GetPressure(CHIO_4_5) / KVUnits::mbar);
            if (kvd->GetRingNumber() == 6)
               kvd->SetPressure(kvps->GetPressure(CHIO_6_7) / KVUnits::mbar);
            if (kvd->GetRingNumber() >= 8 && kvd->GetRingNumber() <= 12)
               kvd->SetPressure(kvps->GetPressure(CHIO_8_12) / KVUnits::mbar);
            if (kvd->GetRingNumber() >= 13 && kvd->GetRingNumber() <= 17)
               kvd->SetPressure(kvps->GetPressure(CHIO_13_17) / KVUnits::mbar);
            if (kvd->GetPressure() == 0.0) {
               kvd->SetDetecting(kFALSE);
            }
            else {
               kvd->SetDetecting(kTRUE);
            }
         }
      }
      next_chio.Reset();
   }
}
 
 
 
 
void KVINDRA::SetDetectorParametersForRun(KVDBRun* dbr, const TString& myname)
{
   // Override base method in order to set ChIo pressures (if defined)
   KVASMultiDetArray::SetDetectorParametersForRun(dbr, myname);
   SetChIoPressures(dbr);
}
 
 
#ifdef WITH_BUILTIN_GRU
 
 
Bool_t KVINDRA::handle_raw_data_event_ebyedat(KVGANILDataReader& rdr)
{
   // Set raw data in detectors/array coming from a GANIL EBYEDAT format
   // acquisition file.
 
   fEbyedatData = kTRUE;
   // loop over fired data parameters
   TIter nxt_frd(&rdr.GetFiredDataParameters());
   KVEBYEDAT_ACQParam* eby_par;
   while ((eby_par = (KVEBYEDAT_ACQParam*)nxt_frd())) {
      handle_ebyedat_raw_data_parameter(eby_par->GetName(), eby_par->GetData());
   }
 
   return kTRUE;
}
 
#endif
 
#ifdef WITH_MFM
 
 
Bool_t KVINDRA::handle_raw_data_event_mfmframe_ebyedat(const MFMEbyedatFrame& f)
{
   // General method for reading raw data in MFM-encapsulated ebyedat format
   // Fills list of hit acquisition parameters.
   // Returns kTRUE if at least one parameter belonging to the array is present.
   //
   // Any unknown parameters in the event (i.e. ones for which no KVACQParam object
   // has been defined) are written in the fReconParameters list with names
   //    "ACQPAR.[array name].[parameter name]"
   //
   // Retrieve CENTRUM timestamp from data if present.
   // It will be added to fReconParameters as a 64-bit value "INDRA.TS" (if != 0)
   // Event number is retrieved and stored as "INDRA.EN" (if != 0)
   // Any parameter which appears as [name] and [name]_UP is an unsigned 32-bit value
   // split into two 16-bit words. We replace the two parameters with a 64-bit
   // value (to hold correctly all unsigned 32-bit values) with [name].
 
   fEbyedatData = kTRUE;
 
   uint16_t val;
   std::string lab;
 
   for (int i = 0; i < f.GetNbItems(); ++i) {
      f.GetDataItem(i, lab, val);
      handle_ebyedat_raw_data_parameter(lab.c_str(), val);
   }
 
   ULong64_t ts = f.GetCENTRUMTimestamp();
   if (ts != 0) fReconParameters.SetValue64bit("INDRA.TS", ts);
   ULong64_t en = f.GetEventNumber();
   if (en != 0) fReconParameters.SetValue64bit("INDRA.EN", en);
   int npars = fReconParameters.GetNpar();
   std::vector<TString> names;
   for (int i = 0; i < npars; ++i) {
      TString name = fReconParameters.GetParameter(i)->GetName();
      if (name.EndsWith("_UP")) names.push_back(name);
   }
   if (names.size()) {
      for (std::vector<TString>::iterator it = names.begin(); it != names.end(); ++it) {
         TString name = (*it);
         name.Remove(name.Index("_UP"), 3);
         TString name_up = (*it);
         ULong64_t par_up = fReconParameters.GetIntValue(name_up);
         ULong64_t par = fReconParameters.GetIntValue(name);
         UInt_t par32 = (par_up << 16) + par;
         fReconParameters.RemoveParameter(name_up);
         fReconParameters.RemoveParameter(name);
         fReconParameters.SetValue64bit(name, par32);
      }
   }
   return kTRUE;
}
 
#endif
 
 
 
void KVINDRA::SetIDCodeForIDTelescope(KVIDTelescope* idt) const
{
   // Set the INDRA-specific general identification code for the given telescope
 
   if (idt->InheritsFrom(KVIDPhoswich::Class())) idt->SetIDCode(IDCodes::ID_PHOSWICH);
   else if (idt->InheritsFrom(KVIDINDRACsI::Class())) idt->SetIDCode(IDCodes::ID_CSI_PSA);
   else if (idt->InheritsFrom(KVIDSiCsI::Class())) idt->SetIDCode(IDCodes::ID_SI_CSI);
   else if (idt->InheritsFrom(KVIDChIoSi::Class())) idt->SetIDCode(IDCodes::ID_CI_SI);
   else if (idt->InheritsFrom(KVIDChIoCsI::Class())) idt->SetIDCode(IDCodes::ID_CI_CSI);
   else if (idt->InheritsFrom(KVIDChIoSi75::Class())) idt->SetIDCode(IDCodes::ID_CI_SI75);
   else if (idt->InheritsFrom(KVIDSi75SiLi::Class())) idt->SetIDCode(IDCodes::ID_SI75_SILI);
   else if (idt->InheritsFrom(KVIDChIoSi75::Class())) idt->SetIDCode(IDCodes::ID_CI_SI75);
   else {
      Error("SetIDCodeForIDTelescope", "Request for telescope name=%s of unknown class=%s",
            idt->GetName(), idt->IsA()->GetName());
   }
}
 
 
 
 
void KVINDRA::SetPinLasersForCsI()
{
   // Sets the KVCsI::fPinLaser member of each CsI detector with the number of the
   // pin laser associated for the stability control of these detectors.
   //
   // We look for a file with the following format:
   //
   // CSI_0101     1
   // CSI_0102     1
   // CSI_0103     1
   // CSI_0104     1
   // etc.
   //
   // i.e. 'name of CsI detector'   'number of pin laser (1-8)'
   // Comment lines must begin with '#'
   //
   // The default name of this file is defined in .kvrootrc by
   //
   // INDRADB.CsIPinCorr:          CsI_PILA.dat
   //
   // Dataset-specific version can be specified:
   //
   // INDRA_e999.INDRADB.CsIPinCorr:    CorrCsIPin_2054.dat
   //
   // This file should be in the directory corresponding to the current dataset,
   // i.e. in $KVROOT/KVFiles/name_of_dataset
 
   ifstream pila_file;
   if (gDataSet->OpenDataSetFile(gDataSet->GetDataSetEnv("INDRADB.CsIPinCorr", ""), pila_file)) {
 
      Info("SetPinLasersForCsI", "Setting correspondance CsI-PinLaser using file %s.",
           gDataSet->GetDataSetEnv("INDRADB.CsIPinCorr", ""));
      // read file, set correspondance
      KVString line;
      line.ReadLine(pila_file);
      while (pila_file.good()) {
         if (!line.BeginsWith("#")) {
            line.Begin(" ");
            KVString detname = line.Next(kTRUE);
            KVCsI* det = (KVCsI*)GetDetector(detname.Data());
            Int_t pila = line.Next(kTRUE).Atoi();
            if (det) {
               det->SetPinLaser(pila);
            }
         }
         line.ReadLine(pila_file);
      }
      pila_file.close();
   }
   else {
      Info("SetPinLasersForCsI", "File %s not found. Correspondance Csi-PinLaser is unknown.",
           gDataSet->GetDataSetEnv("CsIPinCorr", ""));
   }
}
 
 
 
//void KVINDRA::LinkToCodeurs()
//{
 
//   // Link detectors with electronic modules
//   // for the moment only QDC for Si and ChIo are implemented
//   // This information is accessible via KVINDRADetector::GetNumeroCodeur()
//   // To be active one has to put in the dataset directory
//   // a file name Codeurs.dat containing the name of the file for the concerned type
//   // of electronic module
//   // for example see INDRA_e613 dataset
//   // [dataset name].INDRADB.Codeurs:    ...
 
 
//   KVFileReader flist;
//   TString fp;
//   if (!KVBase::SearchKVFile(gDataSet->GetDataSetEnv("INDRADB.Codeurs", ""), fp, gDataSet->GetName())) {
//      Warning("LinkToCodeurs", "Fichier %s, inconnu au bataillon", gDataSet->GetDataSetEnv("INDRADB.Codeurs", ""));
//      return;
//   }
 
//   if (!flist.OpenFileToRead(fp.Data())) {
//      //Error("ReadGainList","Error opening file named %s",fp.Data());
//      return;
//   }
//   Info("LinkToCodeurs()", "Reading correspondance Codeur-Detecteur ...");
 
//   TEnv* env = 0;
//   KVINDRADetector* idet = 0;
//   while (flist.IsOK()) {
//      flist.ReadLine(NULL);
//      KVString file = flist.GetCurrentLine();
//      if (file != "") {
//         if (KVBase::SearchKVFile(file.Data(), fp, gDataSet->GetName())) {
//            env = new TEnv();
//            env->ReadFile(fp.Data(), kEnvAll);
//            TEnvRec* rec = 0;
//            TObjArray* toks = 0;
//            TIter it(env->GetTable());
//            while ((rec = (TEnvRec*)it.Next())) {
//               if (!strcmp(rec->GetName(), "type")) {
//                  Info("LinkToCodeurs", "Module type %s", rec->GetValue());
//               }
//               else {
//                  toks = TString(rec->GetValue()).Tokenize(",");
//                  for (Int_t ii = 0; ii < toks->GetEntries(); ii += 1) {
//                     idet = (KVINDRADetector*)gIndra->GetDetector(((TObjString*)toks->At(ii))->GetString().Data());
//                     if (idet)
//                        idet->SetNumeroCodeur(TString(rec->GetName()).Atoi());
//                  }
//                  delete toks;
//               }
//            }
//            delete env;
//         }
//      }
//   }
 
//   flist.CloseFile();
 
 
//}
 
 
 
 
 
void KVINDRA::GetDetectorEvent(KVDetectorEvent* detev, const TSeqCollection* fired_dets)
{
   // Overrides KVASMultiDetArray::GetDetectorEvent.
   // If the list of fired detectors is given (meaning we are reading raw data)
   // then we check that what we have read is in fact an INDRA event
   // (see KVINDRATriggerInfo::IsINDRAEvent()) : if not, we do not try to find the hit groups.
 
   if (((fired_dets && fired_dets->GetEntries()) || fFiredDetectors.GetEntries())
         && (GetTriggerInfo() && !GetTriggerInfo()->IsINDRAEvent())) return;
   KVASMultiDetArray::GetDetectorEvent(detev, fired_dets);
}
 
 
 
//void KVINDRA::SetGGtoPGConversionFactors()
//{
//   // Sets the parameters for linear conversion of silicon & ChIo coder values
//   // between GG and PG, using the following formula:
//   //
//   //   PG  = alpha + beta*(GG - GG_0) + PG_0
//   //
//   // where GG_0 and PG_0 are respectively GG and PG pedestals
//   //
//   // We look for the file whose name is given by the .kvrootrc variable
//   //    [dataset].INDRADB.GGtoPGFactors:
//   // or by default
//   //    INDRADB.GGtoPGFactors:
//   // and expect to find in it a line for each detector of the form:
//   //    Det_Name   alpha   beta
//   // Comments in the file can be written on lines beginning with the character '#'
 
//   ifstream datfile;
//   if (!gDataSet->OpenDataSetFile(gDataSet->GetDataSetEnv("INDRADB.GGtoPGFactors", ""), datfile)) {
 
//      Info("SetGGtoPGConversionFactors", "Cannot open file with parameters for conversion (%s).",
//           gDataSet->GetDataSetEnv("INDRADB.GGtoPGFactors", ""));
//      return;
//   }
//   else {
//      Info("SetGGtoPGConversionFactors", "Reading parameters from file %s",
//           gDataSet->GetDataSetEnv("INDRADB.GGtoPGFactors", ""));
 
//      Char_t detname[30];
//      Double_t a, b;
//      TString aline;
//      aline.ReadLine(datfile);
//      while (datfile.good()) {
 
//         if (aline[0] != '#') {   //skip comments
 
//            sscanf(aline.Data(), "%s %lf %lf", detname, &a, &b);
//            KVINDRADetector* det = (KVINDRADetector*)GetDetector(detname);
//            if (!det) {
//               //no detector found with cou, mod and type
//               Error("SetGGtoPGConversionFactors", "Unknown detector : %s", detname);
//            }
//            else {
//               det->SetGGtoPGConversionFactors(a, b);
//               //Info("SetGGtoPGConversionFactors", "%s : PG = %f + %f * GG", detname, a, b);
//            }
//         }
//         aline.ReadLine(datfile);
//      }                            //while( datfile.good()
//      datfile.close();
//   }
//}
 
 
 
 
void KVINDRA::CreateROOTGeometry()
{
   // Overrides KVASMultiDetArray::CreateGeoManager in order to use INDRAGeometryBuilder
   // which builds the TGeo representation of INDRA using the Y. Huguet CAO data.
   //
   // The optional arguments (dx,dy,dz) are the half-lengths in centimetres of the "world"/"top" volume
   // into which all the detectors of the array are placed. This should be big enough so that all detectors
   // fit in. The default values of 500 give a "world" which is a cube 1000cmx1000cmx1000cm (with sides
   // going from -500cm to +500cm on each axis).
   //
   // If closegeo=kFALSE we leave the geometry open for other structures to be added.
 
   if (!IsBuilt()) {
      Error("CreateROOTGeometry", "gIndra has to be build first");
      return;
   }
   if (!GetNavigator()) {
      //Error("CreateROOTGeometry","No existing navigator"); return;
      SetNavigator(new KVRangeTableGeoNavigator(gGeoManager, KVMaterial::GetRangeTable()));
      GetNavigator()->SetNameCorrespondanceList("INDRA.names");
   }
 
   // set up shape & matrix pointers in detectors
   Info("CreateROOTGeometry", "Scanning geometry shapes and matrices...");
   KVGeoImport gimp(gGeoManager, KVMaterial::GetRangeTable(), this, kFALSE);
   gimp.SetNameCorrespondanceList("INDRA.names");
   KVNucleusEvent evt;
   KVNucleus* nuc = evt.AddNucleus();
   nuc->SetZAandE(1, 1, 1);
   KVINDRADetector* det;
   TIter next(GetDetectors());
   Int_t nrootgeo = 0;
   while ((det = (KVINDRADetector*)next())) {
      nuc->SetTheta(det->GetTheta());
      nuc->SetPhi(det->GetPhi());
      gimp.SetLastDetector(0);
      gimp.PropagateEvent(&evt);
      if (!(det->ROOTGeo())) {
         Info("CreateROOTGeometry", "Volume checking for %s", det->GetName());
         Double_t theta0 = det->GetTheta();
         Double_t phi0 = det->GetPhi();
         for (Double_t TH = theta0 - 0.5; TH <= theta0 + 0.5; TH += 0.1) {
            for (Double_t PH = phi0 - 10; PH <= phi0 + 10; PH += 1) {
               nuc->SetTheta(TH);
               nuc->SetPhi(PH);
               gimp.SetLastDetector(0);
               gimp.PropagateEvent(&evt);
               if (det->ROOTGeo()) break;
            }
            if (det->ROOTGeo()) break;
         }
      }
      if (!(det->ROOTGeo())) {
         Info("CreateROOTGeometry", "Volume checking failed for : %s", det->GetName());
      }
      // check etalon trajectories (if etalons are present)
      if (det->ROOTGeo() && det->GetRingNumber() > 9) {
         if (GetDetector(Form("SI75_%d", det->GetRingNumber())) || GetDetector(Form("SILI_%d", det->GetRingNumber()))) {
            if ((det->IsCalled("CSI_1002") || det->IsCalled("CSI_1102")
                  || det->IsCalled("CSI_1202") || det->IsCalled("CSI_1304")
                  || det->IsCalled("CSI_1403") || det->IsCalled("CSI_1503")
                  || det->IsCalled("CSI_1602") || det->IsCalled("CSI_1702"))
                  && det->GetNode()->GetNDetsInFront() < 2) {
               Info("CreateROOTGeometry", "Trajectory checking for %s", det->GetName());
               Double_t theta0 = det->GetTheta();
               Double_t phi0 = det->GetPhi();
               for (Double_t TH = theta0 - 0.5; TH <= theta0 + 0.5; TH += 0.1) {
                  for (Double_t PH = phi0 - 10; PH <= phi0 + 10; PH += 1) {
                     nuc->SetTheta(TH);
                     nuc->SetPhi(PH);
                     gimp.SetLastDetector(0);
                     gimp.PropagateEvent(&evt);
                     if (det->GetNode()->GetNDetsInFront() == 2) break;
                  }
                  if (det->GetNode()->GetNDetsInFront() == 2) break;
               }
            }
         }
      }
      nrootgeo += (det->ROOTGeo());
   }
 
   Info("CreateROOTGeometry", "ROOT geometry initialised for %d/%d detectors", nrootgeo, GetDetectors()->GetEntries());
 
   // Set up trajectories
   TIter it(GetDetectors());
   KVDetector* d;
   while ((d = (KVDetector*)it())) d->GetNode()->RehashLists();// make sure detector nodes are correct
   AssociateTrajectoriesAndNodes();
   DeduceGroupsFromTrajectories();
   FillTrajectoryIDTelescopeLists();
   CalculateReconstructionTrajectories();
   GetNavigator()->AbsorbDetectorPaths(&gimp);
}
 
 
 
 
void KVINDRA::SetROOTGeometry(Bool_t on)
{
   // Override base class method
   // If ROOT geometry is requested but has not been built, we create it
 
   if (on) {
      CreateGeoManager();
      INDRAGeometryBuilder igb;
      igb.Build(kFALSE, fCloseGeometryNow);
      if (fCloseGeometryNow) PerformClosedROOTGeometryOperations();
   }
   else {
      KVMultiDetArray::SetROOTGeometry(on);
   }
}
 
 
 
 
KVGroupReconstructor* KVINDRA::GetReconstructorForGroup(const KVGroup* g) const
{
   // Special INDRA group reconstructors:
   //   KVINDRAForwardGroupReconstructor       rings 1-9
   //   KVINDRABackwardGroupReconstructor      rings 10-17
   //   KVINDRAEtalonGroupReconstructor for groups with etalon telescopes
 
   KVGroupReconstructor* gr(nullptr);
   if (GetGroup(g->GetName())) { // make sure group belongs to us
      // etalons ?
      if (g->GetDetectorByType("SILI") || g->GetDetectorByType("SI75"))
         gr = KVGroupReconstructor::Factory("INDRA.etalon");
      else {
         KVINDRADetector* id = (KVINDRADetector*)g->GetDetectors()->First();
         if (id->GetRingNumber() < 10) {
            gr = KVGroupReconstructor::Factory("INDRA.forward");
         }
         else {
            gr = KVGroupReconstructor::Factory("INDRA.backward");
         }
      }
   }
   return gr;
}
 
 
 
 
void KVINDRA::SetReconParametersInEvent(KVReconstructedEvent* e) const
{
   // If "INDRA.EN" parameter has been set, we use it to set the event number
 
   KVASMultiDetArray::SetReconParametersInEvent(e);
   if (GetReconParameters().HasValue64bit("INDRA.EN")) e->SetNumber(GetReconParameters().GetValue64bit("INDRA.EN"));
}
 
 
 
 
void KVINDRA::SetRawDataFromReconEvent(KVNameValueList& l)
{
   // Overrides base method in KVMultiDetArray.
   //
   // If we are reading old reconstructed data with EBYEDAT parameters, we need special
   // treatment to decode the detector name & signal type.
 
   if (!l.GetBoolValue("INDRA.MESYTEC") && !l.GetBoolValue("INDRA.EBYEDAT")) {
      // assume we are reading old recon data with parameters in list such as "ACQPAR.INDRA.SI_0201_PG"
      fMesytecData = kFALSE;
      fEbyedatData = kTRUE;
      prepare_to_handle_new_raw_data();
 
      int N = l.GetNpar();
      for (int i = 0; i < N; ++i) {
         KVNamedParameter* np = l.GetParameter(i);
 
         KVString name(np->GetName());
         if (name.BeginsWith("ACQPAR")) {
            // 2 '.' => 3 values
            int dots = name.GetNValues(".");
            assert(dots == 3); // sanity check
            name.Begin(".");
            name.Next(); // "ACQPAR"
            if (name.Next() != "INDRA") continue; // check name of array - somebody else's parameter ?
            std::string parname(name.Next());
            KVString detname;
            KVString sig_type;
            KVDetector* det = nullptr;
            if (parname.find('_') != kNPOS) {
               // try to split into '[detname]_[partype]'
               detname = parname.substr(0, parname.rfind('_'));
               sig_type = parname.substr(parname.rfind('_') + 1);
               det = GetDetector(detname);
            }
            else {
               sig_type = parname;
            }
            add_and_set_detector_signal(det, detname, np->GetDouble(), sig_type);
         }
      }
   }
   else {
      fMesytecData = l.GetBoolValue("INDRA.MESYTEC");
      fEbyedatData = l.GetBoolValue("INDRA.EBYEDAT");
      KVASMultiDetArray::SetRawDataFromReconEvent(l);
   }
}
 
 
 
 
void KVINDRA::MakeCalibrationTables(KVExpDB* db)
{
   // Override base method in order to read ChIo pressures for each run
   KVMultiDetArray::MakeCalibrationTables(db);
   // if the actual database class inherits from KVINDRADB (historic INDRA data),
   // the ChIo pressures have already been read (& from a file with a different format)
   if (!db->InheritsFrom("KVINDRADB")) ReadChIoPressures(db);
}
 
 
 
 
void KVINDRA::InitialiseRawDataReading(KVRawDataReader* r)
{
   // Call this method just after opening a raw data file in order to perform any
   // necessary initialisations, depending on the type of data
 
   KVASMultiDetArray::InitialiseRawDataReading(r);
 
   fMesytecData = kFALSE;
   fEbyedatData = kFALSE;
#ifdef WITH_MESYTEC
#ifdef WITH_MFM
   if (r->GetDataFormat() == "MFM") {
      TString crate = gDataSet->GetDataSetEnv("MesytecCrateConf", "");
      if (crate != "") {
         Info("InitialiseRawDataReading", "Setting Mesytec crate config for run %d, raw run number=%d",
              GetCurrentRunNumber(), r->GetRunNumberReadFromFile());
 
         TString channel_conf_list = gDataSet->GetDataSetEnv("MesytecChannelsConf", "");
         if (channel_conf_list != "") {
            TEnv conf_list;
            if (conf_list.ReadFile(gDataSet->GetFullPathToDataSetFile(channel_conf_list), kEnvUser) == -1) {
               Fatal("InitialiseRawDataReading",
                     "Could not open file containing names of Mesytec channel config for each run: %s",
                     gDataSet->GetFullPathToDataSetFile(channel_conf_list).Data());
            }
            KVString files = conf_list.GetValue("File", "");
            files.Begin(" ");
            while (!files.End()) {
               KVString next_file = files.Next(kTRUE);
               KVNumberList runlist = conf_list.GetValue(Form("%s.RunList", next_file.Data()), "");
               if (runlist.Contains(r->GetRunNumberReadFromFile())) {
                  // found file for this run
                  Info("InitialiseRawDataReading", "Using file %s for raw run number %d",
                       next_file.Data(), r->GetRunNumberReadFromFile());
                  dynamic_cast<KVMFMDataFileReader*>(r)->InitialiseMesytecConfig(
                     gDataSet->GetFullPathToDataSetFile(crate).Data(),
                     gDataSet->GetFullPathToDataSetFile(next_file).Data()
                  );
               }
            }
         }
         else {
            Fatal("InitialiseRawDataReading",
                  "Name of file containing names of Mesytec channel config for each run should be defined in dataset variable %s.MesytecChannelsConf",
                  GetDataSet().Data());
         }
      }
   }
#endif
#endif
}
 
 
#ifdef WITH_MFM
#ifdef WITH_MESYTEC
Bool_t KVINDRA::handle_raw_data_event_mfmframe_mesytec_mdpp(const MFMMesytecMDPPFrame& f)
{
   // Read a raw data event from a Mesytec MFM Frame.
   //
   // All data is transferred to KVDetectorSignal objects, either associated to detectors of the
   // array (if they can be identified), either associated more globally with the array/event
   // itself. The latter are created as needed and go into the fExtraRawDataSignals list.
 
   if (!fMesytecData) {
      // first time we read data, we tweak the parameter names
      mesytec::module::set_data_type_alias("qdc_long", "TotLight"); // => CsI
      mesytec::module::set_data_type_alias("qdc_short", "R"); // => like old CsI 'R' (rapide) parameter
      mesytec::module::set_data_type_alias("tdc", "T"); // => like old 'marqueur de temps' parameters
      mesytec::module::set_data_type_alias("adc", "ADC"); // => Si or ChIo signals
 
      fMesytecData = kTRUE;
   }
   fReconParameters.SetValue64bit("INDRA.TS", f.GetTimeStamp());
   auto mfmfilereader = dynamic_cast<KVMFMDataFileReader*>(fRawDataReader);
   mfmfilereader->GetMesytecBufferReader().read_event_in_buffer(
      (const uint8_t*)f.GetPointUserData(), f.GetBlobSize(),
#ifdef MESYTEC_DATA_NO_MDPP_NAMESPACE
#ifdef MESYTEC_DATA_BUFFER_READER_CALLBACK_WITH_EVENT_AND_SETUP
   [ = ](const mesytec::event & evt, const mesytec::experimental_setup & setup) {
#else
   [ = ](const mesytec::event & evt) {
#endif
#else
   [ = ](const mesytec::mdpp::event & evt) {
#endif
#ifndef MESYTEC_DATA_BUFFER_READER_CALLBACK_WITH_EVENT_AND_SETUP
      auto& setup = mfmfilereader->GetMesytecBufferReader().get_setup();
#endif
      // loop over module data in event, set data in detectors when possible
#ifdef MESYTEC_DATA_NO_PUBLIC_MEMBERS
      for (auto& mdat : evt.get_module_data()) {
         auto mod_id = mdat.get_module_id();
#else
      for (auto& mdat : evt.modules) {
         auto mod_id = mdat.module_id;
#endif
         auto& current_module = setup.get_module(mod_id);
         if (current_module.is_mdpp_module()) {
            // data for detectors from MDPP module
#ifdef MESYTEC_DATA_NO_PUBLIC_MEMBERS
            for (auto& voie : mdat.get_channel_data()) {
               KVString detname(setup.get_detector(mod_id, voie.get_channel_number()));
               KVString sig_type(current_module.get_data_type_name(voie.get_data_type()));
               Double_t sig_data = voie.get_data();
#else
            for (auto& voie : mdat.data) {
               KVString detname(setup.get_detector(mod_id, voie.channel));
               KVString sig_type(voie.data_type);
               Double_t sig_data = voie.data;
 
#endif
               add_and_set_detector_signal(GetDetector(detname), detname, sig_data, sig_type);
            }
         }
         else if (current_module.is_mvlc_scaler()) {
            // data from 64 bit scalers in MVLC
            //
            // 4 data words of 16 bits (least significant is first word) => 64 bit scaler
#ifdef MESYTEC_DATA_NO_PUBLIC_MEMBERS
            if (mdat.get_channel_data().size() == 4) {
#else
            if (mdat.data.size() == 4) {
#endif
               ULong64_t x = 0;
               int i = 0;
#ifdef MESYTEC_DATA_NO_PUBLIC_MEMBERS
               for (auto& d : mdat.get_channel_data()) x += ((uint64_t)d.get_data_word()) << (16 * (i++));
#else
               for (auto& d : mdat.data) x += ((uint64_t)d.data_word) << (16 * (i++));
#endif
               fReconParameters.SetValue64bit(current_module.name.c_str(), x);
            }
//            else {
//               Warning("handle_raw_data_event_mfmframe_mesytec_mdpp",
//                       "Scaler data %s is corrupt : %d words instead of 4",
//                       current_module.name.c_str(), mdat.data.size());
//            }
         }
      }
   },
   f.GetRevision()
   );
   return kTRUE;
}
#endif
#endif