Anp::FillBasic Class Reference

#include <FillBasic.h>

List of all members.

Public Member Functions
	FillBasic ()
	~FillBasic ()
void	Config (const Registry &reg)
bool	Fill (const NtpStRecord &ntprec, const NtpSREvent &record)
bool	Fill (const NtpStRecord &ntprec, const NtpSRShower &record)
bool	Fill (const NtpStRecord &ntprec, const NtpSRTrack &record)
const Basic	GetBasic () const
const Vertex	GetBegVtx () const
const Vertex	GetEndVtx () const
Private Member Functions
const Vertex	Get (const NtpSRVertex &vtx, bool isfid_) const
bool	Check (const Basic &basic, const NtpSRPlane &plane) const
bool	Fill (Basic &basic, const NtpSRStripPulseHeight &ph) const
const Basic	Get (const NtpStRecord &record, const int *index_array, int nstrip) const
Private Attributes
bool	fCheck
double	fTrackZOffset
Basic	fBasic
Vertex	fBegVtx
Vertex	fEndVtx

---

Detailed Description

Definition at line 27 of file FillBasic.h.

---

Constructor & Destructor Documentation

Anp::FillBasic::FillBasic

(

)

Definition at line 32 of file FillBasic.cxx.

   :fCheck(true),
    fTrackZOffset(0.0392)
{    
}

Anp::FillBasic::~FillBasic

(

)

Definition at line 39 of file FillBasic.cxx.

{
}

---

Member Function Documentation

bool Anp::FillBasic::Check	(	const Basic &	basic,
		const NtpSRPlane &	plane
	)			const [private]

Definition at line 219 of file FillBasic.cxx.

References NtpSRPlane::beg, Anp::Basic::beg_uplane, Anp::Basic::beg_vplane, NtpSRPlane::begu, NtpSRPlane::begv, NtpSRPlane::end, Anp::Basic::end_uplane, Anp::Basic::end_vplane, NtpSRPlane::endu, NtpSRPlane::endv, Msg::kWarning, MSG, NtpSRPlane::n, NtpSRPlane::nu, Anp::Basic::nuplane, NtpSRPlane::nv, and Anp::Basic::nvplane.

Referenced by Fill().

{
   //
   // Check consistency of Basic object, this is useful bug check.
   //

   bool result = true;

   if(plane.nu + plane.nv != plane.n)
   {
      MSG("FillAlg", Msg::kWarning) << "Number of planes does not match" << endl;
      result = false;
   }
   if(basic.nuplane != short(plane.nu))
   {
      MSG("FillAlg", Msg::kWarning) << "Mismatched number of U planes" << endl;
      result = false;
   }
   if(basic.nvplane != short(plane.nv))
   {
      MSG("FillAlg", Msg::kWarning) << "Mismatched number of V planes" << endl;
      result = false;
   }
   
   if(plane.beg < plane.end)
   {
      if(basic.beg_uplane != short(plane.begu))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched beg U plane" << endl;
         result = false;
      }
      if(basic.beg_vplane != short(plane.begv))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched beg V plane" << endl;
         result = false;
      }
      if(basic.end_uplane != short(plane.endu))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched end U plane" << endl;
         result = false;
      }  
      if(basic.end_vplane != short(plane.endv))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched end V plane" << endl;
         result = false;
      }
   }
   else if(plane.beg > plane.end)
   {    
      if(basic.beg_uplane != short(plane.endu))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched beg U plane" << endl;
         result = false;
      }
      if(basic.beg_vplane != short(plane.endv))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched beg V plane" << endl;
         result = false;
      }
      if(basic.end_uplane != short(plane.begu))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched end U plane" << endl;
         result = false;
      }  
      if(basic.end_vplane != short(plane.begv))
      {
         MSG("FillAlg", Msg::kWarning) << "Mismatched end V plane" << endl;
         result = false;
      }
   }
   else
   {
      MSG("FillAlg", Msg::kWarning) << "beg and end plane are equal" << endl;
      result = false;
   }

   return result;
}

