chain2run.C File Reference

Functions
void	chain2run (TChain &c, char *rout="r.root")
bool	GoodEvent (EdbView *v)

Function Documentation

chain2run()

void chain2run	(	TChain &	c,
		char *	rout = "r.root"
	)

{
  // input - chain of files(runs) - output is one run
 
  EdbRun r( rout, "RECREATE" );
  r.GetFile()->SetCompressionLevel(9);
  EdbView *v = r.GetView();
 
  c.SetBranchStatus("clusters",  0  );
  c.SetBranchAddress("frames",   v->GetFramesAddr()   );
  c.SetBranchAddress("headers",  v->GetHeaderAddr()   );
  c.SetBranchAddress("segments", v->GetSegmentsAddr() );
 
  int n=c.GetEntries();
  printf("entries=%d\n",n);
  for(int i=0; i<n; i++) {
    c.GetEntry(i);
 
    if( !GoodEvent(v) ) continue;      // selections
    v->GetHeader()->SetStatus( c.GetTreeNumber() );      // keep memory about each tree (if needed)
    r.AddView(v);
    if( !(i%1000) ) printf("%d (%d\%)\n",i,100*i/n);
  }
  r.Close();
}

GoodEvent()

bool GoodEvent

(

EdbView *

v

)

{
  // do all selections here
 
  if( v->Nsegments() < 0) return 0;
  TObjArray *segs = v->GetSegments();
  int n = v->Nsegments();
  for(int i=0; i<n; i++) {
    s = v->GetSegment(i);
    if( s->GetPuls()<8 )    
      { segs->RemoveAt(i); continue; }
    if( TMath::Sqrt(s->GetTx()*s->GetTx() + s->GetTy()*s->GetTy()) >0.5 ) 
      { segs->RemoveAt(i); continue;}
  }
  segs->Compress();
 
  return true;
}