FillShortEvent.h

Go to the documentation of this file.

#ifndef ANP_FILLSHORTEVENT_H
#define ANP_FILLSHORTEVENT_H

// $Id: FillShortEvent.h,v 1.3 2009/10/15 21:31:51 jyuko Exp $
//
//

// C/C++
#include <string>

// Local
#include "PhysicsNtuple/AlgEvent.h"
#include "PhysicsNtuple/AlgSnarl.h"
#include "PhysicsNtuple/Plot/Data.h"
#include "PhysicsNtuple/ShortMu/StripLink.h"

class TDirectory;
class TH1;
class TH2;

namespace Anp
{
     class FillShortEvent: public AlgSnarl
     {
     private:
     
       typedef StripLink* StripLinkPtr;

       typedef std::map<unsigned short, std::vector<Anp::StripIter> > PlaneMap;
       typedef std::map<unsigned short, std::vector<StripLinkPtr> > PlaneMapList;
       
   public:
      FillShortEvent();
      virtual ~FillShortEvent();
      bool Run(Record &record);
      void Config(const Registry& reg);
      bool Init(const Header& );

      void End( const DataBlock &);
      bool Explore( Event& ievent, const Record& record);
      
   private:
      bool GetPlaneStrips(PlaneMap& plane_mapu, PlaneMap& plane_mapv, const Event ievent, const Record& record);
      //      void GetEventQuantities(PlaneMap plane_m, const Vertex vtx,double& mean, double& offset, double& angl, double & pearson);
      bool GetEventQuantities(PlaneMap plane_m, std::map<int,double>* vect);
      double BVD_BestLine(PlaneMap& map);
      //      void GetMeanLine(PlaneMap map, double& deriv, double& cons, double& sig_deriv, double& sig_cons, double& gNHits, bool best);
      // void GetWellSeparated(PlaneMap plane_m,PlaneMap& result,  double num);
      double GetProbability(StripLink* cur, StripLink* cmp);
      bool   FillLink(StripLink& striplink, double ins,double ext, double ddx, StripIter strip, StripLink* link);
      bool   FillLink(StripLink& striplink, double ins,double ext, double ddx, StripIter strip);

      double fRecoHighE; // highest energy event 
      double fNEndPlanes; // number of planes you iterate back
      double fMinADCHit; // minimum ADC hit 
      double fNPlaneSkip; // number of planes you are allowed to skip
      double fTimeWindow; // timing window
      int fStripWindow;
      int fKeyBase;
      int fMinDLength;
      double fCutOff;
      int fKeyPass;
      int fKeyPass2;
      std::string fDirName;
      TDirectory *fDir;
      bool fDebug;
      bool fErase;  
      bool   fPlot;
      int    fNMiss;
      int    fTooShort;
      int    fBVDmiss;

   };
   
}
   
#endif

---