MUON/Config_acc.C

   1 enum gentype_t {hijing, gun, box, pythia, param, cocktail, fluka, halo, ntuple, scan, doublescan};
   2
   3 gentype_t gentype=param;
   4 Int_t ntracks=1000;
   5
   6 void Config()
   7 {
   8     gSystem->Load("libgeant321");
   9     new     TGeant3("C++ Interface to Geant3");
  10     TGeant3 *geant3 = (TGeant3*)gMC;
  11
  12 //=======================================================================
  13 //  Create the output file
  14
  15 TFile *rootfile = new TFile("galice.root","recreate");
  16 rootfile->SetCompressionLevel(2);
  17 TGeant3 *geant3 = (TGeant3*)gMC;
  18
  19
  20 //=======================================================================
  21 // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
  22  geant3->fGctrak->maxnst=1000000;
  23
  24 geant3->SetTRIG(1);          //Number of events to be processed
  25 geant3->SetSWIT(4,100);
  26 geant3->SetDEBU(0,0,1);
  27 geant3->SetDCAY(1);
  28 geant3->SetPAIR(1);
  29 geant3->SetCOMP(1);
  30 geant3->SetPHOT(1);
  31 geant3->SetPFIS(0);
  32 geant3->SetDRAY(1);
  33 geant3->SetANNI(1);
  34 geant3->SetBREM(1);
  35 geant3->SetMUNU(1);
  36 geant3->SetCKOV(0);
  37 geant3->SetHADR(4); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
  38 geant3->SetLOSS(1);
  39 geant3->SetMULS(1);
  40 geant3->SetRAYL(0);
  41 geant3->SetAUTO(0); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
  42 geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
  43 geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
  44 Float_t cut    = 1.e-1; // 100MeV cut by default
  45 Float_t tofmax = 1.e10;
  46 //              GAM    ELEC   NHAD   CHAD   MUON  EBREM  MUHAB EDEL MUDEL MUPA TOFMAX
  47 geant3->SetCUTS(1.e-4, 1.e-4, 1.e-3, 1.e-4, 1.e-4, cut,  cut,  cut, cut,  cut, 1.e-5);
  48
  49  gAlice->TrackingLimits( 700, 2200);
  50
  51 //
  52 //=======================================================================
  53 // ************* STEERING parameters FOR ALICE SIMULATION **************
  54 // --- Specify event type to be tracked through the ALICE setup
  55 // --- All positions are in cm, angles in degrees, and P and E in GeV
  56
  57  switch(gentype)
  58  {
  59  case gun:
  60 //*********************************************
  61 // Example for Fixed Particle Gun             *
  62 //*********************************************
  63      AliGenFixed *gener = new AliGenFixed(ntracks);
  64      gener->SetMomentum(3);
  65      gener->SetPhiRange(0);
  66      gener->SetThetaRange(0);
  67      gener->SetOrigin(0,0,0);          //vertex position
  68      gener->SetPart(kProton);          //GEANT particle type
  69      break;
  70  case box:
  71 //*********************************************
  72 // Example for Moving Particle Gun            *
  73 //*********************************************
  74      AliGenBox *gener = new AliGenBox(ntracks);
  75      gener->SetMomentumRange(3,4);
  76      gener->SetPhiRange(-360,360);
  77      gener->SetThetaRange(2., 10.);
  78      gener->SetOrigin(25,25,510.5);
  79      gener->SetVertexSmear(kPerTrack);
  80      //vertex position
  81      gener->SetSigma(1.8, 1.8,0);           //Sigma in (X,Y,Z) (cm) on IP position
  82      gener->SetPart(kProton);                    //GEANT particle type
  83      break;
  84  case scan:
  85 //*********************************************
  86 // Scanning on a grid                         *
  87 //*********************************************
  88      AliGenScan *gener = new AliGenScan(-1);
  89      gener->SetMomentumRange(20,20);
  90      gener->SetPhiRange(90,90);
  91      gener->SetThetaRange(0,0);
  92      //vertex position
  93      gener->SetSigma(3,3,0);           //Sigma in (X,Y,Z) (cm) on IP position
  94      gener->SetPart(kMuonMinus);
  95      gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
  96      break;
  97  case doublescan:
  98 //*********************************************
  99 // Scanning on a grid                         *
 100 //*********************************************
 101      AliGenDoubleScan *gener = new AliGenDoubleScan(-1);
 102      gener->SetMomentumRange(4,4);
 103      gener->SetPhiRange(0,360);
 104      gener->SetThetaRange(0,0);
 105      //vertex position
 106      gener->SetSigma(3,3,0);           //Sigma in (X,Y,Z) (cm) on IP position
 107      gener->SetPart(8);
 108      gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
 109      gener->SetDistance(1);
 110
 111      break;
 112
 113  case hijing:
 114      AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
 115      gener->SetMomentumRange(0,999);
 116      gener->SetPhiRange(0,360);
 117      gener->SetThetaRange(0.104,33.52);
 118      gener->SetOrigin(0,0,0);          //vertex position
 119      gener->SetSigma(0,0,5.6);         //Sigma in (X,Y,Z) (cm) on IP position
 120      break;
 121
 122  case pythia:
 123 //********************************************
 124 // Example for Charm  Production with Pythia *
 125 //********************************************
 126
 127      AliGenPythia *gener = new AliGenPythia(ntracks);
 128      gener->SetMomentumRange(0,999);
 129      gener->SetPhiRange(0,360);
 130      gener->SetThetaRange(0., 180.);
 131      gener->SetYRange(-10,10);
 132      gener->SetPtRange(0,100);
 133      gener->SetOrigin(0,0,0);          // vertex position
 134      gener->SetVertexSmear(kPerEvent);
 135      gener->SetSigma(0,0,5.6);         // Sigma in (X,Y,Z) (cm) on IP position
 136      gener->SetProcess(kPyCharm);
 137      gener->SetEnergyCMS(5500.);
 138      break;
 139
 140  case param:
 141 //*******************************************************
 142 // Example for J/psi  Production from  Parameterisation *
 143 //*******************************************************
 144      AliGenParam *gener = new AliGenParam(ntracks, AliGenMUONlib::kUpsilon);
 145      gener->SetMomentumRange(0,999);
 146      gener->SetPtRange(0,100.);
 147      gener->SetPhiRange(0,360);
 148      gener->SetYRange(2.5,4);
 149      gener->SetPtRange(0,50);
 150      gener->SetOrigin(0,0,0);          //vertex position
 151      gener->SetSigma(0,0,0);           //Sigma in (X,Y,Z) (cm) on IP position
 152      gener->SetForceDecay(kDiMuon);
 153      gener->SetCutOnChild(1);
 154      gener->SetChildThetaRange(2,9);
 155      gener->SetTrackingFlag(1);
 156
 157      break;
 158
 159  case fluka:
 160 //*******************************************************
 161 // Example for a FLUKA Boundary Source                  *
 162 //*******************************************************
 163      AliGenFLUKAsource *gener = new AliGenFLUKAsource(-1);
 164      gener->AddFile("$(ALICE_ROOT)/data/all32.root");
 165      rootfile->cd();
 166      gener->SetPartFlag(7);
 167      gener->SetMomentumRange(0,999);
 168      gener->SetPhiRange(0,360);
 169      gener->SetThetaRange(0., 180.);
 170      gener->SetAgeMax(1.e-5);
 171
 172 //  31.7 events
 173      gener->SetFraction(0.0315);
 174      break;
 175
 176  case ntuple:
 177 //*******************************************************
 178 // Example for reading from a external file                  *
 179 //*******************************************************
 180      AliGenExtFile *gener = new AliGenExtFile(-1);
 181      gener->SetFileName("$(ALICE_ROOT)/data/pbpb.root");
 182      gener->SetThetaRange(0.104,33.52);
 183      gener->SetOrigin(0,0,0);          //vertex position
 184      gener->SetSigma(0,0,5.6);         //Sigma in (X,Y,Z) (cm) on IP position
 185      gener->SetVertexSmear(kPerTrack);
 186      gener->SetTrackingFlag(1);
 187      break;
 188
 189  case halo:
 190 //*******************************************************
 191 // Example for Tunnel Halo Source                       *
 192 //*******************************************************
 193      AliGenHalo *gener = new AliGenHalo(ntracks);
 194      gener->SetFileName("/h1/morsch/marsip/marsip5.mu");
 195      break;
 196
 197  case cocktail:
 198 //*******************************************************
 199 // Example for a Cocktail                               *
 200 //*******************************************************
 201      AliGenCocktail *gener = new AliGenCocktail();
 202      gener->SetMomentumRange(0,10);
 203      gener->SetPhiRange(0,360);
 204      gener->SetThetaRange(45.,135);
 205
 206      pions   = new AliGenParam(100, AliGenMUONlib::kPion);
 207      gener->AddGenerator(pions  , "Pions"  , 100);
 208      break;
 209  }
 210
 211 // Activate this line if you want the vertex smearing to happen
 212 // track by track
 213 //
 214
 215 gener->Init();
 216 gAlice->SetField(-2,1);    //Specify maximum magnetic field in Tesla (neg. ==> default field)
 217 Int_t iMAG    =0;
 218 Int_t iABSO   =0;
 219 Int_t iDIPO   =0;
 220 Int_t iSHIL   =0;
 221 Int_t iPIPE   =0;
 222 Int_t iMUON   =1;
 223
 224 //=================== Alice BODY parameters =============================
 225 AliBODY *BODY = new AliBODY("BODY","Alice envelop");
 226
 227 if(iMAG) {
 228 //=================== MAG parameters ============================
 229 // --- Start with Magnet since detector layouts may be depending ---
 230 // --- on the selected Magnet dimensions ---
 231 AliMAG *MAG  = new AliMAG("MAG","Magnet");
 232 }
 233
 234 if(iABSO) {
 235 //=================== ABSO parameters ============================
 236     AliABSO *ABSO  = new AliABSOv1("ABSO","Muon Absorber");
 237 }
 238
 239 if(iDIPO) {
 240 //=================== DIPO parameters ============================
 241
 242     AliDIPO *DIPO  = new AliDIPOv2("DIPO","Dipole version 2");
 243 }
 244
 245
 246 if(iSHIL) {
 247 //=================== SHIL parameters ============================
 248     AliSHIL *SHIL  = new AliSHILv1("SHIL","Shielding");
 249 }
 250
 251
 252 if(iPIPE) {
 253 //=================== PIPE parameters ============================
 254     AliPIPE *PIPE  = new AliPIPEv0("PIPE","Beam Pipe");
 255 }
 256
 257
 258 if(iMUON) {
 259 //=================== MUON parameters ===========================
 260
 261 AliMUON *MUON  = new AliMUONv0("MUON","normal MUON");
 262  MUON->SetAcceptance(kTRUE, 0.1, 12.);
 263  MUON->SetIshunt(0);
 264  MUON->SetMaxStepGas(0.1);
 265  MUON->SetMaxStepAlu(0.1);
 266 }
 267 }