| 1 | } else if(TString(GetType())=="MissingET") {
|
|---|
| 2 |
|
|---|
| 3 | // MissingET as invisible track (like a photon)
|
|---|
| 4 | TEveTrackPropagator *trkProp = new TEveTrackPropagator();
|
|---|
| 5 | trkProp->SetMagField(0., 0., -tk_Bz_);
|
|---|
| 6 | trkProp->SetMaxR(tkRadius_);
|
|---|
| 7 | trkProp->SetMaxZ(tkHalfLength_);
|
|---|
| 8 |
|
|---|
| 9 | TEveTrack *eveTrack;
|
|---|
| 10 |
|
|---|
| 11 | Int_t counter = 0;
|
|---|
| 12 |
|
|---|
| 13 | TIter itMet(branch_);
|
|---|
| 14 | MissingET *MET;
|
|---|
| 15 |
|
|---|
| 16 | // Missing Et
|
|---|
| 17 | while((MET = (MissingET*) itMet.Next())) {
|
|---|
| 18 | TParticle pb(13, 1, 0, 0, 0, 0,
|
|---|
| 19 | (tkRadius_ * MET->MET/maxPt_)*cos(MET->Phi),
|
|---|
| 20 | (tkRadius_ * MET->MET/maxPt_)*sin(MET->Phi),
|
|---|
| 21 | 0., MET->MET,
|
|---|
| 22 | // MET->P4().Pz(), MET->P4().E(),
|
|---|
| 23 | 0.0, 0.0, 0.0, 0.0);
|
|---|
| 24 | eveTrack = new TEveTrack(&pb, counter, trkProp);
|
|---|
| 25 | eveTrack->SetName(Form("MissingEt [%d]", counter++));
|
|---|
| 26 | eveTrack->SetStdTitle();
|
|---|
| 27 |
|
|---|
| 28 | eveTrack->SetRnrPoints(0);
|
|---|
| 29 | eveTrack->SetMarkerColor(kMagenta);
|
|---|
| 30 | eveTrack->SetMarkerStyle(4);
|
|---|
| 31 | eveTrack->SetMarkerSize(2.);
|
|---|
| 32 | eveTrack->SetLineWidth(2.);
|
|---|
| 33 |
|
|---|
| 34 | eveTrack->SetLineStyle(7);
|
|---|
| 35 | data_->AddElement(eveTrack);
|
|---|
| 36 | eveTrack->SetLineColor(GetColor());
|
|---|
| 37 | eveTrack->MakeTrack();
|
|---|
| 38 | }
|
|---|
| 39 |
|
|---|
| 40 | /* old version
|
|---|
| 41 | TIter itMet(branch_);
|
|---|
| 42 | MissingET *MET;
|
|---|
| 43 | TEveArrow *eveMet;
|
|---|
| 44 | // Missing Et
|
|---|
| 45 | while((MET = (MissingET*) itMet.Next())) {
|
|---|
| 46 | eveMet = new TEveArrow((tkRadius_ * MET->MET/maxPt_)*cos(MET->Phi), (tkRadius_ * MET->MET/maxPt_)*sin(MET->Phi), 0., 0., 0., 0.);
|
|---|
| 47 | eveMet->SetMainColor(GetColor());
|
|---|
| 48 | eveMet->SetTubeR(0.04);
|
|---|
| 49 | eveMet->SetConeR(0.08);
|
|---|
| 50 | eveMet->SetConeL(0.10);
|
|---|
| 51 | eveMet->SetPickable(kTRUE);
|
|---|
| 52 | eveMet->SetName("Missing Et");
|
|---|
| 53 | eveMet->SetTitle(Form("Missing Et (%.1f GeV)",MET->MET));
|
|---|
| 54 | data_->AddElement(eveMet);
|
|---|
| 55 | }
|
|---|
| 56 | */
|
|---|
| 57 | }
|
|---|
