Changeset 4072 for trunk/MagicSoft/Mars/mtemp/mifae/library

Timestamp:

05/14/04 16:35:57 (21 years ago)

Author:

rico

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae/library

Files:

• IFAELinkDef.h (modified) (1 diff)
• MSrcPlace.cc (added)
• MSrcPlace.h (added)
• MSrcPosFromFile.cc (modified) (5 diffs)
• MSrcPosFromFile.h (modified) (3 diffs)
• MSrcRotate.cc (modified) (10 diffs)
• MSrcRotate.h (modified) (1 diff)
• Makefile (modified) (1 diff)

trunk/MagicSoft/Mars/mtemp/mifae/library/IFAELinkDef.h

@@ -8 +8 @@
 #pragma link C++ class MPSFFit+;
 #pragma link C++ class MPSFFitCalc+;
+#pragma link C++ class MSrcPlace+;
 #pragma link C++ class MSrcPosFromFile+;
 #pragma link C++ class MSrcRotate+;

trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcPosFromFile.cc

@@ -27 +27 @@
 // MSrcPosFromFile
 //
-// Task to calculate the position of the source as a function of run number
+// Task to set the position of the source as a function of run number according
+// to the positions read from an input file
 //
 //  Output Containers:
@@ -39 +40 @@
 
 #include "MParList.h"
-
 #include "MSrcPosFromFile.h"
 
@@ -53 +53 @@
 
 static const TString gsDefName  = "MSrcPosFromFile";
-static const TString gsDefTitle = "Calculate position of the (off axis) source";
+static const TString gsDefTitle = "Set the position of the (off axis) source according to input file";
 
 // -------------------------------------------------------------------------
 //
-// Default constructor
+// Default constructor. cardpath is the name of the input file (.pos) where the
+// source positions are stored in the format Run# x y (in mm). mode indicates whether
+// to read or to save the positions of the source in/from the internal histogram
 //
 MSrcPosFromFile::MSrcPosFromFile(TString cardpath, OnOffMode_t mode, const char *name, const char *title)
-    : fRawRunHeader(NULL), fSrcPos(NULL), fMode(mode), fSourcePositionFilePath(cardpath)
-{
-    fName  = name  ? name  : gsDefName.Data();
-    fTitle = title ? title : gsDefTitle.Data();
-
-    fLastRun = 0;
-
-// Count the number of runs in the card with the source poistions
-    ReadSourcePositionsFile(kCount);
-
-    fRunList = new Int_t[fNumRuns];
-    fRunSrcPos = new MSrcPosCam[fNumRuns];
-    fRunMap = new TExMap(fNumRuns);
-
-    Float_t cameraSize   = 600; //[mm]
-    Float_t binPrecision =  3;  //[mm] ~ 0.01 deg
-
-    UInt_t nbins =  (UInt_t)(cameraSize*2/binPrecision);
-
-    fHistSrcPos = new TH2F("HistSrcPos","",nbins,-cameraSize,cameraSize,nbins,-cameraSize,cameraSize);
-
-// Read card with the source poistions
-    ReadSourcePositionsFile(kRead);
+  : fRawRunHeader(NULL), fSourcePositionFilePath(cardpath)
+{
+  fName  = name  ? name  : gsDefName.Data();
+  fTitle = title ? title : gsDefTitle.Data();
+  SetMode(mode);
+
+  fLastRun = 0;
+  
+  // Count the number of runs in the card with the source poistions
+  ReadSourcePositionsFile(kCount);
+  
+  fRunList = new Int_t[fNumRuns];
+  fRunSrcPos = new MSrcPosCam[fNumRuns];
+  fRunMap = new TExMap(fNumRuns);
+  
+  // Read card with the source poistions
+  ReadSourcePositionsFile(kRead);
 }
 
 MSrcPosFromFile::~MSrcPosFromFile()
 {
-    delete [] fRunList;
-    delete [] fRunSrcPos;
-    delete fRunMap;
-    delete fHistSrcPos;
+  delete [] fRunList;
+  delete [] fRunSrcPos;
+  delete fRunMap;
 }
 
@@ -97 +92 @@
 Int_t MSrcPosFromFile::PreProcess(MParList *pList)
 {
-
-    fRawRunHeader = (MRawRunHeader*)pList->FindObject(AddSerialNumber("MRawRunHeader"));
-    if (!fRawRunHeader)
+  if(!SearchForSrcPos(pList))
+    return kFALSE;
+  
+  fRawRunHeader = (MRawRunHeader*)pList->FindObject(AddSerialNumber("MRawRunHeader"));
+  if (!fRawRunHeader)
     {
       *fLog << err << AddSerialNumber("MRawRunHeader") << " not found ... aborting" << endl;
-        return kFALSE;
-    }
-
-    fSrcPos = (MSrcPosCam*)pList->FindCreateObj("MSrcPosCam");
-    if (!fSrcPos)
-        return kFALSE;
-
-    return kTRUE;
+      return kFALSE;
+    }
+  
+  return kTRUE;
 }
 
 // --------------------------------------------------------------------------
 //
