KVFAZIAIDSiCsI_FAZIASYM.cpp

Go to the documentation of this file.

//Created by KVClassFactory on Tue Sep  8 16:14:25 2015
//Author: ,,,
 
#include "KVFAZIAIDSiCsI_FAZIASYM.h"
#include "KVDataSet.h"
//#include <cstdio>
//#include <cstdlib>
 
ClassImp(KVFAZIAIDSiCsI_FAZIASYM)
 
// BEGIN_HTML <!--
/* -->
<h2>KVFAZIAIDSiCsI</h2>
<h4>id telescope to manage FAZIA Si-CsI identification</h4>
<!-- */
// --> END_HTML
 
 
 
KVFAZIAIDSiCsI_FAZIASYM::KVFAZIAIDSiCsI_FAZIASYM()
{
   // Default constructor
   SetType("Si-CsI");
   fBelowProton = 0;
   fSiThreshold = 0;
 
}
 
 
 
 
KVFAZIAIDSiCsI_FAZIASYM::~KVFAZIAIDSiCsI_FAZIASYM()
{
   // Destructor
}
 
 
 
 
 
Bool_t KVFAZIAIDSiCsI_FAZIASYM::Identify(KVIdentificationResult* idr, Double_t x, Double_t y)
{
   //Particle identification and code setting using identification grid (class KVIDZAGrid).
 
   idr->SetIDType(GetType());
   idr->IDattempted = kTRUE;
 
   //perform identification
   Double_t si2 = (y < 0. ? GetIDMapY() : y);
   Double_t csi = (x < 0. ? GetIDMapX() : x);
 
   //std::cout << "Inside KVFAZIAIDSiCsI_FAZIASYM::IDentify csi=" << csi << " si2=" << si2 << "\n";
   //test if line below proton and si threshold are there
   //if yes test the position of the point respect to
   //these lines
   //Bool_t OKproton = ((fBelowProton && fBelowProton->TestPoint(csi, si2)) || !fBelowProton);
   //Bool_t OKthreshold = ((fSiThreshold && fSiThreshold->TestPoint(csi, si2)) || !fSiThreshold);
   //if (OKproton && OKthreshold){
 
 
   if (fBelowProton) {
      if (fBelowProton->TestPoint(csi, si2)) idr->deltaEpedestal = KVIdentificationResult::deltaEpedestal_NO;
      else idr->deltaEpedestal = KVIdentificationResult::deltaEpedestal_YES;
   }
   else {
      idr->deltaEpedestal = KVIdentificationResult::deltaEpedestal_UNKNOWN;
   }
 
   if (TheGrid->IsIdentifiable(csi, si2)) {
      // std::cout << "Identifico con TheGrid\n";
      TheGrid->Identify(csi, si2, idr);
      //  cout << "Ho identificato con TheGrid" << endl;
      //   printf ("The Grid: %d %d %f\n", idr->Z, idr->A, idr->PID);
 
      //utilizzo TheGrid_lcp per identificare nuovamente Z=1 e Z=2
      if (TheGrid_lcp->IsIdentifiable(csi, si2) && idr->Z <= 2 && idr->Z >= 0) {
         //std::cout << "Identifico con TheGrid_lcp\n";
         TheGrid_lcp->Identify(csi, si2, idr);
         // cout << "Ho identificato con TheGrid_lcp" << endl;
         //printf ("The Grid_lcp: %d %d %f\n", idr->Z, idr->A, idr->PID);
      }
   }
   else {
      idr->IDOK = kFALSE;
      idr->IDquality = KVIDZAGrid::kICODE8;
   }
 
 
 
   // set general ID code
   idr->IDcode = GetIDCode();
 
   return kTRUE;
 
}
 
 
 
 
 
void KVFAZIAIDSiCsI_FAZIASYM::Initialize()
{
   // Initialisation of telescope before identification.
   // This method MUST be called once before any identification is attempted.
   // Initialisation of grid is performed here.
   // IsReadyForID() will return kTRUE if a grid is associated to this telescope for the current run.
   //std::cout << "Inizialize FAZIASYM SiCsI" << std::endl;
   //TheGrid = (KVIDZAGrid*) GetIDGrid();
   TheGrid = (KVIDZAGrid*) GetListOfIDGrids()->At(0);
   TheGrid_lcp = (KVIDZAGrid*) GetListOfIDGrids()->At(1);
   if (TheGrid) {
      // std::cout << "Inizializzo TheGrid" << std::endl;
      SetHasMassID(TheGrid->HasMassIDCapability());
      TheGrid->Initialize();
      fBelowProton = (KVIDCutContour*)TheGrid->GetCut("PIEDESTAL");//fBelowProton = (KVIDCutLine*)TheGrid->GetCut("PIEDESTAL");
      fSiThreshold = (KVIDCutLine*)TheGrid->GetCut("threshold");
      SetBit(kReadyForID);
   }
   else {
      ResetBit(kReadyForID);
   }
 
   if (TheGrid_lcp) {
      //std::cout << "Inizializzo TheGrid_lcp" << std::endl;
      SetHasMassID(TheGrid_lcp->HasMassIDCapability());
      TheGrid_lcp->Initialize();
      fBelowProton = (KVIDCutContour*)TheGrid->GetCut("PIEDESTAL");//fBelowProton = (KVIDCutLine*)TheGrid->GetCut("PIEDESTAL");
      fSiThreshold = (KVIDCutLine*)TheGrid_lcp->GetCut("threshold");
      SetBit(kReadyForID);
   }
   else {
      ResetBit(kReadyForID);
   }
 
   if (!gDataSet->HasCalibIdentInfos()) {// for filtering simulations
      SetBit(kReadyForID);
      SetHasMassID();// in principle mass identification always possible
   }
}
 
 
 
 
void KVFAZIAIDSiCsI_FAZIASYM::SetIdentificationStatus(KVReconstructedNucleus* n)
{
   // For filtering simulations
   //
   // Z-dependence of A identification:
   //    all ok if Z<=14, decreasing probability for 15<=Z<=18
   //    no A identification for Z>18
 
   n->SetZMeasured();
   fMassIDProb->SetParameters(16.5, .4);
   Bool_t okmass = (n->GetZ() <= 14) || (n->GetZ() < 19 && gRandom->Uniform() < fMassIDProb->Eval(n->GetZ()));
   if (okmass) {
      n->SetAMeasured();
   }
   else
      n->SetZ(n->GetZ());
}