AliGenMUONCocktailpp is used with modified values of the

[u/mrichter/AliRoot.git] / EVGEN / AliGenMUONCocktailpp.cxx

diff --git a/EVGEN/AliGenMUONCocktailpp.cxx b/EVGEN/AliGenMUONCocktailpp.cxx
index 476c3c67432064da7d64ffaec1b69810262d781c..45d0107a9b754f75641cb685dbab4c281eafbbd6 100644 (file)
--- a/EVGEN/AliGenMUONCocktailpp.cxx
+++ b/EVGEN/AliGenMUONCocktailpp.cxx
@@ -107,7 +107,50 @@ void AliGenMUONCocktailpp::CreateCocktail()
     Double_t sigmaupsilonPP = 0.228e-6;
     Double_t sigmaccbar = 11.2e-3;
     Double_t sigmabbbar = 0.51e-3;
-    
+
+    Bool_t showSIGJPSI = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGJPSI")) {
+      sigmajpsi = atof(gSystem->Getenv("COCKTAIL_SIGJPSI"));
+      AliInfo("Cross-section for JPsi set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGJPSI = kFALSE;
+    }
+    Bool_t showSIGPSIP = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGPSIP")) {
+      sigmapsiP = atof(gSystem->Getenv("COCKTAIL_SIGPSIP"));
+      AliInfo("Cross-section for Psi-prime set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGPSIP = kFALSE;
+    }
+    Bool_t showSIGUPSI = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGUPSI")) {
+      sigmaupsilon = atof(gSystem->Getenv("COCKTAIL_SIGUPSI"));
+      AliInfo("Cross-section for Upsilon 1S set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGUPSI = kFALSE;
+    }
+    Bool_t showSIGUPSIP = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGUPSIP")) {
+      sigmaupsilonP = atof(gSystem->Getenv("COCKTAIL_SIGUPSIP"));
+      AliInfo("Cross-section for Upsilon 2S set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGUPSIP = kFALSE;
+    }
+    Bool_t showSIGUPSIPP = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGUPSIPP")) {
+      sigmaupsilonPP = atof(gSystem->Getenv("COCKTAIL_SIGUPSIPP"));
+      AliInfo("Cross-section for Upsilon 3S set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGUPSIPP = kFALSE;
+    }
+    Bool_t showSIGCCBAR = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGCCBAR")) {
+      sigmaccbar = atof(gSystem->Getenv("COCKTAIL_SIGCCBAR"));
+      AliInfo("Cross-section for c-cbar set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGCCBAR = kFALSE;
+    }
+    Bool_t showSIGBBBAR = kTRUE;
+    if (gSystem->Getenv("COCKTAIL_SIGBBBAR")) {
+      sigmabbbar = atof(gSystem->Getenv("COCKTAIL_SIGBBBAR"));
+      AliInfo("Cross-section for b-bbar set from external value");
+      if (!gSystem->Getenv("COCKTAIL_SIGSHOW")) showSIGBBBAR = kFALSE;
+    }
+   
     AliInfo(Form("the parametrised resonances uses the decay mode %d",fDecayModeResonance));
 
 // Generation using CDF scaled distribution for pp@14TeV (from D. Stocco)
@@ -123,8 +166,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
 
     genjpsi->Init();  // generation in 4pi
     ratiojpsi = sigmajpsi / sigmaReaction * genjpsi->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax); // get weight
-    AliInfo(Form("jpsi prod. cross-section in pp(14 TeV) %5.3g b",sigmajpsi));
-    AliInfo(Form("jpsi prod. probability per collision in acceptance %5.3g",ratiojpsi));
+    if (showSIGJPSI) {
+      AliInfo(Form("jpsi prod. cross-section in pp(14 TeV) %5.3g b",sigmajpsi));
+      AliInfo(Form("jpsi prod. probability per collision in acceptance %5.3g",ratiojpsi));
+    }
 // second step: generation in selected kinematical range
     genjpsi->SetPtRange(ptMin, ptMax);  
     genjpsi->SetYRange(yMin, yMax);
@@ -143,8 +188,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
 
     genpsiP->Init();  // generation in 4pi
     ratiopsiP = sigmapsiP / sigmaReaction * genpsiP->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
-    AliInfo(Form("psiP prod. cross-section in pp(14 TeV) %5.3g b",sigmapsiP));
-    AliInfo(Form("psiP prod. probability per collision in acceptance %5.3g",ratiopsiP));
+    if (showSIGPSIP) {
+      AliInfo(Form("psiP prod. cross-section in pp(14 TeV) %5.3g b",sigmapsiP));
+      AliInfo(Form("psiP prod. probability per collision in acceptance %5.3g",ratiopsiP));
+    }
 // second step: generation in selected kinematical range
     genpsiP->SetPtRange(ptMin, ptMax);  
     genpsiP->SetYRange(yMin, yMax);
