jet v2 analysis : give user control over random seed in toy mc event generator ...

diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.cxx b/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.cxx
index 6a62ffa2fb7af4c34ecdcf25aab7de923701f9df..760d7e43c30bad24b8ea3a51cae54fe0a5da4175 100644 (file)
--- a/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.cxx
+++ b/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.cxx
@@ -39,7 +39,7 @@ AliAnalysisTaskJetFlowMC::AliAnalysisTaskJetFlowMC() : AliAnalysisTaskSE("AliAna
     }
 }
 //_____________________________________________________________________________
-AliAnalysisTaskJetFlowMC::AliAnalysisTaskJetFlowMC(const char *name) : AliAnalysisTaskSE(name), fTracksOutName("JetFlowMC"), fTracksInName("PicoTrack"), fTracksIn(0), fTracksOut(0), fReuseTracks(kFALSE), fMult(2200), fCenBin(-1), fCentralityClasses(0), fFuncVn(0), fOutputList(0), fTrackSpectrum(0), fRandomizeEta(kTRUE), fJetSpectrumSF(0), fNoOfSFJets(0), fHistIntV2(0), fHistIntV3(0), fFlowFluctuations(-10), fMaxNumberOfIterations(100), fPsi2(-10), fPsi3(-10), fPrecisionPhi(1e-10), fDetectorType(kVZEROC), fHistSFJetSpectrum(0), fHistSFJetEtaPhi(0) {
+AliAnalysisTaskJetFlowMC::AliAnalysisTaskJetFlowMC(const char *name, Int_t seed) : AliAnalysisTaskSE(name), fTracksOutName("JetFlowMC"), fTracksInName("PicoTrack"), fTracksIn(0), fTracksOut(0), fReuseTracks(kFALSE), fMult(2200), fCenBin(-1), fCentralityClasses(0), fFuncVn(0), fOutputList(0), fTrackSpectrum(0), fRandomizeEta(kTRUE), fJetSpectrumSF(0), fNoOfSFJets(0), fHistIntV2(0), fHistIntV3(0), fFlowFluctuations(-10), fMaxNumberOfIterations(100), fPsi2(-10), fPsi3(-10), fPrecisionPhi(1e-10), fDetectorType(kVZEROC), fHistSFJetSpectrum(0), fHistSFJetEtaPhi(0) {
     // constructor
     DefineInput(0, TChain::Class());
     DefineOutput(1, TList::Class());
@@ -56,6 +56,11 @@ AliAnalysisTaskJetFlowMC::AliAnalysisTaskJetFlowMC(const char *name) : AliAnalys
         fHistToyDeltaPhi[i]             = 0x0;
         fHistToyVn[i]                   = 0x0;
     }
+    // by default construction replace gRandom
+    if(seed > -1) {
+        if(gRandom) delete gRandom;
+        gRandom = new TRandom3(seed);
+    }
 }
 //_____________________________________________________________________________
 AliAnalysisTaskJetFlowMC::~AliAnalysisTaskJetFlowMC()

diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.h b/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.h
index bf5473dbed29c0bc60d4cc919a93df3e55264758..77278a758cff479f6c38df3ecb949db589a1c8a0 100644 (file)
--- a/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.h
+++ b/PWG/FLOW/Tasks/AliAnalysisTaskJetFlowMC.h
@@ -20,7 +20,7 @@ class AliAnalysisTaskJetFlowMC : public AliAnalysisTaskSE
         enum detectorType       {kVZEROA, kVZEROC, kVZEROComb};  // detector that was used
         // constructors, destructor
         AliAnalysisTaskJetFlowMC();
-        AliAnalysisTaskJetFlowMC(const char *name);
+        AliAnalysisTaskJetFlowMC(const char *name, Int_t seed = 0);
         virtual ~AliAnalysisTaskJetFlowMC();
         // mirror image of PickTrackMaker
         void    UserCreateOutputObjects();