StMiniHijing.C

// $Id: StMiniHijing.C,v 1.7 2006/07/24 19:25:26 calderon Exp $
// owner: Manuel Calderon de la Barca Sanchez
//
// what it does: reads .geant.root file from hijing data, produces minimc.root file 
//               runs a chain with the makers: 
//                 StMcEventMaker,StAssociationMaker,
//                 StMiniMcEventMaker
// $Log: StMiniHijing.C,v $
// Revision 1.7  2006/07/24 19:25:26  calderon
// Load EEmcUtil, needed by StMcEvent.
//
// Revision 1.6  2004/03/30 03:16:15  calderon
// Modifications for running in bfc.
//  - Changed to use StiIOInterface (IOMaker in normal mode, TreeMaker in bfc)
//  - Cleaned up Init(), InitRun() to handle the changing file names.
//  - Initialize lots of variables and pointers in constructor.
//  - Delete some pointers in Finish (deleting the TTree causes a seg fault, though.)
//  - Note that currently the StHits in the ITTF chain don't have a usedInFit() flag,
//    so there will be many messages complaining about this.
//  - Removed the mDebug data member, every Maker already has one, so change
//    to use that throughout the package.
//
// Revision 1.5  2003/07/09 01:07:23  calderon
// Addition of FTPC reference multiplicity
// Addition of other multiplicity values for StMiniMcEvent
// Changes to reflect the use of the setters and getters, no longer
// access the data members directly.
//
// Revision 1.4  2002/06/28 22:15:12  calderon
// Changes to deal with seg. faults in the file name handling:
// Conventions:
// StMiniMcMaker looks for the input file from the IO maker to figure out
// if the file has changed.  This is done using TString::Contains() in Make().
// Usually we will run one file at a time, but in order not to break Bum's scheme of being
// able to process several files in one go, this is left as is.  However, for
// embedding, the file name is not enough, in Eric's new scheme there are repeated
// file names.  This is resolved by adding a prefix to the output file name.  However,
// this prefix should not be overwritten, so the current code only replaces the
// string inside the output file name pertaining to the input file name, and leaves
// the prefix of the output file intact.  This was done for embedding looking for
// st_physics, and here is where the problem arose: hijing files begin with a different
// prefix.  To solve this problem, the input file name prefix is now an input parameter
// in the macro.
//
// StMiniEmbed.C and StMiniHijing.C now conform to this convention.  StMiniEmbed.C
// did not change its prototype, because all embedding files have st_phyics as prefix.
// StMiniHijing.C changed its prototype, now it takes as an input argument the prefix,
// but in order not to break Jenn's scripts if she was already using this macro,
// this parameter was added at the end and defaults to "rcf", which is appropriate
// for hijing files reconstructed in rcf.
//
// Revision 1.3  2002/06/07 02:21:59  calderon
// Protection against empty vector in findFirstLastHit
// $Log: StMiniHijing.C,v $
// Revision 1.7  2006/07/24 19:25:26  calderon
// Load EEmcUtil, needed by StMcEvent.
//
// Revision 1.6  2004/03/30 03:16:15  calderon
// Modifications for running in bfc.
//  - Changed to use StiIOInterface (IOMaker in normal mode, TreeMaker in bfc)
//  - Cleaned up Init(), InitRun() to handle the changing file names.
//  - Initialize lots of variables and pointers in constructor.
//  - Delete some pointers in Finish (deleting the TTree causes a seg fault, though.)
//  - Note that currently the StHits in the ITTF chain don't have a usedInFit() flag,
//    so there will be many messages complaining about this.
//  - Removed the mDebug data member, every Maker already has one, so change
//    to use that throughout the package.
//
// Revision 1.5  2003/07/09 01:07:23  calderon
// Addition of FTPC reference multiplicity
// Addition of other multiplicity values for StMiniMcEvent
// Changes to reflect the use of the setters and getters, no longer
// access the data members directly.
//
// Revision 1.4  2002/06/28 22:15:12  calderon
// Changes to deal with seg. faults in the file name handling:
//

class StChain;
StChain *chain=0;

