runEvent.C

TString opt = "in fcsWFF fcsCluster fcsPoint";

void runEvent(Int_t nevents=100, Int_t pedLedPhy=2, Int_t eventDisplay=100, Int_t debug=0,
            const Char_t *file="st_physics_22066021_raw_0000002.event.root"){

    TString f(file);
    int l1 = f.Length();
    int l2 = f.Last('/');
    TString dir(f);  dir.Remove(l2+1,l1-l2);
    TString fne(f);  fne.Remove(0,l2+1);
    TString crun(fne);
    int l3= crun.First('_'); crun.Remove(0,l3+1);
    int l4= crun.First('_'); crun.Remove(0,l4+1);
    int l5 =crun.Length(); int l6= crun.First('_'); crun.Remove(l6,l5-l6);
    int run=crun.Atoi();
    int day=run/1000;    
    //printf("l1=%d l2=%d l3=%d l4=%d l5=%d l6=%d\n",l1,l2,l3,l4,l5,l6);
    //printf("f  =%s\n",f.Data());
    //printf("dir=%s\n",dir.Data());
    //printf("fne=%s\n",fne.Data());
    //printf("run=%s\n",run.Data());
    printf("Run=%d yearday=%d\n",run,day);
    char edout[200];

    gROOT->LoadMacro("bfc.C");  // Load big "full" chain
    bfc(-1,opt,file);     // Setup but do not init
        
    St_db_Maker *dbMk= (St_db_Maker*) chain->GetMaker("db");
    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=(StFcsDbMaker *)chain->GetMaker("fcsDbMkr");
    StFcsDb* fcsDb = (StFcsDb*) chain->GetDataSet("fcsDb");
    fcsDb->setDbAccess(0);
    fcsDb->setDebug(debug);
    fcsDb->forceUniformGain(0.0053*2.0); //tail has 1/2 integral
    fcsDb->forceUniformGainCorrection(1.0);

    StFcsWaveformFitMaker *wff=(StFcsWaveformFitMaker *)chain->GetMaker("StFcsWaveformFitMaker");
    StFcsClusterMaker *clu=(StFcsClusterMaker *)chain->GetMaker("StFcsClusterMaker");
    StFcsPointMaker *poi=(StFcsPointMaker *)chain->GetMaker("StFcsPointMaker");
    wff->setDebug(1);
    wff->setEnergySelect(10);    
    char wffout[100];
    char wffout2[100];
    if(pedLedPhy==1){
        if(run<1106734){
            wff->setCenterTimeBins(207,190,240);
        }else{
            wff->setCenterTimeBins(207+8,190+8,240+8);
        }      
        wff->setPedTimeBins(0,150);
        skip=5;
    }else{
        wff->setCenterTimeBins(50,30,100);
    }
    sprintf(wffout,"%d/%d.wff",day,run);
    wff->setFileName(wffout,10,0);
    sprintf(wffout2,"%d/%d.wfftime.png",day,run);
    wff->setMeasureTime(wffout2);   

    gSystem->Load("StFcsEventDisplay");
    StFcsEventDisplay* fcsed;
    if(pedLedPhy>0 && eventDisplay>0){
        gSystem->Load("StEpdUtil");
        fcsed = new StFcsEventDisplay();
        fcsed->setMaxEvents(eventDisplay);
        sprintf(edout,"%d/%d.eventDisplay.png",day,run);
        fcsed->setFileName(edout);
        fcsed->setFilter(1);
    }

    gSystem->Load("StFcsQaMaker");
    StFcsQaMaker *qaMkr=new StFcsQaMaker("FcsQa");  
    qaMkr->setRun(run);
    if(pedLedPhy==0){
        qaMkr->setMaxTimeBins(1024);
        qaMkr->setMaxAdc(512);      
    }else if(pedLedPhy==1){
        qaMkr->setMaxTimeBins(256);
        qaMkr->setSumTimeBins(210,225);  
        qaMkr->setMaxAdc(2048);      
        qaMkr->setMaxAdcSum(20000);      
    }else if(pedLedPhy==2){
        qaMkr->setMaxTimeBins(160);
        qaMkr->setSumTimeBins(35,60);  
        qaMkr->setMaxAdc(1024);      
        qaMkr->setMaxAdcSum(2000);      
    }

    // Intializet the chain
    chain->Init();

    // Process events
    chain->EventLoop(nevents);    
}