Changeset 232 in Idarraga

Timestamp:

Sep 5, 2011, 5:23:25 PM (13 years ago)

Author:

idarraga

Message:

 

File:

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

allpix/src/AllPixFEI3StandardDigitizer.cc

@@ -188 +188 @@
                 MipTOT=10;
                 MipCharge=40000;
+                //MipCharge=20000;
                 CounterDepth=15;
-                Lv1Unit = 25;
+                Lv1Unit = 25*ns;
                 chipNoise = 125*elec;
                 m_digitIn.thl = 3200*elec;
@@ -443 +444 @@
         //TF1 * f1 = new TF1("expint","exp((x-[0])*(x-[0])/2.*[1])",-100.,100.);
 
+        G4double hitsETotal = 0.;
+
         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);
+                hitsETotal += eHitTotal;
 
                 //under-depletion
@@ -453 +457 @@
                 tempPixel.second = (*hitsCollection)[itr]->GetPixelNbY();
 
-                G4cout << "x : " << tempPixel.first << " ,  y : " << tempPixel.second << G4endl;
+                G4cout << "x : " << tempPixel.first << " ,  y : " << tempPixel.second << ", E = " << eHitTotal/keV << G4endl;
+
 
                 G4double xpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().x();
@@ -566 +571 @@
 
         }
+
+        G4cout << "total = " << hitsETotal/keV << " keV" << G4endl;
+
         // Now create digits.  One per pixel
         map<pair<G4int, G4int>, G4double >::iterator pCItr = pixelsContent.begin();
@@ -573 +581 @@
                 // If the charge in a given pixel is over the threshold
                 double threshold =CLHEP::RandGauss::shoot(m_digitIn.thl,chipNoise);
-                //cout << TString::Format("Threshold= %f", threshold/keV) << endl;
+
+                G4cout << "pixel : " << (*pCItr).first.first << " , " << (*pCItr).first.second
+                                << ", E =  " << ((*pCItr).second)/keV << " keV | thl = "
+                                << TString::Format("Threshold= %f keV", threshold/keV) << G4endl;
                 //if((*pCItr).second > m_digitIn.thl)
                 if((*pCItr).second > threshold)