Changeset 178 in Idarraga for allpix

Timestamp:

Apr 24, 2011, 3:40:54 PM (13 years ago)

Author:

benoit

Message:

Full Field addition, with 1D field

Location:

allpix

Files:

• include/AllPixFEI3StandardDigitizer.hh (modified) (3 diffs)
• src/AllPixFEI3StandardDigitizer.cc (modified) (20 diffs)

allpix/include/AllPixFEI3StandardDigitizer.hh

@@ -56 +56 @@
   G4double SlimEdgeEffect(G4int nX,G4double xpos,G4double eHit);
   G4bool isSlimEdge(G4int nX, G4int nY);
+  void Efield1D(G4double z);
   G4double elec;
 
@@ -70 +71 @@
   G4double b;
   G4double c;
+  G4double Neff;
 
   G4double detectorThickness;
@@ -114 +116 @@
   G4bool doSlimEdge;
 
+  G4double epsilon;
+  G4double echarge;
+  G4int precision;
+
 
   //////////////////////////////////////////////////////

allpix/src/AllPixFEI3StandardDigitizer.cc

@@ -13 +13 @@
 #include "TF1.h"
 #include "CLHEP/Random/RandGauss.h"
+#include "CLHEP/Random/RandFlat.h"
 
 //#include "SystemOfUnits.h"
@@ -42 +43 @@
         //////////////////////////
 
-        biasVoltage=600.0; //[V]
+        biasVoltage=1000.0; //[V]
         Temperature = 300.0;
         detectorThickness = gD->GetSensorZ();
@@ -98 +99 @@
         //////////////////////
 
-        fluence = 1e15/cm2;
+        fluence = 0.0/cm2;
+        //Beta_electrons = 5e-16*cm2/ns;
         Beta_electrons = 5e-16*cm2/ns;
         b=7.94e-3/cm;
@@ -104 +106 @@
         trappingTime = 1.0/(Beta_electrons*fluence);
 
-        double echarge = 1.60217646e-19;
-        double epsilon = 11.8*8.854187817e-12/m;
-
-        double Neff=0;
+         echarge = 1.60217646e-19;
+         epsilon = 11.8*8.854187817e-12/m;
+
+        Neff=0;
         double Neff0=0;
         if(bulkType) Neff0=-1.0/(resistivity*Default_Hole_Mobility*(s/cm2)*echarge*cm3);
         else  Neff0=1.0/(resistivity*Default_Electron_Mobility*(s/cm2)*echarge*cm3);
 
-        if(!bulkType)Neff=Neff*TMath::Exp(-c*fluence)-b*fluence;
+        if(!bulkType)Neff=Neff0*TMath::Exp(-c*fluence)-b*fluence;
         else Neff=Neff0-b*fluence;
 
@@ -118 +120 @@
         depletedDepth=TMath::Sqrt(2*epsilon*biasVoltage/(echarge*TMath::Abs(Neff)));
         if(depletedDepth>detectorThickness)depletedDepth=detectorThickness;
-        electricFieldY = biasVoltage/depletedDepth; // V/um
+        electricFieldZ = biasVoltage/depletedDepth; // V/um
         electricFieldX = 0; // V/um
-        electricFieldZ = 0; // V/um
+        electricFieldY = 0; // V/um
 
         G4cout << "!!!!!!!Radiation Damage Report !!!!!!!!" << endl
@@ -137 +139 @@
 
         doTrapping =true;
-        doFullField = false;
+        doFullField = true;
 
 
@@ -148 +150 @@
         dtIni = 0.01*ns;
 
+        precision = 5;
+
 
         /////////////////////////////
@@ -168 +172 @@
         case 3 :
                 MipTOT=60;
-                MipCharge=21000;
+                MipCharge=22000;
                 CounterDepth=255;
                 Lv1Unit = 25*ns;
@@ -229 +233 @@
 
         case 1 :
-                doSlimEdge=true;
+                doSlimEdge=false;
                 GRShift =250*um;
                 break;
@@ -294 +298 @@
         G4double ytemp=y;
         G4double ztemp=z;
-
+        //cout << "z : " << z/um << endl;
         vector<G4double> step;