void StMiniHijing(Int_t nevents=3,
                  const char* MainFile=
                  "hijing_b0_3/rcf0147_01*geant.root",
                  const char* outDir = "./",
                  int commonHits=3,
                  TString filePrefix="rcf")
{

  // Dynamically link needed shared libs
  //gSystem->Load("StarRoot");
  gSystem->Load("St_base");
  gSystem->Load("StChain");

  gSystem->Load("St_Tables");
  gSystem->Load("StUtilities");
  gSystem->Load("StIOMaker");
  gSystem->Load("StarClassLibrary");
    
  gSystem->Load("StDetectorDbMaker");
  gSystem->Load("StTpcDb");
  gSystem->Load("StEvent");
  gSystem->Load("StEventMaker"); 
  gSystem->Load("StEmcUtil"); 
  gSystem->Load("StEEmcUtil");
  gSystem->Load("StMcEvent");
  gSystem->Load("StMcEventMaker");
  gSystem->Load("StAssociationMaker");
  gSystem->Load("StMcAnalysisMaker");

  gSystem->Load("StMiniMcEvent");
  gSystem->Load("StMiniMcMaker");
//   gSystem->Load("Common");

  chain = new StChain("StChain"); 
  chain->SetDebug();
  
  // Now we add Makers to the chain...
  
  StIOMaker* ioMaker = new StIOMaker("IO","r",MainFile,"bfcTree");
  ioMaker->SetDebug();
  ioMaker->SetIOMode("r");
  ioMaker->SetBranch("*",0,"0");                 //deactivate all branches
  ioMaker->SetBranch("geantBranch",0,"r"); //activate geant Branch
  ioMaker->SetBranch("eventBranch",0,"r"); //activate Event Branch
  // ioMaker->SetBranch("runcoBranch",0,"r"); //activate runco Branch
  
  //     const char *mainDB = "MySQL:Geometry_tpc";
  //     St_db_Maker *dbMk = new St_db_Maker("Geometry",mainDB);
  //     dbMk->SetDebug();
  
  //     const char *calibDB = "MySQL:Calibrations_tpc";
  //     St_db_Maker *calibMk = new St_db_Maker("Calibrations",calibDB);
  //     calibMk->SetDebug();
  
  //     StTpcDbMaker *tpcDbMk = new StTpcDbMaker("tpcDb");
  
  // Note, the title "events" is used in the Association Maker, so don't change it.
//   StEventMaker*       eventReader   = new StEventMaker("events","title");
//   eventReader->doPrintMemoryInfo = kFALSE;

  StMcEventMaker*     mcEventReader = new StMcEventMaker; // Make an instance...
  //     mcEventReader->doPrintMemoryInfo = kFALSE;
  //     mcEventReader->doUseTpc = kTRUE;
  //     mcEventReader->doUseSvt = kTRUE;
  //     mcEventReader->doUseFtpc = kTRUE;
  //     mcEventReader->doUseRich = kTRUE;
  StAssociationMaker* associator    = new StAssociationMaker;
  //associator->doPrintMemoryInfo = kTRUE;

  //
  // mdst of association stuff
  //

  TString filename   = MainFile;
  int fileBeginIndex = filename.Index(filePrefix,0);
  filename.Remove(0,fileBeginIndex);

  
  StMiniMcMaker *krap = new StMiniMcMaker;
  krap->SetDebug(1);
  krap->setGhost();
  krap->setOutDir(outDir);
  krap->setPtCut(0.);
  krap->setFileName(filename);
  krap->setFilePrefix(filePrefix);

  // Define the cuts for the Associations
  
  StMcParameterDB* parameterDB = StMcParameterDB::instance();  
  // TPC
  parameterDB->setXCutTpc(.5); // 5 mm
  parameterDB->setYCutTpc(.5); // 5 mm
  parameterDB->setZCutTpc(.5); // 5 mm
  parameterDB->setReqCommonHitsTpc(commonHits); // Require 3 hits in common for tracks to be associated
  // FTPC
  parameterDB->setRCutFtpc(.3); // 3 mm
  parameterDB->setPhiCutFtpc(5*(3.1415927/180.0)); // 5 degrees
  parameterDB->setReqCommonHitsFtpc(3); // Require 3 hits in common for tracks to be associated
  // SVT
  parameterDB->setXCutSvt(.08); // 800 um
  parameterDB->setYCutSvt(.08); // 800 um
  parameterDB->setZCutSvt(.08); // 800 um
  parameterDB->setReqCommonHitsSvt(1); // Require 1 hits in common for tracks to be associated
  
  
  // now execute the chain member functions
  
  //chain->PrintInfo();
  Int_t initStat = chain->Init(); // This should call the Init() method in ALL makers
  if (initStat) chain->Fatal(initStat, "during Init()");
  
  int istat=0,iev=1;
 EventLoop: if (iev<=nevents && istat!=2) {
   chain->Clear();
   cout << "---------------------- Processing Event : " << iev << " ----------------------" << endl;
   istat = chain->Make(iev); // This should call the Make() method in ALL makers
   if (istat == 2) { cout << "Last  Event Processed. Status = " << istat << endl; }
   if (istat == 3) { cout << "Error Event Processed. Status = " << istat << endl; }
   iev++; goto EventLoop;
 } // Event Loop
 
 chain->Finish();
  
}