]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EMCAL/Config.C
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Correcting sources of warnings
[u/mrichter/AliRoot.git] / EMCAL / Config.C
diff --git a/EMCAL/Config.C b/EMCAL/Config.C
index a4e5a81db380c4ef196e05fcc34fb099ae1afd80..89e0853231a8b90edccf2780e46e540f5cf193fa 100644 (file)
--- a/EMCAL/Config.C
+++ b/EMCAL/Config.C
@@ -6,8 +6,11 @@ enum PprRun_t
     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, 
-    kJetPlusBg,    kGammaPlusBg, kParam_8000_Ecal, 
-    kJets_100,
+    kJetPlusBg,    kGammaPlusBg, 
+    kParam_8000_Ecal, kParam_4000_Ecal, 
+    kJets_50,       kJets_75,      kJets_100,      kJets_200,
+    kGammaJets_50,  kGammaJets_75, kGammaJets_100, kGammaJets_200,
+    kGammaJets_250, kGammaJets_300,
     kGammaGun, kGammaBox
 };
 
@@ -22,10 +25,18 @@ enum PprRad_t
 };
 
 // This part for configuration    
-//static PprRun_t run = kParam_8000_Ecal;
-static PprRun_t run = kJets_100;
-static PprGeo_t geo = kHoles;
-static PprRad_t rad = kNoGluonRadiation;
+static PprRun_t run          = kHijing_jj75;
+//static PprRun_t run = kJets_100;
+static PprGeo_t geo          = kHoles;
+static PprRad_t rad          = kGluonRadiation;
+static Int_t    eventsPerRun = 100;
+static Int_t    simpleJet    = 1;
+if (rad == kGluonRadiation) 
+{
+    simpleJet = 0;
+}
+
+
 // Comment line 
 static TString  comment;
 
@@ -103,6 +114,7 @@ void Config()
     gAlice->SetDebug(1);
 
     AliGenerator*  gener = GeneratorFactory(run);
+    gener->SetPtRange(0.,1.e10);
     gener->Init();
     gener->SetTrackingFlag(1);
     
@@ -134,26 +146,26 @@ void Config()
     field->SetSolenoidField(4.);
     gAlice->SetField(field);    
 
-    Int_t iMAG    = 0;
-    Int_t iITS    = 0;
-    Int_t iTPC    = 0;
-    Int_t iTOF    = 0;
+    Int_t iMAG    = 1;
+    Int_t iITS    = 1;
+    Int_t iTPC    = 1;
+    Int_t iTOF    = 1;
     Int_t iRICH   = 0;
     Int_t iZDC    = 0;
     Int_t iCASTOR = 0;
-    Int_t iTRD    = 0;
-    Int_t iABSO   = 0;
+    Int_t iTRD    = 1;
+    Int_t iABSO   = 1;
     Int_t iDIPO   = 0;
-    Int_t iHALL   = 0;
-    Int_t iFRAME  = 0;
+    Int_t iHALL   = 1;
+    Int_t iFRAME  = 1;
     Int_t iSHIL   = 0;
-    Int_t iPIPE   = 0;
-    Int_t iFMD    = 0;
+    Int_t iPIPE   = 1;
+    Int_t iFMD    = 1;
     Int_t iMUON   = 0;
-    Int_t iPHOS   = 0;
-    Int_t iPMD    = 0;
-    Int_t iSTART  = 0;
-    Int_t iVZERO  = 0;
+    Int_t iPHOS   = 1;
+    Int_t iPMD    = 1;
+    Int_t iSTART  = 1;
+    Int_t iVZERO  = 1;
     Int_t iEMCAL  = 1;    
 // 
 //  enable/disable StepManager()
@@ -464,7 +476,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetMomentum(100.);
        gener->SetPhi(240.);
        gener->SetTheta(91.);
-       gener->SetPart(kElectron);
+       gener->SetPart(kGamma);
      break;
 
     case kGammaBox:
