Subversion Repositories f9daq

Rev

Rev 273 | Rev 319 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 273 Rev 277
Line 75... Line 75...
75 
   if (svar.Contains("GetCharge"  )) fHtype[id]=3;
 75 
   if (svar.Contains("GetCharge"  )) fHtype[id]=3;
76 
   if (svar.Contains("GetPid"     )) fHtype[id]=4;
 76 
   if (svar.Contains("GetPid"     )) fHtype[id]=4;
77 
   if (svar.Contains("GetXMomentum")) fHtype[id]=5;
 77 
   if (svar.Contains("GetXMomentum")) fHtype[id]=5;
78 
   if (svar.Contains("GetYMomentum")) fHtype[id]=6;
 78 
   if (svar.Contains("GetYMomentum")) fHtype[id]=6;
79 
   if (svar.Contains("GetZMomentum")) fHtype[id]=7;
 79 
   if (svar.Contains("GetZMomentum")) fHtype[id]=7;
80 
   if (svar.Contains("GetTranverseMomentum")) fHtype[id]=8;
 80 
   if (svar.Contains("GetTransverseMomentum")) fHtype[id]=8;
81 
   if (svar.Contains("GetTheta"))             fHtype[id]=9;
 81 
   if (svar.Contains("GetTheta"))             fHtype[id]=9;
82 
   if (svar.Contains("GetCosTheta"))          fHtype[id]=10;
 82 
   if (svar.Contains("GetCosTheta"))          fHtype[id]=10;
83 
 
 83 
 
84 
   fH[id]= new TH1F(TString::Format("h%d",id), TString::Format("%s;%s;N",name,axis[fHtype[id]].Data()), nbins, min, max);
 84 
   fH[id]= new TH1F(TString::Format("h%d",id), TString::Format("%s;%s;N",name,axis[fHtype[id]].Data()), nbins, min, max);
85 
   
 85 
   
Line 104... Line 104...
104 
      case 4 : val  = p->pid(); break;
 104 
      case 4 : val  = p->pid(); break;
105 
      case 5 : val  = p->px(); break;
 105 
      case 5 : val  = p->px(); break;
106 
      case 6 : val  = p->py(); break;
 106 
      case 6 : val  = p->py(); break;
107 
      case 7 : val  = p->pz(); break;
 107 
      case 7 : val  = p->pz(); break;
108 
      case 8 : val  = p->GetTransverseMomentum(); break;
 108 
      case 8 : val  = p->GetTransverseMomentum(); break;
109 
      case 9 : val  = (p->GetMomentum()!=0) ? p->pz()/p->GetMomentum() : 0; val = 180*acos(val)/TMath::Pi(); break;
 109 
      case 9 : val  = (p->GetMomentum()!=0) ? p->pz()/p->GetMomentum() : 0; val = 180.0*TMath::ACos(val)/TMath::Pi(); break;
110 
      case 10: val  = (p->GetMomentum()!=0) ? p->pz()/p->GetMomentum() : 0; break;
 110 
      case 10: val  = (p->GetMomentum()!=0) ? p->pz()/p->GetMomentum() : 0; break;
111 
      default: val  = 0 ; break;
 111 
      default: val  = 0 ; break;
112 
   }
 112 
   }
113 
   fH[hid]->Fill(val);
 113 
   fH[hid]->Fill(val);
114 
}  
 114 
}