* adding threshold to CsI and SiLi

NPSimulation/Detectors/MUST2/MUST2Array.cc

@@ -939,11 +939,13 @@ void MUST2Array::ReadSensitive(const G4Event*) {
     for (unsigned int i = 0; i < sizeSiLi; i++) {
       double ESiLi = RandGauss::shoot(SiLiScorer->GetEnergy(i), ResoSiLi);
       vector<unsigned int> level = SiLiScorer->GetLevel(i);
+      if(ESiLi>ThresholdSiLi){
       m_Event->SetSiLiE(level[0], level[1],
                         NPL::EnergyToADC(ESiLi, 0, 250, 8192, 16384));
       double timeSiLi = RandGauss::shoot(SiLiScorer->GetTime(i), ResoTimeMust);
       m_Event->SetSiLiT(level[0], level[1],
                         NPL::EnergyToADC(timeSiLi, 0, 1000, 16384, 8192));
+      }
     }
 
     // CsI //
@@ -954,11 +956,13 @@ void MUST2Array::ReadSensitive(const G4Event*) {
     for (unsigned int i = 0; i < sizeCsI; i++) {
       double ECsI = RandGauss::shoot(CsIScorer->GetEnergy(i), ResoCsI);
       vector<unsigned int> level = CsIScorer->GetLevel(i);
+      if(ECsI>ThresholdCsI){
       m_Event->SetCsIE(level[0], level[1],
                        NPL::EnergyToADC(ECsI, 0, 250, 8192, 16384));
       double timeCsI = RandGauss::shoot(CsIScorer->GetTime(i), ResoTimeMust);
       m_Event->SetCsIT(level[0], level[1],
                        NPL::EnergyToADC(timeCsI, 0, 1000, 16384, 8192));
+      }
     }
   }
 }

NPSimulation/Detectors/MUST2/MUST2Array.hh

@@ -40,7 +40,10 @@ namespace MUST2
    const G4double ResoCsI      = 0.08           ;// = 188 kev of resolution //   Unit is MeV/2.35
    const G4double ResoStrip    = 0.0149         ;// 0.0223 = 52keV of Resolution   //   Unit is MeV/2.35  14.861996
    const G4double TimeOffset   = 500            ;// 500 ns stop
-   const G4double ThresholdSi  = 500 * keV;
+   // Threshold
+   const G4double ThresholdSi   = 500 * keV;
+   const G4double ThresholdSiLi = 500 * keV;
+   const G4double ThresholdCsI  = 1000 * keV;
    // Geometry
    const G4double FaceFront = 11.*cm ;
    const G4double FaceBack = 16.5*cm ;