From 770d1cadd31dde82b6dcb06d942e64ab6b326aa6 Mon Sep 17 00:00:00 2001
From: adrien-matta <a.matta@surrey.ac.uk>
Date: Mon, 13 Jul 2015 17:33:27 -0700
Subject: [PATCH] * Adding S1554 analysis
---
NPAnalysis/S1554/22.si | 181 +++++++
NPAnalysis/S1554/22.sie | 181 +++++++
NPAnalysis/S1554/Analysis.cxx | 161 ++++++
NPAnalysis/S1554/Analysis.h | 87 +++
NPAnalysis/S1554/CMakeLists.txt | 31 ++
NPAnalysis/S1554/RunToTreat.txt | 6 +
NPAnalysis/S1554/SharcCalibration.txt | 496 ++++++++++++++++++
NPAnalysis/S1554/SharcEff.root | Bin 0 -> 5788 bytes
NPAnalysis/S1554/Show.C | 59 +++
.../S1554/macro/Energy/EnergyCalibrator.cxx | 64 +++
10 files changed, 1266 insertions(+)
create mode 100644 NPAnalysis/S1554/22.si
create mode 100644 NPAnalysis/S1554/22.sie
create mode 100644 NPAnalysis/S1554/Analysis.cxx
create mode 100644 NPAnalysis/S1554/Analysis.h
create mode 100644 NPAnalysis/S1554/CMakeLists.txt
create mode 100644 NPAnalysis/S1554/RunToTreat.txt
create mode 100644 NPAnalysis/S1554/SharcCalibration.txt
create mode 100644 NPAnalysis/S1554/SharcEff.root
create mode 100644 NPAnalysis/S1554/Show.C
create mode 100644 NPAnalysis/S1554/macro/Energy/EnergyCalibrator.cxx
diff --git a/NPAnalysis/S1554/22.si b/NPAnalysis/S1554/22.si
new file mode 100644
index 000000000..6d9163232
--- /dev/null
+++ b/NPAnalysis/S1554/22.si
@@ -0,0 +1,181 @@
+ 0.00 8.082999E+01
+ 0.96 8.013461E+01
+ 1.92 7.807824E+01
+ 2.88 7.474845E+01
+ 3.84 7.028513E+01
+ 4.80 6.487207E+01
+ 5.76 5.872597E+01
+ 6.72 5.208399E+01
+ 7.68 4.519057E+01
+ 8.64 3.828463E+01
+ 9.60 3.158781E+01
+ 10.56 2.529468E+01
+ 11.52 1.956526E+01
+ 12.48 1.452010E+01
+ 13.44 1.023819E+01
+ 14.41 6.757372E+00
+ 15.37 4.076943E+00
+ 16.33 2.162190E+00
+ 17.29 9.501766E-01
+ 18.25 3.563639E-01
+ 19.21 2.815332E-01
+ 20.18 6.185614E-01
+ 21.14 1.258682E+00
+ 22.10 2.096945E+00
+ 23.06 3.036686E+00
+ 24.03 3.992911E+00
+ 24.99 4.894570E+00
+ 25.95 5.685780E+00
+ 26.92 6.326100E+00
+ 27.88 6.789981E+00
+ 28.85 7.065583E+00
+ 29.81 7.153095E+00
+ 30.78 7.062749E+00
+ 31.74 6.812668E+00
+ 32.71 6.426687E+00
+ 33.68 5.932258E+00
+ 34.64 5.358531E+00
+ 35.61 4.734672E+00
+ 36.58 4.088454E+00
+ 37.54 3.445147E+00
+ 38.51 2.826715E+00
+ 39.48 2.251286E+00
+ 40.45 1.732889E+00
+ 41.42 1.281428E+00
+ 42.39 9.028334E-01
+ 43.36 5.993811E-01
+ 44.33 3.701125E-01
+ 45.30 2.113377E-01
+ 46.27 1.171760E-01
+ 47.25 8.010797E-02
+ 48.22 9.151219E-02
+ 49.19 1.421664E-01
+ 50.17 2.226974E-01
+ 51.14 3.239696E-01
+ 52.12 4.374054E-01
+ 53.09 5.552348E-01
+ 54.07 6.706767E-01
+ 55.04 7.780529E-01
+ 56.02 8.728439E-01
+ 57.00 9.516905E-01
+ 57.98 1.012351E+00
+ 58.95 1.053624E+00
+ 59.93 1.075241E+00
+ 60.91 1.077740E+00
+ 61.89 1.062333E+00
+ 62.87 1.030760E+00
+ 63.86 9.851529E-01
+ 64.84 9.279005E-01
+ 65.82 8.615252E-01
+ 66.80 7.885734E-01
+ 67.79 7.115213E-01
+ 68.77 6.326960E-01
+ 69.76 5.542144E-01
+ 70.75 4.779374E-01
+ 71.73 4.054407E-01
+ 72.72 3.379991E-01
+ 73.71 2.765836E-01
+ 74.70 2.218697E-01
+ 75.69 1.742540E-01
+ 76.68 1.338789E-01
+ 77.67 1.006615E-01
+ 78.66 7.432630E-02
+ 79.65 5.443989E-02
+ 80.64 4.044575E-02
+ 81.64 3.169825E-02
+ 82.63 2.749481E-02
+ 83.63 2.710534E-02
+ 84.62 2.979840E-02
+ 85.62 3.486372E-02
+ 86.62 4.163078E-02
+ 87.61 4.948352E-02
+ 88.61 5.787110E-02
+ 89.61 6.631512E-02
+ 90.61 7.441344E-02
+ 91.61 8.184095E-02
+ 92.61 8.834787E-02
+ 93.62 9.375581E-02
+ 94.62 9.795218E-02
+ 95.62 1.008834E-01
+ 96.63 1.025469E-01
+ 97.63 1.029835E-01
+ 98.64 1.022685E-01
+ 99.65 1.005040E-01
+ 100.65 9.781104E-02
+ 101.66 9.432246E-02
+ 102.67 9.017676E-02
+ 103.68 8.551262E-02
+ 104.69 8.046432E-02
+ 105.70 7.515820E-02
+ 106.71 6.971000E-02
+ 107.73 6.422301E-02
+ 108.74 5.878703E-02
+ 109.75 5.347803E-02
+ 110.77 4.835837E-02
+ 111.78 4.347747E-02
+ 112.80 3.887291E-02
+ 113.82 3.457166E-02
+ 114.83 3.059159E-02
+ 115.85 2.694294E-02
+ 116.87 2.362977E-02
+ 117.89 2.065139E-02
+ 118.91 1.800357E-02
+ 119.93 1.567961E-02
+ 120.95 1.367119E-02
+ 121.97 1.196908E-02
+ 123.00 1.056353E-02
+ 124.02 9.444573E-03
+ 125.04 8.602090E-03
+ 126.07 8.025798E-03
+ 127.10 7.705073E-03
+ 128.12 7.628729E-03
+ 129.15 7.784749E-03
+ 130.18 8.160007E-03
+ 131.20 8.740017E-03
+ 132.23 9.508727E-03
+ 133.26 1.044837E-02
+ 134.29 1.153943E-02
+ 135.32 1.276065E-02
+ 136.35 1.408918E-02
+ 137.38 1.550076E-02
+ 138.42 1.697010E-02
+ 139.45 1.847118E-02
+ 140.48 1.997771E-02
+ 141.51 2.146364E-02
+ 142.55 2.290360E-02
+ 143.58 2.427345E-02
+ 144.62 2.555071E-02
+ 145.65 2.671506E-02
+ 146.69 2.774869E-02
+ 147.72 2.863666E-02
+ 148.76 2.936713E-02
+ 149.80 2.993156E-02
+ 150.84 3.032475E-02
+ 151.87 3.054490E-02
+ 152.91 3.059347E-02
+ 153.95 3.047504E-02
+ 154.99 3.019707E-02
+ 156.03 2.976962E-02
+ 157.07 2.920499E-02
+ 158.11 2.851733E-02
+ 159.15 2.772229E-02
+ 160.19 2.683658E-02
+ 161.23 2.587755E-02
+ 162.27 2.486286E-02
+ 163.31 2.381009E-02
+ 164.36 2.273641E-02
+ 165.40 2.165836E-02
+ 166.44 2.059152E-02
+ 167.48 1.955044E-02
+ 168.53 1.854842E-02
+ 169.57 1.759748E-02
+ 170.61 1.670831E-02
+ 171.65 1.589026E-02
+ 172.70 1.515139E-02
+ 173.74 1.449850E-02
+ 174.78 1.393724E-02
+ 175.83 1.347220E-02
+ 176.87 1.310696E-02
+ 177.91 1.284422E-02
+ 178.96 1.268585E-02
+ 180.00 1.263293E-02
diff --git a/NPAnalysis/S1554/22.sie b/NPAnalysis/S1554/22.sie
new file mode 100644
index 000000000..213132212
--- /dev/null
+++ b/NPAnalysis/S1554/22.sie
@@ -0,0 +1,181 @@
+ 0.00 1.108964E+02
+ 0.95 1.110523E+02
+ 1.91 1.114927E+02
+ 2.86 1.121379E+02
+ 3.82 1.128632E+02
+ 4.77 1.135110E+02
+ 5.73 1.139048E+02
+ 6.68 1.138653E+02
+ 7.64 1.132258E+02
+ 8.59 1.118459E+02
+ 9.55 1.096229E+02
+ 10.50 1.064992E+02
+ 11.46 1.024653E+02
+ 12.41 9.755997E+01
+ 13.37 9.186575E+01
+ 14.32 8.550242E+01
+ 15.28 7.861815E+01
+ 16.24 7.137960E+01
+ 17.19 6.396181E+01
+ 18.15 5.653859E+01
+ 19.11 4.927403E+01
+ 20.06 4.231555E+01
+ 21.02 3.578859E+01
+ 21.98 2.979331E+01
+ 22.94 2.440281E+01
+ 23.89 1.966312E+01
+ 24.85 1.559439E+01
+ 25.81 1.219322E+01
+ 26.77 9.435687E+00
+ 27.73 7.280947E+00
+ 28.69 5.675011E+00
+ 29.65 4.554579E+00
+ 30.61 3.850707E+00
+ 31.57 3.492168E+00
+ 32.53 3.408406E+00
+ 33.49 3.532023E+00
+ 34.45 3.800745E+00
+ 35.42 4.158881E+00
+ 36.38 4.558288E+00
+ 37.34 4.958884E+00
+ 38.31 5.328754E+00
+ 39.27 5.643927E+00
+ 40.24 5.887884E+00
+ 41.20 6.050856E+00
+ 42.17 6.128991E+00
+ 43.13 6.123433E+00
+ 44.10 6.039375E+00
+ 45.07 5.885112E+00
+ 46.03 5.671148E+00
+ 47.00 5.409363E+00
+ 47.97 5.112272E+00
+ 48.94 4.792380E+00
+ 49.91 4.461646E+00
+ 50.88 4.131056E+00
+ 51.85 3.810298E+00
+ 52.82 3.507552E+00
+ 53.80 3.229356E+00
+ 54.77 2.980570E+00
+ 55.74 2.764411E+00
+ 56.72 2.582539E+00
+ 57.69 2.435203E+00
+ 58.67 2.321405E+00
+ 59.64 2.239105E+00
+ 60.62 2.185416E+00
+ 61.60 2.156815E+00
+ 62.57 2.149341E+00
+ 63.55 2.158783E+00
+ 64.53 2.180848E+00
+ 65.51 2.211315E+00
+ 66.49 2.246159E+00
+ 67.48 2.281662E+00
+ 68.46 2.314493E+00
+ 69.44 2.341775E+00
+ 70.43 2.361118E+00
+ 71.41 2.370647E+00
+ 72.40 2.369000E+00
+ 73.38 2.355315E+00
+ 74.37 2.329203E+00
+ 75.36 2.290707E+00
+ 76.34 2.240255E+00
+ 77.33 2.178606E+00
+ 78.32 2.106789E+00
+ 79.31 2.026050E+00
+ 80.31 1.937788E+00
+ 81.30 1.843508E+00
+ 82.29 1.744760E+00
+ 83.29 1.643104E+00
+ 84.28 1.540061E+00
+ 85.28 1.437087E+00
+ 86.27 1.335536E+00
+ 87.27 1.236642E+00
+ 88.27 1.141505E+00
+ 89.27 1.051069E+00
+ 90.27 9.661266E-01
+ 91.27 8.873059E-01
+ 92.27 8.150780E-01
+ 93.27 7.497588E-01
+ 94.27 6.915178E-01
+ 95.28 6.403877E-01
+ 96.28 5.962774E-01
+ 97.29 5.589851E-01
+ 98.30 5.282139E-01
+ 99.30 5.035865E-01
+ 100.31 4.846610E-01
+ 101.32 4.709457E-01
+ 102.33 4.619132E-01
+ 103.34 4.570140E-01
+ 104.35 4.556883E-01
+ 105.37 4.573767E-01
+ 106.38 4.615300E-01
+ 107.39 4.676168E-01
+ 108.41 4.751303E-01
+ 109.42 4.835938E-01
+ 110.44 4.925649E-01
+ 111.46 5.016383E-01
+ 112.48 5.104484E-01
+ 113.49 5.186702E-01
+ 114.51 5.260203E-01
+ 115.53 5.322562E-01
+ 116.56 5.371761E-01
+ 117.58 5.406175E-01
+ 118.60 5.424558E-01
+ 119.62 5.426026E-01
+ 120.65 5.410034E-01
+ 121.67 5.376358E-01
+ 122.70 5.325068E-01
+ 123.73 5.256507E-01
+ 124.76 5.171268E-01
+ 125.78 5.070166E-01
+ 126.81 4.954218E-01
+ 127.84 4.824615E-01
+ 128.87 4.682707E-01
+ 129.90 4.529970E-01
+ 130.94 4.367990E-01
+ 131.97 4.198441E-01
+ 133.00 4.023061E-01
+ 134.04 3.843630E-01
+ 135.07 3.661953E-01
+ 136.11 3.479834E-01
+ 137.14 3.299062E-01
+ 138.18 3.121384E-01
+ 139.22 2.948494E-01
+ 140.25 2.782006E-01
+ 141.29 2.623445E-01
+ 142.33 2.474222E-01
+ 143.37 2.335628E-01
+ 144.41 2.208811E-01
+ 145.45 2.094772E-01
+ 146.49 1.994351E-01
+ 147.53 1.908216E-01
+ 148.58 1.836864E-01
+ 149.62 1.780610E-01
+ 150.66 1.739592E-01
+ 151.71 1.713764E-01
+ 152.75 1.702908E-01
+ 153.79 1.706629E-01
+ 154.84 1.724373E-01
+ 155.88 1.755426E-01
+ 156.93 1.798933E-01
+ 157.98 1.853908E-01
+ 159.02 1.919246E-01
+ 160.07 1.993742E-01
+ 161.12 2.076104E-01
+ 162.16 2.164973E-01
+ 163.21 2.258939E-01
+ 164.26 2.356557E-01
+ 165.31 2.456370E-01
+ 166.36 2.556920E-01
+ 167.40 2.656768E-01
+ 168.45 2.754513E-01
+ 169.50 2.848803E-01
+ 170.55 2.938353E-01
+ 171.60 3.021957E-01
+ 172.65 3.098502E-01
+ 173.70 3.166977E-01
+ 174.75 3.226486E-01
+ 175.80 3.276257E-01
+ 176.85 3.315647E-01
+ 177.90 3.344148E-01
+ 178.95 3.361396E-01
+ 180.00 3.367169E-01
diff --git a/NPAnalysis/S1554/Analysis.cxx b/NPAnalysis/S1554/Analysis.cxx
new file mode 100644
index 000000000..014f58c8f
--- /dev/null
+++ b/NPAnalysis/S1554/Analysis.cxx
@@ -0,0 +1,161 @@
+/*****************************************************************************
+ * 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 2015 *
+ * 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"NPOptionManager.h"
+////////////////////////////////////////////////////////////////////////////////
+Analysis::Analysis(){
+}
+////////////////////////////////////////////////////////////////////////////////
+Analysis::~Analysis(){
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::Init(){
+ InitOutputBranch();
+ InitInputBranch();
+
+ Sharc = (TSharcPhysics*) m_DetectorManager -> GetDetector("Sharc");
+ LightCD2 = EnergyLoss("proton_CD2.G4table","G4Table",100 );
+ LightSi = EnergyLoss("proton_Si.G4table","G4Table",100);
+ BeamCD2 = EnergyLoss("Si28[0.0]_CD2.G4table","G4Table",100);
+ myReaction = new NPL::Reaction();
+ myReaction->ReadConfigurationFile(NPOptionManager::getInstance()->GetReactionFile());
+ TargetThickness = m_DetectorManager->GetTargetThickness()*micrometer;
+ OriginalBeamEnergy = myReaction->GetBeamEnergy();
+ Rand = TRandom3();
+ DetectorNumber = 0 ;
+ ThetaNormalTarget = 0 ;
+ ThetaM2Surface = 0;
+ Si_E_M2 = 0 ;
+ CsI_E_M2 = 0 ;
+ Energy = 0;
+ E_M2 = 0;
+
+ ThetaSharcSurface = 0;
+ X_Sharc = 0 ;
+ Y_Sharc = 0 ;
+ Z_Sharc = 0 ;
+ Si_E_Sharc = 0 ;
+ E_Sharc = 0;
+ Si_X_Sharc = 0;
+ Si_Y_Sharc = 0;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::TreatEvent(){
+ // Reinitiate calculated variable
+ ReInitValue();
+ double XTarget = 0;
+ double YTarget = 0;
+ TVector3 BeamDirection = TVector3(0,0,1);
+ double BeamEnergy = BeamCD2.Slow(OriginalBeamEnergy,TargetThickness*0.5,0);
+ myReaction->SetBeamEnergy(BeamEnergy);
+ ////////////////////////////////////////////////////////////////////////////
+ ////////////////////////////////////////////////////////////////////////////
+ //////////////////////////// LOOP on Sharc//////////////////
+ if(Sharc->Strip_E.size()>0){
+ /************************************************/
+ // Part 1 : Impact Angle
+ ThetaSharcSurface = 0;
+ ThetaNormalTarget = 0;
+ TVector3 HitDirection = Sharc -> GetPositionOfInteraction(0);
+ ThetaLab = HitDirection.Angle( BeamDirection );
+ ThetaNormalTarget = HitDirection.Angle( TVector3(0,0,1) ) ;
+
+ /************************************************/
+
+ /************************************************/
+ // Part 2 : Impact Energy
+
+ Energy = ELab = 0;
+ if(Sharc->PAD_E[0]>0){
+ Energy = Sharc->PAD_E[0];
+ }
+
+ Energy += Sharc->Strip_E[0];
+ // Target Correction
+ ELab = LightCD2.EvaluateInitialEnergy( Energy ,TargetThickness*0.5, ThetaNormalTarget);
+ /************************************************/
+
+ /************************************************/
+ // Part 3 : Excitation Energy Calculation
+ Ex = myReaction -> ReconstructRelativistic( ELab , ThetaLab );
+ /************************************************/
+
+ /************************************************/
+ // Part 4 : Theta CM Calculation
+ ThetaCM = myReaction -> EnergyLabToThetaCM( ELab , ThetaLab)/deg;
+ ThetaLab=ThetaLab/deg;
+ ThetaLab=Rand.Uniform(ThetaLab-0.5,ThetaLab+0.5);
+ /************************************************/
+ }//end loop Sharc
+
+}
+
+////////////////////////////////////////////////////////////////////////////////
+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");
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::InitInputBranch(){
+}
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::ReInitValue(){
+ Ex = -1000 ;
+ ELab = -1000;
+ ThetaLab = -1000;
+ ThetaCM = -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{
+ public:
+ proxy(){
+ NPL::AnalysisFactory::getInstance()->SetConstructor(Analysis::Construct);
+ }
+};
+
+proxy p;
+}
+
diff --git a/NPAnalysis/S1554/Analysis.h b/NPAnalysis/S1554/Analysis.h
new file mode 100644
index 000000000..5b9ce24e5
--- /dev/null
+++ b/NPAnalysis/S1554/Analysis.h
@@ -0,0 +1,87 @@
+#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"NPVAnalysis.h"
+#include"NPEnergyLoss.h"
+#include"NPReaction.h"
+#include"RootOutput.h"
+#include"RootInput.h"
+#include "TSharcPhysics.h"
+#include "TInitialConditions.h"
+#include <TRandom3.h>
+#include <TVector3.h>
+#include <TMath.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 ELab;
+ double ThetaLab;
+ double ThetaCM;
+ NPL::Reaction* myReaction;
+TInitialConditions* myInit ;
+ // Energy loss table: the G4Table are generated by the simulation
+ EnergyLoss LightCD2;
+ EnergyLoss LightAl;
+ EnergyLoss LightSi;
+ EnergyLoss BeamCD2;
+ TVector3 BeamImpact;
+
+ double TargetThickness ;
+ // Beam Energy
+ double OriginalBeamEnergy ; // AMEV
+ // intermediate variable
+ TRandom3 Rand ;
+ int DetectorNumber ;
+ double ThetaNormalTarget;
+ double ThetaM2Surface ;
+ double Si_E_M2 ;
+ double CsI_E_M2 ;
+ double Energy ;
+ double E_M2 ;
+
+ double ThetaSharcSurface ;
+ double X_Sharc ;
+ double Y_Sharc ;
+ double Z_Sharc ;
+ double Si_E_Sharc ;
+ double E_Sharc ;
+ double Si_X_Sharc ;
+ double Si_Y_Sharc ;
+ TSharcPhysics* Sharc;
+};
+#endif
diff --git a/NPAnalysis/S1554/CMakeLists.txt b/NPAnalysis/S1554/CMakeLists.txt
new file mode 100644
index 000000000..1d6a34417
--- /dev/null
+++ b/NPAnalysis/S1554/CMakeLists.txt
@@ -0,0 +1,31 @@
+cmake_minimum_required (VERSION 2.8)
+#Finding NPTool
+set(NPTOOL "$ENV{NPTOOL}")
+set(NPLIB "${NPTOOL}/NPLib")
+set(NPTOOL_INCLUDE_DIR "${NPLIB}/include")
+set(NPTOOL_LIB_DIR "${NPLIB}/lib")
+
+include("${NPLIB}/FindROOT.cmake")
+
+project (NPAnalysis)
+set(CMAKE_BUILD_TYPE Release)
+# Add root to the link and include directories
+include_directories( ${ROOT_INCLUDE_DIR})
+link_directories( ${ROOT_LIBRARY_DIR})
+include_directories( ${NPTOOL_INCLUDE_DIR})
+link_directories( ${NPTOOL_LIB_DIR})
+
+# Get the compilator flag from root to assure consistancy
+EXEC_PROGRAM(${ROOT_CONFIG_EXECUTABLE}
+ ARGS "--cflags"
+ OUTPUT_VARIABLE root_cflags )
+
+set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} ${root_cflags}")
+
+# If the compiler is Clang, silence the unrecognised flags
+if(${CMAKE_CXX_COMPILER_ID} MATCHES ".*Clang.*")
+ set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Qunused-arguments -undefined dynamic_lookup")
+endif()
+
+add_library(NPAnalysis SHARED Analysis.cxx)
+target_link_libraries(NPAnalysis ${ROOT_LIBRARIES} -L${NPLIB}/lib -lNPCore -lNPPhysics)
diff --git a/NPAnalysis/S1554/RunToTreat.txt b/NPAnalysis/S1554/RunToTreat.txt
new file mode 100644
index 000000000..cd169430b
--- /dev/null
+++ b/NPAnalysis/S1554/RunToTreat.txt
@@ -0,0 +1,6 @@
+TTreeName
+ DataS1554
+RootFileName
+ ~/Desktop/Unpacker/data3507*
+
+
diff --git a/NPAnalysis/S1554/SharcCalibration.txt b/NPAnalysis/S1554/SharcCalibration.txt
new file mode 100644
index 000000000..da9e75e84
--- /dev/null
+++ b/NPAnalysis/S1554/SharcCalibration.txt
@@ -0,0 +1,496 @@
+SHARC_D1_STRIP_FRONT1_E 1 9.71809e-05
+SHARC_D1_STRIP_FRONT2_E 0.164687 2.26585e-05
+SHARC_D1_STRIP_FRONT3_E 0.174147 2.28104e-05
+SHARC_D1_STRIP_FRONT4_E 0.162916 2.26244e-05
+SHARC_D1_STRIP_FRONT5_E 0.163951 2.18268e-05
+SHARC_D1_STRIP_FRONT6_E 0.166109 2.20008e-05
+SHARC_D1_STRIP_FRONT8_E 0.160454 2.15271e-05
+SHARC_D1_STRIP_FRONT9_E 0.165885 2.19531e-05
+SHARC_D1_STRIP_FRONT10_E 0.162686 2.20978e-05
+SHARC_D1_STRIP_FRONT11_E 0.174349 2.19763e-05
+SHARC_D1_STRIP_FRONT12_E 0.160229 2.2092e-05
+SHARC_D1_STRIP_FRONT14_E 0.168189 2.17203e-05
+SHARC_D1_STRIP_FRONT15_E 0.158883 2.13092e-05
+SHARC_D1_STRIP_FRONT16_E 0.145459 2.17219e-05
+SHARC_D1_STRIP_BACK1_E 0.153838 2.60365e-05
+SHARC_D1_STRIP_BACK2_E 0.151547 2.25412e-05
+SHARC_D1_STRIP_BACK3_E 0.159238 2.24302e-05
+SHARC_D1_STRIP_BACK4_E 0.155658 2.24398e-05
+SHARC_D1_STRIP_BACK5_E 0.15404 2.23263e-05
+SHARC_D1_STRIP_BACK6_E 0.161111 2.23206e-05
+SHARC_D1_STRIP_BACK7_E 0.158584 2.24573e-05
+SHARC_D1_STRIP_BACK8_E 0.160069 2.21756e-05
+SHARC_D1_STRIP_BACK9_E 0.154648 2.2053e-05
+SHARC_D1_STRIP_BACK10_E 0.152995 2.22484e-05
+SHARC_D1_STRIP_BACK11_E 0.15407 2.23135e-05
+SHARC_D1_STRIP_BACK12_E 0.150986 2.24622e-05
+SHARC_D1_STRIP_BACK13_E 0.164412 2.23237e-05
+SHARC_D1_STRIP_BACK14_E 0.157589 2.17609e-05
+SHARC_D1_STRIP_BACK15_E 0.154085 2.22351e-05
+SHARC_D1_STRIP_BACK16_E 0.152418 2.2914e-05
+SHARC_D1_STRIP_BACK17_E 0.145236 2.24632e-05
+SHARC_D1_STRIP_BACK18_E 0.149564 2.1789e-05
+SHARC_D1_STRIP_BACK19_E 0.160901 2.21543e-05
+SHARC_D1_STRIP_BACK20_E 0.161678 2.20323e-05
+SHARC_D1_STRIP_BACK21_E 0.159755 2.24461e-05
+SHARC_D1_STRIP_BACK22_E 0.149375 2.22273e-05
+SHARC_D1_STRIP_BACK23_E 0.158633 2.27563e-05
+SHARC_D1_STRIP_BACK24_E 0.146264 2.61009e-05
+SHARC_D2_STRIP_FRONT2_E 0.181992 2.15866e-05
+SHARC_D2_STRIP_FRONT3_E 0.173214 2.14823e-05
+SHARC_D2_STRIP_FRONT4_E 0.174668 2.14928e-05
+SHARC_D2_STRIP_FRONT5_E 0.166655 2.15335e-05
+SHARC_D2_STRIP_FRONT6_E 0.15884 2.16476e-05
+SHARC_D2_STRIP_FRONT7_E 0.168056 2.15086e-05
+SHARC_D2_STRIP_FRONT8_E 0.162027 2.13847e-05
+SHARC_D2_STRIP_FRONT9_E 0.157742 2.12032e-05
+SHARC_D2_STRIP_FRONT10_E 0.153569 2.09853e-05
+SHARC_D2_STRIP_FRONT11_E 0.153675 2.06502e-05
+SHARC_D2_STRIP_FRONT12_E 0.154205 2.1239e-05
+SHARC_D2_STRIP_FRONT13_E 0.147765 2.12137e-05
+SHARC_D2_STRIP_FRONT14_E 0.141683 2.10439e-05
+SHARC_D2_STRIP_FRONT15_E 0.153863 2.12348e-05
+SHARC_D2_STRIP_FRONT16_E 0.154194 2.11789e-05
+SHARC_D2_STRIP_BACK2_E 0.161791 2.16439e-05
+SHARC_D2_STRIP_BACK3_E 0.158998 2.13021e-05
+SHARC_D2_STRIP_BACK4_E 0.155088 2.12472e-05
+SHARC_D2_STRIP_BACK5_E 0.151944 2.14828e-05
+SHARC_D2_STRIP_BACK6_E 0.155744 2.11617e-05
+SHARC_D2_STRIP_BACK7_E 0.159097 2.12061e-05
+SHARC_D2_STRIP_BACK8_E 0.156703 2.11961e-05
+SHARC_D2_STRIP_BACK9_E 0.167073 2.14596e-05
+SHARC_D2_STRIP_BACK10_E 0.15311 2.14144e-05
+SHARC_D2_STRIP_BACK11_E 0.150879 2.17228e-05
+SHARC_D2_STRIP_BACK12_E 0.148134 2.13678e-05
+SHARC_D2_STRIP_BACK13_E 0.150769 2.09629e-05
+SHARC_D2_STRIP_BACK14_E 0.158156 2.15112e-05
+SHARC_D2_STRIP_BACK15_E 0.166186 2.15378e-05
+SHARC_D2_STRIP_BACK16_E 0.164185 2.19222e-05
+SHARC_D2_STRIP_BACK17_E 0.149017 2.15916e-05
+SHARC_D2_STRIP_BACK18_E 0.15531 2.14079e-05
+SHARC_D2_STRIP_BACK19_E 0.158552 2.1928e-05
+SHARC_D2_STRIP_BACK20_E 0.156214 2.16906e-05
+SHARC_D2_STRIP_BACK21_E 0.156121 2.18188e-05
+SHARC_D2_STRIP_BACK22_E 0.161442 2.16562e-05
+SHARC_D2_STRIP_BACK23_E 0.159057 2.10402e-05
+SHARC_D2_STRIP_BACK24_E 0.133271 2.23707e-05
+SHARC_D3_STRIP_FRONT1_E 0.170498 2.17833e-05
+SHARC_D3_STRIP_FRONT2_E 0.177864 2.16078e-05
+SHARC_D3_STRIP_FRONT3_E 0.166926 2.13756e-05
+SHARC_D3_STRIP_FRONT4_E 0.161308 2.14367e-05
+SHARC_D3_STRIP_FRONT5_E 0.165043 2.12113e-05
+SHARC_D3_STRIP_FRONT6_E 0.152747 2.14685e-05
+SHARC_D3_STRIP_FRONT7_E 0.148567 2.15193e-05
+SHARC_D3_STRIP_FRONT8_E 0.146428 2.15751e-05
+SHARC_D3_STRIP_FRONT9_E 0.157735 2.19116e-05
+SHARC_D3_STRIP_FRONT10_E 0.148186 2.1768e-05
+SHARC_D3_STRIP_FRONT11_E 0.143739 2.15243e-05
+SHARC_D3_STRIP_FRONT12_E 0.141352 2.15329e-05
+SHARC_D3_STRIP_FRONT13_E 0.137131 2.10355e-05
+SHARC_D3_STRIP_FRONT14_E 0.131506 2.09502e-05
+SHARC_D3_STRIP_FRONT15_E 0.146444 2.13021e-05
+SHARC_D3_STRIP_FRONT16_E 0.15261 2.10132e-05
+SHARC_D3_STRIP_BACK1_E 0.141785 2.24297e-05
+SHARC_D3_STRIP_BACK2_E 0.169104 2.15242e-05
+SHARC_D3_STRIP_BACK3_E 0.169554 2.12474e-05
+SHARC_D3_STRIP_BACK4_E 0.165779 2.13588e-05
+SHARC_D3_STRIP_BACK5_E 0.169948 2.16378e-05
+SHARC_D3_STRIP_BACK6_E 0.173751 2.09587e-05
+SHARC_D3_STRIP_BACK7_E 0.171863 2.13399e-05
+SHARC_D3_STRIP_BACK8_E 0.168346 2.15308e-05
+SHARC_D3_STRIP_BACK9_E 0.163165 2.11919e-05
+SHARC_D3_STRIP_BACK10_E 0.176397 2.04944e-05
+SHARC_D3_STRIP_BACK11_E 0.174951 2.09201e-05
+SHARC_D3_STRIP_BACK12_E 0.167448 2.0926e-05
+SHARC_D3_STRIP_BACK13_E 0.168121 2.103e-05
+SHARC_D3_STRIP_BACK14_E 0.168596 2.09978e-05
+SHARC_D3_STRIP_BACK15_E 0.175062 2.10524e-05
+SHARC_D3_STRIP_BACK16_E 0.166321 2.13564e-05
+SHARC_D3_STRIP_BACK17_E 0.162009 2.11844e-05
+SHARC_D3_STRIP_BACK18_E 0.171647 2.13428e-05
+SHARC_D3_STRIP_BACK19_E 0.165273 2.09417e-05
+SHARC_D3_STRIP_BACK20_E 0.178293 2.08393e-05
+SHARC_D3_STRIP_BACK21_E 0.16557 2.10342e-05
+SHARC_D3_STRIP_BACK22_E 0.166472 2.13356e-05
+SHARC_D3_STRIP_BACK23_E 0.178184 2.11782e-05
+SHARC_D4_STRIP_FRONT1_E 0.143893 2.15775e-05
+SHARC_D4_STRIP_FRONT2_E 0.151215 2.19168e-05
+SHARC_D4_STRIP_FRONT3_E 0.13701 2.16589e-05
+SHARC_D4_STRIP_FRONT4_E 0.150528 2.19227e-05
+SHARC_D4_STRIP_FRONT5_E 0.139879 2.20799e-05
+SHARC_D4_STRIP_FRONT6_E 0.132977 2.12784e-05
+SHARC_D4_STRIP_FRONT7_E 0.125682 2.16596e-05
+SHARC_D4_STRIP_FRONT12_E 0.117107 2.16223e-05
+SHARC_D4_STRIP_FRONT13_E 0.115119 2.11676e-05
+SHARC_D4_STRIP_FRONT14_E 0.11296 2.16328e-05
+SHARC_D4_STRIP_FRONT15_E 0.113554 2.12594e-05
+SHARC_D4_STRIP_FRONT16_E 0.120279 2.13435e-05
+SHARC_D4_STRIP_BACK1_E 0.116399 2.21512e-05
+SHARC_D4_STRIP_BACK2_E 0.153885 2.28043e-05
+SHARC_D4_STRIP_BACK3_E 0.116021 2.25722e-05
+SHARC_D4_STRIP_BACK4_E 0.11138 2.22597e-05
+SHARC_D4_STRIP_BACK6_E 0.12692 2.28943e-05
+SHARC_D4_STRIP_BACK7_E 0.090648 2.32e-05
+SHARC_D4_STRIP_BACK8_E 0.110853 2.26683e-05
+SHARC_D4_STRIP_BACK9_E 0.106239 2.33143e-05
+SHARC_D4_STRIP_BACK10_E 0.12392 2.30434e-05
+SHARC_D4_STRIP_BACK11_E 0.116171 2.19613e-05
+SHARC_D4_STRIP_BACK12_E 0.115817 2.19587e-05
+SHARC_D4_STRIP_BACK13_E 0.116895 2.20284e-05
+SHARC_D4_STRIP_BACK14_E 0.11774 2.21507e-05
+SHARC_D4_STRIP_BACK15_E 0.108548 2.24632e-05
+SHARC_D4_STRIP_BACK16_E 0.126504 2.24063e-05
+SHARC_D4_STRIP_BACK18_E 0.13103 2.26364e-05
+SHARC_D4_STRIP_BACK21_E 0.116827 2.23377e-05
+SHARC_D4_STRIP_BACK22_E 0.120895 2.24211e-05
+SHARC_D4_STRIP_BACK23_E 0.103585 2.26749e-05
+SHARC_D4_STRIP_BACK24_E 0.128199 2.23913e-05
+SHARC_D5_STRIP_FRONT1_E 0.0184746 2.19244e-05
+SHARC_D5_STRIP_FRONT2_E 0.0220658 2.18118e-05
+SHARC_D5_STRIP_FRONT3_E 0.0258988 2.20747e-05
+SHARC_D5_STRIP_FRONT4_E 0.012098 2.21849e-05
+SHARC_D5_STRIP_FRONT5_E 0.0132119 2.187e-05
+SHARC_D5_STRIP_FRONT6_E 0.0231413 2.2357e-05
+SHARC_D5_STRIP_FRONT8_E 0.0971659 2.16172e-05
+SHARC_D5_STRIP_FRONT9_E 0.04333 2.12817e-05
+SHARC_D5_STRIP_FRONT10_E 0.0507796 2.15983e-05
+SHARC_D5_STRIP_FRONT11_E 0.0243644 2.19276e-05
+SHARC_D5_STRIP_FRONT12_E 0.0259251 2.17859e-05
+SHARC_D5_STRIP_FRONT13_E 0.0335124 2.19481e-05
+SHARC_D5_STRIP_FRONT14_E 0.0399368 2.17655e-05
+SHARC_D5_STRIP_FRONT15_E 0.0184578 2.12683e-05
+SHARC_D5_STRIP_FRONT16_E 0.022659 2.21454e-05
+SHARC_D5_STRIP_BACK11_E 0.0568813 2.07567e-05
+SHARC_D5_STRIP_BACK12_E 0.0585774 2.08176e-05
+SHARC_D5_STRIP_BACK13_E 0.0653997 2.08075e-05
+SHARC_D5_STRIP_BACK14_E 0.0678985 2.10715e-05
+SHARC_D5_STRIP_BACK15_E 0.0595765 2.09586e-05
+SHARC_D5_STRIP_BACK16_E 0.0626579 2.09575e-05
+SHARC_D5_STRIP_BACK17_E 0.0598144 2.08644e-05
+SHARC_D5_STRIP_BACK18_E 0.0615517 2.06992e-05
+SHARC_D5_STRIP_BACK19_E 0.0617771 2.1087e-05
+SHARC_D5_STRIP_BACK20_E 0.0616774 2.11152e-05
+SHARC_D5_STRIP_BACK21_E 0.0536845 2.0387e-05
+SHARC_D5_STRIP_BACK22_E 0.050151 2.09927e-05
+SHARC_D5_STRIP_BACK23_E 0.04769 2.10759e-05
+SHARC_D5_STRIP_BACK24_E -0.0798476 2.1628e-05
+SHARC_D5_STRIP_BACK25_E 0.0446056 2.15827e-05
+SHARC_D5_STRIP_BACK27_E 0.0474283 2.1244e-05
+SHARC_D5_STRIP_BACK28_E 0.0517215 2.14623e-05
+SHARC_D5_STRIP_BACK29_E 0.042902 2.11875e-05
+SHARC_D5_STRIP_BACK30_E 0.0439854 2.14907e-05
+SHARC_D5_STRIP_BACK31_E 0.0498488 2.11159e-05
+SHARC_D5_STRIP_BACK32_E 0.0418072 2.15367e-05
+SHARC_D5_STRIP_BACK33_E 0.0401053 2.12155e-05
+SHARC_D5_STRIP_BACK34_E 0.0405886 2.1192e-05
+SHARC_D5_STRIP_BACK38_E 0.0436775 2.17032e-05
+SHARC_D5_STRIP_BACK39_E 0.0434241 2.12179e-05
+SHARC_D5_STRIP_BACK40_E 0.0491315 2.16629e-05
+SHARC_D5_STRIP_BACK43_E 0.058932 2.0904e-05
+SHARC_D5_STRIP_BACK44_E 0.0541304 2.0784e-05
+SHARC_D6_STRIP_FRONT3_E 0.0267244 2.14665e-05
+SHARC_D6_STRIP_FRONT4_E 0.0194191 2.16453e-05
+SHARC_D6_STRIP_FRONT5_E 0.0214514 2.1681e-05
+SHARC_D6_STRIP_FRONT6_E 0.0185885 2.14626e-05
+SHARC_D6_STRIP_FRONT7_E 0.0107993 2.16074e-05
+SHARC_D6_STRIP_FRONT8_E 0.0220856 2.17046e-05
+SHARC_D6_STRIP_FRONT9_E 3.53799 0.000324051
+SHARC_D6_STRIP_FRONT10_E 0.0147158 2.17499e-05
+SHARC_D6_STRIP_FRONT11_E 0.0180041 2.11129e-05
+SHARC_D6_STRIP_FRONT12_E 0.00928292 2.1948e-05
+SHARC_D6_STRIP_FRONT13_E 0.0130352 2.14338e-05
+SHARC_D6_STRIP_FRONT14_E 0.0126762 2.13344e-05
+SHARC_D6_STRIP_FRONT15_E 0.0130653 2.18171e-05
+SHARC_D6_STRIP_FRONT16_E 0.011682 2.17847e-05
+SHARC_D6_STRIP_FRONT17_E 0.0155902 2.20876e-05
+SHARC_D6_STRIP_FRONT18_E 0.0108346 2.18593e-05
+SHARC_D6_STRIP_FRONT19_E 0.0145107 2.20757e-05
+SHARC_D6_STRIP_FRONT20_E 0.0107507 2.2692e-05
+SHARC_D6_STRIP_FRONT21_E 0.0140636 2.25682e-05
+SHARC_D6_STRIP_FRONT22_E 0.0172886 2.2133e-05
+SHARC_D6_STRIP_FRONT23_E -0.0252202 2.20485e-05
+SHARC_D6_STRIP_BACK1_E -0.169393 2.38599e-05
+SHARC_D6_STRIP_BACK2_E 0.0275962 2.29444e-05
+SHARC_D6_STRIP_BACK3_E 0.075551 2.30798e-05
+SHARC_D6_STRIP_BACK4_E 0.0514595 2.28556e-05
+SHARC_D6_STRIP_BACK5_E 0.0581235 2.25471e-05
+SHARC_D6_STRIP_BACK6_E 0.0435178 2.26392e-05
+SHARC_D6_STRIP_BACK7_E 0.0463227 2.24721e-05
+SHARC_D6_STRIP_BACK8_E 0.0464039 2.26593e-05
+SHARC_D6_STRIP_BACK9_E 0.0464575 2.23355e-05
+SHARC_D6_STRIP_BACK10_E 0.0440967 2.27238e-05
+SHARC_D6_STRIP_BACK11_E 0.0395257 2.25909e-05
+SHARC_D6_STRIP_BACK12_E 0.0222998 2.22406e-05
+SHARC_D6_STRIP_BACK13_E 0.0178924 2.23953e-05
+SHARC_D6_STRIP_BACK14_E 0.0338867 2.26127e-05
+SHARC_D6_STRIP_BACK15_E 0.0373955 2.33028e-05
+SHARC_D6_STRIP_BACK16_E 0.0276929 2.31612e-05
+SHARC_D6_STRIP_BACK17_E 0.0162478 2.29907e-05
+SHARC_D6_STRIP_BACK18_E 0.0220379 2.29427e-05
+SHARC_D6_STRIP_BACK19_E 0.022612 2.36231e-05
+SHARC_D6_STRIP_BACK20_E 0.0186129 2.32864e-05
+SHARC_D6_STRIP_BACK21_E 0.0151438 2.33919e-05
+SHARC_D6_STRIP_BACK23_E 0.0158743 2.27083e-05
+SHARC_D6_STRIP_BACK25_E 0.0381616 2.27827e-05
+SHARC_D6_STRIP_BACK26_E 0.0248942 2.29145e-05
+SHARC_D6_STRIP_BACK27_E 0.0356308 2.26366e-05
+SHARC_D6_STRIP_BACK29_E 0.0233163 2.29887e-05
+SHARC_D6_STRIP_BACK30_E 0.041849 2.30691e-05
+SHARC_D6_STRIP_BACK31_E 0.0369897 2.27207e-05
+SHARC_D6_STRIP_BACK33_E 0.0441671 2.30209e-05
+SHARC_D6_STRIP_BACK34_E 0.0360893 2.23514e-05
+SHARC_D6_STRIP_BACK35_E 0.0339634 2.24387e-05
+SHARC_D6_STRIP_BACK36_E 0.060896 2.26085e-05
+SHARC_D6_STRIP_BACK37_E 0.00208911 3.18512e-05
+SHARC_D6_STRIP_BACK38_E 0.0168685 2.51637e-05
+SHARC_D6_STRIP_BACK39_E 0.0434731 2.29497e-05
+SHARC_D6_STRIP_BACK40_E 0.0385564 2.24943e-05
+SHARC_D6_STRIP_BACK41_E 3.53799 0.000324051
+SHARC_D6_STRIP_BACK42_E 0.0385198 2.24836e-05
+SHARC_D6_STRIP_BACK43_E 0.0285279 2.29171e-05
+SHARC_D6_STRIP_BACK44_E 0.0265627 2.25233e-05
+SHARC_D6_STRIP_BACK45_E 0.0333433 2.2531e-05
+SHARC_D6_STRIP_BACK46_E 0.0200449 2.29807e-05
+SHARC_D6_STRIP_BACK47_E 0.0293445 2.24907e-05
+SHARC_D7_STRIP_FRONT1_E 0.0111471 2.14709e-05
+SHARC_D7_STRIP_FRONT2_E 0.0109067 2.19744e-05
+SHARC_D7_STRIP_FRONT3_E -0.0144983 2.12644e-05
+SHARC_D7_STRIP_FRONT4_E 0.00958447 2.21958e-05
+SHARC_D7_STRIP_FRONT5_E 0.00809943 2.21e-05
+SHARC_D7_STRIP_FRONT6_E 0.0110324 2.22059e-05
+SHARC_D7_STRIP_FRONT7_E 0.016297 2.17193e-05
+SHARC_D7_STRIP_FRONT8_E 0.0133115 2.20289e-05
+SHARC_D7_STRIP_FRONT9_E 0.00415246 2.08826e-05
+SHARC_D7_STRIP_FRONT10_E 0.00902517 2.13638e-05
+SHARC_D7_STRIP_FRONT11_E 0.00770633 2.15624e-05
+SHARC_D7_STRIP_FRONT12_E -0.00164026 2.13971e-05
+SHARC_D7_STRIP_FRONT13_E 0.00403168 2.16003e-05
+SHARC_D7_STRIP_FRONT14_E 0.00881677 2.22311e-05
+SHARC_D7_STRIP_FRONT15_E -0.0334207 2.20094e-05
+SHARC_D7_STRIP_FRONT16_E 0.00373577 2.20422e-05
+SHARC_D7_STRIP_FRONT17_E 0.00897941 2.21518e-05
+SHARC_D7_STRIP_FRONT18_E -0.000181608 2.2362e-05
+SHARC_D7_STRIP_FRONT19_E 0.00775456 2.20313e-05
+SHARC_D7_STRIP_FRONT20_E 0.0155661 2.19076e-05
+SHARC_D7_STRIP_FRONT21_E 0.0168267 2.21411e-05
+SHARC_D7_STRIP_FRONT22_E 0.00881832 2.23677e-05
+SHARC_D7_STRIP_FRONT23_E 0.00365454 2.21892e-05
+SHARC_D7_STRIP_BACK12_E 0.0490787 2.28767e-05
+SHARC_D7_STRIP_BACK13_E 0.0543162 2.27022e-05
+SHARC_D7_STRIP_BACK14_E 0.0422382 2.26245e-05
+SHARC_D7_STRIP_BACK15_E 0.0255048 2.30319e-05
+SHARC_D7_STRIP_BACK16_E -0.0418349 2.28575e-05
+SHARC_D7_STRIP_BACK17_E 0.0449648 2.30128e-05
+SHARC_D7_STRIP_BACK18_E 0.0311507 2.31352e-05
+SHARC_D7_STRIP_BACK19_E 0.0424563 2.29341e-05
+SHARC_D7_STRIP_BACK20_E 0.0417637 2.28612e-05
+SHARC_D7_STRIP_BACK21_E 0.0373395 2.289e-05
+SHARC_D7_STRIP_BACK22_E 0.0326856 2.28419e-05
+SHARC_D7_STRIP_BACK23_E 0.0357087 2.25265e-05
+SHARC_D7_STRIP_BACK24_E 0.0280085 2.29681e-05
+SHARC_D7_STRIP_BACK25_E 0.0342112 2.34278e-05
+SHARC_D7_STRIP_BACK26_E 0.0367721 2.29647e-05
+SHARC_D7_STRIP_BACK27_E 0.0388958 2.35108e-05
+SHARC_D7_STRIP_BACK29_E 0.0297903 2.32002e-05
+SHARC_D7_STRIP_BACK31_E 0.0247006 2.24339e-05
+SHARC_D7_STRIP_BACK32_E 0.0378452 2.2825e-05
+SHARC_D7_STRIP_BACK33_E 0.0397023 2.32915e-05
+SHARC_D7_STRIP_BACK34_E 0.0288076 2.30423e-05
+SHARC_D7_STRIP_BACK36_E 0.0404461 2.2752e-05
+SHARC_D7_STRIP_BACK38_E 0.0233955 2.34949e-05
+SHARC_D7_STRIP_BACK39_E 0.0298339 2.2881e-05
+SHARC_D7_STRIP_BACK40_E 0.0280095 2.32431e-05
+SHARC_D7_STRIP_BACK41_E 0.0241983 2.29269e-05
+SHARC_D7_STRIP_BACK42_E 0.0218632 2.29931e-05
+SHARC_D7_STRIP_BACK43_E 0.0206106 2.30039e-05
+SHARC_D7_STRIP_BACK44_E 0.0136246 2.25895e-05
+SHARC_D7_STRIP_BACK45_E 0.0337216 2.34113e-05
+SHARC_D7_STRIP_BACK46_E 0.0251993 2.32068e-05
+SHARC_D7_STRIP_BACK47_E 0.0306201 2.27172e-05
+SHARC_D8_STRIP_FRONT1_E -0.00938955 2.13663e-05
+SHARC_D8_STRIP_FRONT2_E -0.000148791 2.14337e-05
+SHARC_D8_STRIP_FRONT3_E -0.00615021 2.16178e-05
+SHARC_D8_STRIP_FRONT4_E 0.0023255 2.13416e-05
+SHARC_D8_STRIP_FRONT5_E 0.00514757 2.17759e-05
+SHARC_D8_STRIP_FRONT11_E -0.00672231 2.14599e-05
+SHARC_D8_STRIP_FRONT12_E -0.00684039 2.14697e-05
+SHARC_D8_STRIP_FRONT13_E -0.00810409 2.207e-05
+SHARC_D8_STRIP_FRONT16_E -0.00685472 2.20266e-05
+SHARC_D8_STRIP_FRONT17_E -0.0131768 2.2297e-05
+SHARC_D8_STRIP_FRONT18_E -0.0029741 2.19304e-05
+SHARC_D8_STRIP_FRONT19_E -0.00693736 2.22931e-05
+SHARC_D8_STRIP_FRONT20_E -0.00171549 2.22701e-05
+SHARC_D8_STRIP_FRONT21_E -0.00929376 2.21767e-05
+SHARC_D8_STRIP_FRONT22_E -0.0278538 2.21839e-05
+SHARC_D8_STRIP_FRONT23_E -0.0022602 2.20649e-05
+SHARC_D8_STRIP_BACK2_E 3.53799 0.000324051
+SHARC_D8_STRIP_BACK3_E 0.0980749 2.72316e-05
+SHARC_D8_STRIP_BACK4_E 0.0567821 2.32254e-05
+SHARC_D8_STRIP_BACK5_E 0.0721435 2.81065e-05
+SHARC_D8_STRIP_BACK6_E 0.0426144 2.28552e-05
+SHARC_D8_STRIP_BACK7_E 0.0413196 2.27331e-05
+SHARC_D8_STRIP_BACK8_E 0.0409884 2.2848e-05
+SHARC_D8_STRIP_BACK9_E 0.0418038 2.32477e-05
+SHARC_D8_STRIP_BACK10_E 0.0316771 2.31747e-05
+SHARC_D8_STRIP_BACK11_E 0.0372514 2.26756e-05
+SHARC_D8_STRIP_BACK12_E 0.0370299 2.27828e-05
+SHARC_D8_STRIP_BACK13_E 0.0301217 2.25029e-05
+SHARC_D8_STRIP_BACK14_E 0.0286195 2.25003e-05
+SHARC_D8_STRIP_BACK15_E 0.0299568 2.29574e-05
+SHARC_D8_STRIP_BACK17_E 3.53799 0.000324051
+SHARC_D8_STRIP_BACK18_E 0.000891401 2.39545e-05
+SHARC_D8_STRIP_BACK19_E 0.00753597 2.40669e-05
+SHARC_D8_STRIP_BACK20_E 0.00554838 2.34512e-05
+SHARC_D8_STRIP_BACK21_E 0.0156689 2.32156e-05
+SHARC_D8_STRIP_BACK22_E 0.0160044 2.27685e-05
+SHARC_D8_STRIP_BACK23_E 0.00154922 2.35261e-05
+SHARC_D8_STRIP_BACK24_E 0.00941722 2.3384e-05
+SHARC_D8_STRIP_BACK25_E 0.00684016 2.32117e-05
+SHARC_D8_STRIP_BACK26_E -0.00548732 2.34809e-05
+SHARC_D8_STRIP_BACK27_E 0.0201157 2.34148e-05
+SHARC_D8_STRIP_BACK28_E 0.018198 2.36256e-05
+SHARC_D8_STRIP_BACK29_E 0.0190719 2.33848e-05
+SHARC_D8_STRIP_BACK30_E 0.0292787 2.29745e-05
+SHARC_D8_STRIP_BACK31_E 0.0285674 2.25822e-05
+SHARC_D8_STRIP_BACK32_E 0.0231226 2.28624e-05
+SHARC_D8_STRIP_BACK33_E 0.0386775 2.27137e-05
+SHARC_D8_STRIP_BACK34_E 0.0290607 2.29847e-05
+SHARC_D8_STRIP_BACK35_E 0.0341422 2.28036e-05
+SHARC_D8_STRIP_BACK36_E 3.53799 0.000324051
+SHARC_D8_STRIP_BACK37_E 0.0188344 2.32843e-05
+SHARC_D8_STRIP_BACK38_E 0.0349254 2.26834e-05
+SHARC_D8_STRIP_BACK39_E 0.0352429 2.33297e-05
+SHARC_D8_STRIP_BACK40_E 0.0323803 2.15357e-05
+SHARC_D8_STRIP_BACK41_E 0.0250631 2.35277e-05
+SHARC_D8_STRIP_BACK43_E 0.0255274 2.2294e-05
+SHARC_D8_STRIP_BACK44_E 0.0234287 2.32661e-05
+SHARC_D8_STRIP_BACK45_E 0.0270765 2.30607e-05
+SHARC_D8_STRIP_BACK46_E 0.0161063 2.26809e-05
+SHARC_D8_STRIP_BACK47_E 0.0176102 2.31787e-05
+SHARC_D10_STRIP_FRONT1_E -0.0246137 2.40276e-05
+SHARC_D10_STRIP_FRONT2_E -0.0247509 2.36702e-05
+SHARC_D10_STRIP_FRONT3_E -0.0223287 2.35191e-05
+SHARC_D10_STRIP_FRONT4_E -0.0231089 2.41774e-05
+SHARC_D10_STRIP_FRONT5_E -0.02485 2.39843e-05
+SHARC_D10_STRIP_FRONT6_E -0.0287019 2.39263e-05
+SHARC_D10_STRIP_FRONT7_E -0.0188632 2.28559e-05
+SHARC_D10_STRIP_FRONT8_E -0.0156374 2.30957e-05
+SHARC_D10_STRIP_FRONT9_E -0.0159356 2.30804e-05
+SHARC_D10_STRIP_FRONT10_E -0.0227716 2.33393e-05
+SHARC_D10_STRIP_FRONT11_E -0.0162527 2.31961e-05
+SHARC_D10_STRIP_FRONT12_E -0.0213478 2.37989e-05
+SHARC_D10_STRIP_FRONT13_E -0.0136014 2.36708e-05
+SHARC_D10_STRIP_FRONT14_E -0.0126604 2.43061e-05
+SHARC_D10_STRIP_FRONT15_E -0.00918633 2.43673e-05
+SHARC_D10_STRIP_FRONT16_E -0.023791 2.35819e-05
+SHARC_D10_STRIP_FRONT17_E -0.0294735 2.42918e-05
+SHARC_D10_STRIP_FRONT18_E -0.0162182 2.39197e-05
+SHARC_D10_STRIP_FRONT19_E -0.0223356 2.45194e-05
+SHARC_D10_STRIP_FRONT20_E -0.0177009 2.39589e-05
+SHARC_D10_STRIP_FRONT21_E -0.0227187 2.43008e-05
+SHARC_D10_STRIP_FRONT22_E -0.0134246 2.37881e-05
+SHARC_D10_STRIP_FRONT23_E -0.0226326 2.38859e-05
+SHARC_D10_STRIP_BACK11_E 0.0584119 2.55133e-05
+SHARC_D10_STRIP_BACK12_E 0.0699967 2.23195e-05
+SHARC_D10_STRIP_BACK13_E 0.0682196 2.45787e-05
+SHARC_D10_STRIP_BACK14_E 0.0503693 2.26041e-05
+SHARC_D10_STRIP_BACK15_E 0.0502025 2.27367e-05
+SHARC_D10_STRIP_BACK16_E 0.0586348 2.2878e-05
+SHARC_D10_STRIP_BACK17_E 0.0603656 2.28189e-05
+SHARC_D10_STRIP_BACK18_E 0.0591394 2.27393e-05
+SHARC_D10_STRIP_BACK19_E 0.0549039 2.2804e-05
+SHARC_D10_STRIP_BACK20_E 0.0494821 2.32013e-05
+SHARC_D10_STRIP_BACK21_E 0.0487158 2.20511e-05
+SHARC_D10_STRIP_BACK22_E 0.0405202 2.2122e-05
+SHARC_D10_STRIP_BACK23_E 0.0319081 2.21269e-05
+SHARC_D10_STRIP_BACK24_E 0.0235831 2.25546e-05
+SHARC_D10_STRIP_BACK25_E 0.0264097 2.24409e-05
+SHARC_D10_STRIP_BACK26_E 0.0318306 2.27382e-05
+SHARC_D10_STRIP_BACK27_E 0.0362605 2.25108e-05
+SHARC_D10_STRIP_BACK28_E 0.0275711 2.30077e-05
+SHARC_D10_STRIP_BACK29_E 0.0337957 2.28222e-05
+SHARC_D10_STRIP_BACK30_E 0.0270704 2.33736e-05
+SHARC_D10_STRIP_BACK31_E 0.0360256 2.22405e-05
+SHARC_D10_STRIP_BACK32_E 0.0481137 2.27023e-05
+SHARC_D10_STRIP_BACK33_E 0.0441632 2.26257e-05
+SHARC_D10_STRIP_BACK34_E 0.0344984 2.28046e-05
+SHARC_D10_STRIP_BACK35_E 0.0407187 2.25085e-05
+SHARC_D10_STRIP_BACK36_E 0.0446098 2.28614e-05
+SHARC_D10_STRIP_BACK37_E 0.0369654 2.24067e-05
+SHARC_D10_STRIP_BACK38_E 0.0426523 2.32392e-05
+SHARC_D10_STRIP_BACK39_E 0.042786 2.29164e-05
+SHARC_D10_STRIP_BACK40_E 0.057569 2.23818e-05
+SHARC_D10_STRIP_BACK41_E 0.0505517 2.22661e-05
+SHARC_D10_STRIP_BACK42_E 0.0427701 2.31198e-05
+SHARC_D10_STRIP_BACK43_E 0.0325823 2.28741e-05
+SHARC_D10_STRIP_BACK44_E 0.0382447 2.2782e-05
+SHARC_D10_STRIP_BACK45_E 0.041409 2.28557e-05
+SHARC_D10_STRIP_BACK46_E 0.0428671 2.32039e-05
+SHARC_D10_STRIP_BACK47_E 0.0471655 2.3262e-05
+SHARC_D12_STRIP_FRONT1_E -0.00484915 2.33405e-05
+SHARC_D12_STRIP_FRONT2_E -0.0014674 2.34682e-05
+SHARC_D12_STRIP_FRONT3_E -0.0130518 2.33915e-05
+SHARC_D12_STRIP_FRONT5_E -0.0190662 2.37877e-05
+SHARC_D12_STRIP_FRONT6_E -0.0207672 2.43109e-05
+SHARC_D12_STRIP_FRONT8_E -0.00795163 2.36519e-05
+SHARC_D12_STRIP_FRONT9_E -0.0137762 2.36385e-05
+SHARC_D12_STRIP_FRONT10_E -0.0179864 2.36641e-05
+SHARC_D12_STRIP_FRONT11_E -0.0171115 2.34168e-05
+SHARC_D12_STRIP_FRONT12_E -0.0156885 2.28764e-05
+SHARC_D12_STRIP_FRONT13_E -0.0107117 2.35969e-05
+SHARC_D12_STRIP_FRONT14_E -0.0229854 2.37416e-05
+SHARC_D12_STRIP_FRONT15_E -0.0137189 2.38798e-05
+SHARC_D12_STRIP_FRONT16_E -0.0140456 2.40203e-05
+SHARC_D12_STRIP_FRONT17_E -0.0200381 2.36551e-05
+SHARC_D12_STRIP_FRONT18_E -0.00805789 2.35912e-05
+SHARC_D12_STRIP_FRONT19_E -0.0186952 2.43875e-05
+SHARC_D12_STRIP_FRONT20_E -0.0144919 2.45319e-05
+SHARC_D12_STRIP_FRONT21_E -0.0163102 2.33092e-05
+SHARC_D12_STRIP_FRONT22_E -0.0137205 2.36079e-05
+SHARC_D12_STRIP_FRONT23_E -0.0228041 2.40822e-05
+SHARC_D12_STRIP_BACK1_E 3.53799 0.000324051
+SHARC_D12_STRIP_BACK2_E -0.168899 2.64359e-05
+SHARC_D12_STRIP_BACK3_E 0.0255178 2.53173e-05
+SHARC_D12_STRIP_BACK4_E 0.00579596 2.50901e-05
+SHARC_D12_STRIP_BACK5_E 0.0135761 2.5215e-05
+SHARC_D12_STRIP_BACK6_E 0.00734285 2.55434e-05
+SHARC_D12_STRIP_BACK7_E 0.0193437 2.46671e-05
+SHARC_D12_STRIP_BACK8_E 0.00600789 2.50851e-05
+SHARC_D12_STRIP_BACK9_E 0.00144002 2.51709e-05
+SHARC_D12_STRIP_BACK10_E -0.0122236 2.49197e-05
+SHARC_D12_STRIP_BACK11_E -0.00690914 2.49627e-05
+SHARC_D12_STRIP_BACK12_E -0.0190408 2.50708e-05
+SHARC_D12_STRIP_BACK13_E -0.0136739 2.52216e-05
+SHARC_D12_STRIP_BACK14_E -0.000614762 2.48763e-05
+SHARC_D12_STRIP_BACK15_E -0.0115212 2.50992e-05
+SHARC_D12_STRIP_BACK16_E -0.0113651 2.48327e-05
+SHARC_D12_STRIP_BACK17_E -0.00670837 2.55571e-05
+SHARC_D12_STRIP_BACK18_E 0.013741 2.57646e-05
+SHARC_D12_STRIP_BACK19_E -0.000419331 2.51958e-05
+SHARC_D12_STRIP_BACK20_E -0.00456227 2.91852e-05
+SHARC_D12_STRIP_BACK21_E -0.0226305 2.50211e-05
+SHARC_D12_STRIP_BACK22_E -0.0111732 2.56048e-05
+SHARC_D12_STRIP_BACK23_E -0.00691075 2.51414e-05
+SHARC_D12_STRIP_BACK24_E -0.00912395 2.51848e-05
+SHARC_D12_STRIP_BACK25_E -0.0104703 2.52789e-05
+SHARC_D12_STRIP_BACK26_E 0.00903197 2.41685e-05
+SHARC_D12_STRIP_BACK27_E -0.00597004 2.5631e-05
+SHARC_D12_STRIP_BACK28_E -0.00411177 2.51907e-05
+SHARC_D12_STRIP_BACK29_E -0.00238402 2.55322e-05
+SHARC_D12_STRIP_BACK30_E -0.0019266 2.56927e-05
+SHARC_D12_STRIP_BACK32_E 0.00903099 2.46706e-05
+SHARC_D12_STRIP_BACK33_E -0.000226763 2.46582e-05
+SHARC_D12_STRIP_BACK34_E 0.0165983 2.48672e-05
+SHARC_D12_STRIP_BACK35_E 0.00202338 2.48451e-05
+SHARC_D12_STRIP_BACK36_E 0.0028113 2.47164e-05
+SHARC_D12_STRIP_BACK37_E 0.00321136 2.50483e-05
+SHARC_D12_STRIP_BACK38_E 0.00273738 2.49378e-05
+SHARC_D12_STRIP_BACK39_E 0.00942349 2.46461e-05
+SHARC_D12_STRIP_BACK40_E -0.00127344 2.51202e-05
+SHARC_D12_STRIP_BACK41_E -0.0502734 2.51267e-05
+SHARC_D12_STRIP_BACK42_E -0.00239558 2.45522e-05
+SHARC_D12_STRIP_BACK43_E -0.00143522 2.4463e-05
+SHARC_D12_STRIP_BACK44_E -0.00154555 2.54788e-05
+SHARC_D12_STRIP_BACK45_E -0.00260741 2.50383e-05
+SHARC_D12_STRIP_BACK46_E -0.00107954 2.52285e-05
+SHARC_D12_STRIP_BACK47_E 0.0096826 2.5544e-05
diff --git a/NPAnalysis/S1554/SharcEff.root b/NPAnalysis/S1554/SharcEff.root
new file mode 100644
index 0000000000000000000000000000000000000000..7a23fa69e60c04c900074535a14ae2c71c13133d
GIT binary patch
literal 5788
zcmbW5bxfSWzwQ?(QrvZsMT-}A_k~i56?d1#rMNp3DN?++ZgH35UY6prP~2Txq=)wW
zadU5e$vrufd^0nddFGiv-Z%No<Lcz(4gegU0RR9C0Dw6Axy+x(LeGWpT;ebOdw@;@
z0AS1l5a6n{p3M1p+vy|?T6zgxr_ZbZqviwndyw+GrIj`S=f9Jm-2eb&O?g`fODqi=
zQ&)32D=QAyzqk2+rU2x>s((rV0P-vnf3AOO0D$QF|97RA(tn-NJO1fK{yDGy)xXeG
z<dVm-k+pQUG}pATbT^fG^MB9uRM4K#03GN6Z-nAhR82x#{GlupGpTiYXg>W%1Uw`#
zX<Lk`i{S<+Bx6ug{S%EYAzmD73@r=T2n81vk3=>CRJ#9eMoQAu1dT;X?w0ThAotem
z=W))i$Ze<ZX@&9o$L^<t*{_xD#(bfP>-fo-X7VCjp;EZK($aoWZ(<U&nRn|)MP`<#
zXNN>);7jL+se)dNAm(kM&zy79OpP;%IwMi3R9|e!WK_r_rEu?~b<Mus7jsLg#5#Xm
zAS%cw5Z3k4FKRUZ#5ghLTl%D#%PXbmmpRH`^Rsv-Btb0#Sp6}t3eEKTe1j4xga#5b
z$!vZU^-Jc_%LG^Uj?nOu^>>$2TXuRvO82d**Bd>>nywN-jr2Pn*A=mzz{JwrN%}EY
z?;nB7q%UgiV?=LRFP@aRL8}jLI&m!%LMo3EOAmP_8}rvIV8@ygn(w-KUwjAWzkzdV
z3=XOcSdZC<q7F5{H!t^0;&4QL{7&1c?Flk8ydJ&^Hr{Or8gI>`$K&k3OCNLNlGk(2
zEk22@<;#N*8YuQnHxyas(Q*82MTj(9VKC1iU78+IVgIH)-dg`=u!5_6u*huqPP8Ok
zdVH{E@=p1p67$j*VYtq|flL~K%9xx~j8}nzVF#znQk+fFTx?Bz@g+IcWIko$-fE~)
zX=g4}+0rbn=6m~O0p{WuCx74`tmm)|YT6ODsngOK$^goU4u>*MCY<c`75LLsaTo>3
zg{D_dJ7t;F&W-+}V>s4fQ@gh~*$fw;3+(hhtHxe^=<qDQ;zn1_6dCe32-CTEB4bTy
zL=nPK7bR>vs^Qgo$D?_N-im3NRyu07b)&jgvm}N%)@`p3b!|r*^U!vp!m1sNvDJFn
z4O^<EHEb!;R0~f<X`L#S?yKdtFPOdP6|IV~_)$1GR!hDb#<L|>b5QlTkIK4g)@b<q
z2D@)UhhR4AI*}rZ?^f$=%K>$d9+^XMkYAP98U;*(G`taX(5^64P}9q5z9E2xH|W#S
z@AiNJyKpTT_<?f?)3scaRGM}(7|?rP<ETM8dY=Twvtl`GI6g>RFYzJWPNlyeos=UZ
znJx)6f8oPgJh9Gw&$Xz^vQz4hG@L|f``V}lfx`9vHk}^~a#|XAQW;@bfVsF%vCfkW
z>-WS46heYJaKFTEoO~E&oqmc9uqL|5-QKEv3zV^o?WzmCB`Dfk+zlDsQ)Fjit@bNS
z66iUpas<nbS}B(pT9Bz{30fZ!w}P%si)C`^vD6jbf83t9lCt9+Snc;jx#-4-PQTKY
z7UxOE6rD=YJdK*kxjs$una+VqDNEX&skebwl($&5#lAfBCi;k{&G@lU{PO8s@)vMV
z<fO(cteT5|h(~|xq)T%FCsCnqMp-{j3>P1r9Jmfu2ta?}LMz{CvkVu`N(vmVyPhcY
z+p3&*S%?@J_VT-I3I|F*%~frKrymC%js`$(6Nz~RS}9JIj0v&6+0@u{9;csdr8H>O
zJ>qlgu(cICIYb|x2x&dK<x~azNXjEh<Rb_DirRKddFuxrA-to%t+S}r$>7U7S!|N4
zS@-sr7nq6M+(Ji!H1UMhN`<@Az%6)r22cAW*@y91zDJ6on)T7c{Zm>Sowc#9>pc~Q
zxtAGCsHs{t{G5m2+bjAis1gNE1Fnf>;_%H9f*f03^SmQ11mxProM!zQgV-&i;k#g&
zDc+i;6<C$rNLl+)zhGahqfYU&n`ahYiFFo#O$(BDTz7-HZY($S*Ov_suNp$AX*(m+
zpI6arp^X<$hwe2YfrgXqF2G8s;qlkYI659<iaXy5B&+I_KKx*w3M10xKG1jYGTA1N
z?C$bAj3H9?IWjjRHv3@a0KNj_j~(Y(99>}cTy(G^^quGElz10}MLKLf%_WX}UJqZX
zst9;BAqhP_`}MfOy9J>XV=GS$%FbrAt*+JEVms!xo1p4XcB$;S{AK9PBX+4cqlf5Z
z1XcHjIz|a>8#MY-F6d+nv|XUE;yQJREYmGtqQ{2i5a>--)Vwl@ZqiIgAdV0@esi%H
z6b|v#<V8U}f3F|>pkJ-RcIRt-kkr2I?cZg2noKRe(b#!!qs7CqKAVs50i(yv;x=4%
z|H6xI83*p6Pk46GvUyxnDRNeq_F}^2QF|P!bLZ2n2dq;MS6Oc-7Y)>D^;}>kc(mHa
zqPq5pZ*%Olqh4I=o9+Ns;#@w7O?Kx({|v=Q;iC8Ro}WrPsIOIA&MFFUTB7O@W;)R_
z%y=fVs81@^hqi4knMwSRxIy|rKa?q&4ppuznosJMk4Va_E&AMaq~HdX;wq1XTU2!O
zF0KrNBx>DEppE$Bfe=)Tbn;-(BNr)&-wTWqE}fhv<gA&{pEC|Pkopi4w=x#*mc*X~
zkNAjKIW$Lf)U=)AGAdghf(QLUdGOL6rH{wXdS>kXU%QcaPwXW(+I$rCTrkzIy=3&h
z0xfoyWwbpIwvK?oFIstBA1OhGT4E8o$NLJ+SmKbhcip}2NX*>+`FIYVWzd1%oe_e%
zfKnVFw+`}k%S&&*W6A6KpW04>ymiIvs!Fsg+d{$LCP+>y$$K8#Ghyn57}c7rqOUsu
zI6UuV^z(rYSVL0nk0`j}(6EG4f_#l~pztAeuuU;>l$mvKlrO0)oh^(3gkeaL_2qC>
zA|In#a%KQw3WT*MY=&46?TjmAj5yUGN0Uj1oD<QUzo{vT)_#%S`5m~#g0Msa434CN
zoY8mz^=#u|z$fxdxA#&hHGmfZrR%G-!fM5*U*^j3!D7aeQXwc?9cPiBa+cSR#QA|^
zaExG5)GUHYz+)dB%%)ZKrk^C~D%;&B^J&sOt&h)@_YLuPk9bLNE|x)oQ&>N1J&kED
z8h(f!=;V$iit+;vQ-SMGQeO_>2}Ka;BB9{7S0S)TX56X@)PMx0p6PkQrn=@J{r%T1
zcT|W9!ORQZvO!w@uVKjUI;KD@`7^Yun*fxF_euGlGA}=A&YIq67-y`~qxP>NZACRk
z9i6(<E^AQ^57x;h6PSWeYK_d-@!nfMrGymh`|LzG`S-82-H7<aG6(;|zGs4>v-?L-
zl79*6FAM@KkTq3o-P|!X++8hA9W7mzAXZNJvQ8dm4j*V7Y$5iR7SElV`!hzxVF2vt
zo-t}xZ_Nq#o=BX^_ju6WAzRQkOWk#UvR{dOnO*LUjff|(GMJa+y2K;E+S6+=*UEG7
z^XK^ECRplk$>J_GNogW~geQc4j8w_87q4(=5I7Y-1fx>@4w3A*dRgwbI}Se|_jEjM
zLyo3^xAJU?+_vwE?B@EOEM~ZO%OD?s2tzqv^)X&|+~a|5LD-I!R*|LtOf(&jpkwan
zL71WV-auowJF!^B)qGyS_Sjt$YNt`E)9x7Oo5(MX-`6^2JAY$PQ5UCkGQJZ-<uRAz
z@9qS5w&WIy8&HFLqHOZR-W`KQV{##KViFR@>cJg7)YfWl+xoX1UcC1OPF6naghnO_
zgu6kY(nC(JD+X?7f2n-xVQCUznj5c_24wceeBR|aF{<R<sWiL+^K}y7bu^T(Z1=^I
z-1Al=XHDYLHRje@YaHAB4C*orOth&>mC!v}J&l=00ZOU^J8s2H;M_cX?@6qoPz|Ha
z<~X14u801&RRY_3S3?(V#l^5Bs=}6F8Y12l$8PrNYT0tB_34j-nQfuBaBEmSj>j~&
z<}IDqNOR@rwRuX}9E~|y^&1rEd*oXa;^v>@EiMhsM}z1G=M^k5)h$%zvouxsQ)qU(
zL2M80SfG>uHjCBIX2Io##nx3ezv3tcuI`<(pBN&Do6ZpSistq92_fgGsFaW|ozk5x
zX|fTh(=O3vY8zZ$lPK49%}~+#ZtP!T5w;pDYYzw3i}P<-InNAMo#<2kgsTL+MGq;b
zrzw))%)eaT8z%^~uaExvrH|cRCn=f<{6}4wvB&Kb@eu_R(qtdG2htxZ(OJiCGkwcR
z?bNHz%v_=9Rch8hH>w?Z{)>ciHjgrNc-zgiqr2)lQ9lrsCtSwg%9NPd%H&Tvt*}%A
z2fk=UYmpf#Yl*VP)VVpWcHcFxbxwEeHZf~(&P-oE%GarN({aOdQkgN7eQ-0BC%MJf
z2YjlvlsZM*lp4aYn_*8l*=5U$*`MW50Gu%(YH&EsLWbakzqn{AI*^Fdx*8CeCx1NS
zhT^l$nPiHui+E#0<<R55s2$L$-6RuHFV>DCK-uGA5r>ngGMlEFPx3wDx_(vh%1!2~
zo>Q(5;Pt6w^Dgdi`3fGaQOcDc{8T!Xm&()m*%fLFsyq49<*pOQy?2)Jx}FBu-%HYk
zxN3oEIMu)zZX<-YDCZefi~H5rGX(FPUH)-~%gDX*+%=?jv^mH8^*t*2%X=?K?^E{S
zhpVaMbEnH1#wiJk=s&;nQ~lF0kdmtxyH!wfjR*;R1SiEmG9P^uY|$Xzpl-_H(1L*y
z%V>~him5UR1t?xBe8o`kD(?Lzjk+kE*VucArj8B2-&s@AIS~5Dxgg&?H;Vri-D;m<
zLj5zQs3Cq@LJ<<dRANdTi|DcH2|;qrkBS@ew%3~m0_^3D+I)KaXFg?=nTdQwQ*9Qb
z^epz-7#2U77v0?ujNT~^s1_}R^@#Pve}U|Cq9>5(9C7Wb$XrR)7krF|VQldiN*3ha
zLE>yOWN-QLQ{c-Ln8}EI@JK4)*l{|7SLyhb3Aosrad+?L>DEZWx{v;MO;X&=$GIt4
zqP?HbUs!|I!N7G$*^RR>eGiY>@MDUCH9S*Kz5|afT-z}V>tvtx^a;(lq`0*QE2Hk`
zi_^NF>}W;_SHa2Icsue%rlssEO>iYnLWzg>F(FxuYB=7U7fYyb25`bkrZgF0{oF~p
z42)V2d~Ie3gn_H0_^ja`cjUq=((Nqu(2=S<utL$*IFC_pKm>d({u4pisCvhw)^)Hf
z<yI*j81YW#Zjas$x<{b#n<I7Th(=bWIS+4@^Api#&^zdmr~b|r_B9thlHxo?bPW;S
zO$VAtE|GTP0(-(a!+Df`pKbaIG4^-47Iuu4f^XZ>5pnb6vz5;FyvIvR1)UWn{wfhl
zTNvbX_)lN13a|yOS(|p{A_WQ}>JAl9nXGXr-(8S+INRMM3j3WbJiHMJbAfXAJ?y7p
z2P2`Syc(s<-bqPeSv96dn%@e_iha1xk+YwOmnR_-sIXH=g&9U$-hW7rdWveZW&Tuh
zo}{E{3me9PK@i-GOuqkVHswGY4}bgOGDh1>7R3Ubwwa>KfGy0UEqE%JL@D9z>}qCa
zCYl&2Q^z)@UB)t|9VLhxJ5!&II7j0vdmdx6oaP=Xjroe{JB3At4u_cxQ{7oRsS#}s
zD!%(m4{=Iog-~O>GmF+O+BA`D4e_)t%J^J(ly8ba7>(9@-ecv!)X9t{<{=ckOBOP3
z@Pp)uNqoPo=GW?k)R2XX8R^c!Ri9bU^fYtE28=RA3=hM62h*!swnC(rGc&<JZTG?i
zXH2~n5SuD1t?5le@})bur=t@#Kc0(Utq57|G`fzKMF%Z6)x@Z`0`2C~OpfDJ+km~7
zaOUu)SZt$C<^0W|HoT#uja}Jd0AnMAz{Ie;mOyd-EltV|fed8KZc2zEtWgL!7U~@5
z#6Ic2xM2|Hx@GA_+ch;9`1U^gM;JA?vmyRK0N2~Ak#8SEeY|O+Z8s#f;5GyGb9i)P
z)IY-@OH`bhrE1MWe@1EQgW4AJJBju3Fb&d?@6OqryF=w#g*}YP!KSSjX148+HmL}L
z<ZR2im)p`;2-X4NLR9c5GNcbZ+^U5TI`d3Q|5B{AN#q4(m#|=%F%#0f7M?!`>+gJ?
z_qPn&gkr8&%#Y!Sm)X22Wp%%-Jz!EkYPtbhvveX`-;QiGZ;Z0i0+fd0vxbHDwES~1
zMH$W=Ab&W7c|3+{WyEDM_kD>NN}~b}*lfH-6hbWYd4L_XFs0vg>CTie=X*}89cNt-
z0zR~|ZXBH;LlUmO5Nx49s?Z%=#S4_nQbe{#F``3ydsq;&VkH{)P}v=neOUfFOv?|%
z6=uOCPZd#L3yaq8p)!j!#jT5WpiqrRe7HN;ASEkA5#<>Fa^@*x_SCD7f7whgAVgTO
zEFp9*IJg%6-I)@2P2zkNyq*90RXV6E<B!5S31$A?bjN}x<sOBO(di$e##f&e#q9aI
zBvEp63HWv{)E>X<VuwBdJm>*M*zUzLYYg!<wxaKxiS41?BzTQ=mCDJ~?kl)sM3Uc3
z>EJ)q1f)rdWG%mT%111U&hC6C#I&qFuKZd>!?o~2rq|P2XCu5`hu!oCEe{_*lDYLO
zo9IB~O(1<FF8mH&k;a-uZrhV-J~*hpy0E=Jfugg_!CGBWH#kxv@A}mFGaEtIH#6=S
zb2~!bW#qMsEMCf3xbfq~*u+fN@SwPEI9We?H~b9DrzXJOxBpRY!eEbroG@PLsxByi
zuiXuEeYYu7EXx6w5*&!B5Xl%G73J)waG*!aph~Au>$o0t>YTiY4M*pUf|qTm!J63p
zw`tB)3L;kJ=ZC*(ymA|%l#k_37fR6Eh_tR=2vCG2>3@r*8#NMvGUf{M3wazxGtF$x
zoGVa^i1(ayWMR$O!n!07#M)EW^0v(Z6&!SeWIM+BI_ht!5JJB;86}49e>L6QGa!Y~
zueKU@T*{vS4V1uAz=oT5&Px7)Qh|bdT@ua0d$d7hln*1^lNSP6IU*t56~5E)xY4qa
z0>JmbRaLDCHlGR)YI?a`R`Pn|B}TX;OO=t^qsBoLa}vopeffQ^N>=W9Tk)S|p5QXc
zzlfc{^zsh*@P?wuw+fq<*Qss@1G2`5g2{*_VDnGk6CL6k<RAuY=Ga_EQK7L8`yv3l
z-D<Z0lUI1mLa;!7q0wEo+CKNHXwgz)mP0AGTr7#^8h^JI{&CYKY?nk?g8H&Vgqg9r
z2Def%qoveCYllE)BdW_HyuNSBS~4-Y)dg;~O^wnwc6WhDFT>3#2<sjh$VQ-NO<Id4
zAw#?;kRL5FRgQuFCM=*Vorm0+4L6${egF1sNe&Q{!Mt)UMekAd+nxj<vbOyh-y#13
zS&4s8|1Xezj#B?Y{eK8#pBJ8~|9_pNKFfvvxuO3bw|^n?^AqAfxjlOkXM6E|0{j<U
CVZ`|W
literal 0
HcmV?d00001
diff --git a/NPAnalysis/S1554/Show.C b/NPAnalysis/S1554/Show.C
new file mode 100644
index 000000000..74b2446eb
--- /dev/null
+++ b/NPAnalysis/S1554/Show.C
@@ -0,0 +1,59 @@
+void Show(){
+
+ TFile* file = new TFile("../../Outputs/Analysis/BeamTest.root");
+ TTree* tree = (TTree*) file->FindObjectAny("PhysicsTree");
+
+ TCanvas* c = new TCanvas("Result","Result",600,600);
+ c->Divide(2,2);
+
+ c->cd(1);
+ tree->Draw("Ex:AddBack_DC/1000.>>h(500,0,10,500,-5,10)","ELab>0 && ThetaLab>90","colz");
+
+ c->cd(2);
+ tree->Draw("Ex>>hE(150,-5,10)","ELab>0 && ThetaLab>90");
+ tree->Draw("Ex>>hE2(150,-5,10)","ELab>0 && ThetaLab>90&&Trifoil.Time>50 && Trifoil.Time<65","same");
+ TH1F* hEx2 = (TH1F*) gDirectory->FindObjectAny("hE2");
+ hEx2->SetFillColor(kOrange+7);
+ hEx2->SetLineColor(kOrange+7);
+ c->cd(3);
+ tree->Draw("ThetaCM>>hCM(36,0,180)","ThetaLab>90 && Trifoil.Time>50 && Trifoil.Time<65 && Ex>-0.8 && Ex< 0.2","");
+ //tree->Draw("ThetaCM>>hCMb(36,0,180)","ThetaLab>90 && Ex>-0.8 && Ex< 0.2","same");
+
+ TH1F* hCM = (TH1F*) gDirectory->FindObjectAny("hCM");
+ //TH1F* hCMb = (TH1F*) gDirectory->FindObjectAny("hCMb");
+ // hCMb->SetMarkerColor(kGreen-3);
+// hCMb->Sumw2();
+
+ TFile* effFile = new TFile("SharcEff.root");
+ hCM->Sumw2();
+ TH1F* eff = (TH1F*) effFile->FindObjectAny("hDetecThetaCM");
+ eff->Sumw2();
+ eff->Rebin(5);
+ hCM->Divide(eff);
+ hCM->GetXaxis()->SetRangeUser(0,60);
+ gPad->SetLogy();
+
+ TGraph* g = new TGraph("22.si");
+ g->Draw("c");
+ double normal = hCM->GetBinContent(hCM->FindBin(8))/g->Eval(8);
+ TF1* f = new TF1("f",Form("%f",normal),0,180);
+ hCM->Divide(f);
+ /*double normalb = hCM->GetBinContent(hCM->FindBin(8))/hCMb->GetBinContent(hCMb->FindBin(8));
+ TF1* fb = new TF1("fb",Form("%f",1./normalb),0,180);
+ hCMb->Divide(fb);
+*/
+ c->cd(4);
+ tree->Draw("ThetaCM>>hCM2(36,0,180)","ThetaLab>90 && Trifoil.Time>50 && Trifoil.Time<65 && Ex>4 && Ex< 5","");
+ TH1F* hCM2 = (TH1F*) gDirectory->FindObjectAny("hCM2");
+ hCM2->Sumw2();
+ hCM2->Divide(eff);
+ hCM2->GetXaxis()->SetRangeUser(0,60);
+ gPad->SetLogy();
+
+ TGraph* g2 = new TGraph("22.sie");
+ g2->Draw("c");
+ normal = hCM2->GetBinContent(hCM->FindBin(8))/g2->Eval(8);
+ TF1* f2 = new TF1("f2",Form("%f",normal),0,180);
+ hCM2->Divide(f2);
+
+}
diff --git a/NPAnalysis/S1554/macro/Energy/EnergyCalibrator.cxx b/NPAnalysis/S1554/macro/Energy/EnergyCalibrator.cxx
new file mode 100644
index 000000000..483efc618
--- /dev/null
+++ b/NPAnalysis/S1554/macro/Energy/EnergyCalibrator.cxx
@@ -0,0 +1,64 @@
+void EnergyCalibrator(){
+ TFile* file = new TFile("~/Desktop/nptool/Outputs/Analysis/Sharc_Calib.root");
+ NPL::CalibrationSource* cs1 = new NPL::CalibrationSource();
+ NPL::EnergyLoss* el = new NPL::EnergyLoss("EnergyLossTable/alpha_Al.G4table","G4Table",10);
+
+ ofstream calfile("SharcCalibration.txt");
+ ofstream deadfile("DeadLayer.txt");
+
+ if(!deadfile.is_open() || !calfile.is_open())
+ exit(1);
+
+ vector<double> coeff;
+
+ for(unsigned int det = 0 ; det < 12 ; det++){
+ TH2* h2d = (TH2*) file->FindObjectAny(Form("SHARC%i_STR_FRONT_E_RAW",det+1));
+ if(h2d){
+ for(unsigned int i = 1 ; i < h2d->GetNbinsX() ; i++){
+ TH1D* h = h2d->ProjectionY("px",i,i);
+
+ NPL::CalibrationSource* cs = new NPL::CalibrationSource();
+ cs->Set_ThreeAlphaSource();
+
+ NPL::SiliconCalibrator* mycal = new NPL::SiliconCalibrator();
+ // double dist = mycal->ZeroExtrapolation(h,cs,el,coeff,0,5000,150000,300000);
+ double dist = mycal->SimpleCalibration(h,cs,el,coeff,150000,300000);
+
+ gPad->Update();
+ if(dist>=0){
+ deadfile << Form("SHARC_D%i_STRIP_FRONT%i_DEADLAYER ",det+1,i) << dist <<endl;
+ calfile << Form("SHARC_D%i_STRIP_FRONT%i_E ",det+1,i) << coeff[0] << " " << coeff[1] << endl;
+ }
+
+ delete h;
+ delete mycal;
+ delete cs;
+ }
+
+ h2d = (TH2*) file->FindObjectAny(Form("SHARC%i_STR_BACK_E_RAW",det+1));
+ if(h2d){
+ for(unsigned int i = 1 ; i < h2d->GetNbinsX() ; i++){
+ TH1D* h = h2d->ProjectionY("px",i,i);
+
+ NPL::CalibrationSource* cs = new NPL::CalibrationSource();
+ cs->Set_ThreeAlphaSource();
+
+ NPL::SiliconCalibrator* mycal = new NPL::SiliconCalibrator();
+ // double dist = mycal->ZeroExtrapolation(h,cs,el,coeff,0,5000,150000,300000);
+ double dist = mycal->SimpleCalibration(h,cs,el,coeff,150000,300000);
+
+ gPad->Update();
+ if(dist>=0){
+ deadfile << Form("SHARC_D%i_STRIP_BACK%i_DEADLAYER ",det+1,i) << dist <<endl;
+ calfile << Form("SHARC_D%i_STRIP_BACK%i_E ",det+1,i) << coeff[0] << " " << coeff[1] << endl;
+ }
+
+ delete h;
+ delete mycal;
+ delete cs;
+ }
+ }
+ }
+ }
+
+}
--
GitLab