void Anp::FillBasic::Config

(

const Registry &

reg

)

Definition at line 44 of file FillBasic.cxx.

References fCheck, fTrackZOffset, Registry::Get(), Registry::KeyExists(), and Anp::Read().

Referenced by Anp::FillTrack::Config(), Anp::FillEvent::Config(), and Anp::FillShower::Config().

{
   //
   // Configure self
   //
   Anp::Read(reg, "FillBasicCheck", fCheck);

   reg.Get("TrackZOffset", fTrackZOffset);

   if(reg.KeyExists("PrintConfigBasic"))
   {
      cout << "FillBasic::Config" << endl
           << "   Check = " << fCheck << endl
           << "   TrackZOffset = " << fTrackZOffset << endl;
   }
}

bool Anp::FillBasic::Fill	(	Basic &	basic,
		const NtpSRStripPulseHeight &	ph
	)			const [private]

Definition at line 185 of file FillBasic.cxx.

References Anp::Basic::adcu, Anp::Basic::adcv, count, Msg::kWarning, NtpSRStripPulseHeight::mip, Anp::Basic::mip, MSG, NtpSRPulseHeight::pe, Anp::Basic::pe, NtpSRPulseHeight::raw, NtpSRPulseHeight::sigcor, Anp::Basic::sigcoru, Anp::Basic::sigcorv, NtpSRPulseHeight::siglin, Anp::Basic::siglin, NtpSRStripPulseHeight::sigmap, and Anp::Basic::sigmap.

{
   //
   // Fill Basic with pulse height information and check consistency between standard
   // pulse height class NtpSRStripPulseHeight and pulse height computed by iterating over
   // NtpSRStrip objects that belong to reconstructed object.
   //

   basic.siglin = ph.siglin;
   basic.pe     = ph.pe;
   basic.sigmap = ph.sigmap;
   basic.mip    = ph.mip;

   if(std::fabs(basic.adcu + basic.adcv - ph.raw) > 1.001)
   {
      static unsigned int count = 0;
      MSG("FillAlg", Msg::kWarning) << "Error " << ++count << " in NtpSRStripPulseHeight adc "
                                    << basic.adcu + basic.adcv << " - " << ph.raw << " = "
                                    << basic.adcu + basic.adcv - ph.raw << std::endl;
      return false;
   }
   if(std::fabs(basic.sigcoru + basic.sigcorv - ph.sigcor) > 1.001)
   {
      static unsigned int count = 0;
      MSG("FillAlg", Msg::kWarning) << "Error " << ++count << " in NtpSRStripPulseHeight sigcor " 
                                    << basic.sigcoru + basic.sigcorv << " - " << ph.sigcor << " = "
                                    << basic.sigcoru + basic.sigcorv - ph.sigcor << std::endl;
      return false;
   }

   return true;
}

bool Anp::FillBasic::Fill	(	const NtpStRecord &	ntprec,
		const NtpSRTrack &	record
	)

Definition at line 132 of file FillBasic.cxx.

References Check(), Anp::Vertex::Clear(), Anp::Basic::Clear(), NtpSRTrack::end, fBasic, fBegVtx, fCheck, fEndVtx, Fill(), fTrackZOffset, Get(), VldContext::GetDetector(), RecRecordImp< T >::GetHeader(), VldContext::GetSimFlag(), RecHeader::GetVldContext(), infid(), Msg::kWarning, MSG, NtpSRTrack::nstrip, NtpSRTrack::ph, NtpSRTrack::plane, NtpSRTrack::stp, NtpSRTrack::vtx, NtpSRVertex::x, NtpSRVertex::y, and NtpSRVertex::z.

