@@ -54,7 +54,6 @@ protected:
5454 while (!genOk ) {
5555 genOk = GeneratorPythia8 ::generateEvent ();
5656 if (mGeneratedEvents % mInverseTriggerRatio == 0 ) {
57- std ::cout <<"generating event with injected signals" <<std ::endl ;
5857 bool found = false;
5958 while (!found ) {
6059 mGeneratorParam -> generateEvent ();
@@ -70,7 +69,6 @@ protected:
7069 }
7170 }
7271 mGeneratedEvents ++ ;
73- std ::cout <<"generated events: " <<mGeneratedEvents <<std ::endl ;
7472 return true;
7573 }
7674
@@ -80,9 +78,7 @@ protected:
8078
8179 bool genOk = false;
8280 if ((mGeneratedEvents - 1 ) % mInverseTriggerRatio == 0 ){ // add injected prompt signals to the stack
83- // mGeneratorParam->importParticles();
8481 int originalSize = mParticles .size ();
85- std ::cout <<"adding " <<mGeneratorParam -> getParticles ().size ()<<" particles to the stack" <<std ::endl ;
8682 for (int ipart = 0 ; ipart < mGeneratorParam -> getParticles ().size (); ipart ++ ){
8783 TParticle part = TParticle (mGeneratorParam -> getParticles ().at (ipart ));
8884 if (part .GetFirstMother () >= 0 ) part .SetFirstMother (part .GetFirstMother () + originalSize );
@@ -98,16 +94,11 @@ protected:
9894 }
9995
10096 bool findSignalInAcceptance () {
101- std ::cout <<"loop over" << mGeneratorParam -> getParticles ().size ()<<" particles" <<std ::endl ;
102- for (int pdg : mSignalsPDGs ) {
103- std ::cout <<"signal pdg: " <<pdg <<std ::endl ;
104- }
10597 for (int ipart = 0 ; ipart < mGeneratorParam -> getParticles ().size (); ipart ++ ) {
10698 TParticle part = TParticle (mGeneratorParam -> getParticles ().at (ipart ));
10799 // make sure all signals are in the acceptance
108100 for (int pdg : mSignalsPDGs ) {
109101 if (part .GetPdgCode () == pdg ) {
110- std ::cout <<"found signal with pdg: " <<part .GetPdgCode ()<<", mother: " <<part .GetFirstMother ()<<std ::endl ;
111102 if (part .GetFirstMother () == -1 ) {
112103 if (part .Y () < mHadronRapidityMin || part .Y () > mHadronRapidityMax ) {
113104 return false;
@@ -116,15 +107,13 @@ protected:
116107 }
117108 }
118109 }
119- std ::cout <<"generated signal in acceptance" <<std ::endl ;
120110 return true;
121111 }
122112
123113private :
124114 Generator * mGeneratorParam ;
125115 unsigned long long mGeneratedEvents ;
126116 int mInverseTriggerRatio ;
127- // Pythia8::Pythia pythiaMBgen; // minimum bias event
128117 std ::vector < int > mSignalsPDGs ;
129118 double mHadronRapidityMin ;
130119 double mHadronRapidityMax ;
@@ -140,8 +129,5 @@ FairGenerator*
140129 gen -> addSignalPDGs (443 ); // J/Psi
141130 gen -> addSignalPDGs (100443 ); // Psi(2S)
142131 gen -> addSignalPDGs (11 ); // e
143- // auto seed = (gRandom->TRandom::GetSeed() % 900000000);
144- // gen->readString("Random:setSeed on");
145- // gen->readString("Random:seed " + std::to_string(seed));
146132 return gen ;
147133}
0 commit comments