runMudst.C

void runMudst(char* file="st_cosmic_adc_22326042_raw_0000005.MuDst.root", 
	      int ifile=-1, Int_t nevt=10, char* outdir=".", int readMuDst=1, int debug=0){  
    gROOT->Macro("Load.C");
    gROOT->Macro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
    gSystem->Load("StEventMaker");
    gSystem->Load("StFcsDbMaker");
    gSystem->Load("StFcsRawHitMaker");
    gSystem->Load("StFcsWaveformFitMaker");
    gSystem->Load("StFcsClusterMaker");
    gSystem->Load("libMinuit");
    gSystem->Load("StFcsPointMaker");

    StChain* chain = new StChain("StChain"); chain->SetDEBUG(0);
    StMuDstMaker* muDstMaker = new StMuDstMaker(0, 0, "", file,".", 1000, "MuDst");
    int n=muDstMaker->tree()->GetEntries();
    printf("Found %d entries in Mudst\n",n);
    int start=0, stop=n;
    if(ifile>=0){
	int start=ifile*nevt;
	int stop=(ifile+1)*nevt-1;
	if(n<start) {printf(" No event left. Exiting\n"); return;}
	if(n<stop)  {printf(" Overwriting end event# stop=%d\n",n); stop=n;}
    }else if(nevt>=0 && nevt<n){
	stop=nevt;
    }else if(nevt==-2){
	stop=2000000000; 
    }
    printf("Doing Event=%d to %d\n",start,stop);
    
    St_db_Maker* dbMk = new St_db_Maker("db","MySQL:StarDb","$STAR/StarDb"); 
    if(dbMk){
	dbMk->SetAttr("blacklist", "tpc");
	dbMk->SetAttr("blacklist", "svt");
	dbMk->SetAttr("blacklist", "ssd");
	dbMk->SetAttr("blacklist", "ist");
	dbMk->SetAttr("blacklist", "pxl");
	dbMk->SetAttr("blacklist", "pp2pp");
	dbMk->SetAttr("blacklist", "ftpc");
	dbMk->SetAttr("blacklist", "emc");
	dbMk->SetAttr("blacklist", "eemc");
	dbMk->SetAttr("blacklist", "mtd");
	dbMk->SetAttr("blacklist", "pmd");
	dbMk->SetAttr("blacklist", "tof");
	dbMk->SetAttr("blacklist", "etof");
	dbMk->SetAttr("blacklist", "rhicf");
    }
    
    StFcsDbMaker *fcsDbMkr= new StFcsDbMaker();
    StFcsDb* fcsDb = (StFcsDb*) chain->GetDataSet("fcsDb");
    //fcsDb->setReadGainFromText();
    fcsDb->setReadGainCorrFromText();
    StEventMaker* eventMk = new StEventMaker();
    StFcsRawHitMaker* hit = new StFcsRawHitMaker();  
    hit->setReadMuDst(readMuDst);
    StFcsWaveformFitMaker *wff= new StFcsWaveformFitMaker();
    //wff->setEnergySelect(10);
    wff->SetDebug(debug);
    StFcsClusterMaker *clu= new StFcsClusterMaker();
    StFcsPointMaker *poi= new StFcsPointMaker();
    clu->SetDebug(debug);
    poi->SetDebug(debug);

     gSystem->Load("StVpdCalibMaker");
   StVpdCalibMaker *vpdCalib = new StVpdCalibMaker();
   vpdCalib->setMuDstIn();

    gSystem->Load("StFcsPi0FinderForEcal");
    StFcsPi0FinderForEcal* fcsPi0Finder = new StFcsPi0FinderForEcal();
    fcsPi0Finder->st(1);

    chain->Init();
    chain->EventLoop(start,stop);
    chain->Finish();
    delete chain;
}