EdbViewGen Class Reference

#include <EdbViewGen.h>

Inheritance diagram for EdbViewGen:

Collaboration diagram for EdbViewGen:

Public Member Functions
	EdbViewGen ()
	EdbViewGen (EdbViewDef &vd)
int	GenAlfaGrains ()
int	GenDeltaGrains ()
int	GenFogGrains (int ngr, TObjArray &grains, int side=0)
int	GenFrames (EdbView &v)
int	GenGrainClusters (EdbView &v, EdbCluster &g)
int	GenGrains (EdbView &v)
int	GenSegGrains (EdbSegment &s)
float	GrainPathMip (float lambda=2.)
	~EdbViewGen ()
Public Member Functions inherited from EdbViewDef
	EdbViewDef ()
void	Print ()
void	SetDef ()
	~EdbViewDef ()

Additional Inherited Members
Public Attributes inherited from EdbViewDef
Float_t	eClaSX
Float_t	eClaSY
Float_t	eClaSZ
Float_t	eClaSZvar
Float_t	eDZdead
Float_t	eFogDens
Float_t	eFogGrainArea
Float_t	eGrainArea
Float_t	eGrainSX
Float_t	eGrainSY
Float_t	eGrainSZ
Int_t	eNframes
Float_t	eTXopt
Float_t	eTYopt
Float_t	eX0
Float_t	eX0opt
Float_t	eXmax
Float_t	eXmin
Float_t	eY0
Float_t	eY0opt
Float_t	eYmax
Float_t	eYmin
Float_t	eZdead
Float_t	eZmax
Float_t	eZmin
Float_t	eZxy

Constructor & Destructor Documentation

EdbViewGen() [1/2]

EdbViewGen::EdbViewGen

(

)

EdbViewGen() [2/2]

EdbViewGen::EdbViewGen ( EdbViewDef &  vd )

inline

: EdbViewDef(vd) {}

~EdbViewGen()

EdbViewGen::~EdbViewGen ( )

inline

{}

Member Function Documentation

GenAlfaGrains()

int EdbViewGen::GenAlfaGrains ( )

inline

{return 0;}   // bg from alfa-particles

GenDeltaGrains()

int EdbViewGen::GenDeltaGrains ( )

inline

{return 0;}  // bg from delta-electrons

GenFogGrains()

int EdbViewGen::GenFogGrains	(	int	ngr,
		TObjArray &	grains,
		int	side = 0
	)

{
  float  x,y,z;
  int    area;
  int    n=0;
 
  printf("GenFogGrains: X:(%f %f)  Y:(%f %f) \n", eXmin,eXmax, eYmin,eYmax );
  while( n<ngr ) {
    x      = eXmin + gRandom->Rndm()*(eXmax-eXmin);
    y      = eYmin + gRandom->Rndm()*(eYmax-eYmin);
    z      = eZmin + gRandom->Rndm()*(eZmax-eZmin);
 
    area = (int)(eFogGrainArea); // + gRandom->Poisson(4) - 2);
    grains.Add( new EdbCluster(x,y,z, area, 0, 0, side, -1) );
    n++;
  }
  return n;
 }

GenFrames()

int EdbViewGen::GenFrames

(

EdbView &

v

)

{
   v.SetNframes(eNframes,0);
   float step=(eZmax-eZmin-1.)/(eNframes-1);
   for(int i=0; i<eNframes; i++)
     v.AddFrame( i, eZmin + 0.5 + i*step );
   
   return eNframes;
 }

GenGrainClusters()

int EdbViewGen::GenGrainClusters	(	EdbView &	v,
		EdbCluster &	g
	)

{
  // generate clusters for the input grain g
 
  int ncl=0;
  if( g.GetArea()==0 ) return ncl;
  int nfr= v.GetNframes();
  float z;
  float dz = 0.5 * gRandom->Gaus(eGrainSZ,eGrainSZ/3.);
  //eClaSZ*TMath::Sqrt(g.GetArea())/3;  //ToDo
  float dx,dy;
  EdbCluster *c;
  for(int i=0; i<nfr; i++) {
    z= v.GetFrame(i)->GetZ();
    if( TMath::Abs(z - g.GetZ()) < dz ) {
      c= new EdbCluster(g);
      gRandom->Rannor(dx,dy);
      c->SetZ(z);
      c->SetX( g.X()+ dx*eClaSX );           // TODO: navesti nauku
      c->SetY( g.Y()+ dy*eClaSY );
      c->SetFrame(i);
      v.AddCluster(c);
      ncl++;
    }
  }
  //printf("ncl = %d\n",ncl);
  return ncl;
}

GenGrains()

int EdbViewGen::GenGrains

(

EdbView &

v

)

