1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
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.
22 #include "DPMcommon.h"
25 #include <TParticle.h>
27 #include <TDatabasePDG.h>
28 #include <TParticlePDG.h>
29 #include <TParticleClassPDG.h>
31 #include <TLorentzVector.h>
32 #include <TClonesArray.h>
33 #include "AliRunLoader.h"
34 #include "AliGenDPMjet.h"
35 #include "AliGenDPMjetEventHeader.h"
37 #include "AliDpmJetRndm.h"
38 #include "AliIonPDGCodes.h"
39 #include "AliHeader.h"
44 ClassImp(AliGenDPMjet)
46 //______________________________________________________________________________
47 AliGenDPMjet::AliGenDPMjet()
51 fMaxImpactParam(100.),
66 fTriggerMultiplicity(0),
67 fTriggerMultiplicityEta(0),
68 fTriggerMultiplicityPtMin(0),
71 fFragmentation(kFALSE),
72 fHeader(AliGenDPMjetEventHeader("DPMJET"))
77 AliDpmJetRndm::SetDpmJetRandom(GetRandom());
81 //______________________________________________________________________________
82 AliGenDPMjet::AliGenDPMjet(Int_t npart)
86 fMaxImpactParam(100.),
101 fTriggerMultiplicity(0),
102 fTriggerMultiplicityEta(0),
103 fTriggerMultiplicityPtMin(0),
106 fFragmentation(kFALSE),
107 fHeader(AliGenDPMjetEventHeader("DPMJET"))
111 // Default PbPb collisions at 5. 5 TeV
115 fTitle= "Particle Generator using DPMJET";
119 AliDpmJetRndm::SetDpmJetRandom(GetRandom());
122 AliGenDPMjet::AliGenDPMjet(const AliGenDPMjet &/*Dpmjet*/)
126 fMaxImpactParam(100.),
141 fTriggerMultiplicity(0),
142 fTriggerMultiplicityEta(0),
143 fTriggerMultiplicityPtMin(0),
146 fFragmentation(kFALSE),
151 // Dummy copy constructor
155 //______________________________________________________________________________
156 AliGenDPMjet::~AliGenDPMjet()
160 //______________________________________________________________________________
161 void AliGenDPMjet::Init()
165 if(fEnergyCMS>0. && fBeamEn<0.1) fBeamEn = fEnergyCMS/2;
166 SetMC(new TDPMjet(fProcess, fAProjectile, fZProjectile, fATarget, fZTarget,
167 fBeamEn,fEnergyCMS));
169 fDPMjet=(TDPMjet*) fMCEvGen;
170 if (fAProjectile == 1 && TMath::Abs(fZProjectile == 1)) fDPMjet->SetfIdp(1);
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);
180 fDPMjet->SetbRange(fMinImpactParam, fMaxImpactParam);
181 fDPMjet->SetPi0Decay(fPi0Decay);
182 fDPMjet->SetDecayAll(fDecayAll);
183 fDPMjet->SetFragmentProd(fFragmentation);
188 fDPMjet->Initialize();
192 //______________________________________________________________________________
193 void AliGenDPMjet::Generate()
195 // Generate one event
197 Double_t polar[3] = {0,0,0};
198 Double_t origin[3] = {0,0,0};
202 // converts from mm/c to s
203 const Float_t kconv = 0.001/2.999792458e8;
209 // Set collision vertex position
210 if (fVertexSmear == kPerEvent) Vertex();
214 // Generate one event
215 // --------------------------------------------------------------------------
218 // --------------------------------------------------------------------------
219 fDPMjet->GenerateEvent();
223 fDPMjet->ImportParticles(&fParticles,"All");
227 fGenImpPar = fDPMjet->GetBImpac();
229 if(TMath::Abs(fXingAngleY) > 1.e-10) BeamCrossAngle();
231 Int_t np = fParticles.GetEntriesFast();
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);
239 Int_t statusCode = iparticle->GetStatusCode();
241 // Initial state particle
245 if (fTriggerMultiplicityEta > 0 && TMath::Abs(iparticle->Eta()) > fTriggerMultiplicityEta)
248 if (iparticle->Pt() < fTriggerMultiplicityPtMin)
251 TParticlePDG* pdgPart = iparticle->GetPDG();
252 if (pdgPart && pdgPart->Charge() == 0)
258 if (multiplicity < fTriggerMultiplicity) continue;
259 Printf("Triggered on event with multiplicity of %d >= %d", multiplicity, fTriggerMultiplicity);
263 if(fkTuneForDiff && ( (TMath::Abs(fEnergyCMS - 900) < 1) || (TMath::Abs(fEnergyCMS - 2760) < 1) || (TMath::Abs(fEnergyCMS - 7000) < 1)) ) {
264 if(!CheckDiffraction() ) continue;
269 if (np == 0) continue;
272 Int_t* newPos = new Int_t[np];
273 Int_t* pSelected = new Int_t[np];
275 for (i = 0; i<np; i++) {
280 // First select parent particles
282 for (i = 0; i<np; i++) {
283 TParticle *iparticle = (TParticle *) fParticles.At(i);
285 // Is this a parent particle ?
287 if (Stable(iparticle)) continue;
289 Bool_t selected = kTRUE;
290 Bool_t hasSelectedDaughters = kFALSE;
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) &&
301 hasSelectedDaughters = DaughtersSelection(iparticle);
304 // Put particle on the stack if it is either selected or
305 // it is the mother of at least one seleted particle
307 if (selected || hasSelectedDaughters) {
311 } // particle loop parents
313 // Now select the final state particles
316 for (i=0; i<np; i++) {
317 TParticle *iparticle = (TParticle *) fParticles.At(i);
319 // Is this a final state particle ?
321 if (!Stable(iparticle)) continue;
323 Bool_t selected = kTRUE;
324 kf = iparticle->GetPdgCode();
325 ks = iparticle->GetStatusCode();
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;
333 // --------------------------------------------------------------------------
336 selected = KinematicSelection(iparticle,0)&&SelectFlavor(kf);
337 if (!fSpectators && selected) selected = (ks == 13 || ks == 14);
340 // Put particle on the stack if selected
346 } // particle loop final state
348 // Write particles to stack
350 for (i = 0; i<np; i++) {
351 TParticle * iparticle = (TParticle *) fParticles.At(i);
352 Bool_t hasMother = (iparticle->GetFirstMother()>=0);
355 kf = iparticle->GetPdgCode();
356 ks = iparticle->GetStatusCode();
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]
366 tof = fTime + kconv*iparticle->T();
369 TParticle* mother = 0;
371 imo = iparticle->GetFirstMother();
372 mother = (TParticle *) fParticles.At(imo);
373 imo = (mother->GetPdgCode() != 92 && mother->GetPdgCode() != 99999) ? newPos[imo] : -1;
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);
394 if (jev >= fNpart || fNpart == -1) {
400 SetHighWaterMark(nt);
403 //______________________________________________________________________________
404 Bool_t AliGenDPMjet::DaughtersSelection(TParticle* iparticle)
407 // Looks recursively if one of the daughters has been selected
409 // printf("\n Consider daughters %d:",iparticle->GetPdgCode());
413 Bool_t hasDaughters = (iparticle->GetFirstDaughter() >=0);
414 Bool_t selected = kFALSE;
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;
424 if (DaughtersSelection(jparticle)) {selected=kTRUE; break; }
434 //______________________________________________________________________________
435 Bool_t AliGenDPMjet::SelectFlavor(Int_t pid)
437 // Select flavor of particle
439 // 4: charm and beauty
446 Int_t ifl = TMath::Abs(pid/100);
447 if (ifl > 10) ifl/=10;
448 res = (fFlavor == ifl);
451 // This part if gamma writing is inhibited
453 res = res && (pid != kGamma && pid != kPi0);
458 //______________________________________________________________________________
459 Bool_t AliGenDPMjet::Stable(TParticle* particle)
461 // Return true for a stable particle
463 int st = particle->GetStatusCode();
464 if(st == 1 || st == -1) return kTRUE;
469 //______________________________________________________________________________
470 void AliGenDPMjet::MakeHeader()
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());
479 fHeader.SetCollisions(DTGLCP.ncp, DTGLCP.nct,
480 fDPMjet->GetProjWounded(),fDPMjet->GetTargWounded());
482 if(fProcDiff>0) fHeader.SetProcessType(fProcDiff);
483 else fHeader.SetProcessType(fDPMjet->GetProcessCode());
485 // Bookkeeping for kinematic bias
486 fHeader.SetTrials(fTrials);
488 fHeader.SetPrimaryVertex(fVertex);
489 fHeader.SetInteractionTime(fTime);
490 fHeader.SetNDiffractive(POEVT1.nsd1, POEVT1.nsd2, POEVT1.ndd);
491 // gAlice->SetGenEventHeader(fHeader);
493 fCollisionGeometry = &fHeader;
496 //______________________________________________________________________________
497 /*void AliGenDPMjet::AddHeader(AliGenEventHeader* fHeader)
499 // Add fHeader to container or runloader
501 fContainer->AddHeader(fHeader);
503 AliRunLoader::Instance()->GetHeader()->SetGenEventHeader(fHeader);
508 //______________________________________________________________________________
509 AliGenDPMjet& AliGenDPMjet::operator=(const AliGenDPMjet& /*rhs*/)
511 // Assignment operator
516 //______________________________________________________________________________
517 void AliGenDPMjet::FinishRun()
519 // Print run statistics
520 fDPMjet->Dt_Dtuout();
524 //______________________________________________________________________________
525 Bool_t AliGenDPMjet::CheckDiffraction()
530 Int_t np = fParticles.GetEntriesFast();
538 const Int_t kNstable=20;
539 const Int_t pdgStable[20] = {
542 12, // Electron Neutrino
562 for (Int_t i = 0; i < np; i++) {
563 TParticle * part = (TParticle *) fParticles.At(i);
565 Int_t statusCode = part->GetStatusCode();
567 // Initial state particle
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]) {
582 Double_t y = part->Y();
596 if(iPart1<0 || iPart2<0) return kFALSE;
601 TParticle * part1 = (TParticle *) fParticles.At(iPart1);
602 TParticle * part2 = (TParticle *) fParticles.At(iPart2);
604 Int_t pdg1 = part1->GetPdgCode();
605 Int_t pdg2 = part2->GetPdgCode();
609 if (pdg1 == 2212 && pdg2 == 2212)
617 if((AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.)>0.5) iPart = iPart2;
620 else if (pdg1 == 2212)
622 else if (pdg2 == 2212)
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));
637 Double_t Mmin, Mmax, wSD, wDD, wND;
639 if(!GetWeightsDiffraction(M, Mmin, Mmax, wSD, wDD, wND)) return kFALSE;
640 if(M>-1 && M<Mmin) return kFALSE;
644 Int_t procType=fDPMjet->GetProcessCode();//fPythia->GetMSTI(1);
645 if(procType== 4) return kFALSE;
647 if(procType== 7) proc0=1;
648 if(procType== 5 || procType== 6) proc0=0;
653 else if(proc0==1) proc=1;
655 if(proc==1 && (AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.) > wDD) return kFALSE;
656 if(proc==2 && (AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.) > wND) return kFALSE;
661 if(proc0!=0) fProcDiff = procType;
666 if(wSD<0) AliError("wSD<0 ! \n");
668 if((AliDpmJetRndm::GetDpmJetRandom())->Uniform(0.,1.)> wSD) return kFALSE;
670 // printf("iPart = %d\n", iPart);
672 if(iPart==iPart1) fProcDiff=5;
673 else if(iPart==iPart2) fProcDiff=6;
675 printf("EROOR: iPart!=iPart1 && iPart!=iPart2\n");
682 // -------------------------------------------------------
683 void AliGenDPMjet::SetIonPDGCodes()
685 // Defining PDG codes for the ions
686 AliIonPDGCodes *pdgcodes = new AliIonPDGCodes();
687 pdgcodes->AddParticlesToPdgDataBase();
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)
695 if(TMath::Abs(fEnergyCMS-900)<1 ){
696 const Int_t nbin=400;
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};
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};
839 if(M<Mmin || M>Mmax) return kTRUE;
842 for(Int_t i=1; i<=nbin; i++)
845 // printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);
852 else if(TMath::Abs(fEnergyCMS-2760)<1 ){
854 const Int_t nbin=400;
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};
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};
997 if(M<Mmin || M>Mmax) return kTRUE;
1000 for(Int_t i=1; i<=nbin; i++)
1003 // printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);
1010 else if(TMath::Abs(fEnergyCMS-7000)<1 ){
1012 const Int_t nbin=400;
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};
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};
1155 if(M<Mmin || M>Mmax) return kTRUE;
1158 for(Int_t i=1; i<=nbin; i++)
1161 // printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);
1173 //______________________________________________________________________________