From c6b905e5ec3a2518f38f7501feafd0bf39f318e0 Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Date: Thu, 25 Jul 2024 10:18:21 +0200 Subject: [PATCH 1/2] Remove protection agains Root6 for Running Starlight --- MC/GeneratorConfig.C | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/MC/GeneratorConfig.C b/MC/GeneratorConfig.C index 640ac8df..600ef6cd 100644 --- a/MC/GeneratorConfig.C +++ b/MC/GeneratorConfig.C @@ -1256,10 +1256,10 @@ GeneratorHijing() AliGenerator * GeneratorStarlight(){ -#if ROOT_VERSION_CODE >= ROOT_VERSION(6,0,0) - printf("ERROR: Starlight not yet compatible with ROOT 6!\n"); - return NULL; -#else +//#if ROOT_VERSION_CODE >= ROOT_VERSION(6,0,0) +// printf("ERROR: Starlight not yet compatible with ROOT 6!\n"); +// return NULL; +//#else gSystem->Load("libStarLight.so"); gSystem->Load("libAliStarLight.so"); @@ -1504,7 +1504,7 @@ GeneratorStarlight(){ genCocktail->AddGenerator(genStarLight,"StarLight",1.); genCocktail->AddGenerator(genEvtGen,"EvtGen",1.); return genCocktail; -#endif +//#endif } From 3b3819d604ffdc4005759ba9324b94597b1b3fb7 Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Date: Thu, 25 Jul 2024 11:18:35 +0200 Subject: [PATCH 2/2] combine Hijing and Pythia jets for p-Pb, adapted from MC/CustomGenerators/PWGHF/EPS_HFbjetpPb.C, needed for ALIROOT-8867 --- MC/CustomGenerators/PWGHF/HIJING_HFbjetpPb.C | 34 +++++++++++++------- 1 file changed, 23 insertions(+), 11 deletions(-) diff --git a/MC/CustomGenerators/PWGHF/HIJING_HFbjetpPb.C b/MC/CustomGenerators/PWGHF/HIJING_HFbjetpPb.C index 2c0eb935..684e0604 100644 --- a/MC/CustomGenerators/PWGHF/HIJING_HFbjetpPb.C +++ b/MC/CustomGenerators/PWGHF/HIJING_HFbjetpPb.C @@ -1,19 +1,24 @@ AliGenerator * GeneratorCustom(TString opt = "") { + + // set the xmldoc path using PYTHIA8DATA enviroement var + gSystem->Setenv("PYTHIA8DATA", gSystem->ExpandPathName("$ALICE_ROOT/PYTHIA8/pythia8/xmldoc")); + AliGenCocktail *ctl = (AliGenCocktail*) GeneratorCocktail("HIJING_HF"); Float_t randHF = gRandom->Rndm(); - if(randHF>0.176797){ // add HIJING generator for p-Pb - AliGenerator *hij = GeneratorHijing(); - ctl->AddGenerator(hij, "Hijing", 1.); + if(randHF>0.176797) + { // add HIJING generator for p-Pb + AliGenerator *hijing = GeneratorHijing(); + ctl->AddGenerator(hijing, "HIJING", 1.); } // Int_t process[2] = {kPythia6HeavyProcess_Charm, kPythia6HeavyProcess_Beauty}; - Int_t decay[4] = {kPythia6HeavyDecay_Hadrons, kPythia6HeavyDecay_HadronsWithV0, kPythia6HeavyDecay_Electron, kPythia6HeavyDecay_All}; + Int_t decay[5] = {kPythia6HeavyDecay_Hadrons, kPythia6HeavyDecay_HadronsWithV0, kPythia6HeavyDecay_Electron, kPythia6HeavyDecay_All, kPythia6HeavyDecay_All_bDecaysEvtGen}; - const Char_t *label[2][4] = { - "chadr PYTHIA", "chadr PYTHIA", "cele PYTHIA", "ccbar PYTHIA", - "bchadr PYTHIA", "bchadr PYTHIA", "bele PYTHIA", "bbbar PYTHIA" + const Char_t *label[2][5] = { + "chadr PYTHIA", "chadr PYTHIA", "cele PYTHIA", "ccbar PYTHIA", "ccbar PYTHIA+EvtGen", + "bchadr PYTHIA", "bchadr PYTHIA", "bele PYTHIA", "bbbar PYTHIA", "bbbar PYTHIA+EvtGen" }; Int_t iprocess = -1; @@ -22,15 +27,15 @@ GeneratorCustom(TString opt = "") else if(opt.Contains("bb")) {iprocess = 1; printf("-- HF process = %s -- \n",opt.Data());} else {printf("-- no HF process set = %s -- \n",opt.Data());} - TString optList[5] = {"had", "hv0", "ele", "all", "jetjet"}; + TString optList[6] = {"had", "hv0", "ele", "all", "bEvtGen", "jetjet"}; Int_t idecay = 0; - for (Int_t iopt = 0; iopt < 5; iopt++){ + for (Int_t iopt = 0; iopt < 6; iopt++){ if (opt.Contains(optList[iopt])) idecay = iopt; } AliGenerator *phf; - if(idecay < 4) { + if(idecay < 5) { Float_t randHF = gRandom->Rndm(); Int_t typeHF = -1; if(randHF < 0.9) { @@ -44,11 +49,18 @@ GeneratorCustom(TString opt = "") printf(">>>>> added generator for HF production: %s \n", label[iprocess][2]); } } - else if(idecay == 4) { + else if(idecay == 5) { phf = GeneratorPythia6Jets(kPythia6Tune_Perugia2011); ctl->AddGenerator(phf, "jetjet", 1.); printf(">>>>> added generator for HF production: %s \n", "jetjet"); } + + if(idecay == 4) { //EvtGen for b-decays + AliGenEvtGen *gener = new AliGenEvtGen(); + gener->SetParticleSwitchedOff(AliGenEvtGen::kBeautyPart); + ctl->AddGenerator(gener,"EvtGen",1.); + } + return ctl; }