AliGenMUONCocktailpp is used with modified values of the
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 1 Feb 2009 10:59:34 +0000 (10:59 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 1 Feb 2009 10:59:34 +0000 (10:59 +0000)
different cross-sections, hidden (as much as possible) from the people
who will analyze afterwards the results, in order to see the accuracy
of the physics analysis when unknown cross sections are at the input of
the Monte-Carlo. By default the hardwired values are used, and only if
intentionally some env variables are defined the values will be
changed, without print-out in the log file (as default).
(Bogdan Vulpescu)

EVGEN/AliGenMUONCocktailpp.cxx

index 476c3c6..45d0107 100644 (file)
@@ -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); 
 
 //-------------------------------------------------------------------