rnegrao's blog

rnegrao's home page
Posts

2013

November (1)
October (1)
August (2)
March (1)
February (1)

2012

October (1)

2011

September (1)
July (2)
April (4)

STAR Protected

You must register or login in order to post into this group.

User login

Navigation

Group notifications

This group offers an RSS feed. Or subscribe to these personalized, sitewide feeds:

All posts:
- Feed
- Page

Temporary Codes

Updated on Tue, 2013-11-19 10:16. Originally created by rnegrao on 2011-10-20 12:52.

re starting on node: rcas6410.rcf.bnl.gov
HOME is not defined
Using getent method
HOME is now /star/u/rnegrao
--------------------------------
STAR Unified Meta Scheduler 1.10.4 we are starting in $SCRATCH : /home/tmp/rnegrao/5F1D60488317E4272DC2ED2E7C793620_0
--------------------------------
total 0
  *******************************************
  *                                         *
  *        W E L C O M E  to  R O O T       *
  *                                         *
  *   Version   5.22/00  17 December 2008   *
  *                                         *
  *  You are welcome to visit our Web site  *
  *          http://root.cern.ch            *
  *                                         *
  *******************************************
 
ROOT 5.22/00 (trunk@26997, De 10 2010, 13:19:00 on linux)
 
CINT/ROOT C/C++ Interpreter version 5.16.29, Jan 08, 2008
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.
*** Float Point Exception is OFF ***
 *** Start at Date : Thu Oct 20 13:08:28 2011
QAInfo:You are using STAR_LEVEL : dev, ROOT_LEVEL : 5.22.00 and node : rcas6410.rcf.bnl.gov
 
Processing do.C(9999,"/star/u/rnegrao/jobs/test/sched5F1D60488317E4272DC2ED2E7C793620_0.list","/star/u/rnegrao/jobs/test/st_physics_5F1D60488317E4272DC2ED2E7C793620_0.root")...
strangeArrayNames-arrayNames-__NARRAYS__ = 0
 loading of shared libraries done
##### void StMuDstMaker::setProbabilityPidFile(const char*) ##### ATTENTION: pid table hardwired to /afs/rhic.bnl.gov/star/packages/DEV/StarDb/dEdxModel/PIDTableP01gl.root
MuDst pointer are created!
adc ok!
StInfo: StDbManagerImpl::lookUpServers(): config file name is /afs/rhic.bnl.gov/star/packages/conf/dbLoadBalancerLocalConfig_BNL.xml
There is no 1  /star/u/rnegrao/dbServers.xml
 looking at 2 /afs/rhic.bnl.gov/star/packages/DEV/StDb/servers/dbServers.xml
StInfo: StDbManagerImpl::lookUpServer() line=679
******** Order of Files searched for dbServers *********
  3. /afs/rhic.bnl.gov/star/packages/DEV/StDb/servers/dbServers.xml
********************************************************
 
StInfo: MysqlDb::Connect: Load balancer took 53 ms, will use db15.star.bnl.gov:3316
StInfo: Ssl_cipher = AES128-SHA
StInfo: MysqlDb::Connect(host,user,pw,database,port) line=474 Server Connecting: DB=StarDb  Host=db15.star.bnl.gov:3316
 --> Connection Time=0.238845 sec    MysqlVersion=5.0.77
StInfo: Configuration for BEMC hit reconstruction
StInfo:   Configuration for detector bemc
StInfo:      switch for deducting pedestal     = 1
StInfo:      switch for calibration            = 1
StInfo:      cutoff type                       = 0
StInfo:      cutoff value                      = -1
StInfo:      save only calibrated hits         = 0
StInfo:      save only if status is ok         = 1
StInfo:      save only if crate is ok          = 1
StInfo:      SAVE ALL FLAG (overwrites above)  = 0
StInfo:   Configuration for detector bprs
StInfo:      switch for deducting pedestal     = 1
StInfo:      switch for calibration            = 1
StInfo:      cutoff type                       = 1
StInfo:      cutoff value                      = 1.5
StInfo:      save only calibrated hits         = 0
StInfo:      save only if status is ok         = 1
StInfo:      save only if crate is ok          = 1
StInfo:      SAVE ALL FLAG (overwrites above)  = 0
StInfo:   Configuration for detector bsmde
StInfo:      switch for deducting pedestal     = 1
StInfo:      switch for calibration            = 1
StInfo:      cutoff type                       = 1
StInfo:      cutoff value                      = 1.5
StInfo:      save only calibrated hits         = 0
StInfo:      save only if status is ok         = 1
StInfo:      save only if crate is ok          = 1
StInfo:      SAVE ALL FLAG (overwrites above)  = 0
StInfo:   Configuration for detector bsmde
StInfo:      switch for deducting pedestal     = 1
StInfo:      switch for calibration            = 1
StInfo:      cutoff type                       = 1
StInfo:      cutoff value                      = 1.5
StInfo:      save only calibrated hits         = 0
StInfo:      save only if status is ok         = 1
StInfo:      save only if crate is ok          = 1
StInfo:      SAVE ALL FLAG (overwrites above)  = 0
StInfo:   Configuration for detector bsmdp
StInfo:      switch for deducting pedestal     = 1
StInfo:      switch for calibration            = 1
StInfo:      cutoff type                       = 1
StInfo:      cutoff value                      = 1.5
StInfo:      save only calibrated hits         = 0
StInfo:      save only if status is ok         = 1
StInfo:      save only if crate is ok          = 1
StInfo:      SAVE ALL FLAG (overwrites above)  = 0
StInfo: Getting database tables for the BEMC detector
StInfo: MysqlDb::Connect: Load balancer took 36 ms, will use db15.star.bnl.gov:3316
StInfo: Ssl_cipher = AES128-SHA
StInfo: MysqlDb::Connect(host,user,pw,database,port) line=474 Server Connecting: DB=Conditions_trg  Host=db15.star.bnl.gov:3316
 --> Connection Time=0.0576891 sec    MysqlVersion=5.0.77
