Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages
NueAna

AngClusterFitAna.cxx File Reference

#include <iostream>
#include <algorithm>
#include <TCanvas.h>
#include <TRotation.h>
#include <TFile.h>
#include <TH2F.h>
#include <TH1F.h>
#include <TH3F.h>
#include <TPolyMarker3D.h>
#include <TPolyLine3D.h>
#include "StandardNtuple/NtpStRecord.h"
#include "CandNtupleSR/NtpSRRecord.h"
#include "CandNtupleSR/NtpSREvent.h"
#include "CandNtupleSR/NtpSRTrack.h"
#include "CandNtupleSR/NtpSRStrip.h"
#include "MessageService/MsgService.h"
#include "NueAna/AngClusterFitAna.h"
#include "AnalysisNtuples/ANtpDefaultValue.h"

Go to the source code of this file.

Functions

 CVSID ("$Id: AngClusterFitAna.cxx,v 1.21 2007/03/01 16:38:49 rhatcher Exp $")
Double_t shwLongFunc (Double_t *x, Double_t *par)
Double_t shwTransFunc (Double_t *x, Double_t *par)

Variables

const Float_t kDetectorPitch = 0.05949
const Float_t kStripWidth = 0.041
const Float_t kRadLengthToPlane = 1.54
const Float_t kMolRadDet = 0.0391
const Float_t kMeuToGeV = 1.0

Function Documentation

CVSID (  "$Id: AngClusterFitAna.  cxx,
v 1.21 2007/03/01 16:38:49 rhatcher Exp $" 
) 
 

Id
AngClusterFitAna.cxx,v 1.21 2007/03/01 16:38:49 rhatcher Exp

class AngClusterFitAna

NueAna package

Purpose: Fit a EM shower profile to the 3D hits contained in the primary cluster determined by AngClusterAna..

Author: Alex Sousa <asousa@minos.phy.tufts.edu>

Created on: Fri May 06 2005

Double_t shwLongFunc (  Double_t *  x,
Double_t *  par
)  [static]
 

Definition at line 135 of file AngClusterFitAna.cxx.

References kRadLengthToPlane.

Referenced by AngClusterFitAna::FitShower().

00136 {
00137 
00138     // Longitudinal profile to be fitted to em showers 
00139     //xx=Distance from the vertex in units of X0 
00140     //par[0]=a, par[1]=b, par[2]=E0
00141    
00142     Float_t xx=x[0];
00143     
00144     Double_t lnf = TMath::Log(kRadLengthToPlane*par[2]*par[1])+(par[0]-1)
00145         *TMath::Log(xx*par[1])-par[1]*xx-TMath::LnGamma(par[0]);
00146     
00147     Double_t f = TMath::Exp(lnf);
00148 
00149     return f;
00150 }

Double_t shwTransFunc (  Double_t *  x,
Double_t *  par
)  [static]
 

Definition at line 152 of file AngClusterFitAna.cxx.

Referenced by AngClusterFitAna::FitShower().

00153 {
00154     //Transverse fit profile
00155     //x[0]=Radial distance from shower axis in units of Moliere radius
00156     //par[0]=Lambda1 attenuation length for high energy core
00157     //par[1]=Lambda2 attenuation length for shower halo
00158     //par[2]=Relative weight between lambda1 and lambda2
00159     //par[3]=E0 Energy normalization
00160 
00161     Float_t xx=x[0];
00162     Double_t f=par[3]*(TMath::Exp(-TMath::Sqrt(xx/par[0]))
00163                        +par[2]*TMath::Exp(-xx/par[1]));    
00164     return f;
00165 }

Variable Documentation

const Float_t kDetectorPitch = 0.05949
 

Definition at line 40 of file AngClusterFitAna.cxx.

const Float_t kMeuToGeV = 1.0
 

Definition at line 49 of file AngClusterFitAna.cxx.

const Float_t kMolRadDet = 0.0391
 

Definition at line 46 of file AngClusterFitAna.cxx.

const Float_t kRadLengthToPlane = 1.54
 

Definition at line 44 of file AngClusterFitAna.cxx.

Referenced by AngClusterFitAna::Analyze(), and shwLongFunc().

const Float_t kStripWidth = 0.041
 

Definition at line 42 of file AngClusterFitAna.cxx.

Generated on Mon Feb 15 11:07:52 2010 for loon by  doxygen 1.3.9.1