FillShortVar.h

Go to the documentation of this file.

#ifndef ANP_FILLSHORTVAR_H
#define ANP_FILLSHORTVAR_H


//
// Fill variables for low energy muon identification
//

// C++
#include <map>

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

class VldContext;

namespace Anp
{
  
   class FillShortVar: public AlgSnarl
   {
   public:

     FillShortVar();      
     virtual ~FillShortVar();
     
     bool Run(Record &record);
     
     void Config(const Registry &reg);
     
     void End(const DataBlock &);

   private:
     
     typedef StripLink* StripLinkIter;
     typedef std::map<short,float> DataMap;
     typedef std::map<unsigned short, std::vector<Anp::StripIter> > PlaneMap;
     typedef std::map<unsigned short, std::vector<StripLinkIter> > PlaneMapList;
     
   private:
      
     bool Study(Track &track, const Record &record,DataMap& dmap, string StripMin);
     void Fill(DataMap& dmapout, DataMap dmapin, int base);
     bool GetTrackMap(Track &track, const Record &record, PlaneMap& track_map,std::map<int,float>& other_map,string StripMin) const;
     
     std::vector<Anp::StripIter> Get(const Track &track,
                                     unsigned short plane,
                                     const Record &record, std::map<int,float> other_map, string StripMin) const;
  
     bool GetScatter(PlaneMap track_map, float& pear, float& dpear, int keyval) const;
     double BVD_BestLine(PlaneMap& map) const ;
     double GetProbability(StripLink* cur, StripLink* cmp) const;
     float GetWPos(const short plane, std::map<int,float> pmap) const;
     const VldContext GetVldc(const Header &header) const;
     float TrackStripEnergy(StripIter si, short index, string StripMin) const;
     float OtherStripEnergy(StripIter sp, const VldContext &vldc,
                            const float lpos, string StripMin) const;
     bool IsGoodTrack(TrackIter itrack, Record& record);
     float GetStripChargeCut(string StripMin) const;
     std::vector<Anp::StripIter> FirstPassOtherStrip(const std::vector<StripIter> instrip, const std::vector<StripIter> track_strip, double minT, double maxT) const;
   private:
     int fMinNPlane;
     int fMaxNPlane;
     int fEndNPlane;
     int fKeyBase;     
     double fTimeWindow;
     int fStripWindow;
     double fMinStripADC;
     double fMinStripMIP;
     bool fFillEvent;
     bool fDebug;
     bool fErase;
     
     bool fqeltrack;
     double ftruelowE;
     double frecolowE;
     int fselectsign;
      //
      // Computer variables
      //
      int fNTrackAll;
      int fNTrackAdd;
   };  
}

#endif

---