AlignBC.h

Go to the documentation of this file.

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

#ifndef TRK_ALIGN_BC
#define TRK_ALIGN_BC

#include "JobControl/JobCModuleFactory.h"
#include "Geometry/GChamGeo.h"
#include <TObject.h>
#include <string>
#include <vector>

//......................................................................
class BCAlignNt : public TObject {
public:
  BCAlignNt() {;}

  int   nclust;
  float chi2;
  float par[4];
  float u[3][4];
  float ufit[3][4];
  float uperp[3][4];
  float z1, z2;
  float chi2rm[3][4];
  float resid[3][4];

  ClassDef(BCAlignNt, 1); // Ntuple for BC alignment summary
};

//......................................................................
/// Module which takes BC track segments and determines alignment
/// constants for all planes. 
/// One needs 6 parameters to determine a coordinate system:
/// three 3D points minus 2 equations which normalizes the (1, 0, 0)
/// and (0, 1, 0) vectors, minus 1 equation which ensures their 
/// orthogonality.
/// If length scale is set by wire spacing, then fixed parameters can
/// be chosen as z of BC3, wide-angle planes of BC3 and BC1, and
/// orientation of BC3 in xy-plane.
///
/// This module is currently meant to do "rough" alignment, i.e. w/o
/// using drift time information. Hence, it is unlikely that
/// XZ or YZ rotations or z-positions will be fit for correctly. 
/// Nevertheless, ignoring rotations in xz and yz planes, we have
/// 6 parameters per chamber (4 planes, z, xy-angle), of which we
/// must fix at least 6 in order for minimization to be meaningful.
///
/// @author Andre Lebedev
///
class AlignBC : public JobCModule {
public:
  typedef struct {
    float  fWire[kNBC][kNPlane]; // Wire can be a cluster, e.g. 80.5
    float  fT[kNBC][kNPlane];
    float  fWeight[kNBC][kNPlane];
  } TrkInfo_t;

  AlignBC(const char* version);
  ~AlignBC() {;}

  JobCResult Ana(const EDMEventHandle& evt);

  void NewRun(int run, int subrun);
  void EndRun(int run, int subrun);
  void Update(const CfgConfig& c);
  
private:
  //static double ComputeChi2(TrkInfo_t& ti);
  bool Fit(const TrkInfo_t* trk, BCAlignNt* nt, bool resolve = false);
  bool PreenTrackList(double res[3][4]);
  bool RoughAlign();
  bool FineAlign();
  void FillNt(const char* name, bool resolve);
  void WriteToDB(int run);

private:
  static AlignBC* fsModule;
  static double fZ1, fZ2; ///< Z locations that determine track params
  std::vector<TrkInfo_t*> fTrkInfo;

  int fDoFineAlign;
  int fMaxAlignAttempt;
  int fMinNClust; ///< Minimum number of clusters per track
  int fMinNTrack;
  int fMaxNTrackEvt;
  std::string fTrkFolder; // Folder where tracks are found
  double fOkRoughShift;
  double fOkFineShift;
  double fDriftVel;
  float  fMinTime;
  float  fMaxTime;
  int    fFillNt;
};

#endif // TRK_ALIGN_BC

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

---