Skip to content
Snippets Groups Projects
Commit 5bd13bfe authored by clantz's avatar clantz
Browse files

Changed recorded coordinates to x and z to reflect beamline coordinates

parent e23b7495
No related branches found
No related tags found
No related merge requests found
...@@ -21,9 +21,9 @@ int main(int argc, char *argv[]){ ...@@ -21,9 +21,9 @@ int main(int argc, char *argv[]){
if(name.find(".") != string::npos) name.erase( name.find_last_of(".") ); if(name.find(".") != string::npos) name.erase( name.find_last_of(".") );
TTree *t = (TTree*)f->Get("lightGuide"); TTree *t = (TTree*)f->Get("lightGuide");
double x,y; double x,z;
t->SetBranchAddress("X",&x); t->SetBranchAddress("X",&x);
t->SetBranchAddress("Y",&y); t->SetBranchAddress("Z",&z);
TCanvas *c = new TCanvas("Data","Data",1280,720); TCanvas *c = new TCanvas("Data","Data",1280,720);
TH2D *h = new TH2D("eff","Relative Efficiency",40,-82,82,40,-42.875,42.875); TH2D *h = new TH2D("eff","Relative Efficiency",40,-82,82,40,-42.875,42.875);
...@@ -31,7 +31,7 @@ int main(int argc, char *argv[]){ ...@@ -31,7 +31,7 @@ int main(int argc, char *argv[]){
int nEvents = t->GetEntries(); int nEvents = t->GetEntries();
for(int ev = 0; ev < nEvents; ev++){ for(int ev = 0; ev < nEvents; ev++){
t->GetEntry(ev); t->GetEntry(ev);
h->Fill(y,x); h->Fill(z,x);
} }
gStyle->SetOptStat(0); gStyle->SetOptStat(0);
......
...@@ -75,9 +75,9 @@ void EventAction::EndOfEventAction(const G4Event* event){ ...@@ -75,9 +75,9 @@ void EventAction::EndOfEventAction(const G4Event* event){
// fill ntuple // // fill ntuple //
G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
analysisManager->FillNtupleDColumn(0, aHit->getPos().x() ); analysisManager->FillNtupleDColumn(0, aHit->getPos().x() );
analysisManager->FillNtupleDColumn(1, aHit->getPos().y() ); analysisManager->FillNtupleDColumn(1, aHit->getPos().z() );
analysisManager->FillNtupleDColumn(2, aHit->getHit().x() ); analysisManager->FillNtupleDColumn(2, aHit->getHit().x() );
analysisManager->FillNtupleDColumn(3, aHit->getHit().y() ); analysisManager->FillNtupleDColumn(3, aHit->getHit().z() );
analysisManager->FillNtupleDColumn(4, aHit->getTime() ); analysisManager->FillNtupleDColumn(4, aHit->getTime() );
analysisManager->AddNtupleRow(); analysisManager->AddNtupleRow();
} }
......
...@@ -63,9 +63,9 @@ RunAction::RunAction(G4String fileName) ...@@ -63,9 +63,9 @@ RunAction::RunAction(G4String fileName)
// Creating ntuple // Creating ntuple
analysisManager->CreateNtuple("lightGuide", "pos and momentum"); analysisManager->CreateNtuple("lightGuide", "pos and momentum");
analysisManager->CreateNtupleDColumn("X"); analysisManager->CreateNtupleDColumn("X");
analysisManager->CreateNtupleDColumn("Y"); analysisManager->CreateNtupleDColumn("Z");
analysisManager->CreateNtupleDColumn("hitX"); analysisManager->CreateNtupleDColumn("hitX");
analysisManager->CreateNtupleDColumn("hitY"); analysisManager->CreateNtupleDColumn("hitZ");
analysisManager->CreateNtupleDColumn("time"); analysisManager->CreateNtupleDColumn("time");
analysisManager->FinishNtuple(); analysisManager->FinishNtuple();
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment