#include <HitCalcAna.h>
Inheritance diagram for HitCalcAna:
Public Types | |
| typedef std::deque< Int_t > | DeqInt_t |
| typedef DeqInt_t::iterator | IterDeqInt_t |
| typedef DeqInt_t::size_type | SizeDeqInt_t |
| typedef std::deque< Float_t > | DeqFloat_t |
| typedef DeqFloat_t::iterator | IterDeqFloat_t |
| typedef std::vector< Int_t > | VecInt_t |
| typedef VecInt_t::iterator | IterVecInt_t |
| typedef std::vector< Float_t > | VecFloat_t |
| typedef VecFloat_t::iterator | IterVecFloat_t |
| typedef std::deque< std::deque< Int_t > > | DeqDeqInt_t |
| typedef DeqDeqInt_t::iterator | IterDeqDeqInt_t |
| typedef DeqDeqInt_t::size_type | SizeDeqDeqInt_t |
Public Member Functions | |
| HitCalcAna (HitCalc &hc) | |
| virtual | ~HitCalcAna () |
| void | Analyze (int evtn, RecRecordImp< RecCandHeader > *srobj) |
| Int_t | ComputeHits (RecRecordImp< RecCandHeader > *srobj, NtpSREvent *event) |
| void | Get3DHit (DeqFloat_t &x, DeqFloat_t &y, DeqFloat_t &z, DeqFloat_t &e) |
| void | Reset () |
Public Attributes | |
| DeqFloat_t | fX |
| DeqFloat_t | fY |
| DeqFloat_t | fZ |
| DeqFloat_t | fE |
| IterDeqFloat_t | fHitIndex |
Private Member Functions | |
| void | VarCalc () |
Private Attributes | |
| HitCalc & | fHitCalc |
|
|
Definition at line 44 of file HitCalcAna.h. |
|
|
Definition at line 35 of file HitCalcAna.h. Referenced by ComputeHits(). |
|
|
Definition at line 31 of file HitCalcAna.h. |
|
|
Definition at line 45 of file HitCalcAna.h. |
|
|
Definition at line 36 of file HitCalcAna.h. |
|
|
Definition at line 32 of file HitCalcAna.h. |
|
|
Definition at line 42 of file HitCalcAna.h. |
|
|
Definition at line 39 of file HitCalcAna.h. |
|
|
Definition at line 46 of file HitCalcAna.h. |
|
|
Definition at line 33 of file HitCalcAna.h. |
|
|
Definition at line 41 of file HitCalcAna.h. Referenced by VarCalc(). |
|
|
Definition at line 38 of file HitCalcAna.h. |
|
|
Definition at line 37 of file HitCalcAna.cxx. 00037 : 00038 fHitCalc(hc) 00039 { 00040 00041 00042 }
|
|
|
Definition at line 44 of file HitCalcAna.cxx. 00045 {
00046
00047
00048 }
|
|
||||||||||||
|
Implements NueAnaBase. Definition at line 50 of file HitCalcAna.cxx. References ComputeHits(), SntpHelpers::GetEvent(), RecRecordImp< T >::GetHeader(), RecPhysicsHeader::GetSnarl(), MSG, and VarCalc(). Referenced by NueDisplayModule::Analyze(), and NueRecordAna::FillReco(). 00052 {
00053
00054 if(srobj==0){
00055 return;
00056 }
00057 if(((dynamic_cast<NtpStRecord *>(srobj))==0)&&
00058 ((dynamic_cast<NtpSRRecord *>(srobj))==0)){
00059 return;
00060 }
00061
00062 MSG("HitCalcAna",Msg::kInfo)<<"In HitCalcAna::Analyze"<<endl;
00063 MSG("HitCalcAna",Msg::kInfo)<<"On Snarl "<<srobj->GetHeader().GetSnarl()
00064 <<" event "<<evtn<<endl;
00065
00066 NtpSREvent *event = SntpHelpers::GetEvent(evtn,srobj);
00067 if(!event){
00068 MSG("HitCalcAna",Msg::kError)<<"Couldn't get event "<<evtn
00069 <<" from Snarl "
00070 <<srobj->GetHeader().GetSnarl()<<endl;
00071 return;
00072 }
00073
00074 Int_t hitTotal;
00075
00076 //Match strips into 3D Hits.
00077 hitTotal=ComputeHits(srobj,event);
00078
00079 if (hitTotal > 3){
00080 //Calculate simple hit variables and fill branch.
00081 VarCalc();
00082
00083 }
00084
00085 }
|
|
||||||||||||
|
Definition at line 110 of file HitCalcAna.cxx. References DeqFloat_t, NtpSREvent::end, fE, NtpVtxFinder::FindVertex(), fX, fY, fZ, VldContext::GetDetector(), RecRecordImp< T >::GetHeader(), RecPhysicsHeader::GetSnarl(), SntpHelpers::GetStrip(), SntpHelpers::GetStripIndex(), SntpHelpers::GetTrack(), RecHeader::GetVldContext(), ReleaseType::IsCedar(), MSG, NtpSRPlane::n, NtpSREvent::nstrip, NtpSREvent::ntrack, NtpSRStrip::plane, NtpSRTrack::plane, pow(), NtpVtxFinder::SetTargetEvent(), NtpSRStrip::tpos, NtpSREvent::trk, NtpSREvent::vtx, NtpVtxFinder::VtxX(), NtpVtxFinder::VtxY(), NtpVtxFinder::VtxZ(), NtpSRVertex::x, NtpSRVertex::y, NtpSRStrip::z, and NtpSRVertex::z. Referenced by Analyze(). 00112 {
00113
00114 if(srobj==0){
00115 return 1;
00116 }
00117 if(((dynamic_cast<NtpStRecord *>(srobj))==0)&&
00118 ((dynamic_cast<NtpSRRecord *>(srobj))==0)){
00119 return 1;
00120 }
00121
00122 // Clear 3D Hit Coordinate vectors
00123 fX.clear();
00124 fY.clear();
00125 fZ.clear();
00126 fE.clear();
00127
00128 //Determine detector type from Validity Context
00129 DetectorType::Detector_t fDetectorType =
00130 srobj->GetHeader().GetVldContext().GetDetector();
00131
00132 //Scale sigcor to pe.
00133 // const Float_t kNFScale = 1.0/100.0; // 92.3/65.2;
00134
00135 Float_t totalVisEnerg;
00136
00137 // Auxiliary hit variables
00138 Float_t hitX;
00139 Float_t hitY;
00140 Float_t hitZ;
00141 Float_t hitE;
00142
00143 //Strip vectors
00144 DeqFloat_t stripUV;
00145 DeqFloat_t stripZ;
00146 DeqFloat_t stripE;
00147
00148 DeqFloat_t stripUVTemp;
00149 DeqFloat_t stripZTemp;
00150 DeqFloat_t stripETemp;
00151
00152 Int_t nPlaneSM = 0; // account for U and V even-odd switch
00153 // from FD SM1 to SM2
00154
00155 Int_t nMax = 0; //Total number of Hits found
00156 Int_t nHit = 0; //Number of hits counter
00157
00158 Float_t vertexX = event->vtx.x; //Munits::meters
00159 Float_t vertexY = event->vtx.y; //Munits::meters
00160 Float_t vertexZ = event->vtx.z; //Munits::meters
00161
00162 if(ReleaseType::IsCedar(release)){
00163 NtpSRTrack *primaryTrack = 0;
00164 NtpSRTrack *ntpTrack = 0;
00165 Int_t index = 0;
00166
00167 for(Int_t i = 0; i < event->ntrack; ++i){
00168 index = event->trk[i];
00169 ntpTrack = SntpHelpers::GetTrack(index, srobj);
00170 if(!primaryTrack) primaryTrack = ntpTrack;
00171
00172 if(ntpTrack->plane.n > primaryTrack->plane.n)
00173 primaryTrack = ntpTrack;
00174 }
00175
00176 NtpStRecord* st = dynamic_cast<NtpStRecord *>(srobj);
00177 NtpVtxFinder vtxf;
00178 vtxf.SetTargetEvent(st, event, primaryTrack);
00179 if(vtxf.FindVertex() > 0){
00180 vertexX = vtxf.VtxX();
00181 vertexY = vtxf.VtxY();
00182 vertexZ = vtxf.VtxZ();
00183 }
00184 }
00185
00186 Int_t event_end = event -> plane.end;
00187 Int_t event_start = event -> plane.beg;
00188
00189 if(fDetectorType == DetectorType::kNear){
00190 if (event_end >= 121) {
00191 MSG("HitCalcAna",Msg::kInfo)<< "ND event ends outside"
00192 << " calorimeter region in plane "
00193 << event_end
00194 <<". Truncated hit matching."
00195 << endl;
00196 event_end = 121;
00197 }
00198 if (event_start >= 121) {
00199 MSG("HitCalcAna",Msg::kError)<< "ND event starts outside"
00200 << " calorimeter region in "
00201 << event_start <<" Stopping."
00202 << endl;
00203 return 1;
00204 }
00205 }
00206
00207 totalVisEnerg = 0.0;
00208 nHit=0;
00209
00210 Float_t stpEDiff = 0.0;
00211
00212 //Loop over all planes in the event
00213 for(Int_t nPlane = event_start;
00214 nPlane < event_end; nPlane++){
00215
00216 //Reset the strip vectors for a new plane
00217 stripUVTemp.clear();
00218 stripZTemp.clear();
00219 stripETemp.clear();
00220
00221 stripUV.clear();
00222 stripZ.clear();
00223 stripE.clear();
00224
00225 // Loop over all strips in the event
00226 for(Int_t iStp = 0; iStp < event->nstrip; iStp++){
00227 Int_t index = SntpHelpers::GetStripIndex(iStp,event);
00228 NtpSRStrip *strip = SntpHelpers::GetStrip(index,srobj);
00229
00230 if(!strip){
00231 MSG("HitCalcAna",Msg::kDebug)<< "Couldn't get strip "<< index
00232 << " in snarl "
00233 << srobj->GetHeader().GetSnarl()
00234 << " next"<<endl;
00235 continue;
00236 }
00237 if(!evtstp0mip){
00238 MSG("HitCalcAna",Msg::kError)<<"No mip strip information"<<endl;
00239 continue;
00240 }
00241
00242
00243 float charge = evtstp0mip[index] + evtstp1mip[index];;
00244
00245 totalVisEnerg += charge;
00246
00247 //Eliminate crosstalk strips
00248 if(TMath::Abs(strip->z-vertexZ) < 4. && (charge)>200.0/sigcormeu){
00249
00250 if(strip->plane == nPlane){
00251 stripUVTemp.push_back(strip->tpos);
00252 stripZTemp.push_back(strip ->z);
00253 stripETemp.push_back(charge);
00254 }
00255
00256 if(strip->plane == (nPlane+1)){
00257 stripUV.push_back(strip->tpos);
00258 stripZ.push_back(strip ->z);
00259 stripE.push_back(charge);
00260 }
00261 }
00262 } // iStp
00263
00264
00265 // Loop over all strips and match hits with strips in
00266 // consecutive planes
00267 for(UInt_t nStrip = 0; nStrip < stripZ.size(); nStrip++) {
00268 for(UInt_t nStripTemp = 0; nStripTemp < stripZTemp.size()
00269 ; nStripTemp++){
00270
00271 MSG("HitCalcAna",Msg::kDebug)<< "StripUV(" << nStrip
00272 << ") = "
00273 << stripUV.at(nStrip) << " "
00274 << "StripUVTemp(" << nStripTemp
00275 << ") = "
00276 << stripUVTemp.at(nStripTemp)
00277 << endl;
00278 MSG("HitCalcAna",Msg::kDebug)<< "StripZ(" << nStrip
00279 << ") = "
00280 << stripZ.at(nStrip) << " "
00281 << "StripZTemp(" << nStripTemp
00282 << ") = "
00283 << stripZTemp.at(nStripTemp)
00284 << endl;
00285
00286 MSG("HitCalcAna",Msg::kDebug)<< "StripE(" << nStrip
00287 << ") = "
00288 << stripE.at(nStrip) << " "
00289 << "StripETemp(" << nStripTemp
00290 << ") = "
00291 << stripETemp.at(nStripTemp)
00292 << endl;
00293
00294 //Only match hits if there are two strips in
00295 //consecutive planes without unphysical values.
00296 if(stripUV.at(nStrip) != stripUVTemp.at(nStripTemp)
00297 && (stripUV.at(nStrip) > -99.
00298 && stripUVTemp.at(nStripTemp) > -99.)){
00299
00300 if(fDetectorType == DetectorType::kFar){ //far detector
00301 stpEDiff=10; //sigcormeu (1GeV~25.76 sigcormeu)
00302
00303 if(nPlane >= 0 && nPlane <= 249)
00304 nPlaneSM = nPlane+1; //u - even and v - odd
00305 if(nPlane >= 250 && nPlane <= 485)
00306 nPlaneSM = nPlane; //u - odd and v - even
00307 }
00308
00309 if(fDetectorType == DetectorType::kNear||
00310 fDetectorType==DetectorType::kCalDet){
00311 nPlaneSM = nPlane; // near detector
00312 stpEDiff=10; //sigcormeu
00313 }
00314 // Only match hits if adjacent strips have reasonably
00315 // similar energy deposition.
00316
00317 // cout << "StpEDiff " << stpEDiff << endl;
00318
00319
00320 if(TMath::Abs(stripE.at(nStrip)
00321 -stripETemp.at(nStripTemp)) < stpEDiff ){
00322
00323 MSG("HitCalcAna",Msg::kDebug)<< "Found Hit "
00324 << nHit << endl;
00325
00326 hitX = pow((-1.0),nPlaneSM)*(sqrt(2.0)/2.0)
00327 *(stripUV.at(nStrip)-stripUVTemp.at(nStripTemp));
00328 hitY = (sqrt(2.0)/2.0)
00329 *(stripUV.at(nStrip)+stripUVTemp.at(nStripTemp));
00330 hitZ = (stripZ.at(nStrip)
00331 +stripZTemp.at(nStripTemp))/2.0;
00332 hitE = (stripE.at(nStrip)
00333 +stripETemp.at(nStripTemp))/2.0;
00334
00335 fX.push_back(hitX-vertexX);
00336 fY.push_back(hitY-vertexY);
00337 fZ.push_back(hitZ-vertexZ);
00338 fE.push_back(hitE);
00339
00340 MSG("HitCalcAna",Msg::kDebug)<< "Vertex Pos(X,Y,Z)= ("
00341 << vertexX << ","
00342 << vertexY << ","
00343 << vertexZ << ")"
00344 << endl << endl;
00345
00346 MSG("HitCalcAna",Msg::kDebug)<< "fX["
00347 << nHit <<"]= "
00348 << fX[nHit]
00349 << endl << endl;
00350
00351 MSG("HitCalcAna",Msg::kDebug)<< "fY["
00352 << nHit <<"]= "
00353 << fY[nHit]
00354 << endl << endl;
00355
00356 MSG("HitCalcAna",Msg::kDebug)<< "fZ["
00357 << nHit <<"]= "
00358 << fZ[nHit]
00359 << endl << endl;
00360
00361 MSG("HitCalcAna",Msg::kDebug)<< "fE["
00362 << nHit <<"]= "
00363 << fE[nHit]
00364 << endl << endl;
00365
00366 nHit++;
00367
00368 } else continue;
00369 }
00370 } // nStripTemp
00371 } // nStrip
00372 } // nPlane
00373
00374 nMax = nHit;
00375
00376 MSG("HitCalcAna",Msg::kInfo)<< "Found "
00377 << nMax
00378 << " 3D Hits."
00379 << endl;
00380
00381 return nMax;
00382 }
|
|
||||||||||||||||||||
|
Definition at line 87 of file HitCalcAna.cxx. References fE, fX, fY, fZ, and MSG. Referenced by NueRecordAna::Analyze(), NueDisplayModule::Analyze(), and NueRecordAna::FillReco(). 00091 {
00092
00093 if(fX.size()==0 || fY.size()==0 || fZ.size()==0 || fE.size()==0){
00094 MSG("HitCalcAna",Msg::kWarning)<< "3D Hits not calculated for event,"
00095 << "stopping any further"
00096 << " Hit related processing."
00097 << endl;
00098 return;
00099 }
00100
00101 //Reference the 3D Hit vectors.
00102 x=fX;
00103 y=fY;
00104 z=fZ;
00105 e=fE;
00106
00107 }
|
|
|
Definition at line 385 of file HitCalcAna.cxx. 00386 {
00387
00388
00389 }
|
|
|
Definition at line 391 of file HitCalcAna.cxx. References fE, fHitCalc, HitCalc::fHitFarAngle, HitCalc::fHitFarMomBalance, HitCalc::fHitLongEnergy, HitCalc::fHitLongEnergyRatio, HitCalc::fHitPeakAngle, HitCalc::fHitPeakMomBalance, HitCalc::fHitTotalEnergy, HitCalc::fHitTransCMEnergy, HitCalc::fHitTransCMEnergyRatio, HitCalc::fHitTransEnergy, HitCalc::fHitTransEnergyRatio, HitCalc::fHitTransLongEnergyRatio, fX, fY, fZ, MSG, s(), and VecFloat_t. Referenced by Analyze(). 00392 {
00393
00394 //Variable declaration
00395
00396 Float_t totalHitEnerg;
00397 Float_t arg;
00398 Float_t argTheta;
00399 Float_t eTrans;
00400 Float_t eTransX;
00401 Float_t eTransY;
00402 Float_t sumETransX;
00403 Float_t sumETransY;
00404 Float_t eLong;
00405 Float_t sTrans;
00406
00407 Float_t phi;
00408 Float_t phiMax;
00409 Float_t phiEMax;
00410 Float_t thetaLMax;
00411 Float_t thetaEMax;
00412 Float_t s;
00413 Float_t sE;
00414
00415 Float_t length;
00416 Float_t lengthMax;
00417 Float_t eMax;
00418
00419 Float_t sumOppSide;
00420 Float_t sumSameSide;
00421 Float_t ratioFar;
00422
00423 Float_t sumOppSideE;
00424 Float_t sumSameSideE;
00425 Float_t ratioEnerg;
00426
00427 Int_t nFar;
00428 Int_t nEMax;
00429
00430 Float_t zMax;
00431 Float_t zL;
00432
00433 VecFloat_t theta;
00434 theta.assign(fZ.size(),0.0);
00435
00437
00438 sumETransX=0.0;
00439 sumETransY=0.0;
00440 eLong=0.0;
00441 sTrans=0.0;
00442 length=0.0;
00443 lengthMax=0.0;
00444 eMax=0.0;
00445 phiMax=0.0;
00446 phiEMax=0.0;
00447 thetaLMax=0.0;
00448 thetaEMax=0.0;
00449 nFar=0;
00450 nEMax=0;
00451 zMax=0;
00452 totalHitEnerg=0.0;
00453
00454 Float_t epsilon=0.00001;
00455
00456 //Preliminary Calc of Total Hit Energy and Hit with highest z
00457 for (UInt_t n = 0; n < fZ.size(); n++){
00458 if(fZ.at(n) >= 0.0){ //Get rid of Hits behind the vertex
00459 totalHitEnerg += fE.at(n);
00460 if(fZ.at(n) > zMax) zMax = fZ.at(n);
00461 }
00462 }
00463
00464
00465 for (UInt_t n = 0; n < fZ.size(); n++){ //BIG Loop over all 3D Hits
00466 if(fZ.at(n) >= 0.0){ //Get rid of Hits behind the vertex
00467 if(fZ.at(n) == 0){ fZ.at(n) += epsilon; } //Protect against nans
00468
00469 arg = fX.at(n)*fX.at(n)+fY.at(n)*fY.at(n);
00470 argTheta = TMath::Sqrt(arg)/fZ.at(n);
00471 theta.at(n) = TMath::ATan(argTheta);
00472 phi = TMath::ACos(fX.at(n)/(argTheta*fZ.at(n) + epsilon));
00473
00474 MSG("HitCalcAna",Msg::kDebug)<< "Intermediate VarCalc Monitor "
00475 << endl
00476 << "arg = " << arg << endl
00477 << "argTheta = " << argTheta
00478 << endl
00479 << "fE[" << n << "] = " << fE.at(n)
00480 << endl
00481 << "fX = " << fX.at(n) << endl
00482 << "fY = " << fY.at(n) << endl
00483 << "fZ = " << fZ.at(n) << endl
00484 << "theta = " << theta.at(n)
00485 << endl << endl;
00486
00487 eTransX = fE.at(n)*fX.at(n)*TMath::Cos(theta.at(n))/fZ.at(n);
00488 eTransY = fE.at(n)*fY.at(n)*TMath::Cos(theta.at(n))/fZ.at(n);
00489 sumETransX += eTransX;
00490 sumETransY += eTransY;
00491
00492 sTrans += (eTransX*eTransX+eTransY*eTransY); //Get square of the
00493 //transverse energy
00494
00495 eLong += fE.at(n)*TMath::Cos(theta.at(n)); //Get Longitudinal
00496 //Energy
00497
00498 MSG("HitCalcAna",Msg::kDebug)<< "Intermediate VarCalc Monitor "
00499 << endl
00500 << "sumETransX[" << n << "] = "
00501 << sumETransX << endl
00502 << "sumETransY[" << n << "] = "
00503 << sumETransY << endl
00504 << "eLong = " << eLong
00505 << endl << endl;
00506
00507 //Calculate Farthest Hit
00508 length = TMath::Sqrt(arg+fZ.at(n)*fZ.at(n));
00509 if(length > lengthMax){
00510 if ( fZ.at(n) > 0.8*zMax
00511 && fZ.at(n) <= zMax){ //Try to avoid
00512 //far noise hits
00513 //causing very high
00514 //theta values.
00515 lengthMax = length;
00516 nFar = n;
00517 zL = fZ.at(n);
00518 phiMax = phi;
00519 thetaLMax = theta.at(n);
00520 }
00521 }
00522
00523 //Calculate Highest Energy Hit:
00524 if(fE.at(n) > eMax){
00525 eMax = fE.at(n);
00526 nEMax = n;
00527 phiEMax = phi;
00528 thetaEMax = theta.at(n);
00529 }
00530 }
00531 }
00532
00533 eTrans = TMath::Sqrt(sumETransX*sumETransX
00534 +sumETransY*sumETransY); //Transverse
00535 //energy
00536
00537 sumSameSide = 0.0;
00538 sumOppSide = 0.0;
00539 sumSameSideE = 0.0;
00540 sumOppSideE = 0.0;
00541
00542 for (UInt_t n=0; n < fZ.size(); n++){
00543 if(fZ.at(n) >= 0.0){
00544 //Calculate Ratio Far variable
00545 s = fX.at(n)*TMath::Cos(phiMax)+fY.at(n)*TMath::Sin(phiMax);
00546
00547 if(s >= 0.0) sumSameSide += fE.at(n)*TMath::Sin(theta.at(n));
00548 else sumOppSide += fE.at(n)*TMath::Sin(theta.at(n));
00549
00550 //Calculate Ratio variable for showers.
00551 sE = fX.at(n)*TMath::Cos(phiEMax)+fY.at(n)*TMath::Sin(phiEMax);
00552
00553 if(sE >= 0.0) sumSameSideE += fE.at(n)*TMath::Sin(theta.at(n));
00554 else sumOppSideE += fE.at(n)*TMath::Sin(theta.at(n));
00555
00556 }
00557 }
00558
00559 if (sumSameSide > 0.0){
00560 ratioFar = sumOppSide/sumSameSide;
00561 }else ratioFar = ANtpDefVal::kFloat;
00562
00563 if (sumSameSideE > 0.0){
00564 ratioEnerg = sumOppSideE/sumSameSideE;
00565 }else ratioEnerg = ANtpDefVal::kFloat;
00566
00567 if (ratioFar > 50.0) {ratioFar = ANtpDefVal::kFloat;}
00568 if (ratioEnerg > 50.0) {ratioEnerg = ANtpDefVal::kFloat;}
00569
00570 MSG("HitCalcAna",Msg::kInfo)<< "Calculated Hit variables " << endl
00571 << "eTotal = " << totalHitEnerg
00572 << " eTrans = " << eTrans
00573 << " eLong = " << eLong
00574 << " sTrans = " << sTrans
00575 << " farMomBalance = " << ratioFar
00576 << " peakMomBalance = " << ratioEnerg
00577 << " farAngle = " << thetaLMax
00578 << " peakAngle = " << thetaEMax
00579 << " zMax = " << zMax << endl;
00580
00581 Float_t transLongRatio=ANtpDefVal::kFloat;
00582 if(eLong > 0.0){transLongRatio=eTrans/eLong;}
00583
00584 //Fill the variables
00585 if (totalHitEnerg > 0.0){
00586 fHitCalc.fHitTotalEnergy = totalHitEnerg;
00587 fHitCalc.fHitTransEnergy = eTrans;
00588 fHitCalc.fHitLongEnergy = eLong;
00589 fHitCalc.fHitTransCMEnergy = sTrans;
00590 fHitCalc.fHitTransEnergyRatio = eTrans/totalHitEnerg;
00591 fHitCalc.fHitLongEnergyRatio = eLong/totalHitEnerg;
00592 fHitCalc.fHitTransLongEnergyRatio = transLongRatio;
00593 fHitCalc.fHitTransCMEnergyRatio = sTrans/totalHitEnerg;
00594 fHitCalc.fHitFarMomBalance = ratioFar;
00595 fHitCalc.fHitPeakMomBalance = ratioEnerg;
00596 fHitCalc.fHitFarAngle = thetaLMax;
00597 fHitCalc.fHitPeakAngle = thetaEMax;
00598 }
00599 else {
00600 MSG("HitCalcAna",Msg::kInfo)<< "No hits found in Event" << endl;
00601 fHitCalc.fHitTotalEnergy = ANtpDefVal::kFloat;
00602 fHitCalc.fHitTransEnergy = ANtpDefVal::kFloat;
00603 fHitCalc.fHitLongEnergy = ANtpDefVal::kFloat;
00604 fHitCalc.fHitTransCMEnergy = ANtpDefVal::kFloat;
00605 fHitCalc.fHitTransEnergyRatio = ANtpDefVal::kFloat;
00606 fHitCalc.fHitLongEnergyRatio = ANtpDefVal::kFloat;
00607 fHitCalc.fHitTransLongEnergyRatio = ANtpDefVal::kFloat;
00608 fHitCalc.fHitTransCMEnergyRatio = ANtpDefVal::kFloat;
00609 fHitCalc.fHitFarMomBalance = ANtpDefVal::kFloat;
00610 fHitCalc.fHitPeakMomBalance = ANtpDefVal::kFloat;
00611 fHitCalc.fHitFarAngle = ANtpDefVal::kFloat;
00612 fHitCalc.fHitPeakAngle = ANtpDefVal::kFloat;
00613
00614 }
00615 }
|
|
|
Definition at line 57 of file HitCalcAna.h. Referenced by ComputeHits(), Get3DHit(), and VarCalc(). |
|
|
Definition at line 71 of file HitCalcAna.h. Referenced by VarCalc(). |
|
|
Definition at line 58 of file HitCalcAna.h. |
|
|
Definition at line 54 of file HitCalcAna.h. Referenced by ComputeHits(), Get3DHit(), and VarCalc(). |
|
|
Definition at line 55 of file HitCalcAna.h. Referenced by ComputeHits(), Get3DHit(), and VarCalc(). |
|
|
Definition at line 56 of file HitCalcAna.h. Referenced by ComputeHits(), Get3DHit(), and VarCalc(). |
1.3.9.1