-        while(z>0 && z<detectorThickness){
+        int counter=0;
+        while( ztemp>0 && ztemp<detectorThickness){
 
                 driftTime+=dt;
-                step=RKF5Integration(x,y,z,dt);
+                step=RKF5Integration(xtemp,ytemp,ztemp,dt);
                 xtemp+=step[0];
                 ytemp+=step[1];
                 ztemp+=step[2];
-                SetDt(dt,step[3]);
+                if(GetElectricFieldNorm(xtemp,ytemp,ztemp)==0)ztemp=detectorThickness;
+                dt=SetDt(dt,step[3]);
+                counter++;
+                //cout << "!!!!!! Drifting youhou! !!!!! " << TString::Format("dt: %f z : %f drift : %f Ez : %f",dt/ns,ztemp/um,driftTime/ns,electricFieldZ*cm) << endl;
         }
 
@@ -316 +324 @@
 
 
-void AllPixFEI3StandardDigitizer::ComputeElectricField(G4double /*x*/, G4double /*y*/, G4double /*z*/){}
+void AllPixFEI3StandardDigitizer::ComputeElectricField(G4double /*x*/, G4double /*y*/, G4double z){
+
+                Efield1D(z);
+}
 
 
@@ -327 +338 @@
 
         //ComputeElectricField(x,y,z);
-        G4double parElectricField = electricFieldY;//GetElectricFieldNorm(x,y,z);
+        G4double parElectricField = electricFieldZ;//GetElectricFieldNorm(x,y,z);
         G4double mobilite = Default_Electron_Mobility*
         TMath::Power((1.0/(1.+ TMath::Power((Default_Electron_Mobility*parElectricField)/
@@ -355 +366 @@
         //G4cout << "Hit position with respect to pixel (um) "<< (a_hit->GetPosWithRespectToPixel().z() + detectorThickness/2.0)/um  <<endl;
         //return (a_hit->GetPosWithRespectToPixel().z() + detectorThickness/2.0 ) /(mobility*electricFieldY);
-        double drift = ( CLHEP::RandGauss::shoot((depletedDepth/detectorThickness)*detectorThickness/2,(depletedDepth/detectorThickness)*detectorThickness/4)) /(mobility*electricFieldY);
+        double drift = ( CLHEP::RandGauss::shoot((depletedDepth/detectorThickness)*detectorThickness/2,(depletedDepth/detectorThickness)*detectorThickness/4)) /(mobility*electricFieldZ);
         //cout << "drift: " << mobility*s/cm2 << endl;
         return drift;
@@ -384 +395 @@
 }
 
+
+void AllPixFEI3StandardDigitizer::Efield1D(G4double z){
+
+
+        electricFieldZ=(detectorThickness-z)*biasVoltage/depletedDepth;
+        if(z>depletedDepth){
+                electricFieldZ=0;
+                //cout << "bip" << endl;
+        }
+
+}
+
+
+
+
 /**
  * Digitizer for FEI3Standard
@@ -415 +441 @@
         for(G4int itr  = 0 ; itr < nEntries ; itr++) {
 
+                G4double eHitTotal = elec*CLHEP::RandGauss::shoot((*hitsCollection)[itr]->GetEdep()/elec,TMath::Sqrt((*hitsCollection)[itr]->GetEdep()/elec)*0.118);
+
+                //under-depletion
+                eHitTotal = eHitTotal*depletedDepth/detectorThickness;
+
+                for(G4int nQ  = 0 ; nQ < precision ; nQ++) {
+
+                double eHit = double(eHitTotal)/precision;
                 // Hit Info
 
@@ -427 +461 @@
 
                 G4double ypos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().y();
-                G4double zpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().z();
-
-                G4double eHit = elec*CLHEP::RandGauss::shoot((*hitsCollection)[itr]->GetEdep()/elec,TMath::Sqrt((*hitsCollection)[itr]->GetEdep()/elec)*0.118);
-
-                //under-depletion
-                eHit = eHit*depletedDepth/detectorThickness;
+                //G4double zpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().z();
+                G4double zpos = CLHEP::RandFlat::shoot()*depletedDepth;
+
 
 
@@ -448 +479 @@
                 }
                 else{
+
+                        // Until we get out position
                         vector<G4double> data = ComputeDriftTimeFullField(xpos,ypos,zpos);
                         driftTime = data[3];
@@ -569 +602 @@
         }
 
+
         // Now create digits.  One per pixel
         map<pair<G4int, G4int>, G4double >::iterator pCItr = pixelsContent.begin();
@@ -599 +633 @@
         }
 
+        }
+
         G4int dc_entries = m_digitsCollection->entries();
         if(dc_entries > 0){
@@ -693 +729 @@
 
 //______________________________________________________________________________
-G4double  AllPixFEI3StandardDigitizer::SetDt(G4double Dt, G4double ErreurMoy)
-{
-
+G4double  AllPixFEI3StandardDigitizer::SetDt(G4double dt, G4double ErreurMoy)
+{
+
+        double Dt=dt;
         if(isnan(ErreurMoy)){Dt=tup;}
         else if(ErreurMoy > Target){ Dt*=0.9;}