prod_sof.C

Go to the documentation of this file.

int prod_sof(TString outpre="", Float_t Mom=15.15, int mode=9999, int run=-1, int applyfull=0)
{
        // Parameter file
        TString inParFile = outpre+"_par.root"; 
        
        // PID table with selection thresholds; can be modified by the user
        TString pidParFile = TString(gSystem->Getenv("VMCWORKDIR"))+"/macro/params/all.par";    
                
        // Output file
        TString outFile = outpre+"_sof.root";
        
        // ------------------------------------------------------------------------
        FairLogger::GetLogger()->SetLogToFile(kFALSE);
        
        // -----   Reconstruction run   -------------------------------------------
        FairRunAna *fRun= new FairRunAna();
        fRun->SetInputFile(outpre+"_pid.root");
        fRun->SetOutputFile(outFile);
        fRun->SetGenerateRunInfo(kFALSE);
                
        // *** initialization
        FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
        
        // *** setup parameter database         
        FairParRootFileIo* parIO = new FairParRootFileIo();
        parIO->open(inParFile);
        FairParAsciiFileIo* parIOPid = new FairParAsciiFileIo();
        parIOPid->open(pidParFile.Data(),"in");
        
        rtdb->setFirstInput(parIO);
        rtdb->setSecondInput(parIOPid);
        rtdb->setOutput(parIO);  
        
        //---------------------Create and Set the Field(s)---------- 
        PndMultiField *fField= new PndMultiField("AUTO");
        fRun->SetField(fField);

        //RhoCalculationTools::ForceConstantBz(20.0);
  
        // ***********************
        // *** SoftTriggerTask ***
        // ***********************
        
        int modeshort = mode;
        if (modeshort>1000000) modeshort /= 1000;
        if (modeshort>1000)    modeshort %= 1000;
        
        cout <<" ****** mode:"<<mode<<"  modeshort:"<<modeshort<<endl;
        
        // this file contains the trigger line definitions
        //TString triggercfg   = TString(gSystem->Getenv("VMCWORKDIR"))+"/macro/softrig/triggerlines_10ch.cfg";
        TString triggercfg   = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/triggerlines.cfg";
        
        // this file contains the cut setup for various modes
        //TString selectioncfg = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/selection_10ch_tight.cfg"; 
        TString selectioncfg = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/selection.cfg"; 
        
        PndSoftTriggerTask *stTask = new PndSoftTriggerTask(Mom, mode, run, triggercfg);
        stTask->SetConfigurationFile(selectioncfg);
        
        //stTask->McMatchAllowPhotos(10, 1.0);
                
        stTask->ApplyFullSelection(applyfull);
        
        stTask->SetPi0SignalParams(0.134, 0.0045);
        stTask->SetKs0SignalParams(0.497, 0.0085);
        stTask->SetEtaSignalParams(0.549, 0.0055); // needs to be refined!
        
        // set PID algos
        TString algo = "PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts;PidAlgoMvd";       
        stTask->SetPidAlgoAll(algo);
        
        stTask->SetTagAll(true);                // tag all modes
        //stTask->SetTagMode(modeshort);  // switch single tags on/off; mode number has to match one from config file
        
        if (modeshort==900)     
        {
                stTask->SetQAAll(true);         // ntuple output for all modes
        }
        else
        {
                stTask->SetQAAll(false);
                stTask->SetQAMode(modeshort);
                stTask->SetQAMctOnly();
        }
        
        stTask->SetQAEvent(true);               // event info
        stTask->SetQAMc(true);
        
        stTask->SetGammaMinE(0.10);             // global energy pre-cut for neutrals 
        stTask->SetTrackMinP(0.10);             // global momentum pre-cut for charged  
        stTask->SetInitialPidCut(0.1);  // global PID pre-cut for charged       
        stTask->SetDstMDiffCut(0.1);    // special cut on D*-D mass difference (to reduce comb and output file size)
        
        //stTask->McMatchAllowPhotos(1,0.05);  // MC truth match shall partially ignore photos photons (here max 1 with E<50 MeV)

        fRun->AddTask(stTask);

        // ***********************
        // *** SoftTriggerTask ***
        // ***********************
        
        
        // *** and run analysis
        fRun->Init(); 
        fRun->Run(0,0);
        
        //gObjectTable->Print();
  return 0;
}