now everything works fine, apparently

Author: maxgalli

src/SCECalorimeterSD.cc

@@ -85,6 +85,7 @@ G4bool SCECalorimeterSD::ProcessHits(G4Step* step,
   hit->Add(edep, stepLength);
   hitTotal->Add(edep, stepLength);
   hit->SetPos(pos);
+  hitTotal->SetPos(pos);
 
   return true;
 }

src/SCEEventAction.cc

@@ -131,27 +131,21 @@ void SCEEventAction::EndOfEventAction(const G4Event* event)
                         analysisManager->FillH1(id + 3, gapHit->GetTrackLength());
                 }
 
-                G4cout << "the size of step vector is: " << gapHit->GetPos().size() << G4endl;
-
-                // fill ntuple
+                // fill ntuple and pos histo
                 analysisManager->FillNtupleDColumn(i, id_tup, absoHit->GetEdep());
                 analysisManager->FillNtupleDColumn(i, id_tup + 1, gapHit->GetEdep());
                 analysisManager->FillNtupleDColumn(i, id_tup + 2, absoHit->GetTrackLength());
                 analysisManager->FillNtupleDColumn(i, id_tup + 3, gapHit->GetTrackLength());
-                if (gapHit->GetPos().size() != 0) {
-                        for (auto j=gapHit->GetPos().begin(); j!=gapHit->GetPos().end(); j++) {
-                                G4cout << "Inside j for" << G4endl;
-                                analysisManager->FillH1(id + 4, std::sqrt(std::pow(j->x(),2)+std::pow(j->y(),2)));
-                                analysisManager->FillNtupleDColumn(i, id_tup + 4, std::sqrt(std::pow(j->x(),2)+std::pow(j->y(),2)));
+
+                if (gapHit->GetPos().empty() == false) {
+                        for (int j=0; j<(int)gapHit->GetPos().size(); ++j) {
+                                analysisManager->FillH1(id + 4, std::sqrt(std::pow(gapHit->GetPos()[j].x(),2)+std::pow(gapHit->GetPos()[j].y(),2)));
+                                analysisManager->FillH2(i, gapHit->GetPos()[j].x(), gapHit->GetPos()[j].y());
+                                analysisManager->FillNtupleDColumn(i, id_tup + 4, std::sqrt(std::pow(gapHit->GetPos()[j].x(),2)+std::pow(gapHit->GetPos()[j].y(),2)));
                         }
                         id += 5;
                         analysisManager->AddNtupleRow(i);
                 }
-
-                else {
-                        id += 4;
-                        analysisManager->AddNtupleRow(i);
-                }
         }
 }
 

src/SCERunAction.cc

@@ -41,7 +41,9 @@ SCERunAction::SCERunAction(int nenergy, int nlayers, G4String fabsmaterial, doub
                 analysisManager->CreateH1("Egap" + std::to_string(i),"Edep in gap", 100, 0., nEnergy*GeV);
                 analysisManager->CreateH1("Labs" + std::to_string(i),"trackL in absorber", 100, 0., nEnergy*m);
                 analysisManager->CreateH1("Lgap" + std::to_string(i),"trackL in gap", 100, 0., nEnergy*m);
-                analysisManager->CreateH1("Pgap" + std::to_string(i),"Transverse position in gap", 60, 0., 4*m);
+                analysisManager->CreateH1("Pgap" + std::to_string(i),"Transverse position in gap", 100, 0., 4*m);
+                analysisManager->CreateH2("Pgap2D" + std::to_string(i),"X and Y position", 100, -2.5*m, 2.5*m, 100, -2.5*m, 2.5*m);
+
 
                 // Creating ntuple
                 //
@@ -58,7 +60,8 @@ SCERunAction::SCERunAction(int nenergy, int nlayers, G4String fabsmaterial, doub
         analysisManager->CreateH1("Egap_tot","Edep in gap", 100, 0., nEnergy*GeV);
         analysisManager->CreateH1("Labs_tot","trackL in absorber", 100, 0., nEnergy*m);
         analysisManager->CreateH1("Lgap_tot","trackL in gap", 100, 0., nEnergy*m);
-                analysisManager->CreateH1("Pgap_tot","Transverse position in gap", 100, 0., 4.*m);
+        analysisManager->CreateH1("Pgap_tot","Transverse position in gap", 100, 0., 4.*m);
+        analysisManager->CreateH2("Pgap2D_tot","X and Y position", 100, -2.5*m, 2.5*m, 100, -2.5*m, 2.5*m);
 
         // Creating ntuple for tot
         //