Histo2DActions/ZenVsAzAngleHistoAction/ZenVsAzAngleHistoAction.cpp

/*
 * ZenVsAzAngleHistoAction.cpp
 *
 *  Created on: 28-marth-2010
 *      Author: Vitaly Malakhov
 */

/*! @file ZenVsAzAngleHistoAction.cpp The ZenVsAzAbnleHistoAction class implementation file. */

#include "ZenVsAzAngleHistoAction.h"
#include "TMath.h"

ZenVsAzAngleHistoAction::ZenVsAzAngleHistoAction(const char *actionName, TString outFileBase, TString mode,
    bool outRoot, bool outText, TString title) :
    Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText){

}

void ZenVsAzAngleHistoAction::OnGood(PamLevel2 *event) {
    PamTrack *track;
    track = event->GetTrack(0);
    Float_t Aaxv = TMath::Abs(track->GetTrkTrack()->axv[0])*TMath::DegToRad();
    Float_t Aayv = TMath::Abs(track->GetTrkTrack()->ayv[0])*TMath::DegToRad();
    Float_t PamZenith = TMath::RadToDeg()*asin(sqrt(pow(sin(Aayv), 2) + pow(sin(Aaxv), 2)));

    Float_t axv = -track->GetTrkTrack()->axv[0] * TMath::DegToRad();
    Float_t ayv = -track->GetTrkTrack()->ayv[0] * TMath::DegToRad();
    Float_t angle = atan(sin(TMath::Abs(ayv))/sin(TMath::Abs(axv))) * TMath::RadToDeg();

    Float_t PamAzim =  360. - angle;
    if(axv>=0 && ayv >=0) PamAzim = angle;
    if(axv<0 && ayv >0) PamAzim = 180. - angle;
    if(axv<0 && ayv <0) PamAzim = 180. + angle;

    PamAzim = PamAzim * TMath::DegToRad();
    PamZenith = (180 - PamZenith) * TMath::DegToRad();

    Fill(TMath::RadToDeg()*PamAzim, TMath::RadToDeg()*PamZenith);
}
1	/*
2	* ZenVsAzAngleHistoAction.cpp
3	*
4	* Created on: 28-marth-2010
5	* Author: Vitaly Malakhov
6	*/
7
8	/! @file ZenVsAzAngleHistoAction.cpp The ZenVsAzAbnleHistoAction class implementation file. /
9
10	#include "ZenVsAzAngleHistoAction.h"
11	#include "TMath.h"
12
13	ZenVsAzAngleHistoAction::ZenVsAzAngleHistoAction(const char *actionName, TString outFileBase, TString mode,
14	bool outRoot, bool outText, TString title) :
15	Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText){
16
17	}
18
19	void ZenVsAzAngleHistoAction::OnGood(PamLevel2 *event) {
20	PamTrack *track;
21	track = event->GetTrack(0);
22	Float_t Aaxv = TMath::Abs(track->GetTrkTrack()->axv[0])*TMath::DegToRad();
23	Float_t Aayv = TMath::Abs(track->GetTrkTrack()->ayv[0])*TMath::DegToRad();
24	Float_t PamZenith = TMath::RadToDeg()*asin(sqrt(pow(sin(Aayv), 2) + pow(sin(Aaxv), 2)));
25
26	Float_t axv = -track->GetTrkTrack()->axv[0] * TMath::DegToRad();
27	Float_t ayv = -track->GetTrkTrack()->ayv[0] * TMath::DegToRad();
28	Float_t angle = atan(sin(TMath::Abs(ayv))/sin(TMath::Abs(axv))) * TMath::RadToDeg();
29
30	Float_t PamAzim = 360. - angle;
31	if(axv>=0 && ayv >=0) PamAzim = angle;
32	if(axv<0 && ayv >0) PamAzim = 180. - angle;
33	if(axv<0 && ayv <0) PamAzim = 180. + angle;
34
35	PamAzim = PamAzim * TMath::DegToRad();
36	PamZenith = (180 - PamZenith) * TMath::DegToRad();
37
38	Fill(TMath::RadToDeg()PamAzim, TMath::RadToDeg()PamZenith);
39	}