Changeset 1311 for trunk/MagicSoft/Simulation/Detector/Camera
- Timestamp:
- 04/27/02 11:48:39 (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx
r1260 r1311 21 21 // 22 22 // $RCSfile: camera.cxx,v $ 23 // $Revision: 1.3 8$23 // $Revision: 1.39 $ 24 24 // $Author: blanch $ 25 // $Date: 2002-0 3-18 18:44:29 $25 // $Date: 2002-04-27 10:48:39 $ 26 26 // 27 27 //////////////////////////////////////////////////////////////////////// … … 170 170 static float ct_PixelWidth_corner_2_corner; 171 171 172 //@: ct_PixelWidth_corner_2_corner / 2173 static float ct_PixelWidth_corner_2_corner_half;174 175 172 //@: Number of mirrors 176 173 static int ct_NMirrors = 0; … … 184 181 //@: Number of pixels 185 182 static int ct_NGapPixels; 186 187 //@: ct_Apot = ct_PixelWidth / 2188 static float ct_Apot;189 190 //@: ct_2Apot = 2 * ct_Apot = ct_PixelWidth191 static float ct_2Apot;192 183 193 184 //@: name of the CT definition file to use … … 421 412 float thetaCT, phiCT; //@< parameters of a given shower 422 413 float thetashw, phishw; //@< parameters of a given shower 423 float core D, coreX, coreY; //@< core position414 float coreX, coreY; //@< core position 424 415 float impactD; //@< impact parameter 425 416 float l1, m1, n1; //@< auxiliary variables … … 464 455 465 456 float qTailCut; //@< Tail Cut value 466 int nIslandsCut; //@< Islands Cut value467 int countIslands=FALSE; //@< Will we count the islands?468 457 int anaPixels; 469 458 … … 478 467 int isorttopo[3]; //@< sorting the topologies 479 468 int icontrigger; //@< number of trigger conditions to be analised 480 UShort_t numPix; //@< number of sets of fadc written counts481 469 482 470 float fpixelthres[CAMERA_PIXELS]; //@< Threshold values … … 484 472 TArrayC *fadcValues; //@< the analog Fadc signal for pixels 485 473 486 float plateScale_cm2deg; //@< plate scale (deg/cm)487 474 float degTriggerZone; //@< trigger area in the camera (radius, in deg.) 488 475 … … 631 618 addElecNoise = add_elec_noise(&FADC_noise, &Trigger_noise); 632 619 simulateNSB = get_nsb( &meanNSB, &nphe2NSB ); 633 countIslands = get_islands_cut( &nIslandsCut );634 620 635 621 // get selections on the parameters … … 765 751 read_pixels(&cam); 766 752 767 Int_t Lev0, Lev1 , Lev2;753 Int_t Lev0, Lev1; 768 754 769 755 fadcValues = new TArrayC(FADC_SLICES); … … 1090 1076 } 1091 1077 1092 // set plate scale (deg/cm) and trigger area (deg) 1093 1094 plateScale_cm2deg = ( ct_Type == 0 ) ? (0.244/2.1) : 0.030952381; 1078 // set trigger area (deg) 1095 1079 1096 1080 if ( ! get_trigger_radius( °TriggerZone ) ) … … 1234 1218 fread( (char*)&mcevth, mcevth.mysize(), 1, inputfile ); 1235 1219 1236 // calculate core distance and impact parameter1237 1238 coreD = mcevth.get_core(&coreX, &coreY);1239 1240 1220 // calculate impact parameter (shortest distance betwee the original 1241 1221 // trajectory of the primary (assumed shower-axis) and the … … 1286 1266 // therefore, we calculate the impact parameter as the distance 1287 1267 // between the CT axis and the core position 1288 1268 1269 mcevth.get_core(&coreX, &coreY); 1270 1289 1271 impactD = dist_r_P( 0., 0., 0., l1, m1, n1, coreX, coreY, 0. ); 1290 1272 … … 1456 1438 1457 1439 // Set trigger flags to zero 1458 Lev1= Lev2=0;1440 Lev1=0; 1459 1441 1460 1442 // loop over multiplicity of trigger configuration … … 1468 1450 // loop over topologies 1469 1451 for(itopocount=Trigger_loop_ltop;itopocount<=Trigger_loop_utop;itopocount++){ 1470 Lev1= Lev2=0;1452 Lev1=0; 1471 1453 1472 1454 if(itopocount==0 && imulticount>7) continue; … … 1500 1482 Lev0=0; 1501 1483 } 1502 numPix=0; 1484 1503 1485 for (Int_t ii=0;ii<Lev1;ii++){ 1504 1486 if (Write_McTrig){ … … 1605 1587 Lev0 = (Short_t) Trigger.ZeroLevel() ; 1606 1588 1607 Lev1 = Lev2 =0 ;1589 Lev1 = 0 ; 1608 1590 1609 1591 // … … 1623 1605 } 1624 1606 1625 numPix=0;1626 1607 for(Int_t ii=0;ii<Lev1;ii++){ 1627 1608 // Loop over different level one triggers … … 2219 2200 2220 2201 ct_PixelWidth_corner_2_corner = ct_PixelWidth / cos(RAD(30.0)); 2221 ct_PixelWidth_corner_2_corner_half =2222 ct_PixelWidth_corner_2_corner * 0.50;2223 ct_Apot = ct_PixelWidth / 2;2224 ct_2Apot = ct_Apot * 2.0;2225 2202 2226 2203 log( "read_ct_file", "Pixel width: %f cm\n", ct_PixelWidth ); … … 2516 2493 double dxfirst_pix; 2517 2494 double dyfirst_pix; 2518 double ddxseg1, ddxseg2, ddxseg3, ddxseg4, ddxseg5 , ddxseg6;2519 double ddyseg1, ddyseg2, ddyseg3, ddyseg4, ddyseg5 , ddyseg6;2495 double ddxseg1, ddxseg2, ddxseg3, ddxseg4, ddxseg5; 2496 double ddyseg1, ddyseg2, ddyseg3, ddyseg4, ddyseg5; 2520 2497 2521 2498 … … 2597 2574 ddxseg5 = -ddxseg2; 2598 2575 ddyseg5 = 0.; 2599 ddxseg6 = -ddxseg1;2600 ddyseg6 = ddyseg1;2601 2576 2602 2577 /* to find the position of pixel ipixno take the position of the first pixel in the ring and move */ … … 3378 3353 // 3379 3354 // $Log: not supported by cvs2svn $ 3355 // Revision 1.38 2002/03/18 18:44:29 blanch 3356 // Small modificatin to set the electronic Noise in the MMcTrigHeader class. 3357 // 3380 3358 // Revision 1.37 2002/03/18 16:42:20 blanch 3381 3359 // The data member fProductionSite of the MMcRunHeader has been set to 0, … … 3491 3469 // 3492 3470 // $Log: not supported by cvs2svn $ 3471 // Revision 1.38 2002/03/18 18:44:29 blanch 3472 // Small modificatin to set the electronic Noise in the MMcTrigHeader class. 3473 // 3493 3474 // Revision 1.37 2002/03/18 16:42:20 blanch 3494 3475 // The data member fProductionSite of the MMcRunHeader has been set to 0,
Note:
See TracChangeset
for help on using the changeset viewer.