AlProject56.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: AlProject56.cxx,v 1.2 2006/11/07 21:29:33 lebedev Exp $
//
// lebedev@physics.harvard.edu
////////////////////////////////////////////////////////////////////////

#include "Alignment/AlProject56.h"

#include "Bfield/bfield.h"
#include "Config/CfgConfig.h"
#include "Config/CfgParam.h"
#include "EventDataModel/EDMEventHandle.h"
#include "TrkRBase/TrackSeg.h"
#include "TrkRBase/Wire.h"

MODULE_DECL(AlProject56);
ClassImp(AlProjectNt);

using namespace std;

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

AlProject56::AlProject56(const char* version) :
  JobCModule("AlProject56")
{
  // Register module to be updated with a given config version
  SetCfgVersion(version);
  // Make sure that configuration has been loaded
  CheckInit();

  fNt = new AlProjectNt;
  fTree = new TTree("ptree", "Segment projection tree");
  fTree->Branch("fNt", "AlProjectNt", &fNt);
  
}

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

AlProject56::~AlProject56()
{
  fTree->Write();
  delete fTree;
}

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

JobCResult AlProject56::Ana(const EDMEventHandle& evt)
{
  vector<const TrackSeg*> seg;
  vector<const Wire*> wires;

  try { evt.Reco().Get(fSegFolder.c_str(), seg); }
  catch (...) { }
  if (seg.empty()) return kFailed;

  double x[3];

  for (int cham = 0; cham < kNCham; ++cham) {
    GChamGeo& geo = GMIPPGeo::Instance().Cham(cham);
    
    wires.clear();
    try { evt.Cal().Get(geo.Name(), wires); }
    catch (...) { continue; }
    if (wires.empty()) continue;

    for (unsigned int i = 0; i < seg.size(); ++i) {
      const TrackSeg* trk = seg[i];
      // Skip tracks which have this chamber in the fit
      if (trk->Point(cham)) continue; 
      
      // Predict u in every plane
      for (int j = 0; j < kNPlane; ++j) {
    trk->GetXY(geo.Z(j + 1), x, 0);
    double u;
    geo.XYToU(j + 1, x[0], x[1], u);
    fNt->u[j] = u;
    geo.XYToU(j + 1, trk->dXdZ(), trk->dYdZ(), u);
    fNt->dudz[j] = u;
      }
      for (unsigned int j = 0; j < wires.size(); ++j) {
    const Wire* w = wires[j];

    fNt->cham = cham;
    fNt->plane = w->Plane();;
    fNt->wire = w->WireHit();
    fTree->Fill();
      }
    }
  }
  
  return kPassed;
}

//......................................................................
  
void AlProject56::NewRun(int run, int subrun)
{
  GMIPPGeo::Instance(run);
  config_bfield(run, subrun);
}

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

void AlProject56::Update(const CfgConfig& c)
{
  c("Debug").Get(fDebug);
  c("SegFolder").Get(fSegFolder);

  // Change the flag to signify that our configuration has been read
  fIsInit = true;
}

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

---