Message boards : CMS Application : Could not decay:pi0 with mass:0 will throw event away!
Message board moderation

To post messages, you must log in.

AuthorMessage
Crystal Pellet
Volunteer moderator
Volunteer tester

Send message
Joined: 14 Jan 10
Posts: 1268
Credit: 8,421,616
RAC: 2,139
Message 28081 - Posted: 7 Dec 2016, 10:49:39 UTC

After 35000 events all events throwing away. Job will stop after 10000 rejections.

From batch 161205_162337:ireid_crab_BPH-RunIISummer15GS-00046_BD job 7228
ID: 28081 · Report as offensive     Reply Quote
maeax

Send message
Joined: 2 May 07
Posts: 2071
Credit: 156,091,089
RAC: 103,567
Message 28082 - Posted: 7 Dec 2016, 16:44:12 UTC

ID: 28082 · Report as offensive     Reply Quote
Crystal Pellet
Volunteer moderator
Volunteer tester

Send message
Joined: 14 Jan 10
Posts: 1268
Credit: 8,421,616
RAC: 2,139
Message 28083 - Posted: 7 Dec 2016, 17:21:20 UTC

Maybe the 10,000 rejections can safely be reduced to e.g. 1,000, so the VM can start faster with the next job.
I also had such a case after only 3,000 events from 200,000 events.
The VM was 2 hours useless consuming CPU, before the job was aborted by itself.
ID: 28083 · Report as offensive     Reply Quote
ivan
Volunteer moderator
Project tester
Volunteer developer
Volunteer tester
Project scientist
Avatar

Send message
Joined: 29 Aug 05
Posts: 997
Credit: 6,264,307
RAC: 71
Message 28122 - Posted: 15 Dec 2016, 0:12:47 UTC - in response to Message 28083.  
Last modified: 15 Dec 2016, 0:22:36 UTC

Maybe the 10,000 rejections can safely be reduced to e.g. 1,000, so the VM can start faster with the next job.
I also had such a case after only 3,000 events from 200,000 events.
The VM was 2 hours useless consuming CPU, before the job was aborted by itself.

OK, I need to find out if that's a configurable parameter. Please remind me if I don't respond within a couple of days... :-)
ID: 28122 · Report as offensive     Reply Quote
ivan
Volunteer moderator
Project tester
Volunteer developer
Volunteer tester
Project scientist
Avatar

Send message
Joined: 29 Aug 05
Posts: 997
Credit: 6,264,307
RAC: 71
Message 28131 - Posted: 16 Dec 2016, 11:37:45 UTC - in response to Message 28122.  
Last modified: 16 Dec 2016, 11:43:12 UTC

Maybe the 10,000 rejections can safely be reduced to e.g. 1,000, so the VM can start faster with the next job.
I also had such a case after only 3,000 events from 200,000 events.
The VM was 2 hours useless consuming CPU, before the job was aborted by itself.

OK, I need to find out if that's a configurable parameter. Please remind me if I don't respond within a couple of days... :-)

After a bit of gwgling I found some code at http://svn.cern.ch/guest/evtgen/tags/R00-04-00/src/EvtGen.cpp:
int times=0;
do {
times+=1;
EvtStatus::initRejectFlag();
p->decay();
//ok then finish.
if ( EvtStatus::getRejectFlag()==0 ) {
times=0;
}

...
if ( times==10000) {
report(ERROR,"EvtGen") << "Your event has been rejected 10000 times! << endl;
report(ERROR,"EvtGen") << "Will now abort."<< endl;
::abort();
times=0;
}
} while (times);


So, if that is indeed the code we are using, then there doesn't seem to be an easy way of changing the limit.

[Edit] The exact same code is in the latest release :-( [/Edit]
ID: 28131 · Report as offensive     Reply Quote
Crystal Pellet
Volunteer moderator
Volunteer tester

Send message
Joined: 14 Jan 10
Posts: 1268
Credit: 8,421,616
RAC: 2,139
Message 28132 - Posted: 16 Dec 2016, 13:12:14 UTC - in response to Message 28131.  


...
if ( times==10000) {
report(ERROR,"EvtGen") << "Your event has been rejected 10000 times! << endl;
report(ERROR,"EvtGen") << "Will now abort."<< endl;
::abort();
times=0;
}
} while (times);


So, if that is indeed the code we are using, then there doesn't seem to be an easy way of changing the limit.

Yes, that's exactly the text, I've seen in the Running.log.
The easy way would be skipping 1 zero (in statement and text to report), but of course I've no idea what's the impact elsewhere.
The programmers maybe can tell why 10,000 times is set. I suppose it was rather random chosen. At least a lower value would decrease (almost) idle CPU-time.

btw: yesterday I had again such a task on the dev-project. Thrown away 10,000 times after about the 35,000th event out of 200,000.
ID: 28132 · Report as offensive     Reply Quote
ivan
Volunteer moderator
Project tester
Volunteer developer
Volunteer tester
Project scientist
Avatar

Send message
Joined: 29 Aug 05
Posts: 997
Credit: 6,264,307
RAC: 71
Message 28133 - Posted: 16 Dec 2016, 15:19:19 UTC - in response to Message 28132.  

I found 135 in one batch, so there's a bit over 1%. It's worth noting that this is in EvtGen (EventGenerator) which is separate from Pythia (but obviously used by it). It's included in CMSSW as an external, and as you can see from the code, there's no way we can make a change to it. Short of convincing the developers to make the change...
ID: 28133 · Report as offensive     Reply Quote

Message boards : CMS Application : Could not decay:pi0 with mass:0 will throw event away!


©2024 CERN