Changeset 8425 for trunk/MagicSoft/Mars/mcalib
- Timestamp:
- 04/20/07 12:48:09 (18 years ago)
- Location:
- trunk/MagicSoft/Mars/mcalib
- Files:
-
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mcalib/CalibLinkDef.h
r8408 r8425 20 20 #pragma link C++ class MCalibrationIntensityChargeCam+; 21 21 #pragma link C++ class MCalibrationIntensityBlindCam+; 22 #pragma link C++ class MCalibrationIntensityQECam+;23 22 #pragma link C++ class MCalibrationIntensityRelTimeCam+; 24 23 #pragma link C++ class MCalibrationIntensityTestCam+; -
trunk/MagicSoft/Mars/mcalib/MCalibCalcFromPast.cc
r8417 r8425 59 59 #include "MCalibrationIntensityChargeCam.h" 60 60 #include "MCalibrationIntensityBlindCam.h" 61 #include "MCalibrationIntensityQECam.h"62 61 //#include "MCalibrationIntensityRelTimeCam.h" 63 62 64 63 #include "MBadPixelsCam.h" 65 64 65 #include "MCalibrationQECam.h" 66 #include "MCalibrationQEPix.h" 66 67 #include "MCalibrationChargePix.h" 67 68 #include "MCalibrationChargeCalc.h" … … 126 127 127 128 // 128 // Look for the MBadPixels IntensityCam129 // Look for the MBadPixels Cam 129 130 // 130 131 fBadPixels = (MBadPixelsCam*)pList->FindObject("MBadPixelsCam"); … … 132 133 { 133 134 *fLog << err << "MBadPixelsCam not found... abort." << endl; 135 return kFALSE; 136 } 137 138 // 139 // Look for the MBadPixels Intensity Cam 140 // 141 fQECam = (MCalibrationQECam*)pList->FindObject("MCalibrationQECam"); 142 if (!fQECam) 143 { 144 *fLog << err << "MCalibrationQECam not found... abort." << endl; 134 145 return kFALSE; 135 146 } … … 160 171 161 172 fIntensCharge = (MCalibrationIntensityChargeCam*)pList->FindCreateObj("MCalibrationIntensityChargeCam"); 162 fIntensQE = (MCalibrationIntensityQECam*) pList->FindCreateObj("MCalibrationIntensityQECam");163 164 173 if (!fIntensCharge) 165 174 return kFALSE; 166 if (!fIntensQE)167 return kFALSE;168 175 169 176 MCalibrationChargeCam *chargeinit = (MCalibrationChargeCam*)pList->FindObject("MCalibrationChargeCam"); 170 MCalibrationQECam *qeinit = (MCalibrationQECam*) pList->FindObject("MCalibrationQECam");171 177 172 178 if (chargeinit) … … 174 180 else 175 181 *fLog << "Could not find initial MCalibrationChargeCam, cannot initialize intensity cam" << endl; 176 177 if (qeinit)178 fIntensQE->SetCam(qeinit,0);179 else180 *fLog << "Could not find initial MCalibrationQECam, cannot initialize intensity cam" << endl;181 182 182 183 if (!fChargeCalc) … … 386 387 *fLog << "MBadPixelsCam..."; 387 388 fBadPixels->Clear(); // FIXME:::::: MERGE PreExcl!!!! 389 // IS IT REALLY NECESSARY? 390 391 *fLog << "MCalibrationQECam..."; 392 fQECam->Clear(); 393 fQECam->SetRunNumber(runnumber); 394 // IS IT REALLY NECESSARY? 395 396 // IF SIMPLE ENOUGH, REMOVE THE FUNCTION! 388 397 389 398 if (fIntensCharge) … … 404 413 } 405 414 406 if (fIntensQE)407 {408 *fLog << "MCalibrationQECam...";409 MCalibrationCam *cold = fIntensQE->GetCam();410 fIntensQE->AddToList(Form("MCalibrationQECam%04d",fNumCam),*fGeom);411 fIntensQE->GetCam()->SetRunNumber(runnumber);412 if (cold)413 fIntensQE->DeleteOldCam(cold);414 }415 415 if (fIntensBlind) 416 416 { -
trunk/MagicSoft/Mars/mcalib/MCalibCalcFromPast.h
r8417 r8425 19 19 class MCalibrationIntensityChargeCam; 20 20 class MCalibrationIntensityBlindCam; 21 class MCalibration IntensityQECam;21 class MCalibrationQECam; 22 22 23 23 class MBadPixelsCam; … … 34 34 MCalibrationIntensityChargeCam *fIntensCharge; //! Intensity Charge Cam (to be created) 35 35 MCalibrationIntensityBlindCam *fIntensBlind; //! Intensity Blind Cam (to be created) 36 MCalibrationIntensityQECam *fIntensQE; //! Intensity QE Cam (to be created)37 36 38 MBadPixelsCam *fBadPixels; 37 MCalibrationQECam *fQECam; //! QE Cam 38 MBadPixelsCam *fBadPixels; //! Bad Pixels Cam 39 39 40 40 MCalibrationChargeCalc *fChargeCalc; //! Charge Calibration Task -
trunk/MagicSoft/Mars/mcalib/MCalibColorSteer.cc
r8417 r8425 59 59 #include "MCalibrationIntensityChargeCam.h" 60 60 #include "MCalibrationIntensityBlindCam.h" 61 #include "MCalibrationIntensityQECam.h"62 61 #include "MCalibrationIntensityRelTimeCam.h" 63 62 64 63 #include "MCalibrationPattern.h" 64 #include "MCalibrationQECam.h" 65 65 #include "MCalibrationChargeCam.h" 66 66 #include "MCalibrationChargeCalc.h" … … 156 156 if (pList->FindObject(AddSerialNumber("MHCalibrationChargeCam"))) 157 157 { 158 159 fIntensCharge = (MCalibrationIntensityChargeCam*)pList->FindCreateObj("MCalibrationIntensityChargeCam"); 160 fIntensQE = (MCalibrationIntensityQECam*) pList->FindCreateObj("MCalibrationIntensityQECam"); 161 162 fChargeCalc = (MCalibrationChargeCalc*)tlist->FindObject("MCalibrationChargeCalc"); 163 164 *fLog << inf << "Found MHCalibrationChargeCam ... " << flush; 165 166 if (!fIntensCharge) 158 *fLog << inf << "Found MHCalibrationChargeCam ... " << flush; 159 160 fIntensCharge = (MCalibrationIntensityChargeCam*)pList->FindCreateObj("MCalibrationIntensityChargeCam"); 161 if (!fIntensCharge) 162 return kFALSE; 163 164 fQECam = (MCalibrationQECam*)pList->FindCreateObj("MCalibrationQECam"); 165 if (!fQECam) 166 return kFALSE; 167 168 fChargeCalc = (MCalibrationChargeCalc*)tlist->FindObject("MCalibrationChargeCalc"); 169 if (!fChargeCalc) 167 170 { 168 *fLog << err << "Could not find nor create MCalibrationIntensityChargeCam abort... " << endl; 169 return kFALSE; 170 } 171 172 if (!fIntensQE) 173 { 174 *fLog << err << "Could not find nor create MCalibrationIntensityQECam abort... " << endl; 175 return kFALSE; 176 } 177 178 if (!fChargeCalc) 179 { 180 *fLog << err << "Could not find MCalibrationChargeCalc abort... " << endl; 181 return kFALSE; 171 *fLog << err << "Could not find MCalibrationChargeCalc abort... " << endl; 172 return kFALSE; 182 173 } 183 174 } 184 175 185 176 // 186 177 // Look for the MFillH name "FillRelTimeCam". In case yes, initialize the … … 327 318 328 319 fBad->Clear(); // FIXME:::::: MERGE PreExcl!!!! 320 fQECam->Clear(); 329 321 330 322 if (fIntensCharge) … … 337 329 return kFALSE; 338 330 *fLog << inf << "New MCalibrationChargeCam with name: " << cam->GetName() << endl; 339 }340 if (fIntensQE)341 {342 fIntensQE->AddToList(Form("MCalibrationQECam%s",namep.Data()),*fGeom);343 *fLog << inf << "New MCalibrationQECam with: " << namep << endl;344 331 } 345 332 -
trunk/MagicSoft/Mars/mcalib/MCalibColorSteer.h
r8417 r8425 19 19 class MCalibrationIntensityChargeCam; 20 20 class MCalibrationIntensityBlindCam; 21 class MCalibrationIntensityQECam;22 21 class MCalibrationIntensityRelTimeCam; 22 class MCalibrationQECam; 23 23 class MBadPixelsCam; 24 24 … … 32 32 MCalibrationIntensityChargeCam *fIntensCharge; //! 33 33 MCalibrationIntensityBlindCam *fIntensBlind; //! 34 MCalibrationIntensityQECam *fIntensQE; //!35 34 MCalibrationIntensityRelTimeCam *fIntensRelTime; //! 35 MCalibrationQECam *fQECam; //! 36 36 MBadPixelsCam *fBad; //! 37 37 -
trunk/MagicSoft/Mars/mcalib/MCalibrateData.cc
r8408 r8425 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MCalibrateData.cc,v 1.6 6 2007-04-17 12:39:14tbretz Exp $2 ! $Name: not supported by cvs2svn $:$Id: MCalibrateData.cc,v 1.67 2007-04-20 11:48:08 tbretz Exp $ 3 3 ! -------------------------------------------------------------------------- 4 4 ! … … 118 118 #include "MCalibrationChargePix.h" 119 119 120 #include "MCalibrationIntensityQECam.h"121 120 #include "MCalibrationQECam.h" 122 121 #include "MCalibrationQEPix.h" … … 157 156 MCalibrateData::MCalibrateData(CalibrationMode_t calmode,const char *name, const char *title) 158 157 : fGeomCam(NULL), fBadPixels(NULL), fCalibrations(NULL), fIntensCalib(NULL), 159 fQEs(NULL), f IntensQE(NULL), fSignals(NULL), fCerPhotEvt(NULL), fCalibConstCam(NULL),158 fQEs(NULL), fSignals(NULL), fCerPhotEvt(NULL), fCalibConstCam(NULL), 160 159 fPedestalFlag(kNo), fSignalType(kPhot), fRenormFactor(1.), fScaleFactor(1.) 161 160 { … … 249 248 return kFALSE; 250 249 } 251 252 fIntensQE = (MCalibrationIntensityQECam*)pList->FindObject(AddSerialNumber("MCalibrationIntensityQECam"));253 if (fIntensQE)254 *fLog << inf << "Found MCalibrationIntensityQECam ... " << endl;255 250 256 251 fQEs = (MCalibrationQECam*)pList->FindObject(AddSerialNumber("MCalibrationQECam")); … … 370 365 } 371 366 372 const MCalibrationQECam *qecam = fIntensQE ? (MCalibrationQECam*)fIntensQE->GetCam() : fQEs; 373 if(fCalibrationMode == kBlindPixel && !qecam->IsBlindPixelMethodValid()) 367 if(fCalibrationMode == kBlindPixel && !fQEs->IsBlindPixelMethodValid()) 374 368 { 375 369 *fLog << warn << "Blind pixel calibration method not valid, switching to F-factor method" << endl; … … 377 371 } 378 372 379 if(fCalibrationMode == kPinDiode && ! qecam->IsPINDiodeMethodValid())373 if(fCalibrationMode == kPinDiode && !fQEs->IsPINDiodeMethodValid()) 380 374 { 381 375 *fLog << warn << "PIN diode calibration method not valid, switching to F-factor method" << endl; … … 383 377 } 384 378 385 if(fCalibrationMode == kCombined && ! qecam->IsCombinedMethodValid())379 if(fCalibrationMode == kCombined && !fQEs->IsCombinedMethodValid()) 386 380 { 387 381 *fLog << warn << "Combined calibration method not valid, switching to F-factor method" << endl; … … 520 514 521 515 MCalibrationChargeCam *chargecam = NULL; 522 MCalibrationQECam *qecam = NULL;523 516 if (updatecam) 524 {525 517 chargecam = fCalibrations; 526 qecam = fQEs;527 }528 518 else 529 {530 519 chargecam = fIntensCalib ? (MCalibrationChargeCam*)fIntensCalib->GetCam() : fCalibrations; 531 qecam = fIntensQE ? (MCalibrationQECam*) fIntensQE->GetCam() : fQEs;532 }533 520 534 521 // … … 569 556 calibFFactor = pix.GetMeanFFactorFADC2Phot(); 570 557 571 const MCalibrationQEPix &qe = (MCalibrationQEPix&)(* qecam)[pixidx];558 const MCalibrationQEPix &qe = (MCalibrationQEPix&)(*fQEs)[pixidx]; 572 559 573 560 if (updatecam) -
trunk/MagicSoft/Mars/mcalib/MCalibrateData.h
r8408 r8425 22 22 23 23 class MCalibrationIntensityChargeCam; 24 class MCalibrationIntensityQECam;25 24 26 25 class MCalibrateData : public MTask … … 66 65 MCalibrationIntensityChargeCam *fIntensCalib; //! Calibration constants from interlaced calib. 67 66 MCalibrationQECam *fQEs; //! Quantum efficiencies 68 MCalibrationIntensityQECam *fIntensQE; //! Quantum efficiencies from interlaced calib.69 67 MExtractedSignalCam *fSignals; //! Integrated charge in FADCs counts 70 68 MSignalCam *fCerPhotEvt; //! Cerenkov Photon Event used for calculation -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc
r8423 r8425 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MCalibrationChargeCalc.cc,v 1.17 4 2007-04-20 08:55:18 tbretz Exp $2 ! $Name: not supported by cvs2svn $:$Id: MCalibrationChargeCalc.cc,v 1.175 2007-04-20 11:48:08 tbretz Exp $ 3 3 ! -------------------------------------------------------------------------- 4 4 ! … … 222 222 223 223 #include "MCalibrationIntensityChargeCam.h" 224 #include "MCalibrationIntensityQECam.h"225 224 #include "MCalibrationIntensityBlindCam.h" 226 225 … … 367 366 fCam = NULL; 368 367 fHCam = NULL; 369 fIntensQE = NULL;370 368 fQECam = NULL; 371 369 fIntensBlind = NULL; … … 433 431 } 434 432 433 // Obsolete... it is already set in the constructor isn't it? 435 434 fPulserColor = MCalibrationCam::kNONE; 436 435 … … 444 443 // - MGeomCam 445 444 // - MCalibrationIntensityChargeCam or MCalibrationChargeCam 446 // - MCalibration IntensityQECam or MCalibrationQECam445 // - MCalibrationQECam 447 446 // - MBadPixelsCam 448 447 // … … 476 475 { 477 476 477 // 478 // Geometry 479 // 478 480 fGeom = (MGeomCam*)pList->FindObject("MGeomCam"); 479 481 if (!fGeom) … … 482 484 return kFALSE; 483 485 } 484 486 487 // 488 // Charge histogram 489 // 490 fHCam = (MHCalibrationChargeCam*)pList->FindObject(AddSerialNumber("MHCalibrationChargeCam")); 491 if (!fHCam) 492 { 493 *fLog << err << "Cannot find MHCalibrationChargeCam ... abort." << endl; 494 *fLog << "Maybe you forget to call an MFillH for the MHCalibrationChargeCam before..." << endl; 495 return kFALSE; 496 } 497 498 // 499 // Intensity containers 500 // 485 501 fIntensCam = (MCalibrationIntensityChargeCam*)pList->FindObject(AddSerialNumber("MCalibrationIntensityChargeCam")); 486 502 if (fIntensCam) … … 497 513 } 498 514 499 fHCam = (MHCalibrationChargeCam*)pList->FindObject(AddSerialNumber("MHCalibrationChargeCam"));500 if (!fHCam)501 {502 *fLog << err << "Cannot find MHCalibrationChargeCam ... abort." << endl;503 *fLog << "Maybe you forget to call an MFillH for the MHCalibrationChargeCam before..." << endl;504 return kFALSE;505 }506 507 fIntensQE = (MCalibrationIntensityQECam*)pList->FindObject(AddSerialNumber("MCalibrationIntensityQECam"));508 if (fIntensQE)509 *fLog << inf << "Found MCalibrationIntensityQECam... " << flush;510 else511 {512 fQECam = (MCalibrationQECam*)pList->FindObject(AddSerialNumber("MCalibrationQECam"));513 if (!fQECam)514 {515 *fLog << err << "Cannot find MCalibrationQECam ... abort." << endl;516 *fLog << "Maybe you forget to call an MFillH for the MHCalibrationQECam before..." << endl;517 return kFALSE;518 }519 }520 521 515 fIntensBlind = (MCalibrationIntensityBlindCam*)pList->FindObject(AddSerialNumber("MCalibrationIntensityBlindCam")); 522 516 if (fIntensBlind) … … 528 522 *fLog << inf << "No MCalibrationBlindCam found... no Blind Pixel method!" << endl; 529 523 } 530 531 fHBlindCam = (MHCalibrationChargeBlindCam*)pList->FindObject(AddSerialNumber("MHCalibrationChargeBlindCam")); 532 if (!fHBlindCam) 533 *fLog << inf << "No MHCalibrationChargeBlindCam found... no Blind Pixel method!" << endl; 524 525 // 526 // Required containers 527 // 528 fQECam = (MCalibrationQECam*)pList->FindObject(AddSerialNumber("MCalibrationQECam")); 529 if (!fQECam) 530 { 531 *fLog << err << "Cannot find MCalibrationQECam ... abort." << endl; 532 return kFALSE; 533 } 534 534 535 535 fBadPixels = (MBadPixelsCam*)pList->FindObject(AddSerialNumber("MBadPixelsCam")); … … 546 546 if (!fPINDiode) 547 547 *fLog << inf << "No MCalibrationChargePINDiode found... no PIN Diode method!" << endl; 548 549 MCalibrationQECam *qecam = fIntensQE 550 ? (MCalibrationQECam*) fIntensQE->GetCam() : fQECam; 548 549 fHBlindCam = (MHCalibrationChargeBlindCam*)pList->FindObject(AddSerialNumber("MHCalibrationChargeBlindCam")); 550 if (!fHBlindCam) 551 *fLog << inf << "No MHCalibrationChargeBlindCam found... no Blind Pixel method!" << endl; 552 553 554 // 555 // Further initialization 556 // 551 557 MCalibrationChargeCam *chargecam = fIntensCam 552 558 ? (MCalibrationChargeCam*)fIntensCam->GetCam() : fCam; 553 559 554 UInt_t npixels= fGeom->GetNumPixels();560 const UInt_t npixels = fGeom->GetNumPixels(); 555 561 556 562 for (UInt_t i=0; i<npixels; i++) 557 563 { 558 564 559 565 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*chargecam)[i]; 560 MCalibrationQEPix &pqe = (MCalibrationQEPix&) (*qecam) [i]; 561 MBadPixelsPix &bad = (*fBadPixels)[i]; 562 563 if (bad.IsBad()) 564 { 566 567 if ((*fBadPixels)[i].IsBad()) 568 { 569 MCalibrationQEPix &pqe = (MCalibrationQEPix&)(*fQECam)[i]; 565 570 pix.SetExcluded(); 566 571 pqe.SetExcluded(); 567 572 continue; 568 569 573 } 574 570 575 if (IsDebug()) 571 576 pix.SetDebug(); … … 585 590 586 591 const MCalibrationCam::PulserColor_t col = fCalibPattern->GetPulserColor(); 592 593 // 594 // If pulser color didn't change there is nothing to do 595 // 587 596 const Float_t strength = fCalibPattern->GetPulserStrength(); 588 597 const Float_t strdiff = TMath::Abs(strength-fStrength); 589 590 598 if (col == fPulserColor && strdiff < 0.05 ) 591 599 { … … 593 601 return kTRUE; 594 602 } 595 603 604 // 605 // If there is no color is nothing to do 606 // 596 607 if (col == MCalibrationCam::kNONE) 597 608 return kTRUE; 598 599 609 600 610 // … … 627 637 } 628 638 } 629 639 640 // 641 // Keep pulser color and strength in mind and output some information 642 // 630 643 fPulserColor = col; 631 644 fStrength = strength; 632 645 633 646 *fLog << inf << "Found new colour ... " << flush; 634 647 635 648 switch (col) 636 649 { … … 644 657 *fLog << inf << " with strength: " << strength << endl; 645 658 659 // 660 // Distribute new color to all containers 661 // 646 662 fHCam->SetColor(col); 647 663 if (fHBlindCam) … … 660 676 fPINDiode->SetColor(col); 661 677 678 // 679 // Reset number of processed events 680 // 662 681 fNumProcessed = 0; 663 682 … … 738 757 739 758 MCalibrationBlindCam *blindcam = fIntensBlind ? (MCalibrationBlindCam*) fIntensBlind->GetCam() : fBlindCam; 740 MCalibrationQECam *qecam = fIntensQE ? (MCalibrationQECam*) fIntensQE->GetCam() : fQECam;741 759 MCalibrationChargeCam *chargecam = fIntensCam ? (MCalibrationChargeCam*)fIntensCam->GetCam() : fCam; 742 760 … … 844 862 // Finalize Blind Pixel 845 863 // 846 qecam->SetBlindPixelMethodValid(FinalizeBlindCam());864 fQECam->SetBlindPixelMethodValid(FinalizeBlindCam()); 847 865 848 866 // 849 867 // Finalize PIN Diode 850 868 // 851 qecam->SetBlindPixelMethodValid(FinalizePINDiode());869 fQECam->SetBlindPixelMethodValid(FinalizePINDiode()); 852 870 853 871 // … … 907 925 908 926 chargecam->SetReadyToSave(); 909 qecam ->SetReadyToSave();927 fQECam ->SetReadyToSave(); 910 928 fBadPixels->SetReadyToSave(); 911 929 … … 1720 1738 } 1721 1739 1722 MCalibrationQECam *qecam = fIntensQE1723 ? (MCalibrationQECam*) fIntensQE->GetCam() : fQECam;1724 1740 MCalibrationChargeCam *chargecam = fIntensCam 1725 1741 ? (MCalibrationChargeCam*)fIntensCam->GetCam() : fCam; 1726 1742 1727 1743 MCalibrationChargePix &avpix = (MCalibrationChargePix&)chargecam->GetAverageArea(0); 1728 MCalibrationQEPix &qepix = (MCalibrationQEPix&) qecam->GetAverageArea(0);1744 MCalibrationQEPix &qepix = (MCalibrationQEPix&) fQECam->GetAverageArea(0); 1729 1745 1730 1746 if (IsDebug()) … … 1738 1754 / qepix.GetPMTCollectionEff() 1739 1755 / qepix.GetLightGuidesEff(fPulserColor) 1740 / qecam->GetPlexiglassQE();1756 / fQECam->GetPlexiglassQE(); 1741 1757 1742 1758 const Float_t avphotrelvar = ( IsUseExternalNumPhes() … … 1746 1762 + qepix.GetPMTCollectionEffRelVar() 1747 1763 + qepix.GetLightGuidesEffRelVar(fPulserColor) 1748 + qecam->GetPlexiglassQERelVar();1764 + fQECam->GetPlexiglassQERelVar(); 1749 1765 1750 1766 const UInt_t nareas = fGeom->GetNumAreas(); … … 1772 1788 1773 1789 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*chargecam)[i]; 1774 MCalibrationQEPix &qpix = (MCalibrationQEPix&) (* qecam)[i];1790 MCalibrationQEPix &qpix = (MCalibrationQEPix&) (*fQECam) [i]; 1775 1791 1776 1792 MBadPixelsPix &bad = (*fBadPixels)[i]; … … 1788 1804 qpix.SetFFactorMethodValid( kTRUE , fPulserColor ); 1789 1805 1790 if (!qpix.UpdateFFactorMethod( qecam->GetPlexiglassQE() ))1806 if (!qpix.UpdateFFactorMethod( fQECam->GetPlexiglassQE() )) 1791 1807 *fLog << warn << GetDescriptor() 1792 1808 << ": Cannot update Quantum efficiencies with the F-Factor Method" << endl; … … 1917 1933 1918 1934 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*chargecam)[i]; 1919 MCalibrationQEPix &qpix = (MCalibrationQEPix&) (* qecam)[i];1935 MCalibrationQEPix &qpix = (MCalibrationQEPix&) (*fQECam) [i]; 1920 1936 1921 1937 if ((*fBadPixels)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun)) … … 1953 1969 MCalibrationBlindCam *blindcam = fIntensBlind 1954 1970 ? (MCalibrationBlindCam*) fIntensBlind->GetCam(): fBlindCam; 1955 MCalibrationQECam *qecam = fIntensQE1956 ? (MCalibrationQECam*) fIntensQE->GetCam() : fQECam;1957 1971 MCalibrationChargeCam *chargecam = fIntensCam 1958 1972 ? (MCalibrationChargeCam*)fIntensCam->GetCam() : fCam; … … 1968 1982 1969 1983 const Float_t photons = blindcam->GetFluxInsidePlexiglass() * (*fGeom)[0].GetA() 1970 / qecam->GetPlexiglassQE();1984 / fQECam->GetPlexiglassQE(); 1971 1985 chargecam->SetNumPhotonsBlindPixelMethod(photons); 1972 1986 1973 1987 const Float_t photrelvar = blindcam->GetFluxInsidePlexiglassRelVar() 1974 + qecam->GetPlexiglassQERelVar();1988 + fQECam->GetPlexiglassQERelVar(); 1975 1989 1976 1990 if (photrelvar > 0.) … … 1986 2000 { 1987 2001 1988 MCalibrationQEPix &qepix = (MCalibrationQEPix&) (*qecam)[i];2002 MCalibrationQEPix &qepix = (MCalibrationQEPix&)(*fQECam)[i]; 1989 2003 1990 2004 if (!blindcam || !(blindcam->IsFluxInsidePlexiglassAvailable())) … … 2008 2022 / blindcam->GetFluxInsidePlexiglass() 2009 2023 / geo.GetA() 2010 * qecam->GetPlexiglassQE();2024 * fQECam->GetPlexiglassQE(); 2011 2025 2012 2026 const Float_t qerelvar = blindcam->GetFluxInsidePlexiglassRelVar() 2013 + qecam->GetPlexiglassQERelVar()2027 + fQECam->GetPlexiglassQERelVar() 2014 2028 + pix.GetPheFFactorMethodRelVar(); 2015 2029 … … 2018 2032 qepix.SetBlindPixelMethodValid( kTRUE , fPulserColor ); 2019 2033 2020 if (!qepix.UpdateBlindPixelMethod( qecam->GetPlexiglassQE()))2034 if (!qepix.UpdateBlindPixelMethod( fQECam->GetPlexiglassQE())) 2021 2035 *fLog << warn << GetDescriptor() 2022 2036 << ": Cannot update Quantum efficiencies with the Blind Pixel Method" << endl; … … 2046 2060 const UInt_t npixels = fGeom->GetNumPixels(); 2047 2061 2048 MCalibrationQECam *qecam = fIntensQE2049 ? (MCalibrationQECam*) fIntensQE->GetCam() : fQECam;2050 2062 MCalibrationChargeCam *chargecam = fIntensCam 2051 2063 ? (MCalibrationChargeCam*)fIntensCam->GetCam() : fCam; … … 2061 2073 { 2062 2074 2063 MCalibrationQEPix &qepix = (MCalibrationQEPix&) (*qecam)[i];2075 MCalibrationQEPix &qepix = (MCalibrationQEPix&)(*fQECam)[i]; 2064 2076 2065 2077 if (!fPINDiode) … … 2113 2125 const UInt_t npixels = fGeom->GetNumPixels(); 2114 2126 2115 MCalibrationQECam *qecam = fIntensQE2116 ? (MCalibrationQECam*) fIntensQE->GetCam() : fQECam;2117 2118 2127 for (UInt_t i=0; i<npixels; i++) 2119 2128 { 2120 2129 2121 MCalibrationQEPix &qepix = (MCalibrationQEPix&) (*qecam)[i];2130 MCalibrationQEPix &qepix = (MCalibrationQEPix&)(*fQECam)[i]; 2122 2131 2123 2132 if (!(*fBadPixels)[i].IsUnsuitable (MBadPixelsPix::kUnsuitableRun)) -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.h
r8417 r8425 29 29 class MCalibrationChargeCam; 30 30 class MHCalibrationChargeCam; 31 class MCalibrationIntensityQECam;32 31 class MCalibrationQECam; 33 32 class MGeomCam; … … 93 92 MHCalibrationChargeBlindCam *fHBlindCam; //! Charges histograms of the Blind Pixels 94 93 MCalibrationChargePINDiode *fPINDiode; //! Calibrated Charges of the PIN Diode 95 MCalibrationIntensityQECam *fIntensQE; //! Intensity Calibration Quantum Efficiencies of all pixels96 94 MCalibrationQECam *fQECam; //! Calibrated Quantum Efficiencies of all pixels 97 95 MGeomCam *fGeom; //! Camera geometry -
trunk/MagicSoft/Mars/mcalib/MCalibrationIntensityCam.cc
r8140 r8425 34 34 // per calibration camera result and area index 35 35 // 36 // See also: MCalibrationIntensityChargeCam, MCalibrationIntensityQECam, 37 // MCalibrationIntensityRelTimeCam, 38 // MCalibrationCam, MCalibrationPix, 39 // MCalibrationQECam, MCalibrationQEPix, 40 // MHCalibrationChargePix, MHCalibrationChargeCam 41 // MCalibrationChargeBlindPix, MCalibrationChargePINDiode 36 // See also: MCalibrationIntensityChargeCam, MCalibrationIntensityRelTimeCam, 37 // MCalibrationCam, MCalibrationPix, MHCalibrationChargePix, 38 // MHCalibrationChargeCam, MCalibrationChargeBlindPix, 39 // MCalibrationChargePINDiode 42 40 // 43 41 // Inline Functions: -
trunk/MagicSoft/Mars/mcalib/Makefile
r8408 r8425 42 42 MCalibrationIntensityChargeCam.cc \ 43 43 MCalibrationIntensityBlindCam.cc \ 44 MCalibrationIntensityQECam.cc \45 44 MCalibrationIntensityRelTimeCam.cc \ 46 45 MCalibrationIntensityTestCam.cc \
Note:
See TracChangeset
for help on using the changeset viewer.