ana_jpsi_task_ST.C

Go to the documentation of this file.

// ********************************************
// Example analysis in a task with accessing the OnlineEventFilterInfo
// written by the SoftTriggerTask
// ********************************************

int ana_jpsi_task_ST(TString Fname="test", int nevts=0, double pbarmom = 6.232069)
{
        // *** set this to your output path
        TString WorkDir = TString(gSystem->Getenv("VMCWORKDIR"));
        TString OutPath = WorkDir+"/macro/softrig";     
        TString OutFile = OutPath + "/" + Fname(Fname.Last('/')+1,Fname.Length()); // cut away input path
                
        // *** the output file for FairRunAna
        TString InFile  = Fname;
        if (!InFile.EndsWith(".root")) InFile+="_fast.root";
                                        
        // *** initialization
        FairLogger::GetLogger()->SetLogToFile(kFALSE);

        FairRunAna* fRun = new FairRunAna();
        fRun->SetGenerateRunInfo(kFALSE);
        fRun->SetInputFile(InFile);
        fRun->SetOutputFile(OutFile+"_loose.root");

        // *** take constant field; needed for PocaVtx
        RhoCalculationTools::ForceConstantBz(20.0);
        
        // --------------------------------
        // *** SoftTriggerTask ***
        // --------------------------------
        
        // this file contains the trigger line definitions
        TString triggercfg   = WorkDir+"/softrig/triggerlines.cfg";
        
        // this file contains the cut setup for the 10 channels from soft trigger report
        TString selectioncfg = WorkDir+"/softrig/selection_10ch_loose.cfg"; 
        
        PndSoftTriggerTask *stTask = new PndSoftTriggerTask(pbarmom,0,0,triggercfg);
        stTask->SetConfigurationFile(selectioncfg);
        
        // set parameters for pi0, KS, eta selection
        stTask->SetPi0SignalParams(0.136, 0.0045);
        stTask->SetEtaSignalParams(0.552, 0.009);
        stTask->SetKs0SignalParams(0.497, 0.008);
        
        //TString algo = "PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts"; // FullSim
        TString algo = "PidChargedProbability"; // Fast Sim
        stTask->SetPidAlgoAll(algo);
        
        // Full selection switch
        // 0: detailed selection turned off
        // 1: exclusive mode (modes w/o detailed cut definitions are rejected)
        // 2: open mode (modes w/o detailed cut definitions just have to fulfill their mass window criterion)
        stTask->ApplyFullSelection(1);
        
        stTask->SetTagAll(true);         // tag all modes
        //stTask->SetTagMode(120,true);  // example to switch single tags on/off; mode number has to match one from config file
        
        stTask->SetQAAll(true);          // ntuple output for all modes off
        stTask->SetQAEvent(true);        // ntuple output for event info
        stTask->SetQAMode(201,true);     // example to switch single QA on/off
        
        stTask->SetGammaMinE(0.15);      // global energy pre-cut for neutrals 
        stTask->SetTrackMinP(0.15);      // global momentum pre-cut for charged         
        stTask->SetInitialPidCut(0.1);   // global PID pre-cut for charged      
        
        stTask->McMatchAllowPhotos(1,0.05);  // MC truth match shall partially ignore photos photons (here max 1 with E<50 MeV)
        
        fRun->AddTask(stTask);

        // --------------------------------
        // *** Analysis Task ***
        // --------------------------------

        PndAnaWithTrigger *anaTask = new PndAnaWithTrigger(pbarmom, OutFile+"_ana.root");
        fRun->AddTask(anaTask);

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