// Config file test for MUON spectormeter // Remember to define the directory and option // gAlice->SetConfigFunction("Config('$HOME','box');"); void Config(char directory[100]="", char option[6]="param") { //===================================================================== // Config file for MUON test //===================================================================== // Libraries required by geant321 gSystem->Load("libgeant321.so"); new TGeant3TGeo("C++ Interface to Geant3"); //======================================================================= // Create the output file Text_t filename[100]; sprintf(filename,"%sgalice.root",directory); cout << ">>> Output file is " << filename << endl; cout << ">>> Config.C: Creating Run Loader ..."<Fatal("Config.C","Can not instatiate the Run Loader"); return; } rl->SetCompressionLevel(2); rl->SetNumberOfEventsPerFile(100); gAlice->SetRunLoader(rl); //======================================================================= // For having more debuging messages //AliLog::SetModuleDebugLevel("MUON", 1); //======================================================================= // Set External decayer TVirtualMCDecayer *decayer = new AliDecayerPythia(); decayer->SetForceDecay(kAll); decayer->Init(); gMC->SetExternalDecayer(decayer); //======================================================================= // ******* GEANT STEERING parameters FOR ALICE SIMULATION ******* gMC->SetProcess("DCAY",1); gMC->SetProcess("PAIR",1); gMC->SetProcess("COMP",1); gMC->SetProcess("PHOT",1); gMC->SetProcess("PFIS",0); gMC->SetProcess("DRAY",0); gMC->SetProcess("ANNI",1); gMC->SetProcess("BREM",1); gMC->SetProcess("MUNU",1); gMC->SetProcess("CKOV",1); gMC->SetProcess("HADR",1); gMC->SetProcess("LOSS",2); gMC->SetProcess("MULS",1); gMC->SetProcess("RAYL",1); Float_t cut = 1.e-3; // 1MeV cut by default Float_t tofmax = 1.e10; gMC->SetCut("CUTGAM", cut); gMC->SetCut("CUTELE", cut); gMC->SetCut("CUTNEU", cut); gMC->SetCut("CUTHAD", cut); gMC->SetCut("CUTMUO", cut); gMC->SetCut("BCUTE", cut); gMC->SetCut("BCUTM", cut); gMC->SetCut("DCUTE", cut); gMC->SetCut("DCUTM", cut); gMC->SetCut("PPCUTM", cut); gMC->SetCut("TOFMAX", tofmax); //======================================================================= // Examples of generators. Only option param is sistematically tested if (!strcmp(option,"box")) { AliGenBox * gener = new AliGenBox(1); gener->SetMomentumRange(20.,20.1); gener->SetPhiRange(0., 360.); gener->SetThetaRange(171.000,178.001); gener->SetPart(13); // Muons gener->SetOrigin(0.,0., 0.); //vertex position gener->SetSigma(0.0, 0.0, 0.0); //Sigma in (X,Y,Z) (cm) on IP position } if (!strcmp(option,"gun")) { AliGenFixed *gener = new AliGenFixed(ntracks); gener->SetMomentum(10); gener->SetPhiRange(0.); gener->SetThetaRange(0.); gener->SetOrigin(30,30,1200);//vertex position gener->SetPart(13); //GEANT particle type 13 is muons } if (!strcmp(option,"scan")) { AliGenScan *gener = new AliGenScan(-1); gener->SetMomentumRange(10,10); gener->SetPhiRange(0, 0); gener->SetThetaRange(-180, -180); //vertex position //gener->SetSigma(1,1,0); //Sigma in (X,Y,Z) (cm) on IP position gener->SetPart(kRootino); gener->SetRange(100, -300., 300., 100, -300., 300., 1, 2000, 2000); } if (!strcmp(option,"param")) { AliGenParam *gener = new AliGenParam(1, AliGenMUONlib::kUpsilon); gener->SetMomentumRange(0,999); gener->SetPtRange(0,100.); gener->SetPhiRange(0., 360.); gener->SetCutOnChild(1); gener->SetChildPhiRange(0.,360.); gener->SetChildThetaRange(171.0,178.0); gener->SetOrigin(0,0,0); //vertex position gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position gener->SetForceDecay(kDiMuon); gener->SetTrackingFlag(1); gener->Init(); } if (!strcmp(option,"hijing")) { //Hijing generator from ConfigPPR in macros AliGenHijing *gener = new AliGenHijing(-1); // centre of mass energy gener->SetEnergyCMS(5500.); // reference frame gener->SetReferenceFrame("CMS"); // projectile gener->SetProjectile("A", 208, 82); gener->SetTarget ("A", 208, 82); // tell hijing to keep the full parent child chain gener->KeepFullEvent(); // enable jet quenching gener->SetJetQuenching(1); // enable shadowing gener->SetShadowing(1); // neutral pion and heavy particle decays switched off gener->SetDecaysOff(1); // Don't track spectators gener->SetSpectators(0); // kinematic selection gener->SetSelectAll(0); // impact parameter range gener->SetImpactParameterRange(0., 5.); // 0. - 5. fm corresponds to ~10% most central gener->Init(); } if (!strcmp(option,"muoncocktail")) { // Muon cocktail for PbPb AliGenMUONCocktail * gener = new AliGenMUONCocktail(); gener->SetPtRange(1.,100.); // Transverse momentum range gener->SetPhiRange(0.,360.); // Azimuthal angle range gener->SetYRange(-4.0,-2.5); gener->SetMuonPtCut(0.5); gener->SetMuonThetaCut(171.,178.); gener->SetMuonMultiplicity(2); gener->SetImpactParameterRange(0.,5.); // 10% most centra PbPb collisions gener->SetVertexSmear(kPerTrack); gener->SetOrigin(0,0,0); // Vertex position gener->SetSigma(0,0,0.0); // Sigma in (X,Y,Z) (cm) on IP position gener->Init(); } //============================================================= // Field (L3 0.5 T) AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG); gAlice->SetField(field); //============================================================= //=================== Alice BODY parameters ============================= AliBODY *BODY = new AliBODY("BODY","Alice envelop"); //=================== ABSO parameters ============================ AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber"); //=================== DIPO parameters ============================ AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2"); //================== HALL parameters ============================ AliHALL *HALL = new AliHALL("HALL", "Alice Hall"); //=================== PIPE parameters ============================ AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe"); //=================== SHIL parameters ============================ AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2"); //=================== MUON Subsystem =========================== AliMUON *MUON = new AliMUONv1("MUON", "default"); // The 3 switches below are to be used for the trigger code // their default value is set in AliMUON.h // activate trigger cluster-size (0=default, 1=cluster-size according to AliMUONResponseTriggerV1 // MUON->SetTriggerResponseV1(0); // activate 4/4 trigger coincidence (0=default (coinc 3/4), 1=coinc 4/4) // MUON->SetTriggerCoinc44(0); // activate trigger chamber efficiency by cells (0=default, 1=trigger efficiency according to AliMUONTriggerEfficiencyCells // MUON->SetTriggerEffCells(0); // To get same as above w/o noise-only digits for the tracker do : // MUON->SetDigitizerWithNoise(kKALSE); // // If SetAlign, the detection elements transformations // are taken from the input file and not from the code // MUON->SetAlign("transform.dat"); // To generate and read scaler trigger events in rawdata // MUON->SetTriggerScalerEvent(); // If you want to play with builders, first reset the geometry builder, // and then add yours. // MUON->ResetGeometryBuilder(); // MUON->AddGeometryBuilder(new AliMUONSt1GeometryBuilderV2(MUON)); // MUON->AddGeometryBuilder(new AliMUONSt2GeometryBuilderV2(MUON)); // MUON->AddGeometryBuilder(new AliMUONSlatGeometryBuilder(MUON)); // MUON->AddGeometryBuilder(new AliMUONTriggerGeometryBuilder(MUON)); } Float_t EtaToTheta(Float_t arg){ return (180./TMath::Pi())*2.*atan(exp(-arg)); }