Use CA10RAM to get 10%* Discount.
Order Now
I just need help fixing the code under the example histogram so it can run
INSTRUCTIONS TO CANDIDATESANSWER ALL QUESTIONS
I just need help fixing the code under the example histogram so it can run
// $Id$
//
// File: DEventProcessor_search.cc
// Created: Mon Apr 3 11:38:03 EDT 2006
// Creator: davidl (on Darwin swire-b241.jlab.org 8.4.0 powerpc)
//
#include <map>
using namespace std;
#include "DEventProcessor_search.h"
// Routine used to create our DEventProcessor
extern "C"{
void InitPlugin(JApplication *app){
InitJANAPlugin(app);
app->AddProcessor(new DEventProcessor_search());
}
} // "C"
//------------------
// init
//------------------
jerror_t DEventProcessor_search::init(void)
{
return NOERROR;
}
//------------------
// brun
//------------------
jerror_t DEventProcessor_search::brun(JEventLoop *eventLoop, int32_t runnumber)
{
// Called for every run
// Get the calibration for the run and extract the beam spot position
DGeometry* dgeom = NULL;
DApplication* dapp = dynamic_cast<DApplication*> (eventLoop->GetJApplication());
if(dapp)
dgeom = dapp->GetDGeometry(runnumber);
if(dgeom){
dgeom->GetTargetZ(m_targetZ);
dgeom->GetFCALPosition( m_fcalX, m_fcalY, m_fcalZ );
dgeom->GetCCALPosition( m_ccalX, m_ccalY, m_ccalZ );
}else{
cerr << "No geometry accessible." << endl;
return RESOURCE_UNAVAILABLE;
}
//-----
jana::JCalibration *jcalib = japp->GetJCalibration(runnumber);
std::map<string, float> beam_spot;
jcalib->Get("PHOTON_BEAM/beam_spot", beam_spot);
m_beamSpotX = beam_spot.at("x");
m_beamSpotY = beam_spot.at("y");
cout << "Target Z: " << m_targetZ << endl;
cout << "FCal Z: " << m_fcalZ << endl;
cout << "CCal Z: " << m_ccalZ << endl;
return NOERROR;
}
//------------------
// evnt
//------------------
jerror_t DEventProcessor_search::evnt(JEventLoop *loop, uint64_t eventnumber){
//Called for every event
//Use vectors of data type to get them for each event
vector<const DFCALShower*> fcal; //Vector of DFCalShower pointers
TVector3 fP;
//Event loop will get objects that match the type of the vector
loop->Get(fcal);
//Example of a histogram
Fill1DWeightedHistogram("pi0_search", "example", "nFCalShowers", fcal.size(), 1, "Number of FCAL Showers in an event; Number of FCal Showers; Count", 10, 0, 10);
for(int i=0;i<fcal.size(); i++){
fcal[i]->getEnergy();
fcal[i]->getPosition().X();
fcal[i]->getPosition().Y();
fcal[i]->getPosition().Z();
DVector3 v1 = fcal[i]->getPosition();
v1.SetX(v1.X()- m_beamSpotX);
v1.SetY(v1.Y() - m_beamSpotY);
v1.SetZ(v1.Z() - m_targetZ);
v1.Phi();
v1.Theta();
double p1 = e1;
double px = p1 * sin(v1.Theta()) * cos(v1.Phi());
double py = p1 * sin(v1.Theta()) * cos(v1.Phi());
double pz = p1 * sin(v1.Theta());
v2.SetX(v2.X()-m_beamSpotX);
v2.SetY(v2.Y()-m_beamSpotY):
v2.SetZ(v2.Z()-m_targetZ);
v2.Phi();
v2.Theta();
double p2 = e2;
double px2 = p2 * sin(v2.Theta()) * cos(v2.Phi());
double py2 = p2 * sin(v2.Theta()) * cos(v2.Phi());
double pz2 = p2 * sin(v2.Theta());
v3.SetX(v3.X()-m_beamSpotX);
v3.SetY(v3.Y()-m_beamSpotY):
v3.SetZ(v3.Z()-m_targetZ);
v3.Phi();
v3.Theta();
double p3 = e3;
double px3 = p3 * sin(v3.Theta()) * cos(v3.Phi());
double py3 = p3 * sin(v3.Theta()) * cos(v3.Phi());
double pz3 = p3 * sin(v3.Theta());
DVector3 v4 = fcal[i]->getEnergy();
v4.Sete();
for(int i=0;i<fcal.size();i++){
TLorentzVector sh(px,py,pz,E1);
}
if(sh.E()>0.5){
continue;
}
else{
for(int j=i+2;j<fcal.size(); j++);{
TLorentzVector sh1(px2,py2,pz2,E2);
if(sh1.E()>0.5){
continue;
}
else{
for(int k=i+2;k<fcal.size(); k++);{
TLorentzVector sh2(px3,py3,pz3,E3);
if(sh2.E()>0.5){
continue;
}
else{
TLorentzVector w = sh + sh1 + sh2;
Double_t pho = w.M();
Fill1DWeightedHistogram("pi0_search2","shower2","nFCalShowers2",fcal.size(), 1, "Number of FCal Showers in an event; Number of FCal Showers; Count", 10, 0, 10);
}
}
}
}
}
}
}
return NOERROR;
}
}
//------------------
// erun
//------------------
jerror_t DEventProcessor_search::erun(void)
{
// Any final calculations on histogram (like dividing them)
// should be done here. This may get called more than once.
return NOERROR;
}
//------------------
// fini
//------------------
jerror_t DEventProcessor_search::fini(void)
{
return NOERROR;
}
Attachments:
Related Questions
. Introgramming & Unix Fall 2018, CRN 44882, Oakland University Homework Assignment 6 - Using Arrays and Functions in C
DescriptionIn this final assignment, the students will demonstrate their ability to apply two ma
. The standard path finding involves finding the (shortest) path from an origin to a destination, typically on a map. This is an
Path finding involves finding a path from A to B. Typically we want the path to have certain properties,such as being the shortest or to avoid going t
. Develop a program to emulate a purchase transaction at a retail store. This program will have two classes, a LineItem class and a Transaction class. The LineItem class will represent an individual
Develop a program to emulate a purchase transaction at a retail store. Thisprogram will have two classes, a LineItem class and a Transaction class. Th
. SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of Sea Ports. Here are the classes and their instance variables we wish to define:
1 Project 1 Introduction - the SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of
. Project 2 Introduction - the SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of Sea Ports. Here are the classes and their instance variables we wish to define:
1 Project 2 Introduction - the SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of
