Merge pull request #313 from gluparel/master

Add Config for XiC production in pp collisions with Pythia6
alisw · Sep 24, 2019 · ffe5f26 · ffe5f26
2 parents 4e7fffb + 10a82b4
commit ffe5f26
Show file tree

Hide file tree

Showing 2 changed files with 22 additions and 1 deletion.
diff --git a/MC/CustomDecayTables/XiC/decaytable.dat b/MC/CustomDecayTables/XiC/decaytable.dat
@@ -274,7 +274,7 @@
               1    0    0.224000      2212      -321       211         0         0
               1    0    0.276000      2212      -313         0         0         0     
       4132  Xi_c0             Xi_cbar0          0  0  1     2.4710      0.00000     0.00000  3.36000E-02  0  1	
-              1    0    0.100000       211      3312         0         0         0     
+              1    0    1.000000       211      3312         0         0         0     
        111  pi0                                 0  0  0     0.13498    0.000000    0.00000  2.56000E-05  0  1
               1    0    0.988230        22        22         0         0         0
               1    2    0.011740        22       -11        11         0         0

diff --git a/MC/CustomGenerators/PWGHF/Pythia6_pp_XiC.C b/MC/CustomGenerators/PWGHF/Pythia6_pp_XiC.C
@@ -0,0 +1,21 @@
+AliGenerator *GeneratorCustom(TString opt = "")
+{
+  Int_t iprocess = uidConfig % 4;
+
+  Int_t particleTrigger[4]={4232, -4232, 4132, -4132};
+  Int_t iparticle =0;
+  if (iprocess>1) iparticle=1;
+
+  TString label[2]={"XiCPlus to hadrons", "XiC0 to hadrons"};
+
+  AliGenCocktail *ctl = static_cast<AliGenCocktail*>(GeneratorCocktail("Perugia2011_HF"));			
+  AliGenPythia *pythia = (AliGenPythia*) GeneratorPythia6(kPythia6Tune_Perugia2011);
+  pythia->SetProcess(kPyCharmppMNRwmi);
+  pythia->SetHeavyQuarkYRange(-1.5, 1.5); 
+  pythia->SetTriggerY(1.0);
+  pythia->SetTriggerParticle(particleTrigger[iprocess],1.5,1.,100); //Xic //pdgcode,eta,ptmin,ptmax
+
+  ctl->AddGenerator(pythia, Form("MC %s",label[iparticle].Data()),1.);
+
+  return ctl;
+}