/[PAMELA software]/eventviewer/flight/src/FEVdetector.cpp
ViewVC logotype

Diff of /eventviewer/flight/src/FEVdetector.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.24 by mocchiut, Wed Nov 14 16:19:41 2007 UTC revision 1.25 by mocchiut, Thu Nov 29 15:53:10 2007 UTC
# Line 1448  void FEVdetector::DisplayEvent(){ Line 1448  void FEVdetector::DisplayEvent(){
1448        text->DrawLatex(var.xxvc+0.17*var.sfx,var.yxvc-0.08*var.sfy,"CAS1");        text->DrawLatex(var.xxvc+0.17*var.sfx,var.yxvc-0.08*var.sfy,"CAS1");
1449        text->DrawLatex(var.xyvc-0.192*var.sfx,var.yyvc-0.08*var.sfy,"CAS4");        text->DrawLatex(var.xyvc-0.192*var.sfx,var.yyvc-0.08*var.sfy,"CAS4");
1450        text->DrawLatex(var.xyvc+0.17*var.sfx,var.yyvc-0.08*var.sfy,"CAS3");        text->DrawLatex(var.xyvc+0.17*var.sfx,var.yyvc-0.08*var.sfy,"CAS3");
1451        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.48*var.sfy,"CARD2");        text->DrawLatex(var.xxvc-0.222*var.sfx,var.yxvc+0.48*var.sfy,"CARD2");
1452        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.48*var.sfy,"CARD3");        text->DrawLatex(var.xxvc+0.185*var.sfx,var.yxvc+0.48*var.sfy,"CARD3");
1453        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.48*var.sfy,"CARD1");        text->DrawLatex(var.xyvc-0.190*var.sfx,var.yyvc+0.48*var.sfy,"CARD1");
1454        text->DrawLatex(var.xyvc+0.13*var.sfx,var.yyvc+0.48*var.sfy,"CARD4");        text->DrawLatex(var.xyvc+0.16*var.sfx,var.yyvc+0.48*var.sfy,"CARD4");
1455        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.39*var.sfy,"CAT3");        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.39*var.sfy,"CAT3");
1456        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.39*var.sfy,"CAT4");        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.39*var.sfy,"CAT4");
1457        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.39*var.sfy,"CAT2");        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.39*var.sfy,"CAT2");
# Line 1913  void FEVdetector::ShowTOF(){ Line 1913  void FEVdetector::ShowTOF(){
1913    if ( !var.TOF ) return;    if ( !var.TOF ) return;
1914    //    //
1915    Int_t tt = 0;    Int_t tt = 0;
1916      Int_t rj = 0;
1917    //    //
1918    // force level0, i.e. show patterntrig variable!    // force level0, i.e. show patterntrig variable!
1919    //    //
# Line 2106  void FEVdetector::ShowTOF(){ Line 2107  void FEVdetector::ShowTOF(){
2107      Float_t s12p = 0.055;      Float_t s12p = 0.055;
2108      ocolo = 10;      ocolo = 10;
2109      TPolyLine *ftof12ya[6];      TPolyLine *ftof12ya[6];
2110        rj=6;
2111      for ( Int_t j=0; j<6; j++){      for ( Int_t j=0; j<6; j++){
2112          rj--;
2113        Float_t xc1[4]={ (s12p*j),       s12p*(j),               s12p*(j+1),     (s12p*j)};        Float_t xc1[4]={ (s12p*j),       s12p*(j),               s12p*(j+1),     (s12p*j)};
2114        Float_t yc1[4]={      -ws13,     0.,  0., -ws13};        Float_t yc1[4]={      -ws13,     0.,  0., -ws13};
2115        Float_t xc[4];            Float_t xc[4];    
# Line 2122  void FEVdetector::ShowTOF(){ Line 2125  void FEVdetector::ShowTOF(){
2125        } else {        } else {
2126          colo = 1;          colo = 1;
2127        };        };
2128        ColorMIP(ms12a[j],colo);        ColorMIP(ms12a[rj],colo);
2129        //        //
2130        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2131        if ( mt12[0][j] == 4095. ) colo = noadc;        if ( mt12[0][rj] == 4095. ) colo = noadc;
2132        //        //
2133        ftof12ya[j]->SetFillColor(colo);        ftof12ya[j]->SetFillColor(colo);
2134        ftof12ya[j]->SetLineWidth(1);        ftof12ya[j]->SetLineWidth(1);
# Line 2135  void FEVdetector::ShowTOF(){ Line 2138  void FEVdetector::ShowTOF(){
2138      // B      // B
2139      s12p = 0.055;      s12p = 0.055;
2140      TPolyLine *ftof12yb[6];      TPolyLine *ftof12yb[6];
2141        rj=6;
2142      for ( Int_t j=0; j<6; j++){      for ( Int_t j=0; j<6; j++){
2143          rj--;
2144        Float_t xc1[4]={ (s12p*j),       s12p*(j+1),               s12p*(j+1),     (s12p*j)};        Float_t xc1[4]={ (s12p*j),       s12p*(j+1),               s12p*(j+1),     (s12p*j)};
2145        Float_t yc1[4]={ -ws13,     -ws13,    0., -ws13};        Float_t yc1[4]={ -ws13,     -ws13,    0., -ws13};
2146        Float_t xc[4];            Float_t xc[4];    
# Line 2151  void FEVdetector::ShowTOF(){ Line 2156  void FEVdetector::ShowTOF(){
2156        } else {        } else {
2157          colo = 1;          colo = 1;
2158        };        };
2159        ColorMIP(ms12b[j],colo);        ColorMIP(ms12b[rj],colo);
2160        //        //
2161        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2162        if ( mt12[1][j] == 4095. ) colo = noadc;        if ( mt12[1][rj] == 4095. ) colo = noadc;
2163        //        //
2164        ftof12yb[j]->SetFillColor(colo);        ftof12yb[j]->SetFillColor(colo);
2165        ftof12yb[j]->SetLineWidth(1);        ftof12yb[j]->SetLineWidth(1);
# Line 2188  void FEVdetector::ShowTOF(){ Line 2193  void FEVdetector::ShowTOF(){
2193      Float_t s21p = 0.075;      Float_t s21p = 0.075;
2194      ocolo = 10;      ocolo = 10;
2195      TPolyLine *ftof21ya[2];      TPolyLine *ftof21ya[2];
2196        rj=2;
2197      for ( Int_t j=0; j<2; j++){      for ( Int_t j=0; j<2; j++){
2198          rj--;
2199        Float_t xc1[4]={ s21p*(j-1),       s21p*(j-1),               s21p*j,     s21p*(j-1)};        Float_t xc1[4]={ s21p*(j-1),       s21p*(j-1),               s21p*j,     s21p*(j-1)};
2200        Float_t yc1[4]={ 0.,                ws2,     ws2,    0.};        Float_t yc1[4]={ 0.,                ws2,     ws2,    0.};
2201        Float_t xc[4];            Float_t xc[4];    
# Line 2204  void FEVdetector::ShowTOF(){ Line 2211  void FEVdetector::ShowTOF(){
2211        } else {        } else {
2212          colo = 1;          colo = 1;
2213        };        };
2214        ColorMIP(ms21a[j],colo);        ColorMIP(ms21a[rj],colo);
2215        //        //
2216        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2217        if ( mt21[0][j] == 4095. ) colo = noadc;        if ( mt21[0][rj] == 4095. ) colo = noadc;
2218        //        //
2219        ftof21ya[j]->SetFillColor(colo);        ftof21ya[j]->SetFillColor(colo);
2220        ftof21ya[j]->SetLineWidth(1);        ftof21ya[j]->SetLineWidth(1);
# Line 2217  void FEVdetector::ShowTOF(){ Line 2224  void FEVdetector::ShowTOF(){
2224      // B      // B
2225      s21p = 0.075;      s21p = 0.075;
2226      TPolyLine *ftof21yb[2];      TPolyLine *ftof21yb[2];
2227        rj=2;
2228      for ( Int_t j=0; j<2; j++){      for ( Int_t j=0; j<2; j++){
2229          rj--;
2230        Float_t xc1[4]={ s21p*(j-1),       s21p*j,               s21p*j,     s21p*(j-1)};        Float_t xc1[4]={ s21p*(j-1),       s21p*j,               s21p*j,     s21p*(j-1)};
2231        Float_t yc1[4]={ 0.,         0.,                ws2,         0.};        Float_t yc1[4]={ 0.,         0.,                ws2,         0.};
2232        Float_t xc[4];            Float_t xc[4];    
# Line 2233  void FEVdetector::ShowTOF(){ Line 2242  void FEVdetector::ShowTOF(){
2242        } else {        } else {
2243          colo = 1;          colo = 1;
2244        };        };
2245        ColorMIP(ms21b[j],colo);        ColorMIP(ms21b[rj],colo);
2246        //        //
2247        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2248        if ( mt21[1][j] == 4095. ) colo = noadc;        if ( mt21[1][rj] == 4095. ) colo = noadc;
2249        //        //
2250        ftof21yb[j]->SetFillColor(colo);        ftof21yb[j]->SetFillColor(colo);
2251        ftof21yb[j]->SetLineWidth(1);        ftof21yb[j]->SetLineWidth(1);
# Line 2351  void FEVdetector::ShowTOF(){ Line 2360  void FEVdetector::ShowTOF(){
2360      Float_t s31p = 0.060;      Float_t s31p = 0.060;
2361      ocolo = 10;      ocolo = 10;
2362      TPolyLine *ftof31xa[3];      TPolyLine *ftof31xa[3];
2363        rj=3;
2364      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2365          rj--;
2366        Float_t xc1[4]={ (s31p*j),       (s31p*j),               s31p*(j+1),       (s31p*j)};        Float_t xc1[4]={ (s31p*j),       (s31p*j),               s31p*(j+1),       (s31p*j)};
2367        Float_t yc1[4]={ 0.,                       ws13,     ws13,    0.};        Float_t yc1[4]={ 0.,                       ws13,     ws13,    0.};
2368        Float_t xc[4];            Float_t xc[4];    
# Line 2367  void FEVdetector::ShowTOF(){ Line 2378  void FEVdetector::ShowTOF(){
2378        } else {        } else {
2379          colo = 1;          colo = 1;
2380        };        };
2381        ColorMIP(ms31a[j],colo);        ColorMIP(ms31a[rj],colo);
2382        //        //
2383        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2384        if ( mt31[0][j] == 4095. ) colo = noadc;        if ( mt31[0][rj] == 4095. ) colo = noadc;
2385        //        //
2386        ftof31xa[j]->SetFillColor(colo);        ftof31xa[j]->SetFillColor(colo);
2387        ftof31xa[j]->SetLineWidth(1);        ftof31xa[j]->SetLineWidth(1);
# Line 2380  void FEVdetector::ShowTOF(){ Line 2391  void FEVdetector::ShowTOF(){
2391      // B      // B
2392      s31p = 0.060;      s31p = 0.060;
2393      TPolyLine *ftof31xb[3];      TPolyLine *ftof31xb[3];
2394        rj=3;
2395      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2396          rj--;
2397        Float_t xc1[4]={ (s31p*j),       s31p*(j+1),               s31p*(j+1),     (s31p*j)};        Float_t xc1[4]={ (s31p*j),       s31p*(j+1),               s31p*(j+1),     (s31p*j)};
2398        Float_t yc1[4]={ 0.,                        0.,     ws13,    0.};        Float_t yc1[4]={ 0.,                        0.,     ws13,    0.};
2399        Float_t xc[4];            Float_t xc[4];    
# Line 2396  void FEVdetector::ShowTOF(){ Line 2409  void FEVdetector::ShowTOF(){
2409        } else {        } else {
2410          colo = 1;          colo = 1;
2411        };        };
2412        ColorMIP(ms31b[j],colo);        ColorMIP(ms31b[rj],colo);
2413        //        //
2414        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2415        if ( mt31[1][j] == 4095. ) colo = noadc;        if ( mt31[1][rj] == 4095. ) colo = noadc;
2416        //        //
2417        ftof31xb[j]->SetFillColor(colo);        ftof31xb[j]->SetFillColor(colo);
2418        ftof31xb[j]->SetLineWidth(1);        ftof31xb[j]->SetLineWidth(1);
# Line 2431  void FEVdetector::ShowTOF(){ Line 2444  void FEVdetector::ShowTOF(){
2444      Float_t s32p = 0.050;      Float_t s32p = 0.050;
2445      ocolo = 10;      ocolo = 10;
2446      TPolyLine *ftof32ya[3];      TPolyLine *ftof32ya[3];
2447        rj=3;
2448      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2449          rj--;
2450        Float_t xc1[4]={ (s32p*j),       (s32p*j),               s32p*(j+1),     (s32p*j)};        Float_t xc1[4]={ (s32p*j),       (s32p*j),               s32p*(j+1),     (s32p*j)};
2451        Float_t yc1[4]={     -ws13,     0.,    0., -ws13};        Float_t yc1[4]={     -ws13,     0.,    0., -ws13};
2452        Float_t xc[4];            Float_t xc[4];    
# Line 2447  void FEVdetector::ShowTOF(){ Line 2462  void FEVdetector::ShowTOF(){
2462        } else {        } else {
2463          colo = 1;          colo = 1;
2464        };        };
2465        ColorMIP(ms32a[j],colo);        ColorMIP(ms32a[rj],colo);
2466        //        //
2467        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2468        if ( mt32[0][j] == 4095. ) colo = noadc;        if ( mt32[0][rj] == 4095. ) colo = noadc;
2469        //        //
2470        ftof32ya[j]->SetFillColor(colo);        ftof32ya[j]->SetFillColor(colo);
2471        ftof32ya[j]->SetLineWidth(1);        ftof32ya[j]->SetLineWidth(1);
# Line 2459  void FEVdetector::ShowTOF(){ Line 2474  void FEVdetector::ShowTOF(){
2474      };      };
2475      s32p = 0.050;      s32p = 0.050;
2476      TPolyLine *ftof32yb[3];      TPolyLine *ftof32yb[3];
2477        rj=3;
2478      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2479          rj--;
2480        Float_t xc1[4]={ (s32p*j),       s32p*(j+1),               s32p*(j+1),     (s32p*j)};        Float_t xc1[4]={ (s32p*j),       s32p*(j+1),               s32p*(j+1),     (s32p*j)};
2481        Float_t yc1[4]={ -ws13,   -ws13,    0., -ws13};        Float_t yc1[4]={ -ws13,   -ws13,    0., -ws13};
2482        Float_t xc[4];            Float_t xc[4];    
# Line 2475  void FEVdetector::ShowTOF(){ Line 2492  void FEVdetector::ShowTOF(){
2492        } else {        } else {
2493          colo = 1;          colo = 1;
2494        };        };
2495        ColorMIP(ms32b[j],colo);        ColorMIP(ms32b[rj],colo);
2496        //        //
2497        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2498        if ( mt32[1][j] == 4095. ) colo = noadc;        if ( mt32[1][rj] == 4095. ) colo = noadc;
2499        //        //
2500        ftof32yb[j]->SetFillColor(colo);        ftof32yb[j]->SetFillColor(colo);
2501        ftof32yb[j]->SetLineWidth(1);        ftof32yb[j]->SetLineWidth(1);
# Line 2922  void FEVdetector::ShowTOF(){ Line 2939  void FEVdetector::ShowTOF(){
2939      ftof11y->Draw();      ftof11y->Draw();
2940      TPolyLine *sftof11[8];      TPolyLine *sftof11[8];
2941      tt = 1;      tt = 1;
2942      if ( level.file == -1 ) tt = 8;      rj=1;
2943        if ( level.file == -1 ){
2944          tt = 8;
2945          rj = 8;
2946        };
2947      for ( Int_t j=0; j<tt; j++){      for ( Int_t j=0; j<tt; j++){
2948        if ( (mt11[0][j]<4095. || mt11[1][j]<4095.)){        rj--;
2949          if ( (mt11[0][rj]<4095. || mt11[1][rj]<4095.)){
2950          if ( ocolo == 10 ) ocolo = noadc;          if ( ocolo == 10 ) ocolo = noadc;
2951          Float_t lowp = (xp11[j]-E11[j])/100.;          Float_t lowp = -(xp11[rj]+E11[rj])/100.;
2952          Float_t higp = (xp11[j]+E11[j])/100.;          Float_t higp = -(xp11[rj]-E11[rj])/100.;
2953          if ( lowp < -0.164 ) lowp = -0.164;          if ( lowp < -0.164 ) lowp = -0.164;
2954          if ( higp > 0.164 ) higp = 0.164;          if ( higp > 0.164 ) higp = 0.164;
2955          if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp11[j] != 0.))){          if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp11[rj] != 0.))){
2956            Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};            Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
2957            Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};            Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};
2958            Float_t nxc[5];                Float_t nxc[5];    
# Line 2955  void FEVdetector::ShowTOF(){ Line 2977  void FEVdetector::ShowTOF(){
2977      Float_t s12p = 0.055;      Float_t s12p = 0.055;
2978      ocolo = 10;      ocolo = 10;
2979      TPolyLine *ftof12y[6];      TPolyLine *ftof12y[6];
2980        rj = 6;
2981      for ( Int_t j=0; j<6; j++){      for ( Int_t j=0; j<6; j++){
2982          rj--;
2983        Float_t nxc1[5]={ (s12p*j),       (s12p*j),               s12p*(j+1), s12p*(j+1),    (s12p*j)};        Float_t nxc1[5]={ (s12p*j),       (s12p*j),               s12p*(j+1), s12p*(j+1),    (s12p*j)};
2984        Float_t nyc1[5]={      -ws13,     0.,  0., -ws13, -ws13};        Float_t nyc1[5]={      -ws13,     0.,  0., -ws13, -ws13};
2985        Float_t nxc[5];            Float_t nxc[5];    
# Line 2972  void FEVdetector::ShowTOF(){ Line 2996  void FEVdetector::ShowTOF(){
2996          colo = 1;          colo = 1;
2997        };        };
2998        if ( level.file == -1 ){        if ( level.file == -1 ){
2999          ColorMIP(ms12a[j]+ms12b[j],colo);          ColorMIP(ms12a[rj]+ms12b[rj],colo);
3000        } else {        } else {
3001          ColorTOFMIP(ms12a[j]+ms12b[j],colo);          ColorTOFMIP(ms12a[rj]+ms12b[rj],colo);
3002        };        };
3003        //        //
3004        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
# Line 3041  void FEVdetector::ShowTOF(){ Line 3065  void FEVdetector::ShowTOF(){
3065      Float_t s21p = 0.075;      Float_t s21p = 0.075;
3066      ocolo = 10;      ocolo = 10;
3067      TPolyLine *ftof21y[2];      TPolyLine *ftof21y[2];
3068        rj=2;
3069      for ( Int_t j=0; j<2; j++){      for ( Int_t j=0; j<2; j++){
3070          rj--;
3071        Float_t nxc1[5]={ s21p*(j-1),       s21p*(j-1),               s21p*j,s21p*j,     s21p*(j-1)};        Float_t nxc1[5]={ s21p*(j-1),       s21p*(j-1),               s21p*j,s21p*j,     s21p*(j-1)};
3072        Float_t nyc1[5]={ 0.,                ws2,     ws2,  0.,  0.};        Float_t nyc1[5]={ 0.,                ws2,     ws2,  0.,  0.};
3073        Float_t nxc[5];            Float_t nxc[5];    
# Line 3058  void FEVdetector::ShowTOF(){ Line 3084  void FEVdetector::ShowTOF(){
3084          colo = 1;          colo = 1;
3085        };        };
3086        if ( level.file == -1 ){        if ( level.file == -1 ){
3087          ColorMIP(ms21a[j]+ms21b[j],colo);          ColorMIP(ms21a[rj]+ms21b[rj],colo);
3088        } else {        } else {
3089          ColorTOFMIP(ms21a[j]+ms21b[j],colo);          ColorTOFMIP(ms21a[rj]+ms21b[rj],colo);
3090        };        };
3091        //        //
3092        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
# Line 3177  void FEVdetector::ShowTOF(){ Line 3203  void FEVdetector::ShowTOF(){
3203        ftof22y->Draw();        ftof22y->Draw();
3204        TPolyLine *sftof22[2];        TPolyLine *sftof22[2];
3205        tt = 1;        tt = 1;
3206        if ( level.file == -1 ) tt = 2;        rj = 1;
3207          if ( level.file == -1 ){
3208            tt = 2;
3209            rj = 2;
3210          };
3211        for (Int_t j=0; j<tt; j++){        for (Int_t j=0; j<tt; j++){
3212            rj--;
3213          //    for ( Int_t j=0; j<2; j++){          //    for ( Int_t j=0; j<2; j++){
3214          if ( (mt22[0][j]<4095. || mt22[1][j]<4095.) ){          if ( (mt22[0][rj]<4095. || mt22[1][rj]<4095.) ){
3215            if ( ocolo == 10 ) ocolo = noadc;            if ( ocolo == 10 ) ocolo = noadc;
3216            Float_t lowp = (xp22[j]-E22[j])/100.;            Float_t lowp = -(xp22[rj]+E22[rj])/100.;
3217            Float_t higp = (xp22[j]+E22[j])/100.;            Float_t higp = -(xp22[rj]-E22[rj])/100.;
3218            if ( lowp < -0.074 ) lowp = -0.074;            if ( lowp < -0.074 ) lowp = -0.074;
3219            if ( higp > 0.074 ) higp = 0.074;            if ( higp > 0.074 ) higp = 0.074;
3220            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp22[j] != 0.))){            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp22[rj] != 0.))){
3221              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
3222              Float_t nyc1[5]={ -0.0015,         -0.0015,                -ws2+0.0015,    -ws2+0.0015,    -0.0015};              Float_t nyc1[5]={ -0.0015,         -0.0015,                -ws2+0.0015,    -ws2+0.0015,    -0.0015};
3223              Float_t nxc[5];                  Float_t nxc[5];    
# Line 3206  void FEVdetector::ShowTOF(){ Line 3237  void FEVdetector::ShowTOF(){
3237        };        };
3238      };      };
3239      //      //
3240        //daquiinbasso
3241      //      //
3242      // S31  X-view      // S31  X-view
3243      //      //
# Line 3263  void FEVdetector::ShowTOF(){ Line 3294  void FEVdetector::ShowTOF(){
3294        ftof31y->Draw();        ftof31y->Draw();
3295        TPolyLine *sftof31[3];        TPolyLine *sftof31[3];
3296        tt = 1;        tt = 1;
3297        if ( level.file == -1 ) tt = 3;        rj = 1;
3298          if ( level.file == -1 ){
3299            tt = 3;
3300            rj = 3;
3301          };
3302        for (Int_t j=0; j<tt; j++){        for (Int_t j=0; j<tt; j++){
3303            rj--;
3304          //    for ( Int_t j=0; j<3; j++){          //    for ( Int_t j=0; j<3; j++){
3305          if ( (mt31[0][j]<4095. || mt31[1][j]<4095.)){          if ( (mt31[0][rj]<4095. || mt31[1][rj]<4095.)){
3306            if ( ocolo == 10 ) ocolo = noadc;            if ( ocolo == 10 ) ocolo = noadc;
3307            Float_t lowp = (xp31[j]-E31[j])/100.;            Float_t lowp = -(xp31[rj]+E31[rj])/100.;
3308            Float_t higp = (xp31[j]+E31[j])/100.;            Float_t higp = -(xp31[rj]-E31[rj])/100.;
3309            if ( lowp < -0.074 ) lowp = -0.074;            if ( lowp < -0.074 ) lowp = -0.074;
3310            if ( higp > 0.074 ) higp = 0.074;            if ( higp > 0.074 ) higp = 0.074;
3311            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp31[j] != 0.))){            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp31[rj] != 0.))){
3312              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
3313              Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};              Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};
3314              Float_t nxc[5];                  Float_t nxc[5];    
# Line 3299  void FEVdetector::ShowTOF(){ Line 3335  void FEVdetector::ShowTOF(){
3335      Float_t s32p = 0.050;      Float_t s32p = 0.050;
3336      ocolo = 10;      ocolo = 10;
3337      TPolyLine *ftof32y[3];      TPolyLine *ftof32y[3];
3338        rj = 3;
3339      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
3340          rj--;
3341        Float_t nxc1[5]={ (s32p*j),       (s32p*j),               s32p*(j+1),s32p*(j+1),     (s32p*j)};        Float_t nxc1[5]={ (s32p*j),       (s32p*j),               s32p*(j+1),s32p*(j+1),     (s32p*j)};
3342        Float_t nyc1[5]={     -ws13,     0.,    0., -ws13, -ws13};        Float_t nyc1[5]={     -ws13,     0.,    0., -ws13, -ws13};
3343        Float_t nxc[5];            Float_t nxc[5];    
# Line 3316  void FEVdetector::ShowTOF(){ Line 3354  void FEVdetector::ShowTOF(){
3354          colo = 1;          colo = 1;
3355        };        };
3356        if ( level.file == -1 ){        if ( level.file == -1 ){
3357          ColorMIP(ms32a[j]+ms32b[j],colo);          ColorMIP(ms32a[rj]+ms32b[rj],colo);
3358        } else {        } else {
3359          ColorTOFMIP(ms32a[j]+ms32b[j],colo);          ColorTOFMIP(ms32a[rj]+ms32b[rj],colo);
3360        };        };
3361        //        //
3362        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
# Line 5032  void FEVdetector::ShowTRK(Bool_t upd){ Line 5070  void FEVdetector::ShowTRK(Bool_t upd){
5070            for (Int_t m = 0; m<3; m++){            for (Int_t m = 0; m<3; m++){
5071              if ( trk->signcluster[l][m] != 0. ){              if ( trk->signcluster[l][m] != 0. ){
5072                if ( planeno == 10 ) {                if ( planeno == 10 ) {
                 x = 1024 - trk->addrcluster[l][m];  
               } else {  
5073                  x = trk->addrcluster[l][m];                  x = trk->addrcluster[l][m];
5074                  } else {
5075                    x = 1024 - trk->addrcluster[l][m];
5076                };                };
5077                linea = new TLine(14.+x,0.01,14.+x,0.3);                linea = new TLine(14.+x,0.01,14.+x,0.3);
5078                linea->SetLineWidth(2);                linea->SetLineWidth(2);
# Line 5104  void FEVdetector::ShowTRK(Bool_t upd){ Line 5142  void FEVdetector::ShowTRK(Bool_t upd){
5142          trkpad[planepad]->cd();          trkpad[planepad]->cd();
5143          trkpad[planepad]->Range(-7.05,0.,7.05,8.);          trkpad[planepad]->Range(-7.05,0.,7.05,8.);
5144          for (Int_t img=0; img<2; img++){          for (Int_t img=0; img<2; img++){
5145            x = singlet->coord[img];            x = -singlet->coord[img];
5146            Float_t xsig = singlet->sgnl;            Float_t xsig = singlet->sgnl;
5147            if ( xsig > 8. ) xsig = 8.;            if ( xsig > 8. ) xsig = 8.;
5148            //            //
# Line 5189  void FEVdetector::ShowTRK(Bool_t upd){ Line 5227  void FEVdetector::ShowTRK(Bool_t upd){
5227              planepad = plane * 2;              planepad = plane * 2;
5228              trkpad[planepad]->cd();              trkpad[planepad]->cd();
5229              trkpad[planepad]->Range(-7.05,0.,7.05,8.);              trkpad[planepad]->Range(-7.05,0.,7.05,8.);
5230              x = track->ym[plane];              x = -track->ym[plane];
5231              xsig = track->dedx_y[plane];              xsig = track->dedx_y[plane];
5232              if ( track->YGood(plane) ){              if ( track->YGood(plane) ){
5233                //            if ( x > -100 ){                //            if ( x > -100 ){
# Line 5251  void FEVdetector::ShowTRK(Bool_t upd){ Line 5289  void FEVdetector::ShowTRK(Bool_t upd){
5289            for (Int_t plane = 0; plane<6; plane++){            for (Int_t plane = 0; plane<6; plane++){
5290              //              //
5291              x = track->xm[plane];              x = track->xm[plane];
5292              Float_t y = track->ym[plane];              Float_t y = -track->ym[plane];
5293              Float_t xsig = track->dedx_x[plane];              Float_t xsig = track->dedx_x[plane];
5294              Float_t ysig = track->dedx_y[plane];              Float_t ysig = track->dedx_y[plane];
5295              //    if ( x > -100 && y > -100. ){              //    if ( x > -100 && y > -100. ){
# Line 5373  void FEVdetector::ShowTRK(Bool_t upd){ Line 5411  void FEVdetector::ShowTRK(Bool_t upd){
5411              for (Int_t sing = 0; sing < min(totsx,totsy); sing++){              for (Int_t sing = 0; sing < min(totsx,totsy); sing++){
5412                //                              //              
5413                x = arrax->At(orderx->At(sing));                x = arrax->At(orderx->At(sing));
5414                Float_t y1 = array1->At(ordery->At(sing));                Float_t y1 = -array1->At(ordery->At(sing));
5415                Float_t y2 = array2->At(ordery->At(sing));                Float_t y2 = -array2->At(ordery->At(sing));
5416                Float_t xsig = sigarrax->At(orderx->At(sing));                Float_t xsig = sigarrax->At(orderx->At(sing));
5417                Float_t ysig = sigarray->At(ordery->At(sing));                Float_t ysig = sigarray->At(ordery->At(sing));
5418                if ( x > -100 && y1 > -100. ){                if ( x > -100 && y1 > -100. ){
# Line 5467  void FEVdetector::ShowTRK(Bool_t upd){ Line 5505  void FEVdetector::ShowTRK(Bool_t upd){
5505                  for (Int_t m = 0; m<3; m++){                                          for (Int_t m = 0; m<3; m++){                        
5506                    if ( trk->signcluster[l][m] != 0. ){                    if ( trk->signcluster[l][m] != 0. ){
5507                      if ( planeno == 10 ) {                      if ( planeno == 10 ) {
                       x = 1024 - trk->addrcluster[l][m];  
                     } else {  
5508                        x = trk->addrcluster[l][m];                        x = trk->addrcluster[l][m];
5509                        } else {
5510                          x = 1024 - trk->addrcluster[l][m];
5511                      };                      };
5512                      x += 1024.;                                      x += 1024.;                
5513                      if ( trk->signcluster[l][m]> ysig  ){                      if ( trk->signcluster[l][m]> ysig  ){
5514                        ysig = trk->signcluster[l][m];                        ysig = trk->signcluster[l][m];
5515                        yh1[ncro] = x-1024.;                        yh1[ncro] = x - 1024.;
5516                        yh2[ncro] = x;                        yh2[ncro] = x;
5517                      };                      };
5518                    };                    };
# Line 5502  void FEVdetector::ShowTRK(Bool_t upd){ Line 5540  void FEVdetector::ShowTRK(Bool_t upd){
5540                    for (Int_t m = 0; m<3; m++){                    for (Int_t m = 0; m<3; m++){
5541                      if ( trk->signcluster[e][m] != 0. ){                      if ( trk->signcluster[e][m] != 0. ){
5542                        if ( planeno == 10 ) {                        if ( planeno == 10 ) {
                         x = 1024 - trk->addrcluster[e][m];  
                       } else {  
5543                          x = trk->addrcluster[e][m];                          x = trk->addrcluster[e][m];
5544                          } else {
5545                            x = 1024 - trk->addrcluster[e][m];
5546                        };                        };
5547                        x += 1024.;                                              x += 1024.;                      
5548                        if ( trk->signcluster[e][m]> ysig  ){                        if ( trk->signcluster[e][m]> ysig  ){
5549                          ysig = trk->signcluster[e][m];                          ysig = trk->signcluster[e][m];
5550                          yh1[ncro] = x-1024.;                          yh1[ncro] = x - 1024.;
5551                          yh2[ncro] = x;                          yh2[ncro] = x;
5552                        };                        };
5553                      };                      };
# Line 5582  void FEVdetector::ShowTrack(){ Line 5620  void FEVdetector::ShowTrack(){
5620        };        };
5621        for (Int_t plane = 0; plane < 22; plane++){              for (Int_t plane = 0; plane < 22; plane++){      
5622          calx[plane] = (L2->GetCaloLevel2()->cbar[plane][0]/100.)*var.sfx + var.xxvc;;          calx[plane] = (L2->GetCaloLevel2()->cbar[plane][0]/100.)*var.sfx + var.xxvc;;
5623          caly[plane] = (L2->GetCaloLevel2()->cbar[plane][1]/100.)*var.sfx + var.xyvc;          caly[plane] = (-L2->GetCaloLevel2()->cbar[plane][1]/100.)*var.sfx + var.xyvc;
5624          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;
5625          calzx[plane] = zpiano[plane] + var.yxvc;          calzx[plane] = zpiano[plane] + var.yxvc;
5626          calzy[plane] = zpiano[plane] + var.yyvc;          calzy[plane] = zpiano[plane] + var.yyvc;
# Line 5631  void FEVdetector::ShowTrack(){ Line 5669  void FEVdetector::ShowTrack(){
5669        //        //
5670        for (Int_t plane = 0; plane < 22; plane++){              for (Int_t plane = 0; plane < 22; plane++){      
5671          calx[plane] = (ctra->tbar[plane][0]/100.)*var.sfx + var.xxvc;;          calx[plane] = (ctra->tbar[plane][0]/100.)*var.sfx + var.xxvc;;
5672          caly[plane] = (ctra->tbar[plane][1]/100.)*var.sfx + var.xyvc;          caly[plane] = (-ctra->tbar[plane][1]/100.)*var.sfx + var.xyvc;
5673          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;
5674          calzx[plane] = zpiano[plane] + var.yxvc;          calzx[plane] = zpiano[plane] + var.yxvc;
5675          calzy[plane] = zpiano[plane] + var.yyvc;          calzy[plane] = zpiano[plane] + var.yyvc;
# Line 5723  void FEVdetector::ShowTrack(){ Line 5761  void FEVdetector::ShowTrack(){
5761        if ( !ifail ){        if ( !ifail ){
5762          for ( Int_t e = 0; e < npoint ; e++) {          for ( Int_t e = 0; e < npoint ; e++) {
5763            tx[e] = ((float)tr->x[e]/100.)*var.sfx + var.xxvc;            tx[e] = ((float)tr->x[e]/100.)*var.sfx + var.xxvc;
5764            ty[e] = ((float)tr->y[e]/100.)*var.sfx + var.xyvc;            ty[e] = (-(float)tr->y[e]/100.)*var.sfx + var.xyvc;
5765            if ( e < 50 ) {            if ( e < 50 ) {
5766              ptx1[e] = ((float)tr->y[e]/100.)*var.sfx + var.xcat;              ptx1[e] = (-(float)tr->y[e]/100.)*var.sfx + var.xcat;
5767              pty1[e] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;              pty1[e] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;
5768            } else {            } else {
5769              ptx2[e-50] = ((float)tr->y[e]/100.)*var.sfx + var.xcat;              ptx2[e-50] = (-(float)tr->y[e]/100.)*var.sfx + var.xcat;
5770              pty2[e-50] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;              pty2[e-50] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;
5771            };            };
5772            if ( e > 84 ) {            if ( e > 84 ) {
5773              tx1[e-85] = ((float)tr->x[e]/100.)*var.nds4 + var.xxvc;              tx1[e-85] = ((float)tr->x[e]/100.)*var.nds4 + var.xxvc;
5774              ty1[e-85] = ((float)tr->y[e]/100.)*var.nds4 + var.xyvc;                  ty1[e-85] = (-(float)tr->y[e]/100.)*var.nds4 + var.xyvc;    
5775              tzx1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yxvc;              tzx1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yxvc;
5776              tzy1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yyvc;              tzy1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yyvc;
5777            };            };
# Line 5973  void FEVdetector::ShowCalo(Bool_t upd){ Line 6011  void FEVdetector::ShowCalo(Bool_t upd){
6011                Xview->Draw("box same");                Xview->Draw("box same");
6012              };                                            };                              
6013              if ( l == 1 ) {                                      if ( l == 1 ) {                        
6014                Yview->Fill(n,21-m,1.);                Yview->Fill(95-n,21-m,1.);
6015                pd2->cd();                pd2->cd();
6016                Yview->Draw("box same");                Yview->Draw("box same");
6017              };              };
# Line 6202  void FEVdetector::ShowCalo(Bool_t upd){ Line 6240  void FEVdetector::ShowCalo(Bool_t upd){
6240                    Xview->Draw("box same");                    Xview->Draw("box same");
6241                  };                                                    };                                  
6242                  if ( l == 1 ) {                                      if ( l == 1 ) {                    
6243                    Yview->Fill(n,21-m,1.);                    Yview->Fill(96-n,21-m,1.);
6244                    pd2->cd();                    pd2->cd();
6245                    Yview->Draw("box same");                    Yview->Draw("box same");
6246                  };                  };
# Line 6321  void FEVdetector::ShowCalo(Bool_t upd){ Line 6359  void FEVdetector::ShowCalo(Bool_t upd){
6359          gDirectory->Delete(yvev.str().c_str());          gDirectory->Delete(yvev.str().c_str());
6360          TH2F *Yview = new TH2F(yvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);                            TH2F *Yview = new TH2F(yvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);                  
6361          Yview->SetFillColor(colo);          Yview->SetFillColor(colo);
6362          Yview->Fill(n,21-m,1.);          Yview->Fill(96-n,21-m,1.);
6363          pd2->cd();          pd2->cd();
6364          Yview->Draw("box same");          Yview->Draw("box same");
6365        };        };
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.23