EMCAL/Config.C

   1 enum PprRun_t
   2 {
   3     kTest50,
   4     kParam_8000,   kParam_4000,  kParam_2000,
   5     kHijing_cent1, kHijing_cent2,
   6     kHijing_per1,  kHijing_per2, kHijing_per3, kHijing_per4,  kHijing_per5,
   7     kHijing_jj25,  kHijing_jj50, kHijing_jj75, kHijing_jj100, kHijing_jj125,
   8     kHijing_gj25,  kHijing_gj50, kHijing_gj75, kHijing_gj100, kHijing_gj125,
   9     kJetPlusBg,    kGammaPlusBg,
  10     kParam_8000_Ecal, kParam_4000_Ecal,
  11     kJets_100,
  12     kGammaGun, kGammaBox
  13 };
  14
  15 enum PprGeo_t
  16 {
  17     kHoles, kNoHoles
  18 };
  19
  20 enum PprRad_t
  21 {
  22     kGluonRadiation, kNoGluonRadiation
  23 };
  24
  25 // This part for configuration
  26 //static PprRun_t run = kParam_8000_Ecal;
  27 static PprRun_t run = kHijing_jj100;
  28 static PprGeo_t geo = kHoles;
  29 static PprRad_t rad = kGluonRadiation;
  30 // Comment line
  31 static TString  comment;
  32
  33 void Config()
  34 {
  35
  36     // 7-DEC-2000 09:00
  37     // Switch on Transition adiation simulation. 6/12/00 18:00
  38     // iZDC=1  7/12/00 09:00
  39     // ThetaRange is (0., 180.). It was (0.28,179.72) 7/12/00 09:00
  40     // Theta range given through pseudorapidity limits 22/6/2001
  41
  42     // Set Random Number seed
  43     // gRandom->SetSeed(12345);
  44
  45     new AliGeant3("C++ Interface to Geant3");
  46
  47     if (!gSystem->Getenv("CONFIG_FILE"))
  48     {
  49         TFile  *rootfile = new TFile("galice.root", "recreate");
  50
  51         rootfile->SetCompressionLevel(2);
  52     }
  53
  54     TGeant3 *geant3 = (TGeant3 *) gMC;
  55
  56     //
  57     // Set External decayer
  58     AliDecayer *decayer = new AliDecayerPythia();
  59
  60     decayer->SetForceDecay(kAll);
  61     decayer->Init();
  62     gMC->SetExternalDecayer(decayer);
  63     //
  64     //
  65     //=======================================================================
  66     // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
  67     geant3->SetTRIG(1);         //Number of events to be processed
  68     geant3->SetSWIT(4, 10);
  69     geant3->SetDEBU(0, 0, 1);
  70     geant3->SetDCAY(1);
  71     geant3->SetPAIR(1);
  72     geant3->SetCOMP(1);
  73     geant3->SetPHOT(1);
  74     geant3->SetPFIS(0);
  75     geant3->SetDRAY(0);
  76     geant3->SetANNI(1);
  77     geant3->SetBREM(1);
  78     geant3->SetMUNU(1);
  79     geant3->SetCKOV(1);
  80     geant3->SetHADR(1);         //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
  81     geant3->SetLOSS(2);
  82     geant3->SetMULS(1);
  83     geant3->SetRAYL(1);
  84     geant3->SetAUTO(1);         //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
  85     geant3->SetABAN(0);         //Restore 3.16 behaviour for abandoned tracks
  86     geant3->SetOPTI(2);         //Select optimisation level for GEANT geometry searches (0,1,2)
  87     geant3->SetERAN(5.e-7);
  88
  89     Float_t cut    = 1.e-3;        // 1MeV cut by default
  90     Float_t tofmax = 1.e10;
  91
  92     //             GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
  93     geant3->SetCUTS(cut, cut, cut, cut, cut, cut, cut, cut, cut, cut, tofmax);
  94     gAlice->TrackingLimits(900, 900);
  95
  96     //
  97     //=======================================================================
  98     // ************* STEERING parameters FOR ALICE SIMULATION **************
  99     // --- Specify event type to be tracked through the ALICE setup
 100     // --- All positions are in cm, angles in degrees, and P and E in GeV
 101
 102
 103 // Generator Configuration
 104     gAlice->SetDebug(1);
 105
 106     AliGenerator*  gener = GeneratorFactory(run);
 107     gener->Init();
 108     gener->SetTrackingFlag(1);
 109
 110
 111
 112     if (rad == kGluonRadiation)
 113     {
 114         //coment= comment.Append(" | Gluon Radiation On");
 115
 116     } else {
 117         //coment= comment.Append(" | Gluon Radiation Off");
 118     }
 119
 120     if (geo == kHoles)
 121     {
 122         //coment= comment.Append(" | Holes for PHOS/RICH");
 123
 124     } else {
 125         //coment= comment.Append(" | No holes for PHOS/RICH");
 126     }
 127
 128     printf("\n \n Comment: %s \n \n", (char*) comment);
 129
 130
 131 // Field (L3 0.4 T)
 132
 133     AliMagFCM* field = new AliMagFCM(
 134         "Map2","$(ALICE_ROOT)/data/field01.dat", 2, 1., 10.);
 135     field->SetSolenoidField(4.);
 136     gAlice->SetField(field);
 137
 138     Int_t iMAG    = 0;
 139     Int_t iITS    = 0;
 140     Int_t iTPC    = 0;
 141     Int_t iTOF    = 0;
 142     Int_t iRICH   = 0;
 143     Int_t iZDC    = 0;
 144     Int_t iCASTOR = 0;
 145     Int_t iTRD    = 0;
 146     Int_t iABSO   = 0;
 147     Int_t iDIPO   = 0;
 148     Int_t iHALL   = 0;
 149     Int_t iFRAME  = 0;
 150     Int_t iSHIL   = 0;
 151     Int_t iPIPE   = 0;
 152     Int_t iFMD    = 0;
 153     Int_t iMUON   = 0;
 154     Int_t iPHOS   = 0;
 155     Int_t iPMD    = 0;
 156     Int_t iSTART  = 0;
 157     Int_t iVZERO  = 0;
 158     Int_t iEMCAL  = 1;
 159 //
 160 //  enable/disable StepManager()
 161 //
 162
 163     Int_t   enableABSO   = 0;
 164     Int_t   enableCASTOR = 0;
 165     Int_t   enableDIPO   = 0;
 166     Int_t   enableFMD    = 0;
 167     Int_t   enableFRAME  = 0;
 168     Int_t   enableHALL   = 0;
 169     Int_t   enableITS    = 0;
 170     Int_t   enableMAG    = 0;
 171     Int_t   enableMUON   = 0;
 172     Int_t   enablePHOS   = 1;
 173     Int_t   enablePIPE   = 0;
 174     Int_t   enablePMD    = 0;
 175     Int_t   enableRICH   = 0;
 176     Int_t   enableSHIL   = 0;
 177     Int_t   enableSTART  = 0;
 178     Int_t   enableTOF    = 0;
 179     Int_t   enableTPC    = 0;
 180     Int_t   enableTRD    = 0;
 181     Int_t   enableZDC    = 0;
 182     Int_t   enableEMCAL  = 1;
 183
 184     //=================== Alice BODY parameters =============================
 185     AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
 186
 187
 188     if (iMAG)
 189     {
 190         //=================== MAG parameters ============================
 191         // --- Start with Magnet since detector layouts may be depending ---
 192         // --- on the selected Magnet dimensions ---
 193         AliMAG *MAG = new AliMAG("MAG", "Magnet");
 194     }
 195
 196
 197     if (iABSO)
 198     {
 199         //=================== ABSO parameters ============================
 200         AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
 201     }
 202
 203     if (iDIPO)
 204     {
 205         //=================== DIPO parameters ============================
 206
 207         AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
 208     }
 209
 210     if (iHALL)
 211     {
 212         //=================== HALL parameters ============================
 213
 214         AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
 215     }
 216
 217
 218     if (iFRAME)
 219     {
 220         //=================== FRAME parameters ============================
 221
 222         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
 223         if (geo == kHoles) {
 224             FRAME->SetHoles(1);
 225         } else {
 226             FRAME->SetHoles(0);
 227         }
 228     }
 229
 230     if (iSHIL)
 231     {
 232         //=================== SHIL parameters ============================
 233
 234         AliSHIL *SHIL = new AliSHILv0("SHIL", "Shielding");
 235     }
 236
 237
 238     if (iPIPE)
 239     {
 240         //=================== PIPE parameters ============================
 241
 242         AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
 243     }
 244
 245   if(iITS) {
 246
 247 //=================== ITS parameters ============================
 248     //
 249     // As the innermost detector in ALICE, the Inner Tracking System "impacts" on
 250     // almost all other detectors. This involves the fact that the ITS geometry
 251     // still has several options to be followed in parallel in order to determine
 252     // the best set-up which minimizes the induced background. All the geometries
 253     // available to date are described in the following. Read carefully the comments
 254     // and use the default version (the only one uncommented) unless you are making
 255     // comparisons and you know what you are doing. In this case just uncomment the
 256     // ITS geometry you want to use and run Aliroot.
 257     //
 258     // Detailed geometries:
 259     //
 260     //
 261     //AliITS *ITS  = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
 262     //
 263     //AliITS *ITS  = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
 264     //
 265     AliITSvPPRasymm *ITS  = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
 266     ITS->SetMinorVersion(2);                                     // don't touch this parameter if you're not an ITS developer
 267     ITS->SetReadDet(kFALSE);                                     // don't touch this parameter if you're not an ITS developer
 268     //    ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");  // don't touch this parameter if you're not an ITS developer
 269     ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
 270     ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
 271     ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
 272     ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
 273     ITS->SetRails(1);        // 1 --> rails in ; 0 --> rails out
 274     ITS->SetCoolingFluid(1);   // 1 --> water ; 0 --> freon
 275     //
 276     //AliITSvPPRsymm *ITS  = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
 277     //ITS->SetMinorVersion(2);                                       // don't touch this parameter if you're not an ITS developer
 278     //ITS->SetReadDet(kFALSE);                                       // don't touch this parameter if you're not an ITS developer
 279     //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det"); // don't touch this parameter if you're not an ITS developer
 280     //ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
 281     //ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
 282     //ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
 283     //ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
 284     //ITS->SetRails(1);              // 1 --> rails in ; 0 --> rails out
 285     //ITS->SetCoolingFluid(1);       // 1 --> water ; 0 --> freon
 286     //
 287     //
 288     // Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful
 289     // for reconstruction !):
 290     //
 291     //
 292     //AliITSvPPRcoarseasymm *ITS  = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
 293     //ITS->SetRails(1);                // 1 --> rails in ; 0 --> rails out
 294     //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
 295     //
 296     //AliITS *ITS  = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
 297     //ITS->SetRails(1);                // 1 --> rails in ; 0 --> rails out
 298     //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
 299     //
 300     //
 301     //
 302     // Geant3 <-> EUCLID conversion
 303     // ============================
 304     //
 305     // SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
 306     // media to two ASCII files (called by default ITSgeometry.euc and
 307     // ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
 308     // The default (=0) means that you dont want to use this facility.
 309     //
 310     ITS->SetEUCLID(0);
 311     if (!enableITS) ITS->DisableStepManager();
 312   }
 313
 314
 315     if (iTPC)
 316     {
 317         //============================ TPC parameters ================================
 318         // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
 319         // --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
 320         // --- sectors are specified, any value other than that requires at least one
 321         // --- sector (lower or upper)to be specified!
 322         // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
 323         // ---           sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
 324         // --- SecLows - number of lower sectors specified (up to 6)
 325         // --- SecUps - number of upper sectors specified (up to 12)
 326         // --- Sens - sensitive strips for the Slow Simulator !!!
 327         // --- This does NOT work if all S or L-sectors are specified, i.e.
 328         // --- if SecAL or SecAU < 0
 329         //
 330         //
 331         //-----------------------------------------------------------------------------
 332
 333         //  gROOT->LoadMacro("SetTPCParam.C");
 334         //  AliTPCParam *param = SetTPCParam();
 335         AliTPC *TPC = new AliTPCv2("TPC", "Default");
 336
 337         // All sectors included
 338         TPC->SetSecAL(-1);
 339         TPC->SetSecAU(-1);
 340         if (!enableTPC) TPC->DisableStepManager();
 341     }
 342
 343
 344     if (iRICH)
 345     {
 346         //=================== RICH parameters ===========================
 347         AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
 348         if (!enableRICH) RICH->DisableStepManager();
 349     }
 350
 351
 352     if (iZDC)
 353     {
 354         //=================== ZDC parameters ============================
 355
 356         AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
 357         if (!enableZDC) ZDC->DisableStepManager();
 358     }
 359
 360     if (iCASTOR)
 361     {
 362         //=================== CASTOR parameters ============================
 363
 364         AliCASTOR *CASTOR = new AliCASTORv1("CASTOR", "normal CASTOR");
 365         if (!enableCASTOR) CASTOR->DisableStepManager();
 366     }
 367
 368     if (iTRD)
 369     {
 370         //=================== TRD parameters ============================
 371
 372         AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
 373
 374         // Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
 375         TRD->SetGasMix(1);
 376         if (geo == kHoles) {
 377             // With hole in front of PHOS
 378             TRD->SetPHOShole();
 379             // With hole in front of RICH
 380             TRD->SetRICHhole();
 381         }
 382             // Switch on TR
 383             AliTRDsim *TRDsim = TRD->CreateTR();
 384             if (!enableTRD) TRD->DisableStepManager();
 385     }
 386
 387     if (iFMD)
 388     {
 389         //=================== FMD parameters ============================
 390
 391         AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
 392         FMD->SetRingsSi1(256);
 393         FMD->SetRingsSi2(64);
 394         FMD->SetSectorsSi1(20);
 395         FMD->SetSectorsSi2(24);
 396         if (!enableFMD) FMD->DisableStepManager();
 397    }
 398
 399     if (iMUON)
 400     {
 401         //=================== MUON parameters ===========================
 402
 403         AliMUON *MUON = new AliMUONv1("MUON", "default");
 404         if (!enableMUON) MUON->DisableStepManager();
 405     }
 406     //=================== PHOS parameters ===========================
 407
 408     if (iPHOS)
 409     {
 410         AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
 411         if (!enablePHOS) PHOS->DisableStepManager();
 412     }
 413
 414
 415     if (iPMD)
 416     {
 417         //=================== PMD parameters ============================
 418         AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
 419         PMD->SetPAR(1., 1., 0.8, 0.02);
 420         PMD->SetIN(6., 18., -580., 27., 27.);
 421         PMD->SetGEO(0.0, 0.2, 4.);
 422         PMD->SetPadSize(0.8, 1.0, 1.0, 1.5);
 423         if (!enablePMD) PMD->DisableStepManager();
 424     }
 425
 426     if (iTOF) {
 427         if (geo == kHoles) {
 428         //=================== TOF parameters ============================
 429             AliTOF *TOF = new AliTOFv2FHoles("TOF", "TOF with Holes");
 430         } else {
 431             AliTOF *TOF = new AliTOFv4("TOF", "normal TOF");
 432         }
 433         if (!enableTOF) TOF->DisableStepManager();
 434     }
 435
 436     if (iSTART)
 437     {
 438         //=================== START parameters ============================
 439         AliSTART *START = new AliSTARTv1("START", "START Detector");
 440         if (!enableSTART) START->DisableStepManager();
 441     }
 442
 443     if (iEMCAL && !iRICH)
 444     {
 445         //=================== EMCAL parameters ============================
 446         AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCALArch1b");
 447         if (!enableEMCAL) EMCAL->DisableStepManager();
 448     }
 449 }
 450
 451 Float_t EtaToTheta(Float_t arg){
 452   return (180./TMath::Pi())*2.*atan(exp(-arg));
 453 }
 454
 455
 456
 457 AliGenerator* GeneratorFactory(PprRun_t run) {
 458     Int_t isw = 3;
 459     if (rad == kNoGluonRadiation) isw = 0;
 460
 461
 462     switch (run) {
 463     case kGammaGun:
 464         gener = new AliGenFixed(1);
 465         gener->SetMomentum(100.);
 466         gener->SetPhi(240.);
 467         gener->SetTheta(91.);
 468         gener->SetPart(kElectron);
 469      break;
 470
 471     case kGammaBox:
 472         gener = new AliGenBox(100);
 473         gener->SetMomentumRange(100., 100.1);
 474         gener->SetPhiRange(0,360);
 475         gener->SetThetaRange(45., 135.);
 476         gener->SetPart(kElectron);
 477         break;
 478
 479     case kTest50:
 480         //coment= comment.Append(":HIJINGparam test 50 particles");
 481         gener = new AliGenHIJINGparaBa(50);
 482         gener->SetMomentumRange(0, 999999.);
 483         gener->SetPhiRange(-180., 180.);
 484         // Set pseudorapidity range from -8 to 8.
 485         Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
 486         Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min
 487         gener->SetThetaRange(thmin,thmax);
 488         break;
 489
 490     case kParam_8000:
 491         //coment= comment.Append(":HIJINGparam N=8000");
 492         AliGenHIJINGpara *gener = new AliGenHIJINGpara(86030);
 493         gener->SetMomentumRange(0, 999999.);
 494         gener->SetPhiRange(-180., 180.);
 495         // Set pseudorapidity range from -8 to 8.
 496         Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
 497         Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min
 498         gener->SetThetaRange(thmin,thmax);
 499         break;
 500     case kParam_4000:
 501         //coment= comment.Append("HIJINGparam N=4000");
 502         AliGenHIJINGpara *gener = new AliGenHIJINGpara(43015);
 503         gener->SetMomentumRange(0, 999999.);
 504         gener->SetPhiRange(-180., 180.);
 505         // Set pseudorapidity range from -8 to 8.
 506         Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
 507         Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min
 508         gener->SetThetaRange(thmin,thmax);
 509         break;
 510     case kParam_2000:
 511         //coment= comment.Append("HIJINGparam N=2000");
 512         AliGenHIJINGpara *gener = new AliGenHIJINGpara(21507);
 513         gener->SetMomentumRange(0, 999999.);
 514         gener->SetPhiRange(-180., 180.);
 515         // Set pseudorapidity range from -8 to 8.
 516         Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
 517         Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min
 518         gener->SetThetaRange(thmin,thmax);
 519         break;
 520
 521     case kParam_8000_Ecal:
 522         //coment= comment.Append(":HIJINGparam N=8000 baryons");
 523         AliGenHIJINGparaBa *gener = new AliGenHIJINGparaBa(82534);
 524         gener->SetMomentumRange(0, 999999.);
 525         gener->SetPhiRange(-180., 180.);
 526         // Set pseudorapidity range from -8 to 8.
 527         Float_t thmin = EtaToTheta( 5);   // theta min. <---> eta max
 528         Float_t thmax = EtaToTheta(-5);  // theta max. <---> eta min
 529         gener->SetThetaRange(thmin,thmax);
 530         break;
 531
 532     case kParam_4000_Ecal:
 533         //coment= comment.Append(":HIJINGparam N=8000 baryons");
 534         AliGenHIJINGparaBa *gener = new AliGenHIJINGparaBa(82534/2.);
 535         gener->SetMomentumRange(0, 999999.);
 536         gener->SetPhiRange(-180., 180.);
 537         // Set pseudorapidity range from -8 to 8.
 538         Float_t thmin = EtaToTheta( 5);   // theta min. <---> eta max
 539         Float_t thmax = EtaToTheta(-5);  // theta max. <---> eta min
 540         gener->SetThetaRange(thmin,thmax);
 541         break;
 542 //
 543 //  Hijing Central
 544 //
 545     case kHijing_cent1:
 546         //coment= comment.Append("HIJING cent1");
 547         AliGenHijing *gener = HijingStandard();
 548 // impact parameter range
 549         gener->SetImpactParameterRange(0., 5.);
 550         break;
 551     case kHijing_cent2:
 552         //coment= comment.Append("HIJING cent2");
 553         AliGenHijing *gener = HijingStandard();
 554 // impact parameter range
 555         gener->SetImpactParameterRange(0., 2.);
 556         break;
 557 //
 558 // Hijing Peripheral
 559 //
 560     case kHijing_per1:
 561         //coment= comment.Append("HIJING per1");
 562         AliGenHijing *gener = HijingStandard();
 563 // impact parameter range
 564         gener->SetImpactParameterRange(5., 8.6);
 565         break;
 566     case kHijing_per2:
 567         //coment= comment.Append("HIJING per2");
 568         AliGenHijing *gener = HijingStandard();
 569 // impact parameter range
 570         gener->SetImpactParameterRange(8.6, 11.2);
 571         break;
 572     case kHijing_per3:
 573         //coment= comment.Append("HIJING per3");
 574         AliGenHijing *gener = HijingStandard();
 575 // impact parameter range
 576         gener->SetImpactParameterRange(11.2, 13.2);
 577         break;
 578     case kHijing_per4:
 579         //coment= comment.Append("HIJING per4");
 580         AliGenHijing *gener = HijingStandard();
 581 // impact parameter range
 582         gener->SetImpactParameterRange(13.2, 15.);
 583         break;
 584     case kHijing_per5:
 585         //coment= comment.Append("HIJING per5");
 586         AliGenHijing *gener = HijingStandard();
 587 // impact parameter range
 588         gener->SetImpactParameterRange(15., 100.);
 589         break;
 590 //
 591 //  Jet-Jet
 592 //
 593     case kHijing_jj25:
 594         //coment= comment.Append("HIJING Jet 25 GeV");
 595         AliGenHijing *gener = HijingStandard();
 596 // impact parameter range
 597         gener->SetImpactParameterRange(0., 5.);
 598         // trigger
 599         gener->SetTrigger(1);
 600         gener->SetPtJet(25.);
 601         gener->SetSimpleJets(1);
 602         gener->SetRadiation(isw);
 603         gener->SetJetEtaRange(-0.3,0.3);
 604         gener->SetJetPhiRange(15.,105.);
 605         break;
 606
 607     case kHijing_jj50:
 608         //coment= comment.Append("HIJING Jet 50 GeV");
 609         AliGenHijing *gener = HijingStandard();
 610 // impact parameter range
 611         gener->SetImpactParameterRange(0., 5.);
 612         // trigger
 613         gener->SetTrigger(1);
 614         gener->SetPtJet(50.);
 615         gener->SetSimpleJets(1);
 616         gener->SetRadiation(isw);
 617         gener->SetJetEtaRange(-0.3,0.3);
 618         gener->SetJetPhiRange(15.,105.);
 619         break;
 620
 621     case kHijing_jj75:
 622         //coment= comment.Append("HIJING Jet 75 GeV");
 623         AliGenHijing *gener = HijingStandard();
 624 // impact parameter range
 625         gener->SetImpactParameterRange(0., 5.);
 626         // trigger
 627         gener->SetTrigger(1);
 628         gener->SetPtJet(75.);
 629         gener->SetSimpleJets(1);
 630         gener->SetRadiation(isw);
 631         gener->SetJetEtaRange(-0.3,0.3);
 632         gener->SetJetPhiRange(15.,105.);
 633         break;
 634
 635     case kHijing_jj100:
 636         //coment= comment.Append("HIJING Jet 100 GeV");
 637         AliGenHijing *gener = HijingStandard();
 638 // impact parameter range
 639         gener->SetImpactParameterRange(10., 10.1);
 640         // trigger
 641         gener->SetTrigger(1);
 642         gener->SetPtJet(100.);
 643         gener->SetSimpleJets(1);
 644         gener->SetRadiation(isw);
 645         gener->SetJetEtaRange(-0.3,0.3);
 646         gener->SetJetPhiRange(15.,105.);
 647         break;
 648
 649     case kHijing_jj125:
 650         //coment= comment.Append("HIJING Jet 125 GeV");
 651         AliGenHijing *gener = HijingStandard();
 652 // impact parameter range
 653         gener->SetImpactParameterRange(0., 5.);
 654         // trigger
 655         gener->SetTrigger(1);
 656         gener->SetPtJet(125.);
 657         gener->SetSimpleJets(1);
 658         gener->SetRadiation(isw);
 659         gener->SetJetEtaRange(-0.3,0.3);
 660         gener->SetJetPhiRange(15.,105.);
 661         break;
 662 //
 663 // Gamma-Jet
 664 //
 665     case kHijing_gj25:
 666         //coment= comment.Append("HIJING Gamma 25 GeV");
 667         AliGenHijing *gener = HijingStandard();
 668 // impact parameter range
 669         gener->SetImpactParameterRange(0., 5.);
 670         // trigger
 671         gener->SetTrigger(2);
 672         gener->SetPtJet(25.);
 673         gener->SetSimpleJets(1);
 674         gener->SetRadiation(isw);
 675         gener->SetJetEtaRange(-0.3,0.3);
 676         gener->SetJetPhiRange(15.,105.);
 677         break;
 678
 679     case kHijing_gj50:
 680         //coment= comment.Append("HIJING Gamma 50 GeV");
 681         AliGenHijing *gener = HijingStandard();
 682 // impact parameter range
 683         gener->SetImpactParameterRange(0., 5.);
 684         // trigger
 685         gener->SetTrigger(2);
 686         gener->SetPtJet(50.);
 687         gener->SetSimpleJets(1);
 688         gener->SetRadiation(isw);
 689         gener->SetJetEtaRange(-0.3,0.3);
 690         gener->SetJetPhiRange(15.,105.);
 691         break;
 692
 693     case kHijing_gj75:
 694         //coment= comment.Append("HIJING Gamma 75 GeV");
 695         AliGenHijing *gener = HijingStandard();
 696 // impact parameter range
 697         gener->SetImpactParameterRange(0., 5.);
 698         // trigger
 699         gener->SetTrigger(2);
 700         gener->SetPtJet(75.);
 701         gener->SetSimpleJets(1);
 702         gener->SetRadiation(isw);
 703         gener->SetJetEtaRange(-0.3,0.3);
 704         gener->SetJetPhiRange(15.,105.);
 705         break;
 706
 707     case kHijing_gj100:
 708         //coment= comment.Append("HIJING Gamma 100 GeV");
 709         AliGenHijing *gener = HijingStandard();
 710 // impact parameter range
 711         gener->SetImpactParameterRange(0., 5.);
 712         // trigger
 713         gener->SetTrigger(2);
 714         gener->SetPtJet(100.);
 715         gener->SetSimpleJets(1);
 716         gener->SetRadiation(isw);
 717         gener->SetJetEtaRange(-0.3,0.3);
 718         gener->SetJetPhiRange(15.,105.);
 719         break;
 720
 721     case kHijing_gj125:
 722         //coment= comment.Append("HIJING Gamma 125 GeV");
 723         AliGenHijing *gener = HijingStandard();
 724 // impact parameter range
 725         gener->SetImpactParameterRange(0., 5.);
 726         // trigger
 727         gener->SetTrigger(2);
 728         gener->SetPtJet(125.);
 729         gener->SetSimpleJets(1);
 730         gener->SetRadiation(isw);
 731         gener->SetJetEtaRange(-0.3,0.3);
 732         gener->SetJetPhiRange(15.,105.);
 733         break;
 734     case kJetPlusBg:
 735         AliGenCocktail *gener = new AliGenCocktail();
 736         gener->SetMomentumRange(0, 999999.);
 737         gener->SetPhiRange(-180., 180.);
 738         // Set pseudorapidity range from -8 to 8.
 739         Float_t thmin = EtaToTheta( 5.);   // theta min. <---> eta max
 740         Float_t thmax = EtaToTheta(-5.);  // theta max. <---> eta min
 741         gener->SetThetaRange(thmin,thmax);
 742
 743 //
 744 //      Underlying Event
 745 //
 746 //      AliGenHIJINGparaBa *bg = new AliGenHIJINGparaBa(82534);
 747         AliGenHIJINGparaBa *bg = new AliGenHIJINGparaBa(10);
 748 //
 749 //      Jets from Pythia
 750 //
 751         AliGenPythia *jets = new AliGenPythia(-1);
 752 //   Centre of mass energy
 753         jets->SetEnergyCMS(5500.);
 754 //   Process type
 755         jets->SetProcess(kPyJets);
 756 //   final state kinematic cuts
 757         jets->SetJetEtaRange(-0.3, 0.3);
 758         jets->SetJetPhiRange(15., 105.);
 759 //   Structure function
 760         jets->SetStrucFunc(kGRV_LO_98);
 761 //
 762 //   Pt transfer of the hard scattering
 763         jets->SetPtHard(100.,100.1);
 764 //   Decay type (semielectronic, semimuonic, nodecay)
 765         jets->SetForceDecay(kAll);
 766 //
 767 //      Add all to cockail ...
 768 //
 769         gener->AddGenerator(jets,"Jets",1);
 770         gener->AddGenerator(bg,"Underlying Event", 1);
 771
 772      break;
 773     case kGammaPlusBg:
 774         AliGenCocktail *gener = new AliGenCocktail();
 775         gener->SetMomentumRange(0, 999999.);
 776         gener->SetPhiRange(-180., 180.);
 777         // Set pseudorapidity range from -8 to 8.
 778         Float_t thmin = EtaToTheta( 5.);   // theta min. <---> eta max
 779         Float_t thmax = EtaToTheta(-5.);  // theta max. <---> eta min
 780         gener->SetThetaRange(thmin,thmax);
 781
 782 //
 783 //      Underlying Event
 784 //
 785         AliGenHIJINGparaBa *bg = new AliGenHIJINGparaBa(82534);
 786 //
 787 //      Jets from Pythia
 788 //
 789         AliGenPythia *jets = new AliGenPythia(-1);
 790 //   Centre of mass energy
 791         jets->SetEnergyCMS(5500.);
 792 //   Process type
 793         jets->SetProcess(kPyDirectGamma);
 794 //   final state kinematic cuts
 795         jets->SetJetEtaRange(-0.3, 0.3);
 796         jets->SetJetPhiRange(15., 105.);
 797         jets->SetGammaEtaRange(-0.12, 0.12);
 798         jets->SetGammaPhiRange(220., 320.);
 799 //   Structure function
 800         jets->SetStrucFunc(kGRV_LO_98);
 801 //
 802 //   Pt transfer of the hard scattering
 803         jets->SetPtHard(100.,100.1);
 804 //   Decay type (semielectronic, semimuonic, nodecay)
 805         jets->SetForceDecay(kAll);
 806 //
 807 //      Add all to cockail ...
 808 //
 809         gener->AddGenerator(jets,"Jets",1);
 810         gener->AddGenerator(bg,"Underlying Event", 1);
 811
 812         break;
 813     case kJets_100:
 814 //  100 GeV Jets
 815         AliGenPythia *gener = PythiaJets(100.);
 816         break;
 817     }
 818
 819     return gener;
 820 }
 821
 822 AliGenHijing* HijingStandard()
 823 {
 824     AliGenHijing *gener = new AliGenHijing(-1);
 825 // centre of mass energy
 826     gener->SetEnergyCMS(5500.);
 827 // reference frame
 828     gener->SetReferenceFrame("CMS");
 829 // projectile
 830     gener->SetProjectile("A", 208, 82);
 831     gener->SetTarget    ("A", 208, 82);
 832 // tell hijing to keep the full parent child chain
 833     gener->KeepFullEvent();
 834 // enable jet quenching
 835     gener->SetJetQuenching(1);
 836 // enable shadowing
 837     gener->SetShadowing(1);
 838 // neutral pion and heavy particle decays switched off
 839     gener->SetDecaysOff(1);
 840 // Don't track spectators
 841     gener->SetSpectators(0);
 842 // kinematic selection
 843     gener->SetSelectAll(0);
 844     return gener;
 845 }
 846
 847 AliGenPythia* PythiaJets(Float_t energy)
 848 {
 849     AliGenPythia *gener = new AliGenPythia(-1);
 850 //   Centre of mass energy
 851     gener->SetEnergyCMS(5500.);
 852 //   Process type
 853     gener->SetProcess(kPyJets);
 854 //   final state kinematic cuts
 855     gener->SetJetEtaRange(-0.3, 0.3);
 856     gener->SetJetPhiRange(15., 105.);
 857 //   Structure function
 858     gener->SetStrucFunc(kGRV_LO_98);
 859 //
 860 //   Pt transfer of the hard scattering
 861     gener->SetPtHard(energy, energy+0.1);
 862 //   Decay type (semielectronic, semimuonic, nodecay)
 863     gener->SetForceDecay(kAll);
 864 //
 865     return gener;
 866 }
 867
 868
 869 AliGenPythia* PythiaGamma(Float_t energy)
 870 {
 871     AliGenPythia *gener = new AliGenPythia(-1);
 872 //   Centre of mass energy
 873     gener->SetEnergyCMS(5500.);
 874 //   Process type
 875     gener->SetProcess(kPyDirectGamma);
 876 //   final state kinematic cuts
 877     gener->SetJetEtaRange(-0.3, 0.3);
 878     gener->SetJetPhiRange(15., 105.);
 879     gener->SetGammaEtaRange(-0.12, 0.12);
 880     gener->SetGammaPhiRange(220., 320.);
 881 //   Structure function
 882     gener->SetStrucFunc(kGRV_LO_98);
 883 //
 884 //   Pt transfer of the hard scattering
 885     gener->SetPtHard(energy, energy+0.1);
 886 //   Decay type (semielectronic, semimuonic, nodecay)
 887     gener->SetForceDecay(kAll);
 888 //
 889     return gener;
 890 }
 891
 892
 893
 894
 895
 896
 897
 898
 899
 900
 901