/[PAMELA software]/PamCut/CollectionActions/Histo2DActions/MassVsRigHistoAction/MassVsRigHistoAction.cpp
ViewVC logotype

Diff of /PamCut/CollectionActions/Histo2DActions/MassVsRigHistoAction/MassVsRigHistoAction.cpp

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

revision 1.1 by pam-mep, Mon Apr 5 11:07:16 2010 UTC revision 1.3 by pam-mep, Wed Apr 7 10:31:43 2010 UTC
# Line 11  Line 11 
11    
12  MassVsRigHistoAction::MassVsRigHistoAction(const char *actionName, TString outFileBase, TString mode,  MassVsRigHistoAction::MassVsRigHistoAction(const char *actionName, TString outFileBase, TString mode,
13      bool outRoot, bool outText, TString title, Int_t Z, TString rigBinFile) :      bool outRoot, bool outText, TString title, Int_t Z, TString rigBinFile) :
14      Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText), _Z(Z), _rigBinFile(RigBinFile){      Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText), _Z(Z), _rigBinFile(rigBinFile){
15    
16  }  }
17    
18  void MassVsRigHistoAction::OnGood(PamLevel2 *event) {  void MassVsRigHistoAction::OnGood(PamLevel2 *event) {
19      Float_t betaev = event->GetTrack(0)->GetToFTrack()->beta[12];      Float_t betaev = event->GetTrack(0)->GetToFTrack()->beta[12];
20        if(betaev>1) betaev = 2-betaev;
21      Float_t rigev = 1./event->GetTrack(0)->GetTrkTrack()->GetDeflection();      Float_t rigev = 1./event->GetTrack(0)->GetTrkTrack()->GetDeflection();
22      Float_t Mass = _Z*rigev*sqrt(pow(betaev,-2)-1);      Float_t Mass = _Z*rigev*sqrt(pow(betaev,-2)-1);
23      Fill(rigev, Mass);      Fill(rigev, Mass);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23