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