Pycell configurable.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Apr 2003 20:39:00 +0000 (20:39 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Apr 2003 20:39:00 +0000 (20:39 +0000)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliGenPythia.h

index 258b2e86d73571eb0b692aa70578c8219cffe5e0..0cd83dbfaaeea6ffa0b701facd381913c289642a 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.3  2003/03/31 15:47:16  morsch
+Unhandled cases added.
+
 Revision 1.2  2003/03/27 09:28:38  morsch
 Charm production process configuration for pPb and pp added
 (A. Dainese, N. Carrer)
@@ -298,6 +301,8 @@ AliGenPythia::AliGenPythia(Int_t npart)
     fFragmentation = kTRUE;
     // Default counting mode
     fCountMode = kCountAll;
+    // Pycel
+    SetPycellParameters();
 }
 
 AliGenPythia::AliGenPythia(const AliGenPythia & Pythia)
@@ -311,6 +316,21 @@ AliGenPythia::~AliGenPythia()
 // Destructor
 }
 
+void AliGenPythia::SetPycellParameters(Float_t etamax, Int_t neta, Int_t nphi,
+                                      Float_t thresh, Float_t etseed, Float_t minet, Float_t r)
+{
+// Set pycell parameters
+    fPycellEtaMax    =  etamax;
+    fPycellNEta      =  neta;
+    fPycellNPhi      =  nphi;
+    fPycellThreshold =  thresh;
+    fPycellEtSeed    =  etseed;
+    fPycellMinEtJet  =  minet;
+    fPycellMaxRadius =  r;
+}
+
+
+
 void AliGenPythia::SetEventListRange(Int_t eventFirst, Int_t eventLast)
 {
   // Set a range of event numbers, for which a table
@@ -418,6 +438,23 @@ void AliGenPythia::Init()
        break;
     }
 //
+//
+//  JetFinder for Trigger
+//
+//  Configure detector (EMCAL like)
+//
+       fPythia->SetPARU(51, fPycellEtaMax);
+       fPythia->SetMSTU(51, fPycellNEta);
+       fPythia->SetMSTU(52, fPycellNPhi);
+//
+//  Configure Jet Finder
+//  
+       fPythia->SetPARU(58,  fPycellThreshold);
+       fPythia->SetPARU(52,  fPycellEtSeed);
+       fPythia->SetPARU(53,  fPycellMinEtJet);
+       fPythia->SetPARU(54,  fPycellMaxRadius);
+       fPythia->SetMSTU(54,  2);
+//
 //  This counts the total number of calls to Pyevnt() per run.
     fTrialsRun = 0;
     fQ         = 0.;
@@ -888,19 +925,6 @@ void AliGenPythia::RecJetsUA1(Float_t eCellMin, Float_t eCellSeed, Float_t eMin,
 //  Calls the Pythia jet finding algorithm to find jets in the current event
 //
 //
-//  Configure detector (EMCAL like)
-//
-    fPythia->SetPARU(51,2.);
-    fPythia->SetMSTU(51,Int_t(96 * 2./0.7));
-    fPythia->SetMSTU(52,3 * 144);
-//
-//  Configure Jet Finder
-//  
-    fPythia->SetPARU(58, eCellMin);
-    fPythia->SetPARU(52, eCellSeed);
-    fPythia->SetPARU(53, eMin);
-    fPythia->SetPARU(54, rMax);
-    fPythia->SetMSTU(54, 2);
 //
 //  Save jets
     Int_t n     = fPythia->GetN();
@@ -946,20 +970,7 @@ void  AliGenPythia::GetJets(Int_t& nJets, Int_t& nJetsTrig, Float_t jets[4][10])
 //  Loading jets from common block
 //
     } else {
-//
-//  Configure detector (EMCAL like)
-//
-       fPythia->SetPARU(51,2.);
-       fPythia->SetMSTU(51,Int_t(96 * 2./0.7));
-       fPythia->SetMSTU(52,3 * 144);
-//
-//  Configure Jet Finder
-//  
-       fPythia->SetPARU(58,  0.0);
-       fPythia->SetPARU(52,  4.0);
-       fPythia->SetPARU(53, 10.0);
-       fPythia->SetPARU(54,  1.0);
-       fPythia->SetMSTU(54,  2);
+
 //
 //  Run Jet Finder
        fPythia->Pycell(nJets);
index 34a9e34667bc52083656ae02c6a3194f23bb5001..f9610c571a3ec4b75d7bf859a7235b4b543bf97f 100644 (file)
@@ -90,7 +90,11 @@ class AliGenPythia : public AliGenMC
     // get triggered jets
     void GetJets(Int_t& njets, Int_t& ntrig, Float_t[4][10]);
     void RecJetsUA1(Float_t eCellMin, Float_t eCellSeed, Float_t eMin, Float_t rMin, 
-                   Int_t& njets, Float_t jets[4][50]);    
+                   Int_t& njets, Float_t jets[4][50]);
+    void SetPycellParameters(Float_t etamax = 2., Int_t neta = 274, Int_t nphi = 432,
+                            Float_t thresh = 0., Float_t etseed = 4.,
+                            Float_t minet = 10., Float_t r = 1.);
+    
     void LoadEvent();
     // Getters
     virtual Process_t    GetProcess() {return fProcess;}
@@ -159,7 +163,13 @@ class AliGenPythia : public AliGenMC
     Float_t     fEtaMaxGamma;       // Maximum eta of triggered gamma
     Float_t     fPhiMinGamma;       // Minimum phi of triggered gamma
     Float_t     fPhiMaxGamma;       // Maximum phi of triggered gamma
-
+    Float_t     fPycellEtaMax;      // Max. eta for Pycell 
+    Int_t       fPycellNEta;        // Number of eta bins for Pycell 
+    Int_t       fPycellNPhi;        // Number of phi bins for Pycell
+    Float_t    fPycellThreshold;   // Pycell threshold
+    Float_t    fPycellEtSeed;      // Pycell seed
+    Float_t    fPycellMinEtJet;    // Pycell min. jet et
+    Float_t    fPycellMaxRadius;   // Pycell cone radius
     StackFillOpt_t fStackFillOpt;   // Stack filling with all particles with
                                     // that flavour or only with selected
                                     // parents and their decays