{
   //
   // Fill Anp::Basic and Anp::Vertex data
   //

   fBasic.Clear();
   fBegVtx.Clear();
   fEndVtx.Clear();

   const bool isfid = infid(ntprec.GetHeader().GetVldContext().GetDetector(),
                            ntprec.GetHeader().GetVldContext().GetSimFlag(),
                            record.vtx.x,
                            record.vtx.y,
                            record.vtx.z - fTrackZOffset);
   

   fBasic  = FillBasic::Get(ntprec, record.stp, record.nstrip);
   fBegVtx = FillBasic::Get(record.vtx, isfid);
   fEndVtx = FillBasic::Get(record.end, false);

   if(!Fill(fBasic, record.ph))
   {      
      MSG("FillAlg", Msg::kWarning) << "Fill failed for NtpSREvent \n" << record;
      return false;
   }

   if(fCheck && !Check(fBasic, record.plane))
   {
      MSG("FillAlg", Msg::kWarning) << "NtpSRTrack plane check failed \n" << record.plane;
      return false;
   }

   return true;
}

bool Anp::FillBasic::Fill	(	const NtpStRecord &	ntprec,
		const NtpSRShower &	record
	)

Definition at line 97 of file FillBasic.cxx.

References Check(), Anp::Vertex::Clear(), Anp::Basic::Clear(), fBasic, fBegVtx, fCheck, fEndVtx, Fill(), Get(), VldContext::GetDetector(), RecRecordImp< T >::GetHeader(), VldContext::GetSimFlag(), RecHeader::GetVldContext(), infid(), Msg::kWarning, MSG, NtpSRShower::nstrip, NtpSRShower::ph, NtpSRShower::plane, NtpSRShower::stp, NtpSRShower::vtx, NtpSRVertex::x, NtpSRVertex::y, and NtpSRVertex::z.

{
   //
   // Fill Anp::Basic and Anp::Vertex data
   //

   fBasic.Clear();
   fBegVtx.Clear();
   fEndVtx.Clear();

   const bool isfid = infid(ntprec.GetHeader().GetVldContext().GetDetector(),
                            ntprec.GetHeader().GetVldContext().GetSimFlag(),
                            record.vtx.x,
                            record.vtx.y,
                            record.vtx.z);
   
   fBasic  = FillBasic::Get(ntprec, record.stp, record.nstrip);
   fBegVtx = FillBasic::Get(record.vtx, isfid);

   if(!Fill(fBasic, record.ph))
   {      
      MSG("FillAlg", Msg::kWarning) << "Fill failed for NtpSRShower \n" << record;
      return false;
   }

   if(fCheck && !Check(fBasic, record.plane))
   {
      MSG("FillAlg", Msg::kWarning) << "NtpSRShower plane check failed \n" << record.plane;
      return false;
   }

   return true;
}

bool Anp::FillBasic::Fill	(	const NtpStRecord &	ntprec,
		const NtpSREvent &	record
	)

Definition at line 62 of file FillBasic.cxx.

References Check(), Anp::Vertex::Clear(), Anp::Basic::Clear(), fBasic, fBegVtx, fCheck, fEndVtx, Get(), VldContext::GetDetector(), RecRecordImp< T >::GetHeader(), VldContext::GetSimFlag(), RecHeader::GetVldContext(), infid(), Msg::kWarning, MSG, NtpSREvent::nstrip, NtpSREvent::ph, NtpSREvent::plane, NtpSREvent::stp, NtpSREvent::vtx, NtpSRVertex::x, NtpSRVertex::y, and NtpSRVertex::z.

Referenced by Fill(), Anp::FillShower::Run(), Anp::FillTrack::Run(), and Anp::FillEvent::Run().

{
   //
   // Fill Anp::Basic and Anp::Vertex data
   //

   fBasic.Clear();
   fBegVtx.Clear();
   fEndVtx.Clear();

   const bool isfid = infid(ntprec.GetHeader().GetVldContext().GetDetector(),
                            ntprec.GetHeader().GetVldContext().GetSimFlag(),
                            record.vtx.x,
                            record.vtx.y,
                            record.vtx.z);

   fBasic  = FillBasic::Get(ntprec, record.stp, record.nstrip);
   fBegVtx = FillBasic::Get(record.vtx, isfid);
   
   if(!Fill(fBasic, record.ph))
   {      
      MSG("FillAlg", Msg::kWarning) << "Fill failed for NtpSREvent \n" << record;
      return false;
   }

   if(fCheck && !Check(fBasic, record.plane))
   {
      MSG("FillAlg", Msg::kWarning) << "NtpSREvent plane check failed \n" << record.plane;
      return false;
   }

   return true;
}