@@ -162,8 +209,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
     if (!gMC) genupsilon->SetDecayer(fDecayer);
     genupsilon->Init();  // generation in 4pi
     ratioupsilon = sigmaupsilon / sigmaReaction * genupsilon->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
-    AliInfo(Form("Upsilon 1S prod. cross-section in pp(14 TeV) %5.3g b",sigmaupsilon));
-    AliInfo(Form("Upsilon 1S prod. probability per collision in acceptance %5.3g",ratioupsilon));
+    if (showSIGUPSI) {
+      AliInfo(Form("Upsilon 1S prod. cross-section in pp(14 TeV) %5.3g b",sigmaupsilon));
+      AliInfo(Form("Upsilon 1S prod. probability per collision in acceptance %5.3g",ratioupsilon));
+    }
 // second step: generation in selected kinematical range
     genupsilon->SetPtRange(ptMin, ptMax);  
     genupsilon->SetYRange(yMin, yMax);
@@ -181,8 +230,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
     if (!gMC) genupsilonP->SetDecayer(fDecayer);  
     genupsilonP->Init();  // generation in 4pi
     ratioupsilonP = sigmaupsilonP / sigmaReaction * genupsilonP->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
-    AliInfo(Form("Upsilon 2S prod. cross-section in pp(14 TeV) %5.3g b",sigmaupsilonP));
-    AliInfo(Form("Upsilon 2S prod. probability per collision in acceptance %5.3g",ratioupsilonP));
+    if (showSIGUPSIP) {
+      AliInfo(Form("Upsilon 2S prod. cross-section in pp(14 TeV) %5.3g b",sigmaupsilonP));
+      AliInfo(Form("Upsilon 2S prod. probability per collision in acceptance %5.3g",ratioupsilonP));
+    }
 // second step: generation in the kinematical range
     genupsilonP->SetPtRange(ptMin, ptMax);  
     genupsilonP->SetYRange(yMin, yMax);
@@ -201,8 +252,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
     if (!gMC) genupsilonPP->SetDecayer(fDecayer);  
     genupsilonPP->Init();  // generation in 4pi
     ratioupsilonPP = sigmaupsilonPP / sigmaReaction * genupsilonPP->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
-    AliInfo(Form("Upsilon 3S prod. cross-section in pp(14 TeV) %5.3g b",sigmaupsilonPP));
-    AliInfo(Form("Upsilon 3S prod. probability per collision in acceptance %5.3g",ratioupsilonPP));
+    if (showSIGUPSIPP) {
+      AliInfo(Form("Upsilon 3S prod. cross-section in pp(14 TeV) %5.3g b",sigmaupsilonPP));
+      AliInfo(Form("Upsilon 3S prod. probability per collision in acceptance %5.3g",ratioupsilonPP));
+    }
 // second step: generation in selected kinematical range
     genupsilonPP->SetPtRange(ptMin, ptMax);  
     genupsilonPP->SetYRange(yMin, yMax);
@@ -219,6 +272,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
     ratioccbar = sigmaccbar/sigmaReaction;
     if (!gMC) gencharm->SetDecayer(fDecayer);  
     gencharm->Init();
+    if (showSIGCCBAR) {
+      AliInfo(Form("c-cbar prod. cross-section in pp(14 TeV) %5.3g b",sigmaccbar));
+      AliInfo(Form("c-cbar prod. probability per collision in acceptance %5.3g",ratioccbar));
+    }
     AddGenerator(gencharm,"CorrHFCharm",ratioccbar);
 
 //------------------------------------------------------------------
@@ -230,6 +287,10 @@ void AliGenMUONCocktailpp::CreateCocktail()
     ratiobbbar = sigmabbbar/sigmaReaction;
     if (!gMC) genbeauty->SetDecayer(fDecayer);  
     genbeauty->Init();
+    if (showSIGBBBAR) {
+      AliInfo(Form("b-bbar prod. cross-section in pp(14 TeV) %5.3g b",sigmabbbar));
+      AliInfo(Form("b-bbar prod. probability per collision in acceptance %5.3g",ratiobbbar));
+    }
     AddGenerator(genbeauty,"CorrHFBeauty",ratiobbbar); 
 
 //-------------------------------------------------------------------