Changeset 3374 for trunk/MagicSoft/Mars


Ignore:
Timestamp:
03/01/04 17:18:18 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
16 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r3372 r3374  
    1313     - fixed an include (from MArrivalTime.h to MArrivalTimeCam.h)
    1414
     15
     16
    1517 2004/03/01: Wolfgang Wittek
    1618 
    … …  
    4648   * mhist/MHStarMap.[h,cc]:
    4749     - changes to use MObservatory member function
     50
     51   * manalysis/MCerPhotAnal.cc, manalysis/MCerPhotAnal2.cc,
     52     manalysis/MCerPhotCalc.cc, manalysis/MMcCalibrationUpdate.cc,
     53     manalysis/MMcPedestalCopy.cc, manalysis/MMcPedestalNSBAdd.cc,
     54     manalysis/MPedCalcPedRun.cc, manalysis/MPedPhotCalc.cc,
     55     mcalib/MCalibrationChargeCalc.cc, mcalib/MMcCalibrationCalc.cc
     56     - changed to use IsMonteCarloRun() now
     57
     58   * manalysis/MMcCalibrationUpdate.cc:
     59     - unified output to log-stream
     60     - replaced sqrt by TMath::Sqrt
     61
     62   * mcalib/MCalibrationChargeCalc.cc:
     63     - unified output to log-stream
     64     - removed .Data() from TString where obsolete
     65     - fixed memory leak using gSystem->ExpandPathName
     66
     67   * mpointing/MPointingPosCalc.cc, mraw/MRawFileWrite.cc:
     68     - adde class name to used kRT*/kTT* enums
     69
     70   * mraw/MRawEvtHeader.h:
     71     - made enum data member of class
     72
     73   * mraw/MRawRunHeader.[h,cc]:
     74     - made kMagicNumber and kMaxFormatVersion a static data
     75       member of the class
     76     - added IsMonteCarloRun member function
    4877
    4978

  • trunk/MagicSoft/Mars/manalysis/MCerPhotAnal.cc

    r2942 r3374  
    103103        return kFALSE;
    104104
     105    fPedestals = NULL;
     106
    105107    const MRawRunHeader *runheader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
    106108    if (!runheader)
    107109        *fLog << warn << dbginf << "Warning - cannot check file type, MRawRunHeader not found." << endl;
    108110    else
    109         if (runheader->GetRunType() == kRTMonteCarlo)
    110         {
    111             fPedestals=NULL;
     111        if (runheader->IsMonteCarloRun())
    112112            return kTRUE;
    113         }
    114113
    115114    fPedestals = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");

  • trunk/MagicSoft/Mars/manalysis/MCerPhotAnal2.cc

    r2946 r3374  
    106106        return kFALSE;
    107107
     108    fPedestals=NULL;
     109
    108110    const MRawRunHeader *runheader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
    109111    if (!runheader)
    110112        *fLog << warn << dbginf << "Warning - cannot check file type, MRawRunHeader not found." << endl;
    111113    else
    112         if (runheader->GetRunType() == kRTMonteCarlo)
    113         {
    114             fPedestals=NULL;
     114        if (runheader->IsMonteCarloRun())
    115115            return kTRUE;
    116         }
    117116
    118117    fPedestals = (MPedPhotCam*)pList->FindCreateObj("MPedPhotCam");

  • trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc

    r2942 r3374  
    154154    }
    155155
    156     Bool_t fIsMcFile = runheader->GetRunType() == kRTMonteCarlo;
    157     if (!fIsMcFile)
     156    if (!runheader->IsMonteCarloRun())
    158157        return kTRUE;
    159158

  • trunk/MagicSoft/Mars/manalysis/MMcCalibrationUpdate.cc

    r3344 r3374  
    9999    }
    100100
    101     return  run->GetRunType() == kRTMonteCarlo;
     101    return run->IsMonteCarloRun();
    102102}
    103103
    … …  
    108108Int_t MMcCalibrationUpdate::PreProcess(MParList *pList)
    109109{
    110   fCalCam = (MCalibrationChargeCam*) pList->FindObject(AddSerialNumber("MCalibrationChargeCam"));
    111   if ( !fCalCam )
    112     {
    113       *fLog << inf << dbginf << AddSerialNumber("MCalibrationChargeCam") << " does not exist... Creating." << endl;
    114 
    115       fCalCam = (MCalibrationChargeCam*) pList->FindCreateObj(AddSerialNumber("MCalibrationChargeCam"));
    116       if ( !fCalCam )
    117         {
    118           *fLog << err << dbginf << "Cannot create " << AddSerialNumber("MCalibrationChargeCam") << "... aborting." << endl;
    119           return kFALSE;
    120         }
    121     }
    122   else
    123     {
    124       fFillCalibrationCam = kFALSE;
    125       *fLog << inf << AddSerialNumber("MCalibrationChargeCam") << " already exists... " << endl;
    126     }
    127 
    128   fPedPhotCam = (MPedPhotCam*) pList->FindCreateObj(AddSerialNumber("MPedPhotCam"));
    129   if ( ! fPedPhotCam)
    130     {
    131       *fLog << err << dbginf << "Cannot create " << AddSerialNumber("MPedPhotCam") << "... aborting." << endl;
    132       return kFALSE;
    133     }
    134 
    135   fSignalCam = (MExtractedSignalCam*) pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
    136   if ( ! fSignalCam)
    137     {
    138       *fLog << err << dbginf << "Cannot find " << AddSerialNumber("MExtractedSignalCam") << "... aborting." << endl;
    139       return kFALSE;
    140     }
    141 
    142   return kTRUE;
    143 
     110    fCalCam = (MCalibrationChargeCam*) pList->FindObject(AddSerialNumber("MCalibrationChargeCam"));
     111    if (!fCalCam)
     112    {
     113        fCalCam = (MCalibrationChargeCam*) pList->FindCreateObj(AddSerialNumber("MCalibrationChargeCam"));
     114        if (!fCalCam)
     115            return kFALSE;
     116    }
     117    else
     118    {
     119        fFillCalibrationCam = kFALSE;
     120        *fLog << inf << AddSerialNumber("MCalibrationChargeCam") << " already exists... " << endl;
     121    }
     122
     123    fPedPhotCam = (MPedPhotCam*) pList->FindCreateObj(AddSerialNumber("MPedPhotCam"));
     124    if (!fPedPhotCam)
     125        return kFALSE;
     126
     127    fSignalCam = (MExtractedSignalCam*) pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
     128    if (!fSignalCam)
     129    {
     130        *fLog << err << AddSerialNumber("MExtractedSignalCam") << " not found... aborting." << endl;
     131        return kFALSE;
     132    }
     133
     134    return kTRUE;
    144135}
    145136
    … …  
    156147    //
    157148    if (!CheckRunType(pList))
    158       {
    159         *fLog << inf << "This is no MC file... skipping." << endl;
     149    {
     150        *fLog << inf << "This is no MC file... skipping." << endl;
    160151        return kTRUE;
    161       }
     152    }
    162153       
    163154    //
    164155    // Now check the existence of all necessary containers.
    165156    //
    166 
    167157    fGeom = (MGeomCam*) pList->FindObject(AddSerialNumber("MGeomCam"));
    168     if ( ! fGeom )
    169       {
    170         *fLog << err << dbginf << "Cannot find " << AddSerialNumber("MGeomCam") << "... aborting." << endl;
    171         return kFALSE;
    172       }
     158    if (!fGeom)
     159    {
     160        *fLog << err << AddSerialNumber("MGeomCam") << " not found... aborting." << endl;
     161        return kFALSE;
     162    }
    173163
    174164    fHeaderFadc = (MMcFadcHeader*)pList->FindObject(AddSerialNumber("MMcFadcHeader"));
    175165    if (!fHeaderFadc)
    176166    {
    177       *fLog << err << dbginf << AddSerialNumber("MMcFadcHeader") << " not found... aborting." << endl;
     167        *fLog << err << AddSerialNumber("MMcFadcHeader") << " not found... aborting." << endl;
    178168        return kFALSE;
    179169    }
    … …  
    182172    // Initialize Fadc simulation parameters:
    183173    //
    184     if ( fAmplitude < 0. )
    185       {
     174    if (fAmplitude < 0)
     175    {
    186176        fAmplitude = fHeaderFadc->GetAmplitud();
    187177        fAmplitudeOuter = fHeaderFadc->GetAmplitudOuter();
    188178        fConversionHiLo = fHeaderFadc->GetLow2HighGain();
    189       }
     179    }
    190180    else   // Check that following files have all the same FADC parameters
    191       {
    192         if ( fabs(fHeaderFadc->GetAmplitud()-fAmplitude) > 1.e-6  ||
     181    {
     182        if ( fabs(fHeaderFadc->GetAmplitud()-fAmplitude)           > 1.e-6  ||
    193183             fabs(fHeaderFadc->GetAmplitudOuter()-fAmplitudeOuter) > 1.e-6  ||
    194              fabs(fConversionHiLo-fHeaderFadc->GetLow2HighGain()) > 1.e-6 )
     184             fabs(fConversionHiLo-fHeaderFadc->GetLow2HighGain())  > 1.e-6 )
    195185          {
    196             *fLog << err << endl << endl << dbginf << "Parameters of MMcFadcHeader are not the same for all the read files. Aborting..." << endl << endl;
    197             return kFALSE;
    198           }
    199       }
     186              *fLog << err << "Parameters of MMcFadcHeader are not the same for all files... aborting." << endl;
     187            return kFALSE;
     188          }
     189    }
    200190
    201191    //
    … …  
    205195    // else has to be done in ReInit:
    206196    //
    207     if ( !fFillCalibrationCam )
    208       return kTRUE;
     197    if (!fFillCalibrationCam)
     198        return kTRUE;
    209199
    210200    //
    … …  
    278268        // counts for the RMS per slice:
    279269        //
    280 
    281         const Float_t pedestrms  = sigpix.IsLoGainUsed()?
    282           sqrt((Double_t)(fSignalCam->GetNumUsedLoGainFADCSlices())) *
    283           (fHeaderFadc->GetPedestalRmsLow(i)>0.? fHeaderFadc->GetPedestalRmsLow(i): 0.01)
    284           :
    285           sqrt((Double_t)(fSignalCam->GetNumUsedHiGainFADCSlices())) *
    286           (fHeaderFadc->GetPedestalRmsHigh(i)>0.? fHeaderFadc->GetPedestalRmsHigh(i) : 0.01);
     270        const Double_t used = (Double_t)(sigpix.IsLoGainUsed() ?
     271                                         fSignalCam->GetNumUsedLoGainFADCSlices() :
     272                                         fSignalCam->GetNumUsedHiGainFADCSlices());
     273
     274        const Float_t rms0 = sigpix.IsLoGainUsed() ?
     275            fHeaderFadc->GetPedestalRmsLow(i) :
     276            fHeaderFadc->GetPedestalRmsHigh(i);
     277
     278        const Float_t pedestrms = TMath::Sqrt(used) * (rms0>0 ? rms0 : 0.01);
    287279
    288280        //
    … …  
    297289
    298290        if (sigpix.IsLoGainUsed())
    299           pedpix.Set(adc2phot*hi2lo*pedestmean,
    300                      adc2phot*hi2lo*pedestrms);
     291            pedpix.Set(adc2phot*hi2lo*pedestmean, adc2phot*hi2lo*pedestrms);
    301292        else
    302           pedpix.Set(adc2phot*pedestmean, adc2phot*pedestrms);
     293            pedpix.Set(adc2phot*pedestmean, adc2phot*pedestrms);
    303294
    304295    }

  • trunk/MagicSoft/Mars/manalysis/MMcPedestalCopy.cc

    r3004 r3374  
    8686    }
    8787
    88     return run->GetRunType() == kRTMonteCarlo;
     88    return run->IsMonteCarloRun();
    8989}
    9090

  • trunk/MagicSoft/Mars/manalysis/MMcPedestalNSBAdd.cc

    r2470 r3374  
    106106    }
    107107
    108     return runheader->GetRunType() == kRTMonteCarlo;
     108    return runheader->IsMonteCarloRun();
    109109}
    110110

  • trunk/MagicSoft/Mars/manalysis/MPedCalcPedRun.cc

    r3199 r3374  
    175175    }
    176176    else
    177         if (runheader->GetRunType() == kRTMonteCarlo)
     177        if (runheader->IsMonteCarloRun())
    178178            return kTRUE;
    179179

  • trunk/MagicSoft/Mars/manalysis/MPedPhotCalc.cc

    r2850 r3374  
    109109  const MRawRunHeader *runheader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
    110110  if (!runheader)
    111     {
     111  {
    112112      *fLog << warn << dbginf;
    113113      *fLog << "Warning - cannot check file type, MRawRunHeader not found." << endl;
    114     }
     114  }
    115115  else
    116     if (runheader->GetRunType() == kRTMonteCarlo)
    117       return kTRUE;
     116      if (runheader->IsMonteCarloRun())
     117          return kTRUE;
    118118 
    119119
    120120  // Initialize arrays
    121121  if(fSumx.GetSize()==0)
    122     {
    123         const UShort_t num = fPedestals->GetSize();
     122  {
     123      const UShort_t num = fPedestals->GetSize();
    124124
    125         fSumx.Set(num);
    126         fSumx2.Set(num);
     125      fSumx.Set(num);
     126      fSumx2.Set(num);
    127127
    128         memset(fSumx.GetArray(),  0, sizeof(Float_t)*num);
    129         memset(fSumx2.GetArray(), 0, sizeof(Float_t)*num);
    130     }
     128      memset(fSumx.GetArray(),  0, sizeof(Float_t)*num);
     129      memset(fSumx2.GetArray(), 0, sizeof(Float_t)*num);
     130  }
    131131
    132     return kTRUE;
     132  return kTRUE;
    133133}
    134134
    … …  
    140140Int_t MPedPhotCalc::Process()
    141141{
    142   for(UInt_t i=0;i<fCerPhot->GetNumPixels();i++)
     142    for(UInt_t i=0;i<fCerPhot->GetNumPixels();i++)
    143143    {
    144       const MCerPhotPix &pix = (*fCerPhot)[i];
     144        const MCerPhotPix &pix = (*fCerPhot)[i];
    145145
    146       const Float_t nphot = pix.GetNumPhotons();
    147       const Int_t idx     = pix.GetPixId();
    148      
    149       fSumx[idx]  += nphot;
    150       fSumx2[idx] += nphot*nphot;
     146        const Float_t nphot = pix.GetNumPhotons();
     147        const Int_t idx     = pix.GetPixId();
     148
     149        fSumx[idx]  += nphot;
     150        fSumx2[idx] += nphot*nphot;
    151151    }
    152  
    153   fPedestals->SetReadyToSave();
    154  
    155   return kTRUE;
     152
     153    fPedestals->SetReadyToSave();
     154
     155    return kTRUE;
    156156}
    157157
    … …  
    161161//
    162162Int_t MPedPhotCalc::PostProcess()
    163   {
     163{
    164164    // Compute pedestals and rms from fSumx and fSumx2 arrays
    165165    const Int_t n    = GetNumExecutions();
    … …  
    167167
    168168    for(Int_t i=0; i<npix; i++)
    169       {
     169    {
    170170        const Float_t sum  = fSumx[i];
    171171        const Float_t sum2 = fSumx2[i];
    172        
     172
    173173        const Float_t photped = sum/n;
    174174        const Float_t photrms = TMath::Sqrt((sum2-sum*sum/n)/(n-1.));
    175175
    176         (*fPedestals)[i].Set(photped,photrms); 
    177       }
    178    
     176        (*fPedestals)[i].Set(photped,photrms);
     177    }
     178
    179179    return kTRUE;
    180180}

  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

    r3353 r3374  
    164164    if (!fRawEvt)
    165165    {
    166       *fLog << err << dbginf << "MRawEvtData not found... aborting." << endl;
     166      *fLog << err << "MRawEvtData not found... aborting." << endl;
    167167      return kFALSE;
    168168    }
    … …  
    170170    const MRawRunHeader *runheader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
    171171    if (!runheader)
    172       *fLog << warn << dbginf << "Warning - cannot check file type, MRawRunHeader not found." << endl;
     172      *fLog << warn << "Warning - cannot check file type, MRawRunHeader not found." << endl;
    173173    else
    174       if (runheader->GetRunType() == kRTMonteCarlo)
    175         {
     174      if (runheader->IsMonteCarloRun())
    176175          return kTRUE;
    177         }
    178176   
    179177    fCam = (MCalibrationChargeCam*)pList->FindCreateObj("MCalibrationChargeCam");
    180178    if (!fCam)
    181       {
    182         *fLog << err << dbginf << "MCalibrationChargeCam could not be created ... aborting." << endl;       
    183         return kFALSE;
    184       }
     179        return kFALSE;
    185180
    186181    fPINDiode = (MCalibrationChargePINDiode*)pList->FindCreateObj("MCalibrationChargePINDiode");
    187182    if (!fPINDiode)
    188       {
    189         *fLog << err << dbginf << "MCalibrationChargePINDiode could not be created ... aborting." << endl;       
    190         return kFALSE;
    191       }
     183        return kFALSE;
    192184
    193185    fCam->SetPINDiode(fPINDiode);
    … …  
    195187    fBlindPixel = (MCalibrationChargeBlindPix*)pList->FindCreateObj("MCalibrationChargeBlindPix");
    196188    if (!fBlindPixel)
    197       {
    198         *fLog << err << dbginf << "MCalibrationChargeBlindPix could not be created ... aborting." << endl;       
    199         return kFALSE;
    200       }
     189        return kFALSE;
    201190
    202191    fCam->SetBlindPixel(fBlindPixel);
    203192
    204     fEvtTime      = (MTime*)pList->FindObject("MTime");
     193    fEvtTime = (MTime*)pList->FindObject("MTime");
    205194
    206195    fPedestals = (MPedestalCam*)pList->FindObject("MPedestalCam");
    207196    if (!fPedestals)
    208       {
    209         *fLog << err << dbginf << "Cannot find MPedestalCam ... aborting" << endl;
    210         return kFALSE;
    211       }
    212 
     197    {
     198        *fLog << err << "MPedestalCam not found... aborting" << endl;
     199        return kFALSE;
     200    }
    213201
    214202    fSignals = (MExtractedSignalCam*)pList->FindObject("MExtractedSignalCam");
    215203    if (!fSignals)
    216       {
    217         *fLog << err << dbginf << "Cannot find MExtractedSignalCam ... aborting" << endl;
    218         return kFALSE;
    219       }
     204    {
     205        *fLog << err << "MExtractedSignalCam not found... aborting" << endl;
     206        return kFALSE;
     207    }
    220208
    221209    return kTRUE;
    … …  
    230218Bool_t MCalibrationChargeCalc::ReInit(MParList *pList )
    231219{
    232  
    233220    fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
    234221    if (!fRunHeader)
    235222    {
    236       *fLog << err << dbginf << ": MRawRunHeader not found... aborting." << endl;
     223      *fLog << err << "MRawRunHeader not found... aborting." << endl;
    237224      return kFALSE;
    238225    }
    239 
    240226
    241227    fGeom = (MGeomCam*)pList->FindObject("MGeomCam");
    242228    if (!fGeom)
    243229    {
    244       *fLog << err << GetDescriptor() << ": No MGeomCam found... aborting." << endl;
    245       return kFALSE;
     230        *fLog << err << "No MGeomCam found... aborting." << endl;
     231        return kFALSE;
    246232    }
    247233
    … …  
    255241
    256242    for (UInt_t i=0; i<npixels; i++)
    257       {
    258        
     243    {
    259244        MCalibrationChargePix &pix = (*fCam)[i];
    260245        pix.DefinePixId(i);
    … …  
    264249        pix.SetAbsTimeBordersLoGain(fSignals->GetFirstUsedSliceLoGain(),
    265250                                    fSignals->GetLastUsedSliceLoGain());
    266        
    267      }
     251    }
    268252   
     253    if (fExcludedPixelsFile.IsNull())
     254        return kTRUE;
     255
    269256    //
    270257    // Look for file to exclude pixels from analysis
    271258    //
    272     if (!fExcludedPixelsFile.IsNull())
    273       {
    274        
    275         fExcludedPixelsFile = gSystem->ExpandPathName(fExcludedPixelsFile.Data());
    276        
     259    gSystem->ExpandPathName(fExcludedPixelsFile);
     260
     261    //
     262    // Initialize reading the file
     263    //
     264    ifstream in(fExcludedPixelsFile);
     265    if (!in)
     266    {
     267        *fLog << warn << "Cannot open file '" << fExcludedPixelsFile << "'" << endl;
     268        return kTRUE;
     269    }
     270
     271    *fLog << inf << "Use excluded pixels from file: '" << fExcludedPixelsFile << "'" << endl;
     272
     273    //
     274    // Read the file and count the number of entries
     275    //
     276    UInt_t pixel = 0;
     277
     278    while (++fNumExcludedPixels)
     279    {
     280        in >> pixel;
     281        if (!in)
     282            break;
     283
    277284        //
    278         // Initialize reading the file
     285        // Check for out of range
    279286        //
    280         ifstream in(fExcludedPixelsFile.Data(),ios::in);
    281 
    282         if (in)
    283           {
    284             *fLog << inf << "Use excluded pixels from file: '" << fExcludedPixelsFile.Data() << "'" << endl;
    285             //
    286             // Read the file and count the number of entries
    287             //
    288             UInt_t pixel = 0;
    289            
    290             while (++fNumExcludedPixels)
    291               {
    292                
    293                 in >> pixel;
    294 
    295                 if (!in.good())
    296                   break;
    297                 //
    298                 // Check for out of range
    299                 //
    300                 if (pixel > npixels)
    301                   {
    302                     *fLog << warn << "WARNING: To be excluded pixel: " << pixel
    303                           << " is out of range " << endl;
    304                     continue;
    305                   }
    306                 //
    307                 // Exclude pixel
    308                 //
    309                 MCalibrationChargePix &pix = (*fCam)[pixel];
    310                 pix.SetExcluded();
    311                
    312                 *fLog << GetDescriptor() << inf << ": Exclude Pixel: " << pixel << endl;
    313                
    314               }
    315            
    316             if (--fNumExcludedPixels == 0)
    317               *fLog << warn << "WARNING: File '" << fExcludedPixelsFile.Data()
    318                     << "'" << " is empty " << endl;
    319             else
    320               fCam->SetNumPixelsExcluded(fNumExcludedPixels);
    321            
    322           }
    323         else
    324           *fLog << warn << dbginf << "Cannot open file '" << fExcludedPixelsFile.Data() << "'" << endl;
    325       }
    326    
     287        if (pixel > npixels)
     288        {
     289            *fLog << warn << "WARNING: To be excluded pixel: " << pixel << " is out of range " << endl;
     290            continue;
     291        }
     292        //
     293        // Exclude pixel
     294        //
     295        MCalibrationChargePix &pix = (*fCam)[pixel];
     296        pix.SetExcluded();
     297
     298        *fLog << GetDescriptor() << inf << ": Exclude Pixel: " << pixel << endl;
     299
     300    }
     301
     302    if (--fNumExcludedPixels == 0)
     303        *fLog << warn << "WARNING: File '" << fExcludedPixelsFile << "'" << " is empty " << endl;
     304    else
     305        fCam->SetNumPixelsExcluded(fNumExcludedPixels);
     306
    327307    return kTRUE;
    328308}
    … …  
    513493  return kTRUE;
    514494}
    515 
    516 
    517 
    518 
    519 
    520 
    521 

  • trunk/MagicSoft/Mars/mcalib/MMcCalibrationCalc.cc

    r3265 r3374  
    8787    }
    8888
    89     return  run->GetRunType() == kRTMonteCarlo;
     89    return  run->IsMonteCarloRun();
    9090}
    9191

  • trunk/MagicSoft/Mars/mpointing/MPointingPosCalc.cc

    r2601 r3374  
    8080    switch (fRunType)
    8181    {
    82     case kRTData:
     82    case MRawRunHeader::kRTData:
    8383        fReport = (MReportDrive*)plist->FindObject("MReportDrive");
    8484        if (!fReport)
    … …  
    8989        return kTRUE;
    9090
    91     case kRTMonteCarlo:
     91    case MRawRunHeader::kRTMonteCarlo:
    9292        fMcEvt = (MMcEvt*)plist->FindObject("MMcEvt");
    9393        if (!fMcEvt)
    … …  
    9898        return kTRUE;
    9999
    100     case kRTPedestal:
     100    case MRawRunHeader::kRTPedestal:
    101101        *fLog << err << "Cannot work in a pedestal Run!... aborting." << endl;
    102102        return kFALSE;
    103103
    104     case kRTCalibration:
     104    case MRawRunHeader::kRTCalibration:
    105105        *fLog << err << "Cannot work in a calibration Run!... aborting." << endl;
    106106        return kFALSE;
    … …  
    132132    switch (fRunType)
    133133    {
    134     case kRTData:
     134    case MRawRunHeader::kRTData:
    135135        fPosition->SetLocalPosition(fReport->GetNominalZd(), fReport->GetNominalAz());
    136136        return kTRUE;
    137137
    138     case kRTMonteCarlo:
     138    case MRawRunHeader::kRTMonteCarlo:
    139139        fPosition->SetLocalPosition(fMcEvt->GetTelescopeTheta()*TMath::RadToDeg(), fMcEvt->GetTelescopePhi()*TMath::RadToDeg());
    140140        return kTRUE;

  • trunk/MagicSoft/Mars/mraw/MRawEvtHeader.h

    r3338 r3374  
    1414class MRawRunHeader;
    1515
    16 //
    17 // Trigger Type (TT)
    18 //
    19 enum {
    20     kTTEvent       = 0,
    21     kTTPedestal    = 1,
    22     kTTCalibration = 2
    23 };
    24 
    2516class MRawEvtHeader : public MParContainer
    2617{
     18public:
     19    //
     20    // Trigger Type (TT)
     21    //
     22    enum {
     23        kTTEvent       = 0,
     24        kTTPedestal    = 1,
     25        kTTCalibration = 2
     26    };
     27
    2728private:
    2829    MTime   *fTime;            //! object to store the time in (ReadEvt)

  • trunk/MagicSoft/Mars/mraw/MRawFileWrite.cc

    r3183 r3374  
    247247    switch (type)
    248248    {
    249     case kTTEvent:
     249    case MRawEvtHeader::kTTEvent:
    250250        fTData->Fill();
    251251        return kTRUE;
    252252
    253     case kTTPedestal:
     253    case MRawEvtHeader::kTTPedestal:
    254254        fTPedestal->Fill();
    255255        return kTRUE;
    256256
    257     case kTTCalibration:
     257    case MRawEvtHeader::kTTCalibration:
    258258        fTCalibration->Fill();
    259259        return kTRUE;

  • trunk/MagicSoft/Mars/mraw/MRawRunHeader.cc

    r3183 r3374  
    6161using namespace std;
    6262
     63const UShort_t MRawRunHeader::kMagicNumber      = 0xc0c0;
     64const Byte_t   MRawRunHeader::kMaxFormatVersion =      3;
     65
    6366// --------------------------------------------------------------------------
    6467//

  • trunk/MagicSoft/Mars/mraw/MRawRunHeader.h

    r3183 r3374  
    1414class MArrayS;
    1515
    16 //
    17 // Magic number to detect the magic file type
    18 //
    19 const UShort_t kMagicNumber      = 0xc0c0;
    20 const Byte_t   kMaxFormatVersion =      3;
    21 
    22 //
    23 // enum for the Run Type. Monte Carlo Runs have
    24 // to have a value greater than 255 (>0xff)
    25 //
    26 enum {
    27     kRTData        = 0x0000,
    28     kRTPedestal    = 0x0001,
    29     kRTCalibration = 0x0002,
    30     kRTMonteCarlo  = 0x0100,
    31     kRTNone        = 0xffff
    32 };
    33 
    3416class MRawRunHeader : public MParContainer
    3517{
     18public:
     19    //
     20    // enum for the Run Type. Monte Carlo Runs have
     21    // to have a value greater than 255 (>0xff)
     22    //
     23    enum {
     24        kRTData        = 0x0000,
     25        kRTPedestal    = 0x0001,
     26        kRTCalibration = 0x0002,
     27        kRTMonteCarlo  = 0x0100,
     28        kRTNone        = 0xffff
     29    };
     30
     31    //
     32    // Magic number to detect the magic file type
     33    //
     34    static const UShort_t kMagicNumber;
     35    static const Byte_t   kMaxFormatVersion;
     36
    3637private:
    3738    /* ---- Run Header Informations ---- */
    … …  
    105106    Int_t GetNumTotalBytes() const;
    106107
     108    Bool_t IsMonteCarloRun() const { return fRunType>0x00ff; }
     109
    107110    void Print(Option_t *t=NULL) const;
    108111
Note: See TracChangeset for help on using the changeset viewer.