AlignChamZ.h

Go to the documentation of this file.

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

#ifndef TRK_ALIGN_Z_H
#define TRK_ALIGN_Z_H

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

#include <vector>
#include <string>

//......................................................................
/// @author Andre Lebedev
/// This module can be used to determine alignment in z for DC1-PWC6
/// using field-on data. In principle, this should also work for
/// field-off data. Reconstructed TrkCand's are picked up, and all
/// single-wire clusters are used, except for those belonging to
/// BC's. At the end of the run, TMinuit is used to minimize either
/// the sum of all chi squared (faster and I think more reliable)
/// or the sum of the squares of correlation coefficients between
/// residuals and dx/dz.
///
class AlignChamZ : public JobCModule {
public:
  typedef struct track__ {
    float fWire[kNCham][kNPlane]; ///< Wire number in the cluster
    float fWeight[kNCham][kNPlane]; ///< Weight of the wire in the fit
  } TrkInfo_t;

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

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

  void Update(const CfgConfig& c);

  double ComputeF();

private:
  double Fit(const TrkInfo_t& ti, double resid[kNCham][kNPlane],
         double dxdz[kNCham]);
  void MakeTree();
  double ComputeCorrF();
  double ComputeChi2F();

private:
  std::vector<TrkInfo_t> fTrackList;
  double fLambdaX[kNCham][kNPlane];
  double fLambdaY[kNCham][kNPlane];
  double fZ0;

  int fDebug; // Controls level of debug printing
  std::string fFolder;
  int fLoadAlignment; 
  int fMinNClust;
  std::vector<int> fFixParam;
  int fUseCorr;

};

#endif // TRK_ALIGN_Z_H

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

---