StRoot: macros/obsolete/dst.C Source File

00001 // $Id: dst.C,v 1.5 1999/05/21 15:33:58 kathy Exp $
00002 // $Log: dst.C,v $
00003 // Revision 1.5  1999/05/21 15:33:58  kathy
00004 // made sure Log & Id are in each file and also put in standard comment line with name of owner
00005 //
00006 // Revision 1.4  1999/05/06 03:21:28  fisyak
00007 // synchronize FTPC and TPC slow/fast
00008 //
00009 // Revision 1.3  1999/01/03 20:59:55  fisyak
00010 // Remove St_geom_Maker
00011 //
00012 // Revision 1.2  1998/12/29 19:38:22  fine
00013 // STAR_shapes: test of the brand new St_NodeView class has been added
00014 //
00015 // Revision 1.1  1998/11/01 16:42:29  fisyak
00016 // dst analysis
00017 //=======================================================================
00018 // owner: Yuri Fisyak
00019 // what it does: 
00020 //=======================================================================
00021 //
00022 TBrowser *b = 0;
00023 class StChain;
00024 StChain  *chain=0;
00025 void Load(){
00026   gSystem->Load("St_base");
00027   gSystem->Load("StChain");
00028   gSystem->Load("xdf2root");
00029   gSystem->Load("St_Tables");
00030   gSystem->Load("St_xdfin_Maker");
00031   gSystem->Load("St_TLA_Maker");
00032 }
00033 void dst(const Int_t Nevents=1,
00034 const Char_t *fileinp = "/disk1/star/test/tfs/psc0065_01_40evts_dst.xdf")
00035 {
00036 // Dynamically link some shared libs
00037   if (gClassTable->GetID("StChain") < 0) Load();
00038   St_XDFFile   *xdf_in   = 0;
00039   if (fileinp)  xdf_in   = new St_XDFFile(fileinp,"r");
00040 //TFile      *root_tree= new TFile("auau_central_hijing.tree.root","RECREATE");
00041 // Create the main chain object
00042   StChain *chain = new StChain("dst_ana");
00043   //  St_params_Maker *params = new St_params_Maker("params","run/params");
00044   //  St_geant_Maker     *geant = new St_geant_Maker("geant","run/geant/Run");
00045   if (xdf_in) {
00046     St_xdfin_Maker *xdfin = new St_xdfin_Maker("xdfin");
00047     xdfin->Init_Done(kTRUE);
00048     chain->SetInputXDFile(xdf_in);
00049   }
00050   chain.PrintInfo();
00051 // Init the mai chain and all its makers
00052   int iInit = chain.Init();
00053   if (iInit) chain.Fatal(iInit,"on init");
00054   gBenchmark->Start("dst");
00055   Int_t i=0;
00056   for (Int_t i =1; i <= Nevents; i++){
00057     if (chain.Make(i)) break;
00058     if (i != Nevents) chain.Clear();
00059     printf ("===========================================\n");
00060     printf ("=========================================== Done with Event no. %d\n",i);
00061     printf ("===========================================\n");
00062   }
00063   if (Nevents > 1) {
00064     chain.Finish();
00065     delete xdf_in;
00066     gBenchmark->Print("dst");
00067   }
00068   else b = new TBrowser("dst_ana");
00069 }