void KRATProcess(){
     
       // Input files
//pc        TString inFile1="/data.local2/mheil/fairroot_results/comp_RPC/r3bsim_2x2x2.5m_35m_131Sn_1n_1000AMeV_100keV.root";
        TString inFile1="r3bsim.root";
     //   TString inFile1      = "r3bsim.root";
//        TString parFile="pippo.root";
     //   TString parFile     = "/data.local2/mheil/fairroot_results/land/r3bpar_10m_land_1n_170MeV.root";
        TString parFile     = "r3bpar.root";
       // Output file	      
//pc       TString outFile ="/data.local2/mheil/fairroot_results/comp_RPC/land_digi_2x2x2.5m_35m_131Sn_1n_1000AMeV_100keV.root";
       TString outFile ="output.root";
     //  TString outFile     = "land_digi.root";
     

       // ----  Load libraries   -------------------------------------------------
       gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
       basiclibs();
     
       gSystem->Load("libGenVector");
     
       gSystem->Load("libGeoBase");
       gSystem->Load("libParBase");
     
     
       gSystem->Load("libBase");
       gSystem->Load("libMCStack");
       gSystem->Load("libField");
       gSystem->Load("libGen");
     
       //----  Load R3B specific libraries ---------------------------------------
       gSystem->Load("libR3Bbase");
       gSystem->Load("libR3BGen");
       gSystem->Load("libR3BPassive");
       gSystem->Load("libR3BData");
       gSystem->Load("libR3BLand");
       gSystem->Load("libR3BChimera");
       
       gSystem->Load("libData");
       gSystem->Load("libKratta"); 
       gSystem->Load("libKrattaTasks"); 
     
     
     
     
       // -----   Timer   --------------------------------------------------------
       TStopwatch timer;
       timer.Start();
       // ------------------------------------------------------------------------
     
       // -----   Digitization run   -------------------------------------------
       FairRunAna *fRun= new FairRunAna();
       fRun->SetInputFile(inFile1);
      // fRun->SetInputFile(parFile);
     //  fRun->AddFile(inFile2);
     //  fRun->AddFile(inFile3);
     //  fRun->AddFile(inFile4);
       fRun->SetOutputFile(outFile);
          TFile* file = new TFile("r3bpar.root");
 	        file->Get("FairBaseParSet");
 	        file->Get("FairGeoNode");
 	        file->Get("FairGeoVolume");
 	    
 TGeoManager* gGeoManager = (TGeoManager*) file->Get("FAIRGeom");
     
       // Verbosity Mode level
       // (0=quiet, 1=event level, 2=track level, 3=debug)
       Int_t iVerbose = 1;
       //Connect the Digitization  Task
       KRATTATreat* tKrattaTreat  = new KRATTATreat();
       //Connect the Tracking Task
//pc       R3BNeutronTracker* tracker  = new R3BNeutronTracker();
     
       // -----   beam energy per nucleon  --------------------------------------
     
     //  Double_t beamEnergy=200.;
     //  Double_t beamBeta=0.5676881;
     
     //  Double_t beamEnergy=600.;
     //  Double_t beamBeta=0.7937626;
     
//pc       Double_t beamEnergy=1000.;
//pc       Double_t beamBeta=0.8760237;
     
     
     
//pc       ((R3BNeutronTracker*) tracker)->UseBeam(beamEnergy,beamBeta);
     
       fRun->AddTask(tKrattaTreat);
//pc       fRun->AddTask(tracker);
     
       // Runtime DataBase info
       FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
/*pc       FairParRootFileIo*  parIo1 = new FairParRootFileIo();
       parIo1->open(parFile.Data());
       rtdb->setFirstInput(parIo1);
       rtdb->setOutput(parIo1);
       rtdb->saveOutput();
     pc*/
       // Load the Root Geometry
//       fRun->LoadGeometry();

       // Number of events to process
       //Int_t nEvents = 2500;

       // -----   Intialise and run   --------------------------------------------
       fRun->Init();
       //fRun->Run(0, nEvents);
       fRun->Run();
       
       // -----   Finish   -------------------------------------------------------
       timer.Stop();
       Double_t rtime = timer.RealTime();
       Double_t ctime = timer.CpuTime();
       cout << endl << endl;
       cout << "Macro finished succesfully." << endl;
       cout << "Output file writen:  "    << outFile << endl;
       cout << "Parameter file writen " << parFile << endl;
       cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
       cout << endl;
       // ------------------------------------------------------------------------
     
     
     }