hi vitaliy,
yes, basically, you do a Setup() to create all of your blastlib2
objects, and then insert the following commands at taste into your code
to spice it up:
// check for new datafile and open it
if (orun!=nrun) {
for (int il=0;il<<gOpt->nList();il++) {
if (run_no(il)==nrun) {
gRecon->clear();
gRecon->Open(il);
orun=nrun;
break;
}
}
if (orun!=nrun)
cerr<<"error: missing coda file for run "<<nrun<<endl;
}
// read specific event from coda file
{
gRaw->Goto(nev);
if (gRecon->GetNEvent()!=nev)
cerr<<"error: missing event "<<nev<<" in run "<<nrun<<endl;
}
// you can then use anything found in TBLRaw for your cuts, etc
{
gRaw->WClrhits[0]; //left sector hits
}
give your script the same options that you would have used with lrn,
ie. "+SQL" or whatever. i haven't actually tried this code, but what
could possibly go wrong?!
--chris
vitaliy ziskin wrote:
> Chris,
> Thanks for your help. Here is what I do.
> 1)Chain up flr files
> 2) Go through all entries and find possible neutron candidates.
> At this point I have an event number, runnumber and would like to peak
> into the Raw for this event to get raw wire chamber info. Is this doable?
>
> Vitaliy
>
> Chris Crawford wrote:
>
>> hi vitaliy,
>> what is the code you use to filter DST/flr/whatever into your own
>> ntuple. if you give me a copy of it, i'll sketch how to analyze the
>> CODA file and recover the lost info. i'm guessing is should be <20
>> extra lines of code to do what you want.
>> --chris
>
>
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:31 EST