KVFAZIAIDTelescope.cpp

//Created by KVClassFactory on Mon Feb 17 13:51:39 2014
//Author: John Frankland,,,
 
#include "KVFAZIAIDTelescope.h"
 
TF1* KVFAZIAIDTelescope::fMassIDProb = 0;
 
ClassImp(KVFAZIAIDTelescope)
 
 
// BEGIN_HTML <!--
/* -->
<h2>KVFAZIAIDTelescope</h2>
<h4>Identification for FAZIA array</h4>
<!-- */
// --> END_HTML
 
 
 
 
KVFAZIAIDTelescope::KVFAZIAIDTelescope()
{
   // Default constructor
   if (!fMassIDProb) {
      fMassIDProb = new TF1("FAZIA-MassIDProb", "1./(exp((x-[0])/[1])+1)", 0, 100);
      fMaxZ = 22.5;
      fSigmaZ = .5;
//      fMassIDProb->SetParameters(22.5, .5);
   }
}
 
 
 
 
void KVFAZIAIDTelescope::AddDetector(KVDetector* d)
{
   // Add a detector to the telescope.
   //
   // Detectors must be added in the order they will be hit by impinging particles,
   // with the last detector being the one particles stopped in the telescope will stop in.
   // i.e. dE1, dE2, ..., Eres
   //
   // Update name of telescope to "ID_[label of 1st detector]_[label of 2nd detector]_ ... _[label of last detector]"
 
   if (d) {
      fDetectors.Add(d);
      d->AddIDTelescope(this);
      if (GetSize() > 1) {
         TString name = "ID_";
         KVDetector* det = nullptr;
         for (int i = 1; i <= GetSize(); ++i) {
            det = GetDetector(i);
            name += det->GetLabel();
            name += "_";
         }
         name += Form("%d", det->GetIndex());
         SetName(name);
      }
      else
         SetName(Form("ID_%s_%d", GetDetector(1)->GetLabel(), GetDetector(1)->GetIndex()));
   }
   else {
      Warning("AddDetector", "Called with null pointer");
   }
}
 
 
//const Char_t* KVFAZIAIDTelescope::GetNewName(KVString oldname)
//{
 
//   KVString tmp = "";
//   KVString lab = "";
//   oldname.Begin("_");
//   KVString id = oldname.Next();
//   if (id != "ID") return "unkonwn";
//   if (oldname.End()) return "unkonwn";
//   KVString det1 = oldname.Next();
//   KVString det2 = "";
//   if (!oldname.End()) {
//      det2 = oldname.Next();
//   }
//   KVString newdet1 = KVFAZIADetector::GetNewName(det1.Data());
//   newdet1.Begin("-");
//   KVString labdet1 = newdet1.Next();
//   KVString idxdet1 = newdet1.Next();
 
//   static KVString newname;
//   if (det2 == "") {
//      newname.Form("ID_%s_%s", labdet1.Data(), idxdet1.Data());
//   }
//   else {
//      KVString newdet2 = KVFAZIADetector::GetNewName(det2.Data());
//      newdet2.Begin("-");
//      KVString labdet2 = newdet2.Next();
//      newname.Form("ID_%s_%s_%s", labdet1.Data(), labdet2.Data(), idxdet1.Data());
//   }
 
//   return newname.Data();
//}
 
 
 
void KVFAZIAIDTelescope::SetIdentificationStatus(KVReconstructedNucleus* n)
{
   // For filtering simulations
   //
   // Z-dependence of A identification:
   // fMassIDProb parameters has to set in the Initialize method
 
 
   fMassIDProb->SetParameters(fMaxZ, fSigmaZ);
 
//    Info("SetIdentificationStatus","%s : %lf %lf",ClassName(),fMassIDProb->GetParameter(0),fMassIDProb->GetParameter(1));
 
   n->SetZMeasured();
   Bool_t okmass = (gRandom->Uniform() < fMassIDProb->Eval(n->GetZ())) && CanIdentify(n->GetZ(), n->GetA());
 
   if (okmass) {
      //reset A to the original mass in case of multiple call of this method
      if (n->GetParameters()->HasParameter("OriginalMass")) n->SetA(n->GetParameters()->GetIntValue("OriginalMass"));
      n->SetAMeasured(kTRUE);
   }
   else {
      //save the original mass in the parameter list in case of multiple call of this method
      n->GetParameters()->SetValue("OriginalMass", n->GetA());
      double e = n->GetE();
      n->SetZ(n->GetZ());
      n->SetE(e);
      n->SetAMeasured(kFALSE);
   }
}