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