- Timestamp:
- 12/18/01 17:24:18 (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/manalysis/MMcPedestalCopy.cc
r1118 r1123 26 26 // // 27 27 // MMcPedestalCopy // 28 // This task looks for the ìnformation about FADC pedestals in // 29 // MMcFadcHeader and translates it to the pedestal values in // 30 // MPedestalCam // 28 31 // // 29 32 // Input Containers: // … … 63 66 Bool_t MMcPedestalCopy::PreProcess(MParList *pList) 64 67 { 68 ReInit(pList); 69 70 return kTRUE; 71 } 72 73 Bool_t MMcPedestalCopy::ReInit(MParList *pList) 74 { 65 75 MRawRunHeader *run = (MRawRunHeader*)pList->FindObject("MRawRunHeader"); 66 76 if (run) 77 { 67 78 if (run->GetRunType() != kRTMonteCarlo) 68 79 { … … 70 81 return kSKIP; 71 82 } 83 } 72 84 else 73 85 *fLog << warn << dbginf << "Warning - cannot check file type, MRawRunHeader not found." << endl; 86 74 87 75 88 const MMcFadcHeader *mcped = (MMcFadcHeader*)pList->FindObject("MMcFadcHeader"); … … 94 107 MPedestalPix &pix = (*pedcam)[i]; 95 108 109 // Here one should compute the Pedestal taking into account how 110 // the MC makes the transformation analogic-digital for the FADC. 111 96 112 const Float_t pedest = mcped->GetPedestal(i); 97 const Float_t pedrms = mcped->GetPedestalRms(i);113 const Float_t pedrms = pedest/sqrt(num); 98 114 99 const Float_t sigma = pedest*sqrt(num);115 const Float_t sigma = mcped->GetPedestalRms(i); 100 116 const Float_t sigrms = sigma/sqrt(num*2); 101 117 102 118 pix.SetPedestal(pedest, sigma); 103 119 pix.SetPedestalRms(pedrms, sigrms); 104 120 }
Note:
See TracChangeset
for help on using the changeset viewer.