AlignChamZNoBF.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: AlignChamZNoBF.h,v 1.1 2006/09/25 14:26:09 lebedev Exp $
//
// lebedev@physics.harvard.edu
////////////////////////////////////////////////////////////////////////

#ifndef TRK_ALIGN_NO_B_F_H
#define TRK_ALIGN_NO_B_F_H

#include "JobControl/JobCModuleFactory.h"
#include "Geometry/GChamGeo.h"

#include <string>
#include <vector>

//......................................................................
/// @author Andre Lebedev
/// This module can be used to fit for DC1-PWC6 misalignment in z
/// and rotation in xy-plane using field-off data. 
/// Track information is saved until the 
/// end of the run, and then TMinuit is used to minimize either the
/// sum of chi squared or the sum of square of correlation coefficients 
/// between residuals and dx/dz. The former method is more general
/// and will pick up rotations much better than the latter.
///
class AlignChamZNoBF : public JobCModule {
public:
  typedef struct track__ {
    float fWire[6][kNPlane]; ///< Wire number in the cluster
    float fWeight[6][kNPlane]; ///< Weight of the wire in the fit
  } TrkInfo_t;

  AlignChamZNoBF(const char* version);
  ~AlignChamZNoBF();

  JobCResult Ana(const EDMEventHandle& evt);
  
  void NewRun(int run, int subrun);
  
  void EndRun(int run, int subrun);

  void Update(const CfgConfig& c);

  double ComputeF();

private:
  double Fit(TrkInfo_t& ti, double resid[6][4], float* dydz = 0,
         float x[6][4] = 0, float y[6][4]  = 0);
  void MakeTree(const char* name);
  double ComputeChi2F();
  double ComputeCorrF();

private:
  std::vector<TrkInfo_t> fTrackList;

  int fDebug; // Controls level of debug printing
  std::string fSegFolder;
  int fMinNClust;
  std::vector<int> fFixCham;
  int fUseCorr;

};

#endif // TRK_ALIGN_NO_B_F_H

////////////////////////////////////////////////////////////////////////

---