Insert numerical exchange amplitudes from c3 fits

diff --git a/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.cxx b/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.cxx
index 53dbc878568875a85000a4cf74dfa9137f20d171..9b59a49851aed180349a99796b2d863e24daa32d 100755 (executable)
--- a/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.cxx
+++ b/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.cxx
@@ -579,7 +579,7 @@ AliFourPion &AliFourPion::operator=(const AliFourPion &obj)
   fEvt = obj.fEvt;
   fTempStruct = obj.fTempStruct;
   fRandomNumber = obj.fRandomNumber;
-  fLEGO = fLEGO;
+  fLEGO = obj.fLEGO;
   fMCcase = obj.fMCcase;
   fAODcase = obj.fAODcase;
   fCollisionType = obj.fCollisionType; 
@@ -934,21 +934,21 @@ void AliFourPion::ParInit()
     }
   }
   // Expansion fit parameters: PbPb, pPb, pp
-  Float_t EWParam1Full[3]={0.96827, 0.96827, 0.96827};// PbPb, pPb, pp
-  Float_t EWParam2Full[3]={10.1674, 2.0, 2.0};
-  Float_t EWParam3Full[3]={1.03190e-01, 1.03190e-01, 1.03190e-01};
-  Float_t EWParam4Full[3]={1.84217e-01, 1.84217e-01, 1.84217e-01};
-  Float_t EWParam5Full[3]={0};
+  Float_t EWParam1Full[3]={1.00746, 1.04269, 1.13887};// (lam_3 / 2)^{1/3}
+  Float_t EWParam2Full[3]={1.16820e+01, 2.77982, 2.24686};// Radius
+  Float_t EWParam3Full[3]={2.02382e-01, 2.73547e-01, 3.20273e-01};// kappa3
+  Float_t EWParam4Full[3]={1.43395e-01, 2.41378e-01, 1.49935e-01};// kappa4
+  Float_t EWParam5Full[3]={0};// kappa5
   Float_t EWParam1Point[3][50]={{0}};
   Float_t EWParam2Point[3][50]={{0}};
   Float_t EWParam3Point[3][50]={{0}};
   Float_t EWParam4Point[3][50]={{0}};
   Float_t EWParam5Point[3][50]={{0}};
-  Float_t LGParam1Full[3]={1.3848, 1.3848, 1.3848};
-  Float_t LGParam2Full[3]={24.9985, 4.0, 4.0};
-  Float_t LGParam3Full[3]={1.99598e-01, 1.99598e-01, 1.99598e-01};
-  Float_t LGParam4Full[3]={-7.55834e-02, -7.55834e-02, -7.55834e-02};
-  Float_t LGParam5Full[3]={-8.81810e-03, -8.81810e-03, -8.81810e-03};
+  Float_t LGParam1Full[3]={0.854617, 0.95057, 0.908581};// (lam_3 / 2)^{1/3}
+  Float_t LGParam2Full[3]={1.06967e+01, 2.7007, 1.68574};// Radius
+  Float_t LGParam3Full[3]={-4.44702e-01, -2.24353e-01, -1.92821e-01};// l1
+  Float_t LGParam4Full[3]={-7.98543e-05, 9.16278e-02, 2.57693e-01};// l2
+  Float_t LGParam5Full[3]={0};// l3
   Float_t LGParam1Point[3][50]={{0}};
   Float_t LGParam2Point[3][50]={{0}};
   Float_t LGParam3Point[3][50]={{0}};
@@ -3164,7 +3164,7 @@ void AliFourPion::UserExec(Option_t *)
                    
                    QinvMCGroup4[0] = qinv12MC; QinvMCGroup4[1] = qinv13MC; QinvMCGroup4[2] = qinv14MC;
                    QinvMCGroup4[3] = qinv23MC; QinvMCGroup4[4] = qinv24MC; QinvMCGroup4[5] = qinv34MC;
-                                 
+                         
                  }
                  if(ch1==ch2 && ch1==ch3 && ch1==ch4 && ENsum==6){
                    ((TH2D*)fOutputList->FindObject("DistQinv4pion"))->Fill(1, qinv12); ((TH2D*)fOutputList->FindObject("DistQinv4pion"))->Fill(2, qinv13); 
@@ -3354,6 +3354,7 @@ void AliFourPion::UserExec(Option_t *)
                              T24 = ExchangeAmpFullSource[type-1]->Eval(qinv24) / pow( float(pow(1-G,3) + 3*G*pow(1-G,2)), float(1/3.));
                              T34 = ExchangeAmpFullSource[type-1]->Eval(qinv34) / pow( float(pow(1-G,3) + 3*G*pow(1-G,2)), float(1/3.));
                            }
+                        
                            weightTotal = (1-G*G)*(T12*T12 + T13*T13 + T14*T14 + T23*T23 + T24*T24 + T34*T34);// 2-pion
                            weightTotal += (4*G*pow(1-G,3)+pow(1-G,4))*(pow(T12,2)*pow(T34,2) + pow(T13,2)*pow(T24,2) + pow(T14,2)*pow(T23,2));// 2-pair
                            weightTotal += (6*G*pow(1-G,2) + 2*pow(1-G,3))*(T12*T13*T23);// 3-pion

diff --git a/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.h b/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.h
index bc1089d78965d10897d18791f05e84cd567824ad..733e2a840e47d7adaf24f94c4327e5a736d2550c 100755 (executable)
--- a/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.h
+++ b/PWGCF/FEMTOSCOPY/Chaoticity/AliFourPion.h
@@ -128,8 +128,7 @@ class AliFourPion : public AliAnalysisTaskSE {
   //
   Float_t cubicInterpolate(Float_t[4], Float_t);
   Float_t nCubicInterpolate(Int_t, Float_t*, Float_t[]);
-  //Double_t functionEW(Double_t *x, Double_t *par);
-  //Double_t functionLG(Double_t *x, Double_t *par);
+  
   
   const char* fname;// name of class
   AliAODEvent            *fAOD; //!    // AOD object
@@ -329,22 +328,7 @@ class AliFourPion : public AliAnalysisTaskSE {
   Float_t farrP1[4][4][4];
   Float_t farrP2[4][4][4];
   
-  /* bool LowQPairSwitch_E0E0[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool LowQPairSwitch_E0E1[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool LowQPairSwitch_E0E2[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool LowQPairSwitch_E0E3[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool LowQPairSwitch_E1E2[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool LowQPairSwitch_E1E3[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool LowQPairSwitch_E2E3[kMultLimitPbPb][kMultLimitPbPb];//!
-  //
-  bool NormQPairSwitch_E0E0[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool NormQPairSwitch_E0E1[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool NormQPairSwitch_E0E2[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool NormQPairSwitch_E0E3[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool NormQPairSwitch_E1E2[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool NormQPairSwitch_E1E3[kMultLimitPbPb][kMultLimitPbPb];//!
-  bool NormQPairSwitch_E2E3[kMultLimitPbPb][kMultLimitPbPb];//!
-  */
+ 
   //
   Char_t fDefaultsCharSwitch[kMultLimitPbPb];//!
   TArrayC *fLowQPairSwitch_E0E0[kMultLimitPbPb];//!