AlignTPCDrift2.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: AlignTPCDrift2.h,v 1.1 2007/03/06 15:49:49 lebedev Exp $
//
// lebedev@physics.harvard.edu
////////////////////////////////////////////////////////////////////////

#ifndef ALIGN_TPC_DRIFT2_H
#define ALIGN_TPC_DRIFT2_H

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

class SwimMIPP;

//......................................................................
/// @author lebedev@physics.harvard.edu
/// Class which picks up global tracks and minimizes the sum of chi
/// squared by scaling drift velocity vector components
///
class AlignTPCDrift2 : public JobCModule {
public:
  typedef struct {
    int   fNTPCHit;
    float fPadX[128];
    float fPadZ[128];
    float fDriftT[128];

    int   fNChamHit;
    float fChamZ[36];
    float fU[36];
    float fW[36];
    float fCosA[36];
    float fSinA[36];
    double fPar[5];
    double fZ0;
    SwimMIPP* fSwim;
  } TrkInfo_t;

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

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

  void Update(const CfgConfig& c);

  static void MinFCN(int& /*npar*/, double* /*gin*/, double& f,
            double* par, int /*iflag*/);

private:
  double ComputeTrkChi2(TrkInfo_t& trk);

private:
  static AlignTPCDrift2* gModule;
  
  std::vector<TrkInfo_t> fTrkInfo;

  int fDebug; // Controls level of debug printing
  std::vector<int> fFixParam; 
  std::string fFolder;
  int fMinTrkVtx;
  int fRequireBeamTrk;
  int fMinTPCHits;
  int fMaxTPCHits;
  int fMinChamHits;
  float fMinP;
  float fMinGoF;
  float fTPCSigmaX;
  float fTPCSigmaY;
  float fMinTrkTime;
  float fMaxTrkTime;
  
};

#endif // ALIGN_TPC_DRIFT2_H

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

---