class AtlAppAnalysisTask: public AtlTask

Customised base class for A++ tasks

Description:

Tasks provide a folder-like structure for running analysis jobs, drawing plots etc. For a general description see the parent class TTask. The AtlAPPAnalysiTask class is especially designed for running A++ analysis jobs made of analysis selectors inheriting from AtlSelector (like the AtlExampleAnalysis). The common case is that the same analysis is run several times. First over DATA and afterwards over various Monte Carlo samples (signal and background). Often the MC samples are grouped into different processes (ttbar production, W+jets etc.) With the help of the AtlAppAnalysisTask class this probem can be easily solved.

Execution of analysis tasks is possible on a local machine (interactive running), on a batch system or on the GRID via the PANDA distributed analysis system. Details are described below.

Two things are needed for running analysis tasks:

1. A task class inheriting from this class. Here, the implementation of the InitSelector() routine is mandatory. In this routine the analysis selector is created and configured, i.e. analysis cuts can be set (via the routine SetCut() and analysis tools (see class AtlAnalysisTool) can be added and configured as well (via SetToolCut()). The advantage of this method is that the cuts and configurations given here apply likewise to each of the analyses managed by this task. As an example, the analysis task for the single-top t-channel analysis is shown.
   Atl_tChannelTask.C:

    //____________________________________________________________________
    //
    // Single-top t-channel analysis task
    // 
    //  
    // Author: Oliver Maria Kind 
    // Update: $Id: AtlAppAnalysisTask.cxx,v 1.77 2017/08/08 11:32:07 kaphle Exp $
    // Copyright: 2009 (C) Oliver Maria Kind
    //
    #ifndef SINGLETOP_Atl_tChannelTask
    #define SINGLETOP_Atl_tChannelTask
    #ifndef ATLAS_AtlAppAnalysisTask
    #include 
    #endif
    #ifndef ATLAS_AtlOverlapRemoval
    #include 
    #endif
    #ifndef ATLAS_AtlSgTop_tChannelFinder
    #include 
    #endif
    #include 
    class Atl_tChannelTask : public AtlAppAnalysisTask {
       
    public:
        Atl_tChannelTask(const char* name, const char* title);
        virtual ~Atl_tChannelTask();
        virtual void InitSelector();
       
        ClassDef(Atl_tChannelTask,0) // Single-top t-channel analysis task
    };
    #endif
    #ifndef __CINT__
    ClassImp(Atl_tChannelTask);
    #endif
    //____________________________________________________________________
    Atl_tChannelTask::Atl_tChannelTask(const char* name,
    				   const char* title) :
        AtlAppAnalysisTask(name, title) {
        //
        // Normal constructor
        //
        SetInputTreeName("t_app");
        AddUserEnv("LIBSINGLETOP");
    }
    //____________________________________________________________________
    Atl_tChannelTask::~Atl_tChannelTask() {
        //
        // Default destructor
        //
    }
    //____________________________________________________________________
    void Atl_tChannelTask::InitSelector() {
        //
        // Initialise the t-channel analysis selector
        //
        // Create single-top t-channel analysis selector
        fSelector = new Atl_tChannelAnalysis(fOutputFileName->Data());
        // ============================    
        // Definition of selection cuts  
        // ============================
     
        //
        // Leptons
        // -------
        // Electrons
        SetCut("fSignalElectronAuthor",        "AtlEMShower::kHighPtElectron");// | AtlEMShower::kSoftElectron");
        SetCut("fSignalElectronIsEM",          "AtlEMShower::kElectronTight");
        SetCut("fSignalElectron_Pt_min",       "10.");
        SetCut("fSignalElectron_Pt_max",       "10.e10");
        SetCut("fSignalElectron_Eta_min",      "-2.5");
        SetCut("fSignalElectron_Eta_max",      "2.5");
        SetCut("fSignalElectron_EtCone20_max", "6.");
        ...
        // Add overlap removal procedure
        AtlOverlapRemoval *overlap = new AtlOverlapRemoval("overlap_removal",
     						       "Overlap Removal Single-Top t-Channel");
        SetToolCut("overlap_removal","fFindWorkingPoint", "kFALSE");
        SetToolCut("overlap_removal","fElDeltaRMax",      "0.2");
        ...
        fSelector->AddTool(overlap);
       
        // Add single-top t-channel kinematic fitter
        AtlSgTop_tChannelFinder *kinfit =
    	new AtlSgTop_tChannelFinder("sgtop_finder_kinfitter",
    	    "Single-Top t-Channel Reconstruction by KinFitter");
        SetToolCut("sgtop_finder_kinfitter", "fMode", "AtlSgTop_tChannelFinder::kKinFit");
        ...
        fSelector->AddTool(kinfit);
    }
    

2. A task script which defines the task folder structure. Note that such a script can contain different taks, e.g. analysis and plotting. For example the script corresponding to the single-top task above, given below.
   task_tChannelAnalysis.C:

    // {
        gROOT->LoadMacro("Atl_tChannelTask.C+");
        // ===============================  
        // Luminosities of the MC samples
        // ===============================
        //-----------
        // single top
        //-----------
        Float_t Lumi_singletop_tchannel_enu = (19750/14.46);
        Float_t Lumi_singletop_tchannel_munu = (19999/14.46);
        Float_t Lumi_singletop_tchannel_taunu = (19500/14.46) ;
        ...
       
        // Main task
        TTask *tChannel = new TTask("tChannel", "Single-top t-Channel Analysis");
        //////////////////////////////////////////////////////////////////////
        //
        // Single-top production
        //
        //////////////////////////////////////////////////////////////////////
        TTask *sgtop = new TTask("single-top", "Single-top production");
        tChannel->Add(sgtop);
     
        // =========
        // t-channel
        // =========
        TTask *sgtop_tchannel = new TTask("t_channel", "t-channel");
        sgtop->Add(sgtop_tchannel);
        // e-nu
        Atl_tChannelTask *sgtop_tchannel_enu
    	= new Atl_tChannelTask("sgtop_tchannel_enu", "e-nu decay");
        sgtop_tchannel_enu->AddInputFiles("/dats3/atlas/mandry/user10.RoccoMandrysch.mc09_7TeV.108340.st_tchan_enu_McAtNlo_Jimmy.merge.AOD.e508_s765_s767_r1250_r1260.22.05.2010.app.test_2/user10.RoccoMandrysch.mc09_7TeV.108340.st_tchan_enu_McAtNlo_Jimmy.merge.AOD.e508_s765_s767_r1250_r1260.22.05.2010.app.test_2*.root");
        sgtop_tchannel_enu->SetOutputFile("$LIBSINGLETOP/hist/tChannelAnalysis/tchannel_enu.root");
        sgtop_tchannel->Add(sgtop_tchannel_enu);
     
        // mu-nu
        Atl_tChannelTask *sgtop_tchannel_munu
    	= new Atl_tChannelTask("sgtop_tchannel_munu", "mu-nu decay");
        sgtop_tchannel_munu->AddInputFiles("/dats3/atlas/mandry/user10.RoccoMandrysch.mc09_7TeV.108341.st_tchan_munu_McAtNlo_Jimmy.merge.AOD.e508_s765_s767_r1250_r1260.22.05.2010.app.test/user10.RoccoMandrysch.mc09_7TeV.108341.st_tchan_munu_McAtNlo_Jimmy.merge.AOD.e508_s765_s767_r1250_r1260.22.05.2010.app.test*.root");
        sgtop_tchannel_munu->SetOutputFile("$LIBSINGLETOP/hist/tChannelAnalysis/tchannel_munu.root");
        sgtop_tchannel->Add(sgtop_tchannel_munu);
     
        // tau-nu
        Atl_tChannelTask *sgtop_tchannel_taunu
    	= new Atl_tChannelTask("sgtop_tchannel_taunu", "tau-nu decay");
        sgtop_tchannel_taunu->AddInputFiles("/dats3/atlas/mandry/app_20_02_2010/user09.RoccoMandrysch.ganga.mc08.108342.st_tchan_taunu_McAtNlo_Jimmy.merge.AOD.e405_s495_s520_r808_r838.app.test_20.02.2010.WUPPERTALPROD_SCRATCHDISK/user09.RoccoMandrysch.ganga.mc08.108342.st_tchan_taunu_McAtNlo_Jimmy.merge.AOD.e405_s495_s520_r808_r838.app.test_*.root");
        sgtop_tchannel_taunu->AddInputFiles("/dats3/atlas/mandry/user10.RoccoMandrysch.mc09_7TeV.108342.st_tchan_taunu_McAtNlo_Jimmy.merge.AOD.e508_s765_s767_r1250_r1260.22.05.2010.app.test/user10.RoccoMandrysch.mc09_7TeV.108342.st_tchan_taunu_McAtNlo_Jimmy.merge.AOD.e508_s765_s767_r1250_r1260.22.05.2010.app.test*.root");
        sgtop_tchannel_taunu->SetOutputFile("$LIBSINGLETOP/hist/tChannelAnalysis/tchannel_taunu.root");
        sgtop_tchannel->Add(sgtop_tchannel_taunu);
        // =========
        // s-channel
        // =========
        TTask *sgtop_schannel = new TTask("s_channel", "s-channel");
        sgtop->Add(sgtop_schannel);
     
        // e-nu
        Atl_tChannelTask *sgtop_schannel_enu
    	= new Atl_tChannelTask("sgtop_schannel_enu", "e-nu decay");
        sgtop_schannel_enu->AddInputFiles("/dats3/atlas/mandry/user10.RoccoMandrysch.mc09_7TeV.108343.st_schan_enu_McAtNlo_Jimmy.merge.AOD.e534_s765_s767_r1250_r1260.22.05.2010.app.test/user10.RoccoMandrysch.mc09_7TeV.108343.st_schan_enu_McAtNlo_Jimmy.merge.AOD.e534_s765_s767_r1250_r1260.22.05.2010.app.test*.root");
        sgtop_schannel_enu->SetOutputFile("$LIBSINGLETOP/hist/tChannelAnalysis/schannel_enu.root");
        sgtop_schannel->Add(sgtop_schannel_enu);
        ...  
       
        //////////////////////////////////////////////////////////////////////
        //
        // W + jets production
        //
        //////////////////////////////////////////////////////////////////////
        TTask *WplusJets = new TTask("WplusJets", "W + Jets production");
        tChannel->Add(WplusJets);
     
        // ============
        // e-nu channel
        // ============
        TTask *WplusJets_enu = new TTask("WplusJets_enu", "W+Jets, enu channel");
        WplusJets->Add(WplusJets_enu);
     
        // Np0
        Atl_tChannelTask *WplusJets_enu_Np0
    	= new Atl_tChannelTask("WplusJets_enu_Np0", "W+Jets, enu channel, Np0");
        WplusJets_enu_Np0->AddInputFiles("/dats5/atlas/mandry/user10.RoccoMandrysch.mc09_7TeV.107680.AlpgenJimmyWenuNp0_pt20.merge.AOD.e511_s765_s767_r1250_r1260.22.05.2010.app.test_2/user10.RoccoMandrysch.mc09_7TeV.107680.AlpgenJimmyWenuNp0_pt20.merge.AOD.e511_s765_s767_r1250_r1260.22.05.2010.app.test_2*.root");   
        WplusJets_enu_Np0->SetOutputFile("$LIBSINGLETOP/hist/tChannelAnalysis/Wjets_enu_np0.root");
        WplusJets_enu->Add(WplusJets_enu_Np0);
        ...
        //////////////////////////////////////////////////////////////////////
        //
        // DATA
        //
        //////////////////////////////////////////////////////////////////////
        // L1Calo Trigger Stream (new)
        Atl_tChannelTask *data_l1calo = new Atl_tChannelTask("data-l1-calo", "Data from L1Calo trigger stream");
        data_l1calo->AddInputFiles("/dats2/atlas/mandry/user10.RoccoMandrysch.data10_7TeV.00154813.physics_L1Calo.merge.AOD.f255_m471.22.05.2010.app.test/user10.RoccoMandrysch.data10_7TeV.00154813.physics_L1Calo.merge.AOD.f255_m471.22.05.2010.app.test*.root");
        data_l1calo->SetOutputFile("$LIBSINGLETOP/hist/tChannelAnalysis/data_l1calo.root");
        tChannel->Add(data_l1calo);
        ...  
        //////////////////////////////////////////////////////////////////////
        //
        // Draw all MC-DATA plots
        //
        //////////////////////////////////////////////////////////////////////
        HepDataMCPlotter *pl_sgtop = new HepDataMCPlotter("DATA-MC plotter", "Create all Data-MC plots");
        pl_sgtop->SetWorkingDir("$LIBSINGLETOP/hist/tChannelAnalysis/plots_datamc/");
        pl_sgtop->SetLumiDATA(0.003260);
        // Enable/disable grouping of histograms
        Bool_t group_histograms = kTRUE;
        // ==========
        // Single-top
        // ==========
        //Set grouping of histograms
        pl_sgtop->SetGroupHistograms(group_histograms);
        pl_sgtop->AddMCFolder("sgtop", "sgtop incl", kRed);
        pl_sgtop->AddMCSample("sgtop","$LIBSINGLETOP/hist/tChannelAnalysis/tchannel_enu.root",
                        "MC single-top t-channel e#nu", Lumi_singletop_tchannel_enu, kRed);
        pl_sgtop->AddMCSample("sgtop","$LIBSINGLETOP/hist/tChannelAnalysis/tchannel_munu.root",
                        "MC single-top t-channel #mu#nu", Lumi_singletop_tchannel_munu, kRed-4);
        pl_sgtop->AddMCSample("sgtop","$LIBSINGLETOP/hist/tChannelAnalysis/tchannel_taunu.root",
                        "MC single-top t-channel #tau#nu", Lumi_singletop_tchannel_taunu, kRed-7);
        ...
        tChannel->Add(pl_sgtop);
        // Open task folder in Root's TObject browser
        gROOT->GetListOfTasks()->Add(tChannel);
        gROOT->GetListOfBrowsables()->Add(tChannel);
        new TBrowser;
    }
    

    Atl_tChannelTask::Atl_tChannelTask(const char* name,
                                    const char* title) :
        AtlAppAnalysisTask(name, title) {
        //
        // Normal constructor
        //
        SetInputTreeName("t_app");
        AddUserEnv("LIBSINGLETOP");
    }
 

 Author: Oliver Maria Kind <mailto: kind@mail.desy.de>
 Update: $Id: AtlAppAnalysisTask.cxx,v 1.77 2017/08/08 11:32:07 kaphle Exp $
 Copyright: 2009 (C) Oliver Maria Kind