variable width for LCC

diff --git a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.cxx b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.cxx
index d970f16fc46bdcd8ce0d701b27cc1a1eaf4ea40b..d9f0dd572f974da31eefa210654a60a2a57b764f 100755 (executable)
--- a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.cxx
+++ b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.cxx
@@ -95,7 +95,8 @@ AliAnalysisTaskToyModel::AliAnalysisTaskToyModel()
   fQuandrangularFlowProtons(0.0), fPentangularFlowProtons(0.0),
   fUseDynamicalCorrelations(kFALSE), fDynamicalCorrelationsPercentage(0.1),
   fUseJets(kFALSE), fPtAssoc(0),
-  fUseLCC(kFALSE) {
+  fUseLCC(kFALSE),
+  fSigmaPt(0.1),fSigmaEta(0.5),fSigmaPhi(0.5){
   // Constructor
 
   //======================================================correction
@@ -808,9 +809,9 @@ void AliAnalysisTaskToyModel::Run(Int_t nEvents) {
        Int_t vCharge_LCC = -vCharge;
        
        // Get Kinematics
-       Double_t vPt_LCC  = gRandom->Gaus(vPt,0.1);
-       Double_t vEta_LCC = gRandom->Gaus(vEta,0.5);
-       Double_t vPhi_LCC = gRandom->Gaus(vPhi,0.5);
+       Double_t vPt_LCC  = gRandom->Gaus(vPt,fSigmaPt);
+       Double_t vEta_LCC = gRandom->Gaus(vEta,fSigmaEta);
+       Double_t vPhi_LCC = gRandom->Gaus(vPhi,fSigmaPhi);
 
        if(fUseDebug) 
          Printf("Generated LCC: Charge = %d, eta = %.2f, phi = %.2f, pt = %.2f",vCharge_LCC,vEta_LCC,vPhi_LCC,vPt_LCC);

diff --git a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.h b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.h
index 35d16fd477e466efb6b2fbc5c07a6283f7bf4fad..913d1045b3bd32f4428f14194b830037d9a60610 100755 (executable)
--- a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.h
+++ b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskToyModel.h
@@ -151,8 +151,15 @@ class AliAnalysisTaskToyModel : public TObject {
   void SetUseJets() {fUseJets = kTRUE;}
 
   // Local charge conservation (LCC)
-  void SetUseLCC() {fUseLCC = kTRUE;}
+  void SetUseLCC() {fUseLCC = kTRUE;}//default values for sigma (pt=0.1,eta=0.5,phi=0.5)
 
+  void SetUseLCC(Double_t sigmaPt, Double_t sigmaEta, Double_t sigmaPhi) {
+    fUseLCC   = kTRUE;
+    fSigmaPt  = sigmaPt; 
+    fSigmaEta = sigmaEta; 
+    fSigmaPhi = sigmaPhi; 
+  }
+  
   //============Toy model: List of setters============//
 
  private:
@@ -287,6 +294,9 @@ class AliAnalysisTaskToyModel : public TObject {
   TF1 *fPtAssoc;//pt of associated
 
   Bool_t fUseLCC;//Usage of Local Charge Conservation
+  Double_t fSigmaPt;//sigma for LCC spread in pT
+  Double_t fSigmaEta;//sigma for LCC spread in Eta
+  Double_t fSigmaPhi;//sigma for LCC spread in Phi
 
   AliAnalysisTaskToyModel(const AliAnalysisTaskToyModel&); // not implemented
   AliAnalysisTaskToyModel& operator=(const AliAnalysisTaskToyModel&); // not implemented