Yes the red region is defined by my cut bounds. I looked at the 20 TDCs associated with the PMTs and fit the TDC signal with a Gaussian. I then too the width to be 2*sigma for the Gaussian fit. So I have two arrays of length 20 corresponding to the TDC peak and sigma for each PMT. The difference in the TDC cuts is due to how well the fit fits the TDC signal on each PMT. Some are TDC signals are not as sharp as others. Should I fix them to all have the same width about the TDC peak? <br>
<br><div class="gmail_quote">On Fri, Mar 12, 2010 at 10:55 AM, Brad Sawatzky <<a href="mailto:brads@jlab.org">brads@jlab.org</a>> wrote:<br><blockquote class="gmail_quote">
<div class="im">On Fri, 12 Mar 2010, MATTHEW R POSIK wrote:<br>
<br>
> Yes, that is a good point. Looking at the actual hits there is never<br>
> any above 5 hits for any PMT. But when I do a<br>
> T->MakeSelector("my_selector" ) on the TTree, it outputs a<br>
> my_selector.C and my_selector.h where the .h file lists all variable<br>
> names,types and array sizes. This is what it lists for the DBB.BBcerT*<br>
> stuff:<br>
><br>
> Double_t DBB_BBcerT01[16]; //[Ndata.DBB.BBcerT01]<br>
</div>[ . . . ]<br>
<div class="im">> So there are arrays of length of 16? Am I interpreting this correctly?<br>
<br>
</div>Yup, you're correct.<br>
<br>
[ . . . ]<br>
<div class="im">> I also use this TBranch to set the TDCs to a double array of legth 20.<br>
> Note that if I try to define the length as tdc_values[20][5] or<br>
> something less I get the following error:<br>
><br>
> *** Break *** segmentation violation<br>
</div>[ . . . ]<br>
<div class="im">> Double_t tdc_values[20][20]<br>
> for(pmt=0;pmt<20;pmt++) {<br>
><br>
> TBranch * btdc = chain->GetBranch(Form("DBB.BBcerT%02d",pmt+1) ) ;<br>
> if( btdc ) {<br>
> chain->SetBranchAddress(Form("DBB.BBcerT%02d",pmt+1), &tdc_values[pmt] ,<br>
> &btdc ) ;<br>
<br>
</div>The segfault occurs because you define an array tdc_values[X][Y], and<br>
then ask root to stuff 16 values (the allocated full DBB.BBcerT%02d<br>
array size) into tdc_values[X], which is of length Y. If Y<16 you'll<br>
have problems.<br>
<div class="im"><br>
> I do not think this leads to the same issue as with using the TString<br>
> method, due to the fact that I only loop over the number of hits on<br>
> each PMT for each event:<br>
<br>
> //Loop over PMTs<br>
> for(pmt=0;pmt<20; pmt++) {<br>
> h_hits[pmt]->Fill(tdc_hits[pmt]);<br>
> // Loop over cherenkov tdc hits<br>
> for(j=0;j<tdc_hits[pmt];j++) {<br>
> h_tdc[pmt]->Fill(tdc_values[pmt][j]);<br>
> }<br>
> }<br>
><br>
> So j only goes to the number of hits on that PMT for that event.<br>
> Looking at the TDCs and TDC hits they seem to look good. See plots<br>
> below, red histo is a TDC cut using the Multi-hit TDCs. Does this<br>
> seem correct?<br>
<br>
</div>That approach looks reasonable.<br>
<br>
Why doesn't the red region have the same width for all PMTs though?<br>
Isn't that width just defined by your cut bounds?<br>
<div><div></div><div class="h5"><br>
-- Brad<br>
<br>
--<br>
Brad Sawatzky, PhD <<a href="mailto:brads@jlab.org">brads@jlab.org</a>> -<>- Jefferson Lab / Hall C / C111<br>
Ph: 757-269-5947 -<>- Pager: 757-584-5947 -<>- Fax: 757-269-5235<br>
The most exciting phrase to hear in science, the one that heralds new<br>
discoveries, is not "Eureka!" but "That's funny..." -- Isaac Asimov<br>
_______________________________________________<br>
d2n-analysis-talk mailing list<br>
<a href="mailto:d2n-analysis-talk@jlab.org">d2n-analysis-talk@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/d2n-analysis-talk">https://mailman.jlab.org/mailman/listinfo/d2n-analysis-talk</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Matthew Posik<br>Email: <a href="mailto:posik@temple.edu">posik@temple.edu</a><br>Temple University Physics Dept.<br>Office: BA-319<br>Office #: 215-204-1331 <br>
WebSites:<br>Temple:<br><a href="http://quarks.temple.edu/">http://quarks.temple.edu/</a><br>d2n:<br><a href="http://hallaweb.jlab.org/experiment/E06-014/">http://hallaweb.jlab.org/experiment/E06-014/</a><br><br>