sim_hypGe_extparameters_urqmd.C

Go to the documentation of this file.

// Macro created 5/04/2007 by A.Sanchez
// It creates a geant simulation file for hypGe
int sim_hypGe_extparameters_urqmd(Int_t nEvents, Int_t WhichDetector,Int_t iurqmdFile,Int_t StartEvent,Int_t SecTarget){
        // Load basic libraries
  // If it does not work,  please check the path of the libs and put it         by hands
        gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
  gSystem->Load("libHypGe");  
  gSystem->Load("librazhyp");
  gSystem->Load("libHyp");
  gROOT->LoadMacro("$VMCWORKDIR/macro/hypGe/Marcell/SharedMacros/SharedMacroFunctions.C");
  FairRunSim *fRun = new FairRunSim();
  gRandom->SetSeed();
  TStopwatch timer;
  timer.Start();
  gDebug=0;
  TString urqmdFile = "pbarC_";
        urqmdFile += iurqmdFile;
  TString inFile="/home/steinen/work/urqmd/pbarC_3GeV/";
        inFile += urqmdFile;
        inFile += ".root";
        //Choose geometry
        TString outFile="$SIMDATADIR/Neutron/Geo";                                                      // If no SIMDATADIR, same folder as the macro
        TString GeoFile= GeoFileChooser(WhichDetector);
        outFile += WhichDetector;
        //compose the name of the output file and the simparams file
                outFile += "urqmd_";
                outFile += urqmdFile;
                outFile +="_";
                outFile += nEvents;
                outFile += "Evts_StartEvent_";
                outFile += StartEvent;
                if(SecTarget)
                        outFile += "_SecTar";

        
        TString SimparamsFile;
        SimparamsFile=outFile;
        outFile +=".root";
        SimparamsFile += "__Simparams.root";

        //TString inFile="/d/panda02/urqmd_smm/pbarC_3_GeV.root";
  //TString inFile="/u/asanchez/razhyp_gt12.dat";

        //TString inFile="/data/work/kpha1/bleser/GiBUU-data2b/jobs/carbon0003/Conversion_Class_apr13_carbon0003_allParticles/carbon0003_allParticles_1_32.root"; //GiBUU file 1.5*10⁷ events
  // set the MC version used
  // ------------------------
  
  fRun->SetName("TGeant4");

  fRun->SetOutputFile(outFile);
  // Set Material file Name
  //-----------------------
  
  fRun->SetMaterials("media_pnd_hypGe.geo");
  
  // Create and add detectors
  //-------------------------

  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  fRun->AddModule(Cave); 
        
        if(SecTarget)
        {
                //acc sec. target
                PndHyp *Hyp= new PndHyp("HYP",kFALSE);
                        Hyp->SetAbsorberVol("Absorber"); // absorber layer
                        Hyp->SetSensorVol("Sensor");   // silicon sensor
                        Hyp->SetGeometryFileName("../macro/hyp/Sebastian/SekTarget_open_varAbs4Si5_3Q_HYPbe_1mm_MVD.root");
                        fRun->AddModule(Hyp); 
                
                //add MVD outer barrel
                FairDetector *Mvd = new PndMvdDetector("MVD", kFALSE);
                        Mvd->SetGeometryFileName("../macro/hyp/Sebastian/Mvd-2.2_Simplified_onlyStrip5_z-verschoben550.root");  
                        fRun->AddModule(Mvd);
        }
  PndHypGe *HypGe= new PndHypGe("HYPGE",kTRUE);
                TString nam = gSystem->Getenv("VMCWORKDIR");
                HypGe->SetPathGeo(nam.Data());
                HypGe->SetGeometryFileName(GeoFile);
                fRun->AddModule(HypGe); 

  //  FairDetector *Hyp= new FairHyp("HYP",kTRUE);
  //    Hyp->SetGeometryFileName("HypST_prueba24pipe.geo");
  //    fRun->AddModule(Hyp); 

  //    FairDetector *Drc = new FairDrc("DIRC", kTRUE);
  //    Drc->SetGeometryFileName("dirc.geo");
  //    fRun->AddModule(Drc);

  //  FairModule *Magnet= new FairMagnet("MAGNET");
  //   Magnet->SetGeometryFileName("magnet.geo");
  //   fRun->AddModule(Magnet);

  //  FairDetector *Emc = new FairEmc("EMC",kTRUE);
  //   Emc->SetGeometryFileName("emc_module123.dat"); 
  //   fRun->AddModule(Emc);
  
  // Create and Set Event Generator
  //-------------------------------
  
  FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
  fRun->SetGenerator(primGen);
  
        PndUrqmdSmmGenerator* UrqmdGen = new PndUrqmdSmmGenerator(inFile,StartEvent);
                primGen->SetTarget(-55.,0.);
                primGen->AddGenerator(UrqmdGen);
        

  fRun->SetStoreTraj(kTRUE); // to store particle trajectories 


  //magnetic field: no field when commented put
  PndConstField *fMagField=new PndConstField();
  fMagField->SetField(0.,0.,10.); // values are in kG
  fMagField->SetFieldRegion(-50, 50,-50, 50, -100, 100);
   //values are in cm (xmin,xmax,ymin,ymax,zmin,zmax)
  fRun->SetField(fMagField);
  
  


   fRun->Init();
  
  // Fill the Parameter containers for this run
  //-------------------------------------------
  
   FairRuntimeDb *rtdb=fRun->GetRuntimeDb();
   Bool_t kParameterMerged=kTRUE;
  
   
   FairParRootFileIo* output=new FairParRootFileIo(kParameterMerged);
   output->open(SimparamsFile);
   rtdb->setOutput(output);
   rtdb->saveOutput();
   rtdb->print();


  // Transport nEvents
  // -----------------
  
  fRun->Run(nEvents);
  
  timer.Stop();
  
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
  
  return 0;
}