CGEM BOSS 6.6.5.g
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
KKMC Class Reference

#include <KKMC.h>

+ Inheritance diagram for KKMC:

Public Member Functions

 KKMC (const std::string &name, ISvcLocator *pSvcLocator)
 
StatusCode initialize ()
 
StatusCode execute ()
 
StatusCode finalize ()
 

Detailed Description

Definition at line 18 of file KKMC.h.

Constructor & Destructor Documentation

◆ KKMC()

KKMC::KKMC ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Definition at line 106 of file KKMC.cxx.

106 : Algorithm(name, pSvcLocator)
107{
108 m_numberEvent = 0;
109
110 m_kkseed.clear();
111 m_kkseed.push_back(123456);
112 m_kkseed.push_back(1);
113 m_kkseed.push_back(0);
114
115
116 declareProperty("NumberOfEventPrinted", m_numberEventPrint=100);
117 declareProperty("InitializedSeed", m_kkseed);
118 declareProperty("CMSEnergy", m_cmsEnergy = 3.773);
119 declareProperty("ReadEcmsFromDB", m_ReadEcmsFromDB = false);
120 declareProperty("BeamEnergySpread", m_cmsEnergySpread = 0.0013); //it should be beam energy spread,pingrg-2009-09-24
121 declareProperty("GenerateResonance", m_generateResonance = true);
122 declareProperty("GenerateContinuum", m_generateContinuum = true);
123 declareProperty("GenerateDownQuark", m_generateDownQuark = true);
124 declareProperty("GenerateUpQuark", m_generateUpQuark = true);
125 declareProperty("GenerateStrangeQuark", m_generateStrangeQuark = true);
126 declareProperty("GenerateCharmQuark", m_generateCharmQuark = true);
127 declareProperty("GenerateBottomQuark", m_generateBottomQuark = false);
128 declareProperty("GenerateMuonPair", m_generateMuonPair = true);
129 declareProperty("GenerateTauPair", m_generateTauPair = true);
130 declareProperty("GenerateRho", m_generateRho = true);
131 declareProperty("GenerateOmega", m_generateOmega = true);
132 declareProperty("GeneratePhi", m_generatePhi = true);
133 declareProperty("GenerateJPsi", m_generateJPsi = true);
134 declareProperty("GeneratePsiPrime", m_generatePsiPrime = true);
135 declareProperty("GeneratePsi3770", m_generatePsi3770 = true);
136 declareProperty("GeneratePsi4030", m_generatePsi4030 = true);
137 declareProperty("GeneratePsi4160", m_generatePsi4160 = true);
138 declareProperty("GeneratePsi4415", m_generatePsi4415 = true);
139 declareProperty("GeneratePsi4260", m_generatePsi4260 = true);
140 declareProperty("ThresholdCut", m_DdbarCutPsi3770 = -3.0); //generate DDbar decay, pingrg-2009-10-14
141 declareProperty("TagISR", m_isrtag = false); //Tag ISR photon, false: ID=22, true: ID=-22, pingrg-2010-6-30
142 declareProperty("TagFSR", m_fsrtag = false); //Tag FSR photon, false: ID=22, true: ID=-22, pingrg-2010-6-30
143 declareProperty("ModeIndexExpXS", m_ich = -10); //mode index using the measured cross section
144
145 m_paramRho.clear(); m_paramRho.push_back(0.77457e0); m_paramRho.push_back(147.65e-3); m_paramRho.push_back(6.89e-6);
146 m_paramRh2.clear(); m_paramRh2.push_back(1.465e0); m_paramRh2.push_back(310e-3); m_paramRh2.push_back(0.0e-6);
147 m_paramRh3.clear(); m_paramRh3.push_back(1.700e0); m_paramRh3.push_back(240e-3); m_paramRh3.push_back(0.0e-6);
148 m_paramOme.clear(); m_paramOme.push_back(0.78194e0);m_paramOme.push_back(8.41e-3);m_paramOme.push_back(0.60e-6);
149 m_paramOm2.clear(); m_paramOm2.push_back(1.419e0);m_paramOm2.push_back(174e-3);m_paramOm2.push_back(0.00e-6);
150 m_paramOm3.clear(); m_paramOm3.push_back(1.649e0);m_paramOm3.push_back(220e-3);m_paramOm3.push_back(0.00e-6);
151 m_paramPhi.clear(); m_paramPhi.push_back(1.01942e0);m_paramPhi.push_back(4.46e-3);m_paramPhi.push_back(1.33e-6);
152 m_paramPh2.clear(); m_paramPh2.push_back(1.680e0);m_paramPh2.push_back(150e-3);m_paramPh2.push_back(0.00e-6);
153 m_paramPsi.clear(); m_paramPsi.push_back(3.096916e0);m_paramPsi.push_back(0.0929e-3);m_paramPsi.push_back(5.40e-6);
154 m_paramPs2.clear(); m_paramPs2.push_back(3.686109e0);m_paramPs2.push_back(0.304e-3);m_paramPs2.push_back(2.12e-6);
155 m_paramPs3.clear(); m_paramPs3.push_back(3.77315e0);m_paramPs3.push_back(27.2e-3);m_paramPs3.push_back(0.26e-6);
156 m_paramPs4.clear(); m_paramPs4.push_back(4.039e0);m_paramPs4.push_back(80e-3);m_paramPs4.push_back(0.86e-6);
157 m_paramPs5.clear(); m_paramPs5.push_back(4.153e0);m_paramPs5.push_back(103e-3);m_paramPs5.push_back(0.83e-6);
158 m_paramPs6.clear(); m_paramPs6.push_back(4.421e0);m_paramPs6.push_back(62e-3);m_paramPs6.push_back(0.58e-6);
159 m_paramPs7.clear(); m_paramPs7.push_back(4.263e0);m_paramPs7.push_back(95e-3);m_paramPs7.push_back(0.47e-6);
160 m_paramPs8.clear(); m_paramPs8.push_back(3.872e0);m_paramPs8.push_back(100e-3);m_paramPs8.push_back(0.00e-6);
161 m_paramUps.clear(); m_paramUps.push_back(9.46030e0); m_paramUps.push_back(0.0525e-3); m_paramUps.push_back(1.32e-6);
162 m_paramUp2.clear(); m_paramUp2.push_back(10.02326e0); m_paramUp2.push_back(0.044e-3); m_paramUp2.push_back(0.52e-6);
163 m_paramUp3.clear(); m_paramUp3.push_back(10.3552e0); m_paramUp3.push_back(0.026e-3); m_paramUp3.push_back(0.00e-6);
164 m_paramUp4.clear(); m_paramUp4.push_back(10.580e0); m_paramUp4.push_back(14e-3); m_paramUp4.push_back(0.248e-6);
165 m_paramUp5.clear(); m_paramUp5.push_back(10.865e0); m_paramUp5.push_back(110e-3); m_paramUp5.push_back(0.31e-6);
166 m_paramUp6.clear(); m_paramUp6.push_back(11.019); m_paramUp6.push_back(79e-3); m_paramUp6.push_back(0.13e-6);
167 m_paramZeta.clear(); m_paramZeta.push_back(91.1876e0); m_paramZeta.push_back(2.4952e0); m_paramZeta.push_back(0.08391e0);
168 m_paramW.clear(); m_paramW.push_back(80.43); m_paramW.push_back(2.11e0);
169
170 declareProperty("ResParameterRho", m_paramRho);
171 declareProperty("ResParameterRh2", m_paramRh2);
172 declareProperty("ResParameterRh3", m_paramRh3);
173 declareProperty("ResParameterOme", m_paramOme);
174 declareProperty("ResParameterOm2", m_paramOm2);
175 declareProperty("ResParameterOm3", m_paramOm3);
176 declareProperty("ResParameterPhi", m_paramPhi);
177 declareProperty("ResParameterPh2", m_paramPh2);
178 declareProperty("ResParameterPsi", m_paramPsi);
179 declareProperty("ResParameterPs2", m_paramPs2);
180 declareProperty("ResParameterPs3", m_paramPs3);
181 declareProperty("ResParameterPs4", m_paramPs4);
182 declareProperty("ResParameterPs5", m_paramPs5);
183 declareProperty("ResParameterPs6", m_paramPs6);
184 declareProperty("ResParameterPs7", m_paramPs7);
185 declareProperty("ResParameterPs8", m_paramPs8);
186 declareProperty("ResParameterUps", m_paramUps);
187 declareProperty("ResParameterUp2", m_paramUp2);
188 declareProperty("ResParameterUp3", m_paramUp3);
189 declareProperty("ResParameterUp4", m_paramUp4);
190 declareProperty("ResParameterUp5", m_paramUp5);
191 declareProperty("ResParameterUp6", m_paramUp6);
192 declareProperty("ResParameterZeta", m_paramZeta);
193 declareProperty("ResParameterW", m_paramW);
194
195 // psi(3770) decay
196 declareProperty("Psi3770toNonDDb", m_ps3toNonDDb = 0.0);
197 declareProperty("Psi3770RatioOfD0toDp", m_ps3D0toDp = 0.563);
198 // psi(4030) decay
199 declareProperty("Psi4030toD0D0b", m_ps4toD0D0b = 0.0227);
200 declareProperty("Psi4030toDpDm", m_ps4toDpDm = 0.0167);
201 declareProperty("Psi4030toDsDs", m_ps4toDsDs = 0.0383);
202 declareProperty("Psi4030toD0D0Star", m_ps4toD0D0Star = 0.2952);
203 declareProperty("Psi4030toDpDmStar", m_ps4toDpDmStar = 0.2764);
204 declareProperty("Psi4030toD0StarD0Star", m_ps4toD0StarD0Star=0.2476);
205 declareProperty("Psi4030toDpStarDmStar", m_ps4toDpStarDmStar=0.1041);
206 // psi(4160) decay
207 declareProperty("Psi4160toD0D0b", m_ps5toD0D0b = 0.0190);
208 declareProperty("Psi4160toDpDm", m_ps5toDpDm = 0.0180);
209 declareProperty("Psi4160toDsDs", m_ps5toDsDs = 0.0488);
210 declareProperty("Psi4160toD0D0Star", m_ps5toD0D0Star = 0.1248);
211 declareProperty("Psi4160toDpDmStar", m_ps5toDpDmStar = 0.1240);
212 declareProperty("Psi4160toDsDsStar", m_ps5toDsDsStar = 0.0820);
213 declareProperty("Psi4160toD0StarD0Star", m_ps5toD0StarD0Star=0.3036);
214 declareProperty("Psi4160toDpStarDmStar", m_ps5toDpStarDmStar=0.2838);
215 // interface to EvtGen
216 declareProperty("ParticleDecayThroughEvtGen", m_evtGenDecay = true);
217 declareProperty("RadiationCorrection", m_radiationCorrection = true);
218 //interface set pythia pars
219 m_pypars.clear();
220 declareProperty("setPythiaPars", m_pypars);
221}

