Update Toy Model:
authormiweber <m.weber@cern.ch>
Thu, 9 Oct 2014 14:43:48 +0000 (16:43 +0200)
committermiweber <m.weber@cern.ch>
Thu, 9 Oct 2014 14:43:48 +0000 (16:43 +0200)
- adding option to have a flat eta distribution for randomly produced particles
- default switched OFF

PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.cxx
PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.h

index cdd44ee..be5043a 100755 (executable)
@@ -55,6 +55,7 @@ AliAnalysisTaskToyModel::AliAnalysisTaskToyModel()
   fPtMin(0.0), fPtMax(100.0),
   fEtaMin(-1.0), fEtaMax(1.0),
   fSigmaGaussEta(4.0),
+  fConstantEta(-1.),
   fFixPt(-1.),
   fFixedPositiveRatio(kFALSE),
   fUseAcceptanceParameterization(kFALSE), fAcceptanceParameterization(0),
@@ -675,8 +676,15 @@ void AliAnalysisTaskToyModel::Run(Int_t nEvents) {
       if(fUseDebug) 
        Printf("Generating positive: %d(%d)",iParticleCount+1,nGeneratedPositive);
 
+      // use a constant distribution of particles in eta in a certain range, -fConstantEta - +fConstantEta
+      if(fConstantEta>0){
+       vEta = gRandom->Uniform(-fConstantEta,fConstantEta);
+      }
+
       //Pseudo-rapidity sampled from a Gaussian centered @ 0
-      vEta = gRandom->Gaus(0.0,fSigmaGaussEta);
+      else{
+       vEta = gRandom->Gaus(0.0,fSigmaGaussEta);
+      }
 
       //Fill QA histograms (full phase space)
       fHistEtaTotal->Fill(vEta);
index cb668e0..d3cc332 100755 (executable)
@@ -57,6 +57,9 @@ class AliAnalysisTaskToyModel : public TObject {
   void SetSigmaGaussEta(Double_t sigmaGaussEta){
     fSigmaGaussEta = sigmaGaussEta;
   }
+  void SetConstantEta(Double_t constantEta){
+    fConstantEta = constantEta;
+  }
   void SetFixPt(Double_t fixPt){
     fFixPt = fixPt;
   } 
@@ -225,6 +228,7 @@ class AliAnalysisTaskToyModel : public TObject {
   Double_t fEtaMin; //eta min for acceptance
   Double_t fEtaMax; //eta max for acceptance
   Double_t fSigmaGaussEta; //sigma for the Gaussian distribution of randomly produced particles (default = 4.0)
+  Double_t fConstantEta; //eta value for a constant distribution of particles, if -1. then Gauss is used (default = -1)
   Double_t fFixPt; //fixed pT for unidentified particles (default = -1., i.e. randomly produced pT from fPtSpectraAllCharges)
   Bool_t fFixedPositiveRatio; //fix ratio of produced positive to negative particles (dafult = kFALSE, i.e. randomly produced ratio)