- xy[0] = -1.; xy[1] = -0.5; xy[2] = 1.; xy[3] = 2.5;
- ((TVirtualPad*)l->At(1))->cd();
- if(!GetGraphArray(xy, kTrack, 3, 1)) break;
- return kTRUE;
- case 6: // kTrack phi
- xy[0] = -.3; xy[1] = -5.; xy[2] = .3; xy[3] = 50.;
- if(GetGraph(&xy[0], kTrack , 4)) return kTRUE;
- break;
- case 7: // kTrackIn y
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -.3; xy[1] = -1500.; xy[2] = .3; xy[3] = 5000.;
- pad = ((TVirtualPad*)l->At(0)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selStart=0; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kTrackIn, 0, 1, n, selection)) break;
- pad=((TVirtualPad*)l->At(1)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selStart=fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kTrackIn, 0, 1, n, selection)) break;
- return kTRUE;
- case 8: // kTrackIn y
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -.3; xy[1] = -1500.; xy[2] = .3; xy[3] = 5000.;
- pad = ((TVirtualPad*)l->At(0)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selStart=2*fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kTrackIn, 0, 1, n, selection)) break;
- xy[0] = -.5; xy[1] = -0.5; xy[2] = fgkNresYsegm[fSegmentLevel]-.5; xy[3] = 2.5;
- pad=((TVirtualPad*)l->At(1)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- if(!GetGraphArray(xy, kTrackIn, 1, 1)) break;
- return kTRUE;
- case 9: // kTrackIn z
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -1.; xy[1] = -1000.; xy[2] = 1.; xy[3] = 4000.;
- pad = ((TVirtualPad*)l->At(0)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selection[0]=1;
- if(!GetGraphArray(xy, kTrackIn, 2, 1, 1, selection)) break;
- xy[0] = -1.; xy[1] = -0.5; xy[2] = 1.; xy[3] = 2.5;
- pad = ((TVirtualPad*)l->At(1)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- if(!GetGraphArray(xy, kTrackIn, 3, 1)) break;
- return kTRUE;
- case 10: // kTrackIn phi
- xy[0] = -.3; xy[1] = -5.; xy[2] = .3; xy[3] = 50.;
- if(GetGraph(&xy[0], kTrackIn, 4)) return kTRUE;
- break;
- case 11: // kTrackOut y
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -.3; xy[1] = -50.; xy[2] = .3; xy[3] = 150.;
- pad = ((TVirtualPad*)l->At(0)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selStart=0; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kTrackOut, 0, 1, n, selection)) break;
- pad=((TVirtualPad*)l->At(1)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selStart=fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kTrackOut, 0, 1, n, selection)) break;
- return kTRUE;
- case 12: // kTrackOut y
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -.3; xy[1] = -50.; xy[2] = .3; xy[3] = 150.;
- pad = ((TVirtualPad*)l->At(0)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- selStart=2*fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kTrackOut, 0, 1, n, selection)) break;
- xy[0] = -.5; xy[1] = -0.5; xy[2] = fgkNresYsegm[fSegmentLevel]-.5; xy[3] = 2.5;
- pad=((TVirtualPad*)l->At(1)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- if(!GetGraphArray(xy, kTrackOut, 1, 1)) break;
- return kTRUE;
- case 13: // kTrackOut z
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -1.; xy[1] = -1000.; xy[2] = 1.; xy[3] = 4000.;
- pad = ((TVirtualPad*)l->At(0)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- if(!GetGraphArray(xy, kTrackOut, 2, 1)) break;
- xy[0] = -1.; xy[1] = -0.5; xy[2] = 1.; xy[3] = 2.5;
- pad = ((TVirtualPad*)l->At(1)); pad->cd();
- pad->SetMargin(0.1, 0.1, 0.1, 0.01);
- if(!GetGraphArray(xy, kTrackOut, 3, 1)) break;
- return kTRUE;
- case 14: // kTrackOut phi
- xy[0] = -.3; xy[1] = -5.; xy[2] = .3; xy[3] = 50.;
- if(GetGraph(&xy[0], kTrackOut, 4)) return kTRUE;
- break;
- case 15: // kMCcluster
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.3; xy[1]=-50.; xy[2]=.3; xy[3]=650.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=0; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCcluster, 0, 1, n, selection)) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCcluster, 0, 1, n, selection)) break;
- return kTRUE;
- case 16: // kMCcluster
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.3; xy[1]=-50.; xy[2]=.3; xy[3]=650.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=2*fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCcluster, 0, 1, n, selection)) break;
- ((TVirtualPad*)l->At(1))->cd();
- xy[0]=-.5; xy[1]=-0.5; xy[2]=fgkNresYsegm[fSegmentLevel]-.5; xy[3]=2.5;
- if(!GetGraphArray(xy, kMCcluster, 1, 1)) break;
- return kTRUE;
- case 17: //kMCtracklet [y]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.3; xy[1]=-50.; xy[2]=.3; xy[3] =500.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=0; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtracklet, 0, 1, n, selection)) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtracklet, 0, 1, n, selection)) break;
- return kTRUE;
- case 18: //kMCtracklet [y]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.3; xy[1]=-50.; xy[2]=.3; xy[3] =500.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=2*fgkNresYsegm[fSegmentLevel]/3; for(n=0; n<fgkNresYsegm[fSegmentLevel]/3; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtracklet, 0, 1, n, selection)) break;
- ((TVirtualPad*)l->At(1))->cd();
- xy[0]=-.5; xy[1]=-0.5; xy[2]=fgkNresYsegm[fSegmentLevel]-.5; xy[3]=2.5;
- if(!GetGraphArray(xy, kMCtracklet, 1, 1)) break;
- return kTRUE;
- case 19: //kMCtracklet [z]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-1.; xy[1]=-100.; xy[2]=1.; xy[3] =2500.;
- ((TVirtualPad*)l->At(0))->cd();
- if(!GetGraphArray(xy, kMCtracklet, 2)) break;
- xy[0] = -1.; xy[1] = -0.5; xy[2] = 1.; xy[3] = 2.5;
- ((TVirtualPad*)l->At(1))->cd();
- if(!GetGraphArray(xy, kMCtracklet, 3)) break;
- return kTRUE;
- case 20: //kMCtracklet [phi]
- xy[0]=-.3; xy[1]=-3.; xy[2]=.3; xy[3] =25.;
- if(!GetGraph(&xy[0], kMCtracklet, 4)) break;
- return kTRUE;
- case 21: //kMCtrack [y] ly [0]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-50.; xy[2]=.2; xy[3] =400.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*0.); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer1")) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*0.5); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer1")) break;
- return kTRUE;
- case 22: //kMCtrack [y] ly [1]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-50.; xy[2]=.2; xy[3] =400.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*1.); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer2")) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*1.5); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer2")) break;
- return kTRUE;
- case 23: //kMCtrack [y] ly [2]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-50.; xy[2]=.2; xy[3] =400.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*2.); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer3")) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*2.5); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer3")) break;
- return kTRUE;
- case 24: //kMCtrack [y] ly [3]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-50.; xy[2]=.2; xy[3] =400.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*3.); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer4")) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*3.5); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer4")) break;
- return kTRUE;
- case 25: //kMCtrack [y] ly [4]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-50.; xy[2]=.2; xy[3] =400.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*4.); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer5")) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*4.5); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer5")) break;
- return kTRUE;
- case 26: //kMCtrack [y] ly [5]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-50.; xy[2]=.2; xy[3] =400.;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*5.); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer6")) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=Int_t(fgkNresYsegm[fSegmentLevel]*5.5); for(n=0; n<fgkNresYsegm[fSegmentLevel]/2; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 0, 1, n, selection, "Layer6")) break;
- return kTRUE;
- case 27: //kMCtrack [y pulls]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0] = -.5; xy[1] = -0.5; xy[2] = fgkNresYsegm[fSegmentLevel]-.5; xy[3] = 5.5;
- ((TVirtualPad*)l->At(0))->cd();
- selStart=0; for(n=0; n<6; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 1, 1, n, selection)) break;
- ((TVirtualPad*)l->At(1))->cd();
- selStart=6; for(n=0; n<6; n++) selection[n]=selStart+n;
- if(!GetGraphArray(xy, kMCtrack, 1, 1, n, selection)) break;
- return kTRUE;
- case 28: //kMCtrack [z]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-1.; xy[1]=-1500.; xy[2]=1.; xy[3] =6000.;
- ((TVirtualPad*)l->At(0))->cd();
- if(!GetGraphArray(xy, kMCtrack, 2)) break;
- xy[0] = -1.; xy[1] = -1.5; xy[2] = 1.; xy[3] = 5.;
- ((TVirtualPad*)l->At(1))->cd();
- if(!GetGraphArray(xy, kMCtrack, 3)) break;
- return kTRUE;
- case 29: //kMCtrack [phi/snp]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-.2; xy[1]=-0.5; xy[2]=.2; xy[3] =10.;
- ((TVirtualPad*)l->At(0))->cd();
- if(!GetGraphArray(xy, kMCtrack, 4)) break;
- xy[0] = -.2; xy[1] = -1.5; xy[2] = .2; xy[3] = 5.;
- ((TVirtualPad*)l->At(1))->cd();
- if(!GetGraphArray(xy, kMCtrack, 5)) break;
- return kTRUE;
- case 30: //kMCtrack [theta/tgl]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- xy[0]=-1.; xy[1]=-0.5; xy[2]=1.; xy[3] =5.;
- ((TVirtualPad*)l->At(0))->cd();
- if(!GetGraphArray(xy, kMCtrack, 6)) break;
- xy[0] = -.2; xy[1] = -0.5; xy[2] = .2; xy[3] = 2.5;
- ((TVirtualPad*)l->At(1))->cd();
- if(!GetGraphArray(xy, kMCtrack, 7)) break;
- return kTRUE;
- case 31: //kMCtrack [pt]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- pad = (TVirtualPad*)l->At(0); pad->cd();
- pad->SetMargin(0.125, 0.015, 0.1, 0.015);
- // pi selection
- n=0;
- for(Int_t il(ly0); il<AliTRDgeometry::kNlayer; il+=dly){
- selection[n++] = il*11 + 2; // pi-
- selection[n++] = il*11 + 8; // pi+
- }
- xy[0] = 0.2; xy[1] = -.7; xy[2] = 7.; xy[3] = 4.;
- //xy[0] = 0.2; xy[1] = -1.; xy[2] = 7.; xy[3] = 10.; // SA
- if(!GetGraphArray(xy, kMCtrack, 8, kTRUE, n, selection, "#pi#pm")) break;
- pad->Modified(); pad->Update(); pad->SetLogx();
- pad = (TVirtualPad*)l->At(1); pad->cd();
- pad->SetMargin(0.125, 0.015, 0.1, 0.015);
- // mu selection
- n=0;
- for(Int_t il(ly0); il<AliTRDgeometry::kNlayer; il+=dly){
- selection[n++] = il*11 + 3; // mu-
- selection[n++] = il*11 + 7; // mu+
- }
- if(!GetGraphArray(xy, kMCtrack, 8, kTRUE, n, selection, "#mu#pm")) break;
- pad->Modified(); pad->Update(); pad->SetLogx();
- return kTRUE;
- case 32: //kMCtrack [pt]
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- pad = (TVirtualPad*)l->At(0); pad->cd();
- pad->SetMargin(0.125, 0.015, 0.1, 0.015);
- // p selection
- n=0;
- for(Int_t il(ly0); il<AliTRDgeometry::kNlayer; il+=dly){
- selection[n++] = il*11 + 0; // p bar
- selection[n++] = il*11 + 10; // p
- }
- xy[0] = 0.2; xy[1] = -.7; xy[2] = 7.; xy[3] = 8.;
- //xy[0] = 0.2; xy[1] = -1.; xy[2] = 7.; xy[3] = 10.; // SA
- if(!GetGraphArray(xy, kMCtrack, 8, kTRUE, n, selection, "p&p bar")) break;
- pad->Modified(); pad->Update(); pad->SetLogx();
- pad = (TVirtualPad*)l->At(1); pad->cd();
- pad->SetMargin(0.125, 0.015, 0.1, 0.015);
- // e selection
- n=0;
- for(Int_t il(ly0); il<AliTRDgeometry::kNlayer; il+=dly){
- selection[n++] = il*11 + 4; // e-
- selection[n++] = il*11 + 6; // e+
- }
- xy[0] = 0.2; xy[1] = -1.5; xy[2] = 7.; xy[3] = 12.;
- //xy[0] = 0.2; xy[1] = -1.5; xy[2] = 7.; xy[3] = 14.; // SA
- if(!GetGraphArray(xy, kMCtrack, 8, kTRUE, n, selection, "e#pm")) break;
- pad->Modified(); pad->Update(); pad->SetLogx();
- return kTRUE;
- case 33: //kMCtrack [1/pt] pulls
- xy[0] = 0.; xy[1] = -1.; xy[2] = 2.; xy[3] = 3.5;
- //xy[0] = 0.; xy[1] = -1.; xy[2] = 2.; xy[3] = 4.5; // SA
- gPad->Divide(2, 1, 1.e-5, 1.e-5); l=gPad->GetListOfPrimitives();
- pad = (TVirtualPad*)l->At(0); pad->cd();
- pad->SetMargin(0.125, 0.015, 0.1, 0.015);
- // pi selection
- n=0;
- for(Int_t il(ly0); il<AliTRDgeometry::kNlayer; il+=dly){
- selection[n++] = il*11 + 2; // pi-
- selection[n++] = il*11 + 8; // pi+
- }
- if(!GetGraphArray(xy, kMCtrack, 9, kTRUE, n, selection, "#pi#pm")) break;
- pad = (TVirtualPad*)l->At(1); pad->cd();
- pad->SetMargin(0.125, 0.015, 0.1, 0.015);
- // mu selection
- n=0;
- for(Int_t il(ly0); il<AliTRDgeometry::kNlayer; il+=dly){
- selection[n++] = il*11 + 3; // mu-
- selection[n++] = il*11 + 7; // mu+