--- PamCut/CollectionActions/Histo2DActions/MassVsRigHistoAction/MassVsRigHistoAction.cpp 2010/04/05 11:07:16 1.1.1.1 +++ PamCut/CollectionActions/Histo2DActions/MassVsRigHistoAction/MassVsRigHistoAction.cpp 2010/04/07 10:31:43 1.3 @@ -11,12 +11,13 @@ MassVsRigHistoAction::MassVsRigHistoAction(const char *actionName, TString outFileBase, TString mode, bool outRoot, bool outText, TString title, Int_t Z, TString rigBinFile) : - Histo2DAction (actionName, title, outFileBase, mode, outRoot, outText), _Z(Z), _rigBinFile(RigBinFile){ + Histo2DAction (actionName, title, outFileBase, mode, outRoot, outText), _Z(Z), _rigBinFile(rigBinFile){ } void MassVsRigHistoAction::OnGood(PamLevel2 *event) { Float_t betaev = event->GetTrack(0)->GetToFTrack()->beta[12]; + if(betaev>1) betaev = 2-betaev; Float_t rigev = 1./event->GetTrack(0)->GetTrkTrack()->GetDeflection(); Float_t Mass = _Z*rigev*sqrt(pow(betaev,-2)-1); Fill(rigev, Mass);