// -------------------------------------------------------------------------- // // Macro for standard transport simulation using UrQMD input and GEANT3 // CBM setup with STS only // // V. Friese 14/09/2006 // // -------------------------------------------------------------------------- #include "TRandom3.h" void sts_sim(Int_t nEvents = 1) { // ======================================================================== // Adjust this part according to your requirements // Input file TString inPath = "/d/cbm03/urqmd/auau/25gev/centr/"; TString inFile = inPath + "urqmd.auau.25gev.centr.0000.ftn14"; // Output file TString outFile = Form("data/sts.mc.root",nEvents); // Parameter file TString parFile = Form("data/params.root",nEvents); // Cave geometry TString caveGeom = "cave.geo"; // Target geometry TString targetGeom = "target.geo"; // Beam pipe geometry TString pipeGeom = "pipe_standard.geo"; // Magnet geometry and field map TString magnetGeom = "magnet_active.geo"; TString fieldMap = "FieldActive"; Double_t fieldZ = 50.; // z position of field centre Double_t fieldScale = 1.; // field scaling factor // MVD geometry TString mvdGeom = "mvd_standard.geo"; // STS geometry TString stsGeom = "sts_test_beam_09_1st.geo"; // In general, the following parts need not be touched // ======================================================================== // ---- Debug option ------------------------------------------------- gDebug = 0; // ------------------------------------------------------------------------ // ----- Timer -------------------------------------------------------- TStopwatch timer; timer.Start(); // ------------------------------------------------------------------------ // ---- Load libraries ------------------------------------------------- gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C"); basiclibs(); gSystem->Load("/misc/cbmsoft/cbmdaq/JUN09/sarge32/fileapi/libLmdFile.so"); gSystem->Load("/misc/cbmsoft/cbmdaq/JUN09/sarge32/sw-local/libKnut.so"); gSystem->Load("libGeoBase"); gSystem->Load("libParBase"); gSystem->Load("libBase"); gSystem->Load("libCbmBase"); gSystem->Load("libCbmData"); gSystem->Load("libTestBeam"); gSystem->Load("libField"); gSystem->Load("libGen"); gSystem->Load("libPassive"); gSystem->Load("libMvd"); gSystem->Load("libSts"); // ------------------------------------------------------------------------ // ----- Create simulation run ---------------------------------------- FairRunSim* run = new FairRunSim(); run->SetName("TGeant3"); // Transport engine run->SetOutputFile(outFile); // Output file FairRuntimeDb* rtdb = run->GetRuntimeDb(); // ------------------------------------------------------------------------ // ----- Create media ------------------------------------------------- run->SetMaterials("media.geo"); // Materials // ------------------------------------------------------------------------ // ----- Create geometry ---------------------------------------------- FairModule* cave= new CbmCave("CAVE"); cave->SetGeometryFileName(caveGeom); run->AddModule(cave); /*FairModule* pipe= new CbmPipe("PIPE"); pipe->SetGeometryFileName(pipeGeom); run->AddModule(pipe); FairModule* target= new CbmTarget("Target"); target->SetGeometryFileName(targetGeom); run->AddModule(target); FairModule* magnet= new CbmMagnet("MAGNET"); magnet->SetGeometryFileName(magnetGeom); run->AddModule(magnet); FairDetector* mvd= new CbmMvd("MVD", kTRUE); mvd->SetGeometryFileName(mvdGeom); run->AddModule(mvd);*/ FairDetector* sts= new CbmSts("STS", kTRUE); sts->SetGeometryFileName(stsGeom); run->AddModule(sts); // ------------------------------------------------------------------------ // ----- Create magnetic field ---------------------------------------- /*if ( fieldMap == "FieldActive" || fieldMap == "FieldIron") magField = new CbmFieldMapSym3(fieldMap); else if ( fieldMap == "FieldAlligator" ) magField = new CbmFieldMapSym2(fieldMap); else { cout << "===> ERROR: Field map " << fieldMap << " unknown! " << endl; exit; } magField->SetPosition(0., 0., fieldZ); magField->SetScale(fieldScale); run->SetField(magField);*/ // ------------------------------------------------------------------------ TRandom3* fGen; fGen = new TRandom3(); time_t curtime; time(&curtime); fGen->SetSeed(curtime); Double_t genvert; genvert = fGen->Rndm()*100.; cout<SetThetaRange ( 0., 2.); boxGen->SetPRange (2.3,2.3); boxGen->SetPhiRange (0.,360.); primGen->AddGenerator(boxGen); run->SetGenerator(primGen); /*CbmPrimaryGenerator* primGen = new CbmPrimaryGenerator(); FairIonGenerator *fIongen= new FairIonGenerator(1, 1, 1, 1, 0., 0., 3.5, 0.1, 0.2, 0.); primGen->AddGenerator(fIongen); run->SetGenerator(primGen);*/ // ------------------------------------------------------------------------ run->SetStoreTraj(kTRUE); /* fRun->Init(); CbmTrajFilter *trajFilter = CbmTrajFilter::Instance();// trajFilter->SetStepSizeCut(0.01); // 1 cm // trajFilter->SetStorePrimaries(kTRUE); trajFilter->SetStoreSecondaries(kTRUE); */ // ----- Initialize simulation run ------------------------------------ run->Init(); // ------------------------------------------------------------------------ // ----- Runtime database --------------------------------------------- CbmFieldPar* fieldPar = (CbmFieldPar*) rtdb->getContainer("CbmFieldPar"); //fieldPar->SetParameters(magField); fieldPar->setChanged(); fieldPar->setInputVersion(run->GetRunId(),1); Bool_t kParameterMerged = kTRUE; FairParRootFileIo* parOut = new FairParRootFileIo(kParameterMerged); parOut->open(parFile.Data()); rtdb->setOutput(parOut); rtdb->saveOutput(); rtdb->print(); // ------------------------------------------------------------------------ // ----- Start run ---------------------------------------------------- run->Run(nEvents); // ------------------------------------------------------------------------ // ----- Finish ------------------------------------------------------- timer.Stop(); Double_t rtime = timer.RealTime(); Double_t ctime = timer.CpuTime(); cout << endl << endl; cout << "Macro finished succesfully." << endl; cout << "Output file is " << outFile << endl; cout << "Parameter file is " << parFile << endl; cout << "Real time " << rtime << " s, CPU time " << ctime << "s" << endl << endl; // ------------------------------------------------------------------------ }