FillMuonId.h

Go to the documentation of this file.

#ifndef ANP_FILLMUONID_H
#define ANP_FILLMUONID_H

// $Id: FillMuonId.h,v 1.14 2008/11/11 23:06:36 rustem Exp $
//
// Fill variables for muon identification
//

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

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

class VldContext;

namespace Anp
{
   namespace Fill
   {
      struct AddCharge : public std::unary_function<double, void>
      {
         AddCharge(double value_);
         void operator() (double &x);
         const double value;
      };
   }

   class FillMuonId: public AlgSnarl, public AlgEvent
   {
   public:

      typedef std::vector<double> ChargeVec;
      typedef std::map<short, float> DataMap;

      struct PlaneData
      {
         PlaneData();

         short min_strip;
         short max_strip;

         ChargeVec track_vec;
         ChargeVec other_vec;    

         float tpos;
      };
      
      typedef std::map<short, PlaneData> PMap;
      typedef PMap::const_iterator PIter;

   public:

      FillMuonId();      
      virtual ~FillMuonId();
    
      bool Run(Record &record);

      bool Run(Event &event, const Record &record, const bool pass);

      void Config(const Registry &reg);

      const DataMap Run(const Track &track, const Record &record) const;

      void Fill(const PMap &pmap, DataMap &dmap, short base) const;

      const PMap Select(const Track &track, const Record &record, const std::string &opt = "") const;

      const PMap AddCharge(const PMap &pmap, double calor, double spect) const;      

   private:
      
      void Fill0(PIter beg, PIter end, DataMap &dmap) const;
      void Fill1(PIter beg, PIter end, DataMap &dmap, short base) const;
      void Fill2(PIter beg, PIter end, DataMap &dmap, short base) const;
      void Fill3(PIter beg, PIter end, DataMap &dmap, short base) const;
      void Fill4(PIter beg, PIter end, DataMap &dmap, short base) const;

      const VldContext GetVldc(const Header &header) const;

      float GetLPos(short plane, const PMap &pmap) const;
      float Get(const Strip &strip, const VldContext &vldc, float lpos, short option) const;

   private:

      short fStripWindow;
      short fMinViewNPlane;
      short fKeyBase;
      int   fKeyPass;

      double fTimeWindow;

      double fMinStripAdc;
      double fMinStripMip;

      double fADCShiftMuon;
      double fADCShiftMuonSpect;

      double fADCShiftOther;
      double fADCShiftOtherSpect;

      bool fIsStudy;   // fill additional variables?
      bool fFillDetec; // near only, use all detector strips or only calorimeter?
      bool fFillEvent; // add variable to track and event data vectors
      bool fFillShift; // fill adc shifted variables
      bool fFillShort; // only fill short track information
   };
   
   //
   // Inlined member functions
   //
   inline Fill::AddCharge::AddCharge(const double value_)
      :value(value_)
   {
   }
   inline void Fill::AddCharge::operator() (double &x)
   {
      x += value;
   }
}

#endif

---