runLumiPixel1Digi.C File Reference

Go to the source code of this file.

Functions
int	runLumiPixel1Digi (const int nEvents=10, const int startEvent=0, TString storePath="tmpOutput", std::string misalignment_matrices_path="geo/misalignMatrices-SensorsOnly-100.root", bool use_point_transform_misalignment=false, const int verboseLevel=0, const int pitch=1)

Function Documentation

int runLumiPixel1Digi	(	const int	nEvents = 10,
		const int	startEvent = 0,
		TString	storePath = "tmpOutput",
		std::string	misalignment_matrices_path = "geo/misalignMatrices-SensorsOnly-100.root",
		bool	use_point_transform_misalignment = false,
		const int	verboseLevel = 0,
		const int	pitch = 1
	)

Definition at line 1 of file runLumiPixel1Digi.C.

References ctime, digiparFile, Double_t, fRun, geoPar, gGeoManager, inFile, nEvents, outFile, parFile, parInput1, parInput2, rtdb, rtime, startEvent, storePath, timer, TString, and verboseLevel.

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

        // Input file (MC events)
        TString inFile = storePath + "/Lumi_MC_";
        inFile += startEvent;
        inFile += ".root";
        // Parameter file
        TString parFile = storePath + "/Lumi_Params_";
        parFile += startEvent;
        parFile += ".root";

        // Parameter file
        TString digiparFile = "lumi.digi.par";

        // Output file
        TString outFile = storePath + "/Lumi_digi_";
        outFile += startEvent;
        outFile += ".root";
        std::cout << "DigiFileName: " << outFile.Data() << std::endl;

        // -----   Reconstruction run   -------------------------------------------
        FairRunAna *fRun = new FairRunAna();

        FairFileSource *input_source = new FairFileSource(inFile);
        // input_source->SetEventMeanTime(50);//[ns] 50= interaction rate 20 MHz;
        input_source->SetEventMeanTime(80);  //[ns] 80= interaction rate 12.5 MHz;
        // example: cross-section(DPM el+inlel, 1.5 GeV/c)= 125 mb, L=10^32 -> interaction rate 12.5 MHz
        // input_source->SetEventMeanTime(400);//[ns] 400= interaction rate 5 MHz; //15 GeV
        fRun->SetSource(input_source);
        fRun->SetOutputFile(outFile);

        // -----  Parameter database   --------------------------------------------
        FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
        FairParRootFileIo* parInput1 = new FairParRootFileIo(kTRUE);
        parInput1->open(parFile.Data(), "UPDATE");
        rtdb->setFirstInput(parInput1);
        FairParAsciiFileIo* parInput2 = new FairParAsciiFileIo();
        parInput2->open(digiparFile.Data(), "in");
        rtdb->setSecondInput(parInput2);

        PndLmdHybridHitProducer* lmdStripProd = new PndLmdHybridHitProducer();
        //  PndLmdStripHitProducerDif* ssdStripProd = new PndLmdStripHitProducerDif(); //TEST
        lmdStripProd->SetVerbose(verboseLevel);
        fRun->AddTask(lmdStripProd);

        rtdb->setOutput(parInput1);
        rtdb->print();
        // =====                 End of HitProducers                           =====
        // =========================================================================
        PndSdsGeoPar* geoPar = (PndSdsGeoPar*) (rtdb->getContainer("PndSdsGeoPar"));

        // MialignmentHandler sits here

        if (misalignment_matrices_path != "" && !use_point_transform_misalignment) {
                //load matrices
                cout << "reading misalignMatrices from: " << misalignment_matrices_path << "\n";
                TFile *misalignmentMatrixRootfile = new TFile(misalignment_matrices_path.c_str(), "READ");

                if (misalignmentMatrixRootfile->IsOpen()) {
                        std::map < std::string, TGeoHMatrix > *matrices;

                        gDirectory->GetObject("PndLmdMisalignMatrices", matrices);
                        misalignmentMatrixRootfile->Close();

                        cout << matrices->size() << " matrices successfully read from file.\n";

                        //this call has to be made before fRun->Init();
                        fRun->AddAlignmentMatrices(*matrices);
                }
                else {
                        cout << "file could not be read\n";
                        return 1;
                }
                cout << "starting digi macro\n";
        }

        //we need to Init so the gGeoManager gets populated!
        fRun->Init();

        // -----   Intialise and run   --------------------------------------------
        //fRun->Init();
        fRun->Run(0, nEvents);
        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\a\a" << endl;
        cout << endl;

        // temporary fix to avoid double frees at the destruction of te program for pandaroot/fairroot with root6
        gGeoManager->GetListOfVolumes()->Delete();
        gGeoManager->GetListOfShapes()->Delete();
        delete gGeoManager;

        return 0;
}