StRoot: StEEmcPool/FgtSandbox1/macros/bfc.C Source File

00001 
00002 //                                                                      //
00003 // Macro for running chain with different inputs                        //
00004 // owner:  Yuri Fisyak                                                  //
00005 //                                                                      //
00006 // $Id: bfc.C,v 1.3 2011/04/08 22:35:02 balewski Exp $
00008 class StBFChain;        
00009 class StMessMgr;
00010 #if !defined(__CINT__) || defined(__MAKECINT__)
00011 #include "Stiostream.h"
00012 #include "TSystem.h"
00013 #include "TClassTable.h"
00014 #include "TApplication.h"
00015 #include "TInterpreter.h"
00016 #include "StBFChain.h"
00017 #include "StMessMgr.h"
00018 #else
00019 #endif
00020 #define UseLogger
00021 StBFChain    *chain=0; 
00022 //_____________________________________________________________________
00023 //_________________ Prototypes _______________________________________________
00024 void Usage();
00025 void Load(const Char_t *options="");
00026 TString defChain("y2005h,Test.default.ITTF");
00027 void bfc(Int_t First, Int_t Last,const Char_t *Chain = defChain + ",Display",
00028          const Char_t *infile=0, const Char_t *outfile=0, const Char_t *TreeFile=0);
00029 //       const Char_t *Chain="gstar,y2005h,MakeEvent,trs,sss,svt,ssd,fss,bbcSim,emcY2,tpcI,fcf,ftpc,SvtCL,svtDb,ssdDb,svtIT,ssdIT,ITTF,genvtx,Idst,event,analysis,EventQA,tags,Tree,EvOut,McEvOut,GeantOut,IdTruth,miniMcMk,StarMagField,FieldOn,McAna,Display",//,,NoSimuDb, display, //McQa, 
00030 void bfc(Int_t Last, const Char_t *Chain = defChain,
00031          const Char_t *infile=0, const Char_t *outfile=0, const Char_t *TreeFile=0);
00032          //      const Char_t *Chain="gstar,y2005h,tpcDb,trs,tpc,Physics,Cdst,Kalman,tags,Tree,EvOut,McEvOut,IdTruth,miniMcMk,StarMagField,FieldOn,McAna", // McQA
00033 //_____________________________________________________________________
00034 void Load(const Char_t *options){
00035   cout << "Load system libraries" << endl;
00036   if ( gClassTable->GetID("TGiant3") < 0) { // ! root4star
00037     cout << endl << "Load ";
00038     if (!TString(options).Contains("nodefault",TString::kIgnoreCase) || 
00039          TString(options).Contains("pgf77",TString::kIgnoreCase)) {
00040       const Char_t *pgf77 = "libpgf77VMC";
00041       if (gSystem->DynamicPathName(pgf77,kTRUE) ) {
00042         gSystem->Load(pgf77); cout << " " << pgf77 << " + ";
00043       }
00044     }
00045     if (!TString(options).Contains("nodefault",TString::kIgnoreCase) || 
00046         TString(options).Contains("cern",TString::kIgnoreCase)) {
00047       gSystem->Load("libminicern"); cout << "libminicern";
00048     }
00049     if (!TString(options).Contains("nodefault",TString::kIgnoreCase) || 
00050         TString(options).Contains("mysql",TString::kIgnoreCase)) {
00051       Char_t *mysql = "libmysqlclient";
00052       //
00053       // ATTENTION: The below will FAIL for 64 bits systems (JL 2009/10/22)
00054       //
00055       Char_t *libs[]  = {"", "/usr/mysql/lib/","/usr/lib/", 0}; // "$ROOTSYS/mysql-4.1.20/lib/",
00056       //Char_t *libs[]  = {"/usr/lib/", 0};
00057       Int_t i = 0;
00058       while ((libs[i])) {
00059         TString lib(libs[i]);
00060         lib += mysql;
00061         lib = gSystem->ExpandPathName(lib.Data());
00062         if (gSystem->DynamicPathName(lib,kTRUE)) {
00063           gSystem->Load(lib.Data()); cout << " + " << lib.Data() << endl;
00064           break;
00065         }
00066         i++;
00067       }
00068     }
00069   }
00070   //  if (gClassTable->GetID("TMatrix") < 0) gSystem->Load("StarRoot");// moved to rootlogon.C  StMemStat::PrintMem("load StarRoot");
00071 #ifdef UseLogger
00072   // Look up for the logger option
00073   Bool_t needLogger  = kFALSE;
00074   if (!TString(options).Contains("-logger",TString::kIgnoreCase)) {
00075     needLogger = gSystem->Load("liblog4cxx.so") <= 0;              //  StMemStat::PrintMem("load log4cxx");
00076   }
00077 #endif
00078   gSystem->Load("libSt_base");                                        //  StMemStat::PrintMem("load St_base");
00079 #ifdef UseLogger
00080   if (needLogger) {
00081     gSystem->Load("libStStarLogger.so");
00082     gROOT->ProcessLine("StLoggerManager::StarLoggerInit();");      //  StMemStat::PrintMem("load StStarLogger");
00083   }
00084 #endif
00085   gSystem->Load("libHtml");
00086   gSystem->Load("libStChain");                                        //  StMemStat::PrintMem("load StChain");
00087   gSystem->Load("libStUtilities");                                    //  StMemStat::PrintMem("load StUtilities");
00088   gSystem->Load("libStBFChain");                                      //  StMemStat::PrintMem("load StBFChain");
00089   gSystem->Load("libStChallenger");                                   //  StMemStat::PrintMem("load StChallenger");
00090 }
00091 //_____________________________________________________________________
00092 void bfc(Int_t First, Int_t Last,
00093          const Char_t *Chain,
00094          const Char_t *infile,
00095          const Char_t *outfile,
00096          const Char_t *TreeFile)
00097 { // Chain variable define the chain configuration 
00098   // All symbols are significant (regardless of case)
00099   // "-" sign before requiest means that this option is disallowed
00100   // Chain = "gstar" run GEANT on flight with 10 muons in range |eta| < 1 amd pT = 1GeV/c (default)
00101   // Dynamically link some shared libs
00102   if (gClassTable->GetID("StBFChain") < 0) Load(Chain);
00103   chain = new StBFChain(); cout << "Create chain " << chain->GetName() << endl;
00104   TString tChain(Chain);
00105   chain->cd();
00106   chain->SetDebug(1);
00107   if (Last < -3) return;
00108   chain->SetFlags(Chain);
00109   if (tChain == "" || ! tChain.CompareTo("ittf",TString::kIgnoreCase)) Usage();
00110   chain->Set_IO_Files(infile,outfile);
00111   if (TreeFile) chain->SetTFile(new TFile(TreeFile,"RECREATE"));
00112   gMessMgr->QAInfo() << Form("Process [First=%6i/Last=%6i/Total=%6i] Events",First,Last,Last-First+1) << endm;
00113   if (Last < -2) return;
00114   if (chain->Load() > kStOk) {
00115     gMessMgr->Error() << "Problems with loading of shared library(ies)" << endm;
00116     gSystem->Exit(1);
00117   }
00118   if (Last < -1) return;
00119   if (chain->Instantiate() > kStOk)  { 
00120     gMessMgr->Error() << "Problems with instantiation of Maker(s)" << endm;
00121     gSystem->Exit(1);
00122   }
00123   if (Last < 0) return;
00124   StMaker *dbMk = chain->GetMaker("db");
00125   if (dbMk) dbMk->SetDebug(1);
00126 #if 0
00127   // Insert your maker before "tpc_hits"
00128   Char_t *myMaker = "St_TLA_Maker";
00129   if (gClassTable->GetID(myMaker) < 0) {
00130           gSystem->Load(myMaker);//  TString ts("load "; ts+=myMaker; StMemStat::PrintMem(ts.Data());
00131   }
00132   StMaker *myMk = chain->GetMaker(myMaker);
00133   if (myMk) delete myMk;
00134   myMk = chain->New(myMaker,"before");
00135   if (myMk) {
00136     Char_t *before = "tpc_hits";
00137     StMaker *tclmk = chain->GetMaker(before);
00138     if (tclmk) chain->AddBefore(before,myMk);
00139   }
00140   // Insert your maker after "tpc_hits"
00141   myMk = chain->New(myMaker,"after");
00142   if (myMk) {
00143     Char_t *after = "tpc_hits";
00144     StMaker *tclmk = chain->GetMaker(after);
00145     if (tclmk) chain->AddAfter(after,myMk);
00146   }
00147   // this block is meant as an example ONLY
00148   // The default values are set in StRoot/StPass0CalibMaker/ StTpcT0Maker 
00149   // constructor and are suitable for production. You can change it here
00150   // for test purposes.
00151   if (chain->GetOption("TpcT0")) {
00152     StTpcT0Maker *t0mk = (StTpcT0Maker *) chain->GetMaker("TpcT0");
00153     if (t0mk) t0mk->SetDesiredEntries(10);
00154   }
00155 #endif
00156   {
00157     TDatime t;
00158     gMessMgr->QAInfo() << Form("Run is started at Date/Time %i/%i",t.GetDate(),t.GetTime()) << endm;
00159   }
00160   gMessMgr->QAInfo() << Form("Run on %s in %s",gSystem->HostName(),gSystem->WorkingDirectory()) << endm;
00161   gMessMgr->QAInfo() << Form("with %s", chain->GetCVS()) << endm;
00162   // Init the chain and all its makers
00163   TAttr::SetDebug(0);
00164   chain->SetAttr(".Privilege",0,"*"                );     //All  makers are NOT priviliged
00165   chain->SetAttr(".Privilege",1,"StIOInterFace::*" );     //All IO makers are priviliged
00166   chain->SetAttr(".Privilege",1,"St_geant_Maker::*");     //It is also IO maker
00167   chain->SetAttr(".Privilege",1,"StTpcDbMaker::*");       //It is also TpcDb maker to catch trips
00168   chain->SetAttr(".Privilege",1,"*::tpc_hits"); //May be allowed to act upon excessive events
00169   chain->SetAttr(".Privilege",1,"*::tpx_hits"); //May be allowed to act upon excessive events
00170   chain->SetAttr(".Privilege",1,"StTpcHitMover::*"); //May be allowed to act upon corrupt events
00171   chain->SetAttr(".Privilege",1,"*::tpcChain"); //May pass on messages from sub-makers
00172   chain->SetAttr(".Privilege",1,"StTriggerDataMaker::*"); //TriggerData could reject event based on corrupt triggers
00173 
00174   // add FGT stuff - jan
00175   TString &fnameOut1=chain->GetFileOut();
00176   cout<<"AAA"<<fnameOut1<<endl;
00177   TString fnameOut2=fnameOut1;
00178   fnameOut2.ReplaceAll(".root",".fgt");
00179   //fnameOut2+=".fgt";
00180   cout<<"FGT output core="<<fnameOut2<<endl;
00181 
00182   gSystem->Load("libStEEmcPoolFgtSandbox1");
00183   StFgtSlowSimuMaker  *myMkSM=new StFgtSlowSimuMaker;
00184   StFgtClustFindMaker *myMkCL=new StFgtClustFindMaker;
00185   StFgtClustEvalMaker *myMkEV=new StFgtClustEvalMaker;
00186   // ... position FGT maker in the chain
00187   Char_t *after = "geant"; StMaker *xMk = chain->GetMaker(after);  assert(xMk);
00188   // note , the lines should be in reversed order
00189   chain->AddAfter(after,myMkEV);  //position #3
00190   chain->AddAfter(after,myMkCL); //position #2
00191   chain->AddAfter(after,myMkSM); //position #1
00192 
00193   //...... auxiliary initialization
00194   HList=new  TObjArray;
00195   myMkSM->setHList(HList);  myMkCL->setHList(HList);  myMkEV->setHList(HList);
00196   myMkSM->initFrankModel("StRoot/StEEmcPool/FgtSandbox1/macros/BichselELossProbHighBG.dat");
00197   // myMkSM->forcePerpTracks(); //only  for testing,all tracks will go perp through the GEM gas volume
00198 //  myMkSM->useOnlyDisk(1); //only  for testing, disk # 1...6
00199   myMkSM->useOnlyDisk(0); //All disks
00200   myMkSM->setRadStripGain(1.00,0.03); //relative gain variation: mean=1 & sigma=0.03  
00201   myMkSM->setPhiStripGain(1.00,0.03); //relative gain variation: mean=1.5 & sigma=0.03  
00202   myMkSM->setStripThresh(1.); // a.u. 
00203   myMkSM->setHexGemLatice(0.014, 60.); // (cm) pitch=0.014, (deg) 1st axis angle/deg
00204   myMkSM->setTransDiffusion(0.017); // (cm per 1 cm of path), default=0.017 cm/cm
00205 
00206   myMkCL->setSeedStripThres(30.); // a.u. 
00207   myMkCL->setStripNoiseSigma(12.); // a.u. , default=12=10%
00208 
00209   // chain is ready
00210   chain->ls(3);
00211 
00212 
00213   Int_t iInit = chain->Init();
00214 
00215 
00216 
00217 
00218   if (iInit >=  kStEOF) {chain->FatalErr(iInit,"on init"); return;}
00219   if (Last == 0) return;
00220   StEvtHddr *hd = (StEvtHddr*)chain->GetDataSet("EvtHddr");
00221   if (hd) hd->SetRunNumber(-2); // to be sure that InitRun calls at least once
00222     // skip if any
00223   chain->EventLoop(First,Last,0);
00224   gMessMgr->QAInfo() << "Run completed " << endm;
00225   gSystem->Exec("date");
00226 
00227 
00228   //FGT end....
00229   {
00230     TDatime t;
00231     gMessMgr->QAInfo() << Form("Run is finished at Date/Time %i/%i",
00232                                t.GetDate(),t.GetTime()) << endm;
00233   }
00234   myMkSM->saveHisto(fnameOut2);// one is enough, all histos are merged
00235 
00236 
00237 
00238 
00239 }
00240 //_____________________________________________________________________
00241 void bfc(Int_t Last, 
00242          const Char_t *Chain,
00243          const Char_t *infile, 
00244          const Char_t *outfile, 
00245          const Char_t *TreeFile)
00246 {
00247   bfc(1,Last,Chain,infile,outfile,TreeFile);
00248 }
00249 //____________________________________________________________
00250 void Usage() {
00251   printf ("============= \t U S A G E =============\n");
00252   printf ("bfc(Int_t First,Int_t Last,const Char_t *Chain,const Char_t *infile,const Char_t *outfile,const Char_t *TreeFile)\n");
00253   printf ("bfc(Int_t Last,const Char_t *Chain,const Char_t *infile,const Char_t *outfile,const Char_t *TreeFile)\n");
00254   printf ("bfc(const Char_t *ChainShort,Int_t Last,const Char_t *infile,const Char_t *outfile)\n");
00255   printf ("where\n");
00256   printf (" First     \t- First event to process\t(Default = 1)\n");
00257   printf (" Last      \t- Last  event to process\t(Default = 1)\n");
00258   printf (" Chain     \t- Chain specification   \t(without First &  Last: Default is \"\" which gives this message)\n");
00259   printf ("           \t                        \t with    First || Last: Default is \"gstar tfs\")\n");
00260   printf (" infile    \t- Name of Input file    \t(Default = 0, i.e. use preset file names depending on Chain)\n"); 
00261   printf (" outfile   \t- Name of Output file   \t(Default = 0, i.e. define Output file name from Input one)\n");
00262   printf (" outfile   \t- Name of Tree File     \t(Default = 0, i.e. define Output file name from Input one (tags TNtuple))\n");
00263   printf (" ChainShort\t- Short cut for chain   \t(Default = \"\" -> print out of this message)\n");
00264   gSystem->Exit(1);
00265 }
00266 //_____________________________________________________________________
00267 void bfc(const Char_t *Chain="ittf") {
00268   bfc(-2,Chain);
00269 }