Hi,
I discovered that if you play too much with histogram names you will get
very bad behaviours. For instance
if(l == 2 && i == 0 && thebin == 1){
sprintf(name, "%s%d", "a", 1609);
((TH1D*)gDirectory->Get(name))->Fill(momentumTrk[itrk]);
sprintf(name,"%s%d%s%d","a",1600,"bin",thebin);
((TH1D*)gDirectory->Get(name))->Fill(momentumTrk[itrk]);
}
the two histograms have same binning and same range.
Well, I get two different numbers of entries. The first plot has the
correct number, the second one has an excess. This histogram is not filled
elsewhere. If I comment out the filling I get 0 entries, then no other
filling elsewhere. The extra events have shape proper to a mes
distribution. Other histogram with a similar name are filled with mes.
Very suspicious...
If I change the histogram name things work properly again (same
distribution for both the histograms).
Millions of histograms and tree variables are filled in recoilNtp in this
way. I don't know if I will sleep well in the next days.
Daniele
|