AlignRotB.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: AlignRotB.h,v 1.1 2007/01/09 21:45:45 lebedev Exp $
//
// lebedev@physics.harvard.edu
////////////////////////////////////////////////////////////////////////

#ifndef TRK_ROTB_H
#define TRK_ROTB_H

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

#include <vector>
#include <string>

//......................................................................
/// @author Andre Lebedev
/// This module can be used to determine rotation of B-field components
/// around z-axis. The reason one may want to do this is that in the
/// event that Hall probes were not perfectly well aligned during
/// survey, By will bleed into Bx, which will cause an incorrect
/// kick in y. This module requires TrkCandBuilder with BC tracks
/// pulled into the fit.
///
class AlignRotB : 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;

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

  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();

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 fMinNClust;
  std::vector<int> fFixParam;
};

#endif // ALIGN_ROTB_H

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

---