Resolution_and_acceptance.cxx File Reference

#include <iostream>
#include <sstream>
#include <TApplication.h>
#include <TCanvas.h>
#include <TROOT.h>
#include <TString.h>
#include <TChain.h>
#include <TFile.h>
#include <TClonesArray.h>
#include <TSystem.h>
#include <TH1.h>
#include <TH2.h>
#include <TRotation.h>
#include <TVector3.h>
#include <TMath.h>
#include <TGaxis.h>
#include <TStopwatch.h>
#include <PndMCTrack.h>
#include <PndSdsMCPoint.h>
#include <FairRunAna.h>
#include <FairRootManager.h>
#include <FairGeane.h>
#include <FairRtdbRun.h>
#include <FairRuntimeDb.h>
#include <FairParRootFileIo.h>
#include <PndLmdPerformanceTask.h>
#include <FairLogger.h>
#include <string>

Go to the source code of this file.

Functions
int	main (int nargs, char **args)

Function Documentation

int main	(	int	nargs,
		char **	args
	)

Definition at line 44 of file Resolution_and_acceptance.cxx.

References ctime, DigiFile, Double_t, fRun, Geane, MCFile, outFile, parFile, parInput1, RecoFile, rtdb, rtime, PndLmdPerformanceTask::SetHistFilename(), startEvent, storePath, timer, and TString.

                                 {
        //TApplication myapp("myapp",0,0);
        TString storePath = "./tmpOutput";
        // set the storepath to somewhere else, if provided
        if (nargs == 2) {
                storePath = args[1];
        } else {
                if (nargs == 1)
                        cout << " Hint: you may change the standard path from "
                                        << storePath.Data()
                                        << " to any other folder by providing the path as an argument! "
                                        << endl;
                if (nargs > 2)
                        cout << " Error: too many arguments! " << endl;
        }
        TString startEvent = "0";
        // ========================================================================
        // Input file (MC events)
        TString MCFile = storePath + "/Lumi_MC_";
        MCFile += startEvent;
        MCFile += ".root";
        TString DigiFile = storePath + "/Lumi_digi_";
        DigiFile += startEvent;
        DigiFile += ".root";
        // Digi file
        TString RecoFile = storePath + "/Lumi_reco_";
        RecoFile += startEvent;
        RecoFile += ".root";
        // TCand file
        TString CandFile = storePath + "/Lumi_TCand_";
        CandFile += startEvent;
        CandFile += ".root";
        // Parameter file   << this one
        TString parFile = storePath + "/Lumi_Params_";
        parFile += startEvent;
        parFile += ".root";
        // Track file
        TString TrkFile = storePath + "/Lumi_Track_";
        TrkFile += startEvent;
        TrkFile += ".root";

        // ----  Load libraries   -------------------------------------------------
        gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
        gSystem->Load("libSds");
        gSystem->Load("libSdsReco");
        gSystem->Load("libLmd");
        gSystem->Load("libLmdReco");
        gSystem->Load("libLmdTrk");
        // ------------------------------------------------------------------------
        // Output file
        TString outFile = storePath + "/Lumi_Testing_";
        outFile += startEvent;
        outFile += ".root";

        // Output histogram file
        TString outFilehist = storePath + "/Lumi_performance_histograms_";
        outFilehist += startEvent;
        outFilehist += ".root";

        std::cout << "MCFile  : " << MCFile.Data() << std::endl;
        std::cout << "DigiFile: " << DigiFile.Data() << std::endl;
        std::cout << "RecoFile: " << RecoFile.Data() << std::endl;
        std::cout << "TCandFile: " << CandFile.Data() << std::endl;
        std::cout << "TrackFile: " << TrkFile.Data() << std::endl;
        std::cout << "TestingFile: " << outFile.Data() << std::endl;
        // ---  Now choose concrete engines for the different tasks   -------------
        // ------------------------------------------------------------------------


        // In general, the following parts need not be touched
        // ========================================================================


        // -----   Timer   --------------------------------------------------------
        TStopwatch timer;
        timer.Start();
        // ------------------------------------------------------------------------


        // -----   Reconstruction run   -------------------------------------------
        FairRunAna *fRun = new FairRunAna();
        fRun->SetInputFile(MCFile);
        //fRun->AddFriend(DigiFile);
        //fRun->AddFriend(RecoFile);
        //fRun->AddFriend(CandFile);
        //fRun->AddFriend(TrkFile);
        fRun->SetOutputFile(outFile);
        // ------------------------------------------------------------------------

        FairGeane *Geane = new FairGeane();
        fRun->AddTask(Geane);

        // -----  Parameter database   --------------------------------------------
        FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
        FairParRootFileIo* parInput1 = new FairParRootFileIo(kTRUE);
        parInput1->open(parFile.Data());
        rtdb->setFirstInput(parInput1);

        // =========================================================================
        // ======                  Geane Back-Propagating                     ======
        // =========================================================================

        // -----    SDS hit producer   --------------------------------------------

        // PndEmcMapper *emcMap = PndEmcMapper::Instance(6);
        Double_t fpBeam = 1.5;
        TVector3 IP(0, 0, 0);
        PndLmdPerformanceTask* lmdperformance = new PndLmdPerformanceTask(fpBeam,
                        IP);
        lmdperformance->SetVerbose(0);
        fRun->AddTask(lmdperformance);
        rtdb->setOutput(parInput1);
        lmdperformance->SetHistFilename(outFilehist);
        rtdb->print();
        // =====                        End of Geane                           =====
        // =========================================================================

        // -----   Intialise and run   --------------------------------------------
        fRun->Init();
        //  // PndEmcMapper *emcMap = PndEmcMapper::Instance(6);
        //  PndEmcMapper *emcMap = PndEmcMapper::Instance();
        // //Geane->SetField(fRun->GetField());
        FairRootManager* ioman = FairRootManager::Instance();
        if (!ioman) {
                std::cout << "Error: "
                                << "RootManager not instantiated!" << std::endl;
                return 1;
        }
        fRun->Run(0, ioman->CheckMaxEventNo()); // hmm does not work yet, seems to return 0 but then is runs over all what was intended.
        // ------------------------------------------------------------------------
        rtdb->saveOutput();
        rtdb->print();
        // -----   Finish   -------------------------------------------------------
        timer.Stop();
        Double_t rtime = timer.RealTime();
        Double_t ctime = timer.CpuTime();
        cout << endl << endl;
        cout << "Macro finished succesfully." << endl;
        cout << "Output file is " << outFile << endl;
        cout << "Parameter file is " << parFile << endl;
        cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
        cout << endl;
        // ------------------------------------------------------------------------

        //myapp.Run();
        return 0;
}