@@ -472,7 +484,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetMomentumRange(100., 100.1);
        gener->SetPhiRange(0,360);
        gener->SetThetaRange(45., 135.);
-       gener->SetPart(kElectron);
+       gener->SetPart(kGamma);
        break;
        
     case kTest50:
@@ -516,6 +528,7 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
        gener->SetThetaRange(thmin,thmax);
        break;
+
     case kParam_8000_Ecal:
        //coment= comment.Append(":HIJINGparam N=8000 baryons");
        AliGenHIJINGparaBa *gener = new AliGenHIJINGparaBa(82534);
@@ -526,6 +539,17 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        Float_t thmax = EtaToTheta(-5);  // theta max. <---> eta min 
        gener->SetThetaRange(thmin,thmax);
        break;
+
+    case kParam_4000_Ecal:
+       //coment= comment.Append(":HIJINGparam N=8000 baryons");
+       AliGenHIJINGparaBa *gener = new AliGenHIJINGparaBa(82534/2.);
+       gener->SetMomentumRange(0, 999999.);
+       gener->SetPhiRange(-180., 180.);
+       // Set pseudorapidity range from -8 to 8.
+       Float_t thmin = EtaToTheta( 5);   // theta min. <---> eta max
+       Float_t thmax = EtaToTheta(-5);  // theta max. <---> eta min 
+       gener->SetThetaRange(thmin,thmax);
+       break;
 //
 //  Hijing Central
 //
