DFDCSegment.h

Go to the documentation of this file.

//***********************************************************************
// DFDCSegment.h : definition for a track segment built from pseudopoints
//***********************************************************************

#ifndef DFDCSEGMENT_H
#define DFDCSEGMENT_H

#include <JANA/JObject.h>
using namespace jana;

#include "DFDCHit.h"
#include "DFDCWire.h"
#include "DFDCPseudo.h"

#include <DMatrix.h>
#include <sstream>

// Residuals and arc length
struct fdc_track_t{
  unsigned int hit_id;
  double sign; // Sign of left-right ambiguity resolution for this hit
  double dx,dy; //residuals
  double chi2; // chi2 contribution of this point
  double s; // path length
};


///
/// class DFDCSegment: definition for a track segment in the FDC
/// 
class DFDCSegment : public JObject {
   public :
      JOBJECT_PUBLIC(DFDCSegment);        /// DANA identifier
      
      /// 
      /// DFDCSegment::DFDCSegment():
      /// Default constructor
      ///
      DFDCSegment(){}
      DFDCSegment(const DFDCSegment &segment){
        this->xc=segment.xc;
        this->yc=segment.yc;
        this->rc=segment.rc;
        this->tanl=segment.tanl;
        this->q=segment.q;
        this->D=segment.D;
        this->z_vertex=segment.z_vertex;
        this->phi0=segment.phi0;
        this->Phi1=segment.Phi1;
        this->package=segment.package;
        this->chisq=segment.chisq;
        this->Ndof=segment.Ndof;
      }

      double chisq;
      int Ndof;

      // circle parameters
      double xc,yc,rc;  
      // azimuthal angle of the intersection of the helical path to 
      // the most downstream plane in a given package containing a
      // hit
      double Phi1;                     
      // "vertex" z position and the phi angle there
      double z_vertex,phi0;
      // tangent of the dip angle
      double tanl;
      // charge (+/-1)
      double q;
      // Distance of closest approach to the beam line
      double D;
      // Package number
      unsigned int package;
      
      // List of pseudopoints belonging to this track segment
      vector<const DFDCPseudo *>hits;  

      void toStrings(vector<pair<string,string> > &items)const{
         AddString(items, "xc", "%3.2f", xc);
         AddString(items, "yc", "%3.2f", yc);
         AddString(items, "rc", "%3.2f", rc);
         AddString(items, "phi0(deg)","%3.2f",57.296*phi0);
         AddString(items, "Phi1(rad)", "%3.2f", Phi1);
         AddString(items, "D","%3.2f",D);
         AddString(items, "q","%3.1f",q);
         AddString(items, "Nhits", "%d", hits.size());
      }
};

#endif //DFDCSEGMENT_H