QAmacro_hyp_1.C File Reference

#include "TString.h"

Go to the source code of this file.

Functions
int	QAmacro_hyp_1 (Int_t nEvents=500, TString SimEngine="TGeant4", Bool_t HYP_File=false, Bool_t gam=false)

Function Documentation

int QAmacro_hyp_1	(	Int_t	nEvents = 500,
		TString	SimEngine = "TGeant4",
		Bool_t	HYP_File = false,
		Bool_t	gam = false
	)

Definition at line 21 of file QAmacro_hyp_1.C.

References BeamMomentum, Bool_t, boxGen, Cave, ctime, Double_t, fField, fRun, gam, gDebug, Hyp, kParameterMerged, nEvents, output, primGen, printf(), rtdb, rtime, PndHyp::SetAbsorberVol(), PndHyp::SetGeoVersion(), PndHyp::SetHypGamFEm(), PndHyp::SetHypSDtoFile(), PndHyp::SetListMat(), PndHyp::SetMatbud(), PndHyp::SetSensorVol(), PndHyp::SetTreeFName(), timer, and TString.

{

  //-----User Settings:-----------------------------------------------
  TString  OutputFile     ="hypqasim.root";
  TString  ParOutputfile  ="hypqapar.root";
  TString  MediaFile      ="media_pnd.geo";
  gDebug                  = 0;
  
  Float_t     BeamMomentum   =3.0; // ** change HERE if you run Box generator
  //------------------------------------------------------------------
  TStopwatch timer;
  timer.Start();
  gRandom->SetSeed();
  
  // Create the Simulation run manager--------------------------------
  FairRunSim *fRun = new FairRunSim();
  fRun->SetName(SimEngine.Data() );
  fRun->SetOutputFile(OutputFile.Data());
  //  fRun->SetGenerateRunInfo(kFALSE);
  fRun->SetBeamMom(BeamMomentum);
  fRun->SetMaterials(MediaFile.Data());
  //fRun->SetUseFairLinks(kTRUE);
  FairRuntimeDb *rtdb=fRun->GetRuntimeDb();
  
  
  //---------------------Set Parameter output      ----------
  Bool_t kParameterMerged=kTRUE;
  FairParRootFileIo* output=new FairParRootFileIo(kParameterMerged);
  output->open(ParOutputfile.Data());
  rtdb->setOutput(output);
  
  // Create and add detectors
  
  //-------------------------  CAVE      -----------------
  
  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  fRun->AddModule(Cave);
  
  //-------------------------  HYP     -----------------
  PndHyp *Hyp = new PndHyp("HYP",kTRUE);
  Hyp->SetAbsorberVol("Absorber"); // absorber layer
    Hyp->SetSensorVol("Sensor");   // silicon sensor
    Hyp->SetGeoVersion("List");
    Hyp->SetListMat("HYPboron");
    Hyp->SetListMat("HYPaluminium");
    Hyp->SetListMat("Al97Mg3");
    Hyp->SetListMat("titanium");
    Hyp->SetListMat("HYPcarbon");
    Hyp->SetListMat("siliconinactive");
    Hyp->SetGeometryFileName("TargetSystem_WindowAbsB_insideAbsB.root");
  
    if(HYP_File){
      
      Hyp->SetHypSDtoFile(true,HYP_File);
      
      if(gam){
        Hyp->SetHypGamFEm(gam);
        Hyp->SetTreeFName("hypBupV1T_Decay_gam.root");
      }else{
        Hyp->SetTreeFName("hypBupV1T_Decay.root");
      }
    }
  
    Hyp->SetMatbud(false); //material budget hit Collection production
    
    fRun->AddModule(Hyp);
    
    //gROOT->LoadMacro("$VMCWORKDIR/gconfig/SetFragments.C");
    //FragConfig(fRun);
    
    //fRun->SetUserDecay(kTRUE);
    
    FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
    fRun->SetGenerator(primGen);
    
    FairBoxGenerator* boxGen = new FairBoxGenerator(3312, 1); // 13 = muon; 1 = multipl.
    boxGen->SetPRange(0.1,0.5); // GeV/c
    boxGen->SetPhiRange(0., 360.); // Azimuth angle range [degree]
    boxGen->SetThetaRange(5., 150.); // Polar angle in lab system range [degree]
    boxGen->SetXYZ(0., 0., -55.); // cm
    primGen->AddGenerator(boxGen);

   
  
    PndMultiField *fField= new PndMultiField("FULL");
    fRun->SetField(fField);
    
    // EMC Hit producer
    //-------------------------------
    
    //-------------------------  Initialize the RUN  -----------------
    fRun->Init();
    //-------------------------  Run the Simulation  -----------------
    fRun->Run(nEvents);
    //-------------------------  Save the parameters -----------------
    rtdb->saveOutput();
    //------------------------Print some info and exit----------------
    timer.Stop();
    Double_t rtime = timer.RealTime();
    Double_t ctime = timer.CpuTime();
    printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
    
    cout << " Test passed" << endl;
    cout << " All ok " << endl;
    delete fRun;  
  return 0;
  
};