void create_par_draw_hist()
{
    string outputFilePath = "rich700n_june2019p1_params_mixwls_final.txt";

    // This we need to init digi parameters to draw PMTs
    TString asciiParFile = "";
    TString rootParFile = "/cvmfs/hades.gsi.de/param/sim/mar19/allParam_mar19_sim_run_16000_gen2_22102019.root";
    TString paramSource = "root"; // root, ascii, oracle
    TString paramRelease = "MAR19_v2";

    // create loop obejct and hades
    HLoop loop(kTRUE);

    Int_t refId = -1;
    TString beamtime = "mar19";

    HDst::setupSpectrometer(beamtime,NULL,"start,rich,mdc,tof,rpc,shower,wall");
    HDst::setupParameterSources(paramSource,asciiParFile,rootParFile,paramRelease);


    loop.addMultFiles("/lustre/cbm/users/slebedev/hades/data/digipar/2020.rich_params.dst.e.3.root");
    if(!loop.setInput("-*,+HGeantKine,+HRichHitSim")) {
       cout << "READBACK: ERROR : cannot read input !" << endl;
       exit(1);
    }

    HTaskSet *masterTaskSet = gHades->getTaskSet("all");
    HRich700DigiParCreator *creator = new HRich700DigiParCreator();
    creator->setOutputFilePath(outputFilePath);
    masterTaskSet->add(creator);
    cout << "begin loop.nextEvent(0);" << endl;
    loop.nextEvent(0);
    cout << "end loop.nextEvent(0);" << endl;


    string fileName = "richpar_hist_final.root";
    string outputDir = "results_digipar_final/";

    creator->drawFromFile(fileName, outputDir);
}