Member Function Documentation

◆ execute()

StatusCode KKMC::execute ( )

Definition at line 432 of file KKMC.cxx.

432 {
433 SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
434 int runNo=eventHeader->runNumber();
435 int event=eventHeader->eventNumber();
436 bool newRunFlag=0;
437 if(runNo != 0 && runNo != m_runNo){m_runNo=runNo;newRunFlag = true;}else{newRunFlag=false;}
438
439 //std::cout<<"runNo= "<<runNo<<" m_runNo "<<m_runNo<<" newRunFlag= "<<newRunFlag<<std::endl;
440
441 if(m_ReadEcmsFromDB && newRunFlag){// using cms energy of beam read from database
442 runNo = std::abs(runNo);
443 ReadEb theEb(runNo);
444 double dbEcms = theEb.getEcms();
445 double dbXangle = theEb.getXangle();
446 m_cmsEnergy = dbEcms; //*cos(dbXangle); //here the dbEcms is the enegy in the center-mass system of psi"
447 std::cout<<"The beam energy is "<<dbEcms<<", set for RunNo "<<runNo<<std::endl;
448 xwpar[0] = m_cmsEnergy;
449 KK2F_INITIALIZE(xwpar);
450 }
451 MsgStream log(msgSvc(), name());
452
453 log << MSG::INFO << "KKMC in execute()" << endreq;
454
455
456 HepMC::IO_HEPEVT HepEvtIO;
457
458
459 int KeySkip,iflag;
460 do {
461 HEPEVT_CLEAN();
462 KK2F_MAKE();
463 KK2F_GETKEYSKIP(KeySkip);
464 PSIPP_DDBARCUT(iflag);
465 } while (KeySkip != 0 || iflag ==0);
466
467 int KeyInter;
468 KK2F_GETEVTGENINTERFACE(KeyInter);
469
470 // PSIPP_DDBARCUT(iflag);
471 // if(iflag == 0) {return StatusCode::SUCCESS;}
472
473 if(KeyInter == 0) { // make photos correction
474 int Pos1, Pos2, KFfin, Nhep;
475 HEPEVT_NUMHEP(Nhep);
476 HEPEVT_GETF(Pos1);
477 HEPEVT_GETFBAR(Pos2);
478 HEPEVT_GETKFFIN(KFfin);
479 int Posn = Pos1;
480 if(Pos2 > Posn) Posn = Pos2;
481 Posn = Posn + 1;
482 if(KFfin < 10) Posn = Posn + 1;
483 for(int ip = Posn; ip <= Nhep; ip++) HEPEVT_SETPHOTOSFLAGTRUE(ip);
484 for(int ip = Posn; ip <= Nhep; ip++) PHOTOS(ip);
485 PYHEPC(2);
486 }
487
488 //sleep(5);
489
490 m_numberEvent += 1;
491
492 if(m_numberEvent <= m_numberEventPrint) PYLIST(1);
493 log << MSG::INFO <<" " <<m_numberEvent<<"th event was generated !!" <<endreq;
494
495 PYHEPC(1);
496
497 HepMC::GenEvent* evt = HepEvtIO.read_next_event();
498 evt->set_event_number(m_numberEvent);
499 evt->set_signal_process_id(1);
500 // evt->print();
501
502 // Check if the McCollection already exists
503 SmartDataPtr<McGenEventCol> anMcCol(eventSvc(), "/Event/Gen");
504 if (anMcCol!=0) {
505 // Add event to existing collection
506 MsgStream log(messageService(), name());
507 log << MSG::INFO << "Add McGenEvent to existing collection" << endreq;
508 McGenEvent* mcEvent = new McGenEvent(evt);
509 anMcCol->push_back(mcEvent);
510 } else {
511 // Create Collection and add to the transient store
512 McGenEventCol *mcColl = new McGenEventCol;
513 McGenEvent* mcEvent = new McGenEvent(evt);
514 mcColl->push_back(mcEvent);
515 StatusCode sc = eventSvc()->registerObject("/Event/Gen",mcColl);
516 if (sc != StatusCode::SUCCESS) {
517 log << MSG::ERROR << "Could not register McGenEvent" << endreq;
518 delete mcColl;
519 delete evt;
520 delete mcEvent;
521 return StatusCode::FAILURE;
522 } else {
523 // log << MSG::INFO << "McGenEventCol created and " << npart <<" particles stored in McGenEvent" << endreq;
524 }
525 }
526
527
528
529
530 return StatusCode::SUCCESS;
531
532}
HEPEVT_CLEAN
#define HEPEVT_CLEAN()
Definition: BesBdkRc.cxx:100
runNo
int runNo
Definition: DedxCalibRunByRun.cxx:14
PYHEPC
#define PYHEPC(ICONV)
Definition: KKMC.cxx:53
PSIPP_DDBARCUT
#define PSIPP_DDBARCUT(KEY)
Definition: KKMC.cxx:42
KK2F_GETEVTGENINTERFACE
#define KK2F_GETEVTGENINTERFACE(KEY)
Definition: KKMC.cxx:58
HEPEVT_GETKFFIN
#define HEPEVT_GETKFFIN(KFIN)
Definition: KKMC.cxx:67
PYLIST
#define PYLIST(LIST)
Definition: KKMC.cxx:50
HEPEVT_NUMHEP
#define HEPEVT_NUMHEP(Nhep)
Definition: KKMC.cxx:61
KK2F_MAKE
#define KK2F_MAKE()
Definition: KKMC.cxx:36
KK2F_GETKEYSKIP
#define KK2F_GETKEYSKIP(KEY)
Definition: KKMC.cxx:39
HEPEVT_GETF
#define HEPEVT_GETF(POS)
Definition: KKMC.cxx:63
KK2F_INITIALIZE
#define KK2F_INITIALIZE(XPAR)
Definition: KKMC.cxx:31
HEPEVT_SETPHOTOSFLAGTRUE
#define HEPEVT_SETPHOTOSFLAGTRUE(IP)
Definition: KKMC.cxx:69
HEPEVT_GETFBAR
#define HEPEVT_GETFBAR(POS)
Definition: KKMC.cxx:65
PHOTOS
#define PHOTOS(IP)
Definition: KKMC.cxx:71
McGenEventCol
ObjectVector< McGenEvent > McGenEventCol
Definition: McGenEvent.h:39
msgSvc
IMessageSvc * msgSvc()
Definition: ServiceAccessor.h:13
ReadEb
Definition: ReadEb.h:17

◆ finalize()

StatusCode KKMC::finalize ( )

Definition at line 534 of file KKMC.cxx.

534 {
535
536 MsgStream log(msgSvc(), name());
537
538 log << MSG::INFO << "KKMC in finalize()" << endreq;
539
540 KK2F_FINALIZE();
541 double xSecPb, xErrPb;
542 KK2F_GETXSECMC(xSecPb, xErrPb);
543
544 log << MSG::INFO << "Total MC Xsec = " << xSecPb << " +/- " << xErrPb << endreq;
545 return StatusCode::SUCCESS;
546}
KK2F_GETXSECMC
#define KK2F_GETXSECMC(xsecpb, xerrpb)
Definition: KKMC.cxx:47
KK2F_FINALIZE
#define KK2F_FINALIZE()
Definition: KKMC.cxx:45

◆ initialize()

StatusCode KKMC::initialize ( )

Definition at line 223 of file KKMC.cxx.

223 {
224
225 MsgStream log(msgSvc(), name());
226
227 log << MSG::INFO << "KKMC in initialize()" << endreq;
228
229 //set Bes unified random engine
230 static const bool CREATEIFNOTTHERE(true);
231 StatusCode RndmStatus = service("BesRndmGenSvc", p_BesRndmGenSvc, CREATEIFNOTTHERE);
232 if (!RndmStatus.isSuccess() || 0 == p_BesRndmGenSvc)
233 {
234 log << MSG::ERROR << " Could not initialize Random Number Service" << endreq;
235 return RndmStatus;
236 }
237 CLHEP::HepRandomEngine* engine = p_BesRndmGenSvc->GetEngine("KKMC");
238 KKMCRandom::setRandomEngine(engine);
239 //---
240 if(m_ich == -2 || m_ich>=50 ) {// ISR exclusive default set particle as Psi(4260)
241 m_generatePsi4260 = true;
242 m_generateJPsi = 0;
243 m_generatePsiPrime = 0;
244 m_generatePsi3770 = 0;
245 m_generatePsi4030 = 0;
246 m_generatePsi4160 = 0;
247 m_generatePsi4415 = 0;
248 }
249
250 WHYM_SETDEF(xwpar);
251 xwpar[0] = m_cmsEnergy; // c.m.s energy
252 xwpar[1] = m_cmsEnergySpread; // energy spread of c.m.s.
253 xwpar[3] = 6.0; // fortran output unit
254
255 if(m_generateResonance) // flag indicate to generate Resonance data
256 xwpar[12] = 1.0;
257 else
258 xwpar[12] = 0.0;
259
260 if(m_generateContinuum ) // generate continuum
261 xwpar[3000] = 1.0;
262 else
263 xwpar[3000] = 0.0;
264
265 if(m_generateDownQuark) // d quark production
266 xwpar[400] = 1.0;
267 else
268 xwpar[400] = 0.0;
269
270 if(m_generateUpQuark) // u quark production
271 xwpar[401] = 1.0;
272 else
273 xwpar[401] = 0.0;
274
275 if(m_generateStrangeQuark) // s quark production
276 xwpar[402] = 1.0;
277 else
278 xwpar[402] = 0.0;
279
280 if(m_generateCharmQuark) // c quark production
281 xwpar[403] = 1.0;
282 else
283 xwpar[403] = 0.0;
284
285 if(m_generateBottomQuark) // b quark production
286 xwpar[404] = 1.0;
287 else
288 xwpar[404] = 0.0;
289
290
291 if(m_generateMuonPair) // e+ e- --> mu+ mu-
292 xwpar[412] = 1.0;
293 else
294 xwpar[412] = 0.0;
295
296 if(m_generateTauPair) // e+ e- --> tau+ tau-
297 xwpar[414] = 1.0;
298 else
299 xwpar[414] = 0.0;
300 int keyuds = 0;
301 if(m_generateRho) keyuds |= 1;
302 if(m_generateOmega) keyuds |= 2;
303 if(m_generatePhi) keyuds |= 4;
304
305 int keycharm = 0;
306 if(m_generateJPsi) keycharm |= 1;
307 if(m_generatePsiPrime) keycharm |= 2;
308 if(m_generatePsi3770) keycharm |= 4;
309 if(m_generatePsi4030) keycharm |= 8;
310 if(m_generatePsi4160) keycharm |= 16;
311 if(m_generatePsi4415) keycharm |= 32;
312 if(m_generatePsi4260) keycharm |= 64;
313
314 // resonant parameters
315 int offset = 3100;
316 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramRho[i];
317 offset = offset + 3;
318 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramRh2[i];
319 offset = offset + 3;
320 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramRh3[i];
321 offset = offset + 3;
322 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramOme[i];
323 offset = offset + 3;
324 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramOm2[i];
325 offset = offset + 3;
326 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramOm3[i];
327 offset = offset + 3;
328 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPhi[i];
329 offset = offset + 3;
330 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPh2[i];
331 offset = offset + 3;
332 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPsi[i];
333 offset = offset + 3;
334 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs2[i];
335 offset = offset + 3;
336 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs3[i];
337 offset = offset + 3;
338 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs4[i];
339 offset = offset + 3;
340 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs5[i];
341 offset = offset + 3;
342 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs6[i];
343 offset = offset + 3;
344 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs7[i];
345 offset = offset + 3;
346 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramPs8[i];
347 offset = offset + 3;
348 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramUps[i];
349 offset = offset + 3;
350 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramUp2[i];
351 offset = offset + 3;
352 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramUp3[i];
353 offset = offset + 3;
354 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramUp4[i];
355 offset = offset + 3;
356 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramUp5[i];
357 offset = offset + 3;
358 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramUp6[i];
359 offset = offset + 3;
360 for(int i = 0; i < 3; i++) xwpar[offset+i] = m_paramZeta[i];
361 offset = offset + 3;
362 for(int i = 0; i < 2; i++) xwpar[offset+i] = m_paramW[i];
363 // offset = offset + 2;
364
365 xwpar[3001] = keyuds + 0.0;
366 xwpar[3002] = keycharm + 0.0;
367
368 // psi(3770) decay
369 offset = 3200;
370 xwpar[offset + 0] = m_ps3toNonDDb;
371 xwpar[offset + 1] = m_ps3D0toDp;
372 DDBARMASS.ddbarmassCUT= m_DdbarCutPsi3770;
373 MODEXS.ich = m_ich;
374 //tag ISR or FSR photon
375 if(m_isrtag){PHOTONTAG.isrtag=1;} else {PHOTONTAG.isrtag=0;}
376 if(m_fsrtag){PHOTONTAG.fsrtag=1;} else {PHOTONTAG.fsrtag=0;}
377
378 // psi(4030) decay
379 offset = 3210;
380 xwpar[offset + 0] = m_ps4toD0D0b;
381 xwpar[offset + 1] = m_ps4toDpDm;
382 xwpar[offset + 2] = m_ps4toDsDs;
383 xwpar[offset + 3] = m_ps4toD0D0Star;
384 xwpar[offset + 4] = m_ps4toDpDmStar;
385 xwpar[offset + 5] = m_ps4toD0StarD0Star;
386 xwpar[offset + 6] = m_ps4toDpStarDmStar;
387 // psi(4160) decay
388 offset = 3220;
389 xwpar[offset + 0] = m_ps5toD0D0b;
390 xwpar[offset + 1] = m_ps5toDpDm;
391 xwpar[offset + 2] = m_ps5toDsDs;
392 xwpar[offset + 3] = m_ps5toD0D0Star;
393 xwpar[offset + 4] = m_ps5toDpDmStar;
394 xwpar[offset + 5] = m_ps5toDsDsStar;
395 xwpar[offset + 6] = m_ps5toD0StarD0Star;
396 xwpar[offset + 7] = m_ps5toDpStarDmStar;
397
398 if(!m_radiationCorrection) {
399 xwpar[19] = 0;
400 xwpar[20] = 0;
401 xwpar[26] = 0;
402 }
403
404
405 KK2F_INITIALIZE(xwpar);
406 MY_PYUPD();
407 PYUPDA(1, 22);
408 //for pythia parameter tunning,pingrg-2013-6-9
409 for(int i=0;i<m_pypars.size();i++){
410 pygive_(m_pypars[i].c_str(),strlen(m_pypars[i].c_str()));
411 }
412
413 if(m_evtGenDecay) {
414 KK2F_SETEVTGENINTERFACE(1);
415 TURNOFFTAUDECAY();
416 } else {
417 KK2F_SETEVTGENINTERFACE(0);
418 PHOINI();
419 }
420
421 HepMC::HEPEVT_Wrapper::set_sizeof_real(8);
422 HepMC::HEPEVT_Wrapper::set_sizeof_int(4);
423 HepMC::HEPEVT_Wrapper::set_max_number_entries(4000);
424 // std::cout << "max entries = " << HepMC::HEPEVT_Wrapper::max_number_entries() <<std::endl;
425 // std::cout << "size of real = " << HepMC::HEPEVT_Wrapper::sizeof_real() <<std::endl;
426
427
428 log <<MSG::INFO<< "Finish KKMC initialize()" <<endreq;
429 return StatusCode::SUCCESS;
430}
pygive_
void pygive_(const char *cnfgstr, int length)
MY_PYUPD
#define MY_PYUPD()
Definition: KKMC.cxx:28
KK2F_SETEVTGENINTERFACE
#define KK2F_SETEVTGENINTERFACE(KEY)
Definition: KKMC.cxx:56
PHOINI
#define PHOINI()
Definition: KKMC.cxx:73
WHYM_SETDEF
#define WHYM_SETDEF(XPAR)
Definition: KKMC.cxx:26
DDBARMASS
#define DDBARMASS
Definition: KKMC.cxx:84
PHOTONTAG
#define PHOTONTAG
Definition: KKMC.cxx:90
PYUPDA
#define PYUPDA(MUPDA, LFN)
Definition: KKMC.cxx:79
MODEXS
#define MODEXS
Definition: KKMC.cxx:96
TURNOFFTAUDECAY
#define TURNOFFTAUDECAY()
Definition: KKMC.cxx:76
IBesRndmGenSvc::GetEngine
virtual CLHEP::HepRandomEngine * GetEngine(const std::string &StreamName)=0
Interface to the CLHEP engine.
KKMCRandom::setRandomEngine
static void setRandomEngine(CLHEP::HepRandomEngine *randomEngine)
Definition: KKMCRandom.cxx:32
EvtCyclic3::c_str
char * c_str(Index i)
Definition: EvtCyclic3.cc:252
The documentation for this class was generated from the following files: