StRoot: StSvtSeqAdjMaker/StSvtSeqAdjMaker.h Source File

00001 // $Id: StSvtSeqAdjMaker.h,v 1.26 2007/07/12 20:10:35 fisyak Exp $
00002 // $Log: StSvtSeqAdjMaker.h,v $
00003 // Revision 1.26  2007/07/12 20:10:35  fisyak
00004 // Add forward declaration for ROOT 5.16
00005 //
00006 // Revision 1.25  2005/07/23 03:37:34  perev
00007 // IdTruth + Cleanup
00008 //
00009 // Revision 1.24  2004/03/18 04:05:02  caines
00010 // Remove from global scope variables used in debug mode as they shouldnt be there and caused erratic behaviour, also initialise some variables that valgrind was complaining about - didnt really need to as they are sent back from function which initialises them properly always but doesnt hurt
00011 //
00012 // Revision 1.23  2003/09/10 19:47:36  perev
00013 // ansi corrs
00014 //
00015 // Revision 1.22  2003/09/02 17:59:08  perev
00016 // gcc 3.2 updates + WarnOff
00017 //
00018 // Revision 1.21  2003/07/18 17:15:40  caines
00019 // Fix Pedoffset to be 20 not 10 change variables to int from floats to avoid casting problems, fix that when pedestal goes negative we dont
00020 //
00021 // Revision 1.20  2003/01/21 01:27:44  jeromel
00022 // hfile->write(0 while NULL caused spurious crash.
00023 // Took the oportunity to add GetCVS()
00024 // Some maniaco-compulsive //! alignement fixes ...
00025 //
00026 // Revision 1.19  2002/09/20 19:35:25  caines
00027 // Change building of file name
00028 //
00029 // Revision 1.18  2002/09/19 16:17:49  caines
00030 // Add code to do Juns gain calibration
00031 //
00032 // Revision 1.17  2002/05/09 16:55:08  munhoz
00033 // add reading bad anodes from DB
00034 //
00035 // Revision 1.16  2002/04/22 14:52:54  caines
00036 // Add the .h too
00037 //
00038 // Revision 1.15  2002/01/11 22:49:22  caines
00039 // Fix sequence merging bugs-hopefully
00040 //
00041 // Revision 1.14  2001/12/13 03:08:17  caines
00042 // Can now subtract common mode noise via black anodes 239 and 2
00043 //
00044 // Revision 1.13  2001/10/19 23:31:34  caines
00045 // Correct problem that if anodes were missing didnt do average common mode noise calc
00046 //
00047 // Revision 1.12  2001/09/26 18:42:48  caines
00048 // Fix 2 anode subtraction routines
00049 //
00050 // Revision 1.11  2001/09/16 22:24:26  caines
00051 // Fix for when SVT isnt in every event
00052 //
00053 // Revision 1.10  2001/08/24 20:57:46  caines
00054 // Do common mode noise suppression from first two anodes
00055 //
00056 // Revision 1.9  2001/07/25 14:47:47  caines
00057 // Fix filling histogram only when debug is on
00058 //
00059 // Revision 1.8  2001/07/22 20:31:29  caines
00060 // Better tuning for real data. Common mode noise calc and sub. Avoid overlapping seq. Fill histograms only in debug
00061 //
00062 // Revision 1.7  2001/05/02 02:07:55  caines
00063 // Fix function declaration for Solaris
00064 //
00065 // Revision 1.6  2001/05/01 00:23:58  caines
00066 //  Update h files for use with zsp data
00067 //
00068 // Revision 1.5  2000/11/30 20:45:56  caines
00069 // Dynamically calc prob values, use database
00070 //
00071 // Revision 1.4  2000/08/24 04:23:50  caines
00072 // Improved histograms
00073 //
00074 // Revision 1.3  2000/08/21 12:57:31  caines
00075 // Now opens and reads in ped using CalibMaker
00076 //
00077 // Revision 1.2  2000/07/16 22:32:23  caines
00078 // Now also saves RAW data
00079 //
00080 // Revision 1.1  2000/06/15 20:04:54  caines
00081 // Initial versions of sequence adjusting codes
00082 //
00083 //
00084 #ifndef STAR_StSvtSeqAdj
00085 #define STAR_StSvtSeqAdj
00087 //                                                                      //
00088 // StSvtSeqAdj base class                                               //
00089 //                                                                      //
00091 #ifndef StMaker_H
00092 #include "StMaker.h"
00093 #endif
00094 
00095 #include <Stiostream.h>
00096 #include "Stiostream.h"
00097 #include <string>
00098 
00099 class TH1D;
00100 class TH2F;
00101 class TFile; 
00102 class StSequence;
00103 class StSvtData;
00104 class StSvtHybridData;
00105 class StSvtInverseProducts;
00106 class StSvtPedSub;
00107 class TObjectSet;
00108 class StSvtHybridCollection;
00109 class StSvtBadAnode;
00110 class StSvtProbValues;
00111 class StMCTruth;
00112 
00113 class StSvtSeqAdjMaker : public StMaker
00114 {
00115  public: 
00116   StSvtSeqAdjMaker(const char *name = "SvtSeqAdj");
00117   StSvtSeqAdjMaker(StSvtSeqAdjMaker& analmaker);
00118   ~StSvtSeqAdjMaker();
00119 
00120   virtual Int_t Init();
00121   virtual Int_t InitRun( int runnumber);
00122   virtual Int_t Make();
00123   virtual Int_t Finish();
00124 
00125   Int_t SetSvtData();
00126   Int_t GetSvtRawData();
00127   Int_t GetSvtPedestals();
00128   Int_t GetBadAnodes();
00129   Int_t GetPedOffset(){return mPedOffSet;};
00130   void CommonModeNoiseCalc(int iAnode);
00131   void CommonModeNoiseSub(int iAnode);
00132   void SubtractFirstAnode(int iAnode);
00133   Int_t FindBlackAnodes(); // Find Black anodes on each hybrid
00134   Int_t AdjustSequences1( int iAnode, int Anode); // Find sequences  based on ASICS
00135   Int_t AdjustSequences2(int iAnode, int Anode); //adjust sequences base on LowInvProd
00136   Int_t MergeSequences(StSequence* Seq, int nSeq,StMCTruth *Tru=0); // Merge overlapping sequences
00137   Int_t CreateHist(Int_t tNuOfHyb);
00138   void  MakeHistogramsProb(int index,int Anode);
00139   void  MakeHistogramsAdc(StSvtHybridData* hybridData, int index,int Anode, int Count);
00140   Int_t Reset();
00141 
00142   Int_t SetMinAdcLevels( int MinAdc1,  int MinAbove1, int MinAdc2, int MinAbove2, int PedOffset); // Set the 2 thresholds for a sequence
00143   Int_t SetPedestalFile(const char* PedFile);
00144   Int_t SetLowInvProd(int LowInvProd);// Set the low threshold based on the frequency distribution
00145 
00146   virtual const char *GetCVS() const {
00147     static const char cvs[]="Tag $Name:  $ $Id: StSvtSeqAdjMaker.h,v 1.26 2007/07/12 20:10:35 fisyak Exp $ built "__DATE__" "__TIME__ ; 
00148     return cvs;
00149   }
00150     
00151  protected:
00152   St_ObjectSet* mSvtDataSet;              
00153   StSvtData* mSvtRawData;                 
00154   StSvtData* mSvtAdjData;                 
00155   StSvtHybridData* mHybridRawData ;       
00156   StSvtHybridData* mHybridAdjData ;       
00157 
00158   StSvtHybridCollection* mSvtBadAnodes;   
00159   StSvtHybridCollection* mSvtPedColl;     
00160 
00161   StSvtInverseProducts* mInvProd;         
00162   StSvtPedSub* mSvtPedSub;                
00163   StSvtProbValues* mProbValue;            
00164 
00165   const char* mPedFile;                   
00166 
00167   int* anolist;   
00168   TFile *hfile;   
00169   unsigned long Evt_counts;  
00170 
00171   TH1F* mOcupancyHisto;                   
00172   TH1F* EventOccupancy;                   
00173   TH1D** mInvProdSeqAdj;                  
00174   TH1F** mRawAdc;                         
00175   TH1F** mAdcAfter;                       
00176   TH1F* mCommonModePitch;                 
00177   TH1F* mCommonModeCount;                 
00178   TH2F** mTimeAn;                         
00179   TH1F** mRawPixel;                       
00180   int mNumOfSeq;
00181   int m_n_seq_lo;
00182   int m_n_seq_hi;
00183   int m_thresh_lo;
00184   int m_inv_prod_lo;
00185   int m_thresh_hi;
00186   
00187   int mTotalNumberOfHybrids;
00188   int mPedOffSet;
00189   int mCommonModeNoise[128];
00190   int mCommonModeNoiseAn[128];
00191   int mNAnodes;
00192   int doCommon;
00193   int adcCommon[128];
00194 
00195  private:
00196   string buildFileName(string dir, string fileName, string extention);
00197   string baseName(string s);
00198   ClassDef(StSvtSeqAdjMaker,0)   //virtual base class for Makers
00199 
00200 };
00201 
00202 
00203 #endif