- Timestamp:
- Apr 15, 2011, 11:13:03 AM (13 years ago)
- Location:
- allpix
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
allpix/include/AllPixFEI3StandardDigitizer.hh
r139 r140 128 128 //G4double maxIntegration = 5; // 129 129 130 G4double sigma; 131 132 130 133 131 134 }; -
allpix/src/AllPixFEI3StandardDigitizer.cc
r139 r140 170 170 171 171 172 FEI3173 FEI4, SlimEdge or Conservative174 Omegapix2172 1=FEI3 173 2,3= FEI4, SlimEdge or Conservative 174 4=Omegapix2 175 175 Default = FEI3 176 176 177 177 */ 178 Sensor = "FEI3";178 Sensor = 1; 179 179 180 180 switch (Sensor) { 181 181 182 case "FEI3":182 case 1 : 183 183 doSlimEdge=false; 184 184 GRShift =0*um; 185 case "FEI4-Conservative":185 case 2 : 186 186 doSlimEdge=true; 187 187 GRShift =250*um; 188 case "FEI4-SlimEdge":188 case 3 : 189 189 doSlimEdge=true; 190 190 GRShift =500*um; 191 case "OmegaPix2":191 case 4 : 192 192 doSlimEdge=true; 193 193 GRShift =500*um; … … 367 367 G4double ypos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().y(); 368 368 G4double zpos = (*hitsCollection)[itr]->GetPosWithRespectToPixel().z(); 369 G4double sigma;370 369 G4double driftTime; 371 370 if(doFullField==false){ … … 657 656 if((xpos+pitchX/2. < GRShift))eEdge = eHit*(1-0.5*(xpos+pitchX/2)/(100*um)); 658 657 else{eEdge=eHit;}; 659 if(!isSingleHit){660 Etemp = IntegrateGaussian(xpos/nm,ypos/nm,sigma/nm,(-pitchX/2.0 + i*pitchX)/nm,(-pitchX/2.+(i+1)*pitchX)/nm,(-pitchY/2 + j*pitchY)/nm,(-pitchY/2 + (j+1)*pitchY)/nm, eEdge );661 }662 else663 {664 658 Etemp=eEdge; 665 }666 659 } 667 660 … … 669 662 if((pitchX/2.-xpos < GRShift))eEdge = eHit*(1-0.5*(pitchX/2.-xpos)/(100*um)); 670 663 else{eEdge=eHit;}; 671 672 if(!isSingleHit){ 673 Etemp = IntegrateGaussian(xpos/nm,ypos/nm,sigma/nm,(-pitchX/2.0 + i*pitchX)/nm,(-pitchX/2.+(i+1)*pitchX)/nm,(-pitchY/2 + j*pitchY)/nm,(-pitchY/2 + (j+1)*pitchY)/nm, eEdge ); 674 } 675 else 676 { 677 Etemp=eEdge; 678 } 664 Etemp=eEdge; 679 665 680 666 }
Note: See TracChangeset
for help on using the changeset viewer.