StInfo: loaded a new              bemcPed table with beginTime 2011-04-19 00:00:00 and endTime 2011-05-03 04:00:03
StInfo: loaded a new           bemcStatus table with beginTime 2011-04-19 00:00:00 and endTime 2011-05-06 22:01:53
StInfo: loaded a new            bemcCalib table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new             bemcGain table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new              bprsPed table with beginTime 2008-03-04 10:30:56 and endTime 2037-12-31 12:00:00
StInfo: loaded a new           bprsStatus table with beginTime 2008-12-15 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new            bprsCalib table with beginTime 1999-01-01 00:10:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new             bprsGain table with beginTime 1999-01-01 00:08:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new             bsmdePed table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new          bsmdeStatus table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new           bsmdeCalib table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new            bsmdeGain table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new             bsmdpPed table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new          bsmdpStatus table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new           bsmdpCalib table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new            bsmdpGain table with beginTime 2010-12-20 00:00:00 and endTime 2037-12-31 12:00:00
StInfo: loaded a new    bemcTriggerStatus table with beginTime 2010-05-28 22:00:04 and endTime 2037-12-31 12:00:00
StInfo: loaded a new       bemcTriggerPed table with beginTime 2010-04-29 21:30:03 and endTime 2037-12-31 12:00:00
StInfo: loaded a new       bemcTriggerLUT table with beginTime 2010-05-28 22:00:04 and endTime 2037-12-31 12:00:00
StWarning: Could not make BEMC detector
mMuDst OK!
mMuEvent OK!
mMuEmcCollection OK!
mEmcCollection OK!
mEmcGeometry OK!
9
StWarning: Could not make BEMC detector

Int_t StEmcMaker::Make()
{
	//Getting Makers I need and check if the pointer exists
	//mMuDstMaker = (StMuDstMaker *)GetMaker("StMuDst");
	//if(!mMuDstMaker) cout<<"sey ha!"<<endl;//return kStErr;
 
 
	//Getting StMuDst
	mMuDst = (StMuDst*)GetInputDS("StMuDst");
	//StMuDstMaker *mMuDst = (StMuDstMaker*)GetMaker("MuDst");
	if(mMuDst) cout<<"mMuDst OK!"<<endl;//return kStSkip;
	else cout<<"no mMuDst"<<endl;
 
	//Getting StMuEvent
	mMuEvent = mMuDst->event();
	if(mMuEvent) cout<<"mMuEvent OK!"<<endl;
	else cout<<"no mMuEvent"<<endl;
 
	//Getting muEmcCollection	
	mMuEmcCollection = mMuDst->muEmcCollection();
	if(mMuEmcCollection) cout<<"mMuEmcCollection OK!"<<endl;
	else cout<<"no mMuEmcCollection"<<endl;
 
	//Getting EmcCollection
	mEmcCollection = (StEmcCollection*)mADCMaker->getEmcCollection();
	//mEmcCollection = mMuDst->emcCollection();
	if(mEmcCollection) cout<<"mEmcCollection OK!"<<endl;
	else cout<<"no mEmcCollection"<<endl;
 
	//Getting BEMC Geometry
	mEmcGeometry = StEmcGeom::instance("bemc");
	if(mEmcGeometry) cout<<"mEmcGeometry OK!"<<endl;
	else cout<<"no mEmcGeometry"<<endl;	
 
	cout<<kBarrelEmcTowerId<<endl;
 
	return kStOK;
}

