KVIDZAFromZGrid.h

 
#ifndef __KVIDZAFROMZGRID_H
#define __KVIDZAFROMZGRID_H
 
#include "KVIDZAGrid.h"
#include "KVList.h"
#include "KVUniqueNameList.h"
#include "KVIdentificationResult.h"
#include "KVMultiGaussIsotopeFit.h"
 
 
class interval;
class interval_set;
 
class KVIDZAFromZGrid : public KVIDZAGrid {
 
public:
   enum PIDType {
      kNone,
      kPeackType,
      kIntType
   };
 
protected:
   Int_t  fZmaxInt;
   Bool_t fPIDRange;
   Int_t  fZminInt;
    KVList fTables;
    KVUniqueNameList fFits;
 
   Bool_t fIgnoreMassID;
    Bool_t fHasMassIDRegion;
 
   int is_inside(double pid) const;
 
public:
   KVIDZAFromZGrid();
   ROOT_COPY_CTOR(KVIDZAFromZGrid, KVIDZAGrid)
   ROOT_COPY_ASSIGN_OP(KVIDZAFromZGrid)
 
   void Copy(TObject& obj) const;
 
   virtual void ReadFromAsciiFile(std::ifstream& gridfile);
   virtual void WriteToAsciiFile(std::ofstream& gridfile);
 
   virtual void Identify(Double_t x, Double_t y, KVIdentificationResult*) const;
   virtual double DeduceAfromPID(KVIdentificationResult* idr) const;
   void LoadPIDRanges();
   void ResetPIDRanges();
   void ReloadPIDRanges();
   interval_set* GetIntervalSet(int zint) const;
   KVList* GetIntervalSets()
   {
      return &fTables;
   }
 
   void PrintPIDLimits();
   void ClearPIDIntervals();
   void ExportToGrid();
 
   void Initialize();
 
   Bool_t HasMassIDCapability() const
   {
 
      return (fPIDRange || fTables.GetEntries()) && !fIgnoreMassID;
   }
   void SetOnlyZId(Bool_t /*yes*/ = kTRUE)
   {
 
      KVIDZAGrid::SetOnlyZId(true);
   }
 
   bool MassIdentificationFromMultiGaussFit(KVMultiGaussIsotopeFit*, KVIdentificationResult*) const;
   KVMultiGaussIsotopeFit* GetMultiGaussFit(int z) const
   {
      return (KVMultiGaussIsotopeFit*)fFits.FindObject(Form("MultiGaussIsotopeFit_Z=%d", z));
   }
 
   ClassDef(KVIDZAFromZGrid, 1) //Compute Z and A only from Z lines...
};
 
 
 
class interval: public TNamed {
   int fType;
   int fZ;
   int fA;
 
   double fPID;
   double fPIDMin;
   double fPIDmax;
 
public:
   void Copy(TObject& o) const
   {
      TNamed::Copy(o);
      interval& i = dynamic_cast<interval&>(o);
      i.fType = fType;
      i.fZ = fZ;
      i.fA = fA;
      i.fPID = fPID;
      i.fPIDMin = fPIDMin;
      i.fPIDmax = fPIDmax;
   }
   ROOT_DEF_CTOR(interval, TNamed)
   ROOT_COPY_CTOR(interval, TNamed)
   ROOT_COPY_ASSIGN_OP(interval)
   interval(int zz, int aa, double pid, double pidmin = -1., double pidmax = -1.)
   {
      fZ = zz;
      fA = aa;
      fPID = pid;
      fPIDMin = pidmin;
      fPIDmax = pidmax;
      SetName(Form("%d_%d", zz, aa));
   }
   bool is_inside(double pid)
   {
      if (pid > fPIDMin && pid < fPIDmax) return kTRUE;
      return kFALSE;
   }
   bool is_left_of(double pid)
   {
      return (fPIDmax < pid);
   }
   bool is_right_of(double pid)
   {
      return (fPIDMin > pid);
   }
 
   int    GetA()
   {
      return fA;
   }
   int    GetZ()
   {
      return fZ;
   }
   double GetPID()
   {
      return fPID;
   }
   double GetPIDmin()
   {
      return fPIDMin;
   }
   double GetPIDmax()
   {
      return fPIDmax;
   }
 
   void SetA(int aa)
   {
      fA = aa;
   }
   void SetZ(int zz)
   {
      fZ = zz;
   }
   bool SetPID(double pid)
   {
      fPID = pid;
      if (fPIDMin > 0 && fPIDmax > 0) {
         if (pid >= fPIDMin && pid <= fPIDmax) return true;
         return false;
      }
      return true;
   }
   void SetPIDmin(double pidmin)
   {
      fPIDMin = pidmin;
   }
   void SetPIDmax(double pidmax)
   {
      fPIDmax = pidmax;
   }
 
 
   ClassDef(interval, 1) //
};
 
class interval_set: public TNamed {
public:
   int fType;
   int fZ;
   KVList fIntervals;
   Int_t  fNPIDs;
   TGraph fPIDs;
 
   void Copy(TObject& o) const
   {
      TNamed::Copy(o);
      interval_set& i = dynamic_cast<interval_set&>(o);
      i.fType = fType;
      i.fZ = fZ;
      fIntervals.Copy(i.fIntervals);
      i.fNPIDs = fNPIDs;
    i.fPIDs = fPIDs;
   }
   ROOT_DEF_CTOR(interval_set, TNamed)
   ROOT_COPY_CTOR(interval_set, TNamed)
   ROOT_COPY_ASSIGN_OP(interval_set)
   int GetZ()
   {
      return fZ;
   }
   int GetNPID()
   {
      return fIntervals.GetSize();
   }
   KVList* GetIntervals()
   {
      return &fIntervals;
   }
   TString GetListOfMasses();
 
   interval_set(int zz, int type);
   void   add(int aa, double pid, double pidmin = -1., double pidmax = -1.);
   double eval(KVIdentificationResult* idr);
   bool is_inside(double pid);
   bool is_above(double pid);
 
   ClassDef(interval_set, 1) //
};
 
 
 
 
 
#endif