Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 11:58:25

0001 #ifndef vDriftHistos_H
0002 #define vDriftHistos_H
0003 
0004 #include "TH1.h"
0005 #include "TString.h"
0006 #include "TFile.h"
0007 #include "DTTMax.h"
0008 #include <string>
0009 
0010 // A set of histograms on chamber angle and position
0011 class h4DSegm {
0012 public:
0013   h4DSegm(std::string name_) {
0014     TString N = name_.c_str();
0015     name = name_.c_str();
0016     h4DSegmXPosInCham = new TH1F(N + "_h4DSegmXPosInCham", "4D Segment x position (cm) in Chamber RF", 200, -200, 200);
0017     h4DSegmYPosInCham = new TH1F(N + "_h4DSegmYPosInCham", "4D Segment y position (cm) in Chamber RF", 200, -200, 200);
0018     h4DSegmPhiAngleInCham =
0019         new TH1F(N + "_h4DSegmPhiAngleInCham", "4D Segment phi angle (rad) in Chamber RF", 180, -180, 180);
0020     h4DSegmThetaAngleInCham =
0021         new TH1F(N + "_h4DSegmThetaAngleInCham", "4D Segment theta angle (rad) in Chamber RF", 180, -180, 180);
0022     h4DSegmImpactAngleInCham =
0023         new TH1F(N + "_h4DSegmImpactAngleInCham", "4D Segment impact angle (rad) in Chamber RF", 180, -180, 180);
0024   }
0025   h4DSegm(const TString &name_, TFile *file) {
0026     name = name_;
0027     h4DSegmXPosInCham = (TH1F *)file->Get(name + "_h4DSegmXPosInCham");
0028     h4DSegmYPosInCham = (TH1F *)file->Get(name + "_h4DSegmYPosInCham");
0029     h4DSegmPhiAngleInCham = (TH1F *)file->Get(name + "_h4DSegmPhiAngleInCham");
0030     h4DSegmThetaAngleInCham = (TH1F *)file->Get(name + "_h4DSegmThetaAngleInCham");
0031     h4DSegmImpactAngleInCham = (TH1F *)file->Get(name + "_h4DSegmImpactAngleInCham");
0032   }
0033   ~h4DSegm() {
0034     delete h4DSegmXPosInCham;
0035     delete h4DSegmYPosInCham;
0036     delete h4DSegmPhiAngleInCham;
0037     delete h4DSegmThetaAngleInCham;
0038     delete h4DSegmImpactAngleInCham;
0039   }
0040   void Fill(float x, float y, float phi, float theta, float impact) {
0041     h4DSegmXPosInCham->Fill(x);
0042     h4DSegmYPosInCham->Fill(y);
0043     h4DSegmPhiAngleInCham->Fill(phi);
0044     h4DSegmThetaAngleInCham->Fill(theta);
0045     h4DSegmImpactAngleInCham->Fill(impact);
0046   }
0047   void Fill(float x, float phi) {
0048     h4DSegmXPosInCham->Fill(x);
0049     h4DSegmPhiAngleInCham->Fill(phi);
0050   }
0051   void Write() {
0052     h4DSegmXPosInCham->Write();
0053     h4DSegmYPosInCham->Write();
0054     h4DSegmPhiAngleInCham->Write();
0055     h4DSegmThetaAngleInCham->Write();
0056     h4DSegmImpactAngleInCham->Write();
0057   }
0058 
0059 public:
0060   TH1F *h4DSegmXPosInCham;
0061   TH1F *h4DSegmYPosInCham;
0062   TH1F *h4DSegmPhiAngleInCham;
0063   TH1F *h4DSegmThetaAngleInCham;
0064   TH1F *h4DSegmImpactAngleInCham;
0065 
0066   TString name;
0067 };
0068 
0069 // A set of histograms on SL angle and position
0070 class h2DSegm {
0071 public:
0072   h2DSegm(std::string name_) {
0073     TString N = name_.c_str();
0074     name = name_.c_str();
0075     h2DSegmPosInCham = new TH1F(N + "_h2DSegmPosInCham", "2D Segment position (cm) in Chamber RF", 200, -200, 200);
0076     h2DSegmAngleInCham = new TH1F(N + "_h2DSegmAngleInCham", "2D Segment angle (rad) in Chamber RF", 200, -2, 2);
0077     h2DSegmCosAngleInCham = new TH1F(N + "_h2DSegmCosAngleInCham", "2D Segment cos(angle) in Chamber RF", 200, -2, 2);
0078   }
0079   h2DSegm(const TString &name_, TFile *file) {
0080     name = name_;
0081 
0082     h2DSegmPosInCham = (TH1F *)file->Get(name + "_h2DSegmPosInCham");
0083     h2DSegmAngleInCham = (TH1F *)file->Get(name + "_h2DSegmAngleInCham");
0084     h2DSegmCosAngleInCham = (TH1F *)file->Get(name + "_h2DSegmCosAngleInCham");
0085   }
0086   ~h2DSegm() {
0087     delete h2DSegmPosInCham;
0088     delete h2DSegmAngleInCham;
0089     delete h2DSegmCosAngleInCham;
0090   }
0091   void Fill(float pos, float localAngle) {
0092     h2DSegmPosInCham->Fill(pos);
0093     h2DSegmAngleInCham->Fill(atan(localAngle));
0094     h2DSegmCosAngleInCham->Fill(cos(atan(localAngle)));
0095   }
0096   void Write() {
0097     h2DSegmPosInCham->Write();
0098     h2DSegmAngleInCham->Write();
0099     h2DSegmCosAngleInCham->Write();
0100   }
0101 
0102 public:
0103   TH1F *h2DSegmPosInCham;
0104   TH1F *h2DSegmAngleInCham;
0105   TH1F *h2DSegmCosAngleInCham;
0106 
0107   TString name;
0108 };
0109 
0110 // A set of histograms on SL Tmax
0111 class hTMaxCell {
0112 public:
0113   hTMaxCell(const TString &name_) {
0114     name = name_;
0115 
0116     // book TMax histograms
0117     hTmax123 = new TH1F(name + "_Tmax123", "Tmax123 value", 2000, -1000., 1000.);
0118     hTmax124s72 = new TH1F(name + "_Tmax124_s72", "Tmax124 sigma=sqrt(7/2) value", 2000, -1000., 1000.);
0119     hTmax124s78 = new TH1F(name + "_Tmax124_s78", "Tmax124 sigma=sqrt(7/8) value", 2000, -1000., 1000.);
0120     hTmax134s72 = new TH1F(name + "_Tmax134_s72", "Tmax134 sigma=sqrt(7/2) value", 2000, -1000., 1000.);
0121     hTmax134s78 = new TH1F(name + "_Tmax134_s78", "Tmax134 sigma=sqrt(7/8) value", 2000, -1000., 1000.);
0122     hTmax234 = new TH1F(name + "_Tmax234", "Tmax234 value", 2000, -1000., 1000.);
0123     hTmax_3t0 = new TH1F(name + "_3t0", "Tmax+3*Delta(t0)", 2000, -1000., 1000.);
0124     hTmax_3t0_0 = new TH1F(name + "_3t0_0", "Tmax+3*Delta(t0); 3 hits", 2000, -1000., 1000.);
0125     hTmax_3t0_1 = new TH1F(name + "_3t0_1", "Tmax+3*Delta(t0); one t<5ns", 2000, -1000., 1000.);
0126     hTmax_3t0_2 = new TH1F(name + "_3t0_2", "Tmax+3*Delta(t0); one t<10ns", 2000, -1000., 1000.);
0127     hTmax_3t0_3 = new TH1F(name + "_3t0_3", "Tmax+3*Delta(t0); one t<20ns", 2000, -1000., 1000.);
0128     hTmax_3t0_4 = new TH1F(name + "_3t0_4", "Tmax+3*Delta(t0); one t<50ns", 2000, -1000., 1000.);
0129     hTmax_3t0_5 = new TH1F(name + "_3t0_5", "Tmax+3*Delta(t0); all t>50ns", 2000, -1000., 1000.);
0130     hTmax_2t0 = new TH1F(name + "_2t0", "Tmax+2*Delta(t0)", 2000, -1000., 1000.);
0131     hTmax_2t0_0 = new TH1F(name + "_2t0_0", "Tmax+2*Delta(t0); 3 hits", 2000, -1000., 1000.);
0132     hTmax_2t0_1 = new TH1F(name + "_2t0_1", "Tmax+2*Delta(t0); one t<5ns", 2000, -1000., 1000.);
0133     hTmax_2t0_2 = new TH1F(name + "_2t0_2", "Tmax+2*Delta(t0); one t<10ns", 2000, -1000., 1000.);
0134     hTmax_2t0_3 = new TH1F(name + "_2t0_3", "Tmax+2*Delta(t0); one t<20ns", 2000, -1000., 1000.);
0135     hTmax_2t0_4 = new TH1F(name + "_2t0_4", "Tmax+2*Delta(t0); one t<50ns", 2000, -1000., 1000.);
0136     hTmax_2t0_5 = new TH1F(name + "_2t0_5", "Tmax+2*Delta(t0); all t>50ns", 2000, -1000., 1000.);
0137     hTmax_t0 = new TH1F(name + "_t0", "Tmax+Delta(t0)", 2000, -1000., 1000.);
0138     hTmax_t0_0 = new TH1F(name + "_t0_0", "Tmax+Delta(t0); 3 hits", 2000, -1000., 1000.);
0139     hTmax_t0_1 = new TH1F(name + "_t0_1", "Tmax+Delta(t0); one t<5ns", 2000, -1000., 1000.);
0140     hTmax_t0_2 = new TH1F(name + "_t0_2", "Tmax+Delta(t0); one t<10ns", 2000, -1000., 1000.);
0141     hTmax_t0_3 = new TH1F(name + "_t0_3", "Tmax+Delta(t0); one t<20ns", 2000, -1000., 1000.);
0142     hTmax_t0_4 = new TH1F(name + "_t0_4", "Tmax+Delta(t0); one t<50ns", 2000, -1000., 1000.);
0143     hTmax_t0_5 = new TH1F(name + "_t0_5", "Tmax+Delta(t0); all t>50ns", 2000, -1000., 1000.);
0144     hTmax_0 = new TH1F(name + "_0", "Tmax", 2000, -1000., 1000.);
0145   }
0146 
0147   hTMaxCell(const TString &name_, TFile *file) {
0148     name = name_;
0149     hTmax123 = (TH1F *)file->Get(name + "_Tmax123");
0150     hTmax124s72 = (TH1F *)file->Get(name + "_Tmax124_s72");
0151     hTmax124s78 = (TH1F *)file->Get(name + "_Tmax124_s78");
0152     hTmax134s72 = (TH1F *)file->Get(name + "_Tmax134_s72");
0153     hTmax134s78 = (TH1F *)file->Get(name + "_Tmax134_s78");
0154     hTmax234 = (TH1F *)file->Get(name + "_Tmax234");
0155     hTmax_3t0 = (TH1F *)file->Get(name + "_3t0");
0156     hTmax_3t0_0 = (TH1F *)file->Get(name + "_3t0_0");
0157     hTmax_3t0_1 = (TH1F *)file->Get(name + "_3t0_1");
0158     hTmax_3t0_2 = (TH1F *)file->Get(name + "_3t0_2");
0159     hTmax_3t0_3 = (TH1F *)file->Get(name + "_3t0_3");
0160     hTmax_3t0_4 = (TH1F *)file->Get(name + "_3t0_4");
0161     hTmax_3t0_5 = (TH1F *)file->Get(name + "_3t0_5");
0162     hTmax_2t0 = (TH1F *)file->Get(name + "_2t0");
0163     hTmax_2t0_0 = (TH1F *)file->Get(name + "_2t0_0");
0164     hTmax_2t0_1 = (TH1F *)file->Get(name + "_2t0_1");
0165     hTmax_2t0_2 = (TH1F *)file->Get(name + "_2t0_2");
0166     hTmax_2t0_3 = (TH1F *)file->Get(name + "_2t0_3");
0167     hTmax_2t0_4 = (TH1F *)file->Get(name + "_2t0_4");
0168     hTmax_2t0_5 = (TH1F *)file->Get(name + "_2t0_5");
0169     hTmax_t0 = (TH1F *)file->Get(name + "_t0");
0170     hTmax_t0_1 = (TH1F *)file->Get(name + "_t0_1");
0171     hTmax_t0_2 = (TH1F *)file->Get(name + "_t0_2");
0172     hTmax_t0_3 = (TH1F *)file->Get(name + "_t0_3");
0173     hTmax_t0_4 = (TH1F *)file->Get(name + "_t0_4");
0174     hTmax_t0_5 = (TH1F *)file->Get(name + "_t0_5");
0175     hTmax_0 = (TH1F *)file->Get(name + "_0");
0176   }
0177 
0178   ~hTMaxCell() {
0179     delete hTmax123;
0180     delete hTmax124s72;
0181     delete hTmax124s78;
0182     delete hTmax134s72;
0183     delete hTmax134s78;
0184     delete hTmax234;
0185     delete hTmax_3t0;
0186     delete hTmax_3t0_0;
0187     delete hTmax_3t0_1;
0188     delete hTmax_3t0_2;
0189     delete hTmax_3t0_3;
0190     delete hTmax_3t0_4;
0191     delete hTmax_3t0_5;
0192     delete hTmax_2t0;
0193     delete hTmax_2t0_0;
0194     delete hTmax_2t0_1;
0195     delete hTmax_2t0_2;
0196     delete hTmax_2t0_3;
0197     delete hTmax_2t0_4;
0198     delete hTmax_2t0_5;
0199     delete hTmax_t0;
0200     delete hTmax_t0_0;
0201     delete hTmax_t0_1;
0202     delete hTmax_t0_2;
0203     delete hTmax_t0_3;
0204     delete hTmax_t0_4;
0205     delete hTmax_t0_5;
0206     delete hTmax_0;
0207   }
0208 
0209   void Fill(float tmax123,
0210             float tmax124,
0211             float tmax134,
0212             float tmax234,
0213             dttmaxenums::SigmaFactor s124,  // Give the factor relating the width of the TMax distribution
0214             dttmaxenums::SigmaFactor s134,  // and the cell resolution (for tmax123 and tmax234 is always= sqrt(3/2)).
0215             unsigned t0_123,                // Give the "quantity" of Delta(t0) included in the tmax
0216             unsigned t0_124,                // formula used for Tmax123 or Tma124 or Tma134 or Tma234
0217             unsigned t0_134,
0218             unsigned t0_234,
0219             unsigned hSubGroup  //different t0 hists(at least one hit within a given distance from the wire)
0220   ) {
0221     if (tmax123 > 0.) {
0222       hTmax123->Fill(tmax123);
0223       if (t0_123 == 1) {
0224         hTmax_t0->Fill(tmax123);
0225         switch (hSubGroup) {
0226           case 0:
0227             hTmax_t0_0->Fill(tmax123);
0228             break;
0229           case 1:
0230             hTmax_t0_1->Fill(tmax123);
0231             break;
0232           case 2:
0233             hTmax_t0_2->Fill(tmax123);
0234             break;
0235           case 3:
0236             hTmax_t0_3->Fill(tmax123);
0237             break;
0238           case 4:
0239             hTmax_t0_4->Fill(tmax123);
0240             break;
0241           case 99:
0242             hTmax_t0_5->Fill(tmax123);
0243             break;
0244         }
0245       } else {
0246         hTmax_2t0->Fill(tmax123);
0247         switch (hSubGroup) {
0248           case 0:
0249             hTmax_2t0_0->Fill(tmax123);
0250             break;
0251           case 1:
0252             hTmax_2t0_1->Fill(tmax123);
0253             break;
0254           case 2:
0255             hTmax_2t0_2->Fill(tmax123);
0256             break;
0257           case 3:
0258             hTmax_2t0_3->Fill(tmax123);
0259             break;
0260           case 4:
0261             hTmax_2t0_4->Fill(tmax123);
0262             break;
0263           case 99:
0264             hTmax_2t0_5->Fill(tmax123);
0265             break;
0266         }
0267       }
0268     }
0269     if (tmax124 > 0.) {
0270       (s124 == dttmaxenums::r72) ? hTmax124s72->Fill(tmax124) : hTmax124s78->Fill(tmax124);
0271       if (t0_124 == 0)
0272         hTmax_0->Fill(tmax124);
0273       else if (t0_124 == 1) {
0274         hTmax_t0->Fill(tmax124);
0275         switch (hSubGroup) {
0276           case 0:
0277             hTmax_t0_0->Fill(tmax124);
0278             break;
0279           case 1:
0280             hTmax_t0_1->Fill(tmax124);
0281             break;
0282           case 2:
0283             hTmax_t0_2->Fill(tmax124);
0284             break;
0285           case 3:
0286             hTmax_t0_3->Fill(tmax124);
0287             break;
0288           case 4:
0289             hTmax_t0_4->Fill(tmax124);
0290             break;
0291           case 99:
0292             hTmax_t0_5->Fill(tmax124);
0293             break;
0294         }
0295       } else if (t0_124 == 2) {
0296         hTmax_2t0->Fill(tmax124);
0297         switch (hSubGroup) {
0298           case 0:
0299             hTmax_2t0_0->Fill(tmax124);
0300             break;
0301           case 1:
0302             hTmax_2t0_1->Fill(tmax124);
0303             break;
0304           case 2:
0305             hTmax_2t0_2->Fill(tmax124);
0306             break;
0307           case 3:
0308             hTmax_2t0_3->Fill(tmax124);
0309             break;
0310           case 4:
0311             hTmax_2t0_4->Fill(tmax124);
0312             break;
0313           case 99:
0314             hTmax_2t0_5->Fill(tmax124);
0315             break;
0316         }
0317       } else if (t0_124 == 3) {
0318         hTmax_3t0->Fill(tmax124);
0319         switch (hSubGroup) {
0320           case 0:
0321             hTmax_3t0_0->Fill(tmax124);
0322             break;
0323           case 1:
0324             hTmax_3t0_1->Fill(tmax124);
0325             break;
0326           case 2:
0327             hTmax_3t0_2->Fill(tmax124);
0328             break;
0329           case 3:
0330             hTmax_3t0_3->Fill(tmax124);
0331             break;
0332           case 4:
0333             hTmax_3t0_4->Fill(tmax124);
0334             break;
0335           case 99:
0336             hTmax_3t0_5->Fill(tmax124);
0337             break;
0338         }
0339       }
0340     }
0341     if (tmax134 > 0.) {
0342       (s134 == dttmaxenums::r72) ? hTmax134s72->Fill(tmax134) : hTmax134s78->Fill(tmax134);
0343       if (t0_134 == 0)
0344         hTmax_0->Fill(tmax134);
0345       else if (t0_134 == 1) {
0346         hTmax_t0->Fill(tmax134);
0347         switch (hSubGroup) {
0348           case 0:
0349             hTmax_t0_0->Fill(tmax134);
0350             break;
0351           case 1:
0352             hTmax_t0_1->Fill(tmax134);
0353             break;
0354           case 2:
0355             hTmax_t0_2->Fill(tmax134);
0356             break;
0357           case 3:
0358             hTmax_t0_3->Fill(tmax134);
0359             break;
0360           case 4:
0361             hTmax_t0_4->Fill(tmax134);
0362             break;
0363           case 99:
0364             hTmax_t0_5->Fill(tmax134);
0365             break;
0366         }
0367       } else if (t0_134 == 2) {
0368         hTmax_2t0->Fill(tmax134);
0369         switch (hSubGroup) {
0370           case 0:
0371             hTmax_2t0_0->Fill(tmax134);
0372             break;
0373           case 1:
0374             hTmax_2t0_1->Fill(tmax134);
0375             break;
0376           case 2:
0377             hTmax_2t0_2->Fill(tmax134);
0378             break;
0379           case 3:
0380             hTmax_2t0_3->Fill(tmax134);
0381             break;
0382           case 4:
0383             hTmax_2t0_4->Fill(tmax134);
0384             break;
0385           case 99:
0386             hTmax_2t0_5->Fill(tmax134);
0387             break;
0388         }
0389       } else if (t0_134 == 3) {
0390         hTmax_3t0->Fill(tmax134);
0391         switch (hSubGroup) {
0392           case 0:
0393             hTmax_3t0_0->Fill(tmax134);
0394             break;
0395           case 1:
0396             hTmax_3t0_1->Fill(tmax134);
0397             break;
0398           case 2:
0399             hTmax_3t0_2->Fill(tmax134);
0400             break;
0401           case 3:
0402             hTmax_3t0_3->Fill(tmax134);
0403             break;
0404           case 4:
0405             hTmax_3t0_4->Fill(tmax134);
0406             break;
0407           case 99:
0408             hTmax_3t0_5->Fill(tmax134);
0409             break;
0410         }
0411       }
0412     }
0413     if (tmax234 > 0.) {
0414       hTmax234->Fill(tmax234);
0415       if (t0_234 == 1) {
0416         hTmax_t0->Fill(tmax234);
0417         switch (hSubGroup) {
0418           case 0:
0419             hTmax_t0_0->Fill(tmax234);
0420             break;
0421           case 1:
0422             hTmax_t0_1->Fill(tmax234);
0423             break;
0424           case 2:
0425             hTmax_t0_2->Fill(tmax234);
0426             break;
0427           case 3:
0428             hTmax_t0_3->Fill(tmax234);
0429             break;
0430           case 4:
0431             hTmax_t0_4->Fill(tmax234);
0432             break;
0433           case 99:
0434             hTmax_t0_5->Fill(tmax234);
0435             break;
0436         }
0437       } else {
0438         hTmax_2t0->Fill(tmax234);
0439         switch (hSubGroup) {
0440           case 0:
0441             hTmax_2t0_0->Fill(tmax234);
0442             break;
0443           case 1:
0444             hTmax_2t0_1->Fill(tmax234);
0445             break;
0446           case 2:
0447             hTmax_2t0_2->Fill(tmax234);
0448             break;
0449           case 3:
0450             hTmax_2t0_3->Fill(tmax234);
0451             break;
0452           case 4:
0453             hTmax_2t0_4->Fill(tmax234);
0454             break;
0455           case 99:
0456             hTmax_2t0_5->Fill(tmax234);
0457             break;
0458         }
0459       }
0460     }
0461   }
0462 
0463   void Write() {
0464     // write the Tmax histograms
0465     hTmax123->Write();
0466     hTmax124s72->Write();
0467     hTmax124s78->Write();
0468     hTmax134s72->Write();
0469     hTmax134s78->Write();
0470     hTmax234->Write();
0471     hTmax_3t0->Write();
0472     hTmax_3t0_0->Write();
0473     hTmax_3t0_1->Write();
0474     hTmax_3t0_2->Write();
0475     hTmax_3t0_3->Write();
0476     hTmax_3t0_4->Write();
0477     hTmax_3t0_5->Write();
0478     hTmax_2t0->Write();
0479     hTmax_2t0_0->Write();
0480     hTmax_2t0_1->Write();
0481     hTmax_2t0_2->Write();
0482     hTmax_2t0_3->Write();
0483     hTmax_2t0_4->Write();
0484     hTmax_2t0_5->Write();
0485     hTmax_t0->Write();
0486     hTmax_t0_0->Write();
0487     hTmax_t0_1->Write();
0488     hTmax_t0_2->Write();
0489     hTmax_t0_3->Write();
0490     hTmax_t0_4->Write();
0491     hTmax_t0_5->Write();
0492     hTmax_0->Write();
0493   }
0494 
0495   int GetT0Factor(TH1F *hist) {
0496     unsigned t0 = 999;
0497 
0498     if (hist == hTmax_3t0)
0499       t0 = 3;
0500     else if (hist == hTmax_2t0)
0501       t0 = 2;
0502     else if (hist == hTmax_t0)
0503       t0 = 1;
0504     else if (hist == hTmax_0)
0505       t0 = 0;
0506 
0507     return t0;
0508   }
0509 
0510   TH1F *hTmax123;
0511   TH1F *hTmax124s72;
0512   TH1F *hTmax124s78;
0513   TH1F *hTmax134s72;
0514   TH1F *hTmax134s78;
0515   TH1F *hTmax234;
0516   TH1F *hTmax_3t0;
0517   TH1F *hTmax_3t0_0;
0518   TH1F *hTmax_3t0_1;
0519   TH1F *hTmax_3t0_2;
0520   TH1F *hTmax_3t0_3;
0521   TH1F *hTmax_3t0_4;
0522   TH1F *hTmax_3t0_5;
0523   TH1F *hTmax_2t0;
0524   TH1F *hTmax_2t0_0;
0525   TH1F *hTmax_2t0_1;
0526   TH1F *hTmax_2t0_2;
0527   TH1F *hTmax_2t0_3;
0528   TH1F *hTmax_2t0_4;
0529   TH1F *hTmax_2t0_5;
0530   TH1F *hTmax_t0;
0531   TH1F *hTmax_t0_0;
0532   TH1F *hTmax_t0_1;
0533   TH1F *hTmax_t0_2;
0534   TH1F *hTmax_t0_3;
0535   TH1F *hTmax_t0_4;
0536   TH1F *hTmax_t0_5;
0537   TH1F *hTmax_0;
0538 
0539   TString name;
0540 };
0541 
0542 #endif