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