Changeset 209 in Idarraga

Timestamp:

Jul 14, 2011, 5:49:51 PM (13 years ago)

Author:

benoit

Message:

 

File:

• allpix/src/AllPixFEI3StandardDigitizer.cc (modified) (7 diffs)

allpix/src/AllPixFEI3StandardDigitizer.cc

@@ -43 +43 @@
         //////////////////////////
 
-        biasVoltage=150.0; //[V]
+        biasVoltage=100.0; //[V]
         Temperature = 300.0;
         detectorThickness = gD->GetSensorZ();
-        resistivity=5000;
+        resistivity=9000;
 
         // true = p-type
@@ -150 +150 @@
         dtIni = 0.01*ns;
 
-        precision = 5;
+        precision = 1;
 
 
@@ -166 +166 @@
          */
 
-        FEIX = 3;
+        FEIX = 4;
 
         switch (FEIX) {
@@ -177 +177 @@
                 chipNoise = 300*elec;
                 m_digitIn.thl = 3500*elec;
-                chargeSharingConstant =0.02;
+                chargeSharingConstant =0.0;
                 break;
 
 
         case 4 :
-                MipTOT=5;
-                MipCharge=20000;
+                MipTOT=10;
+                MipCharge=40000;
                 CounterDepth=15;
                 Lv1Unit = 25*ns;
                 chipNoise = 125*elec;
-                m_digitIn.thl = 2000*elec;
+                m_digitIn.thl = 3200*elec;
                 chargeSharingConstant =0.02;
                 break;
@@ -446 +446 @@
                 eHitTotal = eHitTotal*depletedDepth/detectorThickness;
 
+                tempPixel.first  = (*hitsCollection)[itr]->GetPixelNbX();
+                tempPixel.second = (*hitsCollection)[itr]->GetPixelNbY();
+
+                G4double xpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().x();
+                G4double ypos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().y();
+                G4double zpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().z()+detectorThickness/2.;
+                //cout << TString::Format("x,y,z : %5.5f %5.5f %5.5f",xpos/um,ypos/um,zpos/um)<<endl;
+
+
                 for(G4int nQ  = 0 ; nQ < precision ; nQ++) {
 
                 double eHit = double(eHitTotal)/precision;
                 // Hit Info
-
-                tempPixel.first  = (*hitsCollection)[itr]->GetPixelNbX();
-                tempPixel.second = (*hitsCollection)[itr]->GetPixelNbY();
-
-                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();
-                G4double zpos = CLHEP::RandFlat::shoot()*depletedDepth;
+                G4double zpos = CLHEP::RandGauss::shoot(zpos,10*um);
 
 
@@ -620 +621 @@
                         digit->SetPixelCounts(EnergyToTOT((*pCItr).second,threshold));
                         //digit->IncreasePixelCounts(); // Counting mode
+                        cout << "TOT= "<< EnergyToTOT((*pCItr).second,threshold) << endl;
+                        cout << "Energy= "<< (*pCItr).second/elec << endl;
 
                         // MC only //
@@ -748 +751 @@
 {
 
-        G4double dEdt = (MipCharge*elec- threshold)/(MipTOT*Lv1Unit*ns);
-        G4int TOT = TMath::FloorNint((Energy-threshold)/(dEdt*Lv1Unit*ns));
+        G4double dEdt = (MipCharge*elec- threshold)/(MipTOT*Lv1Unit);
+        //G4int TOT = TMath::FloorNint((Energy-threshold)/(dEdt*Lv1Unit*ns));
+        G4int TOT = TMath::FloorNint((Energy-threshold)/(dEdt*Lv1Unit));
+
         if(Energy<threshold) TOT=0;
         if (TOT<0) TOT=0;