123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566 |
- //
- // Created by peter on 8/27/19.
- //
- #include <iostream>
- #include <TStopwatch.h>
- #include <TFile.h>
- #include <TH1F.h>
- #include "CentralityCalculator.h"
- int main(int argc, char** argv)
- {
- TString inFileName, CanvasName;
- if (argc < 5)
- {
- std::cerr << "./calcCent -i INPUTFILE -canv CANVASNAME" << std::endl;
- return 10;
- }
- for (int i = 1; i < argc; i++)
- {
- if (std::string(argv[i]) != "-i" &&
- std::string(argv[i]) != "-canv")
- {
- std::cerr << "\ncalcCent: Unknown parameter: " << argv[i] << std::endl;
- return 11;
- } else {
- if (std::string(argv[i]) == "-i" && i != argc - 1)
- {
- inFileName = argv[++i];
- }
- if (std::string(argv[i]) == "-i" && i == argc - 1)
- {
- std::cerr << "\ncalcCent: Input file name was not specified!" << std::endl;
- return 20;
- }
- if (std::string(argv[i]) == "-canv" && i != argc - 1)
- {
- CanvasName = argv[++i];
- }
- if (std::string(argv[i]) == "-canv" && i == argc - 1)
- {
- std::cerr << "\ncalcCent: Input file name was not specified!" << std::endl;
- return 20;
- }
- }
- }
- TStopwatch timer;
- timer.Start();
- TFile *fi = new TFile(inFileName.Data(),"read");
- TH1F *hist = (TH1F*) fi->Get("hMult");
- CentralityCalculator *cent = new CentralityCalculator();
- cent->SetEnergy(200.);
- cent->SetModel("PHSD");
- cent->SetHistogram(hist);
- cent->Calc();
- cent->DrawHistogram(CanvasName);
- timer.Stop();
- timer.Print();
- return 0;
- }
|