// Silicon r-phi intrinsic resolution analysis driver script
{
  TChain* chain = new TChain("STNTUPLE");

  chain.Add("../4.5.0int4/stnmaker_*_b_svx.root");

  TStnAna x(chain);

  gSystem->CompileMacro("TRPhiWt.cc","k");

  TRPhiWt *wt = new TRPhiWt("RPhiWt","wt","wt2.root");
  int nType=1;
  int nmin[1]={2};
  int nmax[1]={2};
  wt->SetNStripPerClus(nType,nmin,nmax);
  //  wt->SetAnchor(1); // L1 & L3 forced to cluster type 1 (2-strip)
  wt->LoadFile(false);
  wt->PtCut(15.0);
  wt->UseAlignment();

  x.AddModule(wt);

  x.Run();
  
}