KVIDChIoSi.cpp

/***************************************************************************
                          KVIDChIoSi.cpp  -  description
                             -------------------
    begin                : Fri Feb 20 2004
    copyright            : (C) 2004 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 "KVIDChIoSi.h"
#include "KVIDGChIoSi.h"
#include "KVIdentificationResult.h"
 
ClassImp(KVIDChIoSi)
//KVIDChIoSi
//
//Identification in ChIo-Si matrices of INDRA
//
 
//________________________________________________________________________________________//
 
 
 
Bool_t KVIDChIoSi::Identify(KVIdentificationResult* IDR, Double_t x, Double_t y)
{
   //Particle identification and code setting using identification grid KVIDGChIoSi
 
   Bool_t ok = KVINDRAIDTelescope::Identify(IDR, x, y);
 
   //if point lies above Zmax line, we give Zmax as Z of particle (real Z is >= Zmax)
   //general ID code = 5 (Zmin)
   if (IDR->IDquality == KVIDZAGrid::kICODE7) {
      IDR->IDcode = 5;//KVINDRA::ID_STOPPED_IN_FIRST_STAGE;
   }
   //Identified particles with subcode kID_LeftOfBragg are given
   //general ID code 5 (Zmin).
   if (IDR->IDquality == KVIDGChIoSi::k_LeftOfBragg) {
      IDR->IDcode = 5;//KVINDRA::ID_STOPPED_IN_FIRST_STAGE;
   }
   //unidentifiable particles with subcode kID_BelowSeuilSi are given
   //general ID code 5 (Zmin)
   if (IDR->IDquality == KVIDGChIoSi::k_BelowSeuilSi) {
      IDR->IDcode = 5;//KVINDRA::ID_STOPPED_IN_FIRST_STAGE;
      IDR->SetComment("point to identify left of Si threshold line (bruit/arret ChIo?)");
   }
   if (IDR->IDquality == KVIDGChIoSi::k_RightOfEmaxSi)
      IDR->SetComment("point to identify has E_Si > Emax_Si i.e. codeur is saturated. Unidentifiable");
 
   return kTRUE;
}
 
 
 
 
 
void KVIDChIoSi::Initialize()
{
   // Initialisation of telescope before identification.
 
   KVINDRAIDTelescope::Initialize();
   fchio = dynamic_cast<KVChIo*>(GetDetector(1));
}
 
 
 
 
Bool_t KVIDChIoSi::CheckTheoreticalIdentificationThreshold(KVNucleus* ION, Double_t)
{
   // Overrides KVIDTelescope method
   // We lower the calculated energy threshold to be better in agreement with data
   // EINC argument not used as ChIo-Si is always first telescope
 
   Double_t emin = (1.9 / 2.5) * (fchio->GetEIncOfMaxDeltaE(ION->GetZ(), ION->GetA()));
   return (ION->GetEnergy() > emin);
 
}