+// In case we enter a new run, look for the position in the read list
+// If there is no value for that run, take the previous one
 // 
-Bool_t MSrcPosFromFile::ReInit(MParList *pList)
-{
+Int_t MSrcPosFromFile::ComputeNewSrcPosition()
+{
+  const UInt_t run = fRawRunHeader->GetRunNumber();
+  if(run!=fLastRun)
+    {      
+      fLastRun=run;
+      MSrcPosCam* srcpos = (MSrcPosCam*)fRunMap->GetValue(run);
+      
+      Float_t x;
+      Float_t y;
+        
+      if (srcpos)
+        {
+          x = srcpos->GetX();
+          y = srcpos->GetY();
+          
+          GetSrcPosCam()->SetXY(x,y);
+                
+          *fLog << inf << "Source position for run " << run;
+          *fLog << inf << "\tX\t" << setprecision(3) << x;
+          *fLog << inf << "\tY\t" << setprecision(3) << y << endl;
+        }           
+      else
+        {
+          *fLog << inf << "Source position for run " << run << " not found in file. ";
+          *fLog << inf << "Taking previous position: ";
+          *fLog << inf << "\tX\t" << setprecision(3) << x;
+          *fLog << inf << "\tY\t" << setprecision(3) << y << endl;        
+        }
+    }
+  
   return kTRUE;
 }
@@ -122 +146 @@
 // --------------------------------------------------------------------------
 //
+// Read file with table of source positions as function of run
 // 
-Int_t MSrcPosFromFile::Process()
-{
-
-    switch(fMode)
+Int_t MSrcPosFromFile::ReadSourcePositionsFile(UShort_t readmode)
+{
+  
+  ifstream fin(fSourcePositionFilePath);
+  if(!fin)
+    {
+      *fLog << err << "MSrcPosFromFile::ReadSourcePositionsFile. Cannot open file " << fSourcePositionFilePath << endl;
+      return kFALSE;
+    }
+  
+  UInt_t run;
+  Float_t x,y;
+  
+  fNumRuns=0;
+  
+  *fLog << dbg << "MSrcPosFromFile::ReadSourcePositionsFile(" << readmode << ')' << endl;
+  while(1)
+    {
+      fin >> run >> x >> y;
+      if(fin.eof())
+        break;
+      
+      switch(readmode)
         {
-            case kOn:
-            {
-                const UInt_t run = fRawRunHeader->GetRunNumber();
-
-                MSrcPosCam* srcpos = (MSrcPosCam*)fRunMap->GetValue(run);
-                
-                Float_t x;
-                Float_t y;
-                
-                if (srcpos)
-                {
-                    x = srcpos->GetX();
-                    y = srcpos->GetY();
-                    
-                    
-                    if (srcpos && run != fLastRun)
-                    {
-                        fSrcPos->SetXY(x,y);
+        case kCount:
+          {
             
-                        *fLog << inf << "Source position for run " << run;
-                        *fLog << inf << "\tX\t" << setprecision(2) << x;
-                        *fLog << inf << "\tY\t" << setprecision(2) << y << endl;
-                    }
-                    
-                    fLastRun = run;
-                }
-                
-                x = fSrcPos->GetX();
-                y = fSrcPos->GetY();
-                
-                fHistSrcPos->Fill(x,y);
-                break;
-            }
-            case kOff:
-            {
-                Axis_t x;
-                Axis_t y;
-
-                fHistSrcPos->GetRandom2(x,y);
-                fSrcPos->SetXY(x,y);
-
-                break;
-            }
-            default:
-                *fLog << err << "Wrond mode " << fMode << endl;
-                return kFALSE;
+            *fLog << dbg << "Source position for run " << run;
+            *fLog << dbg << "\tX\t" << x << " mm";
+            *fLog << dbg << "\tY\t" << y << " mm" << endl;
+            
+            fNumRuns++;
+            break;
+          }
+        case kRead:
+          {
+            fRunList[fNumRuns] = run;
+            fRunSrcPos[fNumRuns].SetXY(x,y);
+            fRunMap->Add(fRunList[fNumRuns],(Long_t)&(fRunSrcPos[fNumRuns]));
+            fNumRuns++;
+            break;
+          }
+        default:
+          *fLog << err << "Read mode " << readmode << " node defined" << endl;
+          return kFALSE;
         }
-
-    return kTRUE;
-}
-
-
-Int_t MSrcPosFromFile::PostProcess()
-{
-
-    *fLog << dbg << endl;
-    *fLog << dbg << "fHistSrcPos->GetEntries() " << fHistSrcPos->GetEntries() << endl;
-    *fLog << dbg << "fHistSrcPos->ProjectionX()->GetMean() " << fHistSrcPos->ProjectionX()->GetMean() << endl;
-    *fLog << dbg << "fHistSrcPos->ProjectionX()->GetRMS() " << fHistSrcPos->ProjectionX()->GetRMS() << endl;
-    *fLog << dbg << "fHistSrcPos->ProjectionY()->GetMean() " << fHistSrcPos->ProjectionY()->GetMean() << endl;
-    *fLog << dbg << "fHistSrcPos->ProjectionY()->GetRMS() " << fHistSrcPos->ProjectionY()->GetRMS() << endl;
-
-    return kTRUE;
-}
-
-
-Int_t MSrcPosFromFile::ReadSourcePositionsFile(UShort_t readmode)
-{
-
-    ifstream fin(fSourcePositionFilePath);
-    if(!fin)
-    {
-        *fLog << err << "MSrcPosFromFile::ReadSourcePositionsFile. Cannot open file " << fSourcePositionFilePath << endl;
-        return kFALSE;
-    }
-
-    UInt_t run;
-    Float_t x,y;
-
-    fNumRuns=0;
-
-    *fLog << dbg << "MSrcPosFromFile::ReadSourcePositionsFile(" << readmode << ')' << endl;
-    while(1)
-    {
-        fin >> run >> x >> y;
-        if(fin.eof())
-            break;
-
-        switch(readmode)
-        {
-            case kCount:
-            {
-
-                *fLog << dbg << "Source position for run " << run;
-                *fLog << dbg << "\tX\t" << x << " mm";
-                *fLog << dbg << "\tY\t" << y << " mm" << endl;
-
-                fNumRuns++;
-                break;
-            }
-            case kRead:
-            {
-                fRunList[fNumRuns] = run;
-                fRunSrcPos[fNumRuns].SetXY(x,y);
-                fRunMap->Add(fRunList[fNumRuns],(Long_t)&(fRunSrcPos[fNumRuns]));
-                fNumRuns++;
-                break;
-            }
-            default:
-                *fLog << err << "Read mode " << readmode << " node defined" << endl;
-                return kFALSE;
-        }
-    }
-
-    fin.close();
-
-
-    return kTRUE;
-}
+    }
+  
+  fin.close();
+  
+  
+  return kTRUE;
+}

trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcPosFromFile.h

@@ -2 +2 @@
 #define MARS_MSrcPosFromFile
 
-#ifndef MARS_MTask
-#include "MTask.h"
+#ifndef MARS_MSrcPlace
+#include "MSrcPlace.h"
 #endif
 
@@ -10 +10 @@
 #endif
 
-class TH2F;
-
 class MRawRunHeader;
 class MSrcPosCam;
 
-class MSrcPosFromFile : public MTask
+class MSrcPosFromFile : public MSrcPlace
 {
  private:
 
     MRawRunHeader *fRawRunHeader;
-    MSrcPosCam    *fSrcPos;   //  Pointer to the source position
 
     Int_t  fNumRuns;
@@ -29 +26 @@
     TExMap     *fRunMap;   // list of run numbers positions
 
-    TH2F    *fHistSrcPos;
-    UShort_t fMode;
+    TString fSourcePositionFilePath;
 
-    TString fSourcePositionFilePath;
-    Int_t ReadSourcePositionsFile(UShort_t readmode);
-
-    Int_t PreProcess(MParList *plist);
-    Bool_t ReInit(MParList *plist);
-    Int_t Process();
-    Int_t PostProcess();
-
-public:
+    virtual Int_t ReadSourcePositionsFile(UShort_t readmode);
+    virtual Int_t ComputeNewSrcPosition();
     
+    virtual Int_t PreProcess(MParList *plist);
+    
+public:    
     enum ReadMode_t {kCount=0,kRead};
-    enum OnOffMode_t {kOn=0,kOff};
 
     MSrcPosFromFile(TString cardpath, OnOffMode_t mode=kOn, const char *name=NULL, const char *title=NULL);
     ~MSrcPosFromFile();
-    
-    void SetMode(OnOffMode_t mode) {fMode=mode;}
 
     ClassDef(MSrcPosFromFile, 0) // task to calculate the position of the source as a function of the run number 

trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcRotate.cc

@@ -26 +26 @@
 // MSrcRotate
 //
-// Task to rotate the position of the source as a function of Azimuth and 
-// Zenith angles
+// Task to rotate the position of the source as a function of time
 //
 //  Input Containers:
@@ -45 +44 @@
 #include "MRawEvtHeader.h"
 #include "MSrcRotate.h"
-#include "MSrcPosCam.h"
-#include "MDCA.h"
 
 #include "MAstroSky2Local.h"
 #include "MObservatory.h"
+#include "MSrcPosCam.h"
+#include "MDCA.h"
 
 #include "MLog.h"
@@ -59 +58 @@
 
 static const TString gsDefName  = "MSrcRotate";
-static const TString gsDefTitle = "De-rotate position of the source";
+static const TString gsDefTitle = "Rotate position of the source";
 
 // -------------------------------------------------------------------------
@@ -68 +67 @@
 //
 MSrcRotate::MSrcRotate(const char* srcPos, const char* dca, const char *name, const char *title)
-  : fSrcPos(NULL), fDCA(NULL), fRA(0), fDEC(0), fRunNumber(0)
+  : fRA(0), fDEC(0), fRefMJD(0), fRunNumber(0)
 {
     fName  = name  ? name  : gsDefName.Data();
     fTitle = title ? title : gsDefTitle.Data();
 
-    fSrcPosName  = srcPos;
-    fDCAName     = dca;
+    SetSrcPosName(srcPos);
+    SetDCAName(dca);
 }
 
@@ -84 +83 @@
 Int_t MSrcRotate::PreProcess(MParList *pList)
 {
-    // look for/create MSrcPosCam
-  fSrcPos = (MSrcPosCam*)pList->FindObject(AddSerialNumber(fSrcPosName), "MSrcPosCam");
-  if (!fSrcPos)
-    {
-      *fLog << warn << AddSerialNumber(fSrcPosName) << " [MSrcPosCam] not found... creating default container." << endl;
-      fSrcPos = (MSrcPosCam*)pList->FindCreateObj("MSrcPosCam", AddSerialNumber(fSrcPosName));
-      if(!fSrcPos)
-        return kFALSE;
-    }
-
-  // look for/create MDCA
-  fDCA = (MDCA*)pList->FindObject(AddSerialNumber(fDCAName), "MDCA");
-  if (!fDCA)
-    {
-      *fLog << warn << AddSerialNumber(fDCAName) << " [MDCA] not found... creating default container." << endl;
-      fDCA = (MDCA*)pList->FindCreateObj("MDCA", AddSerialNumber(fDCAName));
-      if(!fDCA)
-        return kFALSE;
-    }
-
+  if(!SearchForSrcPos(pList))
+    return kFALSE;
+  
   // look for MRawRunHeader
   fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
@@ -141 +123 @@
 // FIXME: for the time being, this is computed by assuming constant event rate
 //
-Int_t MSrcRotate::Process()
+Int_t MSrcRotate::ComputeNewSrcPosition()
 {  
   if(fRunHeader->GetRunNumber()!=fRunNumber)
     {
       fRunNumber=fRunHeader->GetRunNumber();
-
+      
       // save the number of events, initial and final times
       fNEvts    = fRunHeader->GetNumEvents();
@@ -155 +137 @@
       if((ULong_t)TMath::Nint(fIniTime.GetMjd())!=(ULong_t)TMath::Nint(fFinTime.GetMjd()))
         {
-          *fLog << err << "MSrcRotate::Process Error: Inial and final MJDs are different ("<<fIniTime.GetMjd()<<"!="<<fFinTime.GetMjd()<<")" << endl;
+          *fLog << err << "MSrcRotate::ComputeNewSrcPosition Error: Inial and final MJDs are different ("<<fIniTime.GetMjd()<<"!="<<fFinTime.GetMjd()<<")" << endl;
           return kFALSE;
         }
-
+      
 #ifdef DEBUG
       cout << endl << "********************"  << endl;
@@ -169 +151 @@
     }
   
-
+  
   // Compute the event time
   // FIXME: for the time being, this is computed by assuming constant event rate
@@ -177 +159 @@
   MTime refTime;
   refTime.SetMjd(fRefMJD);
-
+  
   // de-rotate the source position
   const MAstroSky2Local Observation(eventTime, *fObservatory);
   const MAstroSky2Local RefObservation(refTime, *fObservatory);
-
+  
 #ifdef DEBUG
   printf("Run:%d, Event:%d, iniMJD=%15.5f, finMJD=%15.5f, fDeltaT=%15.5f, newMJD=%15.5f, fRefMJD=%15.5f, rotation=%15.5f, ref=%15.5f\n",
@@ -190 +172 @@
   //  cout << "newMJD=" << newMJD << ", fRefMJD="<<fRefMJD<<", rotation="<<Observation.RotationAngle(fRA,fDEC)<<", ref="<<RefObservation.RotationAngle(fRA,fDEC)<< endl;
 #endif
-
-  Double_t rotationAngle = Observation.RotationAngle(fRA,fDEC)-RefObservation.RotationAngle(fRA,fDEC);
+  
+  Double_t rotationAngle = fRefMJD ? 
+    Observation.RotationAngle(fRA,fDEC)-RefObservation.RotationAngle(fRA,fDEC) :
+    Observation.RotationAngle(fRA,fDEC) ;
+  MSrcPosCam* srcpos = GetSrcPosCam();
+  MDCA* dca = GetDCA();
   
   Float_t c = TMath::Cos(rotationAngle);
   Float_t s = TMath::Sin(rotationAngle);
   // perform a rotation of -rotationAngle to move the source back to the "initial" position
-  Float_t newX = c*fSrcPos->GetX()-s*fSrcPos->GetY();
-  Float_t newY = s*fSrcPos->GetX()+c*fSrcPos->GetY();
-
+  Float_t newX = c*srcpos->GetX()-s*srcpos->GetY();
+  Float_t newY = s*srcpos->GetX()+c*srcpos->GetY();
+  
 #ifdef DEBUG
   Double_t rotationAngleComp = fObservatory->RotationAngle(0.1256,2.63);
   cout << "Event " << fEvtHeader->GetDAQEvtNumber() << endl;
-  cout << "newMJD=" << newMJD <<", rotationAngle=" << rotationAngle <<", rotationAngleComp=" << rotationAngleComp << ", oldX="<<fIniSrcPos.GetX()<< ", oldY="<<fIniSrcPos.GetY()<< ", newX="<<newX<< ", newY="<<newY << endl,
+  cout << "newMJD=" << newMJD <<", rotationAngle=" << rotationAngle <<", rotationAngleComp=" << rotationAngleComp << ", oldX="<<fIniSrcPos.GetX()<< ", oldY="<<fIniSrcPos.GetY()<< ", newX="<<newX<< ", newY="<<newY << endl;
 #endif
-
-  fSrcPos->SetX(newX);
-  fSrcPos->SetY(newY);
-  fDCA->SetRefPoint(newX,newY);
+  
+  srcpos->SetX(newX);
+  srcpos->SetY(newY);
+  dca->SetRefPoint(newX,newY);
 
   return kTRUE;

trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcRotate.h

@@ -2 +2 @@
 #define MARS_MSrcRotate
 
-#ifndef MARS_MTask
-#include "MTask.h"
+#ifndef MARS_MSrcPlace
+#include "MSrcPlace.h"
 #endif
 
-#include <TArrayF.h>
 #include "MTime.h"
-#include "MSrcPosCam.h"
 
-class MDCA;
 class MObservatory;
 class MRawEvtHeader;
 class MRawRunHeader;
 
-class MSrcRotate : public MTask
+class MSrcRotate : public MSrcPlace
 {
-private:
-    MSrcPosCam*    fSrcPos;      //  Pointer to the source position
-    MDCA*          fDCA;         //  Pointer to the MDCA object
-    MRawEvtHeader* fEvtHeader;   //  Pointer to the event header
-    MRawRunHeader* fRunHeader;   //  Pointer to the run header
-    MObservatory*  fObservatory; //  Pointer to the MObservatory    
-
-    TString     fSrcPosName;
-    TString     fDCAName;
-
-    Double_t   fRA;          //  [rad] Right ascention 
-    Double_t   fDEC;         //  [rad] Declination
-    Double_t   fRefMJD;      //  [MJ date] reference time for rotation
-    UInt_t     fNEvts;       //  Number of events in file
-    MTime      fIniTime;     //  Run initial time 
-    MTime      fFinTime;     //  Run final time 
-    Double_t   fDeltaT;      //  DeltaT between two events
-    UInt_t     fRunNumber;   //  Current run number
-
-    Int_t PreProcess(MParList *plist);
-    Int_t Process();
-
-public:
-    MSrcRotate(const char* src="MSrcPosCam", const char* dca="MDCA",
-               const char* name=NULL, const char* title=NULL);
-
-    void SetRAandDECandRefMJD(Double_t ra, Double_t dec, Double_t ref) {fRA=ra;fDEC=dec;fRefMJD=ref;}
-
-    ClassDef(MSrcRotate, 0) // task to rotate the position of the source as a function of Azimuth and Zenith angles
+ private:
+  MRawEvtHeader* fEvtHeader;   //  Pointer to the event header
+  MRawRunHeader* fRunHeader;   //  Pointer to the run header
+  MObservatory*  fObservatory; //  Pointer to the MObservatory    
+  
+  Double_t   fRA;          //  [rad] Right ascenssion 
+  Double_t   fDEC;         //  [rad] Declination
+  Double_t   fRefMJD;      //  [MJ date] reference time for rotation
+  UInt_t     fNEvts;       //  Number of events in file
+  MTime      fIniTime;     //  Run initial time 
+  MTime      fFinTime;     //  Run final time 
+  Double_t   fDeltaT;      //  DeltaT between two events
+  UInt_t     fRunNumber;   //  Current run number
+  
+  virtual Int_t PreProcess(MParList *plist);
+  
+ public:
+  MSrcRotate(const char* src="MSrcPosCam", const char* dca="MDCA",
+             const char* name=NULL, const char* title=NULL);
+  
+  void SetRAandDECandRefMJD(Double_t ra, Double_t dec, Double_t ref=0) {fRA=ra;fDEC=dec;fRefMJD=ref;}
+  void SetRAandDEC(Double_t ra, Double_t dec){SetRAandDECandRefMJD(ra,dec);}
+  virtual Int_t ComputeNewSrcPosition();
+  
+  ClassDef(MSrcRotate, 0) // task to rotate the position of the source as a function of Azimuth and Zenith angles
 };
 

trunk/MagicSoft/Mars/mtemp/mifae/library/Makefile

@@ -51 +51 @@
         MPSFFit.cc \
         MPSFFitCalc.cc \
+        MSrcPlace.cc \
         MSrcPosFromFile.cc \
         MSrcRotate.cc \