///////////////////////////////////////////////////////////////////////////////
//
// $Id:  $
//
///////////////////////////////////////////////////////////////////////////////
 
Int_t    usePath = 0;
Int_t    nFile = 0;
TString  thePath;
TString  theFileName;
TString  originalPath;
class    StChain;
StChain  *chain=0;
class StEventDisplayMaker;
StEventDisplayMaker *dsMaker = 0;
TBrowser *b=0;
 
 
const char *dstFile = 0;
const char *xdfFile = 0;
const char *mdcFile = 0;
const char *fileList[] = {dstFile,xdfFile,mdcFile,0};
/*
void do(Int_t nevents=9999)
 
void do(Int_t nevents=9999, 
	    const Char_t *fileIn="/star/data10/reco/pp500_production_2011/ReversedFullField/dev/2011/060/12060048/st_physics_12060048_raw_5020001.MuDst.root",
	    const Char_t *fileOut="test.root")
*/
void do(Int_t nevents,
	   const Char_t *fileIn,
	   const Char_t *fileOut)
 
{ 
    //
    // First load some shared libraries we need
    //
 
    gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
    loadSharedLibraries();
    gSystem->Load("StDetectorDbMaker");
    gSystem->Load("StTpcDb");
    gSystem->Load("StDbUtilities");
    gSystem->Load("StDbLib");
    gSystem->Load("StDbBroker");
    gSystem->Load("St_db_Maker");
    gSystem->Load("libgeometry_Tables");
    gSystem->Load("StDaqLib");
    gSystem->Load("StEmcRawMaker");
    gSystem->Load("StPreEclMaker");
    gSystem->Load("StEpcMaker");
 
    gSystem->Load("St_base");
    gSystem->Load("StUtilities");
    gSystem->Load("StChain");
    gSystem->Load("StBFChain"); 
    gSystem->Load("St_Tables");
 
    gSystem->Load("libgen_Tables");
    gSystem->Load("libsim_Tables");
    gSystem->Load("libglobal_Tables");
    gSystem->Load("geometry");
 
    gSystem->Load("StIOMaker");
    gSystem->Load("StTreeMaker");
    gSystem->Load("StarClassLibrary");
    gSystem->Load("StBichsel");
 
    gSystem->Load("StEvent");
    gSystem->Load("StEmcUtil");
    gSystem->Load("StStrangeMuDstMaker");
    gSystem->Load("StMuDSTMaker");
 
    gSystem->Load("StEmcADCtoEMaker");
 
 
//    gSystem->Load(".sl53_gcc432/obj/StRoot/StAnalysisMaker/StAnalysisMaker.so");
    gSystem->Load("StEmcMaker");
 
    //
    // instantiate the chain
    //
    chain  = new StChain("StChain");
 
    //
    // Input Maker
    //
    //StMuDstMaker *myMuDstMaker = new StMuDstMaker(0,0,"",fileIn,"",10000,"myMuDstMaker");
    StMuDstMaker *myMuDstMaker = new StMuDstMaker(0,0,"",fileIn,"",nevents,"StMuDst");
    if (myMuDstMaker) cout<<"MuDst pointer are created!"<<endl;
 
    StMuDbReader *db  = StMuDbReader::instance(); 
    St_db_Maker *dbMk = new St_db_Maker("StarDb","MySQL:StarDb");
    StEmcADCtoEMaker *adc=new StEmcADCtoEMaker();
    if(adc) cout<<"adc ok!"<<endl;
    else    cout<<"acd not ok!"<<endl;
    StPreEclMaker *ecl=new StPreEclMaker();
 
    //
    // QA Analysis
    //
    //StAnalysisMaker Analysis = new StAnalysisMaker(muDstMaker);
    //StHadEtMaker *Analysis = new StHadEtMaker(muDstMaker,energy);
    //StEmcMaker *Analysis = new StEmcMaker("Analysis");
    StEmcMaker *Analysis = new StEmcMaker("Analysis",adc);
    Analysis->setIO(fileIn,fileOut);
 
    //
    // Initialize chain
    //
 
    Int_t iInit = chain->Init();
    if (iInit) chain->Fatal(iInit,"on init");
 
    //
    // Event loop
    //
 
    int istat=0,i=1;
 EventLoop: if (i <= nevents && istat!=2 && istat!=3 && istat!=4) {
      chain->Clear();
      istat = chain->Make(i);
      if (istat==2) 
	{cout << "Last  event processed. Status = " << istat << endl;}
      if (istat==3) 
	{cout << "Error event processed. Status = " << istat << endl;}
      i++;
      goto EventLoop;
    }
    chain->Finish();
    i--;
}

Groups:

STAR Protected

rnegrao's blog
Login or register to post comments

The STAR experiment

rnegrao's blog

STAR Protected

User login

Navigation

Group notifications

Temporary Codes