const Anp::Basic Anp::FillBasic::Get	(	const NtpStRecord &	record,
		const int *	index_array,
		int	nstrip
	)			const [private]

Definition at line 299 of file FillBasic.cxx.

References Anp::Basic::active_frac, Anp::Basic::adcu, Anp::Basic::adcv, Anp::Basic::beg_uplane, Anp::Basic::beg_vplane, count, Anp::Basic::end_uplane, Anp::Basic::end_vplane, find(), VldContext::GetDetector(), RecRecordImp< T >::GetHeader(), Anp::GetSRStrip(), RecHeader::GetVldContext(), it, Detector::kFar, Detector::kNear, PlaneView::kU, PlaneView::kV, Msg::kWarning, max, Anp::Basic::max_time, min, Anp::Basic::min_time, MSG, Anp::Basic::nuplane, Anp::Basic::nustrip, Anp::Basic::nvplane, Anp::Basic::nvstrip, Anp::Basic::sigcoru, Anp::Basic::sigcorv, and NtpStRecord::stp.

{  
   //
   // Iterate over strips with indexes stored in index_array and
   // compute pulse height and time and plane spans.
   //

   Basic data;

   if(nstrip < 1 || !index_array)
   {
      return data;
   }

   //
   // Initialize variables that use += operator
   //
   data.nustrip = 0;
   data.nvstrip = 0;   
   data.adcu = 0.0;
   data.adcv = 0.0;
   data.sigcoru = 0.0;
   data.sigcorv = 0.0;
   
   const TClonesArray *strip_array = record.stp;
   if(!strip_array)
   {
      MSG("FillAlg", Msg::kWarning) << "NtpStRecord does not have valid strip array." << endl; 
      return data;
   }

   vector<short> uplane_vec, vplane_vec;

   const int entries = strip_array -> GetEntries();

   const Detector::Detector_t detector = record.GetHeader().GetVldContext().GetDetector();

   bool init_time = false;
   for(int i = 0; i < nstrip; ++i)
   {            
      const short index = index_array[i];
      if(index < 0 || index >= entries)
      {
         MSG("FillAlg", Msg::kWarning) << "NtpSRStrip index is out of range: " << index << endl;
         continue; 
      }

      NtpSRStrip *ntpstp = Anp::GetSRStrip(record, index);
      if(!ntpstp)
      {
         continue;
      }
    
      if(index != ntpstp -> index)
      {  
         MSG("FillAlg", Msg::kWarning) << "Mismatched strip and TClonesArray index" << endl;
         continue;
      }

      const short plane = ntpstp -> plane;      

      if(ntpstp -> planeview == PlaneView::kU)
      {
         if(std::find(uplane_vec.begin(), uplane_vec.end(), plane) == uplane_vec.end())
         {
            uplane_vec.push_back(plane);
         }
      }
      else if(ntpstp -> planeview == PlaneView::kV)
      {
         if(std::find(vplane_vec.begin(), vplane_vec.end(), plane) == vplane_vec.end())
         {
            vplane_vec.push_back(plane);
         }
      }
      else
      {
         MSG("FillAlg", Msg::kWarning) << "Unknown planeview " << int(ntpstp -> planeview) << endl;
         continue;
      }

      if(detector == Detector::kNear)
      {
         if(ntpstp -> pmtindex0 > 0)
         {
            MSG("FillAlg", Msg::kWarning) << "In ND East PMT index is positive" << std::endl;
            continue;   
         }

         if(ntpstp -> planeview == PlaneView::kU)
         {
            ++data.nustrip;
            data.adcu += ntpstp -> ph1.raw; 
            data.sigcoru += ntpstp -> ph1.sigcor; 
         }
         else if(ntpstp -> planeview == PlaneView::kV)
         {
            ++data.nvstrip;         
            data.adcv += ntpstp -> ph1.raw;
            data.sigcorv += ntpstp -> ph1.sigcor; 
         }

         if(!init_time)
         {
            data.max_time = ntpstp -> time1;
            data.min_time = ntpstp -> time1;
            init_time = true;
         }
         else
         {
            data.max_time = std::max(ntpstp -> time1, data.max_time);
            data.min_time = std::min(ntpstp -> time1, data.min_time);
         }
      }
      else if(detector == Detector::kFar)      
      {
         if(ntpstp -> planeview == PlaneView::kU)
         {
            ++data.nustrip;

            if(ntpstp -> pmtindex0 > 0 ||
               (ntpstp -> pmtindex0 == 0 && ntpstp -> ph0.raw > 0.0))
            {
               data.adcu += ntpstp -> ph0.raw;
               data.sigcoru += ntpstp -> ph0.sigcor; 
            }
            
            if(ntpstp -> pmtindex1 > 0 ||
               (ntpstp -> pmtindex1 == 0 && ntpstp -> ph1.raw > 0.0))
            {
               data.adcu += ntpstp -> ph1.raw;
               data.sigcoru += ntpstp -> ph1.sigcor; 
            }
         }
         else if(ntpstp -> planeview == PlaneView::kV)
         {
            ++data.nvstrip;

            if(ntpstp -> pmtindex0 > 0 ||
               (ntpstp -> pmtindex0 == 0 && ntpstp -> ph0.raw > 0.0))
            {
               data.adcv += ntpstp -> ph0.raw;
               data.sigcorv += ntpstp -> ph0.sigcor; 
            }

            if(ntpstp -> pmtindex1 > 0 ||
               (ntpstp -> pmtindex1 == 0 && ntpstp -> ph1.raw > 0.0))
            {
               data.adcv += ntpstp -> ph1.raw;
               data.sigcorv += ntpstp -> ph1.sigcor; 
            }
         }


         if(ntpstp -> pmtindex0 > 0 ||
            (ntpstp -> pmtindex0 == 0 && ntpstp -> ph0.raw > 0.0))
         {
            if(!init_time)
            {
               data.max_time = ntpstp -> time0;
               data.min_time = ntpstp -> time0;
               init_time = true;
            }
            else
            {
               data.max_time = std::max(data.max_time, ntpstp -> time0);
               data.min_time = std::min(data.min_time, ntpstp -> time0);
            }
         }
         
         if(ntpstp -> pmtindex1 > 0 ||
            (ntpstp -> pmtindex1 == 0 && ntpstp -> ph1.raw > 0.0))
         {
            if(!init_time)
            {
               data.max_time = ntpstp -> time1;
               data.min_time = ntpstp -> time1;
               init_time = true;
            }
            else
            {
               data.max_time = std::max(data.max_time, ntpstp -> time1);
               data.min_time = std::min(data.min_time, ntpstp -> time1);
            }
         }

         if(!init_time)
         {
            MSG("FillAlg", Msg::kWarning) << "Failed to find initial time" << endl;
            continue;
         }
      }
      else
      {
         MSG("FillAlg", Msg::kWarning) << "Unknown detector type" << std::endl;
      }
   } 

   data.nuplane = uplane_vec.size();
   data.nvplane = vplane_vec.size();

   if(!uplane_vec.empty())
   {
      data.beg_uplane = *std::min_element(uplane_vec.begin(), uplane_vec.end());
      data.end_uplane = *std::max_element(uplane_vec.begin(), uplane_vec.end());
   }

   if(!vplane_vec.empty())
   {
      data.beg_vplane = *std::min_element(vplane_vec.begin(), vplane_vec.end());
      data.end_vplane = *std::max_element(vplane_vec.begin(), vplane_vec.end());
   }

   if(detector == Detector::kNear)
   {
      short min_plane = -1, max_plane = -1, count = 0;

      for(vector<short>::const_iterator it = uplane_vec.begin(); it != uplane_vec.end(); ++it)
      {
         const short plane = *it;
         if(plane > 120)
         {
            break;
         }
         
         ++count;

         if(min_plane < 0 || min_plane > plane)
         {
            min_plane = plane;
         }
         if(max_plane < 0 || max_plane < plane)
         {
            max_plane = plane;
         }
      }
      
      for(vector<short>::const_iterator it = vplane_vec.begin(); it != vplane_vec.end(); ++it)
      {
         const short plane = *it;
         if(plane > 120)
         {
            break;
         }
         
         ++count;        

         if(min_plane < 0 || min_plane > plane)
         {
            min_plane = plane;
         }
         if(max_plane < 0 || max_plane < plane)
         {
            max_plane = plane;
         }
      }
      
      if(max_plane - min_plane > 0)
      {
         data.active_frac = double(count)/double(1 + max_plane - min_plane);
      }
   }
   else
   {
      short min_plane = data.beg_uplane;
      short max_plane = data.end_uplane;

      if(min_plane > data.beg_vplane)
      {
         min_plane = data.beg_vplane;
      }
      if(max_plane < data.end_vplane)
      {
         max_plane = data.end_vplane;
      }
      
      const double nactive = uplane_vec.size() + vplane_vec.size();
      if(max_plane - min_plane > 0)
      {
         data.active_frac = nactive/double(max_plane - min_plane);
      }
   }

   return data;
}

const Anp::Vertex Anp::FillBasic::Get	(	const NtpSRVertex &	vtx,
		bool	isfid_
	)			const [private]

Definition at line 169 of file FillBasic.cxx.

References NtpSRVertex::dcosu, NtpSRVertex::dcosv, NtpSRVertex::dcosz, NtpSRVertex::u, NtpSRVertex::v, and NtpSRVertex::z.

Referenced by Fill().

{
   //
   // This is really just a Vertex constructor, put it here to avoid NtpSRVertex dependency 
   //

   return Vertex(vtx.u,
                 vtx.v,
                 vtx.z,
                 vtx.dcosu,
                 vtx.dcosv,
                 vtx.dcosz,
                 isfid_);
}

const Basic Anp::FillBasic::GetBasic

(

)

const [inline]

Definition at line 72 of file FillBasic.h.

References fBasic.

Referenced by Anp::FillShower::Run(), Anp::FillTrack::Run(), and Anp::FillEvent::Run().

{ return fBasic; }

const Vertex Anp::FillBasic::GetBegVtx

(

)

const [inline]

Definition at line 73 of file FillBasic.h.

References fBegVtx.

Referenced by Anp::FillShower::Run(), Anp::FillTrack::Run(), and Anp::FillEvent::Run().

{ return fBegVtx; }

const Vertex Anp::FillBasic::GetEndVtx

(

void

)

const [inline]

Definition at line 74 of file FillBasic.h.

References fEndVtx.

Referenced by Anp::FillTrack::Run().

{ return fEndVtx; }

---

Member Data Documentation

Basic Anp::FillBasic::fBasic [private]

Definition at line 60 of file FillBasic.h.

Referenced by Fill(), and GetBasic().

Vertex Anp::FillBasic::fBegVtx [private]

Definition at line 61 of file FillBasic.h.

Referenced by Fill(), and GetBegVtx().

bool Anp::FillBasic::fCheck [private]

Definition at line 57 of file FillBasic.h.

Referenced by Config(), and Fill().

Vertex Anp::FillBasic::fEndVtx [private]

Definition at line 62 of file FillBasic.h.

Referenced by Fill(), and GetEndVtx().

double Anp::FillBasic::fTrackZOffset [private]

Definition at line 58 of file FillBasic.h.

Referenced by Config(), and Fill().

---

The documentation for this class was generated from the following files:

• FillBasic.h
• FillBasic.cxx