{
  TObjArray grains;
  int nseg =  v.Nsegments();
  //printf("nsegments = %d\n", nseg);
  EdbSegment   *s = 0;
  int ngrseg=0;
  int ngr;
  for(int iseg=0; iseg<nseg; iseg++) {
    s = v.GetSegment(iseg);
    ngr = GenSegGrains( *s );
    printf("ngr = %d\n", ngr);
    for(int igr=0; igr<ngr; igr++) grains.Add(s->GetElements()->At(igr));
    ngrseg+=ngr;
  }
  //printf("ngrseg = %d\n", ngrseg);
  
  int nfog = (int)(eFogDens*(eXmax-eXmin)*(eYmax-eYmin)*(eZmax-eZmin)/1000.);
  printf("ngrfog = %d\n", nfog);
  GenFogGrains(nfog,grains);
 
  EdbCluster *g;
  ngr = grains.GetEntries();
  printf("%d grains generated\n",ngr);
 
  // remove grains from the dead layer
  int icrem=0;
  for(int i=0; i<ngr; i++) {
    g = (EdbCluster*)grains.At(i);
    if( TMath::Abs(g->GetZ()-eZdead)<eDZdead/2. ) { g->SetArea(0); icrem++;}
  }
  printf("%d grains removed in dead layer\n",icrem);
 
  for(int i=0; i<ngr; i++) {
    g = (EdbCluster*)grains.At(i);
    GenGrainClusters(v,*g);
    //v.AddCluster(new EdbCluster(*g));
  }
  
  for(int iseg=0; iseg<nseg; iseg++) 
    if( v.GetSegment(iseg)->GetElements() )
      v.GetSegment(iseg)->GetElements()->Delete();
 
  return ngr;
}

GenSegGrains()

int EdbViewGen::GenSegGrains

(

EdbSegment &

s

)

{
 // input: segments of the view, assumed that puls of the 
 // segment is correspond to the number of grains to be generated
 // output: array of grains belong to the segment
 //  ToDo: dz only positive?
 
  using namespace TMath;
 
  // could be switched on if serve...
  //
  //   int ngr = s.GetPuls();
  //   TArrayF arr(ngr);
  //   gRandom->RndmArray(ngr, arr.GetArray());
  //   for(int i=0; i<ngr; i++)  printf("%f ",arr[i]);
  //   printf("\n");
  //   TArrayI ind(ngr);
  //   Sort( n, arr.GetArray(), ing.GetArray(), 0 );
 
  float x  = s.GetX0();
  float y  = s.GetY0();
  float z  = s.GetZ0();
  float dz = s.GetDz();
  float zmax=z+dz;
  float  cs = Cos(ATan(Sqrt(s.GetTx()*s.GetTx()+s.GetTy()*s.GetTy())));
  int    area;
  float  step;
  int    ncl=0;
  int    nmax=10000;
 
  float lambda = dz/s.GetPuls();
 
  for(int i=0; i<nmax; i++ ) {
    step = GrainPathMip(lambda)*cs;
    z += step;
    //printf("Zstep: %5.3f \n",step);
    x += s.GetTx()*step;
    y += s.GetTy()*step;
    if( z > zmax  )   break;
    if( x < eXmin )   continue;
    if( x > eXmax )   continue;
    if( y < eYmin )   continue;
    if( y > eYmax )   continue;
    area = (int)(eGrainArea); // + gRandom->Poisson(4) - 2);
   
    double sx,sy,sz, r=gRandom->Gaus(0.,eGrainSX);  // the smearing of the physical grain center position
    gRandom->Sphere(sx,sy,sz,r);
    x+=sx;    y+=sy;    z+=sz;
    if(x<eXmin) continue;
    if(x>eXmax) continue;
    if(y<eYmin) continue;
    if(y<eYmin) continue;
    s.AddElement( new EdbCluster( x+sx, y+sy, z+sz, area, 0, 0, s.GetSide(), s.GetID()) );
    ncl++;
  }
  s.SetPuls(ncl);
  return ncl;
 }

GrainPathMip()

float EdbViewGen::GrainPathMip

(

float

lambda = 2.

)

{
  //=========================================
  //
  // Some emulsion matter properties:
  // density: 2.4 g/cm^3
  //    <A>  = 18.2
  //    <Z>  = 8.9
  //     X0  = 5.5 cm
  // nuclear collision lenght: 33 cm
  //
  // <dE/dX> = 37 KeV/100 microns
  //
  // visible grain density: 30/100 microns
  // grain size: 0.6 micron
  //
  // assuming the ionization energy about 4 eV (to check!)
  // we expect about 100 ionizations/micron
  //
  //=========================================
  
  // lambda    // mean path length
 
  int    nstep=100;  
  double prob0 = 1./nstep;   //prob of grain after single ionization
  float  grsiz=gRandom->Uniform(0.9,1.1);
  float  x=0;
  for(int j=0; j<100000; j++) {
    x+=gRandom->Exp(lambda/nstep);
    if( x>grsiz&&gRandom->Rndm()<prob0 ) {
      break;
    }
  }
  return x;
}

---

The documentation for this class was generated from the following files:

• /home/antonio/fedra_doxygen/src/libEMC/EdbViewGen.h
• /home/antonio/fedra_doxygen/src/libEMC/EdbViewGen.cxx