Here is a little instruction on how to use the RdoDC module.
In your script you should write:
BrDcRdoModule* rdodc = new BrDcRdoModule("T5", "T5");
rdodc->SetCalibFile(Form("calpar/T5CalibRun%d.dat", CalibRunNo));
rdodc->SetOffsetFile(Form("calpar/T5OffsetRun%d.dat", CalibRunNo));
// the calpar directory is the link to PS' directory with calibrations:
//~ufstasze/brahms/brahms_app/ps_app/dc/calib/calpar/
rdodc->SetCalibChan(460);
//it is 460 for T4 and T5, and 290 for T3 ( Pawel, can you correct me if
//I'm wrong )
rdodc->SetVerbose(verboseOption->GetValue());
mainModule->AddModule(rdodc);
//------------------------------------------------------------------------
To read the RdoDC data inside the module is easy:
BrDataTable* RdoDCTable=
outNode->GetDataTable( Form ( "RdoDC %s", GetName() ) );
Int_t NumCleanHits = RdoDCTable->GetEntries();
for(Int_t ih=0;ih<NumCleanHits;ih++)
{
BrDetectorHit* rdohit_p = (BrDetectorHit*)RdoDCTable->At(ih);
Float_t u[2] = {rdohit_p->GetPos()[0], rdohit_p->GetPos()[1] };
//this is the hit distance from centre of the place, there are two
//distances, because of the hit ambiguity
Float_t z = rdohit_p->GetPos()[2];
Int_t im = rdohit_p->GetImod()/100;
Int_t ip = rdohit_p->GetImod()%100;
//im and ip are in fact necessary only to fillin histograms, of course if
//you like it
}
This archive was generated by hypermail 2b30 : Fri Apr 12 2002 - 07:38:41 EDT