AlignChamW0.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: AlignChamW0.h,v 1.2 2006/12/25 06:29:31 lebedev Exp $
// 
// lebedev@physics.harvard.edu
////////////////////////////////////////////////////////////////////////

#ifndef TRK_CHAM_ALIGN_H
#define TRK_CHAM_ALIGN_H

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

class TrackSeg;
class WireClust;
class TH1F;
class TH2F;
class TTree;

//......................................................................

class AlResidNt : public TObject {
public:
  AlResidNt() {;}

  int nclust;
  float u[kNCham][kNPlane];  // U corresponding to actual wire
  float ufit[kNCham][kNPlane]; // U from the overal fit (all chambers in)
  float uperp[kNCham][kNPlane]; // Measure along the wire (all chambers in)
  float resid[kNCham][kNPlane]; // Residual when chamber is not in the fit
  float poq;  // Momentum with all chambers in
  float chi2; // Chi squared, all chambers in
  float poqrm[kNCham][kNPlane]; // Momentum with chamber taken out of the fit
  float chi2rm[kNCham][kNPlane]; // Chi squared, chamber taken out of the fit

  ClassDef(AlResidNt, 2); // Ntuple for residual analysis
};

//......................................................................
/// Module which calculates chamber to chamber alignment. At least
/// 2 chambers should be fixed in order for result to make sense.
/// Result can then be rota
class AlignChamW0 : 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
    float fPoQ;
    float fChi2;
  } TrkInfo_t;

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

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

private:
  bool AdjustBField();
  double Fit(TrkInfo_t& trk, float* p = 0);
  void Fit(AlignChamW0::TrkInfo_t& trk, AlResidNt* nt, int minCham = 0);
  void AnaResid(int i);
  int  PreenTrackList(double res[kNCham][kNPlane], bool withBC);
  bool RoughAlign(bool withBC);

private:
  bool fFieldOn;
  int fNClust[kNCham][kNPlane]; ///< Total number of clusters in each chamber
  double fWire0[kNCham][kNPlane]; ///< Initial wire zero for all planes

  static AlignChamW0* gsModule;
  static int fsCham; ///< Chamber being aligned internally

  std::vector<TrkInfo_t> fTrackList;
  int fNTrackWithBC;
  double fZ0; ///< Z0 -- reference needed to calculate lambda's

  std::string fFolder; ///< Folder where to look for tracks
  double fOkShift;
  int fLoadAlignment;
  int fMaxAlignAttempt; ///< Maximum number of attempts before alignment fails
  int fMinNClust; ///< Minimum number of clusters in Cham 1-6
  int fMinTracks; ///< Minimum number of tracks required to do alignment
  int fMinTracksBC;
  int fUseBC;
  int fAdjustBField;
  double fMinMom;
  double fMaxMom;
};

#endif // TRK_CHAM_ALIGN_H

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

---