Ignore:
Timestamp:
11/04/03 17:08:27 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/manalysis
Files:
6 edited

Legend:

Unmodified
Added
Removed

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

    r2377 r2470  
    5959//
    6060//  Input Containers:
    61 //   MCerPhotEvt[, MBlindPixels]
     61//   MCerPhotEvt
     62//   [MBlindPixels]
    6263//
    6364//  Output Containers:
    64 //   MCerPhotEvt, MBlindPixels
     65//   MCerPhotEvt
     66//   MBlindPixels
    6567//
    6668/////////////////////////////////////////////////////////////////////////////
    … …  
    112114{
    113115    if (TESTBIT(fFlags, kUseBlindPixels))
    114         fPixels = (MBlindPixels*)pList->FindObject("MBlindPixels");
     116        fPixels = (MBlindPixels*)pList->FindObject(AddSerialNumber("MBlindPixels"));
    115117    else
    116         fPixels = (MBlindPixels*)pList->FindCreateObj("MBlindPixels");
     118        fPixels = (MBlindPixels*)pList->FindCreateObj(AddSerialNumber("MBlindPixels"));
    117119    if (!fPixels)
    118120        return kFALSE;
    119121   
    120     fEvt = (MCerPhotEvt*)pList->FindObject("MCerPhotEvt");
     122    fEvt = (MCerPhotEvt*)pList->FindObject(AddSerialNumber("MCerPhotEvt"));
    121123    if (!fEvt)
    122124    {
    … …  
    124126        return kFALSE;
    125127    }
    126     fPed = (MPedestalCam*)pList->FindObject("MPedestalCam");
     128    fPed = (MPedestalCam*)pList->FindObject(AddSerialNumber("MPedestalCam"));
    127129    if (!fEvt)
    128130    {
    … …  
    131133    }
    132134
    133     fGeomCam = (MGeomCam*)pList->FindObject("MGeomCam");
     135    fGeomCam = (MGeomCam*)pList->FindObject(AddSerialNumber("MGeomCam"));
    134136    if (!fGeomCam)
    135137        *fLog << warn << dbginf << "No camera geometry available... can't use interpolation." << endl;

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

    r2424 r2470  
    3535//
    3636//  Input Containers:
    37 //   MRawRunHeader, MRawEvtData, MPedestalCam
     37//   MRawEvtData
     38//   MPedestalCam
     39//   [MRawRunHeader]
    3840//
    3941//  Output Containers:
    … …  
    99101    }
    100102
    101     fRawEvt = (MRawEvtData*)pList->FindObject("MRawEvtData");
     103    fRawEvt = (MRawEvtData*)pList->FindObject(AddSerialNumber("MRawEvtData"));
    102104    if (!fRawEvt)
    103105    {
    … …  
    106108    }
    107109
    108     fPedestals = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
     110    fPedestals = (MPedestalCam*)pList->FindCreateObj(AddSerialNumber("MPedestalCam"));
    109111    if (!fPedestals)
    110112    {
    … …  
    113115    }
    114116
    115     fCerPhotEvt = (MCerPhotEvt*)pList->FindCreateObj("MCerPhotEvt");
     117    fCerPhotEvt = (MCerPhotEvt*)pList->FindCreateObj(AddSerialNumber("MCerPhotEvt"));
    116118    if (!fCerPhotEvt)
    117119        return kFALSE;

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

    r2438 r2470  
    4444//
    4545//  Output Containers:
    46 //   MPedestalCam, MBlindPixels
     46//   [MPedestalCam]
     47//   [MBlindPixels]
    4748//
    4849//////////////////////////////////////////////////////////////////////////////
    … …  
    7980Int_t MGeomApply::PreProcess(MParList *pList)
    8081{
    81     MGeomCam *cam = (MGeomCam*)pList->FindObject("MGeomCam");
     82    MGeomCam *cam = (MGeomCam*)pList->FindObject(AddSerialNumber("MGeomCam"));
    8283    if (cam)
    8384        return kTRUE;
    8485
    85     cam = (MGeomCam*)pList->FindCreateObj(fGeomName, "MGeomCam");
     86    cam = (MGeomCam*)pList->FindCreateObj(AddSerialNumber(fGeomName), "MGeomCam");
    8687
    8788    return cam ? kTRUE : kFALSE;
    … …  
    9798Bool_t MGeomApply::ReInit(MParList *pList)
    9899{
    99     MGeomCam *cam = (MGeomCam*)pList->FindObject("MGeomCam");
     100    MGeomCam *cam = (MGeomCam*)pList->FindObject(AddSerialNumber("MGeomCam"));
    100101    if (!cam)
    101102    {
    … …  
    104105    }
    105106
    106     MPedestalCam *ped = (MPedestalCam*)pList->FindObject("MPedestalCam");
     107    MPedestalCam *ped = (MPedestalCam*)pList->FindObject(AddSerialNumber("MPedestalCam"));
    107108    if (ped)
    108109        ped->InitSize(cam->GetNumPixels());
    109110
    110     MBlindPixels *bnd = (MBlindPixels*)pList->FindObject("MBlindPixels");
     111    MBlindPixels *bnd = (MBlindPixels*)pList->FindObject(AddSerialNumber("MBlindPixels"));
    111112    if (bnd)
    112113        bnd->InitSize(cam->GetNumPixels());

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

    r2454 r2470  
    3333//  Input Containers:
    3434//   MMcFadcHeader
     35//   [MMcRunHeader]
     36//   [MRawRunHeader]
    3537//
    3638//  Output Containers:
    … …  
    104106    // to be done only if this is a MC file.
    105107    //
    106     MPedestalCam *pedcam = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
     108    MMcFadcHeader *fadc = (MMcFadcHeader*)pList->FindObject(AddSerialNumber("MMcFadcHeader"));
     109    if (!fadc)
     110    {
     111        *fLog << err << dbginf << "MMcFadcHeader not found... aborting." << endl;
     112        return kFALSE;
     113    }
     114
     115    MPedestalCam *pedcam = (MPedestalCam*)pList->FindCreateObj(AddSerialNumber("MPedestalCam"));
    107116    if (!pedcam)
    108117        return kFALSE;
    … …  
    111120    if (!mcrun)
    112121        *fLog << warn << dbginf << "MMcRunHeader not found... assuming camera<0.7" << endl;
    113 
    114     MMcFadcHeader *fadc = (MMcFadcHeader*)pList->FindObject("MMcFadcHeader");
    115     if (!fadc)
    116     {
    117         *fLog << err << dbginf << "MMcFadcHeader not found... aborting." << endl;
    118         return kFALSE;
    119     }
    120122
    121123    const int num = pedcam->GetSize();

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

    r2454 r2470  
    2424
    2525/////////////////////////////////////////////////////////////////////////////
    26 //                                                                         //
    27 //  MMcPedestalNSBAdd                                                      //
    28 //  -----------------                                                      //
    29 //                                                                         //
    30 //  This Task adds the contribution of the diffuse NSB to the FADC         //
    31 //  pedestals. We assume that NSB introduces larger fluctuation but does   //
    32 //  not change the mean value.                                             //
    33 //  To be precise we add quadratically to the rms, which is already in     //
    34 //  MPedestalCam, the NSB contribution.                                    //
    35 //  The number of photons from the diffuse NSB follows a poisson           //
    36 //  distribution with expected mean value X, then its standard deviation   //
    37 //  is sqrt(X).                                                            //
    38 //  X is the number of phe per ns so we have to convert in FADC counts per //
    39 //  slice:                                                                 //
    40 //                                                                         //
    41 //  Y = sqrt(X * FADC_time / Number_of_slices * pixel_size)                //
    42 //      * Amp_single_phe_response                                          //
    43 //                                                                         //
    44 //  Input Containers:                                                      //
    45 //   MMcFadcHeader                                                         //
    46 //   MMcRunHeader                                                          //
    47 //   MRawRunHeader                                                         //
    48 //                                                                         //
    49 //  Output Containers:                                                     //
    50 //   MPedestalCam                                                          //
    51 //                                                                         //
     26//
     27//  MMcPedestalNSBAdd
     28//  -----------------
     29//
     30//  This Task adds the contribution of the diffuse NSB to the FADC
     31//  pedestals. We assume that NSB introduces larger fluctuation but does
     32//  not change the mean value.
     33//  To be precise we add quadratically to the rms, which is already in
     34//  MPedestalCam, the NSB contribution.
     35//  The number of photons from the diffuse NSB follows a poisson
     36//  distribution with expected mean value X, then its standard deviation
     37//  is sqrt(X).
     38//  X is the number of phe per ns so we have to convert in FADC counts per
     39//  slice:
     40//
     41//  Y = sqrt(X * FADC_time / Number_of_slices * pixel_size)
     42//      * Amp_single_phe_response
     43//
     44//  Input Containers:
     45//   MMcFadcHeader
     46//   MRawRunHeader
     47//   [MMcRunHeader]
     48//
     49//  Output Containers:
     50//   MPedestalCam
     51//
    5252/////////////////////////////////////////////////////////////////////////////
    5353#include "MMcPedestalNSBAdd.h"
    … …  
    186186    // is correct.
    187187    //
    188     MPedestalCam *pedcam = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
     188    MPedestalCam *pedcam = (MPedestalCam*)pList->FindCreateObj(AddSerialNumber("MPedestalCam"));
    189189    if (!pedcam)
    190190        return kFALSE;
    191191
    192     MMcFadcHeader *fadc = (MMcFadcHeader*)pList->FindObject("MMcFadcHeader");
     192    MMcFadcHeader *fadc = (MMcFadcHeader*)pList->FindObject(AddSerialNumber("MMcFadcHeader"));
    193193    if (!fadc)
    194194    {
    … …  
    197197    }
    198198
    199     MGeomCam *geom = (MGeomCam*)pList->FindObject("MGeomCam");
     199    MGeomCam *geom = (MGeomCam*)pList->FindObject(AddSerialNumber("MGeomCam"));
    200200    if (!geom)
    201201    {

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

    r2209 r2470  
    2525
    2626/////////////////////////////////////////////////////////////////////////////
    27 //                                                                         //
    28 //  MSigmabarCalc                                                          //
    29 //                                                                         //
    30 //  This task calculates Sigmabar using the MSigmabar container and stores //
    31 //  its extremal values together with the corresponding Theta values       //
    32 //  in MSigmabarParam. For the time being, Theta is taken from a Monte     //
    33 //  Carlo container. This is preliminary and to be changed ASAP.           //
    34 //                                                                         //
     27//
     28//  MSigmabarCalc
     29//
     30//  This task calculates Sigmabar using the MSigmabar container and stores
     31//  its extremal values together with the corresponding Theta values       
     32//  in MSigmabarParam. For the time being, Theta is taken from a Monte
     33//  Carlo container. This is preliminary and to be changed ASAP.
     34//
     35//  Input Containers:
     36//   MGeomCam
     37//   MPedestalCam
     38//   MRawRunHeader
     39//   MMcEvt  (FIXME: Must be replaced by a 'real-data' container)
     40//   MCerPhotEvt
     41//
     42//  Output Containers:
     43//   MSigmabar
     44//   MSigmabarParam
     45//
    3546/////////////////////////////////////////////////////////////////////////////
    3647#include "MSigmabarCalc.h"
    … …  
    7283Int_t MSigmabarCalc::PreProcess(MParList *pList)
    7384{
    74     fCam = (MGeomCam*)pList->FindObject("MGeomCam");
     85    fCam = (MGeomCam*)pList->FindObject(AddSerialNumber("MGeomCam"));
    7586    if (!fCam)
    7687    {
    … …  
    7990    }
    8091
    81     fPed = (MPedestalCam*)pList->FindObject("MPedestalCam");
     92    fPed = (MPedestalCam*)pList->FindObject(AddSerialNumber("MPedestalCam"));
    8293    if (!fPed)
    8394    {
    … …  
    94105
    95106    // This is needed for determining min/max Theta
    96     fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
     107    fMcEvt = (MMcEvt*)pList->FindObject(AddSerialNumber("MMcEvt"));
    97108    if (!fMcEvt)
    98109    {
    … …  
    102113    }
    103114
    104     fEvt = (MCerPhotEvt*)pList->FindObject("MCerPhotEvt");
     115    fEvt = (MCerPhotEvt*)pList->FindObject(AddSerialNumber("MCerPhotEvt"));
    105116    if (!fEvt)
    106117      {
    … …  
    109120      }
    110121
    111     fSig = (MSigmabar*)pList->FindCreateObj("MSigmabar");
     122    fSig = (MSigmabar*)pList->FindCreateObj(AddSerialNumber("MSigmabar"));
    112123    if (!fSig)
    113124        return kFALSE;
    114125
    115     fSigParam = (MSigmabarParam*)pList->FindCreateObj("MSigmabarParam");
     126    // FIXME: This might have to change in a
     127    //        'more-than-one-telescope-environment'
     128    fSigParam = (MSigmabarParam*)pList->FindCreateObj(AddSerialNumber("MSigmabarParam"));
    116129    if (!fSigParam)
    117130        return kFALSE;
Note: See TracChangeset for help on using the changeset viewer.