Adding a project for GRIT AGATA at GANIL

Projects/GRIT_AGATA_GANIL/Analysis.cxx

+/*****************************************************************************
+ * Copyright (C) 2009-2014    this file is part of the NPTool Project        *
+ *                                                                           *
+ * For the licensing terms see $NPTOOL/Licence/NPTool_Licence                *
+ * For the list of contributors see $NPTOOL/Licence/Contributors             *
+ *****************************************************************************/
+
+/*****************************************************************************
+ * Original Author: Adrien MATTA  contact address: a.matta@surrey.ac.uk      *
+ *                                                                           *
+ * Creation Date  : march 2025                                               *
+ * Last update    :                                                          *
+ *---------------------------------------------------------------------------*
+ * Decription:                                                               *
+ * Class describing the property of an Analysis object                       *
+ *                                                                           *
+ *---------------------------------------------------------------------------*
+ * Comment:                                                                  *
+ *                                                                           *
+ *                                                                           *
+ *****************************************************************************/
+#include <iostream>
+using namespace std;
+
+#include "Analysis.h"
+#include "NPAnalysisFactory.h"
+#include "NPDetectorManager.h"
+#include "NPFunction.h"
+#include "NPOptionManager.h"
+////////////////////////////////////////////////////////////////////////////////
+Analysis::Analysis() {}
+////////////////////////////////////////////////////////////////////////////////
+Analysis::~Analysis() {}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::Init() {
+  if (NPOptionManager::getInstance()->HasDefinition("simulation")) {
+    cout << "Considering input data as simulation" << endl;
+    simulation = true;
+  }
+  else {
+    cout << "Considering input data as real" << endl;
+    simulation = false;
+  }
+  agata_zShift = 51 * mm;
+
+  // initialize input and output branches
+  if (simulation) {
+    Initial = new TInitialConditions();
+    ReactionConditions = new TReactionConditions();
+  }
+
+  InitOutputBranch();
+  InitInputBranch();
+  // get MUST2 and Gaspard objects
+  GPDTrack = (GaspardTracker*)m_DetectorManager->GetDetector("GaspardTracker");
+
+  // get reaction information
+  reaction.ReadConfigurationFile(NPOptionManager::getInstance()->GetReactionFile());
+  OriginalBeamEnergy = reaction.GetBeamEnergy();
+  // target thickness
+  TargetThickness = m_DetectorManager->GetTargetThickness();
+  string TargetMaterial = m_DetectorManager->GetTargetMaterial();
+  // Cryo target case
+  WindowsThickness = 0;        // m_DetectorManager->GetWindowsThickness();
+  string WindowsMaterial = ""; // m_DetectorManager->GetWindowsMaterial();
+
+  // energy losses
+  string light = NPL::ChangeNameToG4Standard(reaction.GetNucleus3()->GetName());
+  string beam = NPL::ChangeNameToG4Standard(reaction.GetNucleus1()->GetName());
+  LightTarget = NPL::EnergyLoss(light + "_" + TargetMaterial + ".G4table", "G4Table", 100);
+  //  LightAl = NPL::EnergyLoss(light + "_Al.G4table", "G4Table", 100);
+  LightSi = NPL::EnergyLoss(light + "_Si.G4table", "G4Table", 100);
+  BeamCD2 = NPL::EnergyLoss(beam + "_" + TargetMaterial + ".G4table", "G4Table", 100);
+
+  FinalBeamEnergy = BeamCD2.Slow(OriginalBeamEnergy, TargetThickness * 0.5, 0);
+  // FinalBeamEnergy = OriginalBeamEnergy;
+  cout << "Original beam energy: " << OriginalBeamEnergy << " MeV      Mid-target beam energy: " << FinalBeamEnergy
+       << "MeV " << endl;
+  reaction.SetBeamEnergy(FinalBeamEnergy);
+
+  if (WindowsThickness) {
+    cout << "Cryogenic target with windows" << endl;
+    BeamWindow = new NPL::EnergyLoss(beam + "_" + WindowsMaterial + ".G4table", "G4Table", 100);
+    LightWindow = new NPL::EnergyLoss(light + "_" + WindowsMaterial + ".G4table", "G4Table", 100);
+  }
+
+  else {
+    BeamWindow = NULL;
+    LightWindow = NULL;
+  }
+
+  // initialize various parameters
+  Rand = TRandom3();
+  DetectorNumber = 0;
+  ThetaNormalTarget = 0;
+  ThetaM2Surface = 0;
+  ThetaMGSurface = 0;
+  Si_E_M2 = 0;
+  CsI_E_M2 = 0;
+  Energy = 0;
+  ThetaGDSurface = 0;
+  X = 0;
+  Y = 0;
+  Z = 0;
+  dE = 0;
+  BeamDirection = TVector3(0, 0, 1);
+  nbTrack = 0;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::TreatEvent() {
+  // Reinitiate calculated variable
+  ReInitValue();
+  double XTarget, YTarget;
+  TVector3 BeamDirection;
+  XTarget = 0;
+  YTarget = 0;
+  BeamDirection = TVector3(0, 0, 1);
+  OriginalELab = ReactionConditions->GetKineticEnergy(0);
+  OriginalThetaLab = ReactionConditions->GetTheta(0);
+  BeamEnergy = ReactionConditions->GetBeamEnergy();
+  BeamImpact = TVector3(XTarget, YTarget, 0);
+  // determine beam energy for a randomized interaction point in target
+  // 1% FWHM randominastion (E/100)/2.35
+  // reaction.SetBeamEnergy(Rand.Gaus(ReactionConditions->GetIncidentFinalKineticEnergy(),ReactionConditions->GetIncidentFinalKineticEnergy()/235));
+
+  ////////////////////////////////////////////////////////////////////////////
+  //////////////////////////////// GASPARD  //////////////////////////////////
+  ////////////////////////////////////////////////////////////////////////////
+  if ((ELab = GPDTrack->GetEnergyDeposit()) > 0) {
+    TVector3 HitDirection = GPDTrack->GetPositionOfInteraction() - BeamImpact;
+    ELab = GPDTrack->GetEnergyDeposit();
+    ThetaLab = HitDirection.Angle(BeamDirection);
+    ELab = LightTarget.EvaluateInitialEnergy(ELab, TargetThickness * 0.5, 0);
+
+    // Part 3 : Excitation Energy Calculation
+    Ex = reaction.ReconstructRelativistic(ELab, ThetaLab);
+
+    // Part 4 : Theta CM Calculation
+    ThetaCM = reaction.EnergyLabToThetaCM(ELab, ThetaLab) / deg;
+    ThetaLab = ThetaLab / deg;
+  }
+}
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::End() {}
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::InitOutputBranch() {
+  RootOutput::getInstance()->GetTree()->Branch("Ex", &Ex, "Ex/D");
+  RootOutput::getInstance()->GetTree()->Branch("ELab", &ELab, "ELab/D");
+  RootOutput::getInstance()->GetTree()->Branch("ThetaLab", &ThetaLab, "ThetaLab/D");
+  RootOutput::getInstance()->GetTree()->Branch("ThetaCM", &ThetaCM, "ThetaCM/D");
+  RootOutput::getInstance()->GetTree()->Branch("Run", &Run, "Run/I");
+  RootOutput::getInstance()->GetTree()->Branch("X", &X, "X/D");
+  RootOutput::getInstance()->GetTree()->Branch("Y", &Y, "Y/D");
+  RootOutput::getInstance()->GetTree()->Branch("Z", &Z, "Z/D");
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::InitInputBranch() {
+  // RootInput:: getInstance()->GetChain()->SetBranchAddress("GATCONF",&vGATCONF);
+  RootInput::getInstance()->GetChain()->SetBranchStatus("InitialConditions", true);
+  RootInput::getInstance()->GetChain()->SetBranchStatus("fIC_*", true);
+  RootInput::getInstance()->GetChain()->SetBranchAddress("InitialConditions", &Initial);
+  RootInput::getInstance()->GetChain()->SetBranchStatus("ReactionConditions", true);
+  RootInput::getInstance()->GetChain()->SetBranchStatus("fRC_*", true);
+  RootInput::getInstance()->GetChain()->SetBranchAddress("ReactionConditions", &ReactionConditions);
+}
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::ReInitValue() {
+  Ex = -1000;
+  ELab = -1000;
+  BeamEnergy = -1000;
+  ThetaLab = -1000;
+  ThetaCM = -1000;
+  X = -1000;
+  Y = -1000;
+  Z = -1000;
+  dE = -1000;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+//            Construct Method to be pass to the AnalysisFactory              //
+////////////////////////////////////////////////////////////////////////////////
+NPL::VAnalysis* Analysis::Construct() { return (NPL::VAnalysis*)new Analysis(); }
+
+////////////////////////////////////////////////////////////////////////////////
+//            Registering the construct method to the factory                 //
+////////////////////////////////////////////////////////////////////////////////
+extern "C" {
+class proxy_analysis {
+ public:
+  proxy_analysis() { NPL::AnalysisFactory::getInstance()->SetConstructor(Analysis::Construct); }
+};
+
+proxy_analysis p_analysis;
+}
+

Projects/GRIT_AGATA_GANIL/Analysis.h

+#ifndef Analysis_h
+#define Analysis_h
+/*****************************************************************************
+ * Copyright (C) 2009-2014    this file is part of the NPTool Project        *
+ *                                                                           *
+ * For the licensing terms see $NPTOOL/Licence/NPTool_Licence                *
+ * For the list of contributors see $NPTOOL/Licence/Contributors             *
+ *****************************************************************************/
+
+/*****************************************************************************
+ * Original Author: Adrien MATTA  contact address: a.matta@surrey.ac.uk      *
+ *                                                                           *
+ * Creation Date  : march 2025                                               *
+ * Last update    :                                                          *
+ *---------------------------------------------------------------------------*
+ * Decription:                                                               *
+ * Class describing the property of an Analysis object                       *
+ *                                                                           *
+ *---------------------------------------------------------------------------*
+ * Comment:                                                                  *
+ *                                                                           *
+ *                                                                           *
+ *****************************************************************************/
+#include "GaspardTracker.h"
+#include "NPEnergyLoss.h"
+#include "NPReaction.h"
+#include "NPVAnalysis.h"
+#include "RootInput.h"
+#include "RootOutput.h"
+#include "TInitialConditions.h"
+#include "TReactionConditions.h"
+#include <TMath.h>
+#include <TRandom3.h>
+#include <TVector3.h>
+
+class Analysis : public NPL::VAnalysis {
+ public:
+  Analysis();
+  ~Analysis();
+
+ public:
+  void Init();
+  void TreatEvent();
+  void End();
+
+  void InitOutputBranch();
+  void InitInputBranch();
+  void ReInitValue();
+  static NPL::VAnalysis* Construct();
+
+ private:
+  double Ex;
+  double ExNoBeam;
+  double ExNoProton;
+  double EDC;
+  double ELab;
+  double ThetaLab;
+  double ThetaCM;
+  double OriginalELab;
+  double OriginalThetaLab;
+
+  NPL::Reaction reaction;
+  //	Energy loss table: the G4Table are generated by the simulation
+  NPL::EnergyLoss LightTarget;
+  NPL::EnergyLoss LightAl;
+  NPL::EnergyLoss LightSi;
+  NPL::EnergyLoss BeamCD2;
+  NPL::EnergyLoss* BeamWindow;
+  NPL::EnergyLoss* LightWindow;
+
+  double TargetThickness;
+  double WindowsThickness;
+  // Beam Energy
+  double OriginalBeamEnergy; // AMEV
+  double FinalBeamEnergy;
+
+  // intermediate variable
+  TVector3 BeamDirection;
+  TVector3 BeamImpact;
+  TRandom3 Rand;
+  int Run;
+  int DetectorNumber;
+  double ThetaNormalTarget;
+  double ThetaM2Surface;
+  double ThetaMGSurface;
+  double Si_E_M2;
+  double CsI_E_M2;
+  double Energy;
+  double BeamEnergy;
+
+  double ThetaGDSurface;
+  double X;
+  double Y;
+  double Z;
+  // Vamos Branches
+  unsigned long long int LTS;
+
+  // Agata branches
+  double agata_zShift;
+  unsigned long long int TStrack;
+  int nbHits;
+  int nbTrack;
+  float* trackE = new float(100);
+  float* trackX1 = new float(100);
+  float* trackY1 = new float(100);
+  float* trackZ1 = new float(100);
+  float* trackT = new float(100);
+  int* trackCrystalID = new int(100);
+  int nbCores;
+  int* coreId = new int(100);
+  ULong64_t* coreTS = new ULong64_t(100);
+  float* coreE0 = new float(100);
+  //
+  double dE;
+  double dTheta;
+  // Branches and detectors
+  GaspardTracker* GPDTrack;
+  bool simulation;
+  TInitialConditions* Initial;
+  TReactionConditions* ReactionConditions;
+};
+#endif

Projects/GRIT_AGATA_GANIL/detector/grit_agata.detector

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1
+Target
+ THICKNESS= 10 micrometer 
+ RADIUS= 7.5 mm 
+ MATERIAL= CD2
+ ANGLE= 0 deg
+ X= 0 mm 
+ Y= 0 mm 
+ Z= 0 mm
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias THIRDSTAGEFRONT
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias THIRDSTAGEBACK
+  Action= Replace
+  Value= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEFRONT
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEBACK
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEBARREL
+  Action= Replace
+  Value= 1 
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 Annular Back
+GaspardTracker Annular
+ Z= -156.5 mm 
+ RMIN= 16 mm 
+ RMAX= 52 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 Annular Front
+GaspardTracker Annular
+ Z= 156.5 mm 
+ RMIN= 16 mm 
+ RMAX= 52 mm 
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 End-Cap Back
+GaspardTracker Trapezoid
+ X128_Y128= 55.338 -14.346 -146.501 mm
+ X1_Y128= 55.338 14.346 -146.501 mm
+ X128_Y1= 138.519 -48.717 -69.236 mm
+ X1_Y1=  138.518 48.726 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2
+GaspardTracker Trapezoid
+ X128_Y128= 49.215 29.045 -146.501 mm
+ X1_Y128=  28.986 49.274 -146.501 mm
+ X128_Y1= 132.395 63.500 -69.236 mm
+ X1_Y1=  63.492 132.401 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3
+GaspardTracker Trapezoid
+ X128_Y128= 14.263 55.338 -146.501 mm
+ X1_Y128= -14.346 55.338 -146.501 mm
+ X128_Y1=  48.717 138.519 -69.236 mm
+ X1_Y1=  -48.726 138.518 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%4
+GaspardTracker Trapezoid
+ X128_Y128= -29.045 49.215 -146.501 mm
+ X1_Y128= -49.274 28.986 -146.501 mm
+ X128_Y1= -63.500 132.395 -69.236 mm
+ X1_Y1=  -132.401 63.492 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
+GaspardTracker Trapezoid
+ X128_Y128= -55.338 14.346 -146.501 mm
+ X1_Y128= -55.338 -14.346 -146.501 mm
+ X128_Y1= -138.519 48.717 -69.236 mm
+ X1_Y1=  -138.518 -48.726 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6
+GaspardTracker Trapezoid
+ X128_Y128= -49.215 -29.045 -146.501 mm
+ X1_Y128= -28.986 -49.274 -146.501 mm
+ X128_Y1= -132.395 -63.500 -69.236 mm
+ X1_Y1=  -63.492 -132.401 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7
+GaspardTracker Trapezoid
+ X128_Y128= -14.263 -55.338 -146.501 mm
+ X1_Y128= 14.346 -55.338 -146.501 mm
+ X128_Y1= -48.717 -138.519 -69.236 mm
+ X1_Y1=  48.726 -138.518 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8
+GaspardTracker Trapezoid
+ X128_Y128= 29.045 -49.215 -146.501 mm
+ X1_Y128= 49.274 -28.986 -146.501 mm
+ X128_Y1= 63.500 -132.395 -69.236 mm
+ X1_Y1=  132.401 -63.492 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 Barrel
+GaspardTracker Square
+ THETA= 90 deg 
+ PHI= 0 deg
+ R= 143 mm 
+ BETA= 0 0 0 deg 
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2
+GaspardTracker Square
+ THETA= 90 deg 
+ PHI= 45 deg
+ R= 143 mm
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3
+GaspardTracker Square
+ THETA= 90 deg
+ PHI= 90 deg
+ R= 143 mm
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%4
+GaspardTracker Square
+ THETA= 90 deg 
+ PHI= 135 deg 
+ R= 143 mm
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
+GaspardTracker Square
+ THETA= 90 deg
+ PHI= 180 deg 
+ R= 143 mm
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6
+GaspardTracker Square
+ THETA= 90 deg
+ PHI= 225 deg
+ R= 143 mm
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7
+GaspardTracker Square
+ THETA= 90 deg
+ PHI= 270 deg
+ R= 143 mm
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8
+GaspardTracker Square
+ THETA= 90 deg
+ PHI= 315 deg
+ R= 143 mm 
+ BETA= 0 0 0 deg
+ FIRSTSTAGE= 1
+ SECONDSTAGE= 1
+ THIRDSTAGE= 0
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 Front-Cap Back
+GaspardTracker Trapezoid
+ X128_Y128= 55.338 14.346 146.501 mm
+ X1_Y128= 55.338 -14.346 146.501 mm
+ X128_Y1= 138.518 48.726 69.237 mm
+ X1_Y1=  138.519 -48.717 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2
+GaspardTracker Trapezoid
+ X128_Y128= 28.986 49.274 146.501 mm
+ X1_Y128= 49.215 29.045 146.501 mm
+ X128_Y1= 63.492 132.401 69.237 mm
+ X1_Y1=  132.395 63.500 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3
+GaspardTracker Trapezoid
+ X128_Y128= -14.346 55.338 146.501 mm
+ X1_Y128= 14.263 55.338 146.501 mm
+ X128_Y1= -48.726 138.518 69.237 mm
+ X1_Y1=  48.717 138.519 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%4
+GaspardTracker Trapezoid
+ X128_Y128= -49.274 28.986 146.501 mm
+ X1_Y128= -29.045 49.215 146.501 mm
+ X128_Y1= -132.401 63.492 69.237 mm
+ X1_Y1=  -63.500 132.395 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
+GaspardTracker Trapezoid
+ X128_Y128= -55.338 -14.346 146.501 mm
+ X1_Y128= -55.338 14.346 146.501 mm
+ X128_Y1= -138.518 -48.726 69.237 mm
+ X1_Y1=  -138.519 48.717 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6
+GaspardTracker Trapezoid
+ X128_Y128= -28.986 -49.274 146.501 mm
+ X1_Y128= -49.215 -29.045 146.501 mm
+ X128_Y1= -63.492 -132.401 69.237 mm
+ X1_Y1=  -132.395 -63.500 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7
+GaspardTracker Trapezoid
+ X128_Y128= 14.346 -55.338 146.501 mm
+ X1_Y128= -14.263 -55.338 146.501 mm
+ X128_Y1= 48.726 -138.518 69.237 mm
+ X1_Y1=  -48.717 -138.519 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8
+GaspardTracker Trapezoid
+ X128_Y128= 49.274 -28.986 146.501 mm
+ X1_Y128= 29.045 -49.215 146.501 mm
+ X128_Y1= 132.401 -63.492 69.237 mm
+ X1_Y1=  63.500 -132.395 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8

Projects/GRIT_AGATA_GANIL/detector/small_grit_back.detector

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1
+Target
+ THICKNESS= 10 micrometer 
+ RADIUS= 7.5 mm 
+ MATERIAL= CD2
+ ANGLE= 0 deg
+ X= 0 mm 
+ Y= 0 mm 
+ Z= 0 mm
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias THIRDSTAGEFRONT
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias THIRDSTAGEBACK
+  Action= Replace
+  Value= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEFRONT
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEBACK
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEBARREL
+  Action= Replace
+  Value= 1 
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 End-Cap Back
+GaspardTracker Trapezoid
+ X128_Y128= 55.338 -14.346 -146.501 mm
+ X1_Y128= 55.338 14.346 -146.501 mm
+ X128_Y1= 138.519 -48.717 -69.236 mm
+ X1_Y1=  138.518 48.726 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2
+GaspardTracker Trapezoid
+ X128_Y128= 49.215 29.045 -146.501 mm
+ X1_Y128=  28.986 49.274 -146.501 mm
+ X128_Y1= 132.395 63.500 -69.236 mm
+ X1_Y1=  63.492 132.401 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3
+% GaspardTracker Trapezoid
+%  X128_Y128= 14.263 55.338 -146.501 mm
+%  X1_Y128= -14.346 55.338 -146.501 mm
+%  X128_Y1=  48.717 138.519 -69.236 mm
+%  X1_Y1=  -48.726 138.518 -69.237 mm
+%  FIRSTSTAGE= 1
+%  SECONDSTAGE= @SECONDSTAGEBACK 
+%  THIRDSTAGE= @THIRDSTAGEBACK
+%  EnergyResolution= 50
+%  VIS= all
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%4
+GaspardTracker Trapezoid
+ X128_Y128= -29.045 49.215 -146.501 mm
+ X1_Y128= -49.274 28.986 -146.501 mm
+ X128_Y1= -63.500 132.395 -69.236 mm
+ X1_Y1=  -132.401 63.492 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
+GaspardTracker Trapezoid
+ X128_Y128= -55.338 14.346 -146.501 mm
+ X1_Y128= -55.338 -14.346 -146.501 mm
+ X128_Y1= -138.519 48.717 -69.236 mm
+ X1_Y1=  -138.518 -48.726 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6
+GaspardTracker Trapezoid
+ X128_Y128= -49.215 -29.045 -146.501 mm
+ X1_Y128= -28.986 -49.274 -146.501 mm
+ X128_Y1= -132.395 -63.500 -69.236 mm
+ X1_Y1=  -63.492 -132.401 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7
+% GaspardTracker Trapezoid
+%  X128_Y128= -14.263 -55.338 -146.501 mm
+%  X1_Y128= 14.346 -55.338 -146.501 mm
+%  X128_Y1= -48.717 -138.519 -69.236 mm
+%  X1_Y1=  48.726 -138.518 -69.237 mm
+%  FIRSTSTAGE= 1
+%  SECONDSTAGE= @SECONDSTAGEBACK 
+%  THIRDSTAGE= @THIRDSTAGEBACK
+%  EnergyResolution= 50
+%  VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8
+GaspardTracker Trapezoid
+ X128_Y128= 29.045 -49.215 -146.501 mm
+ X1_Y128= 49.274 -28.986 -146.501 mm
+ X128_Y1= 63.500 -132.395 -69.236 mm
+ X1_Y1=  132.401 -63.492 -69.237 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEBACK 
+ THIRDSTAGE= @THIRDSTAGEBACK
+ EnergyResolution= 50
+ VIS= all
+

Projects/GRIT_AGATA_GANIL/detector/small_grit_front.detector

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1
+Target
+ THICKNESS= 10 micrometer 
+ RADIUS= 7.5 mm 
+ MATERIAL= CD2
+ ANGLE= 0 deg
+ X= 0 mm 
+ Y= 0 mm 
+ Z= 0 mm
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias THIRDSTAGEFRONT
+  Action= Replace
+  Value= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias THIRDSTAGEBACK
+  Action= Replace
+  Value= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEFRONT
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEBACK
+  Action= Replace
+  Value= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Alias SECONDSTAGEBARREL
+  Action= Replace
+  Value= 1 
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 Front-Cap Back
+GaspardTracker Trapezoid
+ X128_Y128= 55.338 14.346 146.501 mm
+ X1_Y128= 55.338 -14.346 146.501 mm
+ X128_Y1= 138.518 48.726 69.237 mm
+ X1_Y1=  138.519 -48.717 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2
+GaspardTracker Trapezoid
+ X128_Y128= 28.986 49.274 146.501 mm
+ X1_Y128= 49.215 29.045 146.501 mm
+ X128_Y1= 63.492 132.401 69.237 mm
+ X1_Y1=  132.395 63.500 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3
+% GaspardTracker Trapezoid
+%  X128_Y128= -14.346 55.338 146.501 mm
+%  X1_Y128= 14.263 55.338 146.501 mm
+%  X128_Y1= -48.726 138.518 69.237 mm
+%  X1_Y1=  48.717 138.519 69.236 mm
+%  FIRSTSTAGE= 1
+%  SECONDSTAGE= @SECONDSTAGEFRONT 
+%  THIRDSTAGE= @THIRDSTAGEFRONT
+%  EnergyResolution= 50
+%  VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%4
+GaspardTracker Trapezoid
+ X128_Y128= -49.274 28.986 146.501 mm
+ X1_Y128= -29.045 49.215 146.501 mm
+ X128_Y1= -132.401 63.492 69.237 mm
+ X1_Y1=  -63.500 132.395 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
+GaspardTracker Trapezoid
+ X128_Y128= -55.338 -14.346 146.501 mm
+ X1_Y128= -55.338 14.346 146.501 mm
+ X128_Y1= -138.518 -48.726 69.237 mm
+ X1_Y1=  -138.519 48.717 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6
+GaspardTracker Trapezoid
+ X128_Y128= -28.986 -49.274 146.501 mm
+ X1_Y128= -49.215 -29.045 146.501 mm
+ X128_Y1= -63.492 -132.401 69.237 mm
+ X1_Y1=  -132.395 -63.500 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7
+% GaspardTracker Trapezoid
+%  X128_Y128= 14.346 -55.338 146.501 mm
+%  X1_Y128= -14.263 -55.338 146.501 mm
+%  X128_Y1= 48.726 -138.518 69.237 mm
+%  X1_Y1=  -48.717 -138.519 69.236 mm
+%  FIRSTSTAGE= 1
+%  SECONDSTAGE= @SECONDSTAGEFRONT 
+%  THIRDSTAGE= @THIRDSTAGEFRONT
+%  EnergyResolution= 50
+%  VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8
+GaspardTracker Trapezoid
+ X128_Y128= 49.274 -28.986 146.501 mm
+ X1_Y128= 29.045 -49.215 146.501 mm
+ X128_Y1= 132.401 -63.492 69.237 mm
+ X1_Y1=  63.500 -132.395 69.236 mm
+ FIRSTSTAGE= 1
+ SECONDSTAGE= @SECONDSTAGEFRONT 
+ THIRDSTAGE= @THIRDSTAGEFRONT
+ EnergyResolution= 50
+ VIS= all
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%8

Projects/GRIT_AGATA_GANIL/reaction/132Sndp.reaction

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%  Reaction file for 60Fe(d,p)61Fe reaction  %%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%Beam energy given in MeV ; Excitation in MeV ; emmitance in rad
+Beam
+ Particle= 132Sn
+ Energy= 1320
+ SigmaEnergy= 0 MeV
+ SigmaThetaX= 0.42 deg
+ SigmaPhiY= 0.42 deg
+ SigmaX= 1.7 mm
+ SigmaY= 1.7 mm
+ MeanThetaX= 0.3 deg
+ MeanPhiY= 0 deg
+ MeanX= 1.0 mm
+ MeanY= 2.0 mm
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+TwoBodyReaction
+ Beam= 132Sn
+ Target= 2H
+ Light= 1H
+ Heavy= 133Sn
+ ExcitationEnergyLight= 0.0 MeV
+ ExcitationEnergyHeavy= 0.0 MeV
+ CrossSectionPath= 132Sndp_10A_MeV_2f7_ZR_FRC.lis CS
+ HalfOpenAngleMin= 0 deg
+ HalfOpenAngleMax= 180 deg
+ ShootLight= 1
+ ShootHeavy= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+

Projects/GRIT_AGATA_GANIL/reaction/14Odp.reaction

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%  Reaction file for 60Fe(d,p)61Fe reaction  %%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%Beam energy given in MeV ; Excitation in MeV ; emmitance in rad
+Beam
+ Particle= 14O
+ Energy= 140
+ SigmaEnergy= 0 MeV
+ SigmaThetaX= 0.42 deg
+ SigmaPhiY= 0.42 deg
+ SigmaX= 1.7 mm
+ SigmaY= 1.7 mm
+ MeanThetaX= 0.3 deg
+ MeanPhiY= 0 deg
+ MeanX= 1.0 mm
+ MeanY= 2.0 mm
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+TwoBodyReaction
+ Beam= 14O
+ Target= 2H
+ Light= 1H
+ Heavy= 15O
+ ExcitationEnergyLight= 0.0 MeV
+ ExcitationEnergyHeavy= 0.0 MeV
+ CrossSectionPath= 132Sndp_10A_MeV_2f7_ZR_FRC.lis CS
+ HalfOpenAngleMin= 0 deg
+ HalfOpenAngleMax= 180 deg
+ ShootLight= 1
+ ShootHeavy= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+

Projects/GRIT_AGATA_GANIL/reaction/14Opd.reaction

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%  Reaction file for 60Fe(d,p)61Fe reaction  %%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%Beam energy given in MeV ; Excitation in MeV ; emmitance in rad
+Beam
+ Particle= 14O
+ Energy= 1400
+ SigmaEnergy= 0 MeV
+ SigmaThetaX= 0.42 deg
+ SigmaPhiY= 0.42 deg
+ SigmaX= 1.7 mm
+ SigmaY= 1.7 mm
+ MeanThetaX= 0.3 deg
+ MeanPhiY= 0 deg
+ MeanX= 1.0 mm
+ MeanY= 2.0 mm
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+TwoBodyReaction
+ Beam= 14O
+ Target= 1H
+ Light= 2H
+ Heavy= 13O
+ ExcitationEnergyLight= 0.0 MeV
+ ExcitationEnergyHeavy= 0.0 MeV
+ CrossSectionPath= 132Sndp_10A_MeV_2f7_ZR_FRC.lis CS
+ HalfOpenAngleMin= 0 deg
+ HalfOpenAngleMax= 180 deg
+ ShootLight= 1
+ ShootHeavy= 1
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+