quenching option 3: modified splittings
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 1 Oct 2006 17:53:50 +0000 (17:53 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 1 Oct 2006 17:53:50 +0000 (17:53 +0000)
PYTHIA6/AliGenPythia.cxx

index c7559da..783b434 100644 (file)
@@ -520,6 +520,16 @@ void AliGenPythia::Init()
     if (fQuench) {
        fPythia->InitQuenching(0., 0.1, 0.6e6, 0);
     }
+    if (fQuench == 3) {
+       // Nestor's change of the splittings
+       fPythia->SetPARJ(200, 0.8);
+       fPythia->SetMSTJ(41, 1);  // QCD radiation only
+       fPythia->SetMSTJ(42, 2);  // angular ordering
+       fPythia->SetMSTJ(44, 2);  // option to run alpha_s
+       fPythia->SetMSTJ(47, 0);  // No correction back to hard scattering element
+       fPythia->SetMSTJ(50, 0);  // No coherence in first branching
+       fPythia->SetPARJ(82, 1.); // Cut off for parton showers
+    }
 }
 
 void AliGenPythia::Generate()
@@ -580,7 +590,10 @@ void AliGenPythia::Generate()
            fPythia->Quench();
        } else if (fQuench == 2){
            fPythia->Pyquen(208., 0, 0.);
+       } else if (fQuench == 3) {
+           // Quenching is via multiplicative correction of the splittings
        }
+       
 //
 // Switch hadronisation on
 //