Changeset 177 in Idarraga for allpix

Timestamp:

Apr 24, 2011, 10:16:36 AM (13 years ago)

Author:

benoit

Message:

SLim edges and Radiation damage

Location:

allpix

Files:

• include/AllPixFEI3StandardDigitizer.hh (modified) (1 diff)
• src/AllPixFEI3StandardDigitizer.cc (modified) (17 diffs)

allpix/include/AllPixFEI3StandardDigitizer.hh

@@ -64 +64 @@
   G4double electricFieldY;
   G4double electricFieldZ;
+  G4double depletionVoltage;
+  G4double depletedDepth;
+  G4double resistivity;
+  G4bool   bulkType;
+  G4double b;
+  G4double c;
 
   G4double detectorThickness;

allpix/src/AllPixFEI3StandardDigitizer.cc

@@ -42 +42 @@
         //////////////////////////
 
-        biasVoltage=150.0; //[V]
+        biasVoltage=600.0; //[V]
         Temperature = 300.0;
         detectorThickness = gD->GetSensorZ();
-        electricFieldY = biasVoltage/detectorThickness; // V/um
-        electricFieldX = 0; // V/um
-        electricFieldZ = 0; // V/um
+        resistivity=5000;
+
+        // true = p-type
+        // false = n-type
+        bulkType =false;
+
 
         ////////////////////////////////
@@ -57 +60 @@
         nPixY= gD->GetNPixelsY();
 
-        //////////////////////
-        // Radiation damage //
-        //////////////////////
-
-        fluence = 5e1/cm2;
-        Beta_electrons = 5e-16*cm2/ns;
-        trappingTime = 1.0/(Beta_electrons*fluence);
+        //cout << TString::Format("Lx=%f Ly=%f npixX=%d npixY=%d",pitchX,pitchY,nPixX,nPixY) << endl;
+
 
 
@@ -91 +89 @@
 
 
-        mobility = MobilityElectron(0,0,0);
+
+
+
+
+
+        //////////////////////
+        // Radiation damage //
+        //////////////////////
+
+        fluence = 1e15/cm2;
+        Beta_electrons = 5e-16*cm2/ns;
+        b=7.94e-3/cm;
+        c=3.54e-13*cm2;
+        trappingTime = 1.0/(Beta_electrons*fluence);
+
+        double echarge = 1.60217646e-19;
+        double epsilon = 11.8*8.854187817e-12/m;
+
+        double 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;
+        else Neff=Neff0-b*fluence;
+
+        depletionVoltage=echarge*TMath::Abs(Neff)*detectorThickness*detectorThickness/(2*epsilon);
+        depletedDepth=TMath::Sqrt(2*epsilon*biasVoltage/(echarge*TMath::Abs(Neff)));
+        if(depletedDepth>detectorThickness)depletedDepth=detectorThickness;
+        electricFieldY = biasVoltage/depletedDepth; // V/um
+        electricFieldX = 0; // V/um
+        electricFieldZ = 0; // V/um
+
+        G4cout << "!!!!!!!Radiation Damage Report !!!!!!!!" << endl
+                 << TString::Format("depletionVoltage : %f depleted Depth : %f",depletionVoltage,depletedDepth/um) << endl
+                 << TString::Format("Neff : %e Neff0 : %e ",Neff*cm3,Neff0*cm3);
+        if(bulkType==false && Neff<0) cout << "Bulk inverted" << endl;
+
+
+
+
 
 
@@ -98 +136 @@
         //////////////////////
 
-        doTrapping =false;
+        doTrapping =true;
         doFullField = false;
 
@@ -130 +168 @@
         case 3 :
                 MipTOT=60;
-                MipCharge=20000;
+                MipCharge=21000;
                 CounterDepth=255;
                 Lv1Unit = 25*ns;
                 chipNoise = 300*elec;
-                m_digitIn.thl = 2500*elec;
+                m_digitIn.thl = 3500*elec;
+                chargeSharingConstant =0.02;
+                break;
+
 
         case 4 :
@@ -144 +185 @@
                 m_digitIn.thl = 2000*elec;
                 chargeSharingConstant =0.02;
+                break;
 
 
@@ -154 +196 @@
                 m_digitIn.thl = 800*elec;
                 chargeSharingConstant =0.01;
+                break;
 
 
@@ -165 +208 @@
                 chargeSharingConstant =0.01;
 
-
-
         }
 
@@ -188 +229 @@
 
         case 1 :
-                doSlimEdge=false;
+                doSlimEdge=true;
                 GRShift =250*um;
+                break;
         case 2 :
                 doSlimEdge=true;
                 GRShift =250*um;
+                break;
         case 3 :
                 doSlimEdge=true;
                 GRShift =500*um;
+                break;
         case 4 :
                 doSlimEdge=true;
                 GRShift =500*um;
+                break;
         case 5 :
                 doSlimEdge=true;
                 GRShift =250*um;
+                break;
 
         default :
@@ -207 +253 @@
                 GRShift =0*um;
 
+
+                //fix me !!
+
         }
+
+                mobility = MobilityElectron(0,0,0);
 
 
@@ -281 +332 @@
                         Electron_Saturation_Velocity,Electron_Beta))),1.0/Electron_Beta);
         return mobilite;
-        //cout << mobilite << " " << parElectricField << endl;
+        //cout << "!!!!!!!!!" << mobilite << " " << parElectricField*cm << endl;
 
 }
@@ -304 +355 @@
         //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);
-        return ( detectorThickness) /(mobility*electricFieldY);
-
+        double drift = ( CLHEP::RandGauss::shoot((depletedDepth/detectorThickness)*detectorThickness/2,(depletedDepth/detectorThickness)*detectorThickness/4)) /(mobility*electricFieldY);
+        //cout << "drift: " << mobility*s/cm2 << endl;
+        return drift;
 }
 
@@ -369 +421 @@
 
                 G4double xpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().x();
+
+                //Ugly Hack !!!
+                if(xpos>0)xpos=-pitchX/2+xpos;
+                else xpos=pitchX/2+xpos;
+
                 G4double ypos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().y();
                 G4double zpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().z();
@@ -374 +431 @@
                 G4double eHit = elec*CLHEP::RandGauss::shoot((*hitsCollection)[itr]->GetEdep()/elec,TMath::Sqrt((*hitsCollection)[itr]->GetEdep()/elec)*0.118);
 
-
-
-                if(doSlimEdge){
+                //under-depletion
+                eHit = eHit*depletedDepth/detectorThickness;
+
+
+                //if(doSlimEdge){
                         if (isSlimEdge(tempPixel.first,tempPixel.second)) eHit = SlimEdgeEffect(tempPixel.first,xpos,eHit);
-                }
+                //}
 
 
@@ -386 +445 @@
                          driftTime = ComputeDriftTimeUniformField((*hitsCollection)[itr]);
                          sigma = ComputeDiffusionRMS(driftTime);
+                         //cout << TString::Format("!!!!!!!!! vd/vdep : %f drift time : %f sigma : %f",depletedDepth/detectorThickness,driftTime,sigma) << endl;
                 }
                 else{
@@ -403 +463 @@
                 //G4cout << TString::Format("[MC Truth] Hit Energy=%f, x=%f, y=%f pixel(%i,%i)",eHitTruth/elec,xpos/um,ypos/um,tempPixel.first,tempPixel.second) << endl;
                 //if(fabs(xpos)>=pitchX/2.-10*sigma && fabs(ypos)>pitchY/2.0-10*sigma){
-                if( (fabs(xpos) >= pitchX/2.-3*sigma || fabs(ypos) >=pitchY/2.-3*sigma) ){
+
+                if( (fabs(xpos) >= pitchX/2.-6*sigma || fabs(ypos) >=pitchY/2.-6*sigma) ){
 
 //              if( (fabs(xpos) >= pitchX/2.-3*sigma || fabs(ypos) >=pitchY/2.-3*sigma) ){
@@ -668 +729 @@
         G4double Etemp;
 
+        double xpix= xpos + pitchX/2;
+        //cout << "!!!!!!!!!SLIMEDGES!!!!!!!!!!! " << MipCharge << " " << xpix << endl;
+
+
         if(nX==0)
         {
-                if((xpos+pitchX/2. < GRShift))eEdge = eHit*(1-0.5*(xpos+pitchX/2)/(100*um));
+                if((xpix < GRShift)){
+                        eEdge = eHit*(1-(GRShift-xpix)/(300*um));
+                        //cout << TString::Format("!!!!!!SLIM EDGES!!!!!!!! left nX:%d xpos:%f xpix:%f before:%f after:%f, R:%f",nX,xpos,xpix,eHit,Etemp,eHit*(1-(GRShift-xpix)/(GRShift))) << endl;
+
+                }
                 else{eEdge=eHit;};
                 Etemp=eEdge;
+
                 }
 
         else if(nX==(nPixX-1)){
-                if((pitchX/2.-xpos < GRShift))eEdge = eHit*(1-0.5*(pitchX/2.-xpos)/(100*um));
+                if(xpix > pitchX-GRShift){
+                        eEdge = eHit*(1-(xpix-(pitchX-GRShift))/(300*um));
+                        //cout << TString::Format("!!!!!!SLIM EDGES!!!!!!!! right NX:%d xpos:%f xpix:%f before:%f after:%f R:%f",nX,xpos,xpix,eHit,Etemp,eHit*(1-(xpix-(pitchX-GRShift))/(GRShift))) << endl;
+
+                }
                 else{eEdge=eHit;};
                 Etemp=eEdge;
 
+
         }