Changeset 9936 for trunk/Mars


Ignore:
Timestamp:
09/22/10 09:15:57 (14 years ago)
Author:
tbretz
Message:
Implemented typed 6-8 in MHPhotonEvent
Location:
trunk/Mars
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/Changelog

    r9935 r9936  
    3939     - accept = as a conditional sign
    4040     - make sure the return code of the evaluated filter is corectly returned
     41
     42   * msim/MHPhotonEvent.cc:
     43     - implemented types 6-8
    4144
    4245
  • trunk/Mars/msim/MHPhotonEvent.cc

    r9902 r9936  
    228228
    229229    case 4: // Two times the pixel-0 traversal size
    230     case 5: // The maximum radius (cm)
     230    case 5: // The maximum radius (cm) -- ~ 10 bins / pix
     231    case 6: // The maximum radius (cm) -- ~  5 bins / pix
     232    case 7: // The maximum radius (cm) -- ~  3 bins / pix
     233    case 8: // The maximum radius (cm) -- ~  2 bins / pix
    231234        {
    232235            MGeomCam *c = (MGeomCam*)pList->FindObject("MGeomCam");
     
    237240            }
    238241            // Type 3: Define ~10 bins per pixel
    239             xmax = fType==3 || fType==5 ? c->GetMaxRadius() : 4*(*c)[0].GetT()/10;
    240             num  = fType==3 || fType==5 ? TMath::Nint(10*(2*c->GetMaxRadius())/(2*(*c)[0].GetT())) : 101;
    241 
    242             if (fType==5)
     242            xmax = fType!=4 ? c->GetMaxRadius() : 4*(*c)[0].GetT()/10;
     243
     244            switch (fType)
     245            {
     246            case 3:
     247            case 5:
     248                // ~10 bins per pixel
     249                num = TMath::Nint(10*c->GetMaxRadius()/(*c)[0].GetT());
     250                break;
     251            case 6:
     252                // ~5 bins per pixel
     253                num = TMath::Nint(5*c->GetMaxRadius()/(*c)[0].GetT());
     254                break;
     255            case 7:
     256                // ~3 bins per pixel
     257                num = TMath::Nint(3*c->GetMaxRadius()/(*c)[0].GetT());
     258                break;
     259            case 8:
     260                // ~2 bins per pixel
     261                num = TMath::Nint(2*c->GetMaxRadius()/(*c)[0].GetT());
     262                break;
     263            case 4:
     264                num = 101;
     265                break;
     266            }
     267
     268            if (fType==5 || fType==6 || fType==7)
    243269                xmax /= 10;
    244270
Note: See TracChangeset for help on using the changeset viewer.