@@ -584,7 +608,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(25.);
+       gener->SetPtJet(25.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -597,7 +622,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(50.);
+       gener->SetPtJet(50.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -610,7 +636,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(75.);
+       gener->SetPtJet(75.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -620,10 +647,11 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        //coment= comment.Append("HIJING Jet 100 GeV");
        AliGenHijing *gener = HijingStandard();
 // impact parameter range
-       gener->SetImpactParameterRange(10., 10.1);
+       gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(100.);
+       gener->SetPtJet(100.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -636,7 +664,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(1);
-       gener->SetPtMinJet(125.);
+       gener->SetPtJet(125.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -651,7 +680,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(25.);
+       gener->SetPtJet(25.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -664,7 +694,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(50.);
+       gener->SetPtJet(50.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -677,7 +708,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(75.);
+       gener->SetPtJet(75.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -690,7 +722,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(100.);
+       gener->SetPtJet(100.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -703,7 +736,8 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
        gener->SetImpactParameterRange(0., 5.);
        // trigger
        gener->SetTrigger(2);
-       gener->SetPtMinJet(125.);
+       gener->SetPtJet(125.);
+       gener->SetSimpleJets(simpleJet);
        gener->SetRadiation(isw);
        gener->SetJetEtaRange(-0.3,0.3);
        gener->SetJetPhiRange(15.,105.);   
@@ -743,8 +777,9 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
 //
 //      Add all to cockail ...    
 //
+       gener->AddGenerator(jets,"Jets",1);
        gener->AddGenerator(bg,"Underlying Event", 1);
-//     gener->AddGenerator(jets,"Jets",1);
+
      break;    
     case kGammaPlusBg:
        AliGenCocktail *gener = new AliGenCocktail();
@@ -782,13 +817,104 @@ AliGenerator* GeneratorFactory(PprRun_t run) {
 //
 //      Add all to cockail ...    
 //
-       gener->AddGenerator(bg,"Underlying Event", 1);
        gener->AddGenerator(jets,"Jets",1);
+       gener->AddGenerator(bg,"Underlying Event", 1);
+
        break;
+    case kJets_50:
+//  50 GeV Jets
+        AliGenPythia *gener = PythiaJets(50.);
+        break; 
+    case kJets_75:
+//  75 GeV Jets
+        AliGenPythia *gener = PythiaJets(75.);
+        break; 
     case kJets_100:
 //  100 GeV Jets  
        AliGenPythia *gener = PythiaJets(100.);
-       break;
+       break; 
+    case kJets_200:
+//  200 GeV Jets
+        AliGenPythia *gener = PythiaJets(200.);
+        break;
+    case kGammaJets_50:
+//  50 GeV Jets + Gamma
+        AliGenPythia *gener = PythiaJets(-1);
+        gener->SetEnergyCMS(5500.);
+        gener->SetProcess(kPyDirectGamma);
+        gener->SetJetEtaRange(-0.3,+0.3);
+        gener->SetJetPhiRange(15.,105.);
+        gener->SetGammaEtaRange(-0.12, 0.12);
+        gener->SetGammaPhiRange(220., 320.);
+        gener->SetStrucFunc(kGRV_LO_98);
+        gener->SetPtHard(50.,50.001);
+        gener->SetForceDecay(kAll);
+        break;
+    case kGammaJets_75:
+//  75 GeV Jets + Gamma 
+        AliGenPythia *gener = PythiaJets(-1);
+        gener->SetEnergyCMS(5500.);
+        gener->SetProcess(kPyDirectGamma);
+        gener->SetJetEtaRange(-0.3,+0.3);
+        gener->SetJetPhiRange(15.,105.);
+        gener->SetGammaEtaRange(-0.12, 0.12);
+        gener->SetGammaPhiRange(220., 320.);
+        gener->SetStrucFunc(kGRV_LO_98);
+        gener->SetPtHard(75.,75.001);
+        gener->SetForceDecay(kAll);
+        break; 
+    case kGammaJets_100:
+// 100 GeV Jets + Gamma
+        AliGenPythia *gener = PythiaJets(-1);
+        gener->SetEnergyCMS(5500.);
+        gener->SetProcess(kPyDirectGamma);
+        gener->SetJetEtaRange(-0.3,+0.3);
+        gener->SetJetPhiRange(15.,105.);
+        gener->SetGammaEtaRange(-0.12, 0.12);
+        gener->SetGammaPhiRange(220., 320.);
+        gener->SetStrucFunc(kGRV_LO_98);
+        gener->SetPtHard(100.,100.001);
+        gener->SetForceDecay(kAll);
+        break; 
+    case kGammaJets_200:
+//  200 GeV Jets + Gamma
+        AliGenPythia *gener = PythiaJets(-1);
+        gener->SetEnergyCMS(5500.);
+        gener->SetProcess(kPyDirectGamma);
+        gener->SetJetEtaRange(-0.3,+0.3);
+        gener->SetJetPhiRange(15.,105.);
+        gener->SetGammaEtaRange(-0.12, 0.12);
+        gener->SetGammaPhiRange(220., 320.);
+        gener->SetStrucFunc(kGRV_LO_98);
+        gener->SetPtHard(200.,200.001);
+        gener->SetForceDecay(kAll);
+        break; 
+    case kGammaJets_250:
+//  250 GeV Jets + Gamma
+        AliGenPythia *gener = PythiaJets(-1);
+        gener->SetEnergyCMS(5500.);
+        gener->SetProcess(kPyDirectGamma);
+        gener->SetJetEtaRange(-0.3,+0.3);
+        gener->SetJetPhiRange(15.,105.);
+        gener->SetGammaEtaRange(-0.12, 0.12);
+        gener->SetGammaPhiRange(220., 320.);
+        gener->SetStrucFunc(kGRV_LO_98);
+        gener->SetPtHard(250.,250.001);
+        gener->SetForceDecay(kAll);
+        break; 
+    case kGammaJets_300:
+//  300 GeV Jets + Gamma
+        AliGenPythia *gener = PythiaJets(-1);
+        gener->SetEnergyCMS(5500.);
+        gener->SetProcess(kPyDirectGamma);
+        gener->SetJetEtaRange(-0.3,+0.3);
+        gener->SetJetPhiRange(15.,105.);
+        gener->SetGammaEtaRange(-0.12, 0.12);
+        gener->SetGammaPhiRange(220., 320.);
+        gener->SetStrucFunc(kGRV_LO_98);
+        gener->SetPtHard(300.,300.001);
+        gener->SetForceDecay(kAll);
+        break;
     }
 
     return gener;