00001
00002 #include "TH1.h"
00003 #include "TH2.h"
00004 #include "TProfile.h"
00005 #include "TMath.h"
00006 #include "TRandom.h"
00007
00008 static const char* trackFileName="/afs/fnal.gov/files/data/minos/d03/tobien/track_files/tracks_far_mc_comb_new.out";
00009 static const int nTracks = 20000;
00010 static const int nPlanes = 244;
00011 static const int datPointsMin = 5;
00012 static const int datPointsMax = 200;
00013 static const float residMax = 3.5;
00014 static const float maxTDiff = 15.;
00015 static const int maxIter = 20;
00016
00017
00018 static const double offsetRMS = 0.0;
00019 static const double detectorRMS = 1.;
00020
00021
00022
00023
00024
00025
00026 class track_class {
00027
00028 public:
00029
00030
00031 double dataU[nPlanes][nTracks];
00032 double dataV[nPlanes][nTracks];
00033 double dataUZ[nPlanes];
00034 double dataVZ[nPlanes];
00035
00036 double offsetU[nPlanes];
00037 double offsetV[nPlanes];
00038 double sig[nPlanes];
00039 double avgresU[nPlanes];
00040 double avgresV[nPlanes];
00041 double resU[nPlanes];
00042 double resV[nPlanes];
00043
00044 double stdresU[nPlanes];
00045 double stdresV[nPlanes];
00046 double mavg,bavg,mstd,bstd;
00047
00048 void fittrkU(int n_iter);
00049 void fittrkV(int n_iter);
00050 void datread();
00051 void moddat(int planeOr);
00052 void print_off(int n_iter);
00053 int getZI(int iplane, int tOr);
00054 void offgen();
00055 void offsetread();
00056
00057 track_class();
00058
00059 TH1F *h1, *h2,*h3,*h4,*h5,*h6,*h7,*h20,*h50,*h60;
00060 TH1F *h11, *h12,*h13,*h14,*h15,*h16,*h17,*h21,*h51,*h61;
00061
00062 TH1F *h101, *h102, *h103, *h104, *h105;
00063 TH1F *h111, *h112, *h113, *h114, *h115;
00064 TH1F *h121, *h122, *h123, *h124, *h125;
00065 TH1F *h201, *h202, *h203, *h204, *h205;
00066 TH1F *h211, *h212, *h213, *h214, *h215;
00067 TH1F *h221, *h222, *h223, *h224, *h225;
00068
00069 TH2F *h30, *h31, *h40, *h41, *h70, *h71;;
00070
00071 TProfile *h90, *h80;
00072 TProfile *h91, *h81;
00073 };
00074
