BOSS 7.0.2
BESIII Offline Software System
Loading...
Searching...
No Matches
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 ()
 
 KKMC (const std::string &name, ISvcLocator *pSvcLocator)
 
StatusCode initialize ()
 
StatusCode execute ()
 
StatusCode finalize ()
 

Detailed Description

Definition at line 19 of file Generator/KKMC/KKMC-00-00-55/KKMC/KKMC.h.

Constructor & Destructor Documentation

◆ KKMC() [1/2]

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

Definition at line 104 of file KKMC.cxx.

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

◆ KKMC() [2/2]

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

Member Function Documentation

◆ execute() [1/2]

StatusCode KKMC::execute ( )

Definition at line 446 of file KKMC.cxx.

446 {
447 SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
448 int runNo=eventHeader->runNumber();
449 int event=eventHeader->eventNumber();
450 bool newRunFlag=0;
451 if(runNo != 0 && runNo != m_runNo){m_runNo=runNo;newRunFlag = true;}else{newRunFlag=false;}
452
453 //std::cout<<"runNo= "<<runNo<<" m_runNo "<<m_runNo<<" newRunFlag= "<<newRunFlag<<std::endl;
454
455 /****** Lianjin Wu add ******/
456 if(m_RdMeasuredEcms&& newRunFlag){// using cms energy of beam read from database
457 if(!ecmsSvc->isRunNoValid(runNo)){
458 std::cout<<"ERROR: Can not read the MeasuredEcms, try to turn off the ReadEcmsDB"<<std::endl;
459 return StatusCode::FAILURE;
460 }
461 double dbEcms=ecmsSvc->getEcms(runNo);
462 std::cout<<"INFO: Read the MeasuredEcms"<<std::endl;
463 m_cmsEnergy=dbEcms;
464 xwpar[0] = m_cmsEnergy;
465 KK2F_INITIALIZE(xwpar);
466 }
467 /****** ******/
468 // std::cout<<"The beam energy is "<<m_cmsEnergy<<", set for RunNo "<<runNo<<std::endl;
469
470 MsgStream log(msgSvc(), name());
471
472 log << MSG::INFO << "KKMC in execute()" << endreq;
473
474
475 HepMC::IO_HEPEVT HepEvtIO;
476
477
478 int KeySkip,iflag;
479 do {
480 HEPEVT_CLEAN();
481 KK2F_MAKE();
482 KK2F_GETKEYSKIP(KeySkip);
483 PSIPP_DDBARCUT(iflag);
484 } while (KeySkip != 0 || iflag ==0);
485
486 int KeyInter;
487 KK2F_GETEVTGENINTERFACE(KeyInter);
488
489 // PSIPP_DDBARCUT(iflag);
490 // if(iflag == 0) {return StatusCode::SUCCESS;}
491
492 if(KeyInter == 0) { // make photos correction
493 int Pos1, Pos2, KFfin, Nhep;
494 HEPEVT_NUMHEP(Nhep);
495 HEPEVT_GETF(Pos1);
496 HEPEVT_GETFBAR(Pos2);
497 HEPEVT_GETKFFIN(KFfin);
498 int Posn = Pos1;
499 if(Pos2 > Posn) Posn = Pos2;
500 Posn = Posn + 1;
501 if(KFfin < 10) Posn = Posn + 1;
502 for(int ip = Posn; ip <= Nhep; ip++) HEPEVT_SETPHOTOSFLAGTRUE(ip);
503 for(int ip = Posn; ip <= Nhep; ip++) PHOTOS(ip);
504 PYHEPC(2);
505 }
506
507 //sleep(5);
508
509 m_numberEvent += 1;
510
511 if(m_numberEvent <= m_numberEventPrint) PYLIST(1);
512 log << MSG::INFO <<" " <<m_numberEvent<<"th event was generated !!" <<endreq;
513
514 PYHEPC(1);
515
516 HepMC::GenEvent* evt = HepEvtIO.read_next_event();
517 evt->set_event_number(m_numberEvent);
518 evt->set_signal_process_id(1);
519 // evt->print();
520
521 // Check if the McCollection already exists
522 SmartDataPtr<McGenEventCol> anMcCol(eventSvc(), "/Event/Gen");
523 if (anMcCol!=0) {
524 // Add event to existing collection
525 MsgStream log(messageService(), name());
526 log << MSG::INFO << "Add McGenEvent to existing collection" << endreq;
527 McGenEvent* mcEvent = new McGenEvent(evt);
528 anMcCol->push_back(mcEvent);
529 } else {
530 // Create Collection and add to the transient store
531 McGenEventCol *mcColl = new McGenEventCol;
532 McGenEvent* mcEvent = new McGenEvent(evt);
533 mcColl->push_back(mcEvent);
534 StatusCode sc = eventSvc()->registerObject("/Event/Gen",mcColl);
535 if (sc != StatusCode::SUCCESS) {
536 log << MSG::ERROR << "Could not register McGenEvent" << endreq;
537 delete mcColl;
538 delete evt;
539 delete mcEvent;
540 return StatusCode::FAILURE;
541 } else {
542 // log << MSG::INFO << "McGenEventCol created and " << npart <<" particles stored in McGenEvent" << endreq;
543 }
544 }
545
546
547
548
549 return StatusCode::SUCCESS;
550
551}
#define HEPEVT_CLEAN()
Definition: BesBdkRc.cxx:100
int runNo
Definition: DQA_TO_DB.cxx:12
#define PYHEPC(ICONV)
Definition: KKMC.cxx:51
#define PSIPP_DDBARCUT(KEY)
Definition: KKMC.cxx:40
#define KK2F_GETEVTGENINTERFACE(KEY)
Definition: KKMC.cxx:56
#define HEPEVT_GETKFFIN(KFIN)
Definition: KKMC.cxx:65
#define PYLIST(LIST)
Definition: KKMC.cxx:48
#define HEPEVT_NUMHEP(Nhep)
Definition: KKMC.cxx:59
#define KK2F_MAKE()
Definition: KKMC.cxx:34
#define KK2F_GETKEYSKIP(KEY)
Definition: KKMC.cxx:37
#define HEPEVT_GETF(POS)
Definition: KKMC.cxx:61
#define KK2F_INITIALIZE(XPAR)
Definition: KKMC.cxx:29
#define HEPEVT_SETPHOTOSFLAGTRUE(IP)
Definition: KKMC.cxx:67
#define HEPEVT_GETFBAR(POS)
Definition: KKMC.cxx:63
#define PHOTOS(IP)
Definition: KKMC.cxx:69
virtual bool isRunNoValid(int runNo)=0
virtual double getEcms(int runNo)=0

◆ execute() [2/2]

StatusCode KKMC::execute ( )

◆ finalize() [1/2]

StatusCode KKMC::finalize ( )

Definition at line 553 of file KKMC.cxx.

553 {
554
555 MsgStream log(msgSvc(), name());
556
557 log << MSG::INFO << "KKMC in finalize()" << endreq;
558
560 double xSecPb, xErrPb;
561 KK2F_GETXSECMC(xSecPb, xErrPb);
562
563 log << MSG::INFO << "Total MC Xsec = " << xSecPb << " +/- " << xErrPb << endreq;
564 return StatusCode::SUCCESS;
565}
#define KK2F_GETXSECMC(xsecpb, xerrpb)
Definition: KKMC.cxx:45
#define KK2F_FINALIZE()
Definition: KKMC.cxx:43

◆ finalize() [2/2]

StatusCode KKMC::finalize ( )

◆ initialize() [1/2]

StatusCode KKMC::initialize ( )

Definition at line 225 of file KKMC.cxx.

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

◆ initialize() [2/2]

StatusCode KKMC::initialize ( )

The documentation for this class was generated from the following files: