TDPMjet/AliGenDPMjet.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16
  17 // Generator using DPMJET as an external generator
  18 // The main DPMJET options are accessable for the user through this interface.
  19 // Uses the TDPMjet implementation of TGenerator.
  20
  21 #include <TDPMjet.h>
  22 #include "DPMcommon.h"
  23 #include <TRandom.h>
  24 #include <TArrayI.h>
  25 #include <TParticle.h>
  26 #include <TGraph.h>
  27 #include <TDatabasePDG.h>
  28 #include <TParticlePDG.h>
  29 #include <TParticleClassPDG.h>
  30 #include <TPDGCode.h>
  31 #include <TLorentzVector.h>
  32 #include <TClonesArray.h>
  33 #include "AliRunLoader.h"
  34 #include "AliGenDPMjet.h"
  35 #include "AliGenDPMjetEventHeader.h"
  36 #include "AliRun.h"
  37 #include "AliDpmJetRndm.h"
  38 #include "AliIonPDGCodes.h"
  39 #include "AliHeader.h"
  40 #include "AliStack.h"
  41 #include "AliMC.h"
  42 #include "AliLog.h"
  43
  44 ClassImp(AliGenDPMjet)
  45
  46 //______________________________________________________________________________
  47 AliGenDPMjet::AliGenDPMjet()
  48     :AliGenMC(),
  49      fBeamEn(0.),
  50      fMinImpactParam(0.),
  51      fMaxImpactParam(100.),
  52      fICentr(0),
  53      fSelectAll(0),
  54      fFlavor(0),
  55      fTrials(0),
  56      fSpectators(1),
  57      fSpecn(0),
  58      fSpecp(0),
  59      fDPMjet(0),
  60      fNoGammas(0),
  61      fLHC(0),
  62      fPi0Decay(1),
  63      fDecayAll(0),
  64      fGenImpPar(0.),
  65      fProcess(kDpmMb),
  66      fTriggerMultiplicity(0),
  67      fTriggerMultiplicityEta(0),
  68      fTriggerMultiplicityPtMin(0),
  69      fkTuneForDiff(0),
  70      fProcDiff(0),
  71      fFragmentation(kFALSE),
  72      fHeader(AliGenDPMjetEventHeader("DPMJET"))
  73
  74 {
  75 // Constructor
  76     fEnergyCMS = 5500.;
  77     AliDpmJetRndm::SetDpmJetRandom(GetRandom());
  78 }
  79
  80
  81 //______________________________________________________________________________
  82 AliGenDPMjet::AliGenDPMjet(Int_t npart)
  83     :AliGenMC(npart),
  84      fBeamEn(0.),
  85      fMinImpactParam(0.),
  86      fMaxImpactParam(100.),
  87      fICentr(0),
  88      fSelectAll(0),
  89      fFlavor(0),
  90      fTrials(0),
  91      fSpectators(1),
  92      fSpecn(0),
  93      fSpecp(0),
  94      fDPMjet(0),
  95      fNoGammas(0),
  96      fLHC(0),
  97      fPi0Decay(1),
  98      fDecayAll(0),
  99      fGenImpPar(0.),
 100      fProcess(kDpmMb),
 101      fTriggerMultiplicity(0),
 102      fTriggerMultiplicityEta(0),
 103      fTriggerMultiplicityPtMin(0),
 104      fkTuneForDiff(0),
 105      fProcDiff(0),
 106      fFragmentation(kFALSE),
 107      fHeader(AliGenDPMjetEventHeader("DPMJET"))
 108
 109
 110 {
 111 // Default PbPb collisions at 5. 5 TeV
 112 //
 113     fEnergyCMS = 5500.;
 114     fName = "DPMJET";
 115     fTitle= "Particle Generator using DPMJET";
 116     SetTarget();
 117     SetProjectile();
 118     fVertex.Set(3);
 119     AliDpmJetRndm::SetDpmJetRandom(GetRandom());
 120 }
 121
 122 AliGenDPMjet::AliGenDPMjet(const AliGenDPMjet &/*Dpmjet*/)
 123     :AliGenMC(),
 124      fBeamEn(0.),
 125      fMinImpactParam(0.),
 126      fMaxImpactParam(100.),
 127      fICentr(0),
 128      fSelectAll(0),
 129      fFlavor(0),
 130      fTrials(0),
 131      fSpectators(1),
 132      fSpecn(0),
 133      fSpecp(0),
 134      fDPMjet(0),
 135      fNoGammas(0),
 136      fLHC(0),
 137      fPi0Decay(1),
 138      fDecayAll(0),
 139      fGenImpPar(0.),
 140      fProcess(kDpmMb),
 141      fTriggerMultiplicity(0),
 142      fTriggerMultiplicityEta(0),
 143      fTriggerMultiplicityPtMin(0),
 144      fkTuneForDiff(0),
 145      fProcDiff(0),
 146      fFragmentation(kFALSE),
 147      fHeader(0x0)
 148
 149
 150 {
 151     // Dummy copy constructor
 152     fEnergyCMS = 5500.;
 153 }
 154
 155 //______________________________________________________________________________
 156 AliGenDPMjet::~AliGenDPMjet()
 157 {
 158 // Destructor
 159 }
 160 //______________________________________________________________________________
 161 void AliGenDPMjet::Init()
 162 {
 163 // Initialization
 164
 165     if(fEnergyCMS>0. && fBeamEn<0.1) fBeamEn = fEnergyCMS/2;
 166     SetMC(new TDPMjet(fProcess, fAProjectile, fZProjectile, fATarget, fZTarget,
 167                       fBeamEn,fEnergyCMS));
 168
 169     fDPMjet=(TDPMjet*) fMCEvGen;
 170     if (fAProjectile == 1 && TMath::Abs(fZProjectile == 1)) fDPMjet->SetfIdp(1);
 171
 172     // **** Flag to force central production
 173     // fICentr=1. central production forced
 174     // fICentr<0 && fICentr>-100 -> bmin = fMinImpactParam, bmax = fMaxImpactParam
 175     // fICentr<-99 -> fraction of x-sec. = XSFRAC
 176     // fICentr=-1. -> evaporation/fzc suppressed
 177     // fICentr<-1. -> evaporation/fzc allowed
 178     fDPMjet->SetfFCentr(fICentr);
 179
 180     fDPMjet->SetbRange(fMinImpactParam, fMaxImpactParam);
 181     fDPMjet->SetPi0Decay(fPi0Decay);
 182     fDPMjet->SetDecayAll(fDecayAll);
 183     fDPMjet->SetFragmentProd(fFragmentation);
 184
 185 //
 186 //  Initialize DPMjet
 187 //
 188     fDPMjet->Initialize();
 189 }
 190
 191
 192 //______________________________________________________________________________
 193 void AliGenDPMjet::Generate()
 194 {
 195 // Generate one event
 196
 197   Double_t polar[3]    =   {0,0,0};
 198   Double_t origin[3]   =   {0,0,0};
 199   Double_t p[4]        =   {0};
 200   Float_t tof;
 201
 202 //  converts from mm/c to s
 203   const Float_t kconv = 0.001/2.999792458e8;
 204   Int_t nt  = 0;
 205   Int_t jev = 0;
 206   Int_t kf, ks, imo;
 207   kf = 0;
 208   fTrials = 0;
 209   //  Set collision vertex position
 210   if (fVertexSmear == kPerEvent) Vertex();
 211
 212   while(1)
 213   {
 214 //    Generate one event
 215 // --------------------------------------------------------------------------
 216       fSpecn = 0;
 217       fSpecp = 0;
 218 // --------------------------------------------------------------------------
 219       fDPMjet->GenerateEvent();
 220
 221       fTrials++;
 222
 223       fDPMjet->ImportParticles(&fParticles,"All");
 224       if (fLHC) Boost();
 225
 226       // Temporaneo
 227       fGenImpPar = fDPMjet->GetBImpac();
 228
 229       if(TMath::Abs(fXingAngleY) > 1.e-10) BeamCrossAngle();
 230
 231       Int_t np = fParticles.GetEntriesFast();
 232       //
 233       // Multiplicity Trigger
 234       if (fTriggerMultiplicity > 0) {
 235         Int_t multiplicity = 0;
 236         for (Int_t i = 0; i < np; i++) {
 237           TParticle *  iparticle = (TParticle *) fParticles.At(i);
 238
 239           Int_t statusCode = iparticle->GetStatusCode();
 240
 241           // Initial state particle
 242           if (statusCode != 1)
 243             continue;
 244           // eta cut
 245           if (fTriggerMultiplicityEta > 0 && TMath::Abs(iparticle->Eta()) > fTriggerMultiplicityEta)
 246             continue;
 247           // pt cut
 248           if (iparticle->Pt() < fTriggerMultiplicityPtMin)
 249             continue;
 250
 251           TParticlePDG* pdgPart = iparticle->GetPDG();
 252           if (pdgPart && pdgPart->Charge() == 0)
 253             continue;
 254           ++multiplicity;
 255         }
 256         //
 257         //
 258         if (multiplicity < fTriggerMultiplicity) continue;
 259         Printf("Triggered on event with multiplicity of %d >= %d", multiplicity, fTriggerMultiplicity);
 260       }
 261
 262
 263       if(fkTuneForDiff && ( (TMath::Abs(fEnergyCMS - 900) < 1) || (TMath::Abs(fEnergyCMS - 2760) < 1) || (TMath::Abs(fEnergyCMS - 7000) < 1)) ) {
 264         if(!CheckDiffraction() ) continue;
 265       }
 266
 267
 268       Int_t nc = 0;
 269       if (np == 0) continue;
 270
 271       Int_t i;
 272       Int_t* newPos     = new Int_t[np];
 273       Int_t* pSelected  = new Int_t[np];
 274
 275       for (i = 0; i<np; i++) {
 276           newPos[i]    = i;
 277           pSelected[i] = 0;
 278       }
 279
 280 //      First select parent particles
 281
 282       for (i = 0; i<np; i++) {
 283           TParticle *iparticle = (TParticle *) fParticles.At(i);
 284
 285 // Is this a parent particle ?
 286
 287           if (Stable(iparticle)) continue;
 288
 289           Bool_t  selected             =  kTRUE;
 290           Bool_t  hasSelectedDaughters =  kFALSE;
 291
 292           kf = iparticle->GetPdgCode();
 293           if (kf == 92 || kf == 99999) continue;
 294           ks = iparticle->GetStatusCode();
 295 // No initial state partons
 296           if (ks==21) continue;
 297           if (!fSelectAll) selected = KinematicSelection(iparticle, 0) &&
 298                                SelectFlavor(kf);
 299
 300
 301           hasSelectedDaughters = DaughtersSelection(iparticle);
 302
 303
 304 // Put particle on the stack if it is either selected or
 305 // it is the mother of at least one seleted particle
 306
 307           if (selected || hasSelectedDaughters) {
 308               nc++;
 309               pSelected[i] = 1;
 310           } // selected
 311       } // particle loop parents
 312
 313 // Now select the final state particles
 314
 315
 316       for (i=0; i<np; i++) {
 317           TParticle *iparticle = (TParticle *) fParticles.At(i);
 318
 319 // Is this a final state particle ?
 320
 321           if (!Stable(iparticle)) continue;
 322
 323           Bool_t  selected =  kTRUE;
 324           kf = iparticle->GetPdgCode();
 325           ks = iparticle->GetStatusCode();
 326
 327 // --------------------------------------------------------------------------
 328 // Count spectator neutrons and protons (ks == 13, 14)
 329           if(ks == 13 || ks == 14){
 330               if(kf == kNeutron) fSpecn += 1;
 331               if(kf == kProton)  fSpecp += 1;
 332           }
 333 // --------------------------------------------------------------------------
 334
 335           if (!fSelectAll) {
 336               selected = KinematicSelection(iparticle,0)&&SelectFlavor(kf);
 337               if (!fSpectators && selected) selected = (ks == 13 || ks == 14);
 338           }
 339
 340 // Put particle on the stack if selected
 341
 342           if (selected) {
 343               nc++;
 344               pSelected[i] = 1;
 345           } // selected
 346       } // particle loop final state
 347
 348 // Write particles to stack
 349
 350       for (i = 0; i<np; i++) {
 351           TParticle *  iparticle = (TParticle *) fParticles.At(i);
 352           Bool_t  hasMother   = (iparticle->GetFirstMother()>=0);
 353           if (pSelected[i]) {
 354
 355               kf   = iparticle->GetPdgCode();
 356               ks   = iparticle->GetStatusCode();
 357
 358               p[0] = iparticle->Px();
 359               p[1] = iparticle->Py();
 360               p[2] = iparticle->Pz();
 361               p[3] = iparticle->Energy();
 362               origin[0] = fVertex[0]+iparticle->Vx()/10; // [cm]
 363               origin[1] = fVertex[1]+iparticle->Vy()/10; // [cm]
 364               origin[2] = fVertex[2]+iparticle->Vz()/10; // [cm]
 365
 366               tof = fTime + kconv*iparticle->T();
 367
 368               imo = -1;
 369               TParticle* mother = 0;
 370               if (hasMother) {
 371                   imo = iparticle->GetFirstMother();
 372                   mother = (TParticle *) fParticles.At(imo);
 373                   imo = (mother->GetPdgCode() != 92 && mother->GetPdgCode() != 99999) ? newPos[imo] : -1;
 374               } // if has mother
 375
 376
 377
 378               Bool_t tFlag = (fTrackIt && (ks==1 || ks==-1));
 379               //printf(" AliGemDPMJet->PushTrack: kf %d  ks %d  flag %d\n",kf,ks,tFlag);
 380               if(kf>10000 && (ks==-1 || ks==1000 || ks==1001)) kf += 1000000000;
 381               PushTrack(tFlag, imo, kf,
 382                         p[0], p[1], p[2], p[3],
 383                         origin[0], origin[1], origin[2], tof,
 384                         polar[0], polar[1], polar[2],
 385                         kPNoProcess, nt, 1., ks);
 386               KeepTrack(nt);
 387               newPos[i] = nt;
 388           } // if selected
 389       } // particle loop
 390       delete[] newPos;
 391       delete[] pSelected;
 392       if (nc>0) {
 393           jev += nc;
 394           if (jev >= fNpart || fNpart == -1) {
 395               break;
 396           }
 397       }
 398   } // event loop
 399   MakeHeader();
 400   SetHighWaterMark(nt);
 401 }
 402
 403 //______________________________________________________________________________
 404 Bool_t AliGenDPMjet::DaughtersSelection(TParticle* iparticle)
 405 {
 406 //
 407 // Looks recursively if one of the daughters has been selected
 408 //
 409 //    printf("\n Consider daughters %d:",iparticle->GetPdgCode());
 410     Int_t imin = -1;
 411     Int_t imax = -1;
 412     Int_t i;
 413     Bool_t hasDaughters = (iparticle->GetFirstDaughter() >=0);
 414     Bool_t selected = kFALSE;
 415     if (hasDaughters) {
 416         imin = iparticle->GetFirstDaughter();
 417         imax = iparticle->GetLastDaughter();
 418         for (i = imin; i <= imax; i++){
 419             TParticle *  jparticle = (TParticle *) fParticles.At(i);
 420             Int_t ip = jparticle->GetPdgCode();
 421             if (KinematicSelection(jparticle,0)&&SelectFlavor(ip)) {
 422                 selected=kTRUE; break;
 423             }
 424             if (DaughtersSelection(jparticle)) {selected=kTRUE; break; }
 425         }
 426     } else {
 427         return kFALSE;
 428     }
 429     return selected;
 430 }
 431
 432
 433
 434 //______________________________________________________________________________
 435 Bool_t AliGenDPMjet::SelectFlavor(Int_t pid)
 436 {
 437 // Select flavor of particle
 438 // 0: all
 439 // 4: charm and beauty
 440 // 5: beauty
 441     Bool_t res = 0;
 442
 443     if (fFlavor == 0) {
 444         res = kTRUE;
 445     } else {
 446         Int_t ifl = TMath::Abs(pid/100);
 447         if (ifl > 10) ifl/=10;
 448         res = (fFlavor == ifl);
 449     }
 450 //
 451 //  This part if gamma writing is inhibited
 452     if (fNoGammas)
 453         res = res && (pid != kGamma && pid != kPi0);
 454 //
 455     return res;
 456 }
 457
 458 //______________________________________________________________________________
 459 Bool_t AliGenDPMjet::Stable(TParticle*  particle)
 460 {
 461 // Return true for a stable particle
 462 //
 463     int st = particle->GetStatusCode();
 464     if(st == 1 || st == -1) return kTRUE;
 465     else return kFALSE;
 466
 467 }
 468
 469 //______________________________________________________________________________
 470 void AliGenDPMjet::MakeHeader()
 471 {
 472 // Builds the event header, to be called after each event
 473     fHeader.SetNProduced(fDPMjet->GetNumStablePc());
 474     fHeader.SetImpactParameter(fDPMjet->GetBImpac());
 475     fHeader.SetTotalEnergy(fDPMjet->GetTotEnergy());
 476     fHeader.SetParticipants(fDPMjet->GetProjParticipants(),
 477                             fDPMjet->GetTargParticipants());
 478
 479     fHeader.SetCollisions(DTGLCP.ncp, DTGLCP.nct,
 480         fDPMjet->GetProjWounded(),fDPMjet->GetTargWounded());
 481
 482     if(fProcDiff>0)  fHeader.SetProcessType(fProcDiff);
 483     else fHeader.SetProcessType(fDPMjet->GetProcessCode());
 484
 485     // Bookkeeping for kinematic bias
 486     fHeader.SetTrials(fTrials);
 487     // Event Vertex
 488     fHeader.SetPrimaryVertex(fVertex);
 489     fHeader.SetInteractionTime(fTime);
 490     fHeader.SetNDiffractive(POEVT1.nsd1, POEVT1.nsd2, POEVT1.ndd);
 491 //    gAlice->SetGenEventHeader(fHeader);
 492     AddHeader(&fHeader);
 493     fCollisionGeometry = &fHeader;
 494 }
 495
 496 //______________________________________________________________________________
 497 /*void AliGenDPMjet::AddHeader(AliGenEventHeader* fHeader)
 498 {
 499     // Add fHeader to container or runloader
 500     if (fContainer) {
 501         fContainer->AddHeader(fHeader);
 502     } else {
 503         AliRunLoader::Instance()->GetHeader()->SetGenEventHeader(fHeader);
 504     }
 505 }*/
 506
 507
 508 //______________________________________________________________________________
 509 AliGenDPMjet& AliGenDPMjet::operator=(const  AliGenDPMjet& /*rhs*/)
 510 {
 511 // Assignment operator
 512     return *this;
 513 }
 514
 515
 516 //______________________________________________________________________________
 517 void AliGenDPMjet::FinishRun()
 518 {
 519     // Print run statistics
 520     fDPMjet->Dt_Dtuout();
 521 }
 522
 523
 524 //______________________________________________________________________________
 525 Bool_t AliGenDPMjet::CheckDiffraction()
 526 {
 527
 528   //  printf("AAA\n");
 529
 530    Int_t np = fParticles.GetEntriesFast();
 531
 532    Int_t iPart1=-1;
 533    Int_t iPart2=-1;
 534
 535    Double_t y1 = 1e10;
 536    Double_t y2 = -1e10;
 537
 538   const Int_t kNstable=20;
 539   const Int_t pdgStable[20] = {
 540     22,             // Photon
 541     11,             // Electron
 542     12,             // Electron Neutrino
 543     13,             // Muon
 544     14,             // Muon Neutrino
 545     15,             // Tau
 546     16,             // Tau Neutrino
 547     211,            // Pion
 548     321,            // Kaon
 549     311,            // K0
 550     130,            // K0s
 551     310,            // K0l
 552     2212,           // Proton
 553     2112,           // Neutron
 554     3122,           // Lambda_0
 555     3112,           // Sigma Minus
 556     3222,           // Sigma Plus
 557     3312,           // Xsi Minus
 558     3322,           // Xsi0
 559     3334            // Omega
 560   };
 561
 562      for (Int_t i = 0; i < np; i++) {
 563         TParticle *  part = (TParticle *) fParticles.At(i);
 564
 565         Int_t statusCode = part->GetStatusCode();
 566
 567         // Initial state particle
 568         if (statusCode != 1)
 569           continue;
 570
 571         Int_t pdg = TMath::Abs(part->GetPdgCode());
 572         Bool_t isStable = kFALSE;
 573         for (Int_t i1 = 0; i1 < kNstable; i1++) {
 574           if (pdg == pdgStable[i1]) {
 575             isStable = kTRUE;
 576             break;
 577           }
 578         }
 579         if(!isStable)
 580           continue;
 581
 582         Double_t y = part->Y();
 583
 584         if (y < y1)
 585           {
 586             y1 = y;
 587             iPart1 = i;
 588           }
 589         if (y > y2)
 590         {
 591           y2 = y;
 592           iPart2 = i;
 593         }
 594      }
 595
 596      if(iPart1<0 || iPart2<0) return kFALSE;
 597
 598      y1=TMath::Abs(y1);
 599      y2=TMath::Abs(y2);
 600
 601      TParticle *  part1 = (TParticle *) fParticles.At(iPart1);
 602      TParticle *  part2 = (TParticle *) fParticles.At(iPart2);
 603
 604      Int_t pdg1 = part1->GetPdgCode();
 605      Int_t pdg2 = part2->GetPdgCode();
 606
 607
 608      Int_t iPart = -1;
 609      if (pdg1 == 2212 && pdg2 == 2212)
 610        {
 611          if(y1 > y2)
 612            iPart = iPart1;
 613          else if(y1 < y2)
 614            iPart = iPart2;
 615          else {
 616            iPart = iPart1;
 617            if((AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.)>0.5) iPart = iPart2;
 618          }
 619        }
 620      else if (pdg1 == 2212)
 621        iPart = iPart1;
 622      else if (pdg2 == 2212)
 623        iPart = iPart2;
 624
 625
 626
 627
 628
 629      Double_t M=-1.;
 630      if(iPart>0) {
 631        TParticle *  part = (TParticle *) fParticles.At(iPart);
 632        Double_t E= part->Energy();
 633        Double_t P= part->P();
 634        M= TMath::Sqrt((fEnergyCMS-E-P)*(fEnergyCMS-E+P));
 635      }
 636
 637      Double_t Mmin, Mmax, wSD, wDD, wND;
 638
 639      if(!GetWeightsDiffraction(M, Mmin, Mmax, wSD, wDD, wND)) return kFALSE;
 640      if(M>-1 && M<Mmin) return kFALSE;
 641      if(M>Mmax) M=-1;
 642
 643
 644  Int_t procType=fDPMjet->GetProcessCode();//fPythia->GetMSTI(1);
 645  if(procType== 4) return kFALSE;
 646  Int_t proc0=2;
 647  if(procType== 7) proc0=1;
 648  if(procType== 5 || procType== 6) proc0=0;
 649
 650
 651  Int_t proc=2;
 652  if(M>0) proc=0;
 653  else if(proc0==1) proc=1;
 654
 655  if(proc==1 && (AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.) > wDD) return kFALSE;
 656  if(proc==2 && (AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.) > wND) return kFALSE;
 657
 658
 659
 660     if(proc!=0) {
 661       if(proc0!=0) fProcDiff = procType;
 662       else       fProcDiff = 1;
 663       return kTRUE;
 664     }
 665
 666     if(wSD<0)  AliError("wSD<0 ! \n");
 667
 668     if((AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.)> wSD) return kFALSE;
 669
 670     //    printf("iPart = %d\n", iPart);
 671
 672     if(iPart==iPart1) fProcDiff=5;
 673     else if(iPart==iPart2) fProcDiff=6;
 674     else {
 675       printf("EROOR:  iPart!=iPart1 && iPart!=iPart2\n");
 676
 677     }
 678
 679     return kTRUE;
 680 }
 681
 682 // -------------------------------------------------------
 683 void AliGenDPMjet::SetIonPDGCodes()
 684 {
 685    // Defining PDG codes for the ions
 686    AliIonPDGCodes *pdgcodes = new AliIonPDGCodes();
 687    pdgcodes->AddParticlesToPdgDataBase();
 688 }
 689
 690 // -------------------------------------------------------
 691 Bool_t AliGenDPMjet::GetWeightsDiffraction(Double_t M, Double_t &Mmin, Double_t &Mmax,
 692                                                        Double_t &wSD, Double_t &wDD, Double_t &wND)
 693 {
 694
 695   if(TMath::Abs(fEnergyCMS-900)<1 ){
 696 const Int_t nbin=400;
 697 Double_t bin[]={
 698 1.080000, 1.577300, 2.074600, 2.571900, 3.069200, 3.566500,
 699 4.063800, 4.561100, 5.058400, 5.555700, 6.053000, 6.550300,
 700 7.047600, 7.544900, 8.042200, 8.539500, 9.036800, 9.534100,
 701 10.031400, 10.528700, 11.026000, 11.523300, 12.020600, 12.517900,
 702 13.015200, 13.512500, 14.009800, 14.507100, 15.004400, 15.501700,
 703 15.999000, 16.496300, 16.993600, 17.490900, 17.988200, 18.485500,
 704 18.982800, 19.480100, 19.977400, 20.474700, 20.972000, 21.469300,
 705 21.966600, 22.463900, 22.961200, 23.458500, 23.955800, 24.453100,
 706 24.950400, 25.447700, 25.945000, 26.442300, 26.939600, 27.436900,
 707 27.934200, 28.431500, 28.928800, 29.426100, 29.923400, 30.420700,
 708 30.918000, 31.415300, 31.912600, 32.409900, 32.907200, 33.404500,
 709 33.901800, 34.399100, 34.896400, 35.393700, 35.891000, 36.388300,
 710 36.885600, 37.382900, 37.880200, 38.377500, 38.874800, 39.372100,
 711 39.869400, 40.366700, 40.864000, 41.361300, 41.858600, 42.355900,
 712 42.853200, 43.350500, 43.847800, 44.345100, 44.842400, 45.339700,
 713 45.837000, 46.334300, 46.831600, 47.328900, 47.826200, 48.323500,
 714 48.820800, 49.318100, 49.815400, 50.312700, 50.810000, 51.307300,
 715 51.804600, 52.301900, 52.799200, 53.296500, 53.793800, 54.291100,
 716 54.788400, 55.285700, 55.783000, 56.280300, 56.777600, 57.274900,
 717 57.772200, 58.269500, 58.766800, 59.264100, 59.761400, 60.258700,
 718 60.756000, 61.253300, 61.750600, 62.247900, 62.745200, 63.242500,
 719 63.739800, 64.237100, 64.734400, 65.231700, 65.729000, 66.226300,
 720 66.723600, 67.220900, 67.718200, 68.215500, 68.712800, 69.210100,
 721 69.707400, 70.204700, 70.702000, 71.199300, 71.696600, 72.193900,
 722 72.691200, 73.188500, 73.685800, 74.183100, 74.680400, 75.177700,
 723 75.675000, 76.172300, 76.669600, 77.166900, 77.664200, 78.161500,
 724 78.658800, 79.156100, 79.653400, 80.150700, 80.648000, 81.145300,
 725 81.642600, 82.139900, 82.637200, 83.134500, 83.631800, 84.129100,
 726 84.626400, 85.123700, 85.621000, 86.118300, 86.615600, 87.112900,
 727 87.610200, 88.107500, 88.604800, 89.102100, 89.599400, 90.096700,
 728 90.594000, 91.091300, 91.588600, 92.085900, 92.583200, 93.080500,
 729 93.577800, 94.075100, 94.572400, 95.069700, 95.567000, 96.064300,
 730 96.561600, 97.058900, 97.556200, 98.053500, 98.550800, 99.048100,
 731 99.545400, 100.042700, 100.540000, 101.037300, 101.534600, 102.031900,
 732 102.529200, 103.026500, 103.523800, 104.021100, 104.518400, 105.015700,
 733 105.513000, 106.010300, 106.507600, 107.004900, 107.502200, 107.999500,
 734 108.496800, 108.994100, 109.491400, 109.988700, 110.486000, 110.983300,
 735 111.480600, 111.977900, 112.475200, 112.972500, 113.469800, 113.967100,
 736 114.464400, 114.961700, 115.459000, 115.956300, 116.453600, 116.950900,
 737 117.448200, 117.945500, 118.442800, 118.940100, 119.437400, 119.934700,
 738 120.432000, 120.929300, 121.426600, 121.923900, 122.421200, 122.918500,
 739 123.415800, 123.913100, 124.410400, 124.907700, 125.405000, 125.902300,
 740 126.399600, 126.896900, 127.394200, 127.891500, 128.388800, 128.886100,
 741 129.383400, 129.880700, 130.378000, 130.875300, 131.372600, 131.869900,
 742 132.367200, 132.864500, 133.361800, 133.859100, 134.356400, 134.853700,
 743 135.351000, 135.848300, 136.345600, 136.842900, 137.340200, 137.837500,
 744 138.334800, 138.832100, 139.329400, 139.826700, 140.324000, 140.821300,
 745 141.318600, 141.815900, 142.313200, 142.810500, 143.307800, 143.805100,
 746 144.302400, 144.799700, 145.297000, 145.794300, 146.291600, 146.788900,
 747 147.286200, 147.783500, 148.280800, 148.778100, 149.275400, 149.772700,
 748 150.270000, 150.767300, 151.264600, 151.761900, 152.259200, 152.756500,
 749 153.253800, 153.751100, 154.248400, 154.745700, 155.243000, 155.740300,
 750 156.237600, 156.734900, 157.232200, 157.729500, 158.226800, 158.724100,
 751 159.221400, 159.718700, 160.216000, 160.713300, 161.210600, 161.707900,
 752 162.205200, 162.702500, 163.199800, 163.697100, 164.194400, 164.691700,
 753 165.189000, 165.686300, 166.183600, 166.680900, 167.178200, 167.675500,
 754 168.172800, 168.670100, 169.167400, 169.664700, 170.162000, 170.659300,
 755 171.156600, 171.653900, 172.151200, 172.648500, 173.145800, 173.643100,
 756 174.140400, 174.637700, 175.135000, 175.632300, 176.129600, 176.626900,
 757 177.124200, 177.621500, 178.118800, 178.616100, 179.113400, 179.610700,
 758 180.108000, 180.605300, 181.102600, 181.599900, 182.097200, 182.594500,
 759 183.091800, 183.589100, 184.086400, 184.583700, 185.081000, 185.578300,
 760 186.075600, 186.572900, 187.070200, 187.567500, 188.064800, 188.562100,
 761 189.059400, 189.556700, 190.054000, 190.551300, 191.048600, 191.545900,
 762 192.043200, 192.540500, 193.037800, 193.535100, 194.032400, 194.529700,
 763 195.027000, 195.524300, 196.021600, 196.518900, 197.016200, 197.513500,
 764 198.010800, 198.508100, 199.005400, 199.502700, 200.000000};
 765 Double_t w[]={
 766 1.000000, 0.389079, 0.326612, 0.356363, 0.310726, 0.264037,
 767 0.253810, 0.224655, 0.207990, 0.198149, 0.186803, 0.180209,
 768 0.178928, 0.161772, 0.162314, 0.158409, 0.148379, 0.143557,
 769 0.140357, 0.136272, 0.136957, 0.136606, 0.129376, 0.127671,
 770 0.128406, 0.132239, 0.119311, 0.130227, 0.130233, 0.123312,
 771 0.115257, 0.120105, 0.119733, 0.117116, 0.110096, 0.117509,
 772 0.109149, 0.114372, 0.100322, 0.106227, 0.108696, 0.110352,
 773 0.106748, 0.101475, 0.101837, 0.098585, 0.094433, 0.100148,
 774 0.096505, 0.100155, 0.103526, 0.098161, 0.093647, 0.100121,
 775 0.097426, 0.093414, 0.090241, 0.097228, 0.098098, 0.098340,
 776 0.096779, 0.099265, 0.095101, 0.090342, 0.097469, 0.090663,
 777 0.093422, 0.093103, 0.089648, 0.096923, 0.088079, 0.089360,
 778 0.097772, 0.092932, 0.093554, 0.085784, 0.096395, 0.096304,
 779 0.090183, 0.089255, 0.090265, 0.086262, 0.087044, 0.088965,
 780 0.085825, 0.095073, 0.085933, 0.091005, 0.091707, 0.092428,
 781 0.091689, 0.091224, 0.095256, 0.083407, 0.087983, 0.090320,
 782 0.081580, 0.088077, 0.083478, 0.091309, 0.083734, 0.089906,
 783 0.100155, 0.092728, 0.086542, 0.078091, 0.085261, 0.094302,
 784 0.078063, 0.090070, 0.086566, 0.095020, 0.082124, 0.084791,
 785 0.090624, 0.090236, 0.086344, 0.085706, 0.085913, 0.083107,
 786 0.092557, 0.081144, 0.091254, 0.096139, 0.079458, 0.088124,
 787 0.088777, 0.076652, 0.085168, 0.082326, 0.080435, 0.085022,
 788 0.085693, 0.093957, 0.088134, 0.082347, 0.082139, 0.090980,
 789 0.078397, 0.093595, 0.091562, 0.081789, 0.085026, 0.080868,
 790 0.083455, 0.084590, 0.084124, 0.081486, 0.087242, 0.080928,
 791 0.096392, 0.078324, 0.092093, 0.074267, 0.082925, 0.090717,
 792 0.081467, 0.080435, 0.078825, 0.075843, 0.088940, 0.081451,
 793 0.096047, 0.088102, 0.091710, 0.095208, 0.086160, 0.091305,
 794 0.085406, 0.094432, 0.078227, 0.099870, 0.094140, 0.080130,
 795 0.090707, 0.083268, 0.082222, 0.088767, 0.084477, 0.088069,
 796 0.089382, 0.086164, 0.089123, 0.092799, 0.076710, 0.090727,
 797 0.077097, 0.099905, 0.089733, 0.088101, 0.093705, 0.089215,
 798 0.085110, 0.086032, 0.083719, 0.079693, 0.088116, 0.090519,
 799 0.091150, 0.090855, 0.090547, 0.077773, 0.081914, 0.080864,
 800 0.082935, 0.092952, 0.079390, 0.080255, 0.091123, 0.091331,
 801 0.095160, 0.089343, 0.083353, 0.087445, 0.094036, 0.084719,
 802 0.099665, 0.085104, 0.100912, 0.085958, 0.083972, 0.098284,
 803 0.078318, 0.082042, 0.088007, 0.085469, 0.089984, 0.081181,
 804 0.098850, 0.086409, 0.089070, 0.081579, 0.089622, 0.089396,
 805 0.093922, 0.089472, 0.090806, 0.075034, 0.090346, 0.083871,
 806 0.086931, 0.089207, 0.094425, 0.087830, 0.079537, 0.097316,
 807 0.096513, 0.092264, 0.082211, 0.083841, 0.081861, 0.092591,
 808 0.077785, 0.079646, 0.093721, 0.090735, 0.086910, 0.075837,
 809 0.090729, 0.092800, 0.098704, 0.092441, 0.086404, 0.089344,
 810 0.083650, 0.082569, 0.085753, 0.091430, 0.085460, 0.095210,
 811 0.083226, 0.094882, 0.084856, 0.093350, 0.084579, 0.096744,
 812 0.099211, 0.085142, 0.085499, 0.083395, 0.088352, 0.110418,
 813 0.093788, 0.096140, 0.082758, 0.097145, 0.089170, 0.090720,
 814 0.084708, 0.095927, 0.090330, 0.074239, 0.089532, 0.090514,
 815 0.089823, 0.089709, 0.101840, 0.082676, 0.087157, 0.079221,
 816 0.096460, 0.108192, 0.088904, 0.084510, 0.096624, 0.099242,
 817 0.094470, 0.092473, 0.091745, 0.090439, 0.095316, 0.087963,
 818 0.096400, 0.089479, 0.094880, 0.085964, 0.092775, 0.083200,
 819 0.095133, 0.090079, 0.088828, 0.087600, 0.094123, 0.089135,
 820 0.082617, 0.085109, 0.114091, 0.096331, 0.108465, 0.080318,
 821 0.106576, 0.089671, 0.092023, 0.090722, 0.086603, 0.091788,
 822 0.098375, 0.082712, 0.102681, 0.098869, 0.089051, 0.109972,
 823 0.080440, 0.097860, 0.093074, 0.097028, 0.107826, 0.117152,
 824 0.090968, 0.096790, 0.096725, 0.094641, 0.089535, 0.092561,
 825 0.095828, 0.084320, 0.089942, 0.105476, 0.087495, 0.089805,
 826 0.092238, 0.094141, 0.104390, 0.082958, 0.097449, 0.099594,
 827 0.092640, 0.097332, 0.093223, 0.100183, 0.092511, 0.087035,
 828 0.088741, 0.097856, 0.116682, 0.091732, 0.097753, 0.094283,
 829 0.094235, 0.082016, 0.098370, 0.085676, 0.104529, 0.087319,
 830 0.090289, 0.105314, 0.103634, 0.101261, 0.092764, 0.098217,
 831 0.098939, 0.096071, 0.096071, 0.094027, 0.092713, 0.089542,
 832 0.112293, 0.112293, 0.089531, 0.087752};
 833 wDD = 0.188264;
 834 wND = 0.102252;
 835 wSD = -1;
 836
 837     Mmin = bin[0];
 838     Mmax = bin[nbin];
 839     if(M<Mmin || M>Mmax) return kTRUE;
 840
 841     Int_t ibin=nbin-1;
 842     for(Int_t i=1; i<=nbin; i++)
 843       if(M<=bin[i]) {
 844         ibin=i-1;
 845         //      printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);
 846         break;
 847       }
 848     wSD=w[ibin];
 849     return kTRUE;
 850
 851   }
 852   else if(TMath::Abs(fEnergyCMS-2760)<1 ){
 853
 854 const Int_t nbin=400;
 855 Double_t bin[]={
 856 1.080000, 1.577300, 2.074600, 2.571900, 3.069200, 3.566500,
 857 4.063800, 4.561100, 5.058400, 5.555700, 6.053000, 6.550300,
 858 7.047600, 7.544900, 8.042200, 8.539500, 9.036800, 9.534100,
 859 10.031400, 10.528700, 11.026000, 11.523300, 12.020600, 12.517900,
 860 13.015200, 13.512500, 14.009800, 14.507100, 15.004400, 15.501700,
 861 15.999000, 16.496300, 16.993600, 17.490900, 17.988200, 18.485500,
 862 18.982800, 19.480100, 19.977400, 20.474700, 20.972000, 21.469300,
 863 21.966600, 22.463900, 22.961200, 23.458500, 23.955800, 24.453100,
 864 24.950400, 25.447700, 25.945000, 26.442300, 26.939600, 27.436900,
 865 27.934200, 28.431500, 28.928800, 29.426100, 29.923400, 30.420700,
 866 30.918000, 31.415300, 31.912600, 32.409900, 32.907200, 33.404500,
 867 33.901800, 34.399100, 34.896400, 35.393700, 35.891000, 36.388300,
 868 36.885600, 37.382900, 37.880200, 38.377500, 38.874800, 39.372100,
 869 39.869400, 40.366700, 40.864000, 41.361300, 41.858600, 42.355900,
 870 42.853200, 43.350500, 43.847800, 44.345100, 44.842400, 45.339700,
 871 45.837000, 46.334300, 46.831600, 47.328900, 47.826200, 48.323500,
 872 48.820800, 49.318100, 49.815400, 50.312700, 50.810000, 51.307300,
 873 51.804600, 52.301900, 52.799200, 53.296500, 53.793800, 54.291100,
 874 54.788400, 55.285700, 55.783000, 56.280300, 56.777600, 57.274900,
 875 57.772200, 58.269500, 58.766800, 59.264100, 59.761400, 60.258700,
 876 60.756000, 61.253300, 61.750600, 62.247900, 62.745200, 63.242500,
 877 63.739800, 64.237100, 64.734400, 65.231700, 65.729000, 66.226300,
 878 66.723600, 67.220900, 67.718200, 68.215500, 68.712800, 69.210100,
 879 69.707400, 70.204700, 70.702000, 71.199300, 71.696600, 72.193900,
 880 72.691200, 73.188500, 73.685800, 74.183100, 74.680400, 75.177700,
 881 75.675000, 76.172300, 76.669600, 77.166900, 77.664200, 78.161500,
 882 78.658800, 79.156100, 79.653400, 80.150700, 80.648000, 81.145300,
 883 81.642600, 82.139900, 82.637200, 83.134500, 83.631800, 84.129100,
 884 84.626400, 85.123700, 85.621000, 86.118300, 86.615600, 87.112900,
 885 87.610200, 88.107500, 88.604800, 89.102100, 89.599400, 90.096700,
 886 90.594000, 91.091300, 91.588600, 92.085900, 92.583200, 93.080500,
 887 93.577800, 94.075100, 94.572400, 95.069700, 95.567000, 96.064300,
 888 96.561600, 97.058900, 97.556200, 98.053500, 98.550800, 99.048100,
 889 99.545400, 100.042700, 100.540000, 101.037300, 101.534600, 102.031900,
 890 102.529200, 103.026500, 103.523800, 104.021100, 104.518400, 105.015700,
 891 105.513000, 106.010300, 106.507600, 107.004900, 107.502200, 107.999500,
 892 108.496800, 108.994100, 109.491400, 109.988700, 110.486000, 110.983300,
 893 111.480600, 111.977900, 112.475200, 112.972500, 113.469800, 113.967100,
 894 114.464400, 114.961700, 115.459000, 115.956300, 116.453600, 116.950900,
 895 117.448200, 117.945500, 118.442800, 118.940100, 119.437400, 119.934700,
 896 120.432000, 120.929300, 121.426600, 121.923900, 122.421200, 122.918500,
 897 123.415800, 123.913100, 124.410400, 124.907700, 125.405000, 125.902300,
 898 126.399600, 126.896900, 127.394200, 127.891500, 128.388800, 128.886100,
 899 129.383400, 129.880700, 130.378000, 130.875300, 131.372600, 131.869900,
 900 132.367200, 132.864500, 133.361800, 133.859100, 134.356400, 134.853700,
 901 135.351000, 135.848300, 136.345600, 136.842900, 137.340200, 137.837500,
 902 138.334800, 138.832100, 139.329400, 139.826700, 140.324000, 140.821300,
 903 141.318600, 141.815900, 142.313200, 142.810500, 143.307800, 143.805100,
 904 144.302400, 144.799700, 145.297000, 145.794300, 146.291600, 146.788900,
 905 147.286200, 147.783500, 148.280800, 148.778100, 149.275400, 149.772700,
 906 150.270000, 150.767300, 151.264600, 151.761900, 152.259200, 152.756500,
 907 153.253800, 153.751100, 154.248400, 154.745700, 155.243000, 155.740300,
 908 156.237600, 156.734900, 157.232200, 157.729500, 158.226800, 158.724100,
 909 159.221400, 159.718700, 160.216000, 160.713300, 161.210600, 161.707900,
 910 162.205200, 162.702500, 163.199800, 163.697100, 164.194400, 164.691700,
 911 165.189000, 165.686300, 166.183600, 166.680900, 167.178200, 167.675500,
 912 168.172800, 168.670100, 169.167400, 169.664700, 170.162000, 170.659300,
 913 171.156600, 171.653900, 172.151200, 172.648500, 173.145800, 173.643100,
 914 174.140400, 174.637700, 175.135000, 175.632300, 176.129600, 176.626900,
 915 177.124200, 177.621500, 178.118800, 178.616100, 179.113400, 179.610700,
 916 180.108000, 180.605300, 181.102600, 181.599900, 182.097200, 182.594500,
 917 183.091800, 183.589100, 184.086400, 184.583700, 185.081000, 185.578300,
 918 186.075600, 186.572900, 187.070200, 187.567500, 188.064800, 188.562100,
 919 189.059400, 189.556700, 190.054000, 190.551300, 191.048600, 191.545900,
 920 192.043200, 192.540500, 193.037800, 193.535100, 194.032400, 194.529700,
 921 195.027000, 195.524300, 196.021600, 196.518900, 197.016200, 197.513500,
 922 198.010800, 198.508100, 199.005400, 199.502700, 200.000000};
 923 Double_t w[]={
 924 1.000000, 0.402402, 0.347976, 0.386866, 0.304413, 0.275746,
 925 0.256941, 0.250439, 0.228486, 0.219225, 0.203368, 0.196124,
 926 0.180551, 0.169230, 0.159021, 0.157136, 0.154743, 0.160522,
 927 0.153993, 0.143448, 0.148062, 0.136337, 0.130853, 0.127798,
 928 0.137576, 0.137530, 0.124175, 0.120225, 0.120774, 0.123584,
 929 0.121883, 0.114522, 0.118338, 0.119587, 0.106017, 0.125802,
 930 0.106256, 0.101585, 0.097737, 0.100817, 0.102279, 0.110935,
 931 0.101429, 0.107629, 0.105802, 0.095688, 0.096617, 0.100985,
 932 0.101521, 0.100450, 0.111836, 0.088828, 0.102264, 0.101012,
 933 0.097673, 0.097634, 0.091720, 0.096167, 0.084537, 0.104443,
 934 0.091966, 0.097204, 0.089456, 0.088273, 0.098220, 0.095151,
 935 0.088201, 0.085570, 0.088431, 0.077625, 0.088314, 0.084582,
 936 0.091545, 0.083774, 0.086183, 0.089195, 0.086255, 0.077167,
 937 0.096118, 0.089881, 0.086206, 0.086968, 0.079827, 0.082002,
 938 0.083081, 0.076587, 0.081820, 0.092382, 0.082964, 0.080901,
 939 0.093512, 0.082907, 0.086691, 0.093639, 0.089010, 0.082857,
 940 0.090795, 0.082403, 0.078602, 0.087284, 0.079638, 0.082905,
 941 0.077189, 0.084329, 0.075900, 0.082559, 0.084210, 0.070053,
 942 0.068453, 0.083369, 0.077659, 0.073286, 0.075396, 0.073946,
 943 0.083293, 0.076235, 0.082184, 0.080104, 0.073755, 0.077369,
 944 0.068281, 0.084593, 0.074923, 0.069467, 0.072094, 0.077973,
 945 0.067034, 0.080314, 0.087441, 0.071728, 0.072302, 0.076783,
 946 0.079931, 0.075557, 0.076318, 0.079029, 0.075126, 0.076859,
 947 0.080253, 0.074344, 0.083387, 0.087553, 0.084437, 0.076322,
 948 0.077748, 0.087559, 0.075649, 0.076615, 0.077098, 0.074559,
 949 0.075540, 0.085883, 0.079269, 0.079866, 0.068922, 0.069770,
 950 0.068024, 0.077783, 0.086620, 0.086769, 0.068433, 0.078059,
 951 0.084463, 0.072790, 0.073889, 0.079090, 0.081759, 0.069576,
 952 0.069160, 0.071695, 0.078569, 0.075727, 0.088055, 0.082395,
 953 0.079915, 0.075150, 0.073580, 0.069968, 0.071141, 0.065823,
 954 0.077270, 0.076229, 0.071735, 0.079271, 0.078209, 0.069503,
 955 0.064732, 0.076312, 0.088579, 0.087271, 0.080566, 0.073527,
 956 0.071376, 0.068785, 0.071812, 0.072680, 0.067252, 0.059543,
 957 0.072818, 0.071856, 0.083427, 0.080664, 0.077315, 0.066060,
 958 0.078772, 0.065442, 0.067255, 0.077676, 0.068306, 0.071578,
 959 0.059728, 0.073526, 0.078932, 0.070316, 0.067620, 0.067323,
 960 0.078316, 0.080366, 0.091194, 0.066739, 0.076238, 0.070382,
 961 0.077118, 0.073822, 0.072830, 0.062947, 0.077375, 0.076284,
 962 0.070951, 0.086127, 0.071074, 0.076621, 0.070502, 0.070897,
 963 0.058764, 0.065322, 0.073620, 0.062051, 0.074548, 0.083770,
 964 0.072502, 0.066823, 0.074192, 0.081570, 0.067340, 0.081896,
 965 0.080680, 0.072835, 0.069675, 0.078881, 0.072826, 0.069616,
 966 0.082744, 0.069138, 0.073333, 0.080881, 0.069091, 0.064266,
 967 0.062849, 0.072846, 0.067462, 0.067229, 0.071319, 0.071099,
 968 0.069381, 0.073799, 0.076106, 0.075006, 0.073074, 0.092589,
 969 0.074277, 0.061324, 0.069788, 0.070310, 0.063160, 0.073077,
 970 0.075447, 0.081013, 0.074102, 0.076565, 0.060192, 0.058326,
 971 0.071508, 0.072974, 0.065098, 0.060432, 0.062077, 0.075644,
 972 0.075417, 0.066947, 0.066744, 0.065111, 0.075468, 0.078185,
 973 0.073216, 0.066502, 0.079582, 0.065405, 0.069828, 0.072165,
 974 0.071112, 0.064656, 0.080396, 0.066505, 0.063419, 0.066128,
 975 0.080616, 0.071463, 0.064867, 0.074782, 0.063103, 0.068179,
 976 0.064911, 0.073029, 0.075746, 0.062264, 0.099126, 0.067890,
 977 0.074940, 0.070090, 0.073583, 0.056904, 0.063745, 0.063780,
 978 0.076876, 0.066519, 0.067164, 0.068668, 0.075065, 0.070983,
 979 0.075751, 0.066317, 0.077520, 0.073139, 0.073961, 0.085492,
 980 0.074629, 0.059555, 0.081026, 0.059274, 0.055448, 0.084269,
 981 0.066222, 0.078874, 0.062021, 0.075156, 0.090478, 0.066265,
 982 0.067845, 0.076653, 0.086671, 0.082362, 0.067246, 0.075975,
 983 0.072594, 0.068570, 0.095034, 0.065644, 0.070945, 0.066190,
 984 0.075906, 0.055154, 0.075575, 0.058078, 0.071004, 0.070850,
 985 0.064374, 0.058251, 0.069430, 0.072293, 0.065484, 0.084197,
 986 0.090119, 0.091619, 0.067620, 0.079462, 0.063025, 0.068128,
 987 0.056927, 0.076351, 0.073869, 0.061597, 0.083741, 0.063762,
 988 0.064489, 0.074269, 0.068832, 0.058648, 0.069536, 0.074824,
 989 0.081009, 0.073389, 0.076792, 0.084855, 0.075280, 0.061473,
 990 0.063840, 0.062891, 0.071328, 0.065250};
 991 wDD = 0.077845;
 992 wND = 0.066355;
 993 wSD = -1;
 994
 995     Mmin = bin[0];
 996     Mmax = bin[nbin];
 997     if(M<Mmin || M>Mmax) return kTRUE;
 998
 999     Int_t ibin=nbin-1;
1000     for(Int_t i=1; i<=nbin; i++)
1001       if(M<=bin[i]) {
1002         ibin=i-1;
1003         //      printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);
1004         break;
1005       }
1006     wSD=w[ibin];
1007     return kTRUE;
1008
1009   }
1010   else if(TMath::Abs(fEnergyCMS-7000)<1 ){
1011
1012 const Int_t nbin=400;
1013 Double_t bin[]={
1014 1.080000, 1.577300, 2.074600, 2.571900, 3.069200, 3.566500,
1015 4.063800, 4.561100, 5.058400, 5.555700, 6.053000, 6.550300,
1016 7.047600, 7.544900, 8.042200, 8.539500, 9.036800, 9.534100,
1017 10.031400, 10.528700, 11.026000, 11.523300, 12.020600, 12.517900,
1018 13.015200, 13.512500, 14.009800, 14.507100, 15.004400, 15.501700,
1019 15.999000, 16.496300, 16.993600, 17.490900, 17.988200, 18.485500,
1020 18.982800, 19.480100, 19.977400, 20.474700, 20.972000, 21.469300,
1021 21.966600, 22.463900, 22.961200, 23.458500, 23.955800, 24.453100,
1022 24.950400, 25.447700, 25.945000, 26.442300, 26.939600, 27.436900,
1023 27.934200, 28.431500, 28.928800, 29.426100, 29.923400, 30.420700,
1024 30.918000, 31.415300, 31.912600, 32.409900, 32.907200, 33.404500,
1025 33.901800, 34.399100, 34.896400, 35.393700, 35.891000, 36.388300,
1026 36.885600, 37.382900, 37.880200, 38.377500, 38.874800, 39.372100,
1027 39.869400, 40.366700, 40.864000, 41.361300, 41.858600, 42.355900,
1028 42.853200, 43.350500, 43.847800, 44.345100, 44.842400, 45.339700,
1029 45.837000, 46.334300, 46.831600, 47.328900, 47.826200, 48.323500,
1030 48.820800, 49.318100, 49.815400, 50.312700, 50.810000, 51.307300,
1031 51.804600, 52.301900, 52.799200, 53.296500, 53.793800, 54.291100,
1032 54.788400, 55.285700, 55.783000, 56.280300, 56.777600, 57.274900,
1033 57.772200, 58.269500, 58.766800, 59.264100, 59.761400, 60.258700,
1034 60.756000, 61.253300, 61.750600, 62.247900, 62.745200, 63.242500,
1035 63.739800, 64.237100, 64.734400, 65.231700, 65.729000, 66.226300,
1036 66.723600, 67.220900, 67.718200, 68.215500, 68.712800, 69.210100,
1037 69.707400, 70.204700, 70.702000, 71.199300, 71.696600, 72.193900,
1038 72.691200, 73.188500, 73.685800, 74.183100, 74.680400, 75.177700,
1039 75.675000, 76.172300, 76.669600, 77.166900, 77.664200, 78.161500,
1040 78.658800, 79.156100, 79.653400, 80.150700, 80.648000, 81.145300,
1041 81.642600, 82.139900, 82.637200, 83.134500, 83.631800, 84.129100,
1042 84.626400, 85.123700, 85.621000, 86.118300, 86.615600, 87.112900,
1043 87.610200, 88.107500, 88.604800, 89.102100, 89.599400, 90.096700,
1044 90.594000, 91.091300, 91.588600, 92.085900, 92.583200, 93.080500,
1045 93.577800, 94.075100, 94.572400, 95.069700, 95.567000, 96.064300,
1046 96.561600, 97.058900, 97.556200, 98.053500, 98.550800, 99.048100,
1047 99.545400, 100.042700, 100.540000, 101.037300, 101.534600, 102.031900,
1048 102.529200, 103.026500, 103.523800, 104.021100, 104.518400, 105.015700,
1049 105.513000, 106.010300, 106.507600, 107.004900, 107.502200, 107.999500,
1050 108.496800, 108.994100, 109.491400, 109.988700, 110.486000, 110.983300,
1051 111.480600, 111.977900, 112.475200, 112.972500, 113.469800, 113.967100,
1052 114.464400, 114.961700, 115.459000, 115.956300, 116.453600, 116.950900,
1053 117.448200, 117.945500, 118.442800, 118.940100, 119.437400, 119.934700,
1054 120.432000, 120.929300, 121.426600, 121.923900, 122.421200, 122.918500,
1055 123.415800, 123.913100, 124.410400, 124.907700, 125.405000, 125.902300,
1056 126.399600, 126.896900, 127.394200, 127.891500, 128.388800, 128.886100,
1057 129.383400, 129.880700, 130.378000, 130.875300, 131.372600, 131.869900,
1058 132.367200, 132.864500, 133.361800, 133.859100, 134.356400, 134.853700,
1059 135.351000, 135.848300, 136.345600, 136.842900, 137.340200, 137.837500,
1060 138.334800, 138.832100, 139.329400, 139.826700, 140.324000, 140.821300,
1061 141.318600, 141.815900, 142.313200, 142.810500, 143.307800, 143.805100,
1062 144.302400, 144.799700, 145.297000, 145.794300, 146.291600, 146.788900,
1063 147.286200, 147.783500, 148.280800, 148.778100, 149.275400, 149.772700,
1064 150.270000, 150.767300, 151.264600, 151.761900, 152.259200, 152.756500,
1065 153.253800, 153.751100, 154.248400, 154.745700, 155.243000, 155.740300,
1066 156.237600, 156.734900, 157.232200, 157.729500, 158.226800, 158.724100,
1067 159.221400, 159.718700, 160.216000, 160.713300, 161.210600, 161.707900,
1068 162.205200, 162.702500, 163.199800, 163.697100, 164.194400, 164.691700,
1069 165.189000, 165.686300, 166.183600, 166.680900, 167.178200, 167.675500,
1070 168.172800, 168.670100, 169.167400, 169.664700, 170.162000, 170.659300,
1071 171.156600, 171.653900, 172.151200, 172.648500, 173.145800, 173.643100,
1072 174.140400, 174.637700, 175.135000, 175.632300, 176.129600, 176.626900,
1073 177.124200, 177.621500, 178.118800, 178.616100, 179.113400, 179.610700,
1074 180.108000, 180.605300, 181.102600, 181.599900, 182.097200, 182.594500,
1075 183.091800, 183.589100, 184.086400, 184.583700, 185.081000, 185.578300,
1076 186.075600, 186.572900, 187.070200, 187.567500, 188.064800, 188.562100,
1077 189.059400, 189.556700, 190.054000, 190.551300, 191.048600, 191.545900,
1078 192.043200, 192.540500, 193.037800, 193.535100, 194.032400, 194.529700,
1079 195.027000, 195.524300, 196.021600, 196.518900, 197.016200, 197.513500,
1080 198.010800, 198.508100, 199.005400, 199.502700, 200.000000};
1081 Double_t w[]={
1082 1.000000, 0.526293, 0.446686, 0.437789, 0.366854, 0.333320,
1083 0.291931, 0.266464, 0.253870, 0.248706, 0.232788, 0.220736,
1084 0.209886, 0.202741, 0.188617, 0.182767, 0.178748, 0.169039,
1085 0.175911, 0.169098, 0.171256, 0.146728, 0.144543, 0.159470,
1086 0.153171, 0.151883, 0.144693, 0.136307, 0.140226, 0.135388,
1087 0.141317, 0.151121, 0.131209, 0.144039, 0.124688, 0.128020,
1088 0.119122, 0.121868, 0.122332, 0.119561, 0.115744, 0.102504,
1089 0.114726, 0.109518, 0.115418, 0.114860, 0.110026, 0.107693,
1090 0.103005, 0.115985, 0.108629, 0.105937, 0.101056, 0.101228,
1091 0.113305, 0.110302, 0.104696, 0.107447, 0.099095, 0.107378,
1092 0.103090, 0.111384, 0.090821, 0.109806, 0.093639, 0.096655,
1093 0.110416, 0.104446, 0.098530, 0.095105, 0.106641, 0.101804,
1094 0.091798, 0.094775, 0.104130, 0.088436, 0.119692, 0.099160,
1095 0.086163, 0.089250, 0.101898, 0.091786, 0.087346, 0.095312,
1096 0.107186, 0.085671, 0.093283, 0.091992, 0.114654, 0.086172,
1097 0.084559, 0.097862, 0.079417, 0.094737, 0.089678, 0.105265,
1098 0.099092, 0.087901, 0.086828, 0.084756, 0.077769, 0.084254,
1099 0.083262, 0.092935, 0.088858, 0.093377, 0.083569, 0.084771,
1100 0.084957, 0.084033, 0.095007, 0.076690, 0.087726, 0.097913,
1101 0.084119, 0.076361, 0.084172, 0.091068, 0.094564, 0.091406,
1102 0.088405, 0.090355, 0.085388, 0.084586, 0.085752, 0.084298,
1103 0.091294, 0.081491, 0.075024, 0.082543, 0.088048, 0.074621,
1104 0.078514, 0.077249, 0.078378, 0.092945, 0.075859, 0.077602,
1105 0.074592, 0.078181, 0.081527, 0.080865, 0.078854, 0.078237,
1106 0.085506, 0.089823, 0.072738, 0.078430, 0.077162, 0.081495,
1107 0.090878, 0.099417, 0.100966, 0.087960, 0.085556, 0.081661,
1108 0.078066, 0.089059, 0.073564, 0.092704, 0.073148, 0.098057,
1109 0.067906, 0.079452, 0.090567, 0.082454, 0.077851, 0.079694,
1110 0.091272, 0.078628, 0.096906, 0.098779, 0.088906, 0.070174,
1111 0.083822, 0.084241, 0.093237, 0.071062, 0.075771, 0.096405,
1112 0.098441, 0.086007, 0.069599, 0.078400, 0.083481, 0.081054,
1113 0.085552, 0.069582, 0.071336, 0.073207, 0.087913, 0.080125,
1114 0.075189, 0.067217, 0.073509, 0.099694, 0.080781, 0.073943,
1115 0.062822, 0.076611, 0.087672, 0.077632, 0.063959, 0.077754,
1116 0.084651, 0.097348, 0.063909, 0.055053, 0.087616, 0.074428,
1117 0.101165, 0.078446, 0.070683, 0.071162, 0.091516, 0.069641,
1118 0.070969, 0.069778, 0.083318, 0.066703, 0.074027, 0.067589,
1119 0.070620, 0.081307, 0.094406, 0.076188, 0.061663, 0.077561,
1120 0.076159, 0.071851, 0.074417, 0.076136, 0.069073, 0.075450,
1121 0.070218, 0.089709, 0.079974, 0.082077, 0.076979, 0.070497,
1122 0.076296, 0.087347, 0.073508, 0.088073, 0.067186, 0.059898,
1123 0.085495, 0.100152, 0.079468, 0.093783, 0.082680, 0.077272,
1124 0.081995, 0.068240, 0.071881, 0.073737, 0.072346, 0.079034,
1125 0.078721, 0.067518, 0.068196, 0.081738, 0.082814, 0.082480,
1126 0.069072, 0.066853, 0.081523, 0.073588, 0.082308, 0.091641,
1127 0.072747, 0.084317, 0.072190, 0.079372, 0.067424, 0.077450,
1128 0.062343, 0.090708, 0.065470, 0.086588, 0.071172, 0.066705,
1129 0.070670, 0.070407, 0.096111, 0.066738, 0.081249, 0.072882,
1130 0.075144, 0.060331, 0.074589, 0.076968, 0.085913, 0.072561,
1131 0.064645, 0.078742, 0.075670, 0.065984, 0.080932, 0.069898,
1132 0.065303, 0.096856, 0.057690, 0.065720, 0.066545, 0.068521,
1133 0.068278, 0.069245, 0.086643, 0.063401, 0.070933, 0.070752,
1134 0.066978, 0.058891, 0.070073, 0.078031, 0.082691, 0.101414,
1135 0.075814, 0.072790, 0.057622, 0.093002, 0.084660, 0.079216,
1136 0.070371, 0.070141, 0.076944, 0.067285, 0.078016, 0.077807,
1137 0.066668, 0.066459, 0.059962, 0.062774, 0.083450, 0.064554,
1138 0.067887, 0.064165, 0.072782, 0.067285, 0.052710, 0.096824,
1139 0.071931, 0.064190, 0.074442, 0.082647, 0.055797, 0.078632,
1140 0.061116, 0.063092, 0.049131, 0.074517, 0.069915, 0.079021,
1141 0.088656, 0.101045, 0.090432, 0.076460, 0.067451, 0.071381,
1142 0.069790, 0.077330, 0.077115, 0.084403, 0.065138, 0.056436,
1143 0.088024, 0.069893, 0.055985, 0.089655, 0.062911, 0.075311,
1144 0.086679, 0.093645, 0.068632, 0.064498, 0.057677, 0.081391,
1145 0.062781, 0.075467, 0.061314, 0.073394, 0.084462, 0.068470,
1146 0.071267, 0.060556, 0.072487, 0.063785, 0.079164, 0.070406,
1147 0.073394, 0.063168, 0.066968, 0.064125, 0.063954, 0.072504,
1148 0.058948, 0.057740, 0.083383, 0.067262};
1149 wDD = 0.098024;
1150 wND = 0.053452;
1151 wSD = -1;
1152
1153     Mmin = bin[0];
1154     Mmax = bin[nbin];
1155     if(M<Mmin || M>Mmax) return kTRUE;
1156
1157     Int_t ibin=nbin-1;
1158     for(Int_t i=1; i<=nbin; i++)
1159       if(M<=bin[i]) {
1160         ibin=i-1;
1161         //      printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);
1162         break;
1163       }
1164     wSD=w[ibin];
1165     return kTRUE;
1166
1167   }
1168
1169   return kFALSE;
1170 }
1171
1172
1173 //______________________________________________________________________________