]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - macros/ConfigPPR.C
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Phi range set consistently to 0 ... 360.
[u/mrichter/AliRoot.git] / macros / ConfigPPR.C
diff --git a/macros/ConfigPPR.C b/macros/ConfigPPR.C
index cdf6e2bbbc00cc06ee40d9af13f58ff544c0e705..b87c6cc9c945fb167bd060c4c14fb2f30e34f50d 100644 (file)
--- a/macros/ConfigPPR.C
+++ b/macros/ConfigPPR.C
@@ -4,8 +4,9 @@ enum PprRun_t
     kParam_8000,   kParam_4000,  kParam_2000, 
     kHijing_cent1, kHijing_cent2, 
     kHijing_per1,  kHijing_per2, kHijing_per3, kHijing_per4,  kHijing_per5,
-    kHijing_jj25,  kHijing_jj50, kHijing_jj75, kHijing_jj100, kHijing_jj125, 
-    kHijing_gj25,  kHijing_gj50, kHijing_gj75, kHijing_gj100, kHijing_gj125
+    kHijing_jj25,  kHijing_jj50, kHijing_jj75, kHijing_jj100, kHijing_jj200, 
+    kHijing_gj25,  kHijing_gj50, kHijing_gj75, kHijing_gj100, kHijing_gj200,
+    kHijing_pA
 };
 
 enum PprGeo_t 
@@ -47,7 +48,10 @@ void Config()
     // Set Random Number seed
     // gRandom->SetSeed(12345);
 
-    new AliGeant3("C++ Interface to Geant3");
+    // libraries required by geant321
+    gSystem->Load("libgeant321");
+
+    new TGeant3("C++ Interface to Geant3");
 
     if (!gSystem->Getenv("CONFIG_FILE"))
     {
@@ -167,7 +171,7 @@ void Config()
     Int_t   iTPC    = 1;
     Int_t   iTRD    = 1;
     Int_t   iZDC    = 1;
-    Int_t   iEMCAL  = 0;
+    Int_t   iEMCAL  = 1;
     Int_t   iVZERO  = 1;
     Int_t   iCRT    = 0;
 
@@ -411,7 +415,7 @@ void Config()
     if (iEMCAL)
     {
         //=================== EMCAL parameters ============================
-        AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCALArch1a");
+        AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "G56_2_55_19");
     }
 
      if (iCRT)
@@ -445,7 +449,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        comment = comment.Append(":HIJINGparam test 50 particles");
        AliGenHIJINGpara *gener = new AliGenHIJINGpara(50);
        gener->SetMomentumRange(0, 999999.);
-       gener->SetPhiRange(-180., 180.);
+       gener->SetPhiRange(0., 360.);
        // Set pseudorapidity range from -8 to 8.
        Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
        Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
@@ -455,7 +459,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        comment = comment.Append(":HIJINGparam N=8000");
        AliGenHIJINGpara *gener = new AliGenHIJINGpara(86030);
        gener->SetMomentumRange(0, 999999.);
-       gener->SetPhiRange(-180., 180.);
+       gener->SetPhiRange(0., 360.);
        // Set pseudorapidity range from -8 to 8.
        Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
        Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
@@ -465,7 +469,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        comment = comment.Append("HIJINGparam N=4000");
        AliGenHIJINGpara *gener = new AliGenHIJINGpara(43015);
        gener->SetMomentumRange(0, 999999.);
-       gener->SetPhiRange(-180., 180.);
+       gener->SetPhiRange(0., 360.);
        // Set pseudorapidity range from -8 to 8.
        Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
        Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
@@ -475,7 +479,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        comment = comment.Append("HIJINGparam N=2000");
        AliGenHIJINGpara *gener = new AliGenHIJINGpara(21507);
        gener->SetMomentumRange(0, 999999.);
-       gener->SetPhiRange(-180., 180.);
+       gener->SetPhiRange(0., 360.);
        // Set pseudorapidity range from -8 to 8.
        Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
        Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
@@ -539,10 +543,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(25.);
+       gener->SetPtJet(25.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.3,0.3);
+       gener->SetJetPhiRange(75., 165.);   
        break;
 
     case kHijing_jj50:
@@ -552,10 +557,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(50.);
+       gener->SetPtJet(50.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.3,0.3);
+       gener->SetJetPhiRange(75., 165.);   
        break;
 
     case kHijing_jj75:
@@ -565,10 +571,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(75.);
+       gener->SetPtJet(75.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.3,0.3);
+       gener->SetJetPhiRange(75., 165.);   
        break;
 
     case kHijing_jj100:
@@ -578,23 +585,25 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(100.);
+       gener->SetPtJet(100.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.3,0.3);
+       gener->SetJetPhiRange(75., 165.);   
        break;
 
-    case kHijing_jj125:
-       comment = comment.Append("HIJING Jet 125 GeV");
+    case kHijing_jj200:
+       comment = comment.Append("HIJING Jet 200 GeV");
        AliGenHijing *gener = HijingStandard();
 // impact parameter range
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(125.);
+       gener->SetPtJet(200.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.3,0.3);
+       gener->SetJetPhiRange(75., 165.);   
        break;
 //
 // Gamma-Jet
@@ -606,10 +615,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(25.);
+       gener->SetPtJet(25.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.12, 0.12);
+        gener->SetJetPhiRange(220., 320.);
        break;
 
     case kHijing_gj50:
@@ -619,10 +629,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(50.);
+       gener->SetPtJet(50.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.12, 0.12);
+        gener->SetJetPhiRange(220., 320.);
        break;
 
     case kHijing_gj75:
@@ -632,10 +643,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(75.);
+       gener->SetPtJet(75.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.12, 0.12);
+        gener->SetJetPhiRange(220., 320.);
        break;
 
     case kHijing_gj100:
@@ -645,23 +657,63 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(100.);
+       gener->SetPtJet(100.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.12, 0.12);
+        gener->SetJetPhiRange(220., 320.);
        break;
 
-    case kHijing_gj125:
-       comment = comment.Append("HIJING Gamma 125 GeV");
+    case kHijing_gj200:
+       comment = comment.Append("HIJING Gamma 200 GeV");
        AliGenHijing *gener = HijingStandard();
 // impact parameter range
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(125.);
+       gener->SetPtJet(200.);
        gener->SetRadiation(isw);
-       gener->JetEtaRange(-0.3,0.3);
-       gener->JetPhiRange(15.,105.);   
+       gener->SetSimpleJets(!isw);
+       gener->SetJetEtaRange(-0.12, 0.12);
+        gener->SetJetPhiRange(220., 320.);
+       break;
+    case kHijing_pA:
+       comment = comment.Append("HIJING pA");
+
+       AliGenCocktail *gener  = new AliGenCocktail();
+//     gener->SetTrackingFlag(0);
+       gener->SetOrigin(0, 0, 0);    // vertex position
+       gener->SetSigma(0, 0, 5.3);   // Sigma in (X,Y,Z) (cm) on IP position
+       gener->SetCutVertexZ(1.);     // Truncate at 1 sigma
+       gener->SetVertexSmear(kPerEvent); 
+       AliGenHijing   *hijing = new AliGenHijing(-1);
+// centre of mass energy 
+       hijing->SetEnergyCMS(TMath::Sqrt(82./208.) * 14000.);
+// impact parameter range
+       hijing->SetImpactParameterRange(0., 15.);
+// reference frame
+       hijing->SetReferenceFrame("CMS");
+       hijing->SetBoostLHC(1);
+// projectile
+       hijing->SetProjectile("P", 1, 1);
+       hijing->SetTarget    ("A", 208, 82);
+// tell hijing to keep the full parent child chain
+       hijing->KeepFullEvent();
+// enable jet quenching
+       hijing->SetJetQuenching(0);
+// enable shadowing
+       hijing->SetShadowing(1);
+// Don't track spectators
+       hijing->SetSpectators(0);
+// kinematic selection
+       hijing->SetSelectAll(0);
+//
+       AliGenSlowNucleons*  gray    = new AliGenSlowNucleons(1);
+       AliSlowNucleonModel* model   = new AliSlowNucleonModelExp();
+       gray->SetSlowNucleonModel(model);
+       gray->SetDebug(1);
+       gener->AddGenerator(hijing,"Hijing pPb", 1);
+       gener->AddGenerator(gray,  "Gray Particles",1);
        break;
     }
     return gener;