408{
409 int i;
412
413 {
417 TClonesArray *dc = new TClonesArray("Deriv");
419 float zmin,zmax;
420
422 int nFrame =
t->GetEntries();
423
424 for (int ifr=0; ifr<nFrame; ifr++) {
425
427 if (ifr==0) fpred=
NULL;
429 if (ifr==nFrame-1) fsuc=
NULL;
431
432
433 d->derivate(fpred,
f,fsuc);
434
435
436
437
438
439 new((*dc)[ifr])
Deriv(*
d);
440
441 }
444 if (v->GetNframesTop()) {
447 printf(
"!!!!!!!! dztop=%f, z1=%f, z2=%f\n",
Z1-
Z2,
Z1,
Z2);
450 }
451 }
452 else {
455 printf(
"!!!!!!!! dzbot=%f, z3=%f, z4=%f\n",
Z3-
Z4,
Z3,
Z4);
458 }
459 }
460
461
462 dc->Delete();
465 delete dc;
467 }
469 TClonesArray* frames = v->GetFrames();
470
473 int nframes = v->GetNframes();
474 int nclusters = v->Nclusters();
475
477 float zmax = zmin;
478 for (i=0; i<nframes; i=i+rep) {
480 if (inside)
483 if (zmin>
f->GetZ()) zmin =
f->GetZ();
484 if (zmax<f->GetZ()) zmax =
f->GetZ();
487
488 }
490 for (i=0; i<nclusters; i++) {
491 int j;
492 int flag = false;
493 c = v->GetCluster(i);
494 float xc = c->
GetX();
495 float yc = c->
GetY();
496 if (
In(xc,yc,Xmin,Xmax,Ymin,Ymax)) {
498
500 flag = true;
501 break;
502 }
503 }
504 if (!flag) continue;
509 }
510 }
512
513
514
515
516
517
518
519
520
522
523
524
525}
void d()
Definition: RecDispEX.C:381
FILE * f
Definition: RecDispMC.C:150
int In(float X, float Y, float Xmin, float Xmax, float Ymin, float Ymax)
Definition: bitview.cxx:18
float min(TClonesArray *t)
Definition: bitview.cxx:275
void PreciseBasePosition(float zmin, float zmax, TClonesArray *dc, float &zup, float &zdown)
Definition: bitview.cxx:315
TClonesArray * Smoothing(EdbView *v)
Definition: bitview.cxx:237
float max(TClonesArray *t)
Definition: bitview.cxx:297
Definition: EdbCluster.h:19
Float_t GetX() const
Definition: EdbCluster.h:51
Float_t GetY() const
Definition: EdbCluster.h:52
Float_t GetArea() const
Definition: EdbCluster.h:54
Float_t GetZ() const
Definition: EdbCluster.h:53
Definition: EdbFrame.h:20
Definition: bitmatrix.h:13
float GetZ(void)
Definition: bitmatrix.h:37
void SetBit(long i, long j)
Definition: bitmatrix.cxx:120
float Z3
Definition: bitview.h:31
TBitMatrix * GetLayer(int n)
Definition: bitview.cxx:401
int GetNumberOfLayers(void)
Definition: bitview.h:47
float Z4
Definition: bitview.h:31
float Z1
Definition: bitview.h:31
float Z2
Definition: bitview.h:31
TTree * t
Definition: check_shower.C:4
MFTYPE32 long(MFTYPE MPTYPE *MOCRHOOKFCTPTR)(long HookType
Definition: milocr.h:32
#define NULL
Definition: nidaqmx.h:84