Index: trunk/MagicSoft/Mars/msimcamera/MSimAPD.cc =================================================================== --- trunk/MagicSoft/Mars/msimcamera/MSimAPD.cc (revision 9240) +++ trunk/MagicSoft/Mars/msimcamera/MSimAPD.cc (revision 9240) @@ -0,0 +1,191 @@ +/* ======================================================================== *\ +! +! * +! * This file is part of MARS, the MAGIC Analysis and Reconstruction +! * Software. It is distributed to you in the hope that it can be a useful +! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. +! * It is distributed WITHOUT ANY WARRANTY. +! * +! * Permission to use, copy, modify and distribute this software and its +! * documentation for any purpose is hereby granted without fee, +! * provided that the above copyright notice appear in all copies and +! * that both that copyright notice and this permission notice appear +! * in supporting documentation. It is provided "as is" without express +! * or implied warranty. +! * +! +! +! Author(s): Thomas Bretz 1/2009 +! +! Copyright: Software Development, 2000-2009 +! +! +\* ======================================================================== */ + +////////////////////////////////////////////////////////////////////////////// +// +// MSimAPD +// +// This tasks simulates the individual APDs. Before starting the APD is +// initialized randomly according to the photon rate hitting the APD. Such +// it is assumed that the initial condition of the APD is similar to the real +// one. In this context it is assumed that the events are independent, so +// that the APD is always in the same condition. +// +// For every photon and event the behaviour of the APD is simulated. The +// output is set as weight to the MPhotonData containers. +// +// Remark: +// - The photons MUST be sorted increasing in time. +// - The photon rate used to initialize the APD must match the one used +// to "fill" the random photons. (FIXME: This should be stored somewhere) +// +// Input Containers: +// fNameGeomCam [MGeomCam] +// MPhotonEvent +// MPhotonStatistics +// +// Output Containers: +// MPhotonEvent +// +////////////////////////////////////////////////////////////////////////////// +#include "MSimAPD.h" + +#include +#include + +#include "MLog.h" +#include "MLogManip.h" + +#include "MMath.h" +#include "MParList.h" + +#include "MGeomCam.h" + +#include "MPhotonEvent.h" +#include "MPhotonData.h" + +#include "MAvalanchePhotoDiode.h" + +ClassImp(MSimAPD); + +using namespace std; + +// -------------------------------------------------------------------------- +// +// Default Constructor. +// +MSimAPD::MSimAPD(const char* name, const char *title) +: fGeom(0), fEvt(0), fStat(0) +{ + fName = name ? name : "MSimAPD"; + fTitle = title ? title : " Task to simulate the detection behaviour of APDs"; +} + +// -------------------------------------------------------------------------- +// +// Get the necessary parameter containers +// +Int_t MSimAPD::PreProcess(MParList *pList) +{ + fGeom = (MGeomCam*)pList->FindObject(fNameGeomCam, "MGeomCam"); + if (!fGeom) + { + *fLog << inf << fNameGeomCam << " [MGeomCam] not found..." << endl; + + fGeom = (MGeomCam*)pList->FindCreateObj(fNameGeomCam); + if (!fGeom) + return kFALSE; + } + + fStat = (MPhotonStatistics*)pList->FindObject("MPhotonStatistics"); + if (!fStat) + { + *fLog << err << "MPhotonStatistics not found... aborting." << endl; + return kFALSE; + } + + fEvt = (MPhotonEvent*)pList->FindObject("MPhotonEvent"); + if (!fEvt) + { + *fLog << err << "MPhotonEvent not found... aborting." << endl; + return kFALSE; + } + + return kTRUE; +} + +// -------------------------------------------------------------------------- +// +// Initialize as many APDs as we have pixels in the fGeomCam +// +Bool_t MSimAPD::ReInit(MParList *plist) +{ + if (UInt_t(fAPDs.GetEntriesFast())!=fGeom->GetNumPixels()) + { + fAPDs.Delete(); + + // FIXME: + // * initialize an empty APD and read the APD setup from a file to + // allow different APDs. + // * Make the arguments a data member of MSimAPD + + for (UInt_t i=0; iGetNumPixels(); i++) + //fAPDs.Add(new APD(60, 0.2, 3, 8.75)); + //fAPDs.Add(new APD(60/2, 0.2, 3e-9, 8.75e-9*4)); + fAPDs.Add(new APD(60/2, 0.2, 3, 8.75*4)); + } + + return kTRUE; +} + +// -------------------------------------------------------------------------- +// +// Process all photons through the corresponding APD and set the output +// (weight) accordingly. +// +Int_t MSimAPD::Process() +{ + //const Double_t rate = 40e9; + // FIXME: Where do we get this number from?? + // const Double_t rate = 0.04; + + // Make all APDs look neutral for the first hit by a photon according to the + // average hit rate + const UInt_t npix = fAPDs.GetEntriesFast(); + for (UInt_t idx=0; idx(fAPDs.UncheckedAt(idx))->FillRandom(fFreq, fStat->GetTimeFirst()); + + // Get number of photons + const Int_t num = fEvt->GetNumPhotons(); + + // Loop over all photons + for (Int_t i=0; i(fAPDs.UncheckedAt(idx))->HitRandomCell(t); + + // FIXME: Make a proper simulation of the excess noise!!! + //const Double_t signal = gRandom->Gaus(hits, 0.2*TMath::Sqrt(hits)); + + // Set the weight to the input + ph.SetWeight(hits); + } + + return kTRUE; +} Index: trunk/MagicSoft/Mars/msimcamera/MSimAPD.h =================================================================== --- trunk/MagicSoft/Mars/msimcamera/MSimAPD.h (revision 9240) +++ trunk/MagicSoft/Mars/msimcamera/MSimAPD.h (revision 9240) @@ -0,0 +1,43 @@ +#ifndef MARS_MSimAPD +#define MARS_MSimAPD + +#ifndef MARS_MTask +#include "MTask.h" +#endif + +#ifndef ROOT_TObjArray +#include +#endif + +class MGeomCam; +class MParList; +class MPhotonEvent; +class MPhotonStatistics; + +class MSimAPD : public MTask +{ +private: + MGeomCam *fGeom; //! APD geometry (used to know how many pixels we have) + MPhotonEvent *fEvt; //! Event storing the photon information + MPhotonStatistics *fStat; //! Storing event statistics (needed for the start-time) + + TObjArray fAPDs; //! Array keeping the necessary number of APDs + + TString fNameGeomCam; // Name of the geometry container storing the APD gemeotry + + Double_t fFreq; // Frequency of random phtons which hit the APDs + + // MTask + Int_t PreProcess(MParList *pList); + Bool_t ReInit(MParList *pList); + Int_t Process(); + +public: + MSimAPD(const char *name=NULL, const char *title=NULL); + + void SetNameGeomCam(const char *name="MGeomCam") { fNameGeomCam = name; } + void SetFreq(Float_t f) { fFreq=f; } + + ClassDef(MSimAPD, 0) // Task to simulate the detection behaviour of APDs +}; +#endif