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