@@ -122,16 +122,16 @@ protected:
122122 }
123123
124124 bool findSignalInAcceptance () {
125- std ::cout <<"loop over" << mGeneratorParam -> getParticles ().size ()<<" particles" <<std ::endl ;
125+ // std::cout<<"loop over" << mGeneratorParam->getParticles().size()<<" particles"<<std::endl;
126126 for (int pdg : mSignalsPDGs ) {
127- std ::cout <<"signal pdg: " <<pdg <<std ::endl ;
127+ // std::cout<<"signal pdg: "<<pdg<<std::endl;
128128 }
129129 for (int ipart = 0 ; ipart < mGeneratorParam -> getParticles ().size (); ipart ++ ) {
130130 TParticle part = TParticle (mGeneratorParam -> getParticles ().at (ipart ));
131131 // make sure all signals are in the acceptance
132132 for (int pdg : mSignalsPDGs ) {
133133 if (part .GetPdgCode () == pdg ) {
134- std ::cout <<"found signal with pdg: " <<part .GetPdgCode ()<<", mother: " <<part .GetFirstMother ()<<std ::endl ;
134+ // std::cout<<"found signal with pdg: "<<part.GetPdgCode()<<", mother: "<<part.GetFirstMother()<<std::endl;
135135 if (part .GetFirstMother () == -1 ) {
136136 if (part .Y () < mHadronRapidityMin || part .Y () > mHadronRapidityMax ) {
137137 return false;
@@ -146,19 +146,19 @@ protected:
146146
147147 int getNumSignalsForImpactParameter (double impactParameter ) const
148148 {
149- const double photoproductionImpactParameterCut = 3.5 ; // tmp cut to embed signal > 10% centrality
150- if (impactParameter < photoproductionImpactParameterCut ) {
151- std ::cout <<"impact parameter: " <<impactParameter <<" fm, below photoproduction cut "
152- <<photoproductionImpactParameterCut <<" fm; embedding no signal events" <<std ::endl ;
153- return 0 ;
154- }
149+ // const double photoproductionImpactParameterCut = 3.5; // tmp cut to embed signal > 10% centrality
150+ // if (impactParameter < photoproductionImpactParameterCut) {
151+ // std::cout<<"impact parameter: "<<impactParameter<<" fm, below photoproduction cut "
152+ // <<photoproductionImpactParameterCut<<" fm; embedding no signal events"<<std::endl;
153+ // return 0;
154+ // }
155155 return std ::max (0 , static_cast < int > (std ::lround (5.0 + 0.886202881 * std ::pow (std ::max (0.0 , 17.5 - impactParameter ), 1.7 ))));
156156 }
157157
158158 void appendCurrentSignalParticles (int signalIndex )
159159 {
160160 int originalSize = mParticles .size ();
161- std ::cout <<"adding " <<mGeneratorParam -> getParticles ().size ()<<" particles from signal event " <<signalIndex <<" to the stack" <<std ::endl ;
161+ // std::cout<<"adding "<<mGeneratorParam->getParticles().size()<<" particles from signal event "<<signalIndex<<" to the stack"<<std::endl;
162162 for (size_t ipart = 0 ; ipart < mGeneratorParam -> getParticles ().size (); ipart ++ ){
163163 TParticle part = TParticle (mGeneratorParam -> getParticles ().at (ipart ));
164164 if (part .GetFirstMother () >= 0 ) part .SetFirstMother (part .GetFirstMother () + originalSize );
0 commit comments