FillSnarl.h

Go to the documentation of this file.

#ifndef ANP_FILLSNARL_H
#define ANP_FILLSNARL_H

// $Id: FillSnarl.h,v 1.2 2008/03/31 17:25:35 rustem Exp $
//
// Fill snarl data and summary of snarl strips
//

// C++
#include <algorithm>
#include <vector>

// Local
#include "PhysicsNtuple/AlgSnarl.h"

namespace Anp
{
   class FillSnarl: public AlgSnarl
   {
   public:
      
      FillSnarl();
      virtual ~FillSnarl();
      
      bool Run(Record &record);
      
      void Config(const Registry &reg);

   private:

      bool fData; // fill summary of all strips in snarl
      int  fBase; // base key for snarl data
      int  fJump; // starting value for pulse height threshold scan
      int  fStep; // step for pulse height threshold scan
      int  fStop; // maximum value for pulse height threshold scan

      unsigned int fOver; // maximum number of times strip can appear in same snarl
   };

   //
   // Utility function that counts equal strips
   //
   struct CountPlaneStrip : public std::unary_function<double, void>
   {
      public:
      
      CountPlaneStrip();
      CountPlaneStrip(short plane_, short strip_);
      
      void operator()(const Strip &strip);
      
      unsigned int GetCount() const;
      short GetPlane() const;
      short GetStrip() const;

      private:
      
      short plane;
      short strip;     
      unsigned int count;
   };

   bool operator==(const CountPlaneStrip &lhs, const CountPlaneStrip &rhs);

   //
   // Utility function that count strips below and above threshold
   //
   struct CountAboveBelow : public std::unary_function<double, void>
   {
      public:
      
      CountAboveBelow();
      CountAboveBelow(float threshold_, const std::vector<std::pair<short,short> > &ivec);
      
      void operator()(const Strip &strip);
      
      unsigned int GetNAbove() const;
      unsigned int GetNBelow() const;

      double GetMeanAbove() const;
      double GetMeanBelow() const;
      
      private:
            
      std::vector<std::pair<short,short> > ignorevec;

      float threshold;
      unsigned int nabove;
      unsigned int nbelow;

      double sumabove;
      double sumbelow;